private void Initialize()
{
// Cache the .dbf header
_dbfHeader = new DbaseFileHeader();
using (BinaryReader dbfReader = new BinaryReader(File.OpenRead(Path.Combine(Path.GetDirectoryName(_path), Path.GetFileNameWithoutExtension(_path) + ".dbf"))))
{
_dbfHeader.ReadHeader(dbfReader);
_dbfHeaderOffset = dbfReader.BaseStream.Position;
}
// Need to make one pass over the geometries and pull out the bounding boxes
_spatialIndex = new com.vividsolutions.jts.index.strtree.STRtree();
_extents = new Envelope();
using (BigEndianBinaryReader shpReader = new BigEndianBinaryReader(File.OpenRead(_path)))
using (ShapefileIndexReader shxReader = new ShapefileIndexReader(Path.Combine(Path.GetDirectoryName(_path), Path.GetFileNameWithoutExtension(_path) + ".shx")))
{
// Get the shape type
_type = new ShapefileHeader(shpReader).ShapeType;
while (shxReader.Read())
{
int offset = shxReader.GetOffest();
int length = shxReader.GetLength();
// Move to the start of geometry
shpReader.BaseStream.Position = offset * 2;
double xMin;
double yMin;
double xMax;
double yMax;
int recordNumber = shpReader.ReadIntBE();
int contentLength = shpReader.ReadIntBE();
// Read shape type
int type = shpReader.ReadInt32();
if (type != 1)
{
xMin = shpReader.ReadDouble();
yMin = shpReader.ReadDouble();
xMax = shpReader.ReadDouble();
yMax = shpReader.ReadDouble();
}
else
{
// Point - read x and y
xMin = shpReader.ReadDouble();
yMin = shpReader.ReadDouble();
xMax = yMin;
yMax = yMin;
}
// Build the envelope
Envelope extents = new Envelope(xMin, xMax, yMin, yMax);
// Add to total extents
_extents.expandToInclude(extents);
// Insert the index of the record into the spatial index
_spatialIndex.insert(extents, new ShapefileRecordPointer(recordNumber, offset * 2, contentLength, (int)_dbfHeaderOffset + (_dbfHeader.RecordLength * (recordNumber - 1))));
}
// Build the index once
_spatialIndex.build();
}
}
private void Initialize()
{
// Cache the .dbf header
_dbfHeader = new DbaseFileHeader();
using (BinaryReader dbfReader = new BinaryReader(File.OpenRead(Path.Combine(Path.GetDirectoryName(_path), Path.GetFileNameWithoutExtension(_path) + ".dbf"))))
{
_dbfHeader.ReadHeader(dbfReader);
_dbfHeaderOffset = dbfReader.BaseStream.Position;
}
// Need to make one pass over the geometries and pull out the bounding boxes
_spatialIndex = new com.vividsolutions.jts.index.strtree.STRtree();
_extents = new Envelope();
using (BigEndianBinaryReader shpReader = new BigEndianBinaryReader(File.OpenRead(_path)))
using (ShapefileIndexReader shxReader = new ShapefileIndexReader(Path.Combine(Path.GetDirectoryName(_path), Path.GetFileNameWithoutExtension(_path) + ".shx")))
{
// Get the shape type
_type = new ShapefileHeader(shpReader).ShapeType;
while (shxReader.Read())
{
int offset = shxReader.GetOffest();
int length = shxReader.GetLength();
// Move to the start of geometry
shpReader.BaseStream.Position = offset * 2;
double xMin;
double yMin;
double xMax;
double yMax;
int recordNumber = shpReader.ReadIntBE();
int contentLength = shpReader.ReadIntBE();
// Read shape type
int type = shpReader.ReadInt32();
if (type != 1)
{
xMin = shpReader.ReadDouble();
yMin = shpReader.ReadDouble();
xMax = shpReader.ReadDouble();
yMax = shpReader.ReadDouble();
}
else
{
// Point - read x and y
xMin = shpReader.ReadDouble();
yMin = shpReader.ReadDouble();
xMax = yMin;
yMax = yMin;
}
// Build the envelope
Envelope extents = new Envelope(xMin, xMax, yMin, yMax);
// Add to total extents
_extents.expandToInclude(extents);
// Insert the index of the record into the spatial index
_spatialIndex.insert(extents, new ShapefileRecordPointer(recordNumber, offset * 2, contentLength, (int)_dbfHeaderOffset + (_dbfHeader.RecordLength * (recordNumber - 1))));
}
// Build the index once
_spatialIndex.build();
}
}
