风云4卫星L1数据HDF文件行列如何转化为经纬度？

menseye

各位大神，小白请教个问题：
我利用网上找的代码画出来风云4A的真彩图，但是怎么才能将行列转化为经纬度并且投影呢？
代码附上：
from netCDF4 import Dataset
import matplotlib.pyplot as plt
from skimage import io, data, img_as_float, img_as_ubyte, img_as_uint, img_as_int, exposure
import cv2

hdf_data_path = "D:/Python_total/DATA\FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20210720075336_20210720075753_4000M_V0001.HDF"

nc_obj = Dataset(hdf_data_path)
#print(nc_obj)
#type = nc_obj.variables
#print(type)


B = nc_obj.variables['NOMChannel01'][:]
G = nc_obj.variables['NOMChannel02'][:]
R = nc_obj.variables['NOMChannel03'][:]


# 将数据类型转成int16
B = img_as_int(B)
G = img_as_int(G)
R = img_as_int(R)

# opencv 将rgb三个通道融合
img3 = cv2.merge([R, G, B])

img3 = exposure.adjust_log(img3, inv=True)#调整对比度
img3 = exposure.adjust_gamma(img3, 1.2)  # 图像调暗
T = 2
for i in range(len(img3)):
    for j in range(len(img3[0])):
         r = img3[j][0]
         g = img3[j][1]
         b = img3[j][2]
         img3[j] = (r * 0.6, g, b)
         if r / g > T:
              img3[j] = (g, r * 0.7, b)
#选取特定行列范围的图像
plt.xlim(1000,2747)
plt.ylim(1071,0)

plt.imshow(img3, )
fig =plt.figure(figsize=(15,10))
plt.axis('off')
plt.show()

坎坷

之前有直接转化为经纬度的，就相当于是格点数据，但即使用公式算还是太慢了，就用你这个直接画RGB的方法最快，再琢磨一下框定范围就行了

Lucian

我之前做过风云4号卫星云图的可视化，看我的发的最近一个帖子，好像也是你搜到的类似的代码

这个代码应该有提到一些方法，比如查表，用公式计算之类的，你再多搜索搜索

edwardli

【气Py-47-卫星-哔哩哔哩】 https://b23.tv/LAN7aVx

edwardli

http://bbs.06climate.com/forum.php?mod=viewthread&tid=101507&extra=page%3D2&mobile=no

menseye

Lucian 发表于 2022-4-30 07:15
我之前做过风云4号卫星云图的可视化，看我的发的最近一个帖子，好像也是你搜到的类似的代码

这个代码应 ...

好的，感谢，我再找找

menseye

edwardli 发表于 2022-4-30 10:09
【气Py-47-卫星-哔哩哔哩】 https://b23.tv/LAN7aVx

感谢感谢，正好我最近也在看这个视频

menseye

edwardli 发表于 2022-4-30 10:11
http://bbs.06climate.com/forum.php?mod=viewthread&tid=101507&extra=page%3D2&mobile=no

原来您就是李老师，您做的视频讲的很好

menseye

坎坷 发表于 2022-7-3 17:55
之前有直接转化为经纬度的，就相当于是格点数据，但即使用公式算还是太慢了，就用你这个直接画RGB的方法最 ...

感谢大佬指点，我最近搞到了MATLAB画卫星图的脚本，MATLAB的脚本就是您讲的这个思路！