[PROJ] BoundCRS with BBox

[Even Rouault]

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
> 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/3fccbdb1/attachment-0001.html>