EPSG.io를 통한 proj4 문자열 얻기

EPSG는 European Petroleum Survey Group의 약자로 전세계의 다양한 좌표계를 표준 코드화하여 최대한 많은 정보를 압축하여 제공하는 방식입니다. EPSG의 약자에 Survey, 즉 측량이 있다는 것은 GIS의 근본이 되는 좌표와 그 좌표가 측량으로 측정되고 존재한다는 의미로 측량의 중요성을 나타낸다고 할 수 있습니다.

이러한 EPSG에 대한 정보는 epsg.io라는 사이트를 통해 코드값으로 질의(Query) 함으로써 얻을 수 있습니다 예를 들어, 우리나라의 GRS80 타원체의 UTM-K 좌표계인 EPSG:5179에 대한 정보를 얻기 위해 아래와 같은 url을 호출함으로써, 그 결과를 JSON으로 얻을 수 있습니다.

https://epsg.io/?format=json&q=5179

결과는 아래와 같습니다.

{  
    "status":"ok",
    "number_result":1,
    "results":[  
       {  
          "code":"5179",
          "kind":"CRS-PROJCRS",
          "bbox":[  
             40.27,
             122.71,
             28.6,
             134.28
          ],
          "wkt":"PROJCS[\"Korea 2000 / Unified CS\",GEOGCS[\"Korea 2000\",DATUM[\"Geocentric_datum_of_Korea\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6737\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4737\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",38],PARAMETER[\"central_meridian\",127.5],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",1000000],PARAMETER[\"false_northing\",2000000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AUTHORITY[\"EPSG\",\"5179\"]]",
          "unit":"metre",
          "proj4":"+proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs",
          "name":"Korea 2000 / Unified CS",
          "area":"Republic of Korea (South Korea) - onshore and offshore.",
          "default_trans":15831,
          "trans":[  
             15831
          ],
          "accuracy":1.0
       }
    ]
 }

위의 내용중 상대적으로 의미가 있는 부분은 좌표계의 단위를 나타내는 unit과 해당 좌표계의 국제 명칭인 name과 좌표계가 허용하는 지역의 범위인 bbox와 좌표계 변환을 위한 proj4 문자열인 proj4입니다. 누락된 부분은 타원체간의 변환계수인데요. 위의 경우 proj4의 문자열 값 중 towgs84에 해당되는 부분입니다. 최근 타원체가 GRS80이나 이와 거의 동일한 WGS84가 주로 사용되고, Bessel1841과 같은 타원체는 거의 사용되지 않음으로써 타원체간 변환계수의 중요성은 시간이 갈수록 점점 사라질 것입니다. 그럼에도 아직 우리나라의 GIS 데이터는 Bessel1841 타원체를 근간으로 하는 좌표계가 적용되어 있는 경우가 많습니다. 대표적으로 지적도가 그렇습니다.

여튼, EPSG.io라는 사이트의 서비스를 이용하여 쉽게 좌표계의 제원 정보를 조회할 수 있음으로, 이를 활용하여 좌표계간의 다양성을 해결할 수 있습니다. 예를들어 OpenLayers에는 원하는 좌표계를 등록하여 서로 다른 좌표계간의 좌표 변환이 가능한데, 그 예가 아래와 같습니다.

import proj4 from 'proj4';
import {register} from 'ol/proj/proj4.js';
import {get as getProjection, getTransform} from 'ol/proj.js';

var newProjCode = 'EPSG:5179';

proj4.defs(newProjCode, proj4def);
register(proj4);

var newProj = getProjection(newProjCode);
var fromLonLat = getTransform('EPSG:4326', newProj);

위의 proj4def는 5179 코드로 하여 epsg.io의 서비스를 통해 얻은 proj4의 문자열 값입니다. 최종적으로 formLonLat라는 함수를 얻을 수 있는데, 이함수는 EPSG:4326 좌표계를 EPSG:5179 좌표계로 변환할 수 있는 함수입니다. 참고로 EPSG:4326은 WGS84 타원체의 경위도 좌표계입니다.

답글 남기기 응답 취소