[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