OMI /オーラオゾンデータのダウンロードと描画

ダウンロードリンク:https//cmr.earthdata.nasa.gov/search/concepts/C1266136071-GES_DISC.html

ここに画像の説明を挿入します

バッチ処理の場合は、EarthDataにダウンロードすることをお勧めします。このようなダウンロードリンクは、downthemallプラグインを使用してダウンロードできます。速度は非常に速く、wgetよりもはるかに使いやすくなっています。
ここに画像の説明を挿入します

オゾンの単位

世界では、オゾンの総量を、標準的な大気圧と温度の下での単位面積あたりの等価厚さとして表すのが通例です。通常、単位はドップソン単位(DU)です。1DUは、標準大気圧での厚さ0.01 mm、つまり0.001cmとして定義されます。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

EarthsearchWebサイトからの対応するデータの視覚化。
ここに画像の説明を挿入します
自己描画の効果は、データのローテーションに問題がないことを示しています。
ここに画像の説明を挿入します

おすすめ

転載: blog.csdn.net/wokaowokaowokao12345/article/details/108790199