求助 python极射赤面投影

易烊启迪

大佬们好，我想画一个纬度是40-90，经度0-360的地图，然后我不知道怎么调纬度范围，下面这个图纬度是0-90的北半球图。

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
from netCDF4 import Dataset
import numpy as np
import xarray as xr
from cartopy.util import add_cyclic_point
plt.rcParams['font.sans-serif']=['Arial']#显示中文
plt.rcParams['axes.unicode_minus']=False
#读取数据
data=xr.open_dataset(r'Z:\shuju\o3\nc\OMI0401.nc')
lat = data['lat'].loc[:]
lon = data['lon'][:]
o3=data['ColumnAmountO3'].loc[:,:]
#======绘图=======
fig = plt.figure(dpi=1200)
leftlon, rightlon, lowerlat, upperlat = (0, 360, 0, 90)#这里修改纬度0变成40，图形不变
proj = ccrs.Orthographic(central_longitude=0, central_latitude=90)
data_proj = ccrs.PlateCarree(central_longitude=0)
ax = fig.add_subplot(projection=proj)
ax.set_global()
ax.add_feature(cfeature.LAKES, alpha=0.5) #添加湖泊
ax.coastlines(linewidth=0.5) #添加海岸线
#ax.gridlines(draw_labels=True, x_inline=False, y_inline=False)#去掉地形
lon_formatter = LongitudeFormatter(zero_direction_label=True)
lat_formatter = LatitudeFormatter()
ax.xaxis.set_major_formatter(lon_formatter)
ax.yaxis.set_major_formatter(lat_formatter)
ax.set_title('2020.04.01 OMI',x=1,y=1,fontsize=10)
ax.gridlines(color='black', linestyle=':', xlocs = np.arange(-180,180,45), ylocs = np.arange(0,90,15),draw_labels=True)
cf = ax.contourf(lon, lat, o3, cmap='Spectral_r',levels=np.arange(180,525,15),transform=data_proj)
plt.colorbar(cf, pad=0.042, orientation='horizontal', shrink=0.87, aspect=35)
plt.show()

                               
登录/注册后可看大图

来自苹果APP客户端

名字真难取

你这不是ax.set_global了吗,而且你就设了extent范围，也没写set_extent啊

易烊启迪

名字真难取 发表于 2022-08-23 09:33
你这不是ax.set_global了吗,而且你就设了extent范围，也没写set_extent啊

您好，我应该怎么修改，不太懂

                               
登录/注册后可看大图

来自苹果APP客户端