小程序地圖指定經緯度有偏移的問題

續 小程序讀取照片 EXIF 元信息

為了更直觀地查看照片的 GPS 位置，加了個 map 標簽，傳入 EXIF 中的經緯度。

結果偏成這個樣子：

還以為騰訊地圖有毛病，然后掏出 Google 一番雲雨之后，發現原來是因為不同的地圖用的坐標系不一樣，蘋果地圖用的坐標系是 WGS84，騰訊地圖用的 GCJ02。
所以才有這個偏移。

找到問題了，那就轉換唄：

// https://gist.github.com/jp1017/71bd0976287ce163c11a7cb963b04dd8
// const x_pi = 3.14159265358979324 * 3000.0 / 180.0
// π
const pi = 3.1415926535897932384626
// 長半軸
const a = 6378245.0
// 扁率
const ee = 0.00669342162296594323
const transformlat = (lng, lat) => {
  let ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng))
  ret += (20.0 * Math.sin(6.0 * lng * pi) + 20.0 * Math.sin(2.0 * lng * pi)) * 2.0 / 3.0
  ret += (20.0 * Math.sin(lat * pi) + 40.0 * Math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
  ret += (160.0 * Math.sin(lat / 12.0 * pi) + 320 * Math.sin(lat * pi / 30.0)) * 2.0 / 3.0
  return ret
}
const transformlng = (lng, lat) => {
  let ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng))
  ret += (20.0 * Math.sin(6.0 * lng * pi) + 20.0 * Math.sin(2.0 * lng * pi)) * 2.0 / 3.0
  ret += (20.0 * Math.sin(lng * pi) + 40.0 * Math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
  ret += (150.0 * Math.sin(lng / 12.0 * pi) + 300.0 * Math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
  return ret
}
// 地圖坐標系轉換 蘋果/高德是 wgs84 坐標系，騰訊是 gcj02 坐標系
const wgs84togcj02 = (lng, lat) => {
  let dlat = transformlat(lng - 105.0, lat - 35.0)
  let dlng = transformlng(lng - 105.0, lat - 35.0)
  let radlat = lat / 180.0 * pi
  let magic = Math.sin(radlat)
  magic = 1 - ee * magic * magic
  let sqrtmagic = Math.sqrt(magic)
  dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
  dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * pi)
  const mglat = lat + dlat
  const mglng = lng + dlng
  return { mglng, mglat }
}

const gcj02towgs84 = (lng, lat) => {
  let dlat = transformlat(lng - 105.0, lat - 35.0)
  let dlng = transformlng(lng - 105.0, lat - 35.0)
  let radlat = lat / 180.0 * pi
  let magic = Math.sin(radlat)
  magic = 1 - ee * magic * magic
  let sqrtmagic = Math.sqrt(magic)
  dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
  dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * pi)
  let mglat = lat - dlat
  let mglng = lng - dlng
  return { mglng, mglat }
}

最終效果，iPhone 掃一下這個碼看看：

Android相機用的啥坐標系，又是另一個話題了

A Short Guide To The Chinese Coordinate System