日照分析测绘中高层建筑物三维坐标获取方法.docVIP

日照分析测绘中高层建筑物三维坐标获取方法 　　摘要：对高层建筑高度测量方法进行探讨，提出一种简便可行的测量方法，并对其观测误差进行分析。 　　关键词：日照分析；高层建筑物；误差分析 　　1 引言 　　笔者从事城市测量工作，在平时的日照分析测绘时经常遇到这样的情况，某些高层建筑物及其附属物如电梯间太高，测绘他们的高度时要求的距离就会很远，就超出了平常的免棱镜仪器的距离要求，在日常的工作中，经过试验，笔者找到了一种方法，可以快速获取待测点的三维坐标。 　　2 外业测量步骤 　　首先使用RTK方法获取点A、B的三维坐标，作为已知点为A(XA,YA,HA),B(XB,YB,HB)。（AB应能通视，且都能与P通视） 　　在A点架设全站仪，用B点定向，瞄准待测点P，读出其水平角α，垂直角ι； 　　在B点架设全站仪，用A点定向，瞄准待测点P，读出其水平角β，垂直角κ。 　　 　　3 公式推导及内业计算 　　3.1利用前方交会方法计算P点平面坐标P(XP,YP) 　　 　　计算方法如下： 　　首先通过平面直角坐标反算，可获得AB边的坐标方位角αAB和边长SAB,由坐标方位角αAB和观测角α可推算出坐标方位角αAP，由正弦定理可得AP的边长SAP。由此根据平面直角坐标正算公式，即可求得待定点P的坐标，即 　　 　　当A、B、P按逆时针编号时，αAP=αAB-α，将其代入上式，得 　　 　　顾及XB -XA = ，YB-YA= ，则有 　　 （1） 　　由正弦定理可知： 　　 = 　　将上式带入（1）式，并整理可得 　　 （2） 　　式（2）即为前方交会计算公式，通常称为余切公式，是平面坐标计算的基本公式之一。 　　在此应指出：式（2）是在假定△ABP的点号A(已知点)、B(已知点)、 P(待定点)按逆时针编号的情况下推导出的。若A、B、P按顺时针编号，则相应的余切公式为： 　　 （3） 　　这样，待定点的平面坐标就有了。 　　3.2确定待定点的高程 　　首先由由正弦定理可以确定A、B到P点的水平距离SAP、SBP： 　　由 （式中γ=180°-（α+β））可知 　　SAP=SBP= 　　由三角高程公式可知， 　　HP=HA+ iA+SAP*tanι= HA+iA+ *tanι或 　　HP=HB+ iB+SBP*tanκ= HB+ iB+*tanκ 　　这两个计算公式可相互检核，亦可??平均值作为P点高程。 　　4 误差分析 　　由文献[2]可以知道，平面测角前方交会点位精度估算公式为 　　MP=（4） 　　式中m为测角中误差，一般可根据仪器标称精度和测回数来确定，是个定值。SAB长度可视为已知，没有误差。可以看出，测角交会点位精度主要取决于下式： 　　K= 　　所以，P点的点位精度与K成正比。 　　在实际工作中，使用的全站仪测角标称精度一般是5”,免棱镜测距距离范围在200米左右，所以不能用免棱镜时，SAP和SBP一般要大于200米，这里假设SAP=SBP=300米，SAB一般在100米左右，故可以估计此时的点位精度约为0.1米，满足平面点位误差为0.15米的要求。 　　由于短边三角高程测量精度较高，实际计算中可不考虑其误差。 　　5 实例验证 　　为了验证该方法的可靠性，我们使用Topcon 3005N 型免棱镜全站仪在实地进行了数据采集。同时，为了使计算的三维坐标与实际坐标有可比性，使SAP和SBP均没有超出其测距范围，这样就可以在采集角度ι、α的时候实测P点三维坐标，以此与计算坐标进行对比。 　　已知：XA= 47968.693YA=46134.019 HA=40.167;XB= 47918.445YB= 46184.587 HB= 40.329,可以计算出SAB=71.288 　　观测值 α ι iA β κ iB 　　 11°28′52″ 30°39′33″ 1.613 145°38′52″ 11°43′50″ 1.605 　　计算值 XP YP HP1 HP2 　　 48004.552 46127.159 63.422 63.425 　　实测值 XP YP HP 　　 48004.571 46127.152 63.402 　　差值 -0.019 0.007 0.020 0.023 　　通过以上数据可以判断，这种方法是可行的，但需注意以下几点： 　　1、A、B两点要尽量距离较远，且与待测点组成的三角形是锐角三角形 　　2、定向要准确，仪器高要准确量取 　　3、待测点要清晰可辨 　　参考文献 　　 [1] 潘正风，杨正尧，程效军等 . 数字测图原理与方法[M] . 武汉：武汉大学出版社，2009. 　　[2] 陈立文