MeteoInfoLab脚本示例：CAPLIPSO数据

citrushe

# Add file
fn = 'C:/Users/lenovo/Desktop/data/CAL_LID_L2_333mCLay-ValStage1-V3-30.2014-10-14T04-13-23ZD.hdf'
f = addfile(fn)
# Read data
vname = 'Feature_Classification_Flags'
var = f[vname]
data = var[:,:]
lat = f['latitude'][:,0]
# Extract Feature Type only through bitmask.
data = data & 7
# Subset latitude values for the region of interest (40N to 62N).
lat = lat[3500:4000]
size = lat.shape[0]
data2d = data[3500:4000, 1165:]  # -0.5km to  8.2km
data3d = reshape(data2d, (size, 15, 290))
data = data3d[:,0,:]
# Focus on cloud (=2) data only.
data[data > 2] = 0
data[data < 2] = 0
data[data == 2] = 1
# Generate altitude data according to file specification [1].
alt = zeros(290)
# -0.5km to 8.2km
for j in range (0, 290):
    alt[j] = -0.5 + j*0.03
# Plot
levs = arange(2)
cols = ['w','b']
ls = makesymbolspec('image', levels=levs, colors=cols)
layer = imshow(lat, alt, rot90(data, 3), symbolspec=ls)
colorbar(layer, ticks=['Others','Cloud'])
basename = os.path.basename(fn)
title([basename, 'Feature Type (Bits 1-3) in Feature Classification Flag'])
xlabel('Latitude (degrees north)')
ylabel('Altitude (km)')

MeteoInfo

citrushe 发表于 2016-9-2 16:55
# Add file
fn = 'C:/Users/lenovo/Desktop/data/CAL_LID_L2_333mCLay-ValStage1-V3-30.2014-10-14T04-13- ...

错误提示第14行出错。我们用的数据文件不一样，数据的维长度可能也不同，你看看data的各维长度是多少？运行脚本后，在Console里输入：

>>> data.dims

citrushe

>>> data.dims
[Name: fakeDim92
Min value: 0.0
Max value: 63599.0
Size: 63600
Delta: 1.0, Name: fakeDim93
Min value: 0.0
Max value: 4.0
Size: 5
Delta: 1.0]

显示是这样的

MeteoInfo

citrushe 发表于 2016-9-2 17:28
>>> data.dims
[Name: fakeDim92
Min value: 0.0

数据维数明显是和示例文件不一样的。你需要根据自己数据的情况修改相应代码。

citrushe

MeteoInfo 发表于 2016-9-2 18:39
数据维数明显是和示例文件不一样的。你需要根据自己数据的情况修改相应代码。

好的~多谢王老师

cashia371081

王老师，学习了您的帖子。
希望老师能给出例子中的数据文件。希望增加一个功能：在editor栏选中某行命令时，右键可以增加运行这一行命令的功能。这样就不用copy到console栏去了。
另外MeteoInfoLab现在有没有命令手册。变量查看那一栏现在是不是不能用？data = data & 7 这句命令是什么没看懂。

MeteoInfo

cashia371081 发表于 2016-11-21 09:38
王老师，学习了您的帖子。
希望老师能给出例子中的数据文件。希望增加一个功能：在editor栏选中某行命令时 ...

相关命令可以在这里查阅：http://www.meteothinker.com/docs/meteoinfolab/functions.html，但是不全。也可以在console中输入help(function name)看相关命令的帮助文档。

右键运行选中行命令的功能在以后的版本中会考虑的。

&是位运算的一种，具体你可以在网上搜索一下。

2020Doraemon

sunflower$ 发表于 2016-1-21 10:47
官网有直接给出matlab和IDL 的程序 稍微改改就能用

请问你知道怎么改吗？

突破的渴望

王老师您好，您知道meteoinfo 怎么做calipso某个轨迹点的色比和退偏振比廓线图吧？谢谢您~~

wuyu

老师好厉害，受用匪浅