站点插值作图warning，出图错误

黄小仙儿

用站点数据插值后画雨量分布图，出现警告：
d:\anaconda3\lib\site-packages\scipy\interpolate\rbf.py:257: LinAlgWarning: Ill-conditioned matrix (rcond=1.20757e-21): result may not be accurate.  self.nodes = linalg.solve(self.A, self.di)

把插值后的grid输出了一下

[[    0.             0.             0.         ...  8833.41735458   2553.85585213     0.        ]
[    0.             0.             0.         ...  5394.71074963  10134.48536777     0.        ]
[    0.             0.             0.         ...  7731.83825588  11421.4927206      0.        ] ... [  872.55770493  9633.00927353  2759.9814167  ...     0.      0.             0.        ]
[ 2794.32398796     0.             0.         ...     0.      0.             0.        ]
[ 2810.32857895     0.             0.         ...     0.      0.             0.        ]]
出来的图也不对。

查找了半天也不知道问出在哪里，我把脚本贴出来请大家帮我检查一下，感谢！

proj= ccrs.PlateCarree()  # 简写投影
plt.rcParams['font.sans-serif']=['SimHei']  #显示中文
fig=plt.figure(figsize=(4,2),dpi=500)#添加画布
ax = fig.subplots(1, 1, subplot_kw={'projection': proj})  # 创建子图
extent=[107.5,110.0,33.5,34.8]#限定绘图范围
#隐藏边框线
#plt.axis('off')
lon=df['经度']#读取站点经度
lat=df['纬度']#读取站点纬度
tem=df['小时降水量']#读取站点气温
olon=np.linspace(107.5,110.0,80)#设置网格经度
olat=np.linspace(33.5,34.8,80)#设置网格纬度

olon,olat=np.meshgrid(olon,olat)#网格化
func=Rbf(lon,lat,tem,function='linear')#定义径向基函数插值
tem_new=func(olon,olat)#获得插值后的网格气温
tem_new[tem_new <0 ] = 0
print(tem_new)

clevs = [0.1,10,25,50,100,250,500]
cdict = ['#A9F090','#40B73F','#63B7FF','#0000FE','#FF00FC','#850042']
my_cmap = colors.ListedColormap(cdict)
norm = mpl.colors.BoundaryNorm(clevs,my_cmap.N)
shp_path=r'D:\shp\bb.shp'
a_shapes=list(Reader(shp_path).geometries())
reader = Reader(shp_path)
enshicity = cfeat.ShapelyFeature(reader.geometries(), proj, edgecolor='k', facecolor='none')
ax.set_extent(extent, crs=proj)

contour = plt.contourf(olon,olat, tem_new,clevs,transform = ccrs.PlateCarree(),cmap=my_cmap,norm = norm)

1099221723

rbf的话，如果有两个数据的经纬度一样的话，就会出问题。

betasy

1. 把数据传上来
2. Rbf这个函数是谁的，引用的还是自己写的

黄小仙儿

1099221723 发表于 2022-1-25 16:19
rbf的话，如果有两个数据的经纬度一样的话，就会出问题。

那应该用什么呢？

黄小仙儿

betasy 发表于 2022-1-25 15:46
1. 把数据传上来
2. Rbf这个函数是谁的，引用的还是自己写的

是引用的，scipy的
数据不知怎么传不上去

1099221723

黄小仙儿 发表于 2022-1-25 16:37
那应该用什么呢？

需要对数据提前进行质控，检查数值是否合理，站点是否存在冲突，数据的空间分布是否合理

betasy

黄小仙儿 发表于 2022-1-25 16:42
是引用的，scipy的
数据不知怎么传不上去

dengshunqiang@outlook.com
数据发我邮箱，我看看

nightelfxu

我也出现了这样的问题，有的降水日可以正常出图，有的就插值有问题，请问最后怎么解决了？

edwardli

nightelfxu 发表于 2024-6-29 14:53
我也出现了这样的问题，有的降水日可以正常出图，有的就插值有问题，请问最后怎么解决了？

降水数据中可能有缺测值，比如999999。替换成nan，或者删除该记录就行。

nightelfxu

edwardli 发表于 2024-7-4 15:14
降水数据中可能有缺测值，比如999999。替换成nan，或者删除该记录就行。

李老师，我看了不是缺测的问题，刚刚我发了求助帖，老师有空给指点下