WKT2 with force grid using

[grigoroffbenj]

Hi,
i have a work zone near Calgary that i need to transform NAD83 (old)/ Alberta 3-degree TM reference meridian 114 W (EPSG:3776) with CGVD28 height with GSD95 geoid model TO EPSG:4326 or EPSG:4979 (hellipsoid height).

Before update (Eg. in version 3.3.1) I can't deal with geoid convertion but i can make NTV2 plannimetrics shift with "ca_nrc_ABCSRSV4.tif" file in PROJ_DATA.

I have updated pyproj to 3.4.0 and try to deal with plannimetric and vertical convertion.

I have created a WKT2 of my projection:

COMPOUNDCRS["NAD83 / Alberta 3TM ref merid 114 W + GSD95 height",
	BOUNDCRS[
		SOURCECRS[
			PROJCRS["NAD83 / Alberta 3TM ref merid 114 W",
				BASEGEOGCRS["NAD83",
					DATUM["North American Datum 1983",
						ELLIPSOID["GRS 1980",6378137,298.257222101,
							LENGTHUNIT["metre",1]]],
					PRIMEM["Greenwich",0,
						ANGLEUNIT["degree",0.0174532925199433]],
					ID["EPSG",4269]],
				CONVERSION["Alberta 3-degree TM reference meridian 114 W",
					METHOD["Transverse Mercator",
						ID["EPSG",9807]],
					PARAMETER["Latitude of natural origin",0,
						ANGLEUNIT["degree",0.0174532925199433],
						ID["EPSG",8801]],
					PARAMETER["Longitude of natural origin",-114,
						ANGLEUNIT["degree",0.0174532925199433],
						ID["EPSG",8802]],
					PARAMETER["Scale factor at natural origin",0.9999,
						SCALEUNIT["unity",1],
						ID["EPSG",8805]],
					PARAMETER["False easting",0,
						LENGTHUNIT["metre",1],
						ID["EPSG",8806]],
					PARAMETER["False northing",0,
						LENGTHUNIT["metre",1],
						ID["EPSG",8807]]],
				CS[Cartesian,2],
					AXIS["(E)",east,
						ORDER[1],
						LENGTHUNIT["metre",1]],
					AXIS["(N)",north,
						ORDER[2],
						LENGTHUNIT["metre",1]],
				USAGE[
					SCOPE["Cadastre, survey control and engineering survey in urban areas, typically in all municipalities that previously comprised the 73 Municipal Integrated Surveying and Mapping (MISAM) areas, also known as urban cadastral map areas. For rural areas use UTM."],
					AREA["Canada - Alberta - between 115°30W and 112°30W."],
					BBOX[48.99,-115.5,60,-112.5]],
				ID["EPSG",3776]]],
		TARGETCRS[
			GEOGCRS["WGS 84",
				ENSEMBLE["World Geodetic System 1984 ensemble",
					MEMBER["World Geodetic System 1984 (Transit)"],
					MEMBER["World Geodetic System 1984 (G730)"],
					MEMBER["World Geodetic System 1984 (G873)"],
					MEMBER["World Geodetic System 1984 (G1150)"],
					MEMBER["World Geodetic System 1984 (G1674)"],
					MEMBER["World Geodetic System 1984 (G1762)"],
					MEMBER["World Geodetic System 1984 (G2139)"],
					ELLIPSOID["WGS 84",6378137,298.257223563,
						LENGTHUNIT["metre",1]],
					ENSEMBLEACCURACY[2.0]],
				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]],
				USAGE[
					SCOPE["Horizontal component of 3D system."],
					AREA["World."],
					BBOX[-90,-180,90,180]],
				ID["EPSG",4326]]],
		ABRIDGEDTRANSFORMATION["NAD83 to WGS 84 (8)",
			METHOD["NTv2",
				ID["EPSG",9615]],
			PARAMETERFILE["Latitude and longitude difference file","ca_nrc_ABCSRSV4.tif"],
			USAGE[
				SCOPE["(null/copy) Approximation for medium and low accuracy applications assuming equality between plate-fixed static and earth-fixed dynamic CRSs, ignoring static/dynamic CRS differences."],
				AREA["Canada - Alberta."],
				BBOX[48.99,-120,60,-109.98]],
			ID["EPSG",1702],
			REMARK["Parameter file is from NAD83 to NAD83(CSRS)v4 (3) (code 9244) assuming that NAD83(CSRS)v4 is equivalent to WGS 84 within the accuracy of the transformation. This file AB_CSRS.DAC will need to be renamed to AB_CSRS.gsb to run in some software."]]
	],
	VERTCRS["GSD95",
		VDATUM["NAD83 to NAD83 (EPSG::4269) to GSD95"],
		CS[vertical,1],
		AXIS["gravity-related height (H)",up,LENGTHUNIT["metre",1,ID["EPSG",9001]]],
		ABRIDGEDTRANSFORMATION["GSD95 to WGS 84 ellipsoidal height",
			METHOD["GravityRelatedHeight to Geographic3D"],
			PARAMETERFILE["NAD83 to NAD83 (EPSG::4269) to GSD95","NGSD95.tif"]
		]
	],
	USAGE[
		SCOPE["Cadastre, survey control and engineering survey in urban areas, typically in all municipalities that previously comprised the 73 Municipal Integrated Surveying and Mapping (MISAM) areas, also known as urban cadastral map areas. For rural areas use UTM."],
		AREA["Canada - Alberta - between 115°30W and 112°30W."],
		BBOX[48.99,-115.5,60,-112.5]]
]

