一种基于地理信息应用城市独立坐标系参数获取方法技术

本发明专利技术请求保护一种基于地理信息应用城市独立坐标系参数获取方法，涉及利用高斯投影技术，计算该城市独立坐标系近似参数，从而实现地理信息成果相互转换的方法。本发明专利技术通过获取三个某城市独立坐标系与国家通用坐标系同名点成果，高斯投影正算公式进行变形，简化成一个关于该独立坐标系中央子午线L0的一个一元高次方程，并可利用迭代法解算出中央子午线L0，进而解算该独坐标系其他参数，包括：参考椭球长半径a，中央子午线L0，北方向偏移量X0和东方向偏移量Y0。这样便可以实现该独立坐标系与国家通用坐标系之间的同名点成果转换与应用。本发明专利技术可广泛用于地理空间数据建设，数字城市建设等与地理空间信息坐标转换应用。

【技术实现步骤摘要】

【技术保护点】
一种基于地理信息应用城市独立坐标系参数获取方法，其特征在于：1）利用GPS联测某城市独立坐标系四等以上的控制点，获得至少三个该城市独立坐标系同名点成果P1、P2、P3，在与测得所述城市独立坐标系同名点成果相同的三点上利用GPS与国家通用坐标系联测四等以上控制点，得到国家通用坐标系同名点成果P1′、P2′、P3′；2）将高斯投影正算公式简化为：x=X+NF(l)????????（3）y=NF1(l)?????????（4）式中：F(l)=12tcos2Bl2+124t(5-t2+9η2+4η4)cos4Bl4+1720t(61-58t2+t4)cos6Bl6F1(l)=cosBl+16(1-t2+η2)cos3Bl3+1120(5-18t2+t4+14η2-58η2t2)cos5Bl5其中，经差l=L?L0，以弧度为单位，L0为中央子午线，L为国家通用坐标系同名点成果的经度，B为国家通用坐标系同名点成果的纬度，x，y为城市独立坐标系同名点成果的平面坐标；3）城市独立坐标系参考椭球的扁率e与国家通用坐标系参考椭球的扁率相同，但是两种参考椭球的长半径a不同，在这种情况下城市独立坐标系的经度变化△L=0，纬度变化△B≈0，将P1和P1′、P2和P2′、P3和P3′任意两对分别带入式（3）、式（4）整理可得：x1=X0+X1+N1F(L1?L0)????????（5）y1=Y0+N1F1(L1?L0)??????????（6）x2=X0+X2+N2F(L2?L0)???????（7）y2=Y0+N2F1(L2?L0)?????????（9）x1?x2=X1?X2+N1F(L1?L0)?N2F(L2?L0)????（10）x1-x2y1-y2=X1-X2+N1F(L1-L0)-N2F(L2-L0)N1F1(L1-L0)-N2F1(L2-L0)---(11)式（11）是一个关于L0的一元高次方程，利用高斯迭代算法解出中央子午线L0；将L0代入式（10）求解出椭球长半径a将a，L0分别代入式（5）、式（6）求解出北方偏移量X0和东方偏移量Y0；4）根据步骤3）中获得的城市独立坐标系参数，利用P1、P2、P3和P1′、P2′、P3′的两套坐标系同名点成果，利用Bursa七参数转换模型，计算出坐标转换参数，所述转换参数包括3个平移参数，3个旋转参数和1个比例因子；Bursa七参数转换模型：XTYTZT=ΔXΔYΔZ+0-ZSYSZS0-XS-YSXS0ϵXϵYϵZ+mXSYSZS+XSYSZS---(12)式中，[△X?△Y?△Z]T为3个平移参数，[εX?εZ?εZ]T为3个旋转参数,m为比例因子，[XS?YS?ZS]T，[XT?YT?ZT]T分别为城市独立坐标系和国家通用坐标系同名点成果。FDA0000368761360000014.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：高翔，罗灵军，张泽烈，袁超，蒲德祥，张治清，胡艳，李莉，程宇翔，李仕锋，
申请(专利权)人：重庆市地理信息中心，
类型：发明
国别省市：