我已经回答了这个问题，并将问题的标题从如何使用带有 c++ 的 ogr 库访问多边形层中的顶点更改为如何使用带有 c++的 ogr 库在 ESRI shapefile 中提取几何的顶点。该代码可用于使用 OpenGL 显示形状。
任何使它更好的建议将不胜感激。
例如，您认为我应该使用图层的投影来增强 OpenGL 中的视图吗？如何使用投影并将其应用于 OpenGL 窗口？
到目前为止，我已经尝试了这些类以提取多边形的顶点：
getNumPoints()方法时，我的图层中的所有多边形特征都得到 1。getNumPoints()方法时，我的图层中的所有多边形特征都得到 1。为了在 OpenGL 窗口中绘制多边形，我应该怎么做才能获得多边形 shapefile 中的顶点数和顶点坐标。
我找到了答案：
你应该关心一些关于 ESRI shapefile 的信息。
点层：只能具有wkbpoint类型的特征。因此，您可以定义一个名为MyPoint2D的数据结构，并将该层的坐标存储在具有 MyPoint2D类型成员的向量中。
多点层：只能具有wkbMultipoint类型的特征。因此，您可以将每个特征存储在具有MyPoint2D类型成员的向量中，并且此类向量的集合将保存图层中所有特征的坐标。
折线层：可以具有wkbLineString和wkbMultiLineString类型的特征。所以你可以有一个数据结构MyLine2D，它是一个具有MyPoint2D成员的向量，一个具有 MyLine2D 类型成员的向量将保存 LineFeature，最后层中的特征将存储在具有LineFeature类型成员的向量中。
多边形层：关于多边形，您应该知道每个多边形都由一个外环和几个内环组成（这是我真正的问题，导致我提出这个问题）。而且该层还可以具有 wkbPolygon和wkbMultiPolygon类型的特征。
所以首先我们应该定义这些数据结构：
//data structure for points
typedef struct MyPoint2D
{
double dX;
double dY;
}MyPoint2D;
//Holds Coordinates of Point Shapefile
vector<MyPoint2D>PointLayer;
//data structure for a multipoint feature
typedef struct MultipointFeature
{
vector<MyPoint2D>PointsOfFeature;
}MultipointFeature;
//Holds Coordinates of multiPoint Shapefile
vector<MultipointFeature> MultipointLayer;
//data structure for lines
typedef struct MyLine2D
{
vector<MyPoint2D> LineString;
}MyLine2D;
//data structure for a line feature
typedef struct LineFeature
{
vector<MyLine2D>LinesOfFeature;
}LineFeature;
//Holds Coordinates of Line Shapefile
vector<LineFeature> LineLayer;
//data structure for rings
typedef struct MyRing2D
{
vector<MyPoint2D> RingString;
bool IsClockwised;
}MyRing2D;
//data structure for polygons
typedef struct MyPolygon2D
{
vector<MyRing2D>Polygon;
}MyPolygon2D;
//data structure for a polygon feature
typedef struct PolygonFeature
{
vector<MyPolygon2D>PolygonsOfFeature;
}PolygonFeature;
//Holds Coordinates of Polygon Shapefile
vector<PolygonFeature> PolygonLayer;
//data structure to hold bounding box
typedef struct SBoundingBox
{
float fMaxX;
float fMaxY;
float fMinX;
float fMinY;
}SBoundingBox;
//Bounding Box of Shapefile
SBoundingBox sBoundingBox;
这是整个代码，我用它来提取 ESRI 的 shapefile 中的顶点坐标：
void OpenShapeFile(char* filename)
{
OGRErr error;
OGRDataSource *poDataSource;
poDataSource = OGRSFDriverRegistrar::Open(filename,false);
OGRLayer *poLayer;
poLayer = poDataSource ->GetLayer(0);
OGREnvelope Envelope;
error = poLayer ->GetExtent(&Envelope,true);
sBoundingBox.fMaxX = Envelope.MaxX;
sBoundingBox.fMaxY = Envelope.MaxY;
sBoundingBox.fMinX = Envelope.MinX;
sBoundingBox.fMinY = Envelope.MinY;
OGRwkbGeometryType LayerGeometryType = poLayer ->GetGeomType();
int NumberOfFeatures = poLayer ->GetFeatureCount(true);
poLayer ->ResetReading();
//Point Shapefile
if ( wkbFlatten ( LayerGeometryType ) == wkbPoint )
{
OGRFeature *poFeature;
for ( int i = 0; i < NumberOfFeatures; i++ )
{
poFeature = poLayer ->GetNextFeature();
OGRGeometry *poGeometry;
poGeometry = poFeature ->GetGeometryRef();
if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbPoint )
{
OGRPoint *poPoint = ( OGRPoint * )poGeometry;
MyPoint2D pt;
pt.dX = poPoint ->getX();
pt.dY = poPoint ->getY();
PointLayer.push_back(pt);
}
OGRFeature::DestroyFeature(poFeature);
}
}
//Multipoint Shapefile
if ( wkbFlatten ( LayerGeometryType ) == wkbMultiPoint )
{
OGRFeature *poFeature;
MultipointFeature MultiPoint;
for ( int i = 0; i < NumberOfFeatures; i++ )
{
poFeature = poLayer ->GetNextFeature();
OGRGeometry *poGeometry;
poGeometry = poFeature ->GetGeometryRef();
if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbMultiPoint )
{
OGRMultiPoint *poMultipoint = ( OGRMultiPoint * )poGeometry;
int NumberOfGeometries = poMultipoint ->getNumGeometries();
MultiPoint.PointsOfFeature.resize( NumberOfGeometries );
for ( int j = 0; j < NumberOfGeometries; j++ )
{
OGRGeometry *poPointGeometry = poMultipoint ->getGeometryRef(j);
OGRPoint *poPoint = ( OGRPoint * )poPointGeometry;
MyPoint2D pt;
pt.dX = poPoint ->getX();
pt.dY = poPoint ->getY();
MultiPoint.PointsOfFeature.at(j) = pt;
}
MultipointLayer.push_back(MultiPoint);
}
OGRFeature::DestroyFeature(poFeature);
}
}
//Polyline Shapefile
if ( wkbFlatten ( LayerGeometryType ) == wkbLineString )
{
OGRFeature *poFeature;
LineFeature Polyline;
OGRPoint ptTemp;
for ( int i = 0; i < NumberOfFeatures; i++ )
{
poFeature = poLayer ->GetNextFeature();
OGRGeometry *poGeometry;
poGeometry = poFeature ->GetGeometryRef();
if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbLineString )
{
OGRLineString *poLineString = ( OGRLineString * )poGeometry;
Polyline.LinesOfFeature.resize(1);
int NumberOfVertices = poLineString ->getNumPoints();
Polyline.LinesOfFeature.at(0).LineString.resize(NumberOfVertices);
for ( int k = 0; k < NumberOfVertices; k++ )
{
poLineString ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polyline.LinesOfFeature.at(0).LineString.at(k) = pt;
}
LineLayer.push_back(Polyline);
}
else if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbMultiLineString )
{
OGRMultiLineString *poMultiLineString = ( OGRMultiLineString * )poGeometry;
int NumberOfGeometries = poMultiLineString ->getNumGeometries();
Polyline.LinesOfFeature.resize(NumberOfGeometries);
for ( int j = 0; j < NumberOfGeometries; j++ )
{
OGRGeometry *poLineGeometry = poMultiLineString ->getGeometryRef(j);
OGRLineString *poLineString = ( OGRLineString * )poLineGeometry;
int NumberOfVertices = poLineString ->getNumPoints();
Polyline.LinesOfFeature.at(j).LineString.resize(NumberOfVertices);
for ( int k = 0; k < NumberOfVertices; k++ )
{
poLineString ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polyline.LinesOfFeature.at(j).LineString.at(k) = pt;
}
}
LineLayer.push_back(Polyline);
}
OGRFeature::DestroyFeature(poFeature);
}
}
//Polygon Shapefile
if ( wkbFlatten ( LayerGeometryType ) == wkbPolygon )
{
OGRFeature *poFeature;
PolygonFeature Polygon;
OGRPoint ptTemp;
for ( int i = 0; i < NumberOfFeatures; i++ )
{
poFeature = poLayer ->GetNextFeature();
OGRGeometry *poGeometry;
poGeometry = poFeature ->GetGeometryRef();
if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbPolygon )
{
OGRPolygon *poPolygon = ( OGRPolygon * )poGeometry;
Polygon.PolygonsOfFeature.resize(1);
int NumberOfInnerRings = poPolygon ->getNumInteriorRings();
OGRLinearRing *poExteriorRing = poPolygon ->getExteriorRing();
Polygon.PolygonsOfFeature.at(0).Polygon.resize(NumberOfInnerRings+1);
Polygon.PolygonsOfFeature.at(0).Polygon.at(0).IsClockwised = poExteriorRing ->isClockwise();
int NumberOfExteriorRingVertices = poExteriorRing ->getNumPoints();
Polygon.PolygonsOfFeature.at(0).Polygon.at(0).RingString.resize(NumberOfExteriorRingVertices);
for ( int k = 0; k < NumberOfExteriorRingVertices; k++ )
{
poExteriorRing ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polygon.PolygonsOfFeature.at(0).Polygon.at(0).RingString.at(k) = pt;
}
for ( int h = 1; h <= NumberOfInnerRings; h++ )
{
OGRLinearRing *poInteriorRing = poPolygon ->getInteriorRing(h-1);
Polygon.PolygonsOfFeature.at(0).Polygon.at(h).IsClockwised = poInteriorRing ->isClockwise();
int NumberOfInteriorRingVertices = poInteriorRing ->getNumPoints();
Polygon.PolygonsOfFeature.at(0).Polygon.at(h).RingString.resize(NumberOfInteriorRingVertices);
for ( int k = 0; k < NumberOfInteriorRingVertices; k++ )
{
poInteriorRing ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polygon.PolygonsOfFeature.at(0).Polygon.at(h).RingString.at(k) = pt;
}
}
PolygonLayer.push_back(Polygon);
}
else if ( poGeometry != NULL && wkbFlatten ( poGeometry ->getGeometryType() ) == wkbMultiPolygon )
{
OGRMultiPolygon *poMultiPolygon = ( OGRMultiPolygon * )poGeometry;
int NumberOfGeometries = poMultiPolygon ->getNumGeometries();
Polygon.PolygonsOfFeature.resize(NumberOfGeometries);
for ( int j = 0; j < NumberOfGeometries; j++ )
{
OGRGeometry *poPolygonGeometry = poMultiPolygon ->getGeometryRef(j);
OGRPolygon *poPolygon = ( OGRPolygon * )poPolygonGeometry;
int NumberOfInnerRings = poPolygon ->getNumInteriorRings();
OGRLinearRing *poExteriorRing = poPolygon ->getExteriorRing();
Polygon.PolygonsOfFeature.at(j).Polygon.resize(NumberOfInnerRings+1);
Polygon.PolygonsOfFeature.at(j).Polygon.at(0).IsClockwised = poExteriorRing ->isClockwise();
int NumberOfExteriorRingVertices = poExteriorRing ->getNumPoints();
Polygon.PolygonsOfFeature.at(j).Polygon.at(0).RingString.resize(NumberOfExteriorRingVertices);
for ( int k = 0; k < NumberOfExteriorRingVertices; k++ )
{
poExteriorRing ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polygon.PolygonsOfFeature.at(j).Polygon.at(0).RingString.at(k) = pt;
}
for ( int h = 1; h <= NumberOfInnerRings; h++ )
{
OGRLinearRing *poInteriorRing = poPolygon ->getInteriorRing(h-1);
Polygon.PolygonsOfFeature.at(j).Polygon.at(h).IsClockwised = poInteriorRing ->isClockwise();
int NumberOfInteriorRingVertices = poInteriorRing ->getNumPoints();
Polygon.PolygonsOfFeature.at(j).Polygon.at(h).RingString.resize(NumberOfInteriorRingVertices);
for ( int k = 0; k < NumberOfInteriorRingVertices; k++ )
{
poInteriorRing ->getPoint(k,&ptTemp);
MyPoint2D pt;
pt.dX = ptTemp.getX();
pt.dY = ptTemp.getY();
Polygon.PolygonsOfFeature.at(j).Polygon.at(h).RingString.at(k) = pt;
}
}
}
PolygonLayer.push_back(Polygon);
}
}
OGRFeature::DestroyFeature(poFeature);
}
OGRDataSource::DestroyDataSource(poDataSource);
}
scott_f想要一个 C# 实现，将来可能有人也需要它，所以我将它作为答案发布。我很懒，不想使用嵌套列表和不同的数据结构重写整个事情，Rhino/Grasshopper因为这个实现是插件GH_Structure, IGH_GeometricGoo, IGH_Goo的一部分（是Rhino/GrasshopperAPI 的一部分）。但基本上IGH_Goo在这种情况下是指一个几何类bools，IGH_GeometricGoo它可以包含Polyline和Point3d类，并且GH_Structure是一个多维数据结构（列表列表），它通过.Append()方法“添加”项目，其中第一个参数是要添加的数据第二个是“路径”，因此.Append( item, new GH_Path( 0, 3 ) )将添加item到第 4 个索引处的第一个列表列表中。这可能是个好主意wkbFlatten，以便它在 wkbType 的维度不同时起作用，但我将把它留给读者，因为我的实现不需要它。请注意，layer是作为参数传递的，因此如果您需要文件中的所有数据，则应遍历所有层。
class WKBToRhinoGeoUtilities
{
private string input;
private int layerNumber;
private int EPSGCode;
private int InEPSGCode;
private OSGeo.OGR.Driver driver;
OSGeo.OGR.DataSource dataSource;
private GH_Path path;
private double[] pointList = { 0, 0, 0 };
private long numberOfFeatures;
private int layerCount;
private OSGeo.OGR.Layer layer;
private OSGeo.OGR.wkbGeometryType wkbGeoType;
private int currentProjection;
private OSGeo.OSR.CoordinateTransformation transform;
private OSGeo.OGR.Feature feature;
private OSGeo.OGR.Geometry geo, ring;
public string debug = "Working!";
public WKBToRhinoGeoUtilities( string input, int layerNumber, int EPSGCode,
int InEPSGCode, OSGeo.OGR.DataSource dataSource )
{
this.input = input;
this.layerNumber = layerNumber;
this.EPSGCode = EPSGCode;
this.InEPSGCode = InEPSGCode;
this.dataSource = dataSource;
}
public string Setup()
{
GetDataSourceAttributes();
IdentifyEPSG();
return debug;
}
public GH_Structure<IGH_GeometricGoo> Run( ref GH_Structure<IGH_Goo> cullPattern )
{
switch( wkbGeoType )
{
case OSGeo.OGR.wkbGeometryType.wkbUnknown:
return WKBUnknown( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbPolygon:
return WKBPolygon( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbLineString:
return WKBLineString( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbPoint:
return WKBPoint( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbPoint25D:
return WKBPoint( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbMultiPolygon:
return WKBMultiPolygon( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbMultiLineString:
return WKBMultiLineString( ref cullPattern );
case OSGeo.OGR.wkbGeometryType.wkbGeometryCollection:
return WKBGeometryCollection( ref cullPattern );
default:
return null;
}
}
private void GetDataSourceAttributes()
{
if( dataSource != null )
{
layerCount = dataSource.GetLayerCount();
layer = dataSource.GetLayerByIndex( layerNumber % layerCount );
wkbGeoType = layer.GetGeomType();
numberOfFeatures = unchecked( ( int ) layer.GetFeatureCount( 0 ) );
}
else
{
debug = "Data Source is null!";
}
}
private void IdentifyEPSG()
{
// Source projections, try to get it automatically otherwise resort to user input.
OSGeo.OSR.SpatialReference source = null;
int identifyEPSG = 0;
try
{
source = layer.GetSpatialRef();
identifyEPSG = source.AutoIdentifyEPSG();
}
catch
{
debug = "The file's EPSG code couldn't be found automatically, it will default to EPSG:3857, make sure you input the correct one!";
source = new OSGeo.OSR.SpatialReference( "" );
identifyEPSG = source.ImportFromEPSG( InEPSGCode );
}
currentProjection = Convert.ToInt32( source.GetAttrValue( "AUTHORITY", 1 ) );
// Reprojections.
OSGeo.OSR.SpatialReference target = new OSGeo.OSR.SpatialReference( "" );
if( currentProjection != EPSGCode )
{
target.ImportFromEPSG( EPSGCode );
transform = new OSGeo.OSR.CoordinateTransformation( source, target );
}
}
private GH_Structure<IGH_GeometricGoo> WKBUnknown( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
double[] pointList = { 0, 0, 0 };
int counter = 0;
do
{
path = new GH_Path( counter );
feature = layer.GetNextFeature();
bool pattern = false;
if( feature != null )
{
OSGeo.OGR.Geometry geo = feature.GetGeometryRef();
if( geo != null )
{
if( currentProjection != EPSGCode )
{
geo.Transform( transform );
}
pattern = true;
int pointCount = geo.GetPointCount();
IGH_GeometricGoo geoGoo = null;
Point3d[] tempPointArray = new Point3d[pointCount];
for( int i = 0; i < pointCount; ++i )
{
geo.GetPoint( i, pointList );
tempPointArray[i] = new Point3d( pointList[0], pointList[1], pointList[2] );
}
if( tempPointArray.Length > 2 )
{
Polyline tempPoly = new Polyline( tempPointArray );
geoGoo = GH_Convert.ToGeometricGoo( tempPoly );
geoOutput.Append( geoGoo, path );
}
else if( tempPointArray.Length > 0 )
{
geoGoo = GH_Convert.ToGeometricGoo( tempPointArray[0] );
geoOutput.Append( geoGoo, path );
}
}
}
IGH_Goo patternGoo = GH_Convert.ToGoo( pattern );
cullPattern.Append( patternGoo, path );
counter++;
} while( feature != null );
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBPolygon( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
bool pattern = false;
if( feature != null )
{
geo = feature.GetGeometryRef();
if( geo != null )
{
IGH_GeometricGoo geoGoo = WKBPolygonSingle( geo, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
}
}
else
{
feature = layer.GetNextFeature();
int counter = 0;
while( feature != null )
{
path = new GH_Path( counter );
geo = feature.GetGeometryRef();
bool pattern = false;
if( geo != null )
{
IGH_GeometricGoo geoGoo = WKBPolygonSingle( geo, ref pattern );
geoOutput.Append( geoGoo, path );
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
feature = layer.GetNextFeature();
counter++;
}
}
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBLineString( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
bool pattern = false;
if( feature != null )
{
geo = feature.GetGeometryRef();
if( geo != null )
{
IGH_GeometricGoo geoGoo = WKBLineStringSingle( geo, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
}
}
else
{
int counter = 0;
feature = layer.GetNextFeature();
while( feature != null )
{
path = new GH_Path( counter );
geo = feature.GetGeometryRef();
bool pattern = false;
if( geo != null )
{
IGH_GeometricGoo geoGoo = WKBLineStringSingle( geo, ref pattern );
geoOutput.Append( geoGoo, path );
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
feature = layer.GetNextFeature();
counter++;
}
}
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBPoint( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
bool tempPattern = feature != null ? true : false;
IGH_Goo gooPattern = GH_Convert.ToGoo( tempPattern );
cullPattern.Append( gooPattern, path );
if( feature != null )
{
geo = feature.GetGeometryRef();
IGH_GeometricGoo geoGoo = WKBPointSingle( geo );
geoOutput.Append( geoGoo, path );
}
}
}
else
{
int counter = 0;
feature = layer.GetNextFeature();
while( feature != null )
{
path = new GH_Path( counter );
bool tempPattern = feature != null ? true : false;
IGH_Goo gooPattern = GH_Convert.ToGoo( tempPattern );
cullPattern.Append( gooPattern, path );
geo = feature.GetGeometryRef();
IGH_GeometricGoo geoGoo = WKBPointSingle( geo );
geoOutput.Append( geoGoo, path );
feature = layer.GetNextFeature();
counter++;
}
}
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBMultiPolygon( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
bool pattern = false;
if( feature != null )
{
geo = feature.GetGeometryRef();
for( int j = 0; j < geo.GetGeometryCount(); ++j )
{
OSGeo.OGR.Geometry inGeo = geo.GetGeometryRef(j);
IGH_GeometricGoo geoGoo = WKBPolygonSingle( inGeo, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
}
}
else
{
int counter = 0;
feature = layer.GetNextFeature();
while( feature != null )
{
bool pattern = false;
path = new GH_Path( counter );
geo = feature.GetGeometryRef();
for( int j = 0; j < geo.GetGeometryCount(); ++j )
{
OSGeo.OGR.Geometry inGeo = geo.GetGeometryRef(j);
if( inGeo != null )
{
IGH_GeometricGoo geoGoo = WKBPolygonSingle( inGeo, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
counter++;
feature = layer.GetNextFeature();
}
}
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBMultiLineString( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
bool pattern = false;
if( feature != null )
{
geo = feature.GetGeometryRef();
if( geo != null )
{
for( int j = 0; j < geo.GetGeometryCount(); ++j )
{
ring = geo.GetGeometryRef(j);
IGH_GeometricGoo geoGoo = WKBLineStringSingle( ring, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
}
}
else
{
feature = layer.GetNextFeature();
int counter = 0;
while( feature != null )
{
path = new GH_Path( counter );
geo = feature.GetGeometryRef();
bool pattern = false;
if( geo != null )
{
for( int j = 0; j < geo.GetGeometryCount(); ++j )
{
ring = geo.GetGeometryRef(j);
IGH_GeometricGoo geoGoo = WKBLineStringSingle( ring, ref pattern );
geoOutput.Append( geoGoo, path );
}
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
cullPattern.Append( gooPattern, path );
feature = layer.GetNextFeature();
counter++;
}
}
return geoOutput;
}
private GH_Structure<IGH_GeometricGoo> WKBGeometryCollection( ref GH_Structure<IGH_Goo> cullPattern )
{
GH_Structure<IGH_GeometricGoo> geoOutput = new GH_Structure<IGH_GeometricGoo>();
cullPattern = new GH_Structure<IGH_Goo>();
GH_Structure<IGH_Goo> cullPatternTemp = new GH_Structure<IGH_Goo>();
if ( numberOfFeatures != -1 )
{
for( int i = 0; i < numberOfFeatures; ++i )
{
path = new GH_Path(i);
feature = layer.GetFeature(i);
if( feature != null )
{
geo = feature.GetGeometryRef();
if( geo != null )
{
bool tempPattern = feature != null ? true : false;
IGH_GeometricGoo geoGoo = null;
if( geo != null )
{
for( int j = 0; j < geo.GetGeometryCount(); ++j )
{
path = new GH_Path( j, i );
ring = geo.GetGeometryRef(i);
OSGeo.OGR.wkbGeometryType geoType = ring.GetGeometryType();
switch( geoType )
{
case OSGeo.OGR.wkbGeometryType.wkbPoint:
geoGoo = WKBPointSingle( ring );
break;
case OSGeo.OGR.wkbGeometryType.wkbPoint25D:
geoGoo = WKBPointSingle( ring );
break;
case OSGeo.OGR.wkbGeometryType.wkbPolygon:
geoGoo = WKBPolygonSingle( ring, ref tempPattern );
break;
case OSGeo.OGR.wkbGeometryType.wkbLineString:
geoGoo = WKBLineStringSingle( ring, ref tempPattern );
break;
}
geoOutput.Append( geoGoo, path );
IGH_Goo gooPattern = GH_Convert.ToGoo( tempPattern );
cullPattern.Append( gooPattern, path );
}
}
}
}
}
}
else
{
int counter = 0;
feature = layer.GetNextFeature();
while( feature != null )
{
geo = feature.GetGeometryRef();
bool tempPattern = feature != null ? true : false;
IGH_GeometricGoo geoGoo = null;
if( geo != null )
{
for( int i = 0; i < geo.GetGeometryCount(); ++i )
{
path = new GH_Path( counter, i );
ring = geo.GetGeometryRef(i);
OSGeo.OGR.wkbGeometryType geoType = ring.GetGeometryType();
switch( geoType )
{
case OSGeo.OGR.wkbGeometryType.wkbPoint:
geoGoo = WKBPointSingle( ring );
break;
case OSGeo.OGR.wkbGeometryType.wkbPoint25D:
geoGoo = WKBPointSingle( ring );
break;
case OSGeo.OGR.wkbGeometryType.wkbPolygon:
geoGoo = WKBPolygonSingle( ring, ref tempPattern );
break;
case OSGeo.OGR.wkbGeometryType.wkbLineString:
geoGoo = WKBLineStringSingle( ring, ref tempPattern );
break;
}
geoOutput.Append( geoGoo, path );
IGH_Goo gooPattern = GH_Convert.ToGoo( tempPattern );
cullPattern.Append( gooPattern, path );
}
}
feature = layer.GetNextFeature();
counter++;
}
}
return geoOutput;
}
private IGH_GeometricGoo WKBPointSingle( OSGeo.OGR.Geometry geo )
{
IGH_GeometricGoo geoGoo = null;
if( currentProjection != EPSGCode )
{
geo.Transform( transform );
}
geo.GetPoint( 0, pointList );
Point3d tempPoint = new Point3d( pointList[0], pointList[1], pointList[2] );
geoGoo = GH_Convert.ToGeometricGoo( tempPoint );
return geoGoo;
}
private IGH_GeometricGoo WKBPolygonSingle( OSGeo.OGR.Geometry geo, ref bool pattern )
{
IGH_GeometricGoo geoGoo = null;
ring = geo.GetGeometryRef( 0 );
if( currentProjection != EPSGCode )
{
ring.Transform( transform );
}
int pointCount = ring.GetPointCount();
Point3d[] tempPointArray = new Point3d[pointCount];
for (int j = 0; j < pointCount; ++j)
{
ring.GetPoint( j, pointList );
tempPointArray[j] = new Point3d( pointList[0], pointList[1], pointList[2] );
}
if( tempPointArray.Length > 1 )
{
Polyline polyOut = new Polyline( tempPointArray );
geoGoo = GH_Convert.ToGeometricGoo( polyOut );
pattern = true;
}
return geoGoo;
}
private IGH_GeometricGoo WKBLineStringSingle( OSGeo.OGR.Geometry geo, ref bool pattern )
{
IGH_GeometricGoo geoGoo = null;
if( currentProjection != EPSGCode )
{
geo.Transform( transform );
}
int pointCount = geo.GetPointCount();
Point3d[] tempPointArray = new Point3d[pointCount];
for( int j = 0; j < pointCount; ++j )
{
geo.GetPoint( j, pointList );
tempPointArray[j] = new Point3d( pointList[0], pointList[1], pointList[2] );
}
if( tempPointArray.Length > 1 )
{
Polyline polyOut = new Polyline( tempPointArray );
geoGoo = GH_Convert.ToGeometricGoo( polyOut );
pattern = true;
}
IGH_Goo gooPattern = GH_Convert.ToGoo( pattern );
return geoGoo;
}
}
你会打电话给它（确保你打电话给
GdalConfiguration.ConfigureOgr();
GdalConfiguration.ConfigureGdal();
):
OSGeo.OGR.Driver driver = Utilities.GetExtractorDriver( input, ref debug );
OSGeo.OGR.DataSource dataSource = driver.Open( input, 0 );
WKBToRhinoGeoUtilities wkbToRhino = new WKBToRhinoGeoUtilities( input, layerNumber, EPSGCode, InEPSGCode, dataSource );
debug = wkbToRhino.Setup();
AddRuntimeMessage( GH_RuntimeMessageLevel.Blank, debug );
if( debug != "Driver is null!" || debug != "Data Source is null!" || debug != "File does not exist!" )
{
GH_Structure<IGH_Goo> cullPattern = new GH_Structure<IGH_Goo>();
GH_Structure<IGH_GeometricGoo> geoOutput = wkbToRhino.Run( ref cullPattern );
DA.SetDataTree( 0, geoOutput );
DA.SetDataTree( 1, cullPattern );
}
我有一个Utilities带有这个静态方法的类：
public static OSGeo.OGR.Driver GetExtractorDriver( string inputString, ref string debug )
{
if( File.Exists( inputString ) )
{
debug = "File exists!";
string isRunning = RunGDAL.output;
OSGeo.OGR.Driver driver = null;
string[] fileTypeArray = inputString.Split( '.' );
string fileType = fileTypeArray[fileTypeArray.Length - 1];
if( fileType == "shp" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "ESRI Shapefile" );
}
else if( fileType == "json" || fileType == "geojson" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "GeoJSON" );
}
else if( fileType == "topojson" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "TopoJSON" );
}
else if( fileType == "gpkg" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "GPKG" );
}
else if( fileType == "dxf" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "DXF" );
}
else if( fileType == "kml" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "KML" );
}
else if( fileType == "osm" )
{
driver = OSGeo.OGR.Ogr.GetDriverByName( "OSM" );
}
return driver;
}
else
{
debug = "File does not exist!";
return null;
}
}