如何使用帶有c ++的ogr庫在ESRI shapefile中提取幾何頂點
[英]How to extract vertexes of geometries in ESRI shapefiles using ogr library with c++
問題描述
我已經回答了這個問題，並將問題的標題從如何使用 ogr 庫和 c++ 訪問多邊形層中的頂點更改為如何使用 ogr 庫和 c++提取 ESRI shapefiles 中的幾何頂點。 該代碼可用於使用 OpenGL 顯示形狀。
任何改進它的建議將不勝感激。
例如，您認為我應該使用圖層的投影來增強 OpenGL 中的視圖嗎？ 如何使用投影並將其應用於 OpenGL 窗口？
到目前為止，我已經嘗試過這些類來提取多邊形的頂點：
- OGRPolygon ：它沒有任何方法來獲取頂點數、起點和終點。
- OGRLineString ：當我使用
getNumPoints()方法時，我的圖層中的所有多邊形要素都為 1。 - OGRLinearRing ：當我使用
getNumPoints()方法時，我的圖層中的所有多邊形特征都為 1。
為了在OpenGL窗口中繪制多邊形，我應該怎么做才能獲得多邊形shapefile中的頂點數量和頂點坐標。
2 個解決方案
解決方案1
10 已采納 2013-09-14 21:13:14
我找到了答案：
您應該關心有關 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 形狀文件中頂點的坐標：
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);
}
解決方案2
1 2020-03-22 00:42:46
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 的維度不wkbFlatten它可以工作，但是我將把它留給讀者，因為我的實現不需要它。 請注意， 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;
}
}
使用 GDAL/OGR C++ 將 GeoTiff 的 XY 點重新投影到緯度/經度
[英]Reproject XY points of GeoTiff to Lat/Long using GDAL/OGR C++
有關如何使用簡單的C ++庫（不使用.net）和使用該庫的C ++項目創建.sln的基本說明？
[英]basic instructions on how to create a .sln with simple C++ library (no .net use) and C++ project using that library?
聲明:本站的技術帖子網頁，遵循CC BY-SA 4.0協議，如果您需要轉載，請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.