求帮忙高斯坐标转换经纬度

Posted 2023-05-16

x y
4000151.6933 38407563.1230
4000151.6933 38407263.1230
4000151.6933 38406963.1230
4000151.6933 38406663.1230
4000151.6933 38406383.1230
4000151.6933 38406063.1230
以上几个坐标怎么在google地图上找到啊，我找了高斯坐标转换软件，但是算出来以后输入到googel地图里面是在湖北，而实际应该在山西长治，所以很困惑啊。有没有高手帮忙啊！

B，L，
036:07:38.070270N,112:58:23.352558E,
036:07:37.967265N,112:58:11.356987E,
036:07:37.863926N,112:57:59.361433E,
036:07:37.760254N,112:57:47.365896E,
036:07:37.663192N,112:57:36.170077E,
036:07:37.551909N,112:57:23.374874E,
格式如：36度7分37.551909秒。

经过验证 位置在长治市长治县杜村乡（大致）。 参考技术A 　　高斯坐标转换经纬度：
　　1、首先百度下载高斯坐标转换器。
　　2、需要转换的经纬度输入进去即可。
　　3、点击转换，完成转换。
　　高斯坐标转换器是一款专业坐标转换软件，本软件拥有多种功能，支持高斯坐标正算和高斯坐标反算，可满足学生、办公人员、科研工作者等日常学习工作需求。打开高斯坐标转换器，根据需要选择对应的换算方式，坐标正算是将经纬坐标转化为高斯平面直角坐标，坐标反算是将高斯平面直角坐标转化为经纬坐标，然后选择相应的椭球参数和投影带。高斯坐标即高斯-克吕格坐标系。高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”，又名"等角横切椭圆柱投影”，地球椭球面和平面间正形投影的一种。

高斯克吕格投影，将经纬度转换为投影坐标

1. 预先定义椭球参数

 1         /************************************************************************/
 2         /*                                 84椭球                               */
 3         double a = 6378137.0;
 4         double b = 6356752.3142;
 5         double f = 1 / 298.257223563;
 6         double c = a / (1 - f);
 7         double e2 = 2 * f - f * f;
 8         double ep2 = 1 / ((1 - f) * (1 - f)) - 1;
 9         double dpi = 0.0174532925199432957692;
10         double m0 = a * (1 - e2);
11         double m2 = 1.5 * e2 * m0;
12         double m4 = 1.25 * e2 * m2;
13         double m6 = 7 * e2 * m4 / 6;
14         double m8 = 9 * e2 * m6 / 8;
15         double a0 = m0 + m2 / 2 + 3 * m4 / 8 + 5 * m6 / 16 + 35 * m8 / 128;
16         double a2 = m2 / 2 + m4 / 2 + 15 * m6 / 32 + 7 * m8 / 16;
17         double a4 = m4 / 8 + 3 * m6 / 16 + 7 * m8 / 32;
18         double a6 = m6 / 32 + m8 / 16;
19         double ml = 114;    // 中央经线 根据需要进行设置
20 
21         double _iproj[12] = { a, b, f, c, e2, ep2, dpi, a0, a2, a4, a6, ml };

2. 编写转换函数

 1     int latlon2xy(double &x, double &y)
 2     {
 3         double lat = x * _iproj[6];
 4         double dL = (y - _iproj[11]) * _iproj[6];
 5 
 6         double X = _iproj[7] * lat - _iproj[8] * sin(2 * lat) / 2 + _iproj[9] * sin(4 * lat) / 4 - _iproj[10] * sin(6 * lat) / 6;
 7         double tn = tan(lat);
 8         double tn2 = tn * tn;
 9         double tn4 = tn2 * tn2;
10 
11         double j2 = (1 / pow(1 - _iproj[2], 2) - 1) * pow(cos(lat), 2);
12         double n = _iproj[0] / sqrt(1.0 - _iproj[4] * sin(lat) * sin(lat));
13 
14         double temp[6] = { 0 };
15         temp[0] = n * sin(lat) * cos(lat) * dL * dL / 2;
16         temp[1] = n * sin(lat) * pow(cos(lat), 3) * (5 - tn2 + 9 * j2 + 4 * j2 * j2) * pow(dL, 4) / 24;
17         temp[2] = n * sin(lat) * pow(cos(lat), 5) * (61 - 58 * tn2 + tn4) * pow(dL, 6) / 720;
18         temp[3] = n * cos(lat) * dL;
19         temp[4] = n * pow(cos(lat), 3) * (1 - tn2 + j2) * pow(dL, 3) / 6;
20         temp[5] = n * pow(cos(lat), 5) * (5 - 18 * tn2 + tn4 + 14 * j2 - 58 * tn2 * j2) * pow(dL, 5) / 120;
21 
22         y = X + temp[0] + temp[1] + temp[2];
23         x = temp[3] + temp[4] + temp[5];
24     }