中纬全站仪GSI数据生成CASS9.1绘图坐标文件编程设计

测绘与空间地理信息

ＧＥＯＭＡＴＩＣＳ＆ＳＰＡＴＩＡＬＩＮＦＯＲＭＡＴＩＯＮＴＥＣＨＮＯＬＯＧＹ

Ｖｏｌ.４２ꎬＮｏ.４Ａｐｒ.ꎬ２０１９

中纬全站仪ＧＳＩ数据生成ＣＡＳＳ９.１

绘图坐标文件编程设计

李海灵

(福建省地质工程研究院ꎬ福建福州３５０００１)

摘要:利用全站仪测量碎部点是目前１∶５００数字地籍碎部测量主要的作业方法ꎮ为了节省外业设站、定向时

间ꎬ提高作业效率ꎬ可通过中纬全站仪采集测站点到碎部点的距离与角度等数据(下文称ＧＳＩ数据)ꎬ再用Ｖｉｓｕａｌ达到预期效果ꎮ

Ｂａｓｉｃ６.０编写程序ꎬ根据测站点坐标及ＧＳＩ数据生成ＣＡＳＳ９.１绘图坐标文件ꎮ相关测量实例证明了上述编程可关键词:ＧＳＩ数据ꎻ生成ꎻ坐标文件ꎻ编程设计

中图分类号:Ｐ２２１　　　文献标识码:Ａ　　　文章编号:１６７２－５８６７(２０１９)０４－０１９４－０３

ＰｒｏｇｒａｍｍｉｎｇｆｏｒＧＳＩＤａｔａｏｆＧｅｏｍａｘＥｌｅｃｔｒｏｎｉｃＴｏｔａｌＳｔａｔｉｏｎ

ＧｅｎｅｒａｔｉｎｇＤｒａｗｉｎｇＣｏｏｒｄｉｎａｔｅＦｉｌｅｏｆＣＡＳＳ９.１

(ＦｕｊｉａｎＩｎｓｔｉｔｕｔｅｏｆＧｅｏｌｏｇｉｃａｌＥｎｇｉｎｅｅｒｉｎｇꎬＦｕｚｈｏｕ３５０００１ꎬＣｈｉｎａ)

Ａｂｓｔｒａｃｔ:ＵｓｉｎｇＥｌｅｃｔｒｏｎｉｃＴｏｔａｌＳｔａｔｉｏｎｔｏｍｅａｓｕｒｅｂｒｏｋｅｎｐｏｉｎｔｓｉｓｔｈｅｍａｉｎｍｅｔｈｏｄｏｆｔｈｅ１∶５００ＤｉｇｉｔａｌＣａｄａｓｔｒａｌＤｅｔａｉｌＳｕｒｖｅｙａｔｐｒｅｓｅｎｔ.Ｉｎｏｒｄｅｒｔｏｓａｖｅｔｈｅｔｉｍｅｏｆｂｕｉｌｄｉｎｇｓｔａｔｉｏｎｓａｎｄｏｒｉｅｎｔｉｎｇꎬｉｍｐｒｏｖｅｗｏｒｋｅｆｆｉｃｉｅｎｃｙꎬｉｔｉｓａｇｏｏｄｗａｙｔｏｃｏｌｌｅｃｔｔｈｅｄａｔａｏｆｄｉｓ￣ｄａｔａ.Ｓｏｍｅｒｅｌａｔｅｄｓｕｒｖｅｙｐｒａｃｔｉｃａｌｅｘａｍｐｌｅｓｐｒｏｏｆｔｈａｔｔｈｅａｆｏｒｅｍｅｎｔｉｏｎｅｄｐｒｏｇｒａｍｍｉｎｇｈａｓａｃｈｉｅｖｅｄｔｈｅｅｘｐｅｃｔｅｄｅｆｆｅｃｔ.Ｋｅｙｗｏｒｄｓ:ＧＳＩｄａｔａꎻｇｅｎｅｒａｔｅꎻｃｏｏｒｄｉｎａｔｅｆｉｌｅꎻｐｒｏｇｒａｍｍｉｎｇ

