FY4A卫星边界是多少经纬度？

洗雨

下载了风云四卫星的中国区域云图，但是看描述文件没发现他的最小经纬度和最大经纬度，有哪位知道吗？

   

MeteoInfo

洗雨 发表于 2020-8-28 16:09
# -*- coding: utf-8 -*-
"""
Created on Thu Aug 27 15:06:01 2020

供参考

洗雨

# -*- coding: utf-8 -*-
"""
Created on Thu Aug 27 15:06:01 2020
绘制风云4可见光2km通道图像
@author: micapsadm
"""

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
import h5py  #导入工具包
import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

f = h5py.File(r'K:\风云\FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20190804232334_20190804232750_2000M_V0001.HDF','r')
print(f.keys())
a = f['NOMChannel06'][:]
print((f['NOMChannel01']).shape)
print(a.shape)
print(f['NOMObsColumn'])
print(a.max(),a.min())


box = [70,100,25,40]

xstep,ystep = 5,5
fig = plt.figure(figsize = (8,8))

#选择绘图区域
a = a[500:2000,2000:4000]
print(a)
crs = ccrs.LambertConformal(central_longitude=91.13,central_latitude=29.67)

ax = fig.add_subplot(1,1,1,projection=crs)
# im = ax.imshow(a,cmap='gist_gray')
#可见光云图，vmax设置511,1023等。红外不设置
im = ax.imshow(a,cmap='gist_gray',vmax=511,vmin=0,origin='upper')

plt.show()

墨家大宝

https://blog.csdn.net/modabao/article/details/83834580

                               
登录/注册后可看大图

洗雨

墨家大宝 发表于 2020-8-28 19:28
https://blog.csdn.net/modabao/article/details/83834580

谢谢老师！

洗雨

MeteoInfo 发表于 2020-8-28 17:02
供参考

谢谢老师，不过我按照您的程序运行时出现了 run script...
Traceback (most recent call last):
  File "<iostream>", line 15, in <module>
AttributeError: 'tuple' object has no attribute 'proj'
这个问题应该怎么解决？

洗雨

MeteoInfo 发表于 2020-8-28 17:02
供参考

谢谢老师，不过我按照您的程序运行时出现了以下问题，您看该怎么调整呢？
run script...
Traceback (most recent call last):
  File "<iostream>", line 15, in <module>
AttributeError: 'tuple' object has no attribute 'proj'

MeteoInfo

洗雨 发表于 2020-8-29 11:52
谢谢老师，不过我按照您的程序运行时出现了以下问题，您看该怎么调整呢？
run script...
Traceback (mo ...

是不是版本的问题？用MeteoInfo最新版本