• 서울시 병원 및 약국 인허가 데이터 : http://data.seoul.go.kr/dataList/OA-16179/S/1/datasetView.do

  • 중부원점TM(EPSG:2097) 좌표계 (라고 소개하고 있으나, 실제로는 EPSG:5174 라고 하기도 한다.)
  • 문제1 : 만약 정말로 좌표계가 서울시 열린데이터 광장에서 명시하는 것과 다르다면, 경도의 상수값이 달라진다.
    • 서울시 열린데이터 광장에서 소개한대로 EPSG:2097로 하면 원점 longitude = 127
    • 만약 소개한 것과 달리 EPSG:5174 이면 longitude = 127.0028902777778
  • 문제2: EPSG:2097이든 EPSG:5174든 지구를 오래된 bessel 타원체로 계산한다.
    • 현재 우리의 좌표변환 코드는 GRS80타원체 (네이버, 카카오 등에서 사용) 기준의 공식으로 구현되어 있다. 이에 따라 오차가 생길 수 있다.
  • 그러나 응답 결과 오차가 그렇게 크지 않다.
  • 이 오차를 잡으려면 너무 엄한 곳까지 디깅해야 해서, 일단 여기에서 멈추고 추후 개선을 고려.

• 보정안되었고 오래된 지리원 표준, 실제로는 거의 사용되지 않는다고는 한다.

  • 중부원점(Bessel): 서울 등 중부지역

  • EPSG:2097

    +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs +towgs84=-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43

• 보정되었고 오래된 지루원 표준, 2002년 전까지 지리원과 국가 시스템에서 사용

  • 중부원점(Bessel): KLIS에서 중부지역에 사용중

  • EPSG:5174

    +proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs +towgs84=-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43

//https://www.youtube.com/watch?v=43hR79sRLIk
//위 영상을 그대로 따라갔고, 

export class Coordinates {
  constructor() {}

  coordination(lat: number, long: number): number[]{
    //GRS80 타원체 상수
    //장반경 (미터)
    const a = 6378137;
    //단반경 (미터)
    const b = 6356752.31425;
    //제1이심률 e*e = (a*a - b*b) / a*a = 0.00669438
    const e = 0.00669438;
    //제2이심률 d*d = (a*a - b*b) / b*b = 0.0067395
    const d = 0.0067395;
    //편평률 (a-b)/a = 1/298.257222101
    //중부원점 : log 127 lat 38
    //중부원점 Y축(East) 원점 가산값
    const dY = 200000;
    //중부원점 X축(North) 원점 가산값
    //놀랍게도 서울시 공공데이터 좌푠 Bessel1841 타원체 EPSG:2097
		//때문에 위 유튜브 영상과 달리, 가산값을 500000으로 수정
    const dX = 500000;
    //원점축척계수
    const k = 1;

    let T = Math.tan((lat * Math.PI) / 180) * Math.tan((lat * Math.PI) / 180);
    let C =
      (0.00669438 / (1 - 0.00669438)) *
      Math.cos((lat * Math.PI) / 180) *
      Math.cos((lat * Math.PI) / 180);
    let L = (long * Math.PI) / 180;
    let A = (L - 2.21656815) * Math.cos((lat * Math.PI) / 180);
    let N =
      a /
      Math.sqrt(
        1 -
          e * Math.sin((lat * Math.PI) / 180) * Math.sin((lat * Math.PI) / 180),
      );
    let M =
      a *
      (0.9983242984503243 * ((lat * Math.PI) / 180) -
        0.002514607064228144 * Math.sin((2 * lat * Math.PI) / 180) +
        0.000002639046602129982 * Math.sin((4 * lat * Math.PI) / 180) -
        3.418046101696858e-9 * Math.sin((6 * lat * Math.PI) / 180));
    let Mzero = 4207498.01927;

    /**
     * Y(E)
     *
     */
    let one = ((A * A * A) / 6) * (1 - T + C);

    let two =
      ((A * A * A * A * A) / 120) * (5 - 18 * T + T * T + 72 * C - 58 * d);

    let Y = dY + k * N * (A + one + two);
    console.log('Y:', Y);

    /**
     *
     * X(N)
     */
    let won = ((A * A * A * A) / 24) * (5 - T + 9 * C + 4 * C * C);

    let too =
      ((A * A * A * A * A * A) / 720) *
      (61 - 58 * T + T * T + 600 * C - 330 * d);

    let X = dX + k * (M - Mzero + N * Math.tan((lat * Math.PI) / 180) * ((A * A) / 2 + won + too));
    console.log('X:', X);

    let Xzero = Math.round(X)
    let Yzero = Math.round(Y)
    let Xone = Math.ceil(X - 1000);
    let Xtwo = Math.floor(X + 1000);
    let Yone = Math.ceil(Y - 1000);
    let Ytwo = Math.floor(Y + 1000);
 

    return [Xzero, Xone, Xtwo, Yzero, Yone, Ytwo];
  }
}