解决plt画遥感插值图出现空白块问题(使用interpolate.interp1d进行插值)

出现问题:使用插值画出遥感图后出现空白块问题(如下图所示):

很显然，问题出现在没有值(即值为Nan的问题)，我们首先尝试把所有的值都切换为1，这样保证有个基本值，但是无法解决这个问题。

解决办法:

首先我们查看出现问题的代码:

...
	f = interpolate.interp1d(X1, Y1, kind='cubic')
    x1min = X1.min()
    x1max = X1.max()
    for m in range(0, num2):
         try:
            o3_grid_whole[m][k] = f(m)
            day_grid[m][k] = m + DAY
            height_grid[m][k] = k*0.02
         except:
             o3_grid_whole[m][k] = 1
             day_grid[m][k] = np.nan
             height_grid[m][k] = np.nan
             ...
        

for循环里面的数值处理为对值进行处理，我们需要看到的事o3的值，然后我们将其默认值设置为1，但是画出来的图仍然为大白块。
解决，把代码换成下面这样即可解决:

	f = interpolate.interp1d(X1, Y1, kind='cubic')
    x1min = X1.min()
    x1max = X1.max()
    for m in range(0, num2):
         try:
            o3_grid_whole[m][k] = f(m)
            day_grid[m][k] = m + DAY
            height_grid[m][k] = k*0.02
         except:
             print ("出错的x为:{},x最小值为:{},x最大值为:{}".format(m, x1min, x1max))
             ll += 1
             o3_grid_whole[m][k] = 1
             day_grid[m][k] = m + DAY
             height_grid[m][k] = k*0.02

可以看到修改后的图如下图所示。

问题产生原因:interpolate.interp1d只能进行线段插值而不能实现直线插值。

(1)即我们使用前面已知的数据来构造一条线段，然后我们可以得到一个最小的x_min和x_max。当我们输入后续的值来进行插值的时候，我们只能输入在x_min到x_max之间的值，数学表达为[x_min,x_max]。所以我们在进行插值的时候我们会进行try_catch来避免代码出现报错停止运行。但是这样无法解释会出现白块的问题(第一个代码里面，即使出现问题，我们也会把其值设为1)。
(2)因为我们进行插值画图的时候使用的是散点图，所以我们画的其实是一个个的点，相当于我们需要的是x，y，z(就是图上的值，也可以理解为有值的点)。而我们在try_catch里面，只是把z的值给改了，并没有改变x，y。相当于无法定位到有白块的地方。而白块的地方，因为其xy所对应的z还是没有值的，所以获取不到值，所以得到的还是一个白块

码字不易，点个赞再走吧。