站点数据插值到格点后边界留白的问题

小猴纸

尝试利用MATLAB将站点降水量图进行展示，
开始先是用了sactter散点图，如图1
现在想用等值线表示，插值之后分别做了填色和非填色两种图，出现了大片的空缺，如图2和3
想请教各位如何插值才能保证数据完全覆盖整个省呢？不胜感激

小猴纸

（填色图的标题画图的时候写错了请不要在意）
从散点图能看出来空白的部分其实是应该存在非零数据的，
不知道为什么插值之后就没有了
查了版块也没找到类似的问题和解决方案

小猴纸

附上插值的语句
[X1 Y1]=meshgrid([114.4:0.1:123.1],[33.8:0.1:39.2]);
[X,Y,Z] = griddata(lon,lat,value,X1,Y1,'natural'); %value为读取的站点数据

小猴纸

换了一下插值方式，解决了留白的问题，po一下
X1=linspace(114.4,123.1,100)';
Y1=linspace(33.8,39.2,100);
[X,Y,Z] = griddata(lon,lat,value,X1,Y1,'v4'); %区域内插值
[c,h]=contourf(X1,Y1,Z);


但是出现了新问题，
这样插值之后，精华帖里的掩膜法没有效果了
在省外的图上也有等值线分布
这样情况应该怎么解决呢？

/kfEspresso

gridfit函数
论坛有大神写了个这样函数，可以搜索参考下。