Grads处理风云FY2静止卫星AWX数据总结

wangfei

{:5_213:}{:5_213:}

尉英华

多谢楼主分享

听海落雪

楼主，你好。我想画 FY2G的PRE，出图的时候遇到问题，求助。

ctl如下：
dset F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001_OTG_2017%m2%d2_%h200.AWX
title 2017 FY2G_PRE
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1001 linear 55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 48 linear 00z23JUN2017 1hr
vars 1
pre 0 -1,40,1 FY2G_PRE(m)
endvars

gs如下：

'reinit'

'open F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001.ctl'

'set grads off'
'set grid off'

'set lon 96 106'
'set lat 25 33'

****** 1小时降水 ******
'set rgb 54 255 255 255'  
'set rgb 55 220 250 211'  
'set rgb 56 166 242 143'  
'set rgb 57 61 186 61'   
'set rgb 58 97 184 255'   
'set rgb 59 0 0 225'      
'set rgb 60 250 0 250'   
'set rgb 61 128 0 64'

********* 一页四图 *********

p=1
while(p<=4)

'run page q'p''
'set t 'p''

'set parea 0.5 9.0    1.0 8.0'
'set xlopts 1 5 0.26'
'set ylopts 1 5 0.26'
'set xlab on'
'set ylab on'
'set xlabs 96`3.|98`3.|100`3.|102`3.|104`3.|106`3.`1E'
'set ylabs 25`3.|27`3.|29`3.|31`3.|33`3.`1N'

'set gxout shaded'
'set csmooth on'
'set clevs 0 0.1 1.5 7 15 40 50'
'set ccols 54 55 56 57 58 59 60 61'
'd smth9(pre)'

'set mpdset country1'
'set map 15 1 2'
'draw map'
'set mpdset xnqy_sj'
'set map 15 1 2'
'draw map'
'set mpdset pxdq_dj'
'set map 8 1 5'
'draw map'

'query w2xy 102.2147369385 27.9398202896'
x1=subwrd(result,3)
y1=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.16'
'draw mark 1  'x1' 'y1'  0.25'

'query w2xy 103.4 32.7'
x11=subwrd(result,3)
y11=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.28'
'draw string 'x11' 'y11'  `2  'p+13'BJZ23JUN'

p=p+1
endwhile

************************************************

'printim F:\Data\GrADS\2017\20170624\pic\PRE_2314-2317.png white'
*'enable print F:\Data\GrADS\2017\20170624\pic\PRE_2314-2317.gmf'
*'print'
*'disable print'

aentera

hasugi

谢谢楼主分享

funnybone

思路和分享挺好的，第一张图CLC产品的地理定位从图上看貌似有一点偏移

baggio

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

baggio

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

hqk1105

雪颜 发表于 2016-9-13 15:44
各位吧友注意不要照搬范围，具体看卫星星下点的位置

60°N-60°S， 45°E-165°E

请问如何通过星下点知道经纬度范围？

hqk1105

雪颜 发表于 2016-9-13 15:44
各位吧友注意不要照搬范围，具体看卫星星下点的位置

60°N-60°S， 45°E-165°E

请问经度范围是否是由定位点加减60得到呢？