How can I run a Low Pass Filter over a Digital Elevation Model to Reduce Noise?

6 次查看（过去 30 天）
Jake Collins-May
Jake Collins-May 2017-7-27
编辑： Chad Greene 2017-7-28
I want to run an elevation profile tool over a Digital Elevation Model (DEM) in ArcGIS. The DEM is noisy so the profile returned is jagged. I am looking to run a Low Pass Filter over the DEM in order to "blur" it, reducing the noise, and smoothing the extracted profile.
The DEM is: 17997x8369 pixels 1 Band 32 bit floating point
The code I have been using to try and achieve this is:
I = imread('DEM.tif');
G = fspecial('gaussian',[5 5],2);
Ig = imfilter(I,G,'same');
I2 = im2double(Ig);
imwrite(I2, 'DEM_5_2.tif')
However this simply returns a black box. How would I go about running a low pass filter over this image correctly?
  2 个评论
Chad Greene
Chad Greene 2017-7-28
If the issue is random noise, you may consider a moving median filter like medfilt2, which would remove outliers rather than spreading them out.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Chad Greene
Chad Greene 2017-7-28
编辑：Chad Greene 2017-7-28
It looks like imwrite will not write double to tiff format. Try something like
I2 = im2uint16(Ig);
imwrite(I2, 'DEM_5_2.tif')

类别

在 Help Center 和 File Exchange 中查找有关 Image Category Classification 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by