朗伯72转换为Lat/Lon对
本文关键字:Lon Lat 转换 朗伯 | 更新日期: 2023-09-27 18:09:56
我试图合并两个数据库,以巩固两个客户的网站。然而,客户A一直使用常规的纬度/纬度对进行地理定位,而客户B使用朗伯特72 (X/Y)坐标。
我已经建立了一个应该转换这些坐标的脚本(因为我不确定哪些坐标将在最终合并的数据库中使用,我正在尝试转换它们)。我从这里截取了一些片段:http://zoologie.umh.ac.be/tc/algorithms.aspx
请注意,下面提到的所有坐标都指向比利时的位置。
我正在转换一些坐标,看看计算是否正确,但我得到的坐标似乎是正确的。作为参考,比利时的中心大约是(北50.84323737103243,东4.355735778808594),所以我希望所有的坐标都接近这些值。
我将Lambert 72值(X: 151488250, Y: 170492909)转换为纬度/纬度对,但结果是:(-87.538...., -50.724....),这与期望值相差甚远。如果我转换整个循环(Lambert->LatLon->Lambert,反之亦然),我得到的结果值与我输入的相同,因此我知道我的转换至少是一致的,并且转换是彼此的完美反转。
我也尝试了一些在线转换器工具,他们给了我相同的(-87.538....), -50.724....)结果。
由于多个源产生相同的结果,并且我的转换是彼此的正确反转,我认为计算本身是正确的,但结果值仍然需要进一步转换/偏移?
我认为自己代数学得很好,但地图投影学得一塌糊涂。
有人能解释一下吗?
额外信息- 我希望我把这个贴在了正确的论坛上。我真的不确定把它放在哪里,因为这是地理,数学和编码/转换的混合…
- 上面提到的兰伯特坐标(X: 151488250, Y: 170492909)指向布鲁塞尔的一个位置,所以经纬度的结果应该非常接近(North 50.84323737103243, East 4.355735778808594)。
下面是我的转换函数:
public static Lambert72 LatLon_To_Lambert72(LatLon latlon) { var lat = latlon.Lat; var lng = latlon.Lon; double LongRef = 0.076042943; //=4°21'24"983 double bLamb = 6378388 * (1 - (1 / 297)); double aCarre = Math.Pow(6378388, 2); double eCarre = (aCarre - Math.Pow(bLamb, 2)) / aCarre; double KLamb = 11565915.812935; double nLamb = 0.7716421928; double eLamb = Math.Sqrt(eCarre); double eSur2 = eLamb / 2; //conversion to radians lat = (Math.PI / 180) * lat; lng = (Math.PI / 180) * lng; double eSinLatitude = eLamb * Math.Sin(lat); double TanZDemi = (Math.Tan((Math.PI / 4) - (lat / 2))) * (Math.Pow(((1 + (eSinLatitude)) / (1 - (eSinLatitude))), (eSur2))); double RLamb = KLamb * (Math.Pow((TanZDemi), nLamb)); double Teta = nLamb * (lng - LongRef); double x = 0; double y = 0; x = 150000 + 0.01256 + RLamb * Math.Sin(Teta - 0.000142043); y = 5400000 + 88.4378 - RLamb * Math.Cos(Teta - 0.000142043); return new Lambert72(x, y); } public static LatLon Lambert72_To_LatLon(Lambert72 lb72) { double X = lb72.X; double Y = lb72.Y; double LongRef = 0.076042943; //=4°21'24"983 double nLamb = 0.7716421928; double aCarre = Math.Pow(6378388, 2); double bLamb = 6378388 * (1 - (1 / 297)); double eCarre = (aCarre - Math.Pow(bLamb, 2)) / aCarre; double KLamb = 11565915.812935; double eLamb = Math.Sqrt(eCarre); double eSur2 = eLamb / 2; double Tan1 = (X - 150000.01256) / (5400088.4378 - Y); double Lambda = LongRef + (1 / nLamb) * (0.000142043 + Math.Atan(Tan1)); double RLamb = Math.Sqrt(Math.Pow((X - 150000.01256), 2) + Math.Pow((5400088.4378 - Y), 2)); double TanZDemi = Math.Pow((RLamb / KLamb), (1 / nLamb)); double Lati1 = 2 * Math.Atan(TanZDemi); double eSin = 0; double Mult1 = 0; double Mult2 = 0; double Mult = 0; double LatiN = 0; double Diff = 0; double lat = 0; double lng = 0; do { eSin = eLamb * Math.Sin(Lati1); Mult1 = 1 - eSin; Mult2 = 1 + eSin; Mult = Math.Pow((Mult1 / Mult2), (eLamb / 2)); LatiN = (Math.PI / 2) - (2 * (Math.Atan(TanZDemi * Mult))); Diff = LatiN - Lati1; Lati1 = LatiN; } while (Math.Abs(Diff) > 2.77777E-08); lat = (LatiN * 180) / Math.PI; lng = (Lambda * 180) / Math.PI; return new LatLon(lat, lng); }
我是你在帖子中提到的页面的作者。我不知道你是否解决了你的问题,但你给出的朗伯坐标是不正确的。我认为你必须把它们除以1000。得到x=151488.250和y=170492.909,这是可能的坐标,对应于…布鲁塞尔。在与lat/lng值进行转换时,要小心选择基准