精华内容
下载资源
问答
  • 一个大地坐标与高斯坐标转换工具,能够轻松实现经纬度与现实坐标系的转换。
  • 可以将通用的经纬度转为高斯坐标,可选3种两。 可以将通用的经纬度转为高斯坐标,可选3种两。 可以将通用的经纬度转为高斯坐标,可选3种两。 可以将通用的经纬度转为高斯坐标,可选3种两
  • 本程序是以EXCEL为平台编写的工具,主要功能是进行高斯坐标正反处及换带.
  • 一款高斯坐标转换的软件,简单实用。可以直接实现大地坐标向高斯坐标的转换。
  • 高斯投影转换(正算、反算、换带),支持多行计算,鼠标点击可以查看使用说明。功能很实用。
  • 程序包含两个函数,其中[Gaussian_X,Gaussian_Y]=convert84BLToGauss(longitude,latitude) 函数将经纬度转为o-xy坐标;[longitude,latitude]=convert84GaussToBL(X,Y) 函数将o-xy坐标转为经纬度;选用的6度带宽。
  • 高斯坐标转换小工具,非常好用,可以在不同坐标系间转换
  • 大地坐标高斯/UTM投影计算工具

    万次阅读 2018-05-21 13:00:51
    今天要跟大家分享的是一个坐标投影计算小工具,提供北京54、西安80、WGS84和CGCS2000坐标系进行高斯投影和UTM投影正反算转换计算(平面坐标与大地坐标相互转换),提供地理坐标度分秒格式与十进制度度格式的相互转换...

        今天要跟大家分享的软件是一款坐标投影计算小工具,提供北京54、西安80、WGS84和CGCS2000坐标系进行高斯投影和UTM投影正反算计算(投影坐标与地理坐标相互转换),提供地理坐标度分秒格式与十进制度度格式的相互转换。本软件不提供坐标系之间的相互转换!!!


        坐标系和投影是测绘和GIS专业比较基础但也比较容易弄晕乎的知识点,因此在介绍该软件功能之前,先给大家普及下坐标系和投影的相关知识,讲的不对的地方,欢迎拍砖。

    1、坐标系和投影

    1)坐标系  

        测绘里的坐标系包括大地坐标系(全球的)和地方坐标系(局部的),地方坐标系相对复杂(参数保密),本文仅讨论大地坐标系。我国常用的大地坐标系包括北京54、西安80、WGS84和CGCS2000坐标系。

        大地坐标系是大地测量中以参考椭球面为基准面建立起来的坐标系,大地坐标系分为参心坐标系和地心坐标系2种。参心坐标系是以参考椭球的几何中心为原点的大地坐标系,比如我国的北京54坐标系和西安80坐标系,参心坐标系有大地原点(如西安大地原点);地心坐标系是以地球质心为原点建立的坐标系,比如WGS84和CGCS2000坐标系,地心坐标系没有大地原点。

       大地坐标系包括地理坐标系和投影坐标系。地理坐标系为球面坐标,坐标值是经纬度,包括度分秒格式和各种十进制格式;投影坐标系为平面坐标,坐标值是XY坐标。


    2)投影      

        地理坐标按照一定的方式进行投影计算可以转换成投影坐标,我国常用的投影方式包括高斯投影和UTM投影。UTM投影与高斯投影很相似,唯一差别就是UTM投影的比例系数是0.9996,而高斯投影的比例系数是1。下面以高斯投影为例介绍坐标投影。


        假想有一个椭圆柱面横套在地球椭球体外面,并与某一条子午线(此子午线称为中央子午线)相切,椭圆柱的中心轴通过椭球体中心,然后用一定投影方法,将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上,再将此柱面展开即成为投影面(如上图所示),此投影为高斯投影。

         高斯投影坐标系中的XY轴与ArcGIS软件或AutoCAD软件中的平面坐标系(笛卡尔平面坐标系)中的XY轴刚好是相反的(请看上图右侧的坐标系)。

        高斯投影有严密的计算公式,通过高斯投影正算公式可以将地理坐标转换成投影坐标,通过高斯投影反算公式可以将投影坐标转换成地理坐标,下面要介绍的坐标投影计算软件就是根据这两个公式进行坐标转换计算,具体计算公式请问百度。

        高斯投影计算有个重要的参数就是中央子午线也叫中央经线,因此需要引入6度带和3度带的概念,6度带是自0度经线起每隔经差6度自西向东分带(见下图),依次编号1,2,3,…,60,中央经线的经度用L0表示,它们的关系是:L0=6n-3(n为带号),我国范围内6度带有11个(43~53带);3度带从1.5度经线起每隔经差3度自西向东分带,依次编号1,2,3,…,120,中央经线的经度用L0表示,它们的关系是:L0=3n(n为带号),我国范围内3度带有22个(24~45带)。

        UTM投影通常只采用6度分带,通常也只针对WGS84坐标系。


        设置投影分带的目的是为了在全球范围内统一投影分带标准,但对于坐标投影计算而言,6度带和3度带的目的只是为了计算出中央经线。其实中央经线也可以自定义,比如选取离项目范围最近的整数倍经线作为中央经线,离中央经线距离越小,投影变形越小。

        将地理坐标按照高斯投影计算公式转换成的投影坐标一般包括2种表达,一种是含投影带号的(如:Y=38457136.664, X=2455475.927),另一种是不含投影带号的(如:Y=457136.664, X=2455475.927)。

    3)经验小结    

    • 每个大地坐标系都有两类坐标值,分别是地理坐标(经纬度)和投影坐标(XY坐标);

    • 相同大地坐标系的地理坐标与投影坐标可以通过投影正/反算公式进行相互转换,这种相互转换是严密的,不会产生误差;

    • 本软件的主要功能就是相同大地坐标系中地理坐标与投影坐标之间相互转换,当然ArcGIS软件也有该功能,但要先将坐标值转换成图层后再进行转换操作,步骤相对繁琐。在后面的文章中会单独介绍ArcGIS软件中坐标系和投影的相关知识。

    2、软件功能介绍

        下面介绍的坐标投影计算小工具主要包括2个功能:

    • 提供地理坐标与投影坐标的相互转换(坐标正反投影计算);

    • 提供地理坐标度分秒格式与十进制度格式的相互转换。

    1)坐标正反投影计算

        该功能提供北京54、西安80、WGS84和CGCS2000坐标系进行高斯投影和UTM投影正反算计算。左侧坐标列表为投影坐标,右侧坐标列表为地理坐标,中间的下拉框可以选择坐标系和投影方式,通过设置投影带和带号可以计算出中央经线,也支持手动输入中央经线进行转换(原理不很清楚的,慎用)。


        还可以选择保留小数位数和输出数据列的分割符;右侧地理坐标列表中支持输入“经度-纬度”或者“纬度-经度”的坐标格式,只需要在经纬度下拉框中选择对应的“经纬度”或“纬经度”格式;投影坐标分为“去带号”和“保留带号”2种格式,可以根据需要选择。

    2)地理坐标经纬度格式转换

        该功能提供地理坐标的度分秒格式与十进制度格式之间的相互转换,该计算公式比较简单:

    • 度分秒转十进制度:十进制度=度+分/60+秒/3600;

    • 十进制度转度分秒:度=十进制度取整,分=(十进制度取小数*60)再取整,秒=(十进制度取小数*60)取小数*60。


    3、注意事项

        1)该软件对坐标输入格式有要求,坐标值通过逗号或制表符(Excel中数据列分隔方式)隔开,其他分割方式不支持;另外,度分秒的符号支持搜狗输入法中输出的度分秒符号;

        2)该软件使用时结合文本编辑器和Excel软件进行操作会更方便;

        3)该软件虽然简单易用,但需要对坐标系和投影知识比较熟悉才能正确使用;另外,该软件虽然经过了多年的实际工作验证,但本人不保证转换结果的准确,请验证后使用。

    4、下载地址

        免费下载软件请关注微信公众号。


    展开全文
  • 北京54,西安80,wgs84坐标都是是大地坐标,也就是我们通常所说的经纬度坐标,但是它们基于的椭球体不同,我国当前的基本比例尺地形图都是基于北京54和西安80的,而GPS接受的定位数据是基于WGS84的。 UTM投影坐标 ...

    北京54,西安80wgs84坐标都是是大地坐标,也就是我们通常所说的经纬度坐标,但是它们基于的椭球体不同,我国当前的基本比例尺地形图都是基于北京54和西安80的,而GPS接受的定位数据是基于WGS84的。

    UTM投影坐标

    UTM是一种投影坐标,是将球面经纬度坐标经过投影算法转换成的平面坐标,即通常所说的XY坐标.

    WGS-84坐标系

    WGS-84坐标系(World Geodetic System)是一种国际上采用的地心坐标系。坐标原点为地球质心,其地心空间直角坐标系的Z轴指向国际时间局(BIH)1984.0定义的协议地极(CTP)方向,X轴指向BIH1984.0的协议子午面和CTP赤道的交点,Y轴与Z轴、X轴垂直构成右手坐标系,称为1984年世界大地坐标系。这是一个国际协议地球参考 系统 (ITRS),是目前国际上统一采用的大地坐标系。 GPS 广播星历是以WGS-84坐标系为根据的。
    WGS84坐标系,长轴6378137.000m,短轴6356752.314,扁率1/298.257223563。

    转换参数

    由于采用的椭球基准不一样,并且由于投影的局限性,使的全国各地并不存在一至的转换参数。对于这种转换由于量较大,有条件的话,一般都采用GPS联测已知点,应用GPS 软件 自动完成坐标的转换。当然若条件不许可,且有足够的重合点,也可以进行人工解算。

     

    利用Matlab进行地理坐标和直角坐标相互转换

    Matlab版本:(R2017b)

    这里需要用到的命令有:maps, defaultm, projfwd, projinv

    1,  查看有哪些地图投影方式:

    运行maps查看有哪些投影方式

    >> maps

    MapTools Projections

    地图工具投影
    类名ID字符串
    圆柱形巴尔萨沙特圆柱壳
    圆柱形贝尔曼圆柱贝尔曼
    柱状大黄阿特拉斯MiRA*BSAM
    柱面布劳恩透视圆柱*布劳恩
    圆柱形卡西尼圆柱标准筒
    圆柱形卡西尼圆柱卡西尼
    圆柱形中心圆柱
    圆柱形等面积圆柱壳
    圆柱等距圆柱形EQDCILIN
    圆柱形五边形等值线
    圆柱形正投影GOORO
    圆柱形全景立体图
    圆柱形朗伯圆柱拉姆克林
    圆柱墨卡托圆柱墨卡托
    圆柱形铣床圆柱铣床
    圆柱形板
    圆柱横墨卡托
    Cylindrical Trystan Edwards柱形三角形
    圆柱形通用横向墨卡托(UTM)UTM
    圆柱形润湿圆柱
    假瓢虫
    伪圆柱形夹心碰撞
    假圆柱形岩纹抛物线岩
    Pseudocylindrical Eckert I*ECKTT1
    Pseudocylindrical Eckert II EcKTT2
    伪圆柱埃克特Ⅲ*ECKTT3
    伪圆柱埃克特ⅣECKTT4
    伪圆柱埃克特V*ECKTT5
    伪圆柱埃克特ⅥECKTT6
    伪圆柱形平极抛物线扁钢
    伪圆心平极四次平坦
    伪圆柱平极极正弦平片
    伪圆柱形富尼尔四倍体
    假柱状古鳕
    假圆柱哈塔诺对称等面积哈塔诺
    伪圆柱形KavaRISKY V KAVRSky5
    伪圆柱形KavaRISKY VI KAVRSky6
    Pseudocylindrical Loximuthal *洛西姆
    Pseudocylindrical Modified Sinusoidal(蒂索)*莫迪辛
    Pseudocylindrical Mollweide莫尔维德
    假柱状紫杉醇P5*PtTNS5
    伪圆柱四次自反四次方
    伪圆柱鲁滨孙*鲁滨孙
    伪圆柱正弦正弦曲线
    伪圆柱瓦格纳ⅣWAGELN4
    Pseudocylindrical Winkel·I·温克尔
    圆锥等面积圆锥(阿尔伯斯)-标准EQACONICSTD
    圆锥等面积圆锥(阿尔伯斯)方程
    圆锥等距圆锥标准EQDCONICSTD
    圆锥等距圆锥曲线
    圆锥-朗伯共形圆锥标准LabBurtSTD
    圆锥型朗伯共形圆锥朗伯
    圆锥默多克I圆锥
    圆锥默多克Ⅲ最小误差二次曲线MurdoCH3
    多圆锥多圆锥标准多面体
    多圆锥多圆锥聚变
    多圆锥Van Dr.GruttNⅠVGRIT1
    假性圆锥骨
    伪圆锥沃纳沃纳
    方位共振谐波平均值
    方位等面积方位角
    方位等距方位角
    方位地球仪
    方位角
    方位正投影仪

    CLASS                 NAME                                ID STRING        

    Cylindrical           Balthasart Cylindrical              balthsrt        

    Cylindrical           Behrmann Cylindrical                behrmann        

    Cylindrical           Bolshoi Sovietskii Atlas Mira*      bsam            

    Cylindrical           Braun Perspective Cylindrical*      braun           

    Cylindrical           Cassini Cylindrical-Standard        cassinistd      

    Cylindrical           Cassini Cylindrical                 cassini         

    Cylindrical           Central Cylindrical*                ccylin          

    Cylindrical           Equal Area Cylindrical              eqacylin        

    Cylindrical           Equidistant Cylindrical             eqdcylin        

    Cylindrical           Gall Isographic                     giso            

    Cylindrical           Gall Orthographic                   gortho          

    Cylindrical           Gall Stereographic*                 gstereo         

    Cylindrical           Lambert Cylindrical                 lambcyln        

    Cylindrical           Mercator Cylindrical                mercator        

    Cylindrical           Miller Cylindrical*                 miller          

    Cylindrical           Plate Carree                        pcarree         

    Cylindrical           Transverse Mercator                 tranmerc        

    Cylindrical           Trystan Edwards Cylindrical         trystan         

    Cylindrical           Universal Transverse Mercator (UTM)    utm             

    Cylindrical           Wetch Cylindrical*                  wetch           

    Pseudocylindrical     Apianus II*                         apianus         

    Pseudocylindrical     Collignon                           collig          

    Pseudocylindrical     Craster Parabolic                   craster         

    Pseudocylindrical     Eckert I*                           eckert1         

    Pseudocylindrical     Eckert II                           eckert2         

    Pseudocylindrical     Eckert III*                         eckert3         

    Pseudocylindrical     Eckert IV                           eckert4         

    Pseudocylindrical     Eckert V*                           eckert5         

    Pseudocylindrical     Eckert VI                           eckert6         

    Pseudocylindrical     Flat-Polar Parabolic                flatplrp        

    Pseudocylindrical     Flat-Polar Quartic                  flatplrq        

    Pseudocylindrical     Flat-Polar Sinusoidal               flatplrs        

    Pseudocylindrical     Fournier                            fournier         

    Pseudocylindrical     Goode Homolosine                    goode           

    Pseudocylindrical     Hatano Assymmetrical Equal Area     hatano          

    Pseudocylindrical     Kavraisky V                         kavrsky5        

    Pseudocylindrical     Kavraisky VI                        kavrsky6        

    Pseudocylindrical     Loximuthal*                         loximuth        

    Pseudocylindrical     Modified Sinusoidal (Tissot)*       modsine         

    Pseudocylindrical     Mollweide                           mollweid        

    Pseudocylindrical     Putnins P5*                         putnins5        

    Pseudocylindrical     Quartic Authalic                    quartic         

    Pseudocylindrical     Robinson*                           robinson        

    Pseudocylindrical     Sinusoidal                          sinusoid        

    Pseudocylindrical     Wagner IV                           wagner4         

    Pseudocylindrical     Winkel I*                           winkel          

    Conic                 Equal Area Conic (Albers)-Standard    eqaconicstd     

    Conic                 Equal Area Conic (Albers)           eqaconic        

    Conic                 Equidistant Conic-Standard          eqdconicstd     

    Conic                 Equidistant Conic                   eqdconic        

    Conic                 Lambert Conformal Conic-Standard    lambertstd      

    Conic                 Lambert Conformal Conic             lambert         

    Conic                 Murdoch I Conic*                    murdoch1        

    Conic                 Murdoch III Minimum Error Conic*    murdoch3        

    PolyConic             Polyconic-Standard                  polyconstd      

    PolyConic             Polyconic                           polycon         

    PolyConic             Van Der Grinten I*                  vgrint1         

    PseudoConic           Bonne                               bonne           

    PseudoConic           Werner                              werner          

    Azimuthal             Breusing Harmonic Mean*             breusing        

    Azimuthal             Equal Area Azimuthal (Lambert)      eqaazim         

    Azimuthal             Equidistant Azimuthal*              eqdazim         

    Azimuthal             Globe                               globe            

    Azimuthal             Gnomonic*                           gnomonic        

    Azimuthal             Orthographic*                       ortho           

    Azimuthal             Stereographic                       stereo          

    Azimuthal             Universal Polar Stereographic       ups             

    Azimuthal             Vertical Perspective*               vperspec        

    Pseudoazimuthal       Wiechel Equal Area*                 wiechel         

    Modified Azimuthal    Aitoff*                             aitoff          

    Modified Azimuthal    Briesemeister*                      bries           

    Modified Azimuthal    Hammer*                             hammer         

    * Denotes availability for sphere only

    2,  定义投影结构

    这里对mercator投影结构进行定义,然后对mstruct进行补全

    >>mstruct=defaultm('mercator');

    %定义椭球体长轴,椭率,坐标原点

    >>mstruct.geoid=[ 6378137          ,0.0818191908426215];

    >>mstruct.origin=[0,0,0];

    >>mstruct=defaultm(mstruct);

    3,  进行投影和数据转换

    地理坐标投影到直角坐标:

    >>[x,y] =projfwd(mstruct,lat,lon);

    projfwd - Forward map projection using PROJ.4 map projection library

        This MATLAB function returns the x and y map coordinates from the forward
        projection transformation.

        [x,y] = projfwd(proj,lat,lon)

     

    直角坐标投影到地理坐标:

    >>[lat,lon]=projinv(mstruct,x,y);

    projinv - Inverse map projection using PROJ.4 map projection library

        This MATLAB function returns the latitude and longitude values from the inverse
        projection transformation.

        [lat,lon] = projinv(proj,x,y)

     

    projlist - Map projections supported by projfwd and projinv

        This MATLAB function displays a table of projection names, IDs, and
        availability.

        projlist(listmode)
        S = projlist(listmode)

     

    GPS的输出一般为经纬度坐标,如:北纬(N)41度30分5秒, 东经(E)114度10分30秒。
    UTM(Universal Transverse Mercator)坐标是一个投影坐标系,将地球分为60个区,单位为m。
    二者的转换可以利用商业软件,如Global Mapper等。下面给出了经纬度转UTM的matlab程序。

        转换公式:(只针对北半球,中国地区计算方法)

        Matlab 程序:
       

        运行结果:单位为km
       

    Global Mapper结果对比:

    结果与商业软件在毫米位上保持一致,验证程序的有效性。

    展开全文
  • 高斯坐标经纬度转换器是一款经典的坐标转换工具,支持不同的度投影功能,可以通过输入大地的坐标进行正算反算,得出高斯投影的坐标,非常精准地就能转换出经纬度。软件介绍高斯坐标转换工具最新版是一款功能强劲的...

    高斯坐标经纬度转换器是一款经典的坐标转换工具,支持不同的度带投影功能,可以通过输入大地的坐标进行正算反算,得出高斯投影的坐标,非常精准地就能转换出经纬度。

    软件介绍

    高斯坐标转换工具最新版是一款功能强劲的坐标转换工具,专业实用可以快速的进行各种球面的坐标、平面的坐标之间进行相互的转换,而且也对文件的格式进行支持,可以轻松的即可实现各种投影之间的相互转换,也支持对单独的投影坐标进行转换。

    如果手动算,可以这样做。不过现在一般都用坐标转换工具,可以直接在软件里面设置前、后坐标系统,然后输入转换参数就可以了。

    软件功能

    可以对三度带投影进行快速的计算

    也支持对六度带投影进行计算

    也支持对三度带领带换算

    也支持对六度带领带换算

    也可对三度带,六度带进行计算

    也支持对六度带、三度带进行换算

    高斯坐标转换工具软件特色

    高斯坐标即高斯-克吕格坐标系

    高斯-克吕格(Gauss-Kruger)投影简称高斯投影

    又名等角横切椭圆柱投影

    地球椭球面和平面间正形投影的一种

    根据需要选择对应的换算方式

    坐标正算是将经纬坐标转化为高斯平面直角坐标

    坐标反算是将高斯平面直角坐标转化为经纬坐

    选择相应的椭球参数和投影带

    使用说明

    打开高斯坐标转换器,坐标正算是将经纬坐标转化为高斯平面直角坐标,坐标反算是将高斯平面直角坐标转化为经纬坐标,然后选择相应的椭球参数和投影带。

    展开全文
  • 地理坐标高斯投影坐标互相转换

    千次阅读 2019-12-02 22:17:46
    前言:原作者提供了js版本的坐标转换函数,和三种椭球体参数(北京54、西安80、WGS84),个人在百度上搜索了CGC2000坐标系的椭球参数,添加到函数里,使用java写的转换函数。在ArcMap里使用同一位置的两个不同坐标系...

    前言:原作者提供了js版本的坐标转换函数,和三种椭球体参数(北京54、西安80、WGS84),个人在百度上搜索了CGC2000坐标系的椭球参数,添加到函数里,使用java写的转换函数。在ArcMap里使用同一位置的两个不同坐标系的点进行验证,误差还算可以接受。。。

    转换函数代码:

    public class Convertor {
        private double a;//'椭球体长半轴
        private double b;// '椭球体短半轴
        private double f; //'扁率
        private double e;// '第一偏心率
        private double e1; //'第二偏心率
        private double FE;//'东偏移
        private double FN;//'北偏移
        private double L0;//'中央经度
        private double W0;//'原点纬线
        private double k0;//'比例因子
        /**
         * 幂函数
         * @param e
         * @param i
         * @return
         */
        private double MZ(double e, int i) {
            return Math.pow(e,i);
        }
     
        /**
         * 说明: 用于初始化转换参数
         *
         * @param TuoqiuCanshu    枚举类型,提供北京54、西安80和WGS84三个椭球参数
         * @param CentralMeridian 中央经线
         * @param OriginLatitude  原点纬度,如果是标准的分幅,则该参数是0
         * @param EastOffset      东偏移
         * @param NorthOffset     北偏移
         */
        public Convertor(int TuoqiuCanshu, double CentralMeridian, double OriginLatitude, double EastOffset, double NorthOffset) {
            /**
             * 'Krassovsky (北京54采用) 6378245 6356863.0188
             * 'IAG 75(西安80采用) 6378140 6356755.2882
             * 'WGS 84 6378137 6356752.3142
             * 'CGC 2000 6378137 6356752.31414
             */
            if (TuoqiuCanshu == 0)//北京五四
            {
                a = 6378245;
                b = 6356863.0188;
            } else if (TuoqiuCanshu == 1)// '西安八零
            {
                a = 6378140;
                b = 6356755.2882;
            }
            if (TuoqiuCanshu == 2)//'WGS84
            {
                a = 6378137;
                b = 6356752.3142;
            }
            if (TuoqiuCanshu == 3)//'CGC2000坐标
            {
                a = 6378137;
                b = 6356752.31414;
            }
            f = (a - b) / a;//扁率
            //e = Math.sqrt(1 - MZ((b / a) ,2));//'第一偏心率
            e = Math.sqrt(2 * f - MZ(f, 2));//'第一偏心率
            //eq = Math.sqrt(MZ((a / b) , 2) - 1);//'第二偏心率
            e1 = e / Math.sqrt(1 - MZ(e, 2));//'第二偏心率
            L0 = CentralMeridian;//中央经
            W0 = OriginLatitude;//原点纬线
            k0 = 1;//'比例因子
            FE = EastOffset;//东偏移
            FN = NorthOffset;//北偏移
        }
     
        /**
         * 经纬度坐标转高斯投影坐标
         * @param J 经度
         * @param W 纬度
         * @return
         */
        public double[] JWgetGK(double J, double W) {
            //'给出经纬度坐标,转换为高克投影坐标
            double[] resultP = new double[2];
            double BR = (W - W0) * Math.PI / 180;//纬度弧长
            double lo = (J - L0) * Math.PI / 180; //经差弧度
            double N = a / Math.sqrt(1 - MZ((e * Math.sin(BR)), 2)); //卯酉圈曲率半径
            //求解参数s
            double B0;
            double B2;
            double B4;
            double B6;
            double B8;
            double C = MZ(a, 2) / b;
            B0 = 1 - 3 * MZ(e1, 2) / 4 + 45 * MZ(e1, 4) / 64 - 175 * MZ(e1, 6) / 256 + 11025 * MZ(e1, 8) / 16384;
            B2 = B0 - 1;
            B4 = 15 / 32 * MZ(e1, 4) - 175 / 384 * MZ(e1, 6) + 3675 / 8192 * MZ(e1, 8);
            B6 = 0 - 35 / 96 * MZ(e1, 6) + 735 / 2048 * MZ(e1, 8);
            B8 = 315 / 1024 * MZ(e1, 8);
            double s = C * (B0 * BR + Math.sin(BR) * (B2 * Math.cos(BR) + B4 * MZ((Math.cos(BR)), 3) + B6 * MZ((Math.cos(BR)), 5) + B8 * MZ((Math.cos(BR)), 7)));
            double t = Math.tan(BR);
            double g = e1 * Math.cos(BR);
            double XR = s + MZ(lo, 2) / 2 * N * Math.sin(BR) * Math.cos(BR) + MZ(lo, 4) * N * Math.sin(BR) * MZ((Math.cos(BR)), 3) / 24 * (5 - MZ(t, 2) + 9 * MZ(g, 2) + 4 * MZ(g, 4)) + MZ(lo, 6) * N * Math.sin(BR) * MZ((Math.cos(BR)), 5) * (61 - 58 * MZ(t, 2) + MZ(t, 4)) / 720;
            double YR = lo * N * Math.cos(BR) + MZ(lo, 3) * N / 6 * MZ((Math.cos(BR)), 3) * (1 - MZ(t, 2) + MZ(g, 2)) + MZ(lo, 5) * N / 120 * MZ((Math.cos(BR)), 5) * (5 - 18 * MZ(t, 2) + MZ(t, 4) + 14 * MZ(g, 2) - 58 * MZ(g, 2) * MZ(t, 2));
            resultP[0] = YR + FE;
            resultP[1] = XR + FN;
            return resultP;
        }
     
        /**
         * 高斯投影坐标 转为 经纬度坐标
         * @param X 高斯投影坐标X
         * @param Y 高斯投影坐标Y
         * @return
         */
        public double[] GKgetJW(double X, double Y) {
            //'给出高克投影坐标,转换为经纬度坐标
            double[] resultP =new double[2];
            double El1 = (1 - Math.sqrt(1 - MZ(e, 2))) / (1 + Math.sqrt(1 - MZ(e, 2)));
            double Mf = (Y - FN) / k0;//真实坐标值
            double Q = Mf / (a * (1 - MZ(e, 2) / 4 - 3 * MZ(e, 4) / 64 - 5 * MZ(e, 6) / 256));//角度
            double Bf = Q + (3 * El1 / 2 - 27 * MZ(El1, 3) / 32) * Math.sin(2 * Q) + (21 * MZ(El1, 2) / 16 - 55 * MZ(El1, 4) / 32) * Math.sin(4 * Q) + (151 * MZ(El1, 3) / 96) * Math.sin(6 * Q) + 1097 / 512 * MZ(El1, 4) * Math.sin(8 * Q);
            double Rf = a * (1 - MZ(e, 2)) / Math.sqrt(MZ((1 - MZ((e * Math.sin(Bf)), 2)), 3));
            double Nf = a / Math.sqrt(1 - MZ((e * Math.sin(Bf)), 2));//卯酉圈曲率半径
            double Tf = MZ((Math.tan(Bf)), 2);
            double D = (X - FE) / (k0 * Nf);
            double Cf = MZ(e1, 2) * MZ((Math.cos(Bf)), 2);
            double B = Bf - Nf * Math.tan(Bf) / Rf * (MZ(D, 2) / 2 - (5 + 3 * Tf + 10 * Cf - 9 * Tf * Cf - 4 * MZ(Cf, 2) - 9 * MZ(e1, 2)) * MZ(D, 4) / 24 + (61 + 90 * Tf + 45 * MZ(Tf, 2) - 256 * MZ(e1, 2) - 3 * MZ(Cf, 2)) * MZ(D, 6) / 720);
            double L = L0 * Math.PI / 180 + 1 / Math.cos(Bf) * (D - (1 + 2 * Tf + Cf) * MZ(D, 3) / 6 + (5 - 2 * Cf + 28 * Tf - 3 * MZ(Cf, 2) + 8 * MZ(e1, 2) + 24 * MZ(Tf, 2)) * MZ(D, 5) / 120);
            double Bangle = B * 180 / Math.PI;
            double Langle = L * 180 / Math.PI;
            resultP[0] = Langle;//RW * 180 / Math.PI;
            resultP[1] = Bangle + W0;//RJ * 180 / Math.PI;
            return resultP;
        }
    }

     

    展开全文
  • 因为我们是获取手机的gps坐标,我们起先是使用的甲方提供的坐标转换服务,发现84大地转到深圳独立坐标(深圳高斯平面坐标)时精度误差极大。于是甲方给我们提供了部分控制点信息,让我们自己去完成坐标的转换,不再...
  • 本文整理了一些地理坐标系之间的转换(Java代码)pom依赖com.vividsolutionsjts1.13org.osgeoproj4j0.1.0坐标转换工具类package com.navinfo.datalake.algorithm.geometry;import com.mytest.algorithm.model.Pixel;...
  • 功能:本工具可将百度坐标系的XY(经纬度)坐标转换到WGS84下的XY(经纬度)坐标并生成SHP图层。 注:内附使用说明。
  • C++版本的经纬度坐标转化工具,包含Debug和Release版本的开发库,可以使用Qt和VS2019的集成开发环境进行集成。
  • 【摘 要】 本文针对从事测绘工作者普遍遇到的坐标转换问题,简要介绍ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标转换原理和步骤。【关键词】 ArcGIS 坐标转换 投影变换1 坐标转换简介坐标系统之间的坐标转换...
  • 选择投影坐标系,此处根据坐标大致推算出处于3度的第40度号中心的经度是120度 数据加载完成后,再导出成shp数据,如上图。 此时属性表数据只有X Y Z数据 3 坐标转换 坐标变换,打开ArcToolbox,找到“数据...
  • 坐标批量转换工具

    2011-10-14 13:04:15
    实现直角坐标与各地理坐标坐标的批量转换
  • 高斯坐标正反转换

    2019-06-03 10:10:54
    在工作中遇到坐标转换的问题,给出的坐标是从网上找到的excel表格形式,单位是km,经多次尝试,搞懂了转换方法。 问题分为两种 一,地理坐标转平面坐标(度分秒转米) 这个问题可以利用arcgis来解决,亦可以通过...
  • 地理坐标转换工具合集包括高斯坐标正反算软件,七参数坐标转换工具,坐标换带工具。能转换度分秒转经纬度等等,经纬度转度分秒等,是用于地质坐标转换非常好的一个工具!使用方法高斯坐标正算是将经纬坐标转化为高斯...
  • 提供高斯投影坐标、平面坐标等不同坐标系变换
  • 高斯投影与坐标转换

    千次阅读 2018-11-30 20:26:00
    高斯-克吕格投影特点: 1、中央子午线是直线,其长度不变形;其他子午线是凹向中央子午线的弧线,并以中央子午线为对称轴; 2、赤道线是直线,但有长度变形;其他纬线为凸向赤道的弧线,并以赤道为对称轴; 3、...
  • 1.由于ARCGIS中没有高斯3度投影,所以需要我们自己定义高斯3度投影坐标系(此处以北京为例),并保存 PS:三度投影,中央经线选取时需要根据当地经度就近选择中央经线。高斯3度中央...
  • 投影换带 大地坐标系、投影坐标系、空间直角坐标系的任意互转 程序特点: 支持任意椭球参数 程序内定义了北京54、西安80、WGS84、CGCS2000等常用椭球基准,并且可以采用自定义椭球,只需要输入椭球长轴半径...
  • 【问题描述】:很多情况下,我们知道某一地点的坐标(经纬度),需要计算其在高斯克吕格投影中的号及中央经度。关于该问题,有具体的公式可言,只是计算过程稍微繁琐一些,当然啦,我们可以写程序来解决,谁叫我们...
  • 这在测绘学上就称作高斯投影正反算和高斯坐标换带。 坐标正反算和坐标换带的方法,有很多软件、工具都可以进行计算,包括目前主流的GIS软件、测绘软件,都不在话下。 但王石头认为: 某海外GIS软件庞大臃肿,...
  • 1、常见的地心坐标格式 地心坐标:Y=35607216,X=3024388;...已知条件:一个高斯平面坐标(X,Y)、对应的地理坐标系(西安80、北京54、WGS84等)以及经纬度(L,B)。 计算过程:利用经度L计算出中央...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 9,606
精华内容 3,842
关键字:

高斯坐标换带工具

友情链接: qanlei.zip