ｔａｎｃｅｓａｎｄａｎｇｌｅｓｆｒｏｍａｓｔａｔｉｏｎｔｏａｂｒｏｋｅｎｐｏｉｎｔａｎｄｏｔｈｅｒｒｅｌａｔｅｄｄａｔａ(ｈｅｒｅｉｎａｆｔｅｒｒｅｆｅｒｒｅｄｔｏａｓＧＳＩｄａｔａ)ｂｙＧｅｏｍａｘＥｌｅｃｔｒｏｎｉｃＴｏｔａｌＳｔａｔｉｏｎａｎｄｔｈｅｎｐｒｏｇｒａｍｂｙＶｉｓｕａｌＢａｓｉｃ６.０ꎬｔｏｇｅｎｅｒａｔｅＣＡＳＳ９.１ｄｒａｗｉｎｇｃｏｏｒｄｉｎａｔｅｆｉｌｅｓｂａｓｅｄｏｎｓｉｔｅｃｏｏｒｄｉｎａｔｅｓａｎｄＧＳＩ

ＬＩＨａｉｌｉｎｇ

０　引　言

近３年来ꎬ福建省陆续开展农村不动产测量ꎬ１∶５００数字地籍碎部测量仍主要采用全站仪作业[１]ꎮ全站仪测量碎部点有两种方式:一是使用仪器上装载的程序直接测定碎部点坐标(下文称直接坐标法)ꎻ二是外业采集测站点至碎部点的距离、角度等数据ꎬ内业再由程序计算生成其坐标数据(下文称边角法)ꎮ边角法可以免去外业现场给定测站及定向点坐标ꎬ设站、定向更快捷ꎬ节省时间ꎮ此方法尤其适用分散的农村居民点测量ꎬ实现图根控制(采用ＧＮＳＳＲＴＫ测量)与碎部测量适时进行ꎬ充分利用外业时间ꎮ另外ꎬ如果测站点坐标有更新ꎬ内业就可通过程序适时计算更新碎部点坐标ꎬ提高内业效率ꎮ因此ꎬ边角法比直接坐标法更具优势ꎮ南方ＣＡＳＳ是如今测绘行业内用户量最大、使用最方便快捷的主流地形、地籍成图软件ꎬＣＡＳＳ９.１是目前较新版本ꎬ其符号库符合２００７年版

收稿日期:２０１８－０６－２０

图式(ＧＢ/Ｔ２０２５７.１－２００７)ꎮ但ＣＡＳＳ９.１目前不能直接读取中纬全站仪ＧＳＩ数据ꎬ业内也没有普及ＧＳＩ数据生成ＣＡＳＳ９.１绘图坐标文件的程序ꎬ这就需要编写一个程序解决上述问题ꎮ通过ＶｉｓｕａｌＢａｓｉｃ６.０编程实现上述数据衔接应运而生ꎮ

１　ＧＳＩ数据及ＣＡＳＳ９.１绘图坐标文件解析

中纬全站仪采用ＧＳＩ(串行接口)数据格式结构ꎬ该数据有两种存储格式:８位字符ＧＳＩ－８格式和１６位字符ＧＳＩ－１６格式ꎮ本文针对处理的是ＧＳＩ－８格式ꎮ文件一行记录实例:１１０００３＋０００００１３６２１.３２４＋３５９５９５９０２２.３２４＋０９００６５９０３１...０＋０００６４２５１５１....＋０００８＋０００８７...０＋“５１”表示常数ꎻ“８７”表示棱镜高ꎮ

００００１７００８８....＋００００００００ꎮ其中ꎬ字索引“１１”表示点号ꎻ

ＣＡＳＳ９.１绘图坐标文件(下文简称坐标文件)指的是

“２１”表示水平度盘ꎻ“２２”表示垂直度盘ꎻ“３１”表示斜距ꎻ

作者简介:李海灵(１９８３－)ꎬ男ꎬ福建武平人ꎬ工程师ꎬ学士ꎬ主要从事工程测量方面的应用研究工作ꎮ

