Resample raster using a Gaussian filter

313
5
Jump to solution
09-02-2022 02:50 PM
nik_geo
New Contributor III

I want to resample a satellite image using a Gaussian filter. By resample I mean to increase the pixel size from 15m to 460m. Is there any tool I can use in ArcMap 10.8? 

0 Kudos
1 Solution

Accepted Solutions
nik_geo
New Contributor III

ArcGIS Pro does not have (yet) this functionality. Or at least I didn't find it.

View solution in original post

0 Kudos
5 Replies
DanPatterson
MVP Esteemed Contributor

Resample

Resample—ArcGIS Pro | Documentation

Resample function—ArcGIS Pro | Documentation

should be what you are looking for


... sort of retired...
0 Kudos
nik_geo
New Contributor III

So I tried the Resampling function using a Gaussian blur and Gaussian blur Plus. The issue is that ArcGIS Pro doesn't give me the option to set the the std. Probably I should have mentioned that in my question, but I would like to investigate the resampling results using a Gaussian blur with various std values (from 0.2 to 0.9 with step 0.05). So I am guessing I need to do this programmatically as I couldn't find any GIS software that gives me this option.

0 Kudos
DanPatterson
MVP Esteemed Contributor

Before you reinvent the wheel, have a look at

scipy.ndimage.gaussian_filter — SciPy v1.9.1 Manual

 


... sort of retired...
0 Kudos
nik_geo
New Contributor III

I tested the algorithm and the output is an array with the same dimensions as the input while my goal is to have fewer dimensions. Unless I made a mistake, using the function scipy.ndimage.gaussian_filter I cannot change the pixel size (or the dimensions of the array in python's terms).

To put it differently, I am trying to simulate coarse data (resample panchromatic band) as though they were measured with a coarse PSF.
To do this I need to apply a transfer function (TF; e.g., Gaussian) to the fine data, but with a very large width.
This produces the coarse data. The PSF (point spread function) is the Gaussian filter. 

 

Also, the tool Gauss blur (plus) I think works like this:

From ArcGIS documentation, "Applies a Gaussian convolution to the source raster and calculates pixel value using the distance-weighted value of four nearest pixels from the blurred raster. It is appropriate for removing noise in resampled data and for down-sampling to a larger pixel size."

It simply apply the filter to blur the image then go ahead to resample it (in this case, from their definition, bilinear interpolation)

From ArcGIS "Bilinear Interpolation—Calculates the pixel value using the distance-weighted value of the four nearest pixels. This method is computationally efficient to process".

For a more detailed explanation of what I want to achieve, please refer to the paper 'The effect of point spread function on downscaling continua'. The file size of the pdf is 35mb and I can't upload it. 

0 Kudos
nik_geo
New Contributor III

ArcGIS Pro does not have (yet) this functionality. Or at least I didn't find it.

0 Kudos