OMI/Aura臭氧数据下载及绘图

下载地址：https://cmr.earthdata.nasa.gov/search/concepts/C1266136071-GES_DISC.html

如果批处理，建议在EarthData上下载，这样子的下载链接可以使用downthemall插件进行下载，速度是极快的，比wget好用太多。

臭氧的单位

国际上习惯将臭氧总量表示为相当于再标准大气压和温度下单位面积上的厚度，通常用多普生单位（DU）作为单位。1DU定义为标准大气压下0.01 mm的厚度，也就是0.001 cm。AERONET与OMI中的臭氧单位都是DU。模式反演SBDART中的臭氧单位是cm。

绘图

绘图工具使用的是M_Map工具箱。由于读取的数据是1440*720的矩阵，为了和常规投影形式一致，这里对齐进行90度旋转。旋转后绘图代码如下：
方法1：

image(data)

方法2：

clear all
close all
clc

% info = h5info('C:\Users\Administrator\Desktop\OMI-Aura_L3-OMTO3e_2004m1001_v003-2012m0409t101417.he5','/HDFEOS/GRIDS/OMI Column Amount O3/Data Fields')
data = h5read('K:\OMI-Aura_L3-OMTO3e\OMI-Aura_L3-OMTO3e_2019m1218_v003-2020m0117t132210.he5','/HDFEOS/GRIDS/OMI Column Amount O3/Data Fields/ColumnAmountO3');
data(data==-1267650600228229401496703205376.000000) = NaN;
data = rot90(data);
image(data)

figure
histogram(data)



%使用m_map绘图
m_res = 0.25;
[m_lon, m_lat] = meshgrid(-180:m_res:180-m_res,90-m_res:-m_res:-90);
figure
clf;
m_proj('miller','lat',82);
m_pcolor(m_lon,m_lat,data);

m_coast('color',[0 0 0], 'linewidth', 1.2);
m_grid('tickdir','out','linewi',2);

colormap(m_colmap('jet','step',10));
colorbar

来自Earthsearch网站的对应数据可视化效果。

自己绘制的效果，说明数据旋转的没有问题。