Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I have a mosaic of a 1m LiDar derived DEM. I need to output a subset of the data as a 10m DEM. I am currently using the aggregate tool in ARCGIS 10 to produce a mean value for each new 10m pixel. Any advice on whether this is the best technique for such as task? Is the mean value the best approach with this type of data?

share|improve this question
1  
Can you further explain your workflow on why the data needs to be re-sampled to a lower resolution? There might be a better way to do it than lower the spatial resolution as the first step? –  MLowry Jul 5 '12 at 20:15
add comment

3 Answers

A common mistake (that I've made too) is to down-sample a raster using the resample tool with bilinear interpolation. See this answer for an explanation why this is not good. A raster can be down-sampled in three steps.

  1. The first step might not be required. Reproject the raster to the target extents. Use bilinear interpolation, and keep the output cell size the same as the input resolution (e.g., 1 m). Use the registration point to "snap" the raster corners to the projection. The output extents can be specified in the "Environments", and I suggest specifying the extents with a multiple of 10 m (or whatever resolution). These extents will control where the statistics are determined for the final raster.

  2. Perform Block Statistics (found in Spatial Analyst Tools > Neighborhood). Use a rectangle with 10 cells for both height and width, and choose "MEAN" for a statistic type. Try different shapes and types, if you want. The cell size is the down-sampling ratio.

  3. Since block statistics does not change the raster resolution, the last step is to Resample (found in Data Management Tools > Raster > Raster Processing). Choose 10 m, and use "NEAREST" to pick the block statistic in the centre of the cell.

Steps 2 and 3 can be replaced with's Curtvprice's suggestion to use the Aggregate tool, which will effectively yield the same results using rectangle means.

share|improve this answer
1  
whoa. I've been doing this stuff for 20 years and didn't know that bilinear resampling uses only 4 nearby data points! this definitely needs a higher profile. –  matt wilkie Jul 12 '12 at 20:44
add comment

You can cut a step out of the process by using the Aggregate tool, which avoids having to replicate all those summary value cells, and the Resample.

share|improve this answer
    
Good advice, less steps to perform. –  nadya Dec 4 '12 at 20:47
add comment

All the points mentioned above are important to take note of and I completely agree that bilinear resampling is quite problematic. Although, I am curious why nobody is discussing cubic convolution? The problem with using a block function is that the mean is quite irrelevant when the distribution is non-normal or multimodal, as expected with a lidar derived DEM.

If you have access to the original lidar data, just interpolate the data to the desired resolution using the "Topo to Raster" tool in ArcGIS. If you only have access to the 1m DEM raster, it seems like the best, albeit least efficient, method would be to convert the raster to points and use a thin plate or Bi-cubic spline. This would allow for the resample neighborhood to fit a nonlinear curve to the data.

Alternatively, you could smooth the 1m raster using a Gaussian kernel, approximating the size of your desired resample resolution (10x10), and then a bilinear resample would be much more appropriate. This approach would allow you have direct control over the smoothing parameter and would result in a "locally" normal distribution where the mean becomes relevant as an indicator of central tendency and a linear fit is supported.

share|improve this answer
    
The question owner is not seen for a long time here, but I have the same problem (my lidar rasters are 0.5 m), so I came to the discussion :) Regarding the converting raster to points, for ArcGIS it seems it's much more easy to handle big rasters than millions of points. –  nadya Dec 5 '12 at 20:27
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.