提取FY -4A干涉式大气垂直探测仪GIIRS数据

本文主要是介绍提取FY -4A干涉式大气垂直探测仪GIIRS数据，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

提取FY -4A干涉式大气垂直探测仪GIIRS数据

第一步，注册账号并通过审核

1、在国家气象科学数据中心注册账号；http://data.cma.cn/user/toLogin.html
2、用户类型选择“个人实名注册用户”（或者条件符合，科研实名注册用户），经过1-2个工作日审核；
实名认证

第二步，下载FY-4A GIIRS数据

1、登录风云卫星遥感数据服务网
http://satellite.nsmc.org.cn/PortalSite/Data/Satellite.aspx
2、用户名和密码与“国家气象科学数据中心”注册账号一致；
3、选择FY-4A GIIRS数据，并进行检索；
根据需求进行数据检索
4、进行文件下载，可选中多个至购物车下载，也可单点下载。如果需求文件量不是很大，推荐单点下载。
下载HDF数据

第三步，对HDF文件的读取

可通过HDFView打开也可通过Matlab代码进行数据提取
1、HDFView软件打开此前下载的HDF文件，可以在“Object Attribute Info”以及“General Object Info”面板中了解一些参数信息，如通道数，光谱分辨率等等
HDF文件概述

2、可通过HDFView将数据导出
在这里插入图片描述
3、可通过Matlab的代码读取数据h5read(file_dir,’/ES_NEdRLW’)’;
下面展示一些 内联代码片。

%% 用于提取FY-4A GIIRS的中长波信息
clear
clc
bs1_col = (1650:0.625:2250)';                %中波波数
bs2_col = (700:0.625:1130)';                 %长波波数
wavelen1 = 10000./bs1_col;                   %中波波长
wavelen2 = 10000./bs2_col;                   %长波波长
para1=1E-11;                                 %波数到波长数据的转换的系数
file_dir = 'path\filename';
ES_NEdRLW=h5read(file_dir,'/ES_NEdRLW')';
ES_NEdRMW=h5read(file_dir,'/ES_NEdRMW')';
ES_RealLW=h5read(file_dir,'/ES_RealLW')';
ES_RealMW=h5read(file_dir,'/ES_RealMW')';data1 = ES_NEdRLW.*bs2_col.*bs2_col.*para1;
data2 = ES_NEdRMW.*bs1_col.*bs1_col.*para1;
data3 = ES_RealLW.*bs2_col.*bs2_col.*para1;
data4 = ES_RealMW.*bs1_col.*bs1_col.*para1;

获取实时的风云四号卫星遥感图像

1、进入风云四号卫星天气应用平台（http://rsapp.nsmc.org.cn/geofy/?i=0&isPlay=false&speed=2&sat=fy-4a&pro=geos&type=full_disk&band=28&overlay=&x=-252219.5552927202&y=562329.3184740966&z=1&area=1&ll=0&county=1&duration=36&interval=1&c=false&cp=0.5&st=&et=&ac=&hide=0&s=1）
2、根据需求选择绘制数据
云类型