[PROJ] How to write LAS points in N-E CRS

[Javier Jimenez Shaw]

Hi

I am writing a LAS file using a Northing-Easting CRS (EPSG:2384 for
instance).

I am not sure in which order I should write the point coordinates:
northing-easting-up or easting-northing-up.
* The first option (northing-easting-up) is "consistent" with the
definition of the CRS. However it needs to understand the axes of the CRS
(is this only possible in WKT2?), and operate the values to display them on
the screen. Otherwise you will see a reflected image.
* The second option (easting-northing-up) is much easier to represent in
the screen. And seems to be what many programs are doing/expecting.
However, it is inconsistent with the CRS definition.

I don't find anything about that in LAS specs.

What should I do?

Should I write the WKT string in WKT1 or WKT2 format? if WKT1, with or
without axes?

The same question could apply to ShapeFiles.

projinfo EPSG:2384 -o WKT1:GDAL
WKT1:GDAL string:
PROJCS["Xian 1980 / 3-degree Gauss-Kruger CM 117E",
    GEOGCS["Xian 1980",
        DATUM["Xian_1980",
            SPHEROID["IAG 1975",6378140,298.257,
                AUTHORITY["EPSG","7049"]],
            AUTHORITY["EPSG","6610"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4610"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2384"]]


Thanks

.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__
Entre dos pensamientos racionales
hay infinitos pensamientos irracionales.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20200629/9853f2a6/attachment.html>