求助画图

易烊启迪

大佬们好，求助。我做这个投影图遇到一个问题。就是我这个图的数据经度是0-357.5（间隔2.5），然后纬度是0-90。由于经度不是到360这个数字，所以图形上空白了一部分，缺了这一小部分。就是想请教您们，如何把经度为0的值赋给经度为360的，怎么利用代码把这一部分添加上。数据文件有点大，上传不了。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
def make_map(ax):
    projection = ccrs.Orthographic(central_longitude=-90,central_latitude=90)
    ax.set_global()
    ax.coastlines(linewidth=0.5)
    '''标注坐标轴'''
    '''zero_direction_label=True 有度的标识，False则去掉'''
    lon_formatter = LongitudeFormatter(zero_direction_label=True)
    lat_formatter = LatitudeFormatter()
    '''添加网格线'''
    return ax
data=xr.open_dataset(r'D:\o3 model data\tggo3.cam2.h1.2019-11-01-00000.nc')
time=data['time'][1]
lev=data['lev'][51]
lat = data['lat'][:]
lon = data['lon'][:]
T=data['T'].loc[time,lev,:,:]
data_proj = ccrs.PlateCarree(central_longitude=0)
fig = plt.figure(dpi=1200)
proj = ccrs.Orthographic(central_longitude=0,central_latitude=90)
ax = fig.add_subplot(projection = proj)
ax = make_map(ax)
cf = ax.contourf(lon, lat, T, cmap='Spectral_r',transform=data_proj)
plt.show()

易烊启迪

图是这样的，下方空白了一小部分

                               
登录/注册后可看大图

来自苹果APP客户端

SolaProject

试试这个？
from cartopy.util import add_cyclic_point
T, lon = add_cyclic_point(T, lon)

你是最美七月天

虽然是全球数据，但是采用极射赤面投影在绘制过程中需要插值到0度线，需要采用add_cyclic_point函数，楼上是正解

易烊启迪

SolaProject 发表于 2022-08-21 20:06
试试这个？
from cartopy.util import add_cyclic_point
T, lon = add_cyclic_point(T, lon)

感谢，

                               
登录/注册后可看大图

来自苹果APP客户端

易烊启迪

你是最美七月天 发表于 2022-08-21 22:40
虽然是全球数据，但是采用极射赤面投影在绘制过程中需要插值到0度线，需要采用add_cyclic_point函数，楼上是正解

感谢

                               
登录/注册后可看大图

来自苹果APP客户端

易烊启迪

SolaProject 发表于 2022-08-21 20:06
试试这个？
from cartopy.util import add_cyclic_point
T, lon = add_cyclic_point(T, lon)

请问您知道如果想画40-90纬度的，经度不变0-360，这样区域的地图，可以怎么修改代码吗？

来自苹果APP客户端

SolaProject

易烊启迪 发表于 2022-8-22 20:31
请问您知道如果想画40-90纬度的，经度不变0-360，这样区域的地图，可以怎么修改代码吗？

我也不太了解，不咋画这种图，不过用这个投影应该差不多，算了下差不多这个距离：
projection = ccrs.NearsidePerspective(central_longitude=0.0, central_latitude=90, satellite_height=1945746, false_easting=0)

易烊启迪

SolaProject 发表于 2022-08-24 00:18
我也不太了解，不咋画这种图，不过用这个投影应该差不多，算了下差不多这个距离：
projection = ccrs.NearsidePerspective(central_longitude=0.0, central_latitude=90, satellite_height=1945746, false_easting=0)

好的感谢

                               
登录/注册后可看大图

来自苹果APP客户端