gnudatalanguage / gdl

GDL - GNU Data Language
GNU General Public License v2.0
274 stars 61 forks source link

try to get epsg #1730

Closed sagido closed 7 months ago

sagido commented 7 months ago

auto sr = poDataset->GetSpatialRef(); epsg = sr->GetEPSGGeogCS();

PROJCS["Israeli Grid 05/12",GEOGCS["IG05/12 Intermediate CRS",DATUM["IG05_12_Intermediate_Datum",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","1144"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","6990"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",31.7343936111111],PARAMETER["central_meridian",35.2045169444444],PARAMETER["scale_factor",1.0000067],PARAMETER["false_easting",219529.584],PARAMETER["false_northing",626907.39],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","6991"]]

return epsg=6990 instead of 6991.

How to take the 6991 ?

alaingdl commented 7 months ago

Are you sure this is a GDL issue ? any chance to have more the information ? I cannot play just with those two lines !

GillesDuvert commented 7 months ago

This looks like GDAL commands. Not GDL ! :smile: