<html><head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">Hi Ivan,<br><br>For my application I also had to set the geotransform of the output raster. Just in case you might find it useful, here is the link to my dump_raster function.<br><br><a href="https://github.com/wradlib/wradlib/blob/81e41376088c0ad5551c7b7b7eca4423d79b89ec/wradlib/zonalstats.py#L298">https://github.com/wradlib/wradlib/blob/81e41376088c0ad5551c7b7b7eca4423d79b89ec/wradlib/zonalstats.py#L298</a><br><br>HTH,<br>Kai <br><br><div class="gmail_quote">Am 24. Juli 2019 21:38:34 MESZ schrieb Ivan Lucena <ivan.lucena@outlook.com>:<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>Even, Kai,</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><br>
</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Thank you both for your reply. </div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Adding a band.FlushCache() did not fixed the problem but it was a good point anyway. Thanks Kai.</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
I tried to use the gdal.Rasterize() instead of gdal.RasterizeLayer() without success. There example on autotest varies quit a bit and I could not match one that works.</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
The geojson and geotiff files have the exactly same SRS (EPSG:4326) and they overlap perfectly fine on QGIS. But it seems like the issue might be related to SRS anyway.</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
Running the script from the Spyder IDE I was missing a warning message that now I can see by running "python script.py" on command line:</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace;">Warning 1: Failed to fetch spatial reference on layer MEMORY_LAYER to build transformer, assuming matching coordinate systems.</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
If I understood what Even said, I could fix that issue before calling gdal.RasterizeLayer() or use gdal.Rasterize() with the a transformation operation to fix the problem.</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
I am going to work on that. <span title="">Thank you again so very much,</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
Ivan</div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<br>
</div>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Even Rouault <even.rouault@spatialys.com><br>
<b>Sent:</b> Wednesday, July 24, 2019 10:57 AM<br>
<b>To:</b> gdal-dev@lists.osgeo.org <gdal-dev@lists.osgeo.org><br>
<b>Cc:</b> Ivan Lucena <ivan.lucena@outlook.com><br>
<b>Subject:</b> Re: [gdal-dev] Python RasterizeLayer and gdal_rasterize</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Ivan,<br>
<br>
what comes to mind is a SRS mismatch issue. From what I can see <br>
gdal.RasterizeLayer() takes a pfnTransformer / pTransformArg for that, but I'm <br>
not sure if they are probably mapped to SWIG to something practical.<br>
<br>
You might try instead with the gdal.Rasterize() "utility as a C function" that <br>
is higher level and should take care of that automagically<br>
<br>
usage examples in autotest/utilities/test_gdal_rasterize_lib.py<br>
<br>
Even<br>
</div></span></font></div></blockquote></div><br>-- <br>Diese Nachricht wurde von meinem Android-Gerät mit K-9 Mail gesendet.</body></html>