大地坐标转换成施工坐标公式.doc

大地(高斯平面)坐标系工程坐标系转换大地坐标系--—〉工程坐标系========================待转换点为P,大地坐标为：Xp、Yp工程坐标系原点o：大地坐标：Xo、Yo工程坐标：xo、yo工程坐标系x轴之大地方位角：adX=Xp-XodY=Yp—YoP点转换后之工程坐标为xp、yp:xp=dX＊COS(a）+dY＊SIN（a）+xoyp=—dX*SIN(a)+dY＊COS(a）+yo工程坐标系———>大地坐标系========================待转换点为P,工程坐标为：xp、yp工程坐标系原点o:大地坐标：Xo、Yo工程坐标:xo、yo工程坐标系x轴之大地方位角：adx=xp—xody=yp—yoP点转换后之工程坐标为xp、yp:xp=Xo+dx＊COS（a)-dy＊SIN（a）yp=Yo+dx＊SIN（a)+dy＊COS（a)坐标方位角计算程序置镜点坐标：ZXZY后视点坐标：HXHY方位角：W两点间距离：SLb10←{A，B，C，D}←A〝ZX=〞:B〝ZY=〞：C〝HX=〞：D〝HY=〞：W=tg1（（D-B）÷(C-A)）:(D—B)〉0=〉（C—A）〉0=〉W=W：∟∟(D-B）>0=〉（C-A)<0=>W=W+180：∟∟(D—B)〈0=>（C—A)<0=〉W=W+180:∟∟（D-B）〈0=〉(C—A）〉0=>W=360+W∟∟W=W◢S=√（(D—B)2+（C—A)2)◢Goto0←CASIOfx－4500p坐标计算程序根据坐标计算方位角W＝W＋360△W：“ALF(1～2）＝"L1A“X1＝”：B“Y1＝”:Pol（C“X2”－A，D“Y2"－B：“S＝”▲W＜0直线段坐标计算L1X“X（0）”：Y“Y（0)”：S“S（0)”:A“ALF”L2Lb12L3｛L｝：L“LX"L4M“X（Z)”＝X＋(L－S）cosA▲L5N“Y（Z)"＝Y＋(L－S)sinA▲L6｛B｝：B“B（L）":Q“Q"L7O“X(L）”＝M＋Bcos（A＋Q＋180)▲L8P“Y（L)”＝N＋Bsin(A＋Q＋180）▲L9{C｝：C“B（R）”L10U“X(R)”＝M＋Ccos（A＋Q）▲L11V“Y（R)"＝N＋Csin（A＋Q）▲L12Goto2园曲线段坐标计算L1S“S(0）-Km”：X“X(0)"：Y“Y（0）”：A“ALF”：R“R”：K“K（L＝1,R＝2）”L2Lb12L3{L｝:L“L（X)”L4V＝180／π×（L－S）／R：W＝V／2L5C＝A＋（—1)K×W：D＝2RsinW：F＝A＋（—1）K×VL6M“X（Z)”＝X＋DcosC▲L7N“Y（Z)”＝Y＋DsinC▲L8｛E}：E“B（L)":Q“Q"L9O“X(L)”＝M＋Ecos（F＋Q＋180)▲L10P“Y（L）”＝N＋Esin（F＋Q＋180）▲L11｛G｝:G“B（R)”L12T“X(R）”＝M＋Gcos（F＋Q）▲L13U“Y（R）”＝N＋Gsin(F＋Q)▲L14Goto2正向缓和曲线段坐标计算L1S“ZH—Km”：X“X（ZH)”：Y“Y（ZH）"：A“ALF”：R“R"：H“LS”：K“K(L＝1，R＝2）”L2Lb12L3｛L｝：L“L（X）”L4D＝30（L－S）2／π／R／H:C＝L－S－（L－S）5／90／（R×H)2：B＝A＋D(-1）K:E＝A＋3D（-1）KL5U“X(Z)”＝X＋CcosB▲L6V“Y（Z）”＝Y＋CsinB▲L7{G}：G“B(L)”:Q“Q”L8F“X(L）”＝U＋Gcos（E＋Q＋180）▲L9I“Y（L）”＝V＋Gsin(E＋Q＋180)▲L10｛J｝:J“B（R）”L11M“X（R）”＝U＋Jcos(E＋Q）▲L12N“Y（R)”＝V＋jsin（E＋Q）▲L13Goto2卵形曲线坐标计算X＝1，D＝2）”L1S“Km—YH":E“X（YH）"：F“Y（YH）”：G“ALF”：B“R1”:D“A"：K“K（L＝1，R＝2)”：Q“R1-R2L2Lb12L3{Z｝：Z“L（X)"L4J“L1”＝D2／B:R“RP”＝D2B／(D2＋（-1)Q(Z－S)B）:L“LP”＝D2／RL5M＝（L－J)－(L5－J5)／40／D4＋(L9－J9)／3456／D8L6N＝(L3－J3）／6／D2－（L7－J7）／336／D6＋(L11－J11）／42240／D10L7T＝G－(—1）Q(-1)K×J2×90／D2／πL8X“X（Z）”＝E＋(—1）QMcosT－（—1)KNsinT▲L9Y“Y（Z）”＝F＋（—1)QMsinT＋(—1)KNcosT▲L10A“ALF（P）"＝G＋(-1)K(Z－S)×90×（1／B＋1／R)／πL11{H｝:H“B(L）”：U“Q”L12W“