第４期

李海灵:中纬全站仪ＧＳＩ数据生成ＣＡＳＳ９.１绘图坐标文件编程设计

１９５

含简编码格式的坐标数据文件ꎬ扩展名是“ｄａｔ”ꎮ本文所述生成的坐标文件在原来文件记录“点号”字段上作改进ꎬ将测点点号数字编码(外业所用的)化ꎬ而非流水号ꎬ有助于内业成图中识别当前碎部点是何地物点ꎬ以便作相应图形编辑ꎮ文件一行记录实例(线状地物连线第一点):２０４１ꎬＦ８ꎬ５３３９２９.６９９ꎬ２７６４８５６.１２３ꎬ０ꎮ文件内每一行代表一个点ꎬ各字段分别为:点号(数字编码＋ꎮ地物

测点序号)ꎬ简编码ꎬＹ(东)坐标ꎬＸ(北)坐标ꎬ高程[２]

或线状地物不连线点的点号中不含测点序号ꎬ仅为该点

的外业数字编码ꎮ

２　坐标文件生成程序设计

坐标文件生成主要包含碎部点坐标计算及编码转换ꎮ坐标文件生成流程如图１所示ꎮ

图１　坐标文件生成流程

Ｆｉｇ.１　Ｃｏｏｒｄｉｎａｔｅｆｉｌｅｇｅｎｅｒａｔｉｏｎｐｒｏｃｅｓｓ

２.１　碎部点坐标计算

碎部点坐标计算主要包括碎部点与测站点间的平面坐标增量计算及高差计算ꎮ求得平面坐标增量及高差后ꎬ分别与对应测站点的平面坐标、高程值求和ꎬ即获得待求碎部点坐标ꎮ当碎部点棱镜高为０时ꎬ本程序认为其高程无效“点高程展高程点ꎬ直接将其高程赋值为０ꎮ运用ＣＡＳＳ９.１软件的)ꎬ快捷排除了”工具时ꎬ将不展出高程值为碎部点无效的高０程的点ꎬ提(高如阳台角图面整

２.１.１　洁度ꎮ

依据数学公式平面坐标增量计算

ｓｉｎαＡＢꎬ其中ｙ坐标增量ꎬΔｘＡＢ、Δꎬ可得、水平距离ｙＡＢ:Δ、ＳｘＡＢ＝、αＡＢＳ、坐标方位角分别为直线􀅰ｃｏｓαＡＢꎬΔ[３]ꎮＡＢｙＡＢ＝因为的ｘＳＧＳＩ坐􀅰

标增量、数据包含的字段中只有各碎部点到测站点的斜距ꎬ没有水平距离ꎬ因此ꎬ应对上述平面坐标增量计算式作改进ꎬ

详见①Ｖｉｓｕａｌ计算Ｂａｓｉｃｘ坐标６.０(北坐标定义其函数过程)增量:

(子程序)的代码ꎮ

ｓ１为观测斜距ＰｒｉｖａｔｅＦｕｎｃｔｉｏｎꎬ单位:ｊｓｚｌ＿ｘ(ｓ１ꎬ米ꎻｂ为天顶距ｂꎬａ＿ｄｘꎬꎬ单位ｓｐｊｊ):度、Ａｓ分Ｄｏｕｂｌｅ、秒ꎻ'

ｄｘ为定向边方位角ꎬ单位:度ꎻｓｐｊｊ为水平夹角ꎬ单位:度ａ＿、分、秒

Ｄｉｍｄｅｌ＿ｘꎬｈｄ＿ｂꎬａ＿ｄｑｂꎬｐｉ

ｐｉｈｄ＿ｂ＝３.１４１５９２６５３５８９７９

ａ＿ｄｑｂ＝＝ｚｈｗｄ(ｂ)􀅰(ｐｉ(ａ＿ｄｘ＋ｓｐｊｊ)/􀅰(１８０)

计算

ｐｉ/１８０)'待求边方位角ｄｅｌ＿ｊｓｚｌ＿ｘ＝Ｒｏｕｎｄ(ｓ１􀅰Ｓｉｎ(ｈｄ＿ｂ)􀅰Ｃｏｓ(ａ＿ｄｑｂ)ꎬＥｎｄｘ＝ｄｅｌ＿ｘ３)②Ｐｒｉｖａｔｅ计算Ｆｕｎｃｔｉｏｎ

Ｆｕｎｃｔｉｏｎｙ坐标(东坐标ｊｓｚｌ＿ｙ(ｓ１ꎬ)增量ｂꎬ:Ｄｉｍａ＿ｄｘꎬｓｐｊｊ)ＡｓＤｏｕｂｌｅｐｉｈｄ＿ｂ＝ｄｅｌ＿３.１４１ｙꎬ＝ｚｈｗｄ(ｂ)􀅰(ｐｉ５９２ｈｄ＿ｂꎬ６５３ａ＿ｄｑｂꎬ５８９ｐｉ/７９

计算

ａ＿ｄｑｂ＝(ａ＿ｄｘ＋ｓｐｊｊ)􀅰(１８０)

ｐｉ/１８０)'待求边方位角ｄｅｌ＿ｊｓｚｌ＿ｙ＝Ｒｏｕｎｄ(ｓ１􀅰Ｓｉｎ(ｈｄ＿ｂ)􀅰Ｓｉｎ(ａ＿ｄｑｂ)ꎬ２.１.２　Ｅｎｄｙ高差计算

Ｆｕｎｃｔｉｏｎ

＝ｄｅｌ＿ｙ３)依据公式ｈ＝Ｓ􀅰ｓｉｎα＋(Ｓ􀅰ｃｏｓα)２为经过各种改正后的斜距ꎻα为观测垂直角/２Ｒ＋ꎻｉＲ－ｖ为地球平(式中ꎬＳ

均曲率半径ꎻｉ为设站仪器高ꎻｖ为反射棱镜高)[４]计算高差ꎮＶＢＰｒｉｖａｔｅ子程序代码如下:

