<html><head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div class="moz-cite-prefix">Hi there Richard,</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">The simple array format that your data
is in means that it can be read by GDAL using the ESRI HDR driver:</div>
<div class="moz-cite-prefix"><a class="moz-txt-link-freetext" href="https://gdal.org/drivers/raster/ehdr.html#raster-ehdr">https://gdal.org/drivers/raster/ehdr.html#raster-ehdr</a></div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">This just requires a sidecar file named
<i>clt02.hdr</i>, that contains the information for QGIS to
interpret the binary array. I modified the example in the page
above to create this, and saved it in the same folder as the data:<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">
<pre>ncols 4000</pre>
<pre>nrows 4000</pre>
<pre>cellsize 25</pre>
<pre>xllcorner 1589000.00</pre>
<pre>yllcorner 6135000.00</pre>
<pre>nodata_value 9999.000000</pre>
<pre>nbits 16</pre>
<pre>pixeltype SIGNEDINT</pre>
<pre>byteorder msbfirst</pre>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Load the original binary file into QGIS
and the hdr will be read automatically. It won't have a coordinate
system assigned, but that can be done using the Assign Projection
processing tool within QGIS, which will create a <i>clt02.prj</i>
file in the same folder containing the NZ CRS details.<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">The number of bits can be calculated
from the file size, but I have guessed that the data is a
(standard) signed integer, and indeed the values look reasonable
(30-37).<br>
</div>
<div class="moz-cite-prefix">Using the index file, equivalent .hdr
files could be built for all the NZ tiles<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">The excellent Serval plugin can then be
used to edit the raster values in place.</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Cheers,</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Andy<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">On 01/09/2023 07:28, Richard & Mary
via QGIS-User wrote:<br>
</div>
<blockquote type="cite" cite="mid:fc996004-baaf-19ab-a520-7e6fa2b2e4de@xtra.co.nz">
<div style="padding-bottom: 10px; padding-top: 5px;">
<div style="padding:12px; border:1px solid #8D3970;
background-color:#F7F9FA; color:#8D3970; font-size:14px;
line-height:22px; font-family: Calibri, Arial, Helvetica,
sans-serif;">
<strong>CAUTION:</strong> This e-mail originated outside the
University of Southampton.
</div>
</div>
<div>
<p>Greetings Team</p>
<p>I work as a Radio Engineer in the NZ Public Safety Sector and
a planning tool that I work with uses a land cover and
elevation rasters.</p>
<p>New Zealand is divided up in series of 100km tiles each tile
is 4000 x 4000 representing a 25m square.</p>
<p>The contents of these files are a binary number from 0 to 50
representing individual land cover categories, or 0 to approx
5000 representing the elevation<br>
</p>
<p><br>
</p>
<p>Each of these tiles has a Xmin and Xmax for Easting and a
Ymin and Ymax for Northing on the NZTM2000 grid.</p>
<p><b>I want to be able to display and edit these 25m squares to
fix some errors that affect our radio modelling.</b></p>
<p><br>
</p>
<p>I am using QGIS 3.28.10 and have tried a variety of Raster
imports but get a not supported format.</p>
<p>The file format is know as MSI Planet and is one of the
default formats used for land cover in the radio coverage
modelling tools.</p>
<p>I have attached a sample of the clutter file (this one was
31.25MB before zipping) I have also attached the index which
show where on the NZTM2000 grid this tile resides.</p>
<p>I also have the same files in a format used by EDX which
embeds the East - West North South Limits and the projection,
but these are generated from the initial plain raster.</p>
<p><br>
</p>
<p>I have read some information about using satelite image
stacks to generate land cover information, if I do that i need
to be able to output the plain 100km x 100km block<br>
</p>
<p><br>
</p>
<p>I would be interested to any help you can give me.</p>
<p>Richard Hutchinson<br>
</p>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
QGIS-User mailing list
<a class="moz-txt-link-abbreviated" href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a>
List info: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
Unsubscribe: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
</pre>
</blockquote>
<p><br>
</p>
<pre class="moz-signature" cols="72">--
Andy Harfoot
Skype: gdi_ajph
Teams: <a class="moz-txt-link-abbreviated" href="mailto:ajph@soton.ac.uk">ajph@soton.ac.uk</a>
Phone: +44 2380 590566
GeoData Institute
University of Southampton
Southampton
SO17 1BJ
<a class="moz-txt-link-abbreviated" href="http://www.geodata.soton.ac.uk">www.geodata.soton.ac.uk</a></pre>
</body>
</html>