<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: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=utf-8">
<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:inherit;
panose-1:0 0 0 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:#1F497D;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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-IE" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Hi Andrea,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Please disregard my last email the formula is incorrect.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Apologies..<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Richard<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<table style="FONT-SIZE: 9pt; HEIGHT: 321px; FONT-FAMILY: Arial, sans-serif; WIDTH: 500px; LINE-HEIGHT: 15px" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td style="FONT-SIZE: 9pt; FONT-FAMILY: Arial,sans-serif; WIDTH: 500px; COLOR: #004d44; PADDING-BOTTOM: 6px">
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44"><br>
——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44">Richard McDonnell MSc GIS, FME Certified Professional</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44"><i>FRM Data Management</i></span><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44">Oifig na nOibreacha Poiblí</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">Office of Public Works</span><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #565148">Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36
</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">Jonathan Swift Street, Trim, Co Meath, C15 NX36</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">M +353 87 688 5964 T +353 46 942 2409</span><br>
<a href="https://gov.ie/opw">https://gov.ie/opw</a><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">——<br>
To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonnell@opw.ie<br>
</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/
</span></td>
</tr>
</tbody>
</table>
<table style="FONT-SIZE: 9pt; FONT-FAMILY: Arial, sans-serif; WIDTH: 420px; LINE-HEIGHT: 15px" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td style="FONT-SIZE: 9pt; FONT-FAMILY: Arial,sans-serif; WIDTH: 420px; COLOR: #004d44; PADDING-BOTTOM: 6px">
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44"><br>
——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44">MSc GIS, FME Certified Professional</span><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #004d44">Oifig na nOibreacha Poiblí</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">Office of Public Works</span><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; FONT-WEIGHT: bold; COLOR: #565148">Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36
</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">Jonathan Swift Street, Trim, Co Meath, C15 NX36</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">——</span><br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #565148">M +353 87 688 5964 T +353 46 942 2409</span><br>
<font color="#565148"><a href="https://www.opw.ie">https://https://gov.ie/opw</a></font><br>
<br>
<span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">——<br>
</span><span style="FONT-FAMILY: Arial, sans-serif; COLOR: #004d44">Email Disclaimer:
<a href="https://www.opw.ie/en/disclaimer/">https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/</a></span>
</td>
</tr>
</tbody>
</table>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> QGIS-User <qgis-user-bounces@lists.osgeo.org>
<b>On Behalf Of </b>Richard McDonnell via QGIS-User<br>
<b>Sent:</b> 09 May 2023 09:27<br>
<b>To:</b> andrea antonello <andrea.antonello@gmail.com><br>
<b>Cc:</b> qgis-user@lists.osgeo.org<br>
<b>Subject:</b> Re: [Qgis-user] [pyqgis] create new raster<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Hi Andrea,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">If you are ok substituting the method you are using for GDAL, then you can use the following..<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:10.0pt;font-family:"inherit",serif;border:none windowtext 1.0pt;padding:0cm">gdal_calc.py -A input.tif --outfile=result.tif --calc="(A>=-2000)*(A<=1500)" --NoDataValue=-9999<o:p></o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:10.0pt;font-family:"inherit",serif;border:none windowtext 1.0pt;padding:0cm"><o:p> </o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">You can also implement this in QGIS using GDAL Calculator from the Processing Toolbox, It can also
be ran as a batch in QGIS or built into a model if you are looking to batch it or it’s something that’s regularly done.<o:p></o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">The output type is Float32, as you are using DTM, probably with decimal places a NoDATA value of -9999<o:p></o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Regards,<o:p></o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="vertical-align:baseline"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Richard<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="0" style="width:375.0pt">
<tbody>
<tr>
<td width="500" style="width:375.0pt;padding:0cm 0cm 4.5pt 0cm">
<p class="MsoNormal" style="line-height:11.25pt"><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
——</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
<b>Richard McDonnell MSc GIS, FME Certified Professional</b><br>
<i>FRM Data Management</i><br>
<br>
——<br>
<b>Oifig na nOibreacha Poiblí</b><br>
Office of Public Works<br>
<br>
</span><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36
</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">Jonathan Swift Street, Trim, Co Meath, C15 NX36</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">——</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">M +353 87 688 5964 T +353 46 942 2409</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
<a href="https://gov.ie/opw">https://gov.ie/opw</a><br>
<br>
——<br>
To send me files larger than 30MB, please use the link below <a href="https://filetransfer.opw.ie/filedrop/richard.mcdonnell@opw.ie">
https://filetransfer.opw.ie/filedrop/richard.mcdonnell@opw.ie</a><br>
<br>
Email Disclaimer: <a href="https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/">
https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/</a> <o:p>
</o:p></span></p>
</td>
</tr>
</tbody>
</table>
<p class="MsoNormal"><span style="display:none"><o:p> </o:p></span></p>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="0" style="width:315.0pt">
<tbody>
<tr>
<td width="420" style="width:315.0pt;padding:0cm 0cm 4.5pt 0cm">
<p class="MsoNormal" style="line-height:11.25pt"><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
——</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
<b>MSc GIS, FME Certified Professional</b><br>
<br>
——<br>
<b>Oifig na nOibreacha Poiblí</b><br>
Office of Public Works<br>
<br>
</span><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36
</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">Jonathan Swift Street, Trim, Co Meath, C15 NX36</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">——</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148">M +353 87 688 5964 T +353 46 942 2409</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
</span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#565148"><a href="https://www.opw.ie">https://https://gov.ie/opw</a></span><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:#004D44"><br>
<br>
——<br>
Email Disclaimer: <a href="https://www.opw.ie/en/disclaimer/">https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/</a>
<o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> QGIS-User <<a href="mailto:qgis-user-bounces@lists.osgeo.org">qgis-user-bounces@lists.osgeo.org</a>>
<b>On Behalf Of </b>andrea antonello via QGIS-User<br>
<b>Sent:</b> 09 May 2023 07:28<br>
<b>To:</b> <a href="mailto:qgis-user@lists.osgeo.org">qgis-user@lists.osgeo.org</a><br>
<b>Subject:</b> [Qgis-user] [pyqgis] create new raster<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hello,<o:p></o:p></p>
<div>
<p class="MsoNormal">I am trying to find out the best workflow to create a new raster.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">As an example I take an existing elevation model raster and loop over it to set values between 1500 and 2000 to novalue and write the result to a new raster.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">This is the only way I found to do so:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:green"># create new raster with novalues between 1500 and 2000</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dataType</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dtmLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">dataProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">().</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">dataType</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">1</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">crs</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsCoordinateReferenceSystem</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'EPSG:3003'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">params</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> = {<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'EXTENT'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dtmLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">extent</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(),<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'TARGET_CRS'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">crs</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'PIXEL_SIZE'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dtmLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">rasterUnitsPerPixelX</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(),<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'NUMBER'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">: -</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">9999.0</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:green"># 'OUTPUT_TYPE': dataType,</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'OUTPUT'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsProcessing</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">TEMPORARY_OUTPUT</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">}<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRaster</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">processing</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.run(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'qgis:createconstantrasterlayer'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">params</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)[</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'OUTPUT'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">]<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsRasterLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRaster</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'temp'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#A31515">'gdal'</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">dataProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">()<o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">block</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsRasterBlock</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dataType</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">cols</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">rows</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#AF00DB">for</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">row</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">in</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">range</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">rows</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">):<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#AF00DB">for</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">col</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">in</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">range</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">cols</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">):<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">point</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dtmLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">dataProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">().</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">transformCoordinates</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsPoint</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">col</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">row</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">),
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">transformType</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black"> </span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">value</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">res</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> =
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">dtmLayer</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">dataProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">().</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">sample</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#267F99">QgsPointXY</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">point</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">x</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(),
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">point</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">y</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">()),
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">1</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black"> </span><span style="font-size:10.5pt;font-family:"Courier New";color:#AF00DB">if</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">res</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">and</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">value</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> != -</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">9999.0</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#AF00DB">if</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">value</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> <
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">1000</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">or</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">value</span><span style="font-size:10.5pt;font-family:"Courier New";color:black"> >
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">2000</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black">
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">block</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">setValue</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">row</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">col</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">value</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">setEditable</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">True</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">writeBlock</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">block</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">,
</span><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">band</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">=</span><span style="font-size:10.5pt;font-family:"Courier New";color:#098658">1</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="line-height:14.25pt"><span style="font-size:10.5pt;font-family:"Courier New";color:#001080">newRasterProvider</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">.</span><span style="font-size:10.5pt;font-family:"Courier New";color:#795E26">setEditable</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">(</span><span style="font-size:10.5pt;font-family:"Courier New";color:blue">False</span><span style="font-size:10.5pt;font-family:"Courier New";color:black">)<o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">This code has two main issues:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">1. if I uncomment the line containing OUTPUT_TYPE, I am getting an error about the type passed. But I can't find the right type needed there, it should be the one taken from the original provider.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">2. the resulting raster is scrambled as if there was a shift in the setting of the values. But the QgsRasterBlock seems to be built correctly (rows, cols) and the values set properly (col, row).<br>
3. in the above example, the dtmLayer has an epsg 3033 crs and when loaded manually into QGIS, it is recognized. But when I read the layer's metadata crs with pyQGIS , it is not able to read it and tells me it is invalid. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Has anyone a hint about what I am doing wrong? <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Andrea<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</div>
</body>
</html>