IDL遥感应用入门(7):利用QuikScat数据绘制风速图
发布时间: 2014-01-08
所属分类: IDL遥感应用教程
在上面的几篇文章我我写的是关于利用QuikScat数据,绘制风向图。这次重新利用QuikScat数据中的风速数据,来绘制全球范围内的风速图像,并且通过colorbar来显示相关数据。
QuikScat数据说明
QuikScat提供3个方面数据信息:表面风速,表面风向和降雨信息。
其中数据在0到250属于有效信息,大于250的值均属于无效信息。
计算风速的时候需要乘以0.2,风向需要乘以1.5
更多详细信息可以查看这里
详细代码
;Author:Sailor
;2013-11-08
PRO Course_7
;定义文件路径
MyRootDir='D:\3\'
ImgExtName='.jpg'
;遍历文件夹
filearr = file_search(MyRootDir,'*.gz',count=num);
FOR fileindex=0,num-1,1 DO BEGIN
;打开文件
OPENR , LUN , filearr[fileindex] , /COMPRESS,/GET_LUN
;构造显示标题
filename=FILE_BASENAME(filearr[fileindex]);
filename=STRMID(filename,6,10)
year=STRMID(filename,0,4)+'-'
month=STRMID(filename,4,2)+'-'
day=STRMID(filename,6,2)+''
diplayname=year+month+day
;读入数组
FleArr = BYTARR(1440,720,3)
READU , LUN , FleArr
WdSpd = FleArr(*,*,0)*0.2
;过滤数据
GT250IdxDir = where(WdSpd GT 250*0.2);
WdSpd[GT250IdxDir] = 25.001;
print,max(WdSpd)
;出图&加载色表
SpdImg = IMAGE(WdSpd,RGB_TABLE=39,TITLE=diplayname,POSITION=[0.1,0.1,0.9,1])
grdx= axis('X',LOCATION=[0,0],AXIS_RANGE=[0,360],COORD_TRANSFORM=[0,0.25],MINOR=0,MAJOR=19)
grdy= axis('Y',LOCATION=[0,0],AXIS_RANGE=[-90,90],COORD_TRANSFORM=[-90,0.25],MINOR=0,MAJOR=7)
;COLORBAR
c = COLORBAR(TARGET=SpdImg, ORIENTATION=0,TITLE='Wind Speed (m/s)',POSITION=[0.1,0.1,0.9,0.15])
;SpdImg.save, MyRootDir+diplayname+ImgExtName
;spawn,MyRootDir+diplayname+ImgExtName
;关闭文件释放设备号
CLOSE , LUN
FREE_LUN , LUN
ENDFOR
END
效果如图
这里比较难设置的参数是COORD_TRANSFORM,如果不清楚可以看这里。
说明
文中MyRootDir=’D:\3\’ 是定义的数据路径,这里共享下我用的实验数据 。数据链接:http://pan.baidu.com/s/1i39p5q9 密码:tbvi
相关阅读
声明
1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。
2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。
手机阅读
公众号关注
知识星球
手机阅读
最新GIS干货
私享圈子
请问:我要读取L8 OLI txt文件的几个角度信息,该如何读取指定行或者指定关键字及其所在行的数字呢?
过滤数据为什么要 WdSpd[GT250IdxDir] = 25.001; 而不是0呢?
其中数据在0到250属于有效信息,大于250的值均属于无效信息。
既然是无效值,那直接给0不也可以吗?
给25.001的话不是反而把原来无效单位里的值变成有效的了吗?
0也是有效值,要给0也可以,这里给25.001是因为使用的色表中最大值对应是的白色,这样可视化更容易区分。
原来是这样。不过只有7月1日的色表最大值是25。2,3,4,5,6,7日的数据色表的最大值就不是25了。