无法使用 OGR/GDAL 从形状文件中获取 EPSG

Can't get EPSG from Shapefile with OGR/GDAL

本文关键字：文件 EPSG 获取 GDAL OGR      更新时间：2023-10-16

目前，我正在使用C++和QT的shapefile遮阳板并使用GDAL/OGR库。 我有这种方法来获取我的形状文件的 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 程序中检查这个形状文件，因为 QGIS 和 QGIS 自动检测到他的 EPSG 是 25830。

我的问题是：是否可以使用与我正在执行的操作不同的方法读取投影信息？

我等待你的建议。

多谢。

编辑

这是 .prj 文件的内容：

PROJCS["ETRS89_UTM_zone_30N"，GEOGCS["GCS_ETRS_1989"，DATUM["D_ETRS_1989"，SPHEROID["GRS_1980"，6378137,298.257222101]]，PRIMEM["格林威治"，0]，UNIT["度"，0.017453292519943295]]，投影

["Transverse_Mercator"]，参数["latitude_of_origin"，0]，参数["central_meridian"，-3]，参数["scale_factor"，0.9996]，参数["false_easting"，500000]，参数["false_northing"，0]，单位["米"，1]]

这样的事情应该可以工作：

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

希望对您有所帮助！

相关文章：