无法通过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");

所有这些方法都归还了"".

我在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 [ "格林威治",0],单位[ "度",0.017453292519943295],投影[ "Transverse_Mercator"],PARAMETER [ "latitude_of_origin",0],PARAMETER [ "central_meridian", - 3],PARAMETER [ "scale_factor",0.9996],PARAMETER [ "false_easting",500000],PARAMETER [ "false_northing",0 ],单位[ "表",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(); 

希望能帮助到你！