상황
EPSG:5174 관련해서 어긋남의 문제가 발견되었으며, QGIS, ArcGIS 에서의 결과가 상이하다.
proj의 +towgs84 옵션을 ArcGIS가 해석하지 못하는 듯 하다.
설명
1. 검은색은은 EPSG:5174 전국 DEM 이다.
2. 적색라인은 EPSG:4326 시군구 shapefile 이다.
3. 배경실사는 EPSG:5181 다음(daum.net) 위성 사진이다.
4. 캡쳐된 이미지는, 상기된 세 가지를 EPSG:5181투영하여 시각화 한 결과이다.
이유
EPSG:5174가 보정 안된 좌표계를 사용하고 있다.(참조)
결론
GDAL 유틸리티를 이용하여, EPSG:5174 보정치를 적용한다.
"보정된EPSG:5174" 를 "보정안된EPSG:5174"로 직접 변환 시, 문제를 해결할 수 없었다.
따라서 아래와 같이 "보정된EPSG:5174" 를 "EPSG:4326"으로 변환하고, "보정안된EPSG:5174"로 변환하는 방법을 취한다.
gdalwarp -overwrite -co compress=deflate -co tiled=yes -s_srs "+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" -t_srs "+proj=longlat +datum=WGS84 +no_defs" sk_dem_new_f dem.4326.tiff
gdalwarp -overwrite -co compress=deflate -co tiled=yes -t_srs "+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs" dem.4326.tiff dem.5174.tiff