<div dir="ltr"><div class="gmail_default" style="font-size:small">Maris,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Many thanks for your detailed reply. My LIDAR files are not in LAS format - they are a mixture of ASC and TIF. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I spent a long time learning how to use QGIS and don't want to have to repeat the process with GRASS unless I have to. If there isn't a simple way to get r.contour to work from within the later versions of QGIS, then I'll keep on using the old version as the solution requiring the least effort. From your comments, it would seem that it is how QGIS imports the LIDAR data which has changed and this is why I see the problem I reported. I also realise that QGIS is a global application whereas my work is restricted to the UK using the Ordnance Survey grid so I can't expect a huge resource to look at a narrow application. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Cheers,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Dave</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 15 Feb 2021 at 10:08, Maris Nartiss <<a href="mailto:maris.gis@gmail.com">maris.gis@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Dave,<br>
QGIS hides a bit of GRASS complexity by making a guess for various<br>
parameters. As with any guess – sometimes it works, sometimes it is a<br>
miss (and user has no idea which is the case).<br>
<br>
To get contours out of LAS files:<br>
1) create a location with coordinate system matching one used by LAS<br>
files (be ware – you might need to know it in advance from metadata as<br>
LAS files quite often lack this information);<br>
2) create a mapset for the area of interest (could be whole region or<br>
a single file in case of parallel processing);<br>
3) start GRASS in newly created mapset;<br>
4) set up your computational region (this is most important part!)<br>
with g.region. Don't forget to choose appropriate resolution.<br>
a) if you know the extent in advance (e.g. from a map sheet grid) use that;<br>
b) if you don't know the extent in advance, use actual extent from the<br>
LAS file. I would advocate to use r.in.lidar -s and set the extent<br>
manually with g.region – you can “snap“ your raster to coordinates.<br>
5) import data with r.in.lidar;<br>
6) run r.contour on the map;<br>
7) export with v.out.ogr to Shapefile (#teamshapefile).<br>
<br>
Good luck,<br>
Māris.<br>
<br>
P.S. When you wander into area of 66000 LAS files occupying nice 14T<br>
on your disk, only a few adjustments need to be done + a bit of Python<br>
coding + a bit of cluster management :D<br>
</blockquote></div>