ＤｉｍＦｕｎｃｔｉｏｎｊｓｚｌ＿ｈ(ｓ１ꎬｂꎬｙｇ１ꎬｊｇ)ＡｓＤｏｕｂｌｅｐｉｈｄ＿ｂ＝ｄｅｌ＿ｈꎬ３.１４１５９２ｈｄ＿ｂꎬ６５３ｐｉ

５８９ｄｅｌ＿ｈ＝＝ｚｈｗｄ(ｂ)􀅰(ｐｉＲｏｕｎｄ(ｓ１􀅰Ｃｏｓ(/７９１８０)

(Ｓｉｎ(ｈｄ＿ｂ)))ｈｄ＿ｂ)＋ｙｇ１－ｊｇ＋(ｓ１􀅰

ｊｓｚｌ＿ｈ＾２２.１.３　ＥｎｄＦｕｎｃｔｉｏｎ

＝ｄｅｌ＿ｈ

/(２􀅰６３６９０００)ꎬ３)

起算数据检测是保证输出坐标文件数据准确性很重

碎部点坐标起算数据检测设计

要的一个步骤ꎮ设计思路是根据控制点坐标反算出定向边的平距ＰＪ１、求出测站点与定向点的高差ＧＣ１ꎬ再由外业观测定向点的数据计算定向边的平距ＰＪ２、测站点与定向点的高差ＧＣ２ꎮ当∣ＰＪ２－ＰＪ１∣>５ｃｍ或∣ＧＣ２－ＧＣ１∣>８ｃｍ[５]行代码ꎬ而退出当前过程ꎬ程序弹出起算数据核对提示(Ｅｘｉｔｓｕｂ)ꎮ因此ꎬ并不再往下执ꎬ外业定向后ꎬ

须将定向点作为一碎部点观测并记录ꎮ内业通过本程序生成坐标文件时也将生成一个定向检测文件ꎬ反馈起算数据检测情况ꎮ

２.２　编码转换

