QGIS EPSG proj 확인

421 views
Skip to first unread message

SUNGJIN KANG

unread,
May 27, 2015, 8:09:08 PM5/27/15
to osge...@googlegroups.com
QGIS 2.8.1을 사용하다가 좌표 변환 결과가 영 안 맞는 것 같아, 확인 좀 부탁드립니다.

EPSG : 5173, 5174, 5175, 5176, 5177에 모두 해당할 것 같습니다.

OSGeo 한국어지부 홈페이지에는 아래와 같이 proj가 게시되어 있습니다.

================================================================================================================================

*보정된 서부원점(Bessel) - KLIS에서 서부지역에 사용중

EPSG:5173

+proj=tmerc +lat_0=38 +lon_0=125.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


*보정된 중부원점(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


*보정된 제주원점(Bessel) - KLIS에서 제주지역에 사용중

EPSG:5175

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


*보정된 동부원점(Bessel) - KLIS에서 동부지역에 사용중

EPSG:5176

+proj=tmerc +lat_0=38 +lon_0=129.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


*보정된 동해(울릉)원점(Bessel) - KLIS에서 울릉지역에 사용중

EPSG:5177

+proj=tmerc +lat_0=38 +lon_0=131.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


================================================================================================================================


막상 QGIS에는 아래와 같이 적용되어 있습니다.


================================================================================================================================

EPSG:5173

+proj=tmerc +lat_0=38 +lon_0=125.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs


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


EPSG:5175

+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=550000 +ellps=bessel +units=m +no_defs


EPSG:5176

+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=550000 +ellps=bessel +units=m +no_defs


EPSG:5177

+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=550000 +ellps=bessel +units=m +no_defs


================================================================================================================================

아래 화면은 지적과 TMS Korea (Daum Satellite)를 중첩한 화면입니다.

왼쪽은 OSGeo에 있는 파라미터를 적용한 화면입니다.

아래 화면은 QGIS에 있는 파라미터를 적용한 화면입니다.

고수님들의 조언 부탁드립니다.



BJ Jang

unread,
May 27, 2015, 9:46:37 PM5/27/15
to osge...@googlegroups.com
안녕하세요? 장병진입니다.

먼저 결론부터 알려드리자면, http://OSGeo.kr/17 에 있는 좌표변환이 맞습니다.
QGIS에 있는 한국 좌표계들의 Proj4 인자에는 to_wgs84 인자가 빠져 있어 좌표계 변환만 되고 타원체 변환이 안됩니다.


QGIS, GeoServer 등에 들어있는 좌표계 변환을 위한 Proj4인자들은 EPSG에 등록된 정보를 바탕으로 하고 있습니다.
한국에서 많이 사용하는 좌표계들은 2009년(?) 정도에 국토지리정보원에서 EPSG에 등록했고, 이에 따라 많은 GIS 프로그램들이 한국좌표계를 다룰 수 있게 되었습니다.

그런데 국토지리정보원이 등록한 내용중 타원체 변환 부분의 방식이 좌표계 변환을 위해 가장 많이 사용하는 라이브러리인 Proj4에서 사용하는 좌표계 변환과 다른 문제가 있어 지금과 같은 문제들이 나오는 것입니다.

2006년에 발표된 [국토지리정보원 고시 제 2006-608호]로 국내측지계와 세계측지계 간의 변환, 즉 타원체 변환을 위한 Molodensky-Badekas 모델 인자가 고시되어 있습니다.
이 변환은 10개의 파라미터를 이용하는 방식이고, 이 방식으로 변환하는 인자가 EPSG에 등록되어 있습니다.

그런데 Proj4는 타원체 변환을 위해 Bursa-Wolf 모델을 사용하고 있고, 이 모델은 7개의 파라미터를 사용합니다. 10 파라미터와 7 파라미터를 사용하는 방식은 수학 모델이 달라 서로 호환이 안됩니다.
그래서 QGIS, GeoServer 등 Proj4를 좌표계 변환 라이브러리로 사용하는 프로그램들은 국토지리정보원이 EPSG에 등록한 좌표계변환 인자를 사용할 수 없어서 그냥 비워둬 버린 것입니다.


OSGeo 한국어지부에서도 해당 문제를 윤경담님의 제보로 2013년에 인지하게 되었고, 국토지리정보원에 문의해서 2006년 이전에 사용하던 Bursa-Wolf 모델을 위한 7 파라미터 인자를 받았습니다.
이 때 받은 인자로 http://OSGeo.kr/17 페이지의 내용이 수정되어 지금도 사용하고 있는 것입니다.

상세한 논의 과정은 다음 링크에서 보실 수 있습니다.


QGIS에서 to_wgs84 인자가 정상적으로 들어있는 좌표변환을 사용하시려면,
http://osgeo.kr/146 문서를 참고하시면 됩니다.
최신 QGIS와는 버전차이가 있어 그대로는 안되지만, 참고 자료로는 충분할 것입니다.

GeoServer에서 좌표계 변환을 정상적으로 하려면
data_dir/user_projections/epsg.properties 파일에 한국 좌표계들을 정상적인 값으로 잘 넣어주시면 됩니다.

PostGIS에서 한국 좌표계를 사용하기 하는 방법은
오픈소스 GIS 기초과정의 PostGIS 강의에서 매번 강의되고 있습니다.
윤정환님이 만드신 다음 자료의 26 페이지에 나와 있습니다.


그럼 수고하세요~ 





2015년 5월 28일 오전 9:09, SUNGJIN KANG <sjk...@gaia3d.com>님이 작성:

--
이 메일은 Google 그룹스 'OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter' 그룹에 가입한 분들에게 전송되는 메시지입니다.
이 그룹에서 탈퇴하고 더 이상 이메일을 받지 않으려면 osgeo-kr+u...@googlegroups.com에 이메일을 보내세요.
이 그룹에 게시하려면 osge...@googlegroups.com에 이메일을 보내세요.
http://groups.google.com/group/osgeo-kr에서 이 그룹을 방문하세요.
더 많은 옵션을 보려면 https://groups.google.com/d/optout을(를) 방문하세요.



--

Open Source GIS Technical Manager / e-Cartographer
장병 팀장 | BJ Jang

가이아쓰리디()
서울시 금천구 가산동 60-19 SJ테크노빌 702호
tel
: 02-3397-3475 | fax: 02-3397-3278
e-mailjang...@gmail.com | mobile: 010-3390-7056

Reply all
Reply to author
Forward
0 new messages