发明名称 三维化工园区建立方法
摘要 本发明公开了一种三维化工园区建立方法,该三维化工园区建立方法根据待扫描目标对象特征采用不同类型的三维激光扫描仪获取其三维信息,较好的解决了扫描精度与扫描速度之间的矛盾,且本发明的三维化工园区建立方法利用标靶的WGS-84标准坐标与三维激光扫描仪测得坐标之间的转换关系方便地将三维激光扫描仪获取的三维数据转换至与最新的2000坐标系可以通用的WGS-84标准坐标系,使得获得的最终数据可以直接用于工程、管理,同时可以和现行的GIS技术标准进行对接,具有较大的应用价值。
申请公布号 CN103791887A 申请公布日期 2014.05.14
申请号 CN201410033787.1 申请日期 2014.01.23
申请人 成都九洲电子信息系统股份有限公司 发明人 邢保振
分类号 G01C7/00(2006.01)I;G01B11/24(2006.01)I 主分类号 G01C7/00(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 周永宏
主权项 1.一种三维化工园区建立方法,其特征在于:包括以下步骤:S1、在化工园区内布设全站仪标靶:在待测建筑物内表面及内部设备、管线上安置标靶或棱镜,设布设的标靶的数量为n;S2、测定基准站坐标和标靶坐标:在化工园区内选择空间开阔的任意位置为基准站位置,并基于WGS-84标准坐标系测量基准站的坐标;基准站测定后,用全站仪测量得到每个标靶的WGS-84坐标,并将测量结果标记为{t<sub>i</sub>},i=1,2…n;S3、根据待扫描目标对象特征使用不同类型的三维激光扫描仪获取化工园区的三维数据{y<sub>kj</sub>},其中,k为大于零的自然数,表示待扫描目标对象的编号,j为大于零的自然数,表示所获取三维数据的个数;将获取的三维数据{y<sub>kj</sub>}中标靶的坐标标记为{z<sub>kh</sub>},h=1,2…m,5≤m≤n,表示标靶的编号;S4、坐标转换:S41、在全站仪测得的标靶的坐标{t<sub>i</sub>}中取坐标为{z<sub>kh</sub>}的标靶对应的坐标,并标记为{x<sub>kh</sub>};S42、将{z<sub>kh</sub>}与{x<sub>kh</sub>}进行配准,获取二者间的转换关系,设:x<sub>kh</sub>=Rz<sub>kh</sub>+T      (1)要配准这两个点集,目标就是要找到合适的旋转矩阵R和平移矩阵T使得待配准点集的欧氏距离和最小:<maths num="0001"><![CDATA[<math><mrow><msup><mi>&Sigma;</mi><mn>2</mn></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>kh</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>Rz</mi><mi>kh</mi></msub><mo>+</mo><mi>T</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>假设{x<sub>kh</sub>}与{z<sub>kh</sub>}的重心分别为p、q,则有:<maths num="0002"><![CDATA[<math><mrow><mi>p</mi><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>x</mi><mi>kh</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>q</mi><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>z</mi><mi>kh</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>定义p<sub>kh</sub>=x<sub>kh</sub>-p,q<sub>kh</sub>=z<sub>kh</sub>-q,则式(1)的总误差公式变为:<maths num="0004"><![CDATA[<math><mrow><msup><mi>&Sigma;</mi><mn>2</mn></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mi>kh</mi></msub><mo>-</mo><msub><mi>Rq</mi><mi>kh</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>计算矩阵H:<maths num="0005"><![CDATA[<math><mrow><mi>H</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>q</mi><mi>kh</mi></msub><msup><msub><mi>p</mi><mi>kh</mi></msub><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>对矩阵H作SVD分解:H=UDV<sup>T</sup>       (7)上式的解X=VU<sup>T</sup>即为所求旋转矩阵R;反带入求得平移矩阵TT=p-Rq       (8)S43、根据S42所得的转换关系,将三维激光扫描仪获取化工园区的三维数据{ykj}(y<sub>kj</sub>≠z<sub>kh</sub>)转换至WGS-84标准坐标系。
地址 610041 四川省成都市高新区天府大道中段765号