无法使用 OGR/GDAL 从 Shapefile 获取 EPSG

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");

所有这些方法return "".

我在 gis 程序中检查了这个 shapefile，因为 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["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",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

Can't get EPSG from Shapefile with OGR/GDAL

c++

qt

gdal

ogr