Python完美白化

偷影子的人

楼主，我想画在兰伯特投影的图上，但是我改了程序之后，maskout按你说的改了，然后主程序里改了clip那里，还有 m=Basemap(projection = 'laea', lat_0 = 33, lon_0 = 102., lat_ts = 33, width = 7000000, height = 5000000, resolution = 'l'),
但我这样写了之后，画出图来是空白的。
是什么原因呢？我看网上有人说“数据的投影方式也要改”，我对这个不是很明白，能问一下您具体应该怎么做吗？

平流层的萝卜

偷影子的人 发表于 2020-4-14 14:46
楼主，我想画在兰伯特投影的图上，但是我改了程序之后，maskout按你说的改了，然后主程序里改了clip那里， ...

数据的投影方式这一块，你是否指 xx,yy=m(xx,yy) 这里？确保这个要写上去。

偷影子的人

平流层的萝卜 发表于 2020-4-14 19:03
数据的投影方式这一块，你是否指 xx,yy=m(xx,yy) 这里？确保这个要写上去。

我加上了，还是不行。我把代码给您看吧，太不好意思了。

1. data=ncdata.variables['pres'][0,:,:]
   
2. 
   
3. lat=ncdata.variables['lat'][:]
   
4. lon=ncdata.variables['lon'][:]
   
5. lon1,lon2=lon[0],lon[-1]
   
6. lat1,lat2=lat[-1],lat[0]
   
7. nx=data.shape[1];ny=data.shape[0]
   
8. 
   
9. m = Basemap(projection = 'laea', lat_0 = 33, lon_0 = 102., lat_ts = 33, width = 7000000, height = 5000000, resolution = 'l')
   
10. 
    
11. 
    
12. 
    
13. xx,yy=np.meshgrid(lon,lat)
    
14. xx,yy=m(xx,yy)
    
15. 
    
16. m.readshapefile('E:/research/study/python/maskout/E2/country1','whatevername',color='gray')
    
17. 
    
18. minval,maxval=int(numpy.amin(data)),int(numpy.amax(data))+1
    
19. cs = m.contourf(xx,yy,data,range(minval,maxval),cmap= plt.cm.get_cmap('jet'))
    
20. bar=m.colorbar(cs)
    
21. bar.set_ticks(range(minval-1,maxval,40))
    
22. bar.set_ticklabels(range(minval-1,maxval,40))
    
23. clip=maskout.shp2clip(cs,ax,m,'country1','China')

复制代码

Adiligirl

mark一下，希望以后会用到

笺疏

发生local variable 'clip' referenced before assignment错误怎么办啊

一大碗年糕

感谢大佬

平流层的萝卜

偷影子的人 发表于 2020-4-15 09:28
我加上了，还是不行。我把代码给您看吧，太不好意思了。

你好，出的图是空白的，我觉得可能还有这两个问题：1. ax的设置不对，就是说你没有画在ax里边。2. basemap里关于这个投影方式的参数设置有问题。建议你两个都再复查一下

偷影子的人

平流层的萝卜 发表于 2020-4-19 16:45
你好，出的图是空白的，我觉得可能还有这两个问题：1. ax的设置不对，就是说你没有画在ax里边。2. basema ...

好的。我现在画出来了，谢谢楼主！

平流层的萝卜

偷影子的人 发表于 2020-4-20 15:11
好的。我现在画出来了，谢谢楼主！

最后是什么问题？

一日的天光

感谢分享，燃眉之急！