如何解决无法使用 cartopy 显示纬度和经度或海岸线
我正在使用 netcdf 文件中的 GOES-16 卫星图像。我想裁剪文件,使其只显示波多黎各和轮廓。我可以让它裁剪,但纬度和经度不会显示,波多黎各的海岸线也不会显示。下面是我的示例代码。
import netCDF4 as nc
from netCDF4 import Dataset
import cartopy as ccrs
import cartopy.mpl.ticker as cticker
import matplotlib
import matplotlib.pyplot as plt
import matplotlib.ticker as mticker
#from mpl_toolkits.basemap import Basemap
import metpy
import numpy as np
import cv2
from cv2 import remap
import xarray
fn = 'OR_ABI-L2-ACMC-M6_G16_s20211730001172_e20211730003545_c20211730004314.nc'
ds = xarray.open_dataset(fn)
nc = nc.Dataset(fn)
variable = 'BCM'
bcm = ds['BCM'][:].data
bcm = np.clip(bcm,1)
geo_extent = nc.variables['geospatial_lat_lon_extent']
print(geo_extent)
#set latitude and longitude values from Metadata
min_lon = float(geo_extent.geospatial_westbound_longitude)
max_lon = float(geo_extent.geospatial_eastbound_longitude)
min_lat = float(geo_extent.geospatial_southbound_latitude)
max_lat = float(geo_extent.geospatial_northbound_latitude)
dat = ds.metpy.parse_cf('BCM')
geos = dat.metpy.cartopy_crs
x = dat.x
y = dat.y
print(x.min())
fig = plt.figure(figsize=(12,15))
ax = fig.add_subplot(1,1,projection=geos)
ax.set_extent([-68.07,-64.65,16.7,19.79],crs=geos)
ax.imshow(bcm,origin='upper',extent=(min_lon,max_lon,min_lat,max_lat),transform=geos)
ax.coastlines(resolution='10m',color='black',linewidth=1)
ax.add_feature(ccrs.cartopy.feature.COASTLINE)
plt.title('Puerto Rico',loc='left',fontweight='bold',fontsize=15)
plt.show()
这是地图的图片。 [1]:https://i.stack.imgur.com/Y0qbY.png
解决方法
我想通了。我使用了错误的语法。
dat = ds.metpy.parse_cf('BCM')
geos = dat.metpy.cartopy_crs
x = dat.x
y = dat.y
pc = ccrs.crs.PlateCarree()
fig = plt.figure(figsize=(12,15))
ax = fig.add_subplot(1,1,projection = pc)
ax.set_extent([-68.07,-64.65,16.7,19.79],crs = pc)
ax.imshow(bcm,origin='upper',extent=(x.min(),x.max(),y.min(),y.max()),transform=geos,interpolation = 'none')
ax.coastlines(resolution='10m',color='black',linewidth=1)
ax.add_feature(ccrs.cartopy.feature.COASTLINE)
ax.set_xticks((-68.07,-64.65),minor = 'True',crs= pc)
ax.set_yticks((16.7,19.79),minor= 'True',crs= pc)
plt.title('Puerto Rico Clear-Sky Mask',loc='center',fontweight='bold',fontsize=15)
plt.show()
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。