汪帮主卫星遥感空间分享 http://blog.sciencenet.cn/u/zjwang 湖北鄂州人,2011年博士毕业于吉林大学。现从事卫星遥感技术工作。

博文

The function for reading MOPITT CO Level3 HDF data

已有 5912 次阅读 2010-11-14 13:53 |个人分类:计算程序|系统分类:科研笔记

maps/200303/MOP03-Web-20030314-0.png

MOPITT is an instrument flying on NASA's EOS Terra spacecraft, measuring tropospheric carbon monoxide (CO) on the global scale. MOPITT measurements enable scientists to analyze the distribution, transport, sources and sinks of CO, a trace gas produced by methane oxidation, fossil fuel consumption and biomass burning. MOPITT has been operational since March 2000.

The introduction on MOPITT and the figure are from http://www.acd.ucar.edu/mopitt/, and also, you can download the data from this website.

The following function illustrates how to read MOPITT CO Level 3 data by using MATLAB

function mop_co = wzjReadMOPL3(fname)
   hinf = hdfinfo(fname);
   vgMOP03 = hinf.Vgroup(1);
   vgDts = vgMOP03.Vgroup(1);
   lat = hdfread(vgDts.SDS(1));
   lon = hdfread(vgDts.SDS(2));
   pre = hdfread(vgDts.SDS(3));
   mrpd = hdfread(vgDts.SDS(23));
   mrpn = hdfread(vgDts.SDS(24));
   tcd = hdfread(vgDts.SDS(35));
   tcn = hdfread(vgDts.SDS(36));
   mopco.lat = lat;
   mopco.lon = lon;
   [mopco.longrid, mopco.latgrid] = meshgrid(lon, lat);
   mopco.pre = pre;
   mopco.mrpd = mrpd;
   mopco.mrpn = mrpn;
   tcd(tcd<0)=0;
   mopco.tcd = tcd;
   tcn(tcn<0)=0;
   mopco.tcn = tcn;  
   mop_co = mopco;
end

 



https://blog.sciencenet.cn/blog-43777-383638.html

上一篇:How to call SCIATRAN in your program?
下一篇:The function for reading AIRS CO2 Level3 HDF data
收藏 IP: 202.98.17.*| 热度|

0

发表评论 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-12-26 22:54

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部