국가공간정보 포털에서 "연속지적도_서울_영등포구" 자료를 받아 Postgresql에 저장 하여 사용 하던 도중 배경지도와 이격이 발생하는 부분을 발견 하였다.
해당 좌표계는 Bessel EPSG:5174 이다.
Postgresql 에서 해당 좌표계를 조회해 보자.
SELECT * FROM SPATIAL_REF_SYS WHERE SRID = 5174;
srtext 컬럼값은
"PROJCS["Korean 1985 / Modified Central Belt",GEOGCS["Korean 1985",DATUM["Korean_Datum_1985",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6162"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01 (...)"
proj4text 컬럼값은
"+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs "
http://epsg.io/ 에서 확인해 보면 해당 좌표계에 문제는 없어 보였다.
이격이 발생한 이유는 이곳에서 찾을 수 있었다.
평행이동량(m) | 회전량(″) | 축척변화(ppm) | ||||
Δx | Δy | Δz | Rx | Ry | Rz | λ |
-115.80 | +474.99 | +674.11 | -1.16 | +2.31 | +1.63 | +6.43 |
+towgs84(타원체 변환 계수) 값을 추가하여 해당 좌표계를 추가 해야 한다.
국내 좌표계 towgs84 파라미터 값 계산에 대한 정보는 이곳을 참조하면 된다.https://groups.google.com/forum/#!topic/osgeo-kr/0oMKjzyLPW4타워체 변환 계수 값을 추가 하여 좌표계를 새로 생성해 보자.우선 사용할 srid 값과 srtext, proj4text 값을 살펴보면
+towgs84 추가 srtext
"PROJCS["Korean 1985 / Modified Central Belt",GEOGCS["Korean 1985",DATUM["Korean_Datum_1985",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43],AUTHORITY["EPSG","6162"]],PRIMEM["Gree (...)"+tpwgs84 추가 proj4text
좌표계를 생성 하고 해당 좌표계에 대하여 geometry 정보도 업데이트 시켜야 한다.
"SELECT UpdateGeometrySRID('table_name', 'geometry column', 900916)"
이후 새로 생성한 EPSG:900916 사용, 배경지도와 이격이 사라짐을 확인 하였다.
'꿀팁 - PostGIS' 카테고리의 다른 글
PostGIS 에서 Geometry 오류 수정 (3) | 2019.03.26 |
---|---|
수치지형도 등고선을 이용한 DEM 제작 방법 (2) | 2019.03.22 |
PostGIS에서 Hausdorff 거리를 이용하여 건물통합 건물과 도로명주소 건물 비교 (0) | 2019.03.13 |
PostGIS에서 Voronoi Daiagram을 이용한 경계 분할 (0) | 2019.03.12 |
PostGIS에서 공간 테이블 좌표변환 (0) | 2019.01.09 |