Thursday 29 December 2016

Creating DEM from point data in CSV files using ArcGIS Desktop?


I have a lot of csv files which contain point data containing height, aspect and slope measurements at a resolution of 12,5 meters.


I need to create a dem from this data and have ArcGIS 10 with spatial analyst and 3d Analyst. I did this before with 9.3 but I am not sure if my methodology is correct.




  1. Create point shapefiles from the csv file

  2. interpolate to raster (here I never know whether to use spline, kriging or some other method)?


Last time I did something wrong and rounded all the values up to 1m even though the data was to cm accuracy so this time I want to do it right!



Answer



Because these are regular points, you effectively already have a DEM; it's just in a different format than ArcGIS likes. This makes two different strategies available to you:




  1. Convert the data into a format ArcGIS can handle. One way is to set up a raster extent and cellsize that (i) cover your DEM and (ii) situate each point near the middle of its cell. Import your DEM as a point layer and convert that directly, with no interpolation, into a grid format. (Each cell acquires the value of the unique point lying within it.) From now on, interpolation will be automatically performed whenever you resample the DEM, when aggregating, resizing, or reprojecting it. This gives you access to nearest-neighbor (not recommended), bilinear, and cubic convolution methods. You can create three separate grids this way: one for elevation, one for slope, one for aspect. For this to work, though, you must be using exactly the same coordinate system (including the projection) used to create the DEM in the first place, so that the points are truly evenly spaced and parallel to the coordinate axes.





  2. Pretend the points are "irregular" and interpolate only the elevations using any method you can (splines, IDW, kriging, etc). Derive new slope and aspect DEMs from the interpolated elevations. Unless you are using a Topo2DEM approach and have additional information you can supply--such as stream beds, elevation contours, peaks, or such material--then it is likely this method will merely produce meaningless artifacts that reflect the interpolation method as much as they do the original data. It will also require much more computation than the first method (kriging may even be impossible because of this).




No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...