map_tile/tile/
coordinate.rs

1pub fn wgs842gcj02(lng: f64, lat: f64) -> (f64, f64) {
2    const SEMI_MAJOR_AXIS: f64 = 6378245.0;
3    const EE: f64 = 0.00669342162296594323;
4    const PI: f64 = std::f64::consts::PI;
5
6    let gcj02_x: f64;
7    let gcj02_y: f64;
8
9    let lng2: f64 = lng - 105.0;
10    let lat2: f64 = lat - 35.0;
11    let mut dlat: f64 = -100.0
12        + 2.0 * lng2
13        + 3.0 * lat2
14        + 0.2 * lat2 * lat2
15        + 0.1 * lng2 * lat2
16        + 0.2 * lng2.abs().sqrt()
17        + (20.0 * (6.0 * lng2 * PI).sin() + 20.0 * (2.0 * lng2 * PI).sin()) * 2.0 / 3.0
18        + (20.0 * (lat2 * PI).sin() + 40.0 * (lat2 / 3.0 * PI).sin()) * 2.0 / 3.0
19        + (160.0 * (lat2 / 12.0 * PI).sin() + 320.0 * (lat2 * PI / 30.0).sin()) * 2.0 / 3.0;
20    let mut dlng: f64 = 300.0
21        + lng2
22        + 2.0 * lat2
23        + 0.1 * lng2 * lng2
24        + 0.1 * lng2 * lat2
25        + 0.1 * (lng2.abs()).sqrt()
26        + (20.0 * (6.0 * lng2 * PI).sin() + 20.0 * (2.0 * lng2 * PI).sin()) * 2.0 / 3.0
27        + (20.0 * (lng2 * PI).sin() + 40.0 * (lng2 / 3.0 * PI).sin()) * 2.0 / 3.0
28        + (150.0 * (lng2 / 12.0 * PI).sin() + 300.0 * (lng2 / 30.0 * PI).sin()) * 2.0 / 3.0;
29
30    let radlat = lat / 180.0 * PI;
31    let mut magic = radlat.sin();
32    magic = 1.0 - EE * magic * magic;
33    let sqrtmagic = magic.sqrt();
34    dlat = (dlat * 180.0) / ((SEMI_MAJOR_AXIS * (1.0 - EE)) / (magic * sqrtmagic) * PI);
35    dlng = (dlng * 180.0) / (SEMI_MAJOR_AXIS / sqrtmagic * radlat.cos() * PI);
36    gcj02_y = lat + dlat;
37    gcj02_x = lng + dlng;
38    (gcj02_x, gcj02_y)
39}
40
41#[cfg(test)]
42mod tests {
43    use super::*;
44    #[test]
45    fn test_wgs842gcj02() {
46        assert_eq!(
47            wgs842gcj02(107.55892137719799, 48.048944178597864),
48            (107.56503861907252, 48.05044502209457)
49        );
50    }
51}
map_tile/tile/coordinate.rs

map_tile/tile/
coordinate.rs
