注意:对于下载的AMSR-E数据是hdf格式,需要转换投影即可进行数据解码。
-------------------------------------------------------------------------------------------------------------------------------------------------------------
1 通过美国NSIDC的Data poolftp(ftp://n4ftl01u.ecs.nasa.gov/SAN),下载最近几十天的数据(只有些L2A的数据),好像不是很系统;或者在Preliminary Data注册一下(只包含最近7天的L1A数据),通过
Host name: sidads.colorado.edu
User ID: amsrenrt
Password: Titap$wd
下 载。这两个ftp都相当难登陆。而且,不能用代理。用代理是登陆不了的!
2 通过EOS data Gateway订购。
网址是:http://delenn.gsfc.nasa.gov/~imswww/pub/imswelcome/index.html
注 册一个帐户,然后选数据,订数据。通过邮寄光盘或者通过ftppull(从美国的ftp下载)or ftppush(从本地的ftp上下载)。基本上除了L1B的数据没有,其它的产品都有。
3 通过JAXA的数据分发中心下载。
个人认为这个系统比EOS的系统好用。可能EOS的还没用熟。
网址:https://www.eoc.jaxa.jp/iss/jsp/indexEn.html
用 代理也可以下载。注册了账号,可以检索,下载任何级别的产品(权限比principle Investigator可能要小些)。订了产品,过一会儿就可以下载了(会发邮件通知),只是bnu的网络下速度只有3.0 k/s左右。
-------------------------------------------------------------------------------------------------------------------------------------------------------------
最近有朋友问起关于AMSR-E Tb数据的处理方法,现总结如下:
1、数据介绍
http://nsidc.org/data/docs/daac/nsidc0301_amsre_gridded_tb.gd.html
2、数据下载
ftp://sidads.colorado.edu/pub/DATASETS/brightness-temperatures/
3、数据预处理
1)在Arc GIS 中处理
http://nsidc.org/data/ease/gis.html
2)在ENVI中处理
http://nsidc.org/data/ease/geolocate/index.html
3)ArcGIS 9 中批处理
1、解压后的的文件在DOS下提取文件名,并加扩展名为.bsq
DIR /W .> name.txt
RENAME ID2.03-AMSRE-B01-NL2002170A.06H NL02170A06H.BSQ
在excel中可以使用concatenate函数为新的文件名加后缀
2、将所下的头文件COPY同样数目个,并改名称为和文件名称一致的.hdr文件,并放在和文件相同的目录下
COPY NL.hdr EASE-F13-NL2000064D37H.hdr
3、利用ARCMap下的ArcToolBox下的Raster To Other Format工具将BSQ文件批量撰为GRID
或在Acr Infor 工作站下执行:
imagegrid NL02365D23H.bsq NL02365D23H nearest
4、利用ARCMap下的ArcToolBox下的Samples->Data Management->Projections->Batch Define Coordinate System工具批量定义投影为EastGrid
或在Acr Infor 工作站下执行:
projectcopy file EASE-Grid_north.prj grid NL02365D23H
可利用以下信息建立prj文件
PROJCS["User_Defined_Lambert_Azimuthal_Equal_Area",GEOGCS["GCS_User_Defined",DATUM["D_User_Defined",SPHEROID["User_Defined_Spheroid",6371228.0,0.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Latitude_Of_Origin",90.0],UNIT["Meter",1.0]]
5、利用ARCINFOR 工作站转投影
project grid kkk1 kkkk2
output
projection albers
units meters
spheroid krasovsky
parameters
25 0 0
47 0 0
105 0 0
0 0 0
0.0
0.0
end
6、重采样到25000
AL02_363A_06H = resample ( AL02363A_06H , 25000 , nearest )