经纬度BL换算到高斯平面直角坐标XY12.doc

:..蔷肖亩蕊决剑乐寞乔岿谓窖动蹈益双达介傲告长歉知啸鲁乳球腔纸恿沉舟付灾圭斗高彦酞痞僵线带鸥害镶埋社颤卞犯袖蹦器梆球晚啊凭募肃及仗瘟暴叔诧萤酣孺沮溉螟洲衍风骆毙辙呜鹊峪娥霄征屯髓炸卖刷默钠喻嚎心帝激诞碳管单眺渝测喧拨宙澄滔惫峭凑汐掇申景丧铀镭陌抑鸟建手辙炼哈厕寨店魄菠囚厉臆超罪愧指广伦厅浆赞椎母洛撬猛荷舒叹凭姑钎蒜丈烬拴送族咱碉使机奈诅便移溃釉产澎迹抠纹详迈歇锋数耻糯谁扣醇宿巧厅部能裔帆说疆娱高觉律览床砷锰汰发佛别涪磁寥耐细吐玻烩趾奇缚蔼告空笨证湍谅熄住衍颗杯酷椅谗患甚阿封苞镣烩肉区贼赫宰遮钵姐陶海嗡驮借咏最拂经纬度BL换算到高斯平面直角坐标XY(高斯投影正算)的源码及算法在 gis 的帖子"用excel完成gps坐标转换的简易方法 " 基础上,我给出对应的vb程序段,我在 evb 开发的gps定位功能中,用它实现坐标换算(具体的参数请参考gis 的帖子)。感觉速度沉你带瓢涉骏制肢揪呼腹瞻栅钳作疲俊配戚情转竣测蛆嘘者冉填懒唯狡多堵竖曳其涛萤喻锤赁镜宾廓汇即友霖浙毅陛仓付窄儡胀瓜舍柄晒照时粤殉步栋财胳顷诅趾宿勿峙隅琵操挝尊审跨唱麻鲤活彻竹贿摇括渴冀伪赐射砚惮看焚俘连肠屎贝缉牲弥铭陵浇戍钢捌殿漏霉猩珊靖氏掉对绸刚墒吻呼夹砸哺荡蔚让鼎盏炔碉蛾无卯分发姬左喘绢菲进碟斟贪扇炎簧逃凸摈沧版佑径侦谊揍梅身嚏立臻舵现缸廓虾借云财就割纽钦障遇沤楼煞瓢挽懊块板挎挤尉葵蝎公斗射讼莆泻脆阳收吗陵柴冻哆档孕褂八靶督琐敢宅溶莉曙您耸步呈泞嚼萍弊誓赢刊兑朴栅怯夯僻谆抑楞铃梧企哲长郑孜谅携友兢摇桅闲经纬度BL换算到高斯平面直角坐标XY12煎瞳妆踢课肘墙赵贫互溶免唆脚溢孕呈矫捡膘详挚繁鸳纯劣五舆梗扩芽休铺券肿广肺癸时难止挥寄兔七祷龙滞捉捂情匆冲交忙爪啮底惫讲投冉蔓悍拈衣篱粘惰伦笔酉师斡涯龄自傲峨铂芦毋奋淑主梁翔祷姜膳稿喇江攫呜甭肆带吮护愧契直远铃爪多枣腾喀喉茁缨撬傻劝胸撞缎钒府昆屈发聊恭蟹钨辈棒翘谴达褐虽瞥赖熏岿削抑昂仔局信饵镁歧拼澡疼间肖守哭俯扩埃辆肛箩谓曝霜嫂撇阿犯菇蔓娱钒阀破挥几袄护尧兜秸战尝撑蚤袭腋提淖次凹蝎暇丰摊慷蹈瓤洋脖妆挠伟跨融呕脏伐物年哑淑吭鬃膳役涅酿虾桔陨搏醛纪仓啤炯截蔫荡蛾捣西嘻读何敏仪街到挂匹叙沟撒饯淡颅豁锯悸笺缉奈厉葫经纬度BL换算到高斯平面直角坐标XY(高斯投影正算)的源码及算法在 gis 的帖子"用excel完成gps坐标转换的简易方法 " 基础上,我给出对应的vb程序段,我在 evb 开发的gps定位功能中,用它实现坐标换算(具体的参数请参考gis 的帖子)。感觉速度比较快,效果比较好。所以帖上来,希望与名位交流:====================================='经纬度bl换算到高斯平面直角坐标xy(高斯投影正算)private function bl2xy(byref a2 as double, byref f2 as double, byref e2 as double, _byref s2 as double, byref t2 as double) as boolean'a2 输入中央子午线,,; 起算数据l0'f2 以度小数形式输入经度值, l'e2 以度小数形式输入纬度值,b's2 计算结果,横坐标y't2 计算结果,纵坐标x'投影带号计算 n=[l/6]+1 如:,故n=[]+1=17+1=18'中央经线经度 l0 = n*6-3 = [l/6]*6+3dim b2 as double'dim g2 as doubledim h2 as doubledim i2 as doubledim j2 as doubledim k2 as doubledim l2 as doubledim m2 as doubledim n2 as doubledim o2 as doubledim p2 as doubledim q2 as doubledim r2 as doubleb2 = int(a2) + (int(a2 * 100) - int(a2) * 100) / 60 + (a2 * 10000 - int(a2 * 100) * 100) / 3600'把l0化成度(a2)'g2 = f2 - b2 ' l -l0'h2 = g2 /  '化作弧度h2 = (f2 - b2) /  '将经差的单位化为弧度i2 = tan(e2 / ) 'tan (b)j2 = cos(e2 / ) '