如何使用 ogr 库和 c++ 提取 ESRI 形状文件中几何的折点
How to extract vertexes of geometries in ESRI shapefiles using ogr library with c++
我已经回答了这个问题,并将问题的标题从如何使用 ogr 库和 c++ 访问面图层中的折点更改为如何使用 ogr 库和 c++ 提取 ESRI 形状文件中几何的顶点。该代码可用于使用 OpenGL 显示形状。
任何使其更好的建议将不胜感激。
例如,您是否认为我应该使用图层投影来增强OpenGL中的视图?如何使用投影并将其应用于 OpenGL 窗口?
到目前为止,我已经尝试了这些类来提取多边形的顶点:
- OGRPolygon :它没有任何方法来获取顶点数,起点和终点。
- OGRLineString :当我使用
getNumPoints()
方法时,我得到 1图层中的面要素。 - OGRLinearRing :当我使用
getNumPoints()
方法时,我得到 1图层中的面要素。
我应该怎么做才能在多边形形状文件中获取顶点的数量和顶点坐标,以便在 OpenGL 窗口中绘制多边形。
我找到了答案:
您应该关心有关 ESRI 形状文件的一些信息。
-
点图层 :只能具有 wkbpoint 类型的要素。因此,您可以定义一个名为 MyPoint2D 的数据结构并存储具有类型成员的向量中此图层的坐标我的点2D。
多点 图层:只能具有 wkb多点类型的要素。因此,您可以将每个要素存储在具有 MyPoint2D 类型成员的向量中,并且此类向量的集合将保存图层中所有要素的坐标。
折线图层:可以具有 wkbLineString 和 wkbMultiLineString 类型的特征。因此,您可以拥有一个数据结构 MyLine2D,它是一个具有 MyPoint2D 成员的向量,一个具有 MyLine2D 类型的成员的向量将保存线要素,最后是图层中的要素将存储在具有线要素类型成员的矢量中。
多边形图层 :关于多边形,您应该知道每个多边形都由一个外环和几个内部环组成戒指(这是我真正的问题,导致我问问题)。并且该图层也可以具有类型的要素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 形状文件中折点坐标的整个代码:
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/Grasshopper
API的一部分)。但基本上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;
}
}
- .cpp和.h文件中的模板专用化声明
- 为什么两个不同的未命名名称空间可以共存于一个cpp文件中
- 文本文件中的单词链表
- CMake-按正确顺序将项目与C运行时对象文件链接
- 使用新行和不使用新行读取文件
- 在C++程序中输入的文本文件将不起作用,除非文本被复制和粘贴
- 挂起和取消挂起一个文件DLL
- 如何确定我已使用非编码文件到达 EOF?
- 命名空间中具有.h和.cpp文件的类
- 如何使用ndk-build.cmd构建Android.so文件
- 从包含m行的文件中提取n行,必要时(惰性地)重复该文件
- 读取文件并输入到矢量中
- 在C++中查找文件
- c++库的公共头文件中应该包含什么
- 用c++从输入文件中读取另一行
- Cppcheck生成xml转储文件
- 读取文件的最后一行并输入到链接列表时出错
- 无法编译 rtmidi 测试 cmidiin.cpp 文件, 非法指令
- 无法使用 OGR/GDAL 从形状文件中获取 EPSG
- 如何使用 ogr 库和 c++ 提取 ESRI 形状文件中几何的折点