GrADS实用绘图--图形显示
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
x, y 包含所有的点 z, t 取第一个点 1.GrADS最擅长画2-D图 2. 全球数据的缺省维数♥♥♥
• 改变维数环境
‘set ‘set ‘set ‘set x 1 10’ y 1 10’ z 6’ t 15’ ‘set lon 120 180’ ‘set lat -20 20’ ‘set lev 500’ ‘set time 00z01jan2007’
3
LASG动力论坛版
四. 图形显示(D)1
本节内容提要
• 基本概念
– 数据文件 – .ctl文件 – 维数环境
— 格点绘图
• 内建图形样式及其定制
– 1-D序列图 – 2-D平面图 – 定制
• 自由图形 • 动画*
4
LASG动力论坛版
四. 图形显示(D)1
基本概念
– 数据文件 – 数据描述文件*.ctl
• 绘图区范围设置
ga-> set parea 1 10.5 1 8
9
LASG动力论坛版
四. 图形显示(D)1
— 实现一页多图
• 将画布分成几个小的虚页
panels_demo.gs ~ 各虚页相对独立, 不能跨虚页作图写字符
• 在一个虚页上设置多个绘图区
pages.gs
• 设置绘图窗口尺寸
‘set xsize x y’
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
• 画布/ 虚页/ 绘图区大小
— 格点绘图
命令顺序 set display draw print
– 图形类型
• 等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
• • • • • 标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
• 动画
– 维数环境 x | y | z | t 中有3维变化
• 站点图
13
LASG动力论坛版
四. 图形显示(D)1
set gxout 图形类型
bar line 直方图 折线图
— 设置图形类型
contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 输出统计数据 二进制屏幕输出(set prnopts) 输出二进制文件(set fwrite)
• • • • •
— 格点绘图
气象要素场按照 一定格式存储 说明数据文件的内容
沿纬线方向的点数 (x方向维数) 沿经线方向的点数 (y方向维数) 沿垂直方向的点数 (z方向维数) 有多少个时间点 (t方向维数) 有多少气象要素/ 每个气象要素的层次
构成 “维数环境”♥♥♥
5
LASG动力论坛版
四. 图形显示(D)1
8
LASG动力论坛版
四. 图形显示(D)1
• “画布” — 横放 / 纵放
$ grads -l 或 $ grads –p
— 设置布局
• 虚拟页面尺寸(单位: 英寸!
ga-> set vpage 0 11 0 8.5 ga-> set vpage 0 8.5 0 11
1inch=2.54cm) grads –l grads –p 画布为Letter纸 大小, 不可超出! 请预留坐标和标 题的位置!
1
LASG动力论坛版
GrADS实用绘图
第二讲 图形显示1
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
2
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1— 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2— 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
set strmden [1~10]5
矢量长度比例
set arrscl inch value
1杠=10 knots!
四. 图形显示(D)1
一页多图范例
pages.gs
panels_demo.gs
(注意区别)
12
LASG动力论坛版
四. 图形显示(D)1
• 设置图形类型 • 1-D图 • 2-D图
— 设置图形类型
set gxout …
– 维数环境 x | y | z | t 只有1维变化 – 单站的探空曲线, 台风年频数序列… – 维数环境 x | y | z | t 中有2维变化 – 水平分布 / 垂直剖面 / 经向剖面 / 纬向剖面 / 任意剖面 时间演变 /…
display u;v;mag(u,v)
风羽图 set gxout barb 风羽绘制方式
set hempref auto ~ ~ nhem ~ ~ shem
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
箭头大小比例
set arrowhead ±inch
流线图 set gxout stream 流线密度
line1<line2: #c1 line1>line2: #c2
误差线图 set gxout errbar d line1;line2 颜 色
set ccolor #color
标记颜色
set ccolor #color
标记大小
set digsize #inch
线 宽
set cthick [1~10]
errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
(李晓峰, 2005)
20
LASG动力论坛版
四. 图形显示(D)1
网格填值图 set gxout grid 小数位数
set dignum number
— 2-D单变量图2
马赛克图 set gxout grfill 颜 色
set ccolor rainbow ~ ~ revrain ~ ~ #color
scatter
set xlint int set ylint int
18
LASG动力论坛版
四. 图形显示(D)1
等值线平滑
set cterp on/ off
on/ off/ linear
~ csmooth
— 2-D单变量图1
等值线图 set gxout contour ~ ~ shaded
等值线颜色
(不对图形起作用)
draw title \CRU Global 1961-1990 Monthly Surface Temperature 10
LASG动力论坛版
set vpage 0 11 0 8.5 set parea 1 10.5 1 7
四. 图形显示(D)1
— 2-D单变量图范例
draw xlab lon draw ylab lat
等值线标记间隔
set clskip #number
线 型
set cstyle #style
等值线标记属性
set clopts
#color thick size
等值线上/下限/省略
set cmin value
set cmax value set black vlo vhi
19
LASG动力论坛版
四. 图形显示(D)1
21
LASG动力论坛版
四. 图形显示(D)1
— 2-D单变量图范例
(张丽, 2006)
22
LASG动力论坛版
四. 图形显示(D)1
矢量图 set gxout vector 矢量比例图例
set arrlab on/ off
— 2-D多变量图
颜 色
set ccolor rainbow ~ ~ revrailled bar linefill
— 1-D序列图示例
line
line
outline bar linefill
‘set t 3 9’ ‘define red=lat-lat’ ‘set t 1 12’ … ‘d red’
set xlevs x1 x2 x3… set ylevs y1 y2 y3…
set ccolor rainbow ~ ~ revrain ~ ~ #color
等值线标记方式
set clab on/ off/ forced/ string/ auto
等值线间隔
set cint int
特征线值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
线 宽
set cthick [1~10]
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
14
LASG动力论坛版
四. 图形显示(D)1
柱间距
set bargap percent
— 1-D单变量序列图
颜 色
set ccolor #color
基 线
set barbase top ~ ~ bottom ~ ~ value
set xlabs b1|b2|b3|… set ylabs b1|b2|b3|…
16
LASG动力论坛版
四. 图形显示(D)1
— 1-D双变量序列图
散点图 set gxout scatter d var1;var2 双线填色图 set gxout linefill d line1;line2 颜 色 set lfcols #c1 #c2
线 型
set cstyle #style
标记大小
set digsiz #inch
set xlopts #color thick size set ylopts #color thick size
LASG动力论坛版
set xlpos offset side set ylpos offset side
set mpdset lowres set mproj latlon set map auto set map 1 1 1 draw map
run cbarn.gs
set grads on set grads off
Set timelab on Set timelab off
11
LASG动力论坛版
‘set csmooth on’ ‘set gxout shaded’ ‘set black 760 780’ ‘d var’ ‘set gxout contour’ ‘set cint 5’ ‘d var’ ‘set clevs 0’ ‘set ccolor 1’ ‘d lat’
— 2-D单变量图范例
直方图 set gxout bar
外观属性
set baropts filled ~ ~ outline
颜 色
set ccolor #color
缺测点
set missconn on ~ ~ off
线 宽
set cthick [1~10]
折线图 set gxout line
标记类型
set cmark #mark
字符大小
set digsize #inch
特征值马赛克图 set gxout fgrid 指定特征值/颜色
set fgvals v1 #c1 v2 #c2 v3 #c3 …
格 线
set gridln auto ~ ~ off ~ ~ #color
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
标记类型
set cmark #marker
线 型
set cstyle #style
17
LASG动力论坛版
四. 图形显示(D)1
errbar
— 1-D双变量图示例
set xaxis strt end int set yaxis strt end int set vrange vlo vhi set vrange2 vlo vhi
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件 在实际空间上的结构
高度
经度
UNIT=20 的数据文 件在实际 空间上的 结构
高度
纬 度
纬 度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
• 网格坐标系 x, y, z, t ~ 如同数组下标 • 实际坐标系 lon, lat, lev, time • 缺省维数环境
.ctl文件相当于翻译官——
— 格点绘图
i 经度 j 纬度 k 高度
Write(10,rec=1) ( ( ( temp (i, j, k), i=1,18), j=1,9), k=1,5) Write(20,rec=1) ( ( ( temp (i, j, k), i=1,5), j=1,9), k=1,18)
• 改变维数环境
‘set ‘set ‘set ‘set x 1 10’ y 1 10’ z 6’ t 15’ ‘set lon 120 180’ ‘set lat -20 20’ ‘set lev 500’ ‘set time 00z01jan2007’
3
LASG动力论坛版
四. 图形显示(D)1
本节内容提要
• 基本概念
– 数据文件 – .ctl文件 – 维数环境
— 格点绘图
• 内建图形样式及其定制
– 1-D序列图 – 2-D平面图 – 定制
• 自由图形 • 动画*
4
LASG动力论坛版
四. 图形显示(D)1
基本概念
– 数据文件 – 数据描述文件*.ctl
• 绘图区范围设置
ga-> set parea 1 10.5 1 8
9
LASG动力论坛版
四. 图形显示(D)1
— 实现一页多图
• 将画布分成几个小的虚页
panels_demo.gs ~ 各虚页相对独立, 不能跨虚页作图写字符
• 在一个虚页上设置多个绘图区
pages.gs
• 设置绘图窗口尺寸
‘set xsize x y’
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
• 画布/ 虚页/ 绘图区大小
— 格点绘图
命令顺序 set display draw print
– 图形类型
• 等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
• • • • • 标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
• 动画
– 维数环境 x | y | z | t 中有3维变化
• 站点图
13
LASG动力论坛版
四. 图形显示(D)1
set gxout 图形类型
bar line 直方图 折线图
— 设置图形类型
contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 输出统计数据 二进制屏幕输出(set prnopts) 输出二进制文件(set fwrite)
• • • • •
— 格点绘图
气象要素场按照 一定格式存储 说明数据文件的内容
沿纬线方向的点数 (x方向维数) 沿经线方向的点数 (y方向维数) 沿垂直方向的点数 (z方向维数) 有多少个时间点 (t方向维数) 有多少气象要素/ 每个气象要素的层次
构成 “维数环境”♥♥♥
5
LASG动力论坛版
四. 图形显示(D)1
8
LASG动力论坛版
四. 图形显示(D)1
• “画布” — 横放 / 纵放
$ grads -l 或 $ grads –p
— 设置布局
• 虚拟页面尺寸(单位: 英寸!
ga-> set vpage 0 11 0 8.5 ga-> set vpage 0 8.5 0 11
1inch=2.54cm) grads –l grads –p 画布为Letter纸 大小, 不可超出! 请预留坐标和标 题的位置!
1
LASG动力论坛版
GrADS实用绘图
第二讲 图形显示1
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
2
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1— 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2— 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
set strmden [1~10]5
矢量长度比例
set arrscl inch value
1杠=10 knots!
四. 图形显示(D)1
一页多图范例
pages.gs
panels_demo.gs
(注意区别)
12
LASG动力论坛版
四. 图形显示(D)1
• 设置图形类型 • 1-D图 • 2-D图
— 设置图形类型
set gxout …
– 维数环境 x | y | z | t 只有1维变化 – 单站的探空曲线, 台风年频数序列… – 维数环境 x | y | z | t 中有2维变化 – 水平分布 / 垂直剖面 / 经向剖面 / 纬向剖面 / 任意剖面 时间演变 /…
display u;v;mag(u,v)
风羽图 set gxout barb 风羽绘制方式
set hempref auto ~ ~ nhem ~ ~ shem
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
箭头大小比例
set arrowhead ±inch
流线图 set gxout stream 流线密度
line1<line2: #c1 line1>line2: #c2
误差线图 set gxout errbar d line1;line2 颜 色
set ccolor #color
标记颜色
set ccolor #color
标记大小
set digsize #inch
线 宽
set cthick [1~10]
errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
(李晓峰, 2005)
20
LASG动力论坛版
四. 图形显示(D)1
网格填值图 set gxout grid 小数位数
set dignum number
— 2-D单变量图2
马赛克图 set gxout grfill 颜 色
set ccolor rainbow ~ ~ revrain ~ ~ #color
scatter
set xlint int set ylint int
18
LASG动力论坛版
四. 图形显示(D)1
等值线平滑
set cterp on/ off
on/ off/ linear
~ csmooth
— 2-D单变量图1
等值线图 set gxout contour ~ ~ shaded
等值线颜色
(不对图形起作用)
draw title \CRU Global 1961-1990 Monthly Surface Temperature 10
LASG动力论坛版
set vpage 0 11 0 8.5 set parea 1 10.5 1 7
四. 图形显示(D)1
— 2-D单变量图范例
draw xlab lon draw ylab lat
等值线标记间隔
set clskip #number
线 型
set cstyle #style
等值线标记属性
set clopts
#color thick size
等值线上/下限/省略
set cmin value
set cmax value set black vlo vhi
19
LASG动力论坛版
四. 图形显示(D)1
21
LASG动力论坛版
四. 图形显示(D)1
— 2-D单变量图范例
(张丽, 2006)
22
LASG动力论坛版
四. 图形显示(D)1
矢量图 set gxout vector 矢量比例图例
set arrlab on/ off
— 2-D多变量图
颜 色
set ccolor rainbow ~ ~ revrailled bar linefill
— 1-D序列图示例
line
line
outline bar linefill
‘set t 3 9’ ‘define red=lat-lat’ ‘set t 1 12’ … ‘d red’
set xlevs x1 x2 x3… set ylevs y1 y2 y3…
set ccolor rainbow ~ ~ revrain ~ ~ #color
等值线标记方式
set clab on/ off/ forced/ string/ auto
等值线间隔
set cint int
特征线值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
线 宽
set cthick [1~10]
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
14
LASG动力论坛版
四. 图形显示(D)1
柱间距
set bargap percent
— 1-D单变量序列图
颜 色
set ccolor #color
基 线
set barbase top ~ ~ bottom ~ ~ value
set xlabs b1|b2|b3|… set ylabs b1|b2|b3|…
16
LASG动力论坛版
四. 图形显示(D)1
— 1-D双变量序列图
散点图 set gxout scatter d var1;var2 双线填色图 set gxout linefill d line1;line2 颜 色 set lfcols #c1 #c2
线 型
set cstyle #style
标记大小
set digsiz #inch
set xlopts #color thick size set ylopts #color thick size
LASG动力论坛版
set xlpos offset side set ylpos offset side
set mpdset lowres set mproj latlon set map auto set map 1 1 1 draw map
run cbarn.gs
set grads on set grads off
Set timelab on Set timelab off
11
LASG动力论坛版
‘set csmooth on’ ‘set gxout shaded’ ‘set black 760 780’ ‘d var’ ‘set gxout contour’ ‘set cint 5’ ‘d var’ ‘set clevs 0’ ‘set ccolor 1’ ‘d lat’
— 2-D单变量图范例
直方图 set gxout bar
外观属性
set baropts filled ~ ~ outline
颜 色
set ccolor #color
缺测点
set missconn on ~ ~ off
线 宽
set cthick [1~10]
折线图 set gxout line
标记类型
set cmark #mark
字符大小
set digsize #inch
特征值马赛克图 set gxout fgrid 指定特征值/颜色
set fgvals v1 #c1 v2 #c2 v3 #c3 …
格 线
set gridln auto ~ ~ off ~ ~ #color
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
标记类型
set cmark #marker
线 型
set cstyle #style
17
LASG动力论坛版
四. 图形显示(D)1
errbar
— 1-D双变量图示例
set xaxis strt end int set yaxis strt end int set vrange vlo vhi set vrange2 vlo vhi
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件 在实际空间上的结构
高度
经度
UNIT=20 的数据文 件在实际 空间上的 结构
高度
纬 度
纬 度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
• 网格坐标系 x, y, z, t ~ 如同数组下标 • 实际坐标系 lon, lat, lev, time • 缺省维数环境
.ctl文件相当于翻译官——
— 格点绘图
i 经度 j 纬度 k 高度
Write(10,rec=1) ( ( ( temp (i, j, k), i=1,18), j=1,9), k=1,5) Write(20,rec=1) ( ( ( temp (i, j, k), i=1,5), j=1,9), k=1,18)