This function performs a 2D gaussian-weighted moving-window averaging filter on gridded datasets. It's designed to simplify the process of converting resolution to pixels and figuring out what sigma value to use. Click on the light bulb icon to the right on this page to see examples of use.
Chad Greene (2020). filt2 2D geospatial data filter (https://www.mathworks.com/matlabcentral/fileexchange/61003-filt2-2d-geospatial-data-filter), MATLAB Central File Exchange. Retrieved .
thank you !! Then, I have a problem. We know geospatial data gridded in longitude and latitude, they have different resolution in the two direction, because of the sphere of earth. So, can we change the code for a better use?
Why is the wavelength defined at e-folding scale of power instead of at half-power point?
Why is the cut-off length defined at e-folding scale instead of half-power point?
Very nice job!
how to filter the spatial field with 15km*20km(X-direction and Y-direction,respectively) in the example?
Added bandpass and bandstop filters.
Inspired by: Generate spatial data, moving_average v3.1 (Mar 2008), smooth2a, Noise-Power Spectrum, Generating Gaussian Filter 2D Matrix, filter1, Radially averaged surface roughness/topography power spectrum (PSD), scatstat2 2D local statistics, ndnanfilter.m