Download link: https://cmr.earthdata.nasa.gov/search/concepts/C1266136071-GES_DISC.html
If batch processing, it is recommended to download on EarthData, such download links can be downloaded using downthemall plug-in, the speed is extremely fast, much easier to use than wget.
Unit of ozone
It is customary in the world to express the total amount of ozone as the equivalent thickness per unit area under standard atmospheric pressure and temperature, usually in Doppson units (DU) as the unit. 1DU is defined as the thickness of 0.01 mm under standard atmospheric pressure, which is 0.001 cm. The ozone unit in AERONET and OMI is DU. The unit of ozone in model inversion SBDART is cm.
Drawing
The drawing tool uses the M_Map toolbox. Since the data read is a 1440*720 matrix, in order to be consistent with the conventional projection form, the alignment is rotated by 90 degrees here. The drawing code after rotation is as follows:
Method 1:
image(data)
Method 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
Visualization of corresponding data from Earthsearch website.
The effect of self-drawing shows that there is no problem with data rotation.
