MeteoInfoLab脚本示例:AMSR-E卫星数据投影

本文演示如何读取AMSR-E Land3数据,将其从Cylindrical_Equal_Area投影转换为等经纬度投影,并展示数据可视化。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

AMSR-E(http://nsidc.org/data/amsre/index.html)数据中的Land3数据是HDF-EOS4格式,投影是Cylindrical_Equal_Area。这里示例读取数据并投影至等经纬度投影。

脚本程序:

#Add data file
fn = 'AMSR_E_L3_DailyLand_V06_20091231.hdf'
f = addfile(os.path.join('D:/Temp/hdf', fn))
#vname = 'D_Soil_Moisture'
vname = 'A_TB36.5H_(Res_1)'
data = f[vname][::-1,:]
data[data==0] = -9999
data.fill_value = -9999
data = data * 0.1
yn = data.dimlen(0)
xn = data.dimlen(1)
#Project data
toproj = projinfo()    #longlat projection
x = linspace(-180, 180, xn)
y = linspace(-90, 90, yn)
lon, lat = meshgrid(x, y)
data = data.project(lon, lat, toproj=toproj)
#Plot
axesm(tickfontsize=12)
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer, edgecolor='k')
levs = arange(140, 320, 10)
layer = imshowm(x, y, data, levs)
colorbar(layer, fontsize=12)
xticks(arange(-180, 181, 30))
yticks(arange(-90, 91, 30))
title([fn, vname])

转载于:https://www.cnblogs.com/yaqiang/p/5096002.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值