modis数据_气象编程 | 如何用NCL处理风云4A/MODIS卫星数据？

 欢迎订阅微信公众号:『气象学家』

阅读建议

读取和处理了两种FY-4A和MODIS卫星数据，进行相关产品的绘图，插值为不同分辨率经纬度格点数据并保存为nc格式文件。抛砖引玉，不做更深入的分析，有任何错误欢迎大家批评指正。后文附数据脚本获取方式。

内容目录

00.前言介绍

01.中国区域FY4_AGRI_L2 OLR原始数据Lambert投影绘图和脚本

02.中国区域FY4_AGRI_L2 OLR原始数据经纬度网格投影绘图和脚本

03.中国区域FY4_AGRI_L2 OLR插值数据经纬度网格投影绘图和脚本

04.东亚区域MOD06_L2 Cloud_Top_Temperature原始数据绘图和脚本

05.东亚区域MOD06_L2 Cloud_Top_Temperature插值数据绘图和脚本

06.数据和脚本获取

07.参考

00.前言介绍

工具：NCL、相关地图底图包配料：风云4A数据FY4_AGRI_L2 OLR产品、MOD06_L2产品方法：ESMF_regrid等成品：中国/东亚区域图、经纬度格点数据(e.g. ,  0.1°*0.1°)

01.中国区域FY4_AGRI_L2 OLR原始数据Lambert投影绘图和脚本

FY4A_read_plot_latlon2d_Lambert.ncl

02.中国区域FY4_AGRI_L2 OLR原始数据经纬度网格投影绘图和脚本

; =============================================================================
; Email : Zhpfu.atm@gmail.com

03.中国区域FY4_AGRI_L2 OLR插值数据经纬度网格投影绘图和脚本

FY4A_read_plot_latlon2d_to_rectilinear_grid.ncl

04.东亚区域MOD06_L2 Cloud_Top_Temperature原始数据绘图和脚本

MODIS_Full_Area_read_plot_latlon2d_to_rectilinear_grid.ncl

05.东亚区域MOD06_L2 Cloud_Top_Temperature插值数据绘图和脚本

; =============================================================================
; Author: Gavin | Affiliation: NJU
; Email : Zhpfu.atm@gmail.com
; Last modified:    2019-07-19 01:15
; Filename:        MODIS_Full_Area_regrid_latlon2d_to_rectilinear_grid.ncl
; Description: Transform xy kilometer to lat lon grid;
;              Output data with NetCDF; 
; To Be Determined(TBD):  
;               
; =============================================================================
procedure print_elapsed_time(start_time,title)
local end_time
begin

  end_time = get_cpu_time()
  print("======================================================================")
  print(title + " elapsed time = " + (end_time-start_time) + " CPU seconds.")
  print("======================================================================")
end

; =============================================================================
; ===================================MAIN======================================
; =============================================================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

  start_time = get_cpu_time()

  fpath = "/Users/zhpfu/Downloads/DATA_MAC/MODIS/Full_Area/"
  fname = "MOD06_L2.