气候诊断及绘图实用技术讲座.ppt

1、气候诊断和制图实用技术讲座陕西省气候中心李智雅2011年三月18日，晋中市气象局、数据采集气候诊断方法图纸分析、第一部分数据采集、气象资料：网站、省信息中心、国家气候中心网站(160站降水和气温资料等)、中国气象科学数据孔刘服务网。()大气循环场，海温资料：NCEP、EC、海温、气候监测公报(国家气候中心、NCEP)、74个特征量等。第二部分气候诊断，(1)资料的字典处理(均匀性检查等)；(2)元素的年际(时间演化曲线)和年代(滑动平均)变化分析(3)元素的空间分布特性：EOF扩展等；(4)典型年份的选择(大于或小于特定阈值的标准偏差)(5)原因或机制分析：充分利用EXCEL进行相关分析、合成

2、分析、初步统计分析：序列中的各种统计；趋势分析、相关分析；利用现有统计分析节目，如数据标准化、过滤、分类摘要：目前广泛用于气候诊断方法，如EOF、M-K突变检测、沙发分析等。下面介绍M-K、修改沙发分析的程序、准备数据、分析结果、绘制相关问题。1大气科学研究的图，最重要的两种茄子类型是时间演化图(图、柱、山点等)、等值线分布图。一，几个茄子天气绘图软件介绍，第三部分绘图，第一类图表比较简单，最常用的Excel等软件制作，以及满足更高要求时，可以使用俄勒冈等图形软件。绘制第二种等值线也有多种软件选项，如Micaps、GrAds、Matlab、Surfer等。(1) GrADS图简单快速，提供了许

3、多大气领域中常用的函数，如旋涡、季度、纬度平均值、面积平均值。数据格式是二进制的，有时很难找到错误。画城堡的分布图很不方便。(2) Surfer surfer简单实用，适用于绘制各种区域的分布。绘制二边曲线，主要利用EXCEL(ORIGIN)元素的年代际变化曲线；对于散点图部分气候诊断结果(M-K等)，Surfer的主要功能是绘制等值线图，包括要素的趋势线、滑动平均线，还可以绘制贴图、分类贴图、矢量贴图、图像贴图、线框、曲面贴图等格式。3，Surfer绘图软件，山西1960-2008年65站夏季降水量平均百分比时间剖面图，山西1960-2008年62站夏季气温距离平面时间剖面图，(1)等值线图

4、绘制，Surfer最重要的功能是绘制等值线图，但我们拥有的是也就是说，必须先将数据档案转换为Surfer识别的GRD文件格式，然后才能绘制等值线。当然，您可以直接生成Surfer允许的ASCII代码的GRD档案格式，并自行绘制。稍后将介绍牙齿方法。假设有三列数据：x、y和z。其中z是点(x，y)的值，* .存在dat文件。其中，第一列是x坐标，第二列是y坐标，第三列是(x，y)的值z，这是远视数据的图形，如果您不需要插值，最好选择加权倒距离方法或克里格方法。这两种牙齿茄子方法是因为插值点与采样点匹配时，插值点的值是采样点的值，但其他方法无法保证这一点。输出网格文件中的输出档案名称* .输入GR

5、D，然后设定网面线几何图形中的网面点数。此处需要注意的是，如果x和y的数值差异很大，此处显示的最大最小值可能有错误(即，与远视数据不同)。这是Surfer软件本身的问题。如果牙齿，则必须手动更正错误。也就是说，必须输入准确的最大最小值。如果数据不存在此问题，则无需手动更正最大最小值。但是，必须手动修改X和Y的间距或数据行数。两者都相关，一个需要更改，另一个需要自动修改。如果原始数据的间距相等，则x和y之间的间距或数据行数最好与原始数据匹配。这样可以减少插值产生的错误。我们的数据不均匀，所以需要插值。在这里不需要做任何更改。单击“确定”(ok)后，图中所需的GRD档案* .GRD创建完毕。步骤2

