<div dir="ltr"><div><div><div><div>Thanks Ben,<br><br></div>I will give that a try. I am finding that the Interpolation algorithm is taking forever on a small data set.<br><br></div>Doesn't the Postgres PDAL extension handle point data like this? Would there be a benefit of putting the data in PostGIS?<br><br></div>Thanks,<br><br></div>Calvin<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Apr 20, 2018 at 1:10 PM, Benjamin Ducke <span dir="ltr"><<a href="mailto:benducke@fastmail.fm" target="_blank">benducke@fastmail.fm</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Calvin,<br>
<span class=""><br>
On 20/04/18 18:46, Richard McDonnell wrote:<br>
> Calvin,<br>
> <br>
> Under Processing Toolbox, Type interpolation into the search bar, the<br>
> tool you want can be found there.<br>
<br>
</span>None of the typical GIS interpolation tools will be<br>
of much help here, since they have all been optimized<br>
to work with relatively sparse data. They will be<br>
too inefficient for high-resolution geomagnetic data.<br>
<br>
You data is already dense, with some gaps. So you need<br>
something that is fast and efficient for small<br>
interpolation radii, and will not produce too many<br>
artefacts.<br>
<br>
If you can use the GRASS7 plug-in from within QGIS,<br>
than this is a good processing path:<br>
<br>
1. Rasterize your data points, using a cell size<br>
small enough so that you don't have too much information<br>
loss (i.e. not too many occurrences of more than one data<br>
point per cell). This might take a little trial and error<br>
to find good settings. You should end up with a raster<br>
where you have gaps of max. 1-3 empty cells between cells<br>
with measurement data.<br>
<br>
Now you need to fill these relatively small gaps.<br>
<br>
2. Use r.fill.stats from current GRASS 7.4:<br>
<br>
<a href="https://grass.osgeo.org/grass74/manuals/r.fill.stats.html" rel="noreferrer" target="_blank">https://grass.osgeo.org/<wbr>grass74/manuals/r.fill.stats.<wbr>html</a><br>
<br>
... because that has been designed exactly for the<br>
purpose. Set the interpolation radius as small as<br>
possible, so that it just manages to fill the gaps.<br>
With a small radius, r.fill.stats is very fast and<br>
memory efficient, even on very large raster datasets.<br>
<br>
Note that r.fill.stats is part of GRASS 7 only since<br>
version 7.4. If your QGIS ships with an earlier<br>
version, then you might have to update your GRASS 7<br>
binaries manually.<br>
<br>
Best,<br>
<br>
Ben<br>
<span class=""><br>
> Regards,<br>
> <br>
> Richard.<br>
> <br>
> On 20/04/2018 16:42, C Hamilton wrote:<br>
>> I don't know how many of you have worked with archaeological<br>
>> magnetormetry data, but I have been asked to help with magnetometry<br>
>> data of the ancient Mound Builders in North America. I am looking for<br>
>> some suggestion on how to best turn the irregular magnetic data points<br>
>> into an image.<br>
>><br>
>> I see an Interpolation plugin that perhaps is the solution, but I also<br>
>> see that it is not part of QGIS 3. Is this because it has not been<br>
>> ported yet, or is not going to be porter, or is there a better solution?<br>
>><br>
>> I don't know if any of the Lidar tools would work for this or not.<br>
>> LASTools is a commercial product and the locked version only works<br>
>> with small data sets and adds noise unless you purchase a license. At<br>
>> least that is how it was the last time I used it.<br>
>><br>
>> Thanks for any suggestions.<br>
>><br>
>> Calvin<br>
>><br>
>><br>
>> ______________________________<wbr>_________________<br>
>> Qgis-user mailing list<br>
>> <a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><br>
>> List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a><br>
>> Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a><br>
> <br>
> -- <br>
</span>> *_Richard McDonnell_*<br>
> *GIS Specialist PgD GIS AssocSCSI*<br>
> *OPW FRM Data Management*<br>
> *52 Stephens Green, Dublin 2.*<br>
> TEL: 01 6476543<br>
> <br>
> <<a href="http://www.opw.ie/" rel="noreferrer" target="_blank">http://www.opw.ie/</a>><br>
> <br>
> <br>
> <br>
> Email Disclaimer: <a href="http://www.opw.ie/en/disclaimer/" rel="noreferrer" target="_blank">http://www.opw.ie/en/<wbr>disclaimer/</a><br>
<span class="">> <br>
> <br>
> ______________________________<wbr>_________________<br>
> Qgis-user mailing list<br>
> <a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><br>
> List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a><br>
> Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a><br>
> <br>
<br>
<br>
<br>
-- <br>
</span>Dr. Benjamin Ducke<br>
Deutsches Archäologisches Institut (DAI)<br>
Zentrale Berlin, IT-Referat<br>
* Projekt "Stunde Null" *<br>
<div class="HOEnZb"><div class="h5">______________________________<wbr>_________________<br>
Qgis-user mailing list<br>
<a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><br>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a><br>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-user</a></div></div></blockquote></div><br></div>