编码转换就是将外业测点编码(本文中指ＥＰＳＷ测图编码)转换为ＣＡＳＳ９.１房屋ＥＰＳＷ、９２１表示未加固陡坎测图编码由３软件能识别的相应简编码位数字组成ꎮ其优点是纯数字ꎬ例如:２０４ꎮ

ꎬ观测记录表示混合

时不用进行输入切换ꎬ能便捷输入仪器ꎬ节省外业时间ꎮ因此ꎬ笔者力荐外业使用纯数字编码ꎮ其缺点是初学者不容易记住ＣＡＳＳ９.１ꎮ

简编码有１—３位ꎬ第一位是英文字母ꎬ大小

１９６

测绘与空间地理信息　

　　　　　　　　　　　　　　　２０１９年

写等价ꎬ后面是范围为０—９９的数字ꎬ无意义的０可以省略ꎬ例如:Ａ和Ａ００等价、Ｆ１和Ｆ０１等价ꎮ字母与数字组合的编码在外业不便输入仪器ꎬ影响作业效率ꎮ除了描述实体属性的简编码外ꎬ还有描述线状地物连接关系的连接码ꎬ常用的如“＋”ꎬ意指本点与上一点相连ꎬ连线依测点顺序进行ꎮ

因此ꎬ编码转换的关键是通过ＣＡＳＳ９.１的简编码定义文件(ＪＣＯＤＥ.ＤＥＦ)及符号定义文件(ＷＯＲＫ.ＤＥＦ)找出与外业测点编码相对应的ＣＡＳＳ９.１简编码ꎮ例如:２０４、９２１对应的ＣＡＳＳ９.１简编码为Ｆ８、Ｋ０ꎮ

３　程序界面设计

程序界面内容主要有７项:１)控制点文件(即ＫＮＷ文件ＤＡＴ)定向类型单选按钮文件所在路径指定“打开”按钮ꎻ２)ＧＳＩꎻ７)ꎻ４)“支站坐标提取并转生成文件所在路径指定”按钮ꎻ５)“退出ꎻ３)ＫＮＷ”按钮生成的格式复

ꎻ６)选框ꎮ

在界面上设计ＤＡＴ文件“打开”按钮是为了程序使用者能在第一时间浏览生成结果ꎮ定向类型包含“置零”及“未置零”两项ꎮ“置零”指外业照准定向点时将水平度盘设为零ꎬ此项被选中时ꎬ程序直接采用ＧＳＩ中水平度盘读数进行计算碎部点方位角ꎮ“未置零”选项是为了外业定向忘记将水平度盘设为零的情况而考虑的ꎬ此项被选中时ꎬ程序则将每个碎部点的水平度盘读数减去相应定向点的水平度盘读数ꎬ获得水平夹角ꎬ再计算相应碎部点方位角ꎮ支站坐标提取并转ＫＮＷ格式(ＮＡＳＥＷ软件输出的控制点成果文件格式)选项是为了通过程序提取支站坐标ꎬ并以ＫＮＷ文件存储ꎬ便于后续使用支站坐标ꎮ４　坐标文件生成实例

５００２０１７年３月至４月ꎬ笔者在福建惠安两村庄进行１∶

据来源于此测区数字地籍测量ꎬ以下中纬全站仪ＧＳＩ数据及控制点数１)ＧＳＩ由于篇幅关系数据

ꎮ

ꎬ选取了其中一天外业采集的部分ＧＳＩ

数据作说明ꎬ如图２所示ꎮ

图２　ＧＳＩ数据Ｆｉｇ.２　ＧＳＩｄａｔａ

Ｇ１８２)控制点坐标文件数据

Ｇ１９１０７２７６４９１７.０１８Ｇ２０１０７１０７２７６４９４２.６２５５３３８１０.８８１２７６４９８６.９３７５３３７７６.０６９２６.５０１１５３３７２６.３８４２６.３０１２６.５５１１１

Ｇ２１Ｇ２２１０７Ｇ２３１０７２７６５０４２.９２５５３３６７０.２４９２７.２１８１Ｇ２５１０７２７６５０８４.２８９２７６５１１６.３４０５３３６３９.６８１５３３６８６.０１６２８.４７５２７.１７６１Ｇ２６１０７Ｇ２７１０７２７６５１７４.７４８５３３７６４.５５２２４.１９７１１Ｇ２８１０７２７６５１２２.５８７２７６５０７２.３８８５３３８０９.２９７５３３７８４.７２５２２.９６８２５.０９１１１

