<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-2">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Webdings;
panose-1:5 3 1 2 1 5 9 6 7 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello all,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am trying to reproject elevation raster with gdalwarp in python Gdal API using a transformation defined by -ct string
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">kwargs = {<o:p></o:p></p>
<p class="MsoNormal"> "format": 'GTiff', <o:p></o:p></p>
<p class="MsoNormal"> "resampleAlg": "cubic",<o:p></o:p></p>
<p class="MsoNormal"> "multithread":True <o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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 <o:p></o:p></p>
<p class="MsoNormal" style="margin-left:36.0pt;text-indent:36.0pt">+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" })
<o:p></o:p></p>
<p class="MsoNormal">kwargs.update({"srcSRS":"EPSG:25834","dstSRS":"EPSG:5514"})<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">warp = gdal.Warp(out_raster,in_raster,**kwargs)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Is there a way, how to force vertical transformation?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="mso-ligatures:none">Best regards,<o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="mso-ligatures:none">Ing. Michal Girmala</span></b><span style="mso-ligatures:none"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-ligatures:none"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-ligatures:none">+421 908 133 793 </span><a href="mailto:girmala@photomap.sk"><span style="mso-ligatures:none">girmala@photomap.sk</span></a><span style="color:black;mso-ligatures:none"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-ligatures:none">+421 55 685 4396 </span><a href="https://photomap.photomap.sk/owa/redir.aspx?C=TkX_pwO3Cdy8-2srzfDp5ctFQ587EpI-9faJi0uprt0fCrd688PWCA..&URL=http%3a%2f%2fwww.photomap.sk%2f" target="_blank"><span style="mso-ligatures:none">www.photomap.sk</span></a><span style="color:#0563C1;mso-ligatures:none"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-ligatures:none">_ _<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-ligatures:none"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:Webdings;color:green;mso-ligatures:none">P</span><span style="color:green;mso-ligatures:none"> </span><a href="https://photomap.photomap.sk/owa/redir.aspx?C=tvo5neA_o2nG-tVUhKJQSd_SbCbRRCf3YD7l8laI8Q4fCrd688PWCA..&URL=http%3a%2f%2fthinkbeforeprinting.org%2f" target="_blank"><span style="color:green;mso-ligatures:none">Pred
vytlačením tohto mailu prosím zvážte dopad na životné prostredie.</span></a><span style="color:green;mso-ligatures:none"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>