MODIS积雪面积产品介绍、下载及matlab读取~~

逸儿~~

        
        MODIS（MODerate-resolution Imaging Spectroradiometer）是Terra和Aqua卫星上搭载的主要传感器之一，它的通道范围很广，两颗卫星相互配合每1-2天可重复观测整个地球表面，获取一次全球观测数据。MODIS标准数据产品包括：陆地标准数据产品、大气标准数据产品和海洋标准数据产品等三种主要标准数据产品类型，其中MOD10为陆地2、3级标准数据产品，内容为雪覆盖，有逐日和合成的积雪资料两种。MOD10A1为逐日积雪产品，MOD10A2为MOD10A1的8日合成的积雪产品，分辨率均为500 m，投影方式为等面正弦投影，其中MOD10A2主要用于高纬地区，为了避免每日积雪产品中云层的影响。MOD10C1也为逐日积雪产品，分辨率为0.05°*0.05°，投影方式为经纬度投影，其8日合成积雪产品为MOD10C2。下面主要介绍MOD10C1~~

-------------------------------------------
   MOD10C1是格点资料，不用转换投影方式，其主要数据项包括：

DAY_CMG_Snow_Cover:是每天全球的雪覆盖范围，雪覆盖范围的观察值是以占整个格网的百分比表示的，有效范围是0-100%；

DAY_CMG_Confidence_Index:记录了陆地（Land）在观察值中的百分比。百分比越高，雪范围（Snow extent）的可信度越大，云的遮挡能够减少可信指数；

DAY_CMG_Cloud_Obscured:表明每天里有多少陆地land一直被云所遮挡；

Snow_Spatial_QA:表明整个格网中的质量评价。

----------------------------------------------------------------------

MOD10C1资料下载：

首先登陆数据下载网址： http://reverb.echo.nasa.gov/reve ... C61468236-NSIDC_ECS   

然后输入数据产品名称、选择起始时间，搜索：

记住一定要加入“小车”，然后去“小车”里查看~~

然后点击“Order~~

点击“submit Order”~

填写联系信息，会收到邮件~

邮件中有ftp地址，将ftp地址（红色框框）复制到FTP下载工具（本帖中用的是CuteFTP,这个是收费的，推荐filezilla，免费的，也很好用）中的Host address,下图

用户名是Anonymous,密码是邮箱地址，然后Connect~~

选择类型，然后把要下的hdf文件用shift+鼠标选中，拖到左边电脑里的某个文件夹，然后当所有的都finish时，就搞定啦！！！！！！！

---------------------------------------------------

现在来看数据提取咯~

我下了三年的全球数据，从中提取九个格点每年10月1日至次年5月31日的资料：

原始数据名称是这样的：

为了便于循环读取hdf文件，将文件名进行了批量修改，批量修改文件名可参考网址：http://wenku.baidu.com/view/a2db842327284b73f3425054.html  非常简单快捷神奇！！文件名修改后如下：

文件名含义：如A2011274，代表2011年第274天~~

-------------------------------------------------------------------------

matlab循环读取hdf文件的程序如下：

S='A2011';%%读取2011年10月1日至12月31日9个格点的数据~~

l=0;

s11=999*ones(92,9);%%%2011年10.1-12.31是92天

for day=274:365

    D=int2str(day);

    l=l+1;

    file=strcat(S,D,'.hdf')

    %Day_CMG_Snow_Cover = hdfread(file, '/MOD_CMG_Snow_5km/Data Fields/Day_CMG_Snow_Cover', 'Index', {[1  1],[1  1],[3600  7200]});

    %Day_CMG_Confidence_Index = hdfread(file, '/MOD_CMG_Snow_5km/Data Fields/Day_CMG_Confidence_Index', 'Index', {[1  1],[1  1],[3600  7200]});

    %Day_CMG_Cloud_Obscured = hdfread(file, '/MOD_CMG_Snow_5km/Data Fields/Day_CMG_Cloud_Obscured', 'Index', {[1  1],[1  1],[3600  7200]});

    Snow_Spatial_QA = hdfread(file, '/MOD_CMG_Snow_5km/Data Fields/Snow_Spatial_QA', 'Index', {[1  1],[1  1],[3600  7200]});

    for i=698:700   %%i为提取格点698/699/700

     m=i-697;

       for j=5550:5552   %%j为提取格点5550/5551/5552

           k=j-5549;

           z=(m-1)*3+k;

           %s11(l,z)=Day_CMG_Snow_Cover(i,j);

           %s11(l,z)=Day_CMG_Confidence_Index(i,j);

           %s11(l,z)=Day_CMG_Cloud_Obscured(i,j);

           s11(l,z)= Snow_Spatial_QA (i,j);

       end

   end

end

-------------------------------------------------------

搞定！！撒花~~~~~~~现在开始做人生中最幸福的事儿之一 看着动漫吃西瓜~~~~~~~~~~~

mxgang

66666   可以学习到很多东西啊！！！赞一个

xiao豆豆

{:5_235:}{:5_235:}必须要赞一个，好东西呀！！！！

逸儿~~

xiao豆豆 发表于 2016-5-12 21:51
必须要赞一个，好东西呀！！！！

谢谢谢谢~{:eb500:}

逸儿~~

嘿嘿，谢谢~{:eb511:}

贫道敬孔

多谢楼主分享{:5_257:}

落木萧萧ice

谢谢分享！

bwbj

楼主这里的简直太全面了，赞一个

曦尘

谢谢分享

aas314

感谢楼主  赞