<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:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;
        mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
pre
        {mso-style-priority:99;
        mso-style-link:"HTML Preformatted Char";
        margin:0cm;
        font-size:10.0pt;
        font-family:"Courier New";}
span.HTMLPreformattedChar
        {mso-style-name:"HTML Preformatted Char";
        mso-style-priority:99;
        mso-style-link:"HTML Preformatted";
        font-family:Consolas;
        mso-fareast-language:EN-US;}
span.EmailStyle23
        {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-AU" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">I just installed the latest update to Qgis (3.32.1-3) and noticed there was a GDAL update so I tried clipping an ecw to mask and it appears to work now.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<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="mso-fareast-language:EN-AU">From:</span></b><span lang="EN-US" style="mso-fareast-language:EN-AU"> Michael Dufty
<br>
<b>Sent:</b> Thursday, July 27, 2023 12:53 PM<br>
<b>To:</b> qgis-user@lists.osgeo.org<br>
<b>Subject:</b> RE: [Qgis-user] Clipping ecw to mask not working<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">To respond to Kirk and Andrea’s comments.<o:p></o:p></p>
<p class="MsoNormal">Yes it works fine on geotiff, including geotiff created be exporting the ecw using qgis.  Not always a desirable workaround when dealing with very big ecws.<o:p></o:p></p>
<p class="MsoNormal">It is on Windows 10.  <o:p></o:p></p>
<p class="MsoNormal">GDAL 3.7.1-1 from OSGEO installer<o:p></o:p></p>
<p class="MsoNormal">The ecw displays fine in QGIS.  Output layer is set to Byte to match the source
<o:p></o:p></p>
<p class="MsoNormal">Would be good to know if it works for anyone else.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I thought I would try with a QGIS/GDAL generated ecw  to see if that helps, but don’t seem to be able to generate an ecw either.<o:p></o:p></p>
<p class="MsoNormal">If I use the export layer command to generate one it says it ran successfully but no file is generated.<o:p></o:p></p>
<p class="MsoNormal">If I use translate(convert format) it does give error messages as below.  Not sure if it is the same issue.<o:p></o:p></p>
<p class="MsoNormal">Is GDAL now using the newer ecw sdk that doesn’t support unlicensed writing?  Might need to go back to Manifold for these things.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">QGIS version: 3.32.1-Lima</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">QGIS code revision: 4a078ecf63</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">Qt version: 5.15.3</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">Python version: 3.9.5</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">GDAL version: 3.7.1</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">GEOS version: 3.12.0-CAPI-1.18.0</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">PROJ version: Rel. 9.2.1, June 1st, 2023</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#777777;mso-fareast-language:EN-AU">PDAL version: 2.5.5 (git-version: a7569c)</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">Algorithm started at: 2023-07-27T12:43:23<o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="mso-fareast-language:EN-AU">Algorithm 'Translate (convert format)' starting…</span></b><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">Input parameters:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Courier New";mso-fareast-language:EN-AU">{ 'COPY_SUBDATASETS' : False, 'DATA_TYPE' : 0, 'EXTRA' : '', 'INPUT' : 'W:/Catalyst/GIS Common/Spatial Data/Aerial Imagery/croptest.tif', 'NODATA' : None, 'OPTIONS' : '',
 'OUTPUT' : 'W:/Catalyst/GIS Common/Spatial Data/Aerial Imagery/croptest.ecw', 'TARGET_CRS' : None }</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">GDAL command:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Courier New";mso-fareast-language:EN-AU">gdal_translate -of ECW "W:\\Catalyst\\GIS Common\\Spatial Data\\Aerial Imagery\\croptest.tif" "W:/Catalyst/GIS Common/Spatial Data/Aerial Imagery/croptest.ecw"</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">GDAL command output:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:red;mso-fareast-language:EN-AU">ERROR 1: None of ECW_ENCODE_KEY and ECW_ENCODE_COMPANY were provided.
<br>
Both are required.</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:red;mso-fareast-language:EN-AU">Process returned error code 1</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">Execution completed in 0.10 seconds<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">Results:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Courier New";mso-fareast-language:EN-AU">{'OUTPUT': 'W:/Catalyst/GIS Common/Spatial Data/Aerial Imagery/croptest.ecw'}</span><span style="mso-fareast-language:EN-AU"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-AU">Loading resulting layers<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:red;mso-fareast-language:EN-AU">The following layers were not correctly generated.<br>
• W:/Catalyst/GIS Common/Spatial Data/Aerial Imagery/croptest.ecw<br>
You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.</span><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal"><b><span style="color:#1F497D;mso-fareast-language:EN-GB">Michael Dufty</span></b><span style="color:#1F497D;mso-fareast-language:EN-GB"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<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="mso-fareast-language:EN-AU">From:</span></b><span lang="EN-US" style="mso-fareast-language:EN-AU"> QGIS-User <<a href="mailto:qgis-user-bounces@lists.osgeo.org">qgis-user-bounces@lists.osgeo.org</a>>
<b>On Behalf Of </b>Kirk Schmidt via QGIS-User<br>
<b>Sent:</b> Thursday, July 27, 2023 12:00 AM<br>
<b>To:</b> <a href="mailto:qgis-user@lists.osgeo.org">qgis-user@lists.osgeo.org</a><br>
<b>Subject:</b> Re: [Qgis-user] Clipping ecw to mask not working<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p>Hi Michael:<o:p></o:p></p>
<p>Enhance Compression Wavelet (ECW) is a proprietary compression format.  I would try exporting your image as a geotiff to see if it is being read correctly and does not have internal write protection on it.<o:p></o:p></p>
<p>Kirk Schmidt<o:p></o:p></p>
<div>
<p class="MsoNormal">On 7/26/2023 2:02 AM, Michael Dufty via QGIS-User wrote:<o:p></o:p></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<p class="MsoNormal">I'm having issues trying to clip ecw images to a mask layer.<o:p></o:p></p>
<p class="MsoNormal">Using raster-extraction-clip raster to mask layer.<o:p></o:p></p>
<p class="MsoNormal">Default parameters.<o:p></o:p></p>
<p class="MsoNormal">The tool runs extremely quickly (successfullyin 0.12 seconds), but produces a raster entirely black (0,0,0 for all pixels).  With the correct extent of the mask layer.<o:p></o:p></p>
<p class="MsoNormal">If I tick the box to create an alpha output band it produces an entirely transparent result.<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">Source and mask layers both in EPSG 28350 <o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">Anyone else having this issue?  I cannot see any bug reports. 
<o:p></o:p></p>
<p class="MsoNormal">Tried in 328.8, 3.28.9, 3.32.0 and 3.32.1 with same result.<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">Seems to work fine if input raster is a tif.<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="mso-fareast-language:EN-AU"><o:p> </o:p></span></p>
<pre>_______________________________________________<o:p></o:p></pre>
<pre>QGIS-User mailing list<o:p></o:p></pre>
<pre><a href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a><o:p></o:p></pre>
<pre>List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre>
<pre>Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre>
</blockquote>
<pre>-- <o:p></o:p></pre>
<pre>Kirk Schmidt, MScF, BScF, RPF<o:p></o:p></pre>
<pre>General Manager<o:p></o:p></pre>
<pre>Nortek Resource Solutions Inc.<o:p></o:p></pre>
<pre>26 Church Road<o:p></o:p></pre>
<pre>Sutherlands River, NS<o:p></o:p></pre>
<pre>B0K 1W0<o:p></o:p></pre>
<pre>Tel (902) 922.3607<o:p></o:p></pre>
<pre>Email: <a href="mailto:kirk@nortekresources.com">kirk@nortekresources.com</a><o:p></o:p></pre>
<pre>Web: <a href="http://www.nortekresources.com">www.nortekresources.com</a><o:p></o:p></pre>
</div>
</body>
</html>