Ｇ３０Ｇ２９１０７１０７２７６５０４９.５０８２７６５０１４.３６５５３３８１２.０７６５３３７８６.６５４２４.５５８２５.７０６１１

３)１０７２７６４９９７.９８５５３３８０９.３０９２５.４４５１运行本程序通过程序生成的坐标文件数据

ꎬ获得与上文列出的ＧＳＩ数据对应坐标文件数据如下１３６ꎬＳ１９ꎬ５３３６３９.６８０ꎬ２７６５０８４.２８８ꎬ２８.４７７:

２０４ꎬꎬ５３３６９１.６１３ꎬ２７６５１０３.５６６ꎬ０

２２０ꎬꎬ５３３７０１.７３６ꎬ２７６５１２５.６３０ꎬ０２２０ꎬꎬ５３３７２０.７５４ꎬ２７６５１３９.４５３ꎬ０２０４１ꎬＦ８ꎬ５３３６８８.２７０ꎬ２７６５１２５.０２８ꎬ０２０４２ꎬ１３６ꎬＳ１９ꎬ５３３６３９.６８３ꎬ２７６５０８４.２８６ꎬ２８.４６４＋ꎬ５３３６９１.５０３ꎬ２７６５１２６.８７４ꎬ０２５１１ꎬＭ４ꎬ５３３６５４.６４１ꎬ２７６５０７０.３２４ꎬ０２５１２ꎬ６１１ꎬꎬ５３３６６２.５０２ꎬ２７６５０５２.４９９ꎬ０＋ꎬ５３３６５６.５１７ꎬ２７６５０７１.５３１ꎬ０６０３ꎬꎬ５３３６６２.８５９ꎬ２７６５０５２.３６１ꎬ０２００１ꎬＦ０ꎬ５３３６７３.９５２ꎬ２７６５０４７.６５９ꎬ０２００２ꎬ２４８ꎬꎬ５３３６６９.１６７ꎬ２７６５０５２.９６４ꎬ０＋ꎬ５３３６６９.３５２ꎬ２７６５０５３.０８７ꎬ０

２４８１ꎬＭ１ꎬ５３３６６７.９５０ꎬ２７６５０５５.２０５ꎬ０２４８２ꎬＺ１ꎬＣ１１ꎬ５３３６４８.２０１ꎬ２７６５０７４.５４８ꎬ２８.３１３＋ꎬ５３３６６７.５７７ꎬ２７６５０５４.８５１ꎬ０２２９ꎬꎬ５３３６７３.１８９ꎬ２７６５０４７.０３８ꎬ０２５１ꎬꎬ５３３６６８.５４８ꎬ２７６５０５２.２５０ꎬ０

２２９ꎬꎬ５３３６６５.８０３ꎬ２７６５０５５.９１４ꎬ０２００ꎬꎬ５３３６６６.２０２ꎬ２７６５０５６.６１４ꎬ０

１１４１ꎬＷ０ꎬ５３３６６０.４７９ꎬ２７６５０６３.６２４ꎬ０１１４２ꎬ２２９ꎬꎬ５３３６８３.０９９ꎬ２７６５０５５.３４６ꎬ０

＋ꎬ５３３６６０.９７２ꎬ２７６５０６３.６６０ꎬ０Ｚ２ꎬＣ１１ꎬ５３３６５７.４６４ꎬ２７６５０６０.３２４ꎬ２７.８９２５４７ꎬꎬ５３３６６５.６０２ꎬ２７６５０５５.２４２ꎬ０５４７ꎬꎬ５３３６６２.８２０ꎬ２７６５０５１.７８２ꎬ０２０４１ꎬＦ８ꎬ５３３６９７.３７３ꎬ２７６５０４２.８０３ꎬ０２０４２ꎬ２０４１ꎬＦ８ꎬ５３３６９７.６３６ꎬ２７６５０４２.２１５ꎬ０＋ꎬ５３３６９７.７４５ꎬ２７６５０４２.４２８ꎬ０２０４２ꎬ２０４３ꎬ＋ꎬ５３３６９４.１４０ꎬ２７６５０３８.６９２ꎬ０２０４４ꎬ＋４)＋ꎬ５３３６９６.６７８ꎬ２７６５０３６.１６９ꎬ０ꎬ５３３６９４.３７３ꎬ２７６５０３３.８４６ꎬ０

