Scipy中griddata（）函数解决等高线不闭合的办法

 如以上两图所示，我们在应用matplotlib中的contour()绘制等值线的时候经常会因为边缘数据的缺少，而导致绘制出来的等高线的边缘并不闭合，然而我们有的时候要应用等值线对相关数据进行分析的时候，希望边缘是闭合的的等值线。怎么办呢？

我们先看一下我解决后的等值线如下：

 

 这是两幅相对应的图片，经过处理后的数据边缘部分是闭合的，要做到如上图一样，只需要在我们进行插值处理的时候添加一个参数即可。

scipy中的griddatah函数用来进行插值处理，函数形式如下：

scipy.interpolate.griddata(points, values, xi, method='linear', fill_value=nan, rescale=False)[source]¶

通常我们会进行如下处理：

zz=griddata(points=point,values=s,xi=(xx,yy),method="cubic")

 给定点和要插值的范围以及插值的方法，这时候注意参数fill_value，它默认是nan也就是空值，当进行插值处理时，默认在数值范围外的赋予nan，那么这一部分在会出的时候是不参与绘制的，也就是会出现我们上面出现的等高线不闭合的现象。

如果我们将fill_value的值设置称为参数values值得最小值，那么绘图区域就会铺满整个图面，等高线也就闭合了，代码如下，特做记录。

zz=griddata(points=point,values=s,xi=(xx,yy),method="cubic", fill_value=s.min())