This page hosts a formal specification of Shapefile main file using Kaitai Struct. This specification can be automatically translated into a variety of programming languages to get a parsing library.
All parsing code for Python generated by Kaitai Struct depends on the Python runtime library. You have to install it before you can parse data.
The Python runtime library can be installed from PyPI:
python3 -m pip install kaitaistruct
Parse a local file and get structure in memory:
data = ShapefileMain.from_file("path/to/local/file.shp")
Or parse structure from a bytes:
from kaitaistruct import KaitaiStream, BytesIO
raw = b"\x00\x01\x02..."
data = ShapefileMain(KaitaiStream(BytesIO(raw)))
After that, one can get various attributes from the structure by invoking getter methods like:
data.records # => the size of this section of the file in bytes must equal (header.file_length * 2) - 100
# This is a generated file! Please edit source .ksy file and use kaitai-struct-compiler to rebuild
import kaitaistruct
from kaitaistruct import KaitaiStruct, KaitaiStream, BytesIO
from enum import Enum
if getattr(kaitaistruct, 'API_VERSION', (0, 9)) < (0, 9):
raise Exception("Incompatible Kaitai Struct Python API: 0.9 or later is required, but you have %s" % (kaitaistruct.__version__))
class ShapefileMain(KaitaiStruct):
class ShapeType(Enum):
null_shape = 0
point = 1
poly_line = 3
polygon = 5
multi_point = 8
point_z = 11
poly_line_z = 13
polygon_z = 15
multi_point_z = 18
point_m = 21
poly_line_m = 23
polygon_m = 25
multi_point_m = 28
multi_patch = 31
class PartType(Enum):
triangle_strip = 0
triangle_fan = 1
outer_ring = 2
inner_ring = 3
first_ring = 4
ring = 5
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.header = ShapefileMain.FileHeader(self._io, self, self._root)
self.records = []
i = 0
while not self._io.is_eof():
self.records.append(ShapefileMain.Record(self._io, self, self._root))
i += 1
class MultiPointM(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_points = self._io.read_s4le()
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class BoundingBoxXYZM(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.x = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.y = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.z = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m = ShapefileMain.BoundsMinMax(self._io, self, self._root)
class Point(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.x = self._io.read_f8le()
self.y = self._io.read_f8le()
class Polygon(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
class BoundsMinMax(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.min = self._io.read_f8le()
self.max = self._io.read_f8le()
class PolyLine(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
class MultiPointZ(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_points = self._io.read_s4le()
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.z_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.z_values = []
for i in range(self.number_of_points):
self.z_values.append(self._io.read_f8le())
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class PolyLineZ(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.z_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.z_values = []
for i in range(self.number_of_points):
self.z_values.append(self._io.read_f8le())
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class PolygonZ(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.z_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.z_values = []
for i in range(self.number_of_points):
self.z_values.append(self._io.read_f8le())
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class BoundingBoxXY(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.x = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.y = ShapefileMain.BoundsMinMax(self._io, self, self._root)
class PointM(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.x = self._io.read_f8le()
self.y = self._io.read_f8le()
self.m = self._io.read_f8le()
class PolygonM(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class RecordHeader(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.record_number = self._io.read_s4be()
self.content_length = self._io.read_s4be()
class MultiPoint(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_points = self._io.read_s4le()
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
class FileHeader(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.file_code = self._io.read_bytes(4)
if not self.file_code == b"\x00\x00\x27\x0A":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x27\x0A", self.file_code, self._io, u"/types/file_header/seq/0")
self.unused_field_1 = self._io.read_bytes(4)
if not self.unused_field_1 == b"\x00\x00\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x00\x00", self.unused_field_1, self._io, u"/types/file_header/seq/1")
self.unused_field_2 = self._io.read_bytes(4)
if not self.unused_field_2 == b"\x00\x00\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x00\x00", self.unused_field_2, self._io, u"/types/file_header/seq/2")
self.unused_field_3 = self._io.read_bytes(4)
if not self.unused_field_3 == b"\x00\x00\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x00\x00", self.unused_field_3, self._io, u"/types/file_header/seq/3")
self.unused_field_4 = self._io.read_bytes(4)
if not self.unused_field_4 == b"\x00\x00\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x00\x00", self.unused_field_4, self._io, u"/types/file_header/seq/4")
self.unused_field_5 = self._io.read_bytes(4)
if not self.unused_field_5 == b"\x00\x00\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\x00\x00\x00\x00", self.unused_field_5, self._io, u"/types/file_header/seq/5")
self.file_length = self._io.read_s4be()
self.version = self._io.read_bytes(4)
if not self.version == b"\xE8\x03\x00\x00":
raise kaitaistruct.ValidationNotEqualError(b"\xE8\x03\x00\x00", self.version, self._io, u"/types/file_header/seq/7")
self.shape_type = KaitaiStream.resolve_enum(ShapefileMain.ShapeType, self._io.read_s4le())
self.bounding_box = ShapefileMain.BoundingBoxXYZM(self._io, self, self._root)
class PointZ(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.x = self._io.read_f8le()
self.y = self._io.read_f8le()
self.z = self._io.read_f8le()
self.m = self._io.read_f8le()
class Record(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.header = ShapefileMain.RecordHeader(self._io, self, self._root)
self.contents = ShapefileMain.RecordContents(self._io, self, self._root)
class RecordContents(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.shape_type = KaitaiStream.resolve_enum(ShapefileMain.ShapeType, self._io.read_s4le())
if self.shape_type != ShapefileMain.ShapeType.null_shape:
_on = self.shape_type
if _on == ShapefileMain.ShapeType.poly_line_z:
self.shape_parameters = ShapefileMain.PolyLineZ(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.multi_patch:
self.shape_parameters = ShapefileMain.MultiPatch(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.poly_line_m:
self.shape_parameters = ShapefileMain.PolyLineM(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.polygon:
self.shape_parameters = ShapefileMain.Polygon(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.polygon_z:
self.shape_parameters = ShapefileMain.PolygonZ(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.point_z:
self.shape_parameters = ShapefileMain.PointZ(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.poly_line:
self.shape_parameters = ShapefileMain.PolyLine(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.point_m:
self.shape_parameters = ShapefileMain.PointM(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.polygon_m:
self.shape_parameters = ShapefileMain.PolygonM(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.multi_point:
self.shape_parameters = ShapefileMain.MultiPoint(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.point:
self.shape_parameters = ShapefileMain.Point(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.multi_point_m:
self.shape_parameters = ShapefileMain.MultiPointM(self._io, self, self._root)
elif _on == ShapefileMain.ShapeType.multi_point_z:
self.shape_parameters = ShapefileMain.MultiPointZ(self._io, self, self._root)
class MultiPatch(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.part_types = []
for i in range(self.number_of_parts):
self.part_types.append(KaitaiStream.resolve_enum(ShapefileMain.PartType, self._io.read_s4le()))
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.z_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.z_values = []
for i in range(self.number_of_points):
self.z_values.append(self._io.read_f8le())
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())
class PolyLineM(KaitaiStruct):
def __init__(self, _io, _parent=None, _root=None):
self._io = _io
self._parent = _parent
self._root = _root if _root else self
self._read()
def _read(self):
self.bounding_box = ShapefileMain.BoundingBoxXY(self._io, self, self._root)
self.number_of_parts = self._io.read_s4le()
self.number_of_points = self._io.read_s4le()
self.parts = []
for i in range(self.number_of_parts):
self.parts.append(self._io.read_s4le())
self.points = []
for i in range(self.number_of_points):
self.points.append(ShapefileMain.Point(self._io, self, self._root))
self.m_range = ShapefileMain.BoundsMinMax(self._io, self, self._root)
self.m_values = []
for i in range(self.number_of_points):
self.m_values.append(self._io.read_f8le())