[PROJ] BoundCRS with BBox

[Olivier PENA]

Done : https://github.com/OSGeo/PROJ/issues/2813

Thanks


________________________________
De : Even Rouault <even.rouault at spatialys.com>
Envoyé : mercredi 18 août 2021 13:59
À : Olivier PENA <OPena at isagri.fr>; proj at lists.osgeo.org <proj at lists.osgeo.org>
Objet : Re: [PROJ] BoundCRS with BBox


Olivier,

hum, good finding. This is actually an oversight of the implementation. The WKT2 grammar (http://docs.opengeospatial.org/is/18-010r7/18-010r7.html#149) does indeed allow usage/scope/extent/id/remark for a bound CRS but our C++ classes do not allow for that now. Could you file a ticket in https://github.com/OSGeo/PROJ/issues/new about that ?

Even

Le 18/08/2021 à 13:38, Olivier PENA a écrit :

Hi

I’m looking for a way to set area of use for a BoundCRS I define with WKT like this one :



BOUNDCRS[

    SOURCECRS[

        PROJCRS["Monte Mario / Italy zone 2",

            BASEGEOGCRS["Monte Mario",

                DATUM["Monte Mario",

                    ELLIPSOID["International 1924",6378388,297,

                        LENGTHUNIT["metre",1]]],

                PRIMEM["Greenwich",0,

                    ANGLEUNIT["degree",0.0174532925199433]],

                ID["EPSG",4265]],

            CONVERSION["unnamed",

                METHOD["Transverse Mercator",

                    ID["EPSG",9807]],

                PARAMETER["Latitude of natural origin",0,

                    ANGLEUNIT["degree",0.0174532925199433],

                    ID["EPSG",8801]],

                PARAMETER["Longitude of natural origin",15,

                    ANGLEUNIT["degree",0.0174532925199433],

                    ID["EPSG",8802]],

                PARAMETER["Scale factor at natural origin",0.9996,

                    SCALEUNIT["unity",1],

                    ID["EPSG",8805]],

                PARAMETER["False easting",2520000,

                    LENGTHUNIT["metre",1],

                    ID["EPSG",8806]],

                PARAMETER["False northing",0,

                    LENGTHUNIT["metre",1],

                    ID["EPSG",8807]]],

            CS[Cartesian,2],

                AXIS["x",east,

                    ORDER[1],

                    LENGTHUNIT["metre",1]],

                AXIS["y",north,

                    ORDER[2],

                    LENGTHUNIT["metre",1]],

            ID["EPSG",3004]]],

    TARGETCRS[

        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["latitude",north,

                    ORDER[1],

                    ANGLEUNIT["degree",0.0174532925199433]],

                AXIS["longitude",east,

                    ORDER[2],

                    ANGLEUNIT["degree",0.0174532925199433]],

            ID["EPSG",4326]]],

    ABRIDGEDTRANSFORMATION["Transformation from Monte Mario to WGS84",

        METHOD["Position Vector transformation (geog2D domain)",

            ID["EPSG",9606]],

        PARAMETER["X-axis translation",-50.2,

            ID["EPSG",8605]],

        PARAMETER["Y-axis translation",-50.4,

            ID["EPSG",8606]],

        PARAMETER["Z-axis translation",84.8,

            ID["EPSG",8607]],

        PARAMETER["X-axis rotation",-0.69,

            ID["EPSG",8608]],

        PARAMETER["Y-axis rotation",-2.012,

            ID["EPSG",8609]],

        PARAMETER["Z-axis rotation",0.459,

            ID["EPSG",8610]],

        PARAMETER["Scale difference",0.99997192,

            ID["EPSG",8611]]],

    USAGE[

        SCOPE["unknown"],

        AREA["Italy - Sicily onshore"],

        BBOX[36.59,12.36,38.35,15.71]]

]



I would be happy if proj_get_area_of_use() on the crs return the Sicily BBox.



I ran a debug session and it seems that USAGE node is not taken account in buildBoundCRS() in io.cpp. Same behavior in the json parser.



Is this normal behavior ?

Thanks





_______________________________________________
PROJ mailing list
PROJ at lists.osgeo.org<mailto:PROJ at lists.osgeo.org>
https://lists.osgeo.org/mailman/listinfo/proj


--
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20210818/49ddeda0/attachment.html>