由于篇幅关系定向检测文件数据

ꎬ从本程序生成的该实例定向检测文件中ꎬ选取部分数据以示检测情况ꎬ差:０.００２Ｇ２３－ｍꎮ

Ｇ２２边长检测较差:０.００１具体如下ｍꎬＧ２２:

高程检测较(下转第１９９页)

第４期

杨爱玲等:洪水风险图编制技术研究

１９９

４.４　洪水影响和损失分析

料和调查洪痕资料来率定参数文件ꎮ控制下游出口断面水位相差在０.２０ｍ以内ꎬ洪峰流量相差在５％以内率定糙率参数成果ꎮ

二维模型利用土地利用图、航拍图ꎬ并结合现场调查ꎬ考虑分区内的地形、地貌、植被状况ꎬ根据«洪水风险图编制技术细则»«水力计算手册»等资料ꎬ对不同下垫面赋予不同的糙率值ꎬ并将结果以坐标加糙率的文本格式导入模型ꎬ通过插值得到各网格糙率值ꎮ

根据防洪保护区、蓄滞洪区、城市各方案洪水分析得到的最大淹没范围、最大淹没水深、淹没历时等要素ꎬ结合淹没区社会经济情况ꎬ综合分析评估洪水影响程度(包括淹没范围内、不同淹没水深区域内的人口、资产统计分析等ꎬ并评估洪水损失)ꎮ洪水影响分析包括淹没区的人口与社会经济指标计算及洪水灾害损失评估ꎮ对洪水风险区进行影响分析ꎬ主要是统计不同量级洪水各级水深淹没区域内的经济和社会指标ꎬ从而在一定程度上反映洪水的危害程度ꎮ

影响分析主要考虑淹没范围内不同淹没水深区域(外洪<０.５ｍꎬ０.５—１.０ｍꎬ１.０—２.０ｍꎬ２.０—３.０ｍꎬ>３ｍꎬ内涝０.０５—０.３ｍꎬ０.３—０.５ｍꎬ０.５—１.０ｍꎬ１.０—２.０ｍꎬ２.０—３.０ｍꎬ>３ｍ)内的受淹面积、受淹耕地面积、受灾人口总数、受影响交通线路及重要设施、受影响行政区域以及ＧＤＰ等指标ꎮ

持ꎮ洪水风险图在行政主管部门制定防洪规划、拟订防洪减灾非工程措施、部署防汛抢险应急保障措施等工作中起到重要的作用ꎬ在土地规划管理和便于水资源管理部门依法行政方面起到了很好的技术保障作用ꎬ并且可以增强全民的防洪减灾意识ꎮ

通过项目实施ꎬ逐步规范和统一洪水风险图编制技术体系和技术要求ꎬ提高洪水风险图编制技术水平和技术能力ꎻ建立洪水风险图应用管理制度ꎬ初步确立洪水风险图的法律效力ꎬ使洪水风险管理的实施能够有据可循ꎬ使加大洪水风险和防洪压力的各种社会行为得到控制ꎻ推进洪水风险图在防汛指挥、洪水风险区土地管理、增强全民水患意识和洪水影响评价等领域的实际应用ꎬ为防洪规划与应急决策、灾情评估、居民避险、土地利用开发、灾害保险、公共减灾对策以及灾害教育与宣传等工作提供技术支撑ꎬ进一步提高相关部门的洪水风险管理能力和水平ꎮ

参考文献:

[１]　中国水利水电科学研究院.ＳＬ４８３—２０１０洪水风险图编

