摘要：本文结合航测工作实际生产案例，利用lisp程序语言，编写程序，实现对Z坐标进行整提调整功能。

关键词：lisp程序 Z坐标 调整

中图分类号： O4-34文献标识码：A 文章编号：

0 前言

与传统方法获取的地形图相比，航测立体测图获取的地形图最大的特点是每个节点均有真实的高程坐标，即Z值。目前对航测数据Z值的应用比较多，如进行三维建模、日照分析等。但是由于高程基准的不同，有时需要将不同基准下的航测立体测图高程进行进行转换，即同时对所有地物节点坐标的Z值增加一个数值或减少一个数值。如在ADS80航空摄影中，由于生产的需要，需要对高程系在WGS 84和1985国家高程系之间进行转换，这就需要对测图数据所有节点高程值进行处理，整体增加和减少一个固定值。

1 航测图的节点Z值坐标存放特点

目前航测图一般由立体测图软件的格式转换为Auto CAD软件的DWG格式。主要分为点、线、注记几种情形。以线状地物为例，在航测数据中，线状地物节点的高程值不完全一样，不是DWG数据常规LWPolyline的存储方式，而是3d Polyline形式存在。

航测立体测图数据线状地物节点坐标存储方式如下：

POLYLINELayer: 3100

Space: Model space

Handle = 3344

Open space

VERTEXLayer: 3100

Space: Model space

Handle = 3345

Space

at point, X=66442.174Y=4325491.568Z= 14.736

VERTEXLayer: 3100

Space: Model space

Handle = 3346

Space

at point, X=66444.294Y=4325492.224Z= 14.629

VERTEXLayer: 3100

Space: Model space

Handle = 3348

而一般DWG数据格式线状地物节点坐标存储方式如下：

LWPOLYLINELayer: 81100

Space: Model space

Color: BYLAYERLinetype: CONTINUOUS

Handle = 2F8

Open

Constant width 0.000

at pointX=35000.000Y=4304378.933Z=1.000

at pointX=34999.882Y=4304378.702Z=1.000

从上面可以看到，航测立体测图数据的DWG文件与一般的DWG线状数据主要有两点不一样：航测立体测图数据存储文件每个Z坐标不相同，即其存储的是该线状地物的每一个节点高程值，为三维多义线；而一般DWG数据Z值坐标相同，为一个值，是轻量多义线(LWPolyline)，即通常所说的二维半线；其次，航测立体测图数据是以节点为单位，每个节点为一个Handle，而一般DWG数据为一个实体(点、线、文本)为一个Handle。对于这两点区别是后期程序实现中需要注意的问题。而如何获取3d Polyline数据节点Z值坐标并对Z值坐标进行操作成为该程序的重点。

2 程序实现

程序设计基本思想为首先提取所有数据(点、线、文本)的X、Y、Z坐标。在提取过程中，注意在Auto CAD坐标定义中的X、Y坐标与测量中的X、Y是不同的，前者是数学坐标系，X表示横坐标，Y表示纵坐标。

在程序设计中，首先提取节点坐标组，定义为coor，在coor中依次包含X、Y、Z值，分别以coor_x、coor_y、coor_z定义。程序代码如下：

(setq ent_lst (entget nam))

(setq coor (cdr (assoc 10 ent_lst)))

(setq coor_x (car coor))

(setq coor_y (cadr coor))

(setq coor_z (caddr coor))

然后定义一个高差值h，作为改变高程的变量值，同时提示“请输入高差”，程序代码如下：

(setq h_new (getreal "\n请输入高差:"))

对调整高程后的坐标值定义一个新的坐标组，定义为coor_new，将坐标组代替原坐标组。程序代码如下：

(setq coor_new (list 10coor_x coor_y (+ h_new coor_z)))

(entmod (subst coor_new (assoc 10 ent_lst) ent_lst))

将以上更新坐标的程序作为一个子程序在后续循环中调用。

在实际作业过程由于有大量的点、线、注记，因此对上述过程需要对选中的地物逐个更改，需要对上述过程进行循环，具体程序代码如下：

(setq s0 (ssget ' ((0 . "polyline,INSERT,TEXT"))))

(setq h_new (getreal "\n请输入高差:"))

(if (= h_new nil)

(setq h_new 0)

);if

(repeat (sslength s0)

(setq nam0 (ssname s0 i))

(setq type_ent (cdr (assoc 0 (entget nam0))))