python里面利用cartopy绘制北极极射投影问题

汐年

下面贴上我的代码，还有在上面添加刻度也是个问题？
import matplotlib.pyplot as plt
import matplotlib as mpl
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
import scipy.io as sio
import numpy as np
import matplotlib.colors as mcolors
import matplotlib.ticker as mticker
import matplotlib.path as mpath
from cartopy.util import add_cyclic_point

#load数据，并且进行数据处理
data=sio.loadmat('F:\python_study\sj\data1013')
print data.keys()
percentage_jja=data['percentage_jja']*100
x=data['x']
y=data['y']
tansform_jja=percentage_jja.transpose()
#cycle_jja,cycle_x=add_cyclic_point(percentage_jja,coord=x)
c_x,c_y=np.meshgrid(x,y)

#绘图
#方法一
data_crs = ccrs.PlateCarree()
plt.figure(figsize=(8,8))
ax = plt.axes(projection=ccrs.NorthPolarStereo())#在这里必须设置相对应的经纬度的中心位置，如果不设置的话是无法出图的！
ax.set_extent([-180, 180, 60, 90], crs=ccrs.PlateCarree())
cont=ax.contourf(c_x,c_y,tansform_jja,transform=data_crs)#这两个命令就可以转换投影方式，对于原始数据来说都是PlataCarree
ax.coastlines()

汐年

https://scitools.org.uk/cartopy/ ... light=understanding，这个网页上有介绍projection和transform的，我是按照上面的来的

对方的凤飞飞

看看这个帖子,其实官网有北极投影,http://bbs.06climate.com/forum.p ... &extra=page%3D1,画上标签可以参考这个帖子,完美解决你的问题,现在,我想请教你一个问题,怎么在南极画上经纬度信息呢,我之前研究过,经纬度标记得不对!!!!!!!!!!!!!!

素素是个小可爱

ax.set_extent([-180, 180, 60, 90], crs=ccrs.PlateCarree())这一句改成ax.set_extent([-180, 180, 60, 90], crs=ccrs.NorthPolarStereo())试试呢

汐年

素素是个小可爱 发表于 2019-12-7 21:16
ax.set_extent([-180, 180, 60, 90], crs=ccrs.PlateCarree())这一句改成ax.set_extent([-180, 180, 60, 90 ...

也不行哎。。。。

汐年

对方的凤飞飞 发表于 2019-12-7 13:20
看看这个帖子,其实官网有北极投影,http://bbs.06climate.com/forum.php?mod=viewthread&tid=33601&extra=pa ...

他这个帖子我之前也研究了好久，但是我不知道哪里的问题，一开始的投影我都画不出来，再往后就更出不来图了。

对方的凤飞飞

https://scitools.org.uk/cartopy/ ... lr-gallery-waves-py