안녕하세요? 이번 글에서는 법정구역정보를 확보하는 방법을 정리해 보겠습니다.
법정구역정보는 '국가공간정보포털 오픈마켓'에서 제공하고 있습니다.
국가공간정보포털 오픈마켓: 법정구역정보 | http://data.nsdi.go.kr/dataset/20171128ds00001
로그인 없이 바로 다운로드 받으실 수 있는데요, 압축을 해제하시면 시군구, 읍면동, 리 단위 SHP 파일이 들어 있습니다.
[1] 시군구(SIG)
[2] 읍면동(EMD)
[3] 리 단위(LIO)
위 레이어와 웹 지도를 중첩해볼까요?! '플러그인 > 플러그인 관리 및 설치'에서 'QuickMapServices' 플러그인을 설치해 보겠습니다.
상단 메뉴에서 '웹 > QuickMapServies > OSM > OSM Standard'를 선택해 보겠습니다.
자, 아래와 같이 법정구역정보와 웹 지도가 중첩되었습니다. 한 번 확대해볼까요?!
아래와 같이 법정구역정보와 웹 지도가 약 400m 이격되는 것을 볼 수 있습니다.
원인 파악을 위해 법정구역정보 레이어의 좌표계를 확인해 보겠습니다. 레이어 명을 우클릭한 후 'Set Layer CRS'를 선택합니다.
현재 정의된 레이어 좌표계는 EPSG:5174인데요, Proj4를 보시면 Bessel1841 타원체(+ellps=bessel)인 것을 확인하실 수 있습니다.
이 좌표계는 현재 Bessel 1841과 WGS84 사이의 타원체 변환계수(+towgs94)가 정의되어 있지 않아 위치 오차가 발생함을 확인합니다.
그렇다면 타원체 변환 계수를 포함하여 EPSG:5174 좌표계를 재정의해야 할 텐데요,
관련해서는 OSGeo 한국어지부 신상희 님께서 작성해주신 아래 글을 참조하시면 되겠습니다.
QGIS에서 타원체 변환 계수를 포함시켜 좌표계 설정하는 방법
출처: http://www.osgeo.kr/146 [OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter]
위 글에서 소개된 방법으로 계속 진행해 보겠습니다. '설정 > 사용자 지정 투영체'를 선택하면, '사용자 지정 좌표계 정의' 창이 실행됩니다.
'기존 좌표계에서 파라미터 복사'를 클릭하고,
'5174'를 검색하여 EPSG:5174를 선택해 줍니다.
아래와 같이 명칭은 '새 좌표계', 파라미터는 'EPSG:5174' 좌표계 정보가 추가되는데요,
파라미터에 Proj4에서 타원체 변환계수를 정의하는 '+wgs84' 정보를 추가해 주시면 됩니다.
Bessel1841과 WGS84 사이 타원체 변환계수는 '+towgs84=-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43'입니다.
명칭은 '5174WithWGSPara'로 변경해 주겠습니다.
타원체 변환계수값과 관련해서는 OSGeo 한국어지부 장병진 님이 작성하신 아래 글을 참조하시면 되겠습니다.
[OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter]
출처: http://www.osgeo.kr/17 [OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter]
이제 법정구역정보 레이어 좌표계를 아래와 같이 사용자 지정 좌표계로 변경해 보겠습니다.
이제 오차 없이 법정구역정보와 웹 지도가 제대로 중첩된 것을 보실 수 있습니다.
