来源：https://www.mathworks.com/matlabcentral/fileexchange/61003-filt2-2d-geospatial-data-filter

这个作者其实是基于

ndnanfilter.m - File Exchange - MATLAB Central (mathworks.cn)

这个工具包写的，里面用了fspecial函数构造高斯滤波窗口（Gaussian window）

ndnanfilter.m默认采用矩形窗口，在这里可以构建任意常用的滤波窗口

f = fspecial('gaussian',2*ceil(2.6*sigma)+1,sigma); % sigma是频率/波数 2*pi*sigma 是周期/波长

这里2*ceil(2.6*sigma)+1可以认为是影响范围

% https://stackoverflow.com/questions/60153406/standard-deviation-of-a-discrete-gaussian-filter

这里的作者用的是2*ceil(3*sigma)+1

本来想构建Lanczos二维滤波器，只需要修改fspecial一行即可

实现中遇到些问题，赶时间先采用高斯滤波吧

Lanczos二维滤波器：

https://en.wikipedia.org/wiki/Lanczos_resampling

https://www.mathworks.com/matlabcentral/answers/236122-how-to-do-lanczos-filtering-of-3-d-geospatial-data