温度等值线图的自动化绘制研究.doc

温度等值线图的自动化绘制研究 1 前言　　　　近年发展起来的大气数值模式，基于大气动力学和热力学基础，能够较好的模拟不同地表状况对大气环境的影响状况，通过对比分析，我们就能够准确定量评价某地表类型对大气温度环境的影响程度，从而可为绿地建设规划和改造提供指导，但在建立城市小尺度气象数值模式用于城市地表温度场的模拟过程中，同时还需要进行实际地表温度要素的观测，并绘制其等值线图表，通过对比分析以调整模式内部参数值，从而使模拟值更加接近真实。　　为了获取了局部区域不同下垫面的地表温度分布图，通过对比分析，选择了Surfer 软件进行绘制，Surfer 软件作为Golden Software 公司研制的基于Windows 系列操作系统的插值绘图软件，该软件可处理离散点资料、格点资料，可以绘制等值线图、三维立体图等，也可添加研究区域的部分地理信息、在资料点上显示标注等[1,2]。　　从 Surfer 7.0 开始，软件就增加了Automation 技术，该Automation 技术为ActiveXAutomation 技术，提供了近60 种对象[3,4]，Surfer 7.0 的二次开发平台Scripter[5]可通过进行深入快速开发，其核心编程语言The Sax Basic Language 与Visual Basic（VB）编程语言兼容，还可通过VB 编程语言将其强大的功能嵌入到用户开发的系统中。　　针对小尺度气象模式应用时对高精度格点地表温度实测值的验证需要，即尽量直观与客观地给出实际地表温度的空间分布，以天坛公园为例，本研究，基于公园内不同下垫面的实地观测所获取的地表温度数据，通过Surfer7.0 提供的Scripter 脚本编辑器，应用Surfer7.0提供的Automation 技术，选用克立格法[6]插值方法进行了观测数据的空间插值与区域地表温度等值线图的制作，最终获取了局部区域不同下垫面的相对符合实际的地表温度分布图，基本能够满足城市小尺度气象模式的验证需要，对相关离散数据的空间插值方法也有一定的借鉴意义。2 数据的准备　　　　2.1 地表温度数据的获取　　关于地表温度数据的实地观测，选择在北京天坛公园内分别为纯草地、纯铺装、纯水泥路面、混合下垫面特征的多个样区中进行（图1），每一样区进行多点采样，主要使用2 台Raytek（Raytek Corporation，Germany）标准型红外测温仪进行，同步温度数据验证使用美国FLIR SYSTEMS 公司生产PM695 型热像仪，观测日期选择在2006 年8 月份，采用定点观测的方法。　　另外，以 Arc view 为平台，添加天坛遥感影像（具有北京地方坐标系）作为底图，结合周围的特征，人工确定样区位置并生成点状矢量主题，随后修改其属性表，附加测定的地表温度值。相关的坐标信息通过MapInfo Professional 的坐标提取工具，提出点坐标，然后转换为“xls”格式，最终生成符合要求的离散点文本数据格式。　　　　2.2 离点数据格式转换　　Surfer 绘制等值线的数据有特殊的格式要求，首先应将离散点等数据文件转换成Surfer认识的grd 文件格式，能够转换成grd 文件格式的离散点数据的格式一般为文本格式，它包含的数据项有资料点的位置坐标（如经纬度、WGS84 坐标、北京地方坐标）及属性项（如地表温度），其中属性项可以有多项，但资料点的位置坐标是唯一的[6,7]。　　2.3 离散数据插值　　当准备好数据文件后，关键的工作便是将离散点资料插值成规则的网格点数据，即可利用Surfer Automation 对象所属方法Grid data 来完成插值，具体可通过Grid data 方法内部参数的设定选择一种插值方法，Surfer 软件提供了许多数学插值模型，如距离倒数加权法、最小曲率法、克立格法等。利用这些数学模型，可以很方便地进行插值计算，其中较为精确的插值方法是克立格法（Kriging）[8]，该方法可以反映空间场的各向异性，并且充分利用数据点之间的空间相关性；可以自动识别样品点的空间分布，插值结果在某些情况下更加符合实际。　　　　　　3 利用Scripter 编辑器实现图形自动生成与显示　　　　Surfer 软件提供了脚本编辑器Scripter，与Surfer 在同一程序组下，提供了一系列编写、修改、调试脚本的功能，包括对象浏览器、User Dialog 编辑器、断点、立即窗口等。一个脚本由若干命令组成，脚本运行时命令被逐条执行，其命令格式类似于VB 编程语言，注释命令由英文单引号标明。　　　　3.1 Surfer 的ActiveX Automation 技术　　ActiveX Automation是微软公司建立的标准