制导则[Ｓ].北京:中国水利水电出版社ꎬ２０１１.

[２]　张硕辅ꎬ薛光达ꎬ曾务书.湖南省洪水风险分析防洪风险

图编制及其应用[Ｊ].湖南水利水电ꎬ２００１(１):２２－２４.[３]　贾俊杰ꎬ于得万ꎬ丁曼.松花江防洪保护区洪水风险图编

制研究[Ｊ].东北水利水电ꎬ２０１５ꎬ３３(１１):５３－５５.

[４]　孙玉贤ꎬ郭晓萌ꎬ张玉杰.河南省２０１３年度洪水风险图

编制过程与经验总结[Ｊ].河南水利与南水北调ꎬ２０１５(１４):１３０－１３１.

[５]　崔洪波ꎬ周墨ꎬ李井杰.基于ＧＩＳ空间分析的应急洪水风

险图制作[Ｊ].测绘与空间地理信息ꎬ２０１５ꎬ３８(１１):１０６－１０８.

５　结束语

利用洪水风险图编制项目中得到的成果数据ꎬ可在实时洪水风险图淹没展示等方面展开成果的应用ꎬ并且利用洪水风险图中的基础地形数据ꎬ水文数据、结合实际影像提取数据ꎬ在水灾过程中开展灾前、灾中、灾后的损失评估ꎬ结合防汛预案的对策ꎬ启动应急调度指挥提供支(上接第１９６页)

　　Ｇ２１－Ｇ２２边长检测较差:－０.００３ｍꎬＧ２２高程检测较差:－０.０１１ｍꎮ

Ｚ１－Ｇ２１边长检测较差:０.０００ｍꎬＧ２１高程检测较差:－０.００７ｍꎮ

５)ＣＡＳＳ９.１展绘图

运行ＣＡＳＳ９.１软件ꎬ将上述生成的坐标文件展绘成图(尚未编辑的初始图)ꎬ如图３所示ꎮ

[编辑:任亚茹]

的坐标文件ꎬ便可实现自动绘图ꎬ获得初始图ꎬ接着进行图形编辑ꎬ这样也节省内业成图时间ꎮ全站仪品牌很多ꎬ每个品牌的测量数据格式也不一样ꎬ中纬全站仪ＧＳＩ数据格式具有一定的代表性ꎬ其生成坐标文件的编程方式对其他品牌的全站仪数据生成坐标文件具有一定的参考意义ꎮ

参考文献:

[１]　詹长根ꎬ唐祥云ꎬ刘丽.地籍测量学[Ｍ].第２版.武汉:武

汉大学出版社ꎬ２００５.

[２]　宋志辉ꎬ王建军ꎬ唐弢ꎬ等.测图数据转换为ＣＡＳＳ数据成

图的研究及实现[Ｊ].工程地球物理学报ꎬ２０１１ꎬ８(４):５１０－５１２.

[３]　李金如ꎬ毛启麟ꎬ章学信.地形测量学[Ｍ].第２版.北京:

地质出版社ꎬ１９９８.

[４]　国家测绘局测绘标准化研究所.ＧＢ/Ｔ１２８９８－２００９国家

三、四等水准测量规范[Ｓ].北京:中国标准出版社ꎬ２００９.

[５]　朱留华ꎬ郑凌志ꎬ冯文利ꎬ等.ＴＤ/Ｔ１００１－２０１２地籍调查

规程[Ｓ].北京:中国标准出版社ꎬ２０１２.

图３　ＣＡＳＳ９.１展绘图(部分)

Ｆｉｇ.３　ＣＡＳＳ９.１ｅｘｈｉｂｉｔｉｏｎｄｒａｗｉｎｇ(ｐａｒｔｉａｌ)

５　结束语

通过上述实例证明ꎬ给本程序提供控制点坐标数据及外业ＧＳＩ数据ꎬ就可生成坐标文件ꎮ外业应用有编码采点ꎬ内业使用ＣＡＳＳ９.１软件的“简码识别”工具打开生成

[编辑:张　曦]