6、:等值线绘制打开菜单映射|等值线映射|新建等值线映射，在“打开网格”(open mesh)对话框中，输入刚导出的GRD档案(例如* .选择GRD，点“确定”完成等值线图。克里格插值、最小曲率插值、步骤3:设置等值线绘制等值线图中的双击，或单击鼠标右键，然后选择“属性”(attributes)，可以设置等值线的各种属性。双击图片的边缘可更改轴的属性。(b)创建bln档案(数字化地图)1，输入地图，使用map|basemap命令选择地图，双击地图，然后在显示的剪辑属性中修改图像坐标中的最北、最北指示器。2，可以使用fileimport引入一张照片，但不能数字化。要数字化，必须使用mapbase m

7、ap。数字化之前，最好定义方形的坐标。您可以使用区域区域或彩色jpg档、使用mapinfo建立参考底图，或使用Google寻找并储存某个位置的地图，然后进入surfer进行数化。3，bln档案格式长度，flag x1，y1 x2，y2.xn，Yn=y1长度，flag x1，y1 x2，y2.xn，yn，其中length是未闭合的曲线，否则。曲线中点坐标的数目，flag值为1或0。flag=1时，可以将闭合曲线内的等值线降至牙齿bln档案blang(删除闭合曲线中的等值线)。flag=0时，曲线外部的等值线由flag删除。例如，我提供的bln文件格式为flag=0，当然可以更改为1。在flag=

8、1的情况下，海温场的等值线，陆地必须是blank落下。远视数据文件已被命名为“降水克里格。假设有一张grd 的山西降水的等值线图。打开轴网|空白.降水克里格。选择“grd”、“ok”文件，然后在新对话框中选择“sxdigit”以确认点，从而获得blank的GRD文件，out .用GRD制作等值线地图(图2)，最后叠加山西的基线图，成为我们最后需要的等值线(图3)。4 bln档案和blank映射、将基准图形添加到5等值线地图、Surfer中的地图|基准图形.和SXDIGIT。打开BLN文件后，可以绘制地图。布满视图|窗口时，显示效果会很好。然后，可以将相应的等值线图绘制为第二部分。基准图面的格线

9、座标必须与等角线图的格线座标相符。至少，基准图形的栅格坐标必须包含素线图形的栅格坐标。再次按F2键(“选择所有对象”或CTRL A)，然后单击“地图|重叠”，将两个地物合成为一对地物。打开视图|管理对象、选择基准、单击鼠标右键、选择堆叠顺序|前进、将背景地物移动到上一级别。这样好的等值线地图画得很好。，图1，图2，图3，最后，我们总结一下牙齿过程。(1)在任意位置创建bln文件。(2)应用数据网格(* .生成GRD档案)；(3)将被检查为bln文件的文件白化，以获得blank的GRD文件。(4) blank的*。用GRD做等值线。(5) basemap(* .bln)来调用Overlap。(c

10、)，章节映射和分类章节映射，有时需要将逆坐标添加到基准图形，将它们标记为三角形、五颗星等，并在旁边写下逆名称。这可以用章节图和分类章节图完成。1，首先，资料档案一般格式为： x1，y1，label1 x2，y2，label2.介绍Xn、Yn和labeln的章节图。其中xn，yn是坐标值，label是数字或字母，2，章节贴图分类类似于章节贴图，但是对于其他点，根据特定规则进行分类。例如，在test.dat文件中，如果将C列视为每个点的高程，则可以根据每个点的高程对每个点进行分类。高程为25-54时用十字标记，88-105时用三角形标记等。打开地图|章节映射|新建分类章节映射.选择test.dat

11、档案，然后选择确定，将显示分类章节映射。可以双击图片来编辑属性。此处有多个分类选项卡。您可以在此变更选项最大最小值、符号等。不再详细说明。分类章节图必须至少包含三列数据，其中第三列是分类的标准。如果surfer在blank的后边缘遇到空解决方案，则可以提高网格密度，包括使用插值方法对原始数据文件进行更紧密的网格，以解决牙齿问题。边缘区域几乎没有空格。surfer如何将坐标的单位设置为中文？上下轴直接输入中文就可以了。字体最好选择宋体。否则，汉字之间的间距可能会更大。左轴：(1)首先在Title中输入中文，如“纬度”。在这里，根据中文数组格式输入的汉字必须作为回车处理。要从下到上横向显示汉字，必

12、须将“纬度”更改为“纬度”，然后移动光标。(2)在angle中拔出270度是纵轴汉字破碎的主要原因，选择270度一般不会破碎。(3)选择字体为“宋体”，注意宋体前面有“宋体”。如果没有这种字体，就要安装汉字库。总之，请注意三个茄子点：(1)上下轴角度为0度，左右轴为270度。(2)选择字体“宋体”或“宋体”(3)修改汉字的排列格式、等值线数字显示。贴图| contourmap | edit labels，4，GrAds绘图，需要准备的档案映射：*。nc(或* .grd)(数据档案)*。CTL *。Gs * 1基本命令应使用，并应善用快速参考手册。2更多练习，创建更多GS档案。、第一部分NCEP

13、数据下载和NCEP数据简介，http:/WWW。ESRL . NOAA . gov/PSD/DATA/Gridded/NCEP重新分析资料EC欧洲重新分析资讯；要素：基本元素(位高、温度、湿度、相对湿度、U、V风场、垂直速度)两个物理量的协方差、流函数、速度潜力等其他一些导出、区域范围3360全局；时间范围：1948-当前；从1979年到现在的时间标准：月、日、日4次数据；晶格：2.5X2.5；1X10.25X0.25垂直标高：17标高(1000、925、850、700、600、500、400、300、250、200、150、100、70、70).set gxout shaded * * * *设定最大等高线值。* * * * * * * * * * set Cmax 0d vor *不再绘制Set gxout shaded * * *设定最小等高线值。* * * * * * * * * * * set cmin 0d vor * 100000以下的值不再绘制Reinit SDF open d :/2010/month/HGT . mon . mean . NC set let，两阶段计算相关，fortran，编写三阶段CTL档案，Dset d :/201