This is my python code snippet.

# get CRS object
crsIn = crsIn = pyproj.CRS.from_wkt(<myWKT2String>).to_3d()
crsOut = pyproj.CRS("EPSG:4979").to_3d()
# create a transformGroup for this conversion
transform = pyproj.transformer.TransformerGroup(crsIn, crsOut, always_xy=True, allow_ballpark=True)

Problem description

The TransformerGroup have ballpack with my GSD95 grid file (in fact it's not in PROJ_DATA V1.11) so,

• Downloaded it from CA_NRCAN web site
  https://webapp.csrs-scrs.nrcan-rncan.gc.ca/geod/process/download-helper.php?file_id=GSD95
• I putted it in my PROJ_DATA folder.

But i still have ballpack.

{
  "0": {
    "accuracy": -1,
    "definition": "proj=pipeline step inv proj=tmerc lat_0=0 lon_0=-114 k=0.9999 x_0=0 y_0=0 ellps=GRS80 step proj=hgridshift grids=ca_nrc_ABCSRSV4.tif step proj=unitconvert xy_in=rad xy_out=deg",
    "description": "Inverse of Alberta 3-degree TM reference meridian 114 W + NAD83 to WGS 84 (8) + Transformation from GSD95 to WGS 84 (ballpark vertical transformation, without ellipsoid height to vertical height correction) + axis order change (2D)",
    "proj4": "+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=-114 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +step +proj=hgridshift +grids=ca_nrc_ABCSRSV4.tif +step +proj=unitconvert +xy_in=rad +xy_out=deg"
  }
}

What I'm doing wrong ? WKT2 definition ? or ?
I see in this post #1114 (comment) a different WKT2 definition.

Environment (from conda list & pyproj.datadir.get_user_data_dir()):

affine                    2.3.1              pyhd8ed1ab_0    conda-forge
aiofiles                  22.1.0             pyhd8ed1ab_0    conda-forge
aiohttp                   3.8.3            py38h91455d4_0    conda-forge
aiosignal                 1.2.0              pyhd8ed1ab_0    conda-forge
anyio                     3.6.2              pyhd8ed1ab_0    conda-forge
argon2-cffi               21.3.0             pyhd8ed1ab_0    conda-forge
argon2-cffi-bindings      21.2.0           py38h294d835_2    conda-forge
async-timeout             4.0.2              pyhd8ed1ab_0    conda-forge
attrs                     22.1.0             pyh71513ae_1    conda-forge
bcrypt                    3.2.2            py38h294d835_0    conda-forge
blosc                     1.21.1               h74325e0_3    conda-forge
boost-cpp                 1.78.0               h9f4b32c_1    conda-forge
branca                    0.5.0              pyhd8ed1ab_0    conda-forge
brotli                    1.0.9                h8ffe710_7    conda-forge
brotli-bin                1.0.9                h8ffe710_7    conda-forge
brotlipy                  0.7.0           py38h294d835_1004    conda-forge
bzip2                     1.0.8                h8ffe710_4    conda-forge
ca-certificates           2022.9.24            h5b45459_0    conda-forge
cairo                     1.16.0            hd694305_1014    conda-forge
certifi                   2022.9.24          pyhd8ed1ab_0    conda-forge
cffi                      1.15.1           py38h57701bc_1    conda-forge
cfitsio                   4.1.0                h5a969a9_0    conda-forge
charset-normalizer        2.1.1              pyhd8ed1ab_0    conda-forge
click                     8.1.3            py38haa244fe_0    conda-forge
click-plugins             1.1.1                      py_0    conda-forge
cligj                     0.7.2              pyhd8ed1ab_1    conda-forge
colorama                  0.4.5              pyhd8ed1ab_0    conda-forge
contourpy                 1.0.5            py38hb1fd069_0    conda-forge
cryptography              38.0.2           py38h95f5157_1    conda-forge
curl                      7.85.0               heaf79c2_0    conda-forge
cycler                    0.11.0             pyhd8ed1ab_0    conda-forge
ecdsa                     0.18.0             pyhd8ed1ab_1    conda-forge
exifread                  3.0.0              pyhd8ed1ab_0    conda-forge
expat                     2.4.9                h1537add_0    conda-forge
fastapi                   0.85.1             pyhd8ed1ab_0    conda-forge
fiona                     1.8.22           py38h6f5d430_0    conda-forge
flit-core                 3.7.1              pyhd8ed1ab_0    conda-forge
folium                    0.13.0             pyhd8ed1ab_0    conda-forge
font-ttf-dejavu-sans-mono 2.37                 hab24e00_0    conda-forge
font-ttf-inconsolata      3.000                h77eed37_0    conda-forge
font-ttf-source-code-pro  2.038                h77eed37_0    conda-forge
font-ttf-ubuntu           0.83                 hab24e00_0    conda-forge
fontconfig                2.14.1               hbde0cde_0    conda-forge
fonts-conda-ecosystem     1                             0    conda-forge
fonts-conda-forge         1                             0    conda-forge
fonttools                 4.38.0           py38h91455d4_0    conda-forge
freetype                  2.12.1               h546665d_0    conda-forge
freexl                    1.0.6                h67ca5e6_1    conda-forge
frozenlist                1.3.1            py38h294d835_0    conda-forge
gdal                      3.5.2            py38h5a6f081_5    conda-forge
geopandas                 0.12.0             pyhd8ed1ab_0    conda-forge
geopandas-base            0.12.0             pyha770c72_0    conda-forge
geos                      3.11.0               h39d44d4_0    conda-forge
geotiff                   1.7.1                h4ffd875_4    conda-forge
gettext                   0.21.1               h5728263_0    conda-forge
gmpy2                     2.1.2            py38hb63f22f_0    conda-forge
h11                       0.14.0             pyhd8ed1ab_0    conda-forge
hdf4                      4.2.15               h0e5069d_4    conda-forge
hdf5                      1.12.2          nompi_h57737ce_100    conda-forge
icu                       70.1                 h0e60522_0    conda-forge
idna                      3.4                pyhd8ed1ab_0    conda-forge
intel-openmp              2022.1.0          h57928b3_3787    conda-forge
jinja2                    3.1.2              pyhd8ed1ab_1    conda-forge
joblib                    1.2.0              pyhd8ed1ab_0    conda-forge
jpeg                      9e                   h8ffe710_2    conda-forge
kealib                    1.4.15               hdf81f3a_1    conda-forge
kiwisolver                1.4.4            py38hbd9d945_0    conda-forge
krb5                      1.19.3               hc8ab02b_0    conda-forge
lcms2                     2.12                 h2a16943_0    conda-forge
lerc                      4.0.0                h63175ca_0    conda-forge
libblas                   3.9.0              16_win64_mkl    conda-forge
libbrotlicommon           1.0.9                h8ffe710_7    conda-forge
libbrotlidec              1.0.9                h8ffe710_7    conda-forge
libbrotlienc              1.0.9                h8ffe710_7    conda-forge
libcblas                  3.9.0              16_win64_mkl    conda-forge
libcurl                   7.85.0               heaf79c2_0    conda-forge
libdeflate                1.14                 hcfcfb64_0    conda-forge
libffi                    3.4.2                h8ffe710_5    conda-forge
libgdal                   3.5.2                h283e77c_5    conda-forge
libglib                   2.74.0               h79619a9_0    conda-forge
libiconv                  1.17                 h8ffe710_0    conda-forge
libkml                    1.3.0             hf2ab4e4_1015    conda-forge
liblapack                 3.9.0              16_win64_mkl    conda-forge
libnetcdf                 4.8.1           nompi_h85765be_104    conda-forge
libpng                    1.6.38               h19919ed_0    conda-forge
libpq                     14.5                 h1ea2d34_0    conda-forge
librttopo                 1.1.0               h2842628_11    conda-forge
libspatialindex           1.9.3                h39d44d4_4    conda-forge
libspatialite             5.0.1               hd560d62_21    conda-forge
libsqlite                 3.39.4               hcfcfb64_0    conda-forge
libssh2                   1.10.0               h9a1e1f7_3    conda-forge
libtiff                   4.4.0                h8e97e67_4    conda-forge
libwebp-base              1.2.4                h8ffe710_0    conda-forge
libxcb                    1.13              hcd874cb_1004    conda-forge
libxml2                   2.10.3               hc3477c8_0    conda-forge
libzip                    1.9.2                h519de47_1    conda-forge
libzlib                   1.2.13               hcfcfb64_4    conda-forge
lz4-c                     1.9.3                h8ffe710_1    conda-forge
m2w64-gcc-libgfortran     5.3.0                         6    conda-forge
m2w64-gcc-libs            5.3.0                         7    conda-forge
m2w64-gcc-libs-core       5.3.0                         7    conda-forge
m2w64-gmp                 6.1.0                         2    conda-forge
m2w64-libwinpthread-git   5.0.0.4634.697f757               2    conda-forge
mapclassify               2.4.3              pyhd8ed1ab_0    conda-forge
markupsafe                2.1.1            py38h294d835_1    conda-forge
matplotlib-base           3.6.1            py38h528a6c7_0    conda-forge
mkl                       2022.1.0           h6a75c08_874    conda-forge
mpc                       1.2.1                h54e1faf_0    conda-forge
mpfr                      4.1.0                h8d14728_1    conda-forge
mpir                      3.0.0             he025d50_1002    conda-forge
msys2-conda-epoch         20160418                      1    conda-forge
multidict                 6.0.2            py38h294d835_1    conda-forge
munch                     2.5.0                      py_0    conda-forge
munkres                   1.1.4              pyh9f0ad1d_0    conda-forge
nest-asyncio              1.5.6              pyhd8ed1ab_0    conda-forge
networkx                  2.8.7              pyhd8ed1ab_0    conda-forge
numpy                     1.23.4           py38h90ce339_0    conda-forge
openjpeg                  2.5.0                hc9384bd_1    conda-forge
openssl                   3.0.5                hcfcfb64_2    conda-forge
packaging                 21.3               pyhd8ed1ab_0    conda-forge
pandas                    1.5.1            py38h5846ac1_0    conda-forge
passlib                   1.7.4              pyh9f0ad1d_0    conda-forge
pcre                      8.45                 h0e60522_0    conda-forge
pcre2                     10.37                hdfff0fc_1    conda-forge
pillow                    9.2.0            py38h37aa274_2    conda-forge
pip                       22.3               pyhd8ed1ab_0    conda-forge
pixman                    0.40.0               h8ffe710_0    conda-forge
poppler                   22.10.0              ha6c1112_0    conda-forge
poppler-data              0.4.11               hd8ed1ab_0    conda-forge
postgresql                14.5                 he353ca9_0    conda-forge
proj                      9.1.0                h3863b3b_0    conda-forge
proj-data                 1.11                 hd8ed1ab_0    conda-forge
pthread-stubs             0.4               hcd874cb_1001    conda-forge
pyasn1                    0.4.8                      py_0    conda-forge
pycparser                 2.21               pyhd8ed1ab_0    conda-forge
pydantic                  1.10.2           py38h91455d4_0    conda-forge
pyopenssl                 22.1.0             pyhd8ed1ab_0    conda-forge
pyparsing                 3.0.9              pyhd8ed1ab_0    conda-forge
pyproj                    3.4.0            py38hf350bb6_1    conda-forge
pysocks                   1.7.1              pyh0701188_6    conda-forge
python                    3.8.13          hcf16a7b_0_cpython    conda-forge
python-dateutil           2.8.2              pyhd8ed1ab_0    conda-forge
python-jose               3.3.0              pyh6c4a22f_1    conda-forge
python-multipart          0.0.5                      py_0    conda-forge
python_abi                3.8                      2_cp38    conda-forge
pytz                      2022.5             pyhd8ed1ab_0    conda-forge
rasterio                  1.3.3            py38h01d9f53_1    conda-forge
requests                  2.28.1             pyhd8ed1ab_1    conda-forge
rsa                       4.9                pyhd8ed1ab_0    conda-forge
rtree                     1.0.1            py38h8b54edf_0    conda-forge
scikit-learn              1.1.2            py38hc27f28a_0    conda-forge
scipy                     1.9.3            py38h0f6ee2a_0    conda-forge
setuptools                65.5.0             pyhd8ed1ab_0    conda-forge
shapely                   1.8.5            py38h7869501_0    conda-forge
six                       1.16.0             pyh6c4a22f_0    conda-forge
snappy                    1.1.9                h82413e6_1    conda-forge
sniffio                   1.3.0              pyhd8ed1ab_0    conda-forge
snuggs                    1.4.7                      py_0    conda-forge
spark-parser              1.8.9                    pypi_0    pypi
sqlite                    3.39.4               hcfcfb64_0    conda-forge
starlette                 0.20.4             pyhd8ed1ab_1    conda-forge
tbb                       2021.6.0             h91493d7_0    conda-forge
threadpoolctl             3.1.0              pyh8a188c0_0    conda-forge
tiledb                    2.11.3               h3132609_1    conda-forge
tk                        8.6.12               h8ffe710_0    conda-forge
typing-extensions         4.4.0                hd8ed1ab_0    conda-forge
typing_extensions         4.4.0              pyha770c72_0    conda-forge
ucrt                      10.0.20348.0         h57928b3_0    conda-forge
uncompyle6                3.8.0                    pypi_0    pypi
unicodedata2              14.0.0           py38h294d835_1    conda-forge
urllib3                   1.26.11            pyhd8ed1ab_0    conda-forge
uvicorn                   0.19.0           py38haa244fe_0    conda-forge
vc                        14.2                 hac3ee0b_8    conda-forge
vs2015_runtime            14.29.30139          h890b9b1_8    conda-forge
wheel                     0.37.1             pyhd8ed1ab_0    conda-forge
win_inet_pton             1.1.0            py38haa244fe_4    conda-forge
xdis                      6.0.4                    pypi_0    pypi
xerces-c                  3.2.4                h63175ca_1    conda-forge
xorg-libxau               1.0.9                hcd874cb_0    conda-forge
xorg-libxdmcp             1.1.3                hcd874cb_0    conda-forge
xyzservices               2022.9.0           pyhd8ed1ab_0    conda-forge
xz                        5.2.6                h8d14728_0    conda-forge
yarl                      1.7.2            py38h294d835_2    conda-forge
zlib                      1.2.13               hcfcfb64_4    conda-forge
zstd                      1.5.2                h7755175_4    conda-forge

and PROJ_DATA

>>>import pyproj.datadir
>>>pyproj.datadir.get_user_data_dir()
'C:\\Users\\benjamin\\miniconda3\\envs\\API\\share\\proj'

[snowman2]

I am curious whether you have PROJ_NETWORK enabled or disabled when running this?

[grigoroffbenj]

I have tried with this too but as i said, geoid model GSD95 os not in proj_data.
If i try with this option, and an empty proj-data folder the plannimetric grid ca_nrc_ABCSRSV4.tif id downloaded and used but not geoid.