Question

此时此刻，我正在使用C ++和QT的shapefile遮阳板并使用GDAL / OGR库。我有这种方法来获取我的shapefile的EPSG：

OGRLayer layer = dataset->GetLayer(0);
OGRSpatialReference *spatialRef = layer->GetSpatialRef();

有了这个，我得到了EPSG编号：

atoi(spatialRef->GetAuthorityCode(NULL));

这项工作在我的所有形状文件中都不错。在这种情况下，该方法始终返回null。

我尝试使用：

spatialRef->GetAuthorityCode("PROJCS");
spatialRef->GetAuthorityCode("GEOGCS");
spatialRef->GetAuthorityName("GEOGCS");

所有这些方法都会返回""。

我在gis程序中检查这个shapefile，因为QGIS和QGIS自动检测到他的EPSG是25830。

我的问题是：投影信息是否可以采用与我所做的不同的方法进行处理？

我等你的建议。

非常感谢。

修改

这是.prj文件的内容：

PROJCS [＆＃34; ETRS89_UTM_zone_30N＆＃34;，GEOGCS [＆＃34; GCS_ETRS_1989＆＃34;，DATUM [＆＃34; D_ETRS_1989＆＃34;，SPHEROID [＆＃34; GRS_1980＆＃34;，6378137,298.257222101 ]，PRIMEM [＆＃34;格林威治＆＃34;，0]，单位[＆＃34;度＆＃34;，0.017453292519943295]，投影[＆＃34; Transverse_Mercator＆＃34]，PARAMETER [＆＃34; latitude_of_origin＆＃34;，0]，PARAMETER [＆＃34; central_meridian＆＃34;， - 3]，PARAMETER [＆＃34; scale_factor＆＃34;，0.9996]，PARAMETER [＆＃34; false_easting＆＃34;，500000] ，PARAMETER [＆＃34; false_northing＆＃34;，0]，单位[＆＃34;表＆＃34;，1]]

Answer 1

这样的事情应该有效：

OGRLayer  * layer = dataset->GetLayer(0);
layer->ResetReading();
OGRFeature * feat= layer->GetNextFeature(); 
OGRGeometry * geom = feat->GetGeometryRef(); 
OGRSpatialReference * spatRef = geom->getSpatialReference(); 
int EPSG =  spatRef->GetEPSGGeogCS();

希望它有所帮助！

无法通过OGR / GDAL从Shapefile获取EPSG

1 个答案: