[gdal-dev] Convert GEOCS

[jratike80]

Hi,

There is a command, gdaltransform
https://gdal.org/programs/gdaltransform.html. However, it seems that both
coordinate systems mean the same for Proj library and no conversion in
coordinates happen. That you can test with projinfo that is a Proj utility.

projinfo -s epsg:4326 -t epsg:4674
Candidate operations found: 1
-------------------------------------
Operation No. 1:

INVERSE(EPSG):15894, Inverse of SIRGAS 2000 to WGS 84 (1), 1.0 m, Latin
America - SIRGAS 2000 by country

PROJ string:
+proj=noop

WKT2:2019 string:
COORDINATEOPERATION["Inverse of SIRGAS 2000 to WGS 84 (1)",
    SOURCECRS[
        GEOGCRS["WGS 84",
            DATUM["World Geodetic System 1984",
                ELLIPSOID["WGS 84",6378137,298.257223563,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",4326]]],
    TARGETCRS[
        GEOGCRS["SIRGAS 2000",
            DATUM["Sistema de Referencia Geocentrico para las AmericaS
2000",
                ELLIPSOID["GRS 1980",6378137,298.257222101,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",4674]]],
    METHOD["Geocentric translations (geog2D domain)",
        ID["EPSG",9603]],
    PARAMETER["X-axis translation",0,
        LENGTHUNIT["metre",1],
        ID["EPSG",8605]],
    PARAMETER["Y-axis translation",0,
        LENGTHUNIT["metre",1],
        ID["EPSG",8606]],
    PARAMETER["Z-axis translation",0,
        LENGTHUNIT["metre",1],
        ID["EPSG",8607]],
    OPERATIONACCURACY[1.0],
    USAGE[
        SCOPE["Accuracy 1m."],
        AREA["Latin America - SIRGAS 2000 by country"],
        BBOX[-59.87,-122.19,32.72,-25.28]],
    ID["INVERSE(EPSG)",15894]]


-Jukka Rahkonen-



Moises Silva wrote
> Hello everyone,
> 
>     Is there any command in GDAL that is possible to convert the GEOCS?
>        Source: GCS_SIRGAS_2000
>        Target:  GCS_WGS_1984?
> 
> Thanks a lot,
> Moises
> 
> _______________________________________________
> gdal-dev mailing list

> gdal-dev at .osgeo

> https://lists.osgeo.org/mailman/listinfo/gdal-dev





--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html