抱歉,您的浏览器无法访问本站

本页面需要浏览器支持(启用)JavaScript


了解详情 >

前言

在网上新寻找到的一款工具,目前绘制FY-4A/FY-3D/Sentinel3/
参考朝曦dawn & 官方Doc

  • 支持多种合成产品:雾;沙尘;真彩色等

Install

为了确保安装正确,请使用conda; 下载 Miniconda

$ conda create -c conda-forge -n my_satpy_env python satpy
$ conda activate my_satpy_env

Start

绘制全圆盘真彩色图像。准备一份FY-4A的4KM全圆盘L1数据。

import os, glob
from satpy.scene import Scene
filenames = glob.glob('C://Users/Moe/Desktop/FY4A-_AGRI*4000M_V0001.HDF')
scn = Scene(filenames, reader='agri_l1')
composite = 'true_color'
scn.load([composite])
scn.save_dataset(composite, \filename='{sensor}_{name}.png')

单通道绘制示例

使用Satpy自带的scn.crop框定区域,与使用投影相比使用范围会受到限。

channel = "C02"
filenames = glob.glob('C:/Users/Administrator/Desktop/FY4A*')
scn = Scene(filenames, reader='agri_l1')
scn.load([channel])
lonmin, latmin, lonmax, latmax = 115, 32, 130, 47
scn_c = scn.crop(ll_bbox=(lonmin, latmin, lonmax, latmax))
sat_w_tb_vis = scn_c[Channel].values
area_def = scn_c[Channel].attrs['area']
crs = area_def.to_cartopy_crs()
fig, ax = plt.subplots(subplot_kw=dict(projection=ccrs.PlateCarree()))
ax.set_extent([lonmin, lonmax, latmin, latmax], crs=ccrs.PlateCarree())
img = plt.imshow(sat_w_tb_vis, origin='upper', cmap=plt.cm.binary_r, extent=crs.bounds, transform=crs)
plt.savefig('4A_AGRI.png', dpi=300, bbox_inches='tight')

拓展阅读

评论