[gdal-dev] Forcing -vshift in Gdal Python API

[David Baker]

Michal,
This may be becuase of an issue with the conda-forge install of pyproj.  The grid shift files do not get installed.
Please see this thread for details:
NAD27 to NAD83 Conversions · pyproj4 pyproj · Discussion #1409 
  
|  
|   
|   
|   |    |

   |

  |
|  
|    |  
NAD27 to NAD83 Conversions · pyproj4 pyproj · Discussion #1409
 
In previous versions of pyproj, I was able to convert points between different datums. Specifically, NAD27 to NAD83 with both geographic and Cartesian coordinates. Is there a setting that needs to ...
  |   |

  |

  |

  


David


Sent from Yahoo Mail


On pyproj4 pyproj6, 2024, 6:58 AM, Michal Girmala via gdal-dev <gdal-dev at lists.osgeo.org> wrote:

 <!--#yiv3625194537 filtered {}#yiv3625194537 filtered {}#yiv3625194537 filtered {}#yiv3625194537 p.yiv3625194537MsoNormal, #yiv3625194537 li.yiv3625194537MsoNormal, #yiv3625194537 div.yiv3625194537MsoNormal {margin:0cm;font-size:11.0pt;font-family:"Calibri", sans-serif;}#yiv3625194537 a:link, #yiv3625194537 span.yiv3625194537MsoHyperlink {color:#0563C1;text-decoration:underline;}#yiv3625194537 span.yiv3625194537EmailStyle19 {font-family:"Calibri", sans-serif;color:windowtext;}#yiv3625194537 .yiv3625194537MsoChpDefault {font-size:10.0pt;}#yiv3625194537 filtered {}#yiv3625194537 div.yiv3625194537WordSection1 {}-->
Hello all,
 
  
 
I am trying to reproject elevation raster with gdalwarp in python Gdal API using a transformation defined by -ct string
 
  
 
kwargs = {
 
        "format": 'GTiff', 
 
        "resampleAlg": "cubic",
 
        "multithread":True 
 
}
 
  
 
kwargs.update({"coordinateOperation":"+proj=pipeline +step +inv +proj=utm +zone=34 +ellps=GRS80 +step +proj=push +v_3 +step +proj=cart +ellps=GRS80 +step +inv +proj=helmert +x=485.021 +y=169.465 +z=483.839 +rx=-7.786342 +ry=-4.397554 +rz=-4.102655 +s=0 +convention=coordinate_frame +step +inv +proj=cart 
 
+ellps=bessel +step +proj=pop +v_3 +step +proj=hgridshift +grids=.\\pipeline\\Slovakia_JTSK03_to_JTSK.gsb +step +proj=vgridshift +grids=.\\pipeline\\Slovakia_ETRS89h_to_Baltic1957.gtx +step +proj=krovak +lat_0=49.5 +lon_0=24.8333333333333 +alpha=30.2881397527778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel" })
 
kwargs.update({"srcSRS":"EPSG:25834","dstSRS":"EPSG:5514"})
 
  
 
warp = gdal.Warp(out_raster,in_raster,**kwargs)
 
  
 
  
 
The horizontal reprojection works fine but I can´t force the Gdal to use an elevation shift grid file, because my CRS does not have a vertical system defined.
 
  
 
In the Gdal utility, I can easily do this by using -vshift parameter, but I don´t know how to do this with Python API since osgeo.gdal.WarpOptions does not contain this parameter.
 
  
 
Is there a way, how to force vertical transformation?
 
  
 
  
 
Thank you.
 
  
 
Best regards,
 
Ing. Michal Girmala
 
 
 
+421 908 133 793             girmala at photomap.sk
 
+421 55 685 4396             www.photomap.sk
 
_ _
 
 
 
P Pred vytlačením tohto mailu prosím zvážte dopad na životné prostredie.
 
  
 _______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240726/da85d64f/attachment.htm>