Python完美白化

letaotao22

谢谢，学习了

cashia371081

厉害，厉害！！

气象小小小白

lovechang1314 发表于 2016-3-22 17:02
简单测试了以下，可以用。
有个问题，projection='cyl' 应该是不需要转换坐标的吧？
如果用'stere'这个时 ...

散点图怎么白化

Jareau

请问为什么所有的地方都被覆盖掉了
还有错误/home/jareau/miniconda3/envs/123/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py:782: MatplotlibDeprecationWarning: Passing the minor parameter of set_ticks() positionally is deprecated since Matplotlib 3.2; the parameter will become keyword-only two minor releases later.
  return super(GeoAxes, self).set_xticks(xticks, minor)/home/jareau/miniconda3/envs/123/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py:829: MatplotlibDeprecationWarning: Passing the minor parameter of set_ticks() positionally is deprecated since Matplotlib 3.2; the parameter will become keyword-only two minor releases later.  return super(GeoAxes, self).set_yticks(yticks, minor)

Out[3]:

<matplotlib.colorbar.Colorbar at 0x7fc74c125b50>


/home/jareau/miniconda3/envs/123/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py:388: MatplotlibDeprecationWarning: The 'inframe' parameter of draw() was deprecated in Matplotlib 3.3 and will be removed two minor releases later. Use Axes.redraw_in_frame() instead. If any parameter follows 'inframe', they should be passed as keyword, not positionally.  inframe=inframe)

1. from scipy.interpolate import  Rbf
   
2. import matplotlib.pyplot as plt
   
3. import cartopy.crs as ccrs
   
4. import cartopy.feature as cfeature
   
5. from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
   
6. import cartopy.mpl.ticker as cticker
   
7. import cartopy.io.shapereader as shpreader
   
8. import maskout
   
9. fig_2=plt.figure(figsize=(15,10))
   
10. proj = ccrs.PlateCarree(central_longitude=115)
    
11. leftlon,rightlon,lowerlat,upperlat=(70,137.5,15,52.5)
    
12. lon_formatter = cticker.LongitudeFormatter(zero_direction_label=False)
    
13. lat_formatter = cticker.LatitudeFormatter()
    
14. f= fig_2.add_axes([0.1,0.1,0.9,0.9],projection = proj)
    
15. f.xaxis.set_major_formatter(lon_formatter)
    
16. f.yaxis.set_major_formatter(lat_formatter)
    
17. f.set_extent([leftlon, rightlon, lowerlat, upperlat], crs=ccrs.PlateCarree())
    
18. f.add_feature(cfeature.COASTLINE.with_scale('50m'))
    
19. f.add_feature(cfeature.LAKES, alpha=0.5)
    
20. f.set_xticks(np.arange(leftlon,rightlon+5,5), crs=ccrs.PlateCarree())
    
21. f.set_yticks(np.arange(lowerlat,upperlat+5,5), crs=ccrs.PlateCarree())
    
22. china = shpreader.Reader('country1.dbf').geometries()
    
23. f.add_geometries(china,ccrs.PlateCarree(),facecolor='none', edgecolor='black',zorder = 1)
    
24. f.set_title('EU&tem',loc='left',fontsize =20)
    
25. c2=f.contourf(slon,slat,r2_new[:],levels=np.arange(-1,1.1,0.1), zorder=0,cmap=plt.cm.coolwarm,transform=ccrs.PlateCarree())
    
26. c2_=f.contourf(slon,slat,p2_new,[np.min(p2_new),0.05,np.max(p2_new)],zorder=1,hatches=['.',None],colors="none",transform=ccrs.PlateCarree())
    
27. clip=maskout.shp2clip(c2,f,'country1','China')
    
28. cbposition=fig_2.add_axes([0.1,0.05,0.9,0.07])
    
29. fig_2.colorbar(c2,cax=cbposition,orientation='horizontal',format='%.1f')

复制代码

Jareau

图是这样的

Jareau

Jareau 发表于 2021-5-14 13:01
图是这样的

自行解决了，把中心纬度设置在0度就可以了

hypon

马克Mark一下

AFFan

请问我想mark中国区，但是画图的时候只画[37-47N,75-125E]范围，然后画出来就会是这样，应该怎么改比较好阿T T

代码就只用了一行
clip = shp2clip(plot2,ax2,r'e:\\python\\package\\mask\\country1.shp','China')
C:\Users\NEU\Desktop\20210619214405

AFFan

AFFan 发表于 2021-6-19 21:46
请问我想mark中国区，但是画图的时候只画[37-47N,75-125E]范围，然后画出来就会是这样，应该怎么改比较好阿 ...

害 计算的时候把范围限制在要画图的范围内就好了，感恩

空城不空

Jareau 发表于 2021-5-21 21:49
自行解决了，把中心纬度设置在0度就可以了

请问是改这里吗“proj = ccrs.PlateCarree(central_longitude=115)”，就把115改成0就可以了是吗