如何解决python:scipy.interpolate.griddata 导致脚本结束
我正在尝试使用 scipy interpolate griddata 来插入经纬度数据点的常规网格。但是,当我尝试运行脚本时,我打印了“压缩”行,然后暂停,然后脚本结束而不打印任何其他内容或绘制等高线图。没有显示错误消息。我究竟做错了什么。请参阅下面的相关代码部分。
general_file = r'data\l3bx\netcdf\r0024\ir2_20160815_114117_232_l3bx_v10.nc'
file_data = Dataset(general_file) #Reading in the data
m = Basemap(projection='moll',resolution=None,lat_0=0,lon_0=0)
axis1 = file_data.variables['axis1'][:].data #Extracting values neeed
axis2 = file_data.variables['axis2'][:].data
rad = file_data.variables['radiance'][0].data
lon = file_data.variables['lon'][:][0].data
lat = file_data.variables['lat'][:][0].data
xi = np.arange(-180,190,10)
yi = np.arange(-90,100,10)
grid_x,grid_y = np.meshgrid(xi,yi)
rad = rad.flatten()
lon = lon.flatten()
lat = lat.flatten()
points = np.array(list(zip(lon,lat)))
print('Zipped')
grid_z = interpolate(points,rad,grid_x,grid_y)
plt.contourf(grid_x,grid_y,grid_z,25,cmap='Greys_r')
plt.colorbar()
plt.show()
print('End')
def interpolate(points,grid_y):
grid_z = griddata((lon,lat),(grid_x,grid_y),method='nearest')
print('Interpolated')
return grid_z
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。