[PROJ] PROJ 6.1.0RC1

Bas Couwenberg sebastic at xs4all.nl
Mon May 6 01:29:53 PDT 2019


On 2019-05-06 10:08, Kristian Evers wrote:
> From the previous thread it seems to be a difficult problem without a 
> clear
> Solution. Looking at the test output though it seems as if the vertical 
> shift
> isn't  applied. Whether that is because of the issues discussed 
> previously or
> it is a problem in finding the naptrans2008.gtx grid on the system I 
> can't
> tell from the log file.
> 
> The log file is also reporting that the epsg init file fails to load:
> 
> 	pj_open_lib(epsg): call fopen(/usr/share/proj/epsg) - failed

This is expected, because the epsg init file no longer exist for PROJ 6.

The minimal epsg init file is also moved out of the way to prevent it 
from being used, this file was used with prior versions that only 
supporting a single path in PROJ_LIB.

> So, is this simply an incorrect setting of PROJ_LIB or is it still the 
> same
> not so simple issue as reported previously?

It's probably the not so simple issue.

IIRC epsg:4258 is interpreted as a Geographic 2D CRS which likely causes 
the 3D component to be ignored.

The projinfo output is not something I can make sense of, it is as 
follows:

$ PROJ_LIB=.:/usr/share/proj projinfo -s EPSG:4258 -t rdnap:rdnap
target CRS: parsing of user string failed: crs not found

$ PROJ_LIB=.:/usr/share/proj projinfo -s EPSG:4258 -t "+proj=sterea 
+lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 
+x_0=155000 +y_0=463000 +ellps=bessel +nadgrids=rdtrans2008.gsb 
+geoidgrids=naptrans2008.gtx +units=m +no_defs +type=crs"
Candidate operations found: 1
-------------------------------------
Operation n°1:

unknown id, unknown + unknown to WGS84 + Inverse of ETRS89 to WGS 84 (1) 
+ Transformation from unknown to ETRS89 (ballpark vertical 
transformation, without ellipsoid height to vertical height correction), 
unknown accuracy, Europe - ETRS89, has ballpark transformation

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg +xy_out=rad +step +inv +proj=hgridshift 
+grids=rdtrans2008.gsb +step +proj=sterea +lat_0=52.1561605555556 
+lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 
+ellps=bessel

WKT2_2018 string:
COORDINATEOPERATION["unknown + unknown to WGS84 + Inverse of ETRS89 to 
WGS 84 (1) + Transformation from unknown to ETRS89 (ballpark vertical 
transformation, without ellipsoid height to vertical height 
correction)",
     SOURCECRS[
         GEOGCRS["ETRS89",
             DATUM["European Terrestrial Reference System 1989",
                 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",4258]]],
     TARGETCRS[
         COMPOUNDCRS["unknown",
             BOUNDCRS[
                 SOURCECRS[
                     PROJCRS["unknown",
                         BASEGEOGCRS["unknown",
                             DATUM["Unknown based on Bessel 1841 
ellipsoid",
                                 ELLIPSOID["Bessel 
1841",6377397.155,299.1528128,
                                     LENGTHUNIT["metre",1,
                                         ID["EPSG",9001]]]],
                             PRIMEM["Greenwich",0,
                                 ANGLEUNIT["degree",0.0174532925199433],
                                 ID["EPSG",8901]]],
                         CONVERSION["unknown",
                             METHOD["Oblique Stereographic",
                                 ID["EPSG",9809]],
                             PARAMETER["Latitude of natural 
origin",52.1561605555556,
                                 ANGLEUNIT["degree",0.0174532925199433],
                                 ID["EPSG",8801]],
                             PARAMETER["Longitude of natural 
origin",5.38763888888889,
                                 ANGLEUNIT["degree",0.0174532925199433],
                                 ID["EPSG",8802]],
                             PARAMETER["Scale factor at natural 
origin",0.9999079,
                                 SCALEUNIT["unity",1],
                                 ID["EPSG",8805]],
                             PARAMETER["False easting",155000,
                                 LENGTHUNIT["metre",1],
                                 ID["EPSG",8806]],
                             PARAMETER["False northing",463000,
                                 LENGTHUNIT["metre",1],
                                 ID["EPSG",8807]]],
                         CS[Cartesian,2],
                             AXIS["(E)",east,
                                 ORDER[1],
                                 LENGTHUNIT["metre",1,
                                     ID["EPSG",9001]]],
                             AXIS["(N)",north,
                                 ORDER[2],
                                 LENGTHUNIT["metre",1,
                                     ID["EPSG",9001]]]]],
                 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["unknown to WGS84",
                     METHOD["NTv2",
                         ID["EPSG",9615]],
                     PARAMETERFILE["Latitude and longitude difference 
file","rdtrans2008.gsb",
                         ID["EPSG",8656]]]],
             BOUNDCRS[
                 SOURCECRS[
                     VERTCRS["unknown",
                         VDATUM["unknown"],
                         CS[vertical,1],
                             AXIS["gravity-related height (H)",up,
                                 LENGTHUNIT["metre",1,
                                     ID["EPSG",9001]]]]],
                 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,3],
                             AXIS["latitude",north,
                                 ORDER[1],
                                 ANGLEUNIT["degree",0.0174532925199433]],
                             AXIS["longitude",east,
                                 ORDER[2],
                                 ANGLEUNIT["degree",0.0174532925199433]],
                             AXIS["ellipsoidal height",up,
                                 ORDER[3],
                                 LENGTHUNIT["metre",1]],
                         ID["EPSG",4979]]],
                 ABRIDGEDTRANSFORMATION["unknown to WGS84 ellipsoidal 
height",
                     METHOD["GravityRelatedHeight to Geographic3D"],
                     PARAMETERFILE["Geoid (height correction) model 
file","naptrans2008.gtx",
                         ID["EPSG",8666]]]]]],
     METHOD["PROJ-based operation method (approximate): +proj=pipeline 
+step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg 
+xy_out=rad +step +inv +proj=hgridshift +grids=rdtrans2008.gsb +step 
+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel"],
     USAGE[
         SCOPE["unknown"],
         AREA["Europe - ETRS89"],
         BBOX[32.88,-16.1,84.17,40.18]]]

Kind Regards,

Bas


More information about the PROJ mailing list