NCL读取FY-4A LMIE资料尝试

最近对FY4A的闪电资料很感兴趣，看了小流星大大的这篇用python读取闪电数据的博客后，想用NCL来读取（毕竟还在NCL入门中……），但是后面的处理可能不如python方便，鉴于本人的python水平，还是先把自己熟悉的搞定吧…… 一次读取FY4A雷电数据（LMI）的过程总结_野生的气象小流星的博客-CSDN博客https://blog.csdn.net/weixin_42372313/article/details/114694684?spm=1001.2014.3001.5502

废话不多说，对照着小流星大大的文章，开始读数据吧~

1、资料获取

 之前在气象数据网上下载过FY4A的资料，但是最多只能下载最近3个月的，也就是21年冬季的数据，但是我想研究夏季的闪电，所以找某网友要了，其他还有什么可找到历史资料的途径吗？

2、数据结构

由官方说明得到的信息：

 1）此资料为“事件”，是由一个区域内的闪电聚类而来，推测区域大小与FY4A的分辨率有关（7.8km），这种处理方式可能与LIM是光学仪器的特性有关，但是具体的原理还没有搞明白

2）NC文件

3）变量名及属性（头文件）：主要关注闪电发生的位置（lat,lon），event radiance（ER）和footprint（EFP）

3、简单画一下闪电事件发生的位置：（参考官网scatter plot脚本）

​
begin
;
f = addfile("FY4A-_LMI---xxxxxxxxxxxxxxxxxxxxxx.NC","r")
lon = f->LON 
lat = f->LAT 
wks = gsn_open_wks("png", "FY4A-LIME-map")

;plot settings
res = True 
;res@gsnAddCyclic = False
res@gsnFrame    = False                             ; don't advance yet
res@gsnMaximize = True
res@pmTickMarkDisplayMode="Always"   ;add the marker of angle 
res@tiMainString = "Event Locations"
res@mpMinLatF = min(lat) - 5
res@mpMaxLatF = max(lat) + 5
res@mpMinLonF = min(lon) - 5
res@mpMaxLonF = max(lon) + 5
res@mpDataBaseVersion = "Ncarg4_1"
res@mpDataSetName="Earth..4"
res@mpOutlineOn=True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)
res@mpGeophysicalLineThicknessF=2.0 
res@mpNationalLineThicknessF=2.0

 map = gsn_csm_map(wks,res)                ; Create the map plot, don't draw
 
 mkres               = True
 mkres@gsMarkerIndex = 2
 mkres@gsMarkerColor = "Brown"
 gsn_add_polymarker(wks,map,lon,lat,mkres)

  draw(map)
  frame(wks)

end

​

 感觉图片不是特别清楚hh

@野生的气象小流星  最后感谢小流星大大的教程，对数据处理小白很有帮助！