GrADS绘图的各种设置

MeteorologicalGraphicsTutorialGrADS绘图的各种设置从维数环境到图形输出的全参数配置教程Contents教程目录GrADS绘图系统的完整学习路径，从基础概念到脚本编程的八章核心内容。01GrADS基础概念与工作流程02数据文件准备与描述文件编写03维数环境设置04图形类型与要素配置05坐标轴与标注设置06地图投影与区域控制07变量、函数与表达式08脚本编程与图形输出CHAPTER01GrADS基础概念与工作流程理解网格数据结构、文件体系与软件操作逻辑OverviewGrADS软件概述与核心能力GrADS（GridAnalysisandDisplaySystem）是由NASACOLA实验室开发的气象数据可视化系统，以命令行驱动、脚本可编程的方式，为气象和海洋领域提供专业的网格数据分析与图形展示能力，是全球气象科研与业务的主流工具之一。气象科研数据分析工作场景多格式数据支持GRIB·NetCDF·HDF支持GRIB、NetCDF、HDF等多种气象标准数据格式，可通过数据描述文件（.ctl）灵活定义数据结构和维度映射精细图形控制200+set命令提供超过200条set命令参数，覆盖等值线、颜色填充、矢量场、剖面图等全部气象绘图类型，支持逐像素级图形控制脚本自动化.gs脚本内置GrADSScript描述语言（.gs脚本），支持变量定义、条件判断、循环结构和外部函数调用，实现自动化批量绘图全球广泛应用1988年至今自1988年发布以来持续迭代至2.x版本，在全球气象科研院所、高校大气科学专业和数值预报业务中心广泛部署使用FILESYSTEMGrADS常用文件类型体系GrADS的运行依赖一套完整的文件体系：数据文件存储网格观测或模式输出，数据描述文件定义维度映射，脚本文件实现自动化，图形输出保存可视化成果。DATAFILES数据文件存储网格点数据或站点观测数据，支持GRIB、NetCDF、IEEE二进制等多种格式.grb·.nc·.datCONTROLFILES数据描述文件纯文本格式，定义DSET路径、UNDEF缺测值、TITLE标题、四维坐标及变量列表.ctlSCRIPTFILES脚本文件GrADSScript描述语言编写，支持变量、函数、循环和条件判断，用于批量处理与自动绘图.gsGRAPHICOUTPUT图形输出文件printim输出PNG/JPEG位图，或enableprint与gxprint输出PostScript矢量图，用于论文插图和业务产品PNG·PSGrADS标准工作流程GrADS绘图的各种设置GrADS的操作遵循"启动→打开数据→设置环境→绘图显示→输出保存"的五步标准流程。其中"设置环境"环节是绘图质量控制的核心，通过set命令在进入display命令之前预先配置好维数范围、图形类型、颜色方案等全部参数，确保一次性输出符合要求的科学图形。启动系统启动GrADS后进入命令行交互窗口，默认为landscape横向模式；可通过'grads-blc'启动无GUI批处理模式，适合服务器端批量运行-blc打开数据使用'openxxx.ctl'打开数据描述文件，若需叠加多套数据可用'open-e'追加打开，系统自动分配文件编号供后续变量引用.ctl设置环境通过一系列set命令预设维数环境、图形类型、颜色表、坐标范围等绘图参数，所有set设置持续生效直到被新命令覆盖或用'clear'清除set绘图输出使用'display变量名'（简写'd'）触发绘图，GrADS根据当前set环境自动渲染图形；可用'printim'导出指定分辨率图像displayChapter02数据文件准备与描述文件编写掌握.ctl描述文件的结构与关键字段配置数据描述文件数据描述文件（.ctl）核心结构.ctl描述文件是GrADS读取数据的唯一入口，它以纯文本形式定义了数据文件路径、缺测值标记、四维空间网格（经度/纬度/高度/时间）的划分方式以及变量列表。DSET·数据路径与选项DSET指定数据文件路径（支持^符号表示与ctl同目录），OPTIONS可附加yrev/zrev/byteswapped等标记处理特殊数据排列方式DSET/OPTIONSUNDEF·缺测值定义UNDEF定义缺测值标识（常用-9.99E33），GrADS在绘图时自动跳过该值，避免缺测区域被错误渲染为有效数据-9.99E33XDEF/YDEF·四维网格XDEF/YDEF以LEVELS方式逐一列出经度/纬度值或以LINEAR方式指定起始值和等间距格距，ZDEF和TDEF同理定义垂直层次和时间轴LEVELS/LINEARVARS·变量列表VARS区块列出所有变量，格式为"变量名层数单位描述"，层数为0表示二维场、正整数表示三维场，GrADS据此确定每个变量的数据维度0=2D·N=3DDATAFORMATACCESS常见格点数据格式的接入方式GrADS支持GRIB、NetCDF、HDF和自定义二进制等多种格点数据格式，不同格式的接入路径和配置要点各有差异。GRIB/GRIB2需先用gribmap-ixxx.ctl生成.idx索引文件，.ctl中用DTYPEgrib标记并通过INDEX指向索引，支持多变量多层次自动解析。gribmapNetCDFGrADS2.0+支持sdfopenxxx.nc直接打开，系统自动读取CF约定的维度与变量元数据，省去手动编写ctl的步骤。sdfopenIEEE二进制需在ctl中用OPTIONS指定big_endian或little_endian字节序，DTYPE无特殊标记，是最基础的数据接入方式。OPTIONS多文件拼接通过THEMINIT和OPTIONStemplate标记，可用时间模板自动拼接多个时次的独立数据文件为连续时间序列。templateChapter03维数环境设置精确控制绘图的经度、纬度、高度和时间范围DIMENSIONENVIRONMENT维数环境的核心概念与基础命令维数环境是GrADS绘图前必须预先定义的数据选取范围，通过setlon/lat/lev/time四个命令分别设定经度、纬度、垂直层次和时间轴的取值，维度设为单一值时取切片绘制平面图，设为范围时沿该维度展开绘制剖面图。setlon定义东西方向数据范围，支持360°制和-180°~180°制两种经度表示70°E—140°Esetlat定义南北方向数据范围，常配合中国区域绘图使用15°N—55°Nsetlev选取单一等压面绘制平面图，或设范围绘制垂直剖面500hPasettime指定具体时次，或设范围绘制时间—经度/纬度剖面图HovmöllerDimensionEnvironment维数环境的进阶操作与常见陷阱维数环境的"固定"与"释放"机制决定了GrADS在每次display时如何选取数据维度。正确管理维数状态是避免绘图错误的关键。01维度固定模式setlon116将经度锁定为单一值，适合绘制经度–时间剖面或纬度–高度剖面，固定维度不出现在图形坐标轴中。02维度释放模式setlon70140将经度恢复为范围，GrADS在此范围内自动插值并填充整个绘图区域，释放维度对应图形的坐标轴。03querydims命令实时查询当前四个维度的设定状态（固定/释放/具体数值），是调试绘图问题的必备工具，建议每次display前执行确认。04常见陷阱连续绘制多张图时若未重置setlev，后续图形会继承前一张的层次设定——应养成显式设置全部四维的编程习惯。CHAPTER04图形类型与要素配置掌握等值线、填色图、矢量场、流线图等核心绘图类型及其参数调节SETGXOUT图形输出类型（setgxout）全面解析setgxout命令是GrADS绘图类型选择的总开关，决定了display命令以何种图形形式渲染数据。从等值线到颜色填充、从风矢量到流线图，GrADS提供了十余种图形输出类型，每种类型对应一组专属的美化参数。正确选择gxout类型是绘制专业气象图形的第一步，也是后续所有参数调优的基础。contour绘制数据场的等值线，可配合setccolor/cstyle/cthick/cmark调节线条颜色、虚实、粗细和标记点，适合展示气压场和高度场ccolor·cthickshaded/shaded2以色阶填充等值线间隔区域，shaded2支持平滑过渡效果，配合setclevs和setccols定义色标分级与配色方案clevs·ccolsvector用箭头展示二维矢量场（如风场），箭头长度代表风速、方向代表风向，通过setarrsiz/arrmrg控制箭头尺寸和间距arrsiz·arrmrgstream绘制连续流线展示矢量场方向特征，比vector更直观表现气流走向，通过setstleng/stden控制流线长度和密度stleng·stdenCONTOURPARAMETERS等值线参数精细设置等值线图是气象分析中最核心的图形类型，其外观由一组set参数联合控制：clevs定义等值线数值分布，ccols定义对应颜色，cstyle定义线型，cthick定义线宽，clab控制标注方式。这些参数的精确配置直接决定了等值线图的专业性和可读性。setclevs值列表setclevs-30-20-100102030setcint逐一指定等值线数值，如setclevs-30-20-100102030；或用setcint按等间距自动生成，两者不可同时生效。合理设置数值间隔能突出关键气象特征。clevs·cint·数值间隔setccols颜色列表rgb2551280与clevs一一对应指定每条等值线颜色，支持1-15号标准色和RGB自定义色（如rgb2551280），颜色数量应比clevs多1。科学的配色方案能直观反映数据梯度变化。1-15·RGB·色标映射setcstyle线型1为实线（默认）、2为长虚线、3为短虚线、4为点划线，可对每条等值线单独指定或统一设置，常用于区分正负值区域。线型选择需兼顾打印输出与屏幕显示效果。实线·虚线·点划线setclab标注控制on显示等值线数值标注（默认）、off关闭标注、auto自动选择标注位置、forced强制在每个网格间隔标注，影响图形整洁度。标注密度需根据等值线疏密程度灵活调整。on·off·auto·forcedGRADS色标系统颜色填充图与色标配置颜色填充图通过色阶展示数据场空间分布，核心逻辑：clevs定义断点、ccols指定颜色、cbcols/cbflags控制色标条显示。setclevs定义色标断点如setclevs-30-20-100102030创建7个色级区间，断点数量决定色阶精细程度，通常15–25个色级较为合适15–25级setccols配色方案支持1–15号标准色（1白2红3绿4蓝等），更推荐rgbRGB自定义精确色值，可导入NCL色表实现专业渐变配色RGB自定义shadedvsshaded2shaded按网格单元填充（阶梯效果），shaded2在等值线间做平滑插值（渐变效果），后者视觉效果更细腻平滑插值cbset控制色标条cbseton/off开关色标条，cbcols设置颜色，配合setxlopts/ylopts调整位置和字体，是出图美化最后一步出图收尾WindFieldVisualization矢量场与流线图设置风场可视化是气象绘图的核心需求之一，GrADS通过vector和stream两种图形类型实现，两者均可与颜色填充图叠加使用，形成学术论文中最常见的专业气象图形式。01Vector箭头参数setarrsiz调节箭头尺寸（默认0.05），setarrmrg控制箭头间距防止过密重叠，setarrlabon显示右下角参考箭头及对应风速标注arrsiz·arrmrg02Stream流线参数setstleng控制流线追踪长度（默认0.25英寸），setstden设置流线种子点密度（默认5），值越大流线越密集，展示细节更丰富stleng·stden03Barb风向杆模式setgxoutbarb绘制标准气象风向杆，短划线5kt、长划线10kt、三角旗50kt，是地面天气图和高空图的传统风场表示方式5·10·50kt04复合叠加绘图先dmag(u,v)画风速填色图，再setgxoutvector+du;v叠加风矢量，利用双层display实现填色底图与矢量场的专业叠加效果填色+矢量Chapter05坐标轴与标注设置精细控制坐标轴样式、刻度标注、标题文字和图例配置GrADS·绘图设置坐标轴范围与刻度样式设置坐标轴的精确控制是专业气象图形的基本功。GrADS通过三组参数配置让坐标系统清晰规范，满足学术期刊对插图坐标标注的严格要求。setxaxis显式定义X轴范围和刻度间距。如setxaxis7014010在东经70–140°范围内每隔10°标注一个刻度值。70°—140°×10°setyaxis同理定义Y轴。如setyaxis15555在北纬15–55°范围内每5°标注，适合中国区域地图的标准纬度设置。15°—55°×5°setxlopts/ylopts控制坐标轴标注的字体颜色、线宽和大小。设为off可完全隐藏标注，常用于多图拼接时去除冗余坐标。Color·Width·Sizesetgrid在绘图区叠加经纬网格线。如setgridon33用绿色虚线绘制网格，帮助读者精确定位数据点的地理坐标。on/off·Color·StyleGRADSTEXT&ANNOTATION标题、图注与文字标注设置图形的标题、图注和自定义文字标注是信息传达的重要组成部分。合理配置这些文字元素能让图形自解释性更强，是制作可独立发表的学术插图的关键环节。居中放置drawtitle在图形顶部居中放置标题文字，字体和大小由当前setstrsiz和setfont决定，支持多行标题与自动换行处理setfontsetstrsiz坐标定位drawstring在指定页面坐标位置放置文字标注，采用归一化坐标（0-11×0-8.5英寸），适合添加区域名称或数据说明归一化坐标区域标注0.06–0.10insetstrsiz控制后续所有文字标注的水平和垂直方向大小（英寸），默认0.08×0.07，论文插图建议0.06–0.10水平尺寸垂直尺寸底部控制setannot控制图形底部默认标注行的颜色和线宽，设为off可隐藏；setdbfileoff关闭数据浏览器标注颜色控制线宽设置CHAPTER06地图投影与区域控制配置地图投影方式、海岸线、国界省界等底图要素MAPPROJECTION地图投影方式选择与配置投影决定球面到平面的映射方式，选对投影并精确配参是专业气象制图的基础。LATLON等经纬度投影默认投影方式，经纬线呈直线正交网格，适合中低纬度区域展示，高纬度地区面积变形显著正交网格LAMBERT兰伯特等角圆锥投影形状保真度高的圆锥投影，适合中纬度大陆区域，通过setmpvals设置两条标准纬线控制变形标准纬线ROBINSON罗宾森投影折衷型全球投影，面积和形状变形均衡，适合全球气候分布图展示，视觉效果自然美观全球折衷ORTHOGRAPHIC正射投影模拟从太空看地球的球面透视效果，仅展示半球范围，适合台风路径、极涡分析等球面视觉场景半球透视MapBasemapConfiguration地图底图要素精细控制地图底图的海岸线、政治边界和河流等要素是气象图形的地理参照框架，其精度和样式直接影响图形的专业性和信息完整性。setmpdset<精度>选择底图数据分辨率：lowres（默认，约1°精度）/mres（约0.5°）/hires（约0.1°），论文出版推荐使用hires获得清晰的海岸线hiressetmap<颜色><线型><线宽>控制海岸线和大湖边界的显示样式，如setmap115用白色实线中等粗细绘制，在深色填色图上形成清晰对比115setmpdrawon/off控制是否在绘图时自动叠加地图底图轮廓，设为off可绘制纯数据图（如极坐标图或不需地理参照的剖面图）on/off自定义地图数据通过setmpdset加载高精度GSHHS海岸线或自定义行政边界文件，补充内置数据中缺失的省界、岛屿和特殊区域边界GSHHSCHAPTER07变量、函数与表达式利用内置函数和表达式在绘图时实时计算派生物理量GrADS绘图设置变量引用与基本运算表达式GrADS支持在display命令中直接编写数学表达式，实现数据引用、单位换算和简单物理量计算的"即时绘图"，省去预处理步骤，让绘图流程更加高效灵活。单文件引用直接使用变量名如dt显示温度场、du;v同时显示U/V风分量，GrADS自动按当前维数环境提取对应数据切片。变量名需与数据文件中的定义一致，支持标准气象变量缩写。dt多文件引用.文件编号区分来源如dt.2-t.1计算两个文件温度场的差值，适合模式对比分析和异常场计算。文件编号对应open命令的顺序，从1开始计数。t.2−t.1算术运算四则运算与函数嵌套支持+−*/和括号嵌套，如d(t-273.15)将开尔文转摄氏度，dsqrt(u*u+v*v)直接计算风速。内置数学函数包括sqrt、log、exp、abs等。+−×÷()条件绘图maskout条件屏蔽如dmaskout(t,t-273.15)只在温度大于0℃的区域显示数据，用于突出特定物理意义的区域。第一个参数为待显示变量，第二个为条件表达式。maskout()FunctionsLibrary内置科学函数与气象专用函数GrADS内置完整函数库，可在display命令中直接调用。气象专用函数让复杂物理量诊断一行命令即可完成，空间统计函数支持区域平均和积分运算。数学基础函数abs()绝对值、sqrt()开方、log()/log10()对数、pow()幂运算、sin()/cos()三角函数，可在表达式中自由嵌套组合使用sqrt·log·pow气象诊断函数rh(t,td)相对湿度、tlcl(t,td)抬升凝结高度、hcurl(u,v)垂直涡度、hdivg(u,v)水平散度，一行命令完成复杂物理量诊断rh·hcurl·hdivg空间统计函数ave()区域平均、sum()累积和、mag(u,v)矢量模（风速），用于提取区域特征值；const()缺测值处理，maskout()掩码屏蔽ave·sum·magCHAPTER08脚本编程与图形输出用GrADSScript实现自动化批量绘图与高质量图形导出SCRIPTINGFUNDAMENTALSGrADSScript描述语言基础语法GrADSScript（.gs）是一种轻量级脚本语言，支持变量赋值、条件判断、循环结构和字符串操作，是将重复性绘图操作自动化的核心工具。变量与赋值var=valueresult=subwrd(line,2)var=value直接赋值，result=subwrd(line,2)提取字符串中的第2个词，变量无需声明类型，支持整数、浮点和字符串自动推断AutoTypeInference条件判断if/else/endifwhile(rc!=0)if/else/endif结构支持条件分支，如判断文件是否存在后决定是否跳过绘图，while(rc!=0)循环检测命令执行状态ConditionalBranching循环结构i=1;while(i<=24);...;i=i+1;endwhilei=1;while(i<=24);...;i=i+1;endwhile实现24时次遍历绘图，循环体内可动态拼接时间字符串切换维数环境，实现全自动批量出图BatchProcessing命令执行与交互rc=resultpromptreadfilenamelinerc=result捕获上一条命令的返回值，prompt实现用户交互输入，readfilenameline从外部文件逐行读入参数配置CommandInterfaceOUTPUTMETHODS图形输出方式与分辨率控制GrADS提供位图（printim直出PNG/JPEG）与矢量（PS/gxprint输出EPS）两条路径，分辨率与格式选择直接影响最终呈现质量。printim位图输出指定像素分辨率输出PNG格式，支持白底与透明底参数，日常分析推荐800×600起步PNG/JPEGPostScript矢量输出enableprint开启记录，执行print逐页打印，disableprint关闭，可转EPS/PDF用于排版EPS/PDFgxprint高质量输出GrADS2.x直接输出EPS彩色矢量图，支持CMYK色彩模式，期刊投稿最推荐方式CMYK批处理模式输出以无GUI模式启动，配合脚本循环输出数百张图形，适合后处理与业务产品自动生成-blcGRADS·多子图排列多图布局与子图排列设置多图布局是学术论文中对比分析的常见需求，GrADS通过setparea精确定义每张子图在页面上的位置和尺寸，配合clear命令逐张绘制实现灵活的复合图形排列。区域定义setparea以英寸为单位精确定义绘图区域边界，控制图形在页面中的位置与尺寸。setparealrbt排列策略依次执行clear→setparea→set维数→display循环，精确计算坐标实现网格排列。clear→set→display虚拟页面setvpage定义虚拟页面区域，适合将不同比例的子图放在同一输出文件中。setvpageoff协调设置多子图共享色标时统一clevs，共享坐标轴时仅保留最外侧标注以避免重复。clevs/xlopts/yloptsAnimation&TimeSeries动画显示与时间序列设置通过setloop控制播放模式并逐帧输出PNG序列，可借助ffmpeg合成高质量动画视频展示气象场时间演变。setloop模式设置动画循环方式：foreward正向、backward反向、foreback来回振荡，配合display命令自动遍历所有已设时次，实现数据场动态回放。foreward·backward·foreback交互式动画控制图形窗口提供播放/暂停/快进/快退按钮和帧率调节滑块，实时预览数据场的时间演变，适合探索性分析与参数调试。GUIControls逐帧PNG输出在.gs脚本中用while循环遍历每个时次，依次执行settime→display→printim输出单