regular spaced points… interpolation madness…

The geodata department of the city of Berlin offers a great portfolio of free geodata for everyone to use. One dataset is the result of some LIDAR measurements and is offered as a txt with x,y and d values. This is commonly known as xyz data but should not be mistaken as a simple whatever-delimited text file. This dataset is regular spaced and therefore can be threatened easily with QGIS. But let me first show you some interpolation results which can be produced as well.

Interpolational madness

If you want to create raster data from a point shapefile or a list of coordinates you normally think in the direction of TINs, IDW, kriging, and other more or less advanced interpolation techniques. QGIS offers TIN and IDW out of the box. As I was afraid to select the whole 2.5GB point shapefile as input, I’ve tried with a small subset of points. But first: what does the data look like?

xyd data from LIDAR@Berlin. Notice the right white column … 35.000.000 points…

Looking at it I was frightened… 35.000.000 data points. How can I convert it to a raster using my Mac and not this beautiful Tianhe-2? This is what it looks ones it is loaded into QGIS:

I was sure that interpolating this amount of data was a problem. My first try was done with a small subset. The interpolation methods IDW (inverse distance weighted) and TIN (triangulated irregular network) are placed under Raster->Interpolation after you’ve added the Interpolation plugin which is shipped with QGIS. As this workflow is not correct I’ll just show you some results for this small subset:

The IDW looked best to me, but it was also very lame in computing. And thinking about it once more: I don’t need to interpolate!

Grids are your friend

So I’ve searched the web and gis.stackexchange.com was once again my friend. Once you’ve added your points to QGIS go to Raster->Conversion->Rasterize an fill in the parameters you think are right:

parameters for raster creation

And the results for the subset were created fast and they look much better than the IDW method which took much longer:

I’ve learned: keep it simple to get best results 😉

4
Leave a Reply

avatar
3 Comment threads
1 Thread replies
0 Followers
 
Most reacted comment
Hottest comment thread
2 Comment authors
Jakobtbsprs Recent comment authors

This site uses Akismet to reduce spam. Learn how your comment data is processed.

  Subscribe  
newest oldest most voted
Notify of
tbsprs
Guest
tbsprs

I must admit I cannot follow you. Can you please explain the steps you took from the raw file (list of points) to a shapefile suitable for being used with the Rasterize plug-in? Also please reference the specific data source – I am not sure if you actually used the “ATKIS® Basis-DGM (2m-Rasterweite)” dataset: http://fbinter.stadt-berlin.de/fb/index.jsp?loginkey=showMap&mapId=k_dgm5@senstadt

tbsprs
Guest
tbsprs

PING

Jakob
Guest

@disqus_Sq3X3rzFfa:disqus: maybe you can give a short answer on @tbsprs:disqus question?

trackback

[…] Click to view slideshow. […]