흑백과 디지털

 내가 다시 GIS 에 대해 관심을 가진 이유는 많은 데이터가 공개되었기 때문이다. 아래 사이트에서 데이터들을 얻을 수 있다. 

좌표계의 경우 제가 임의로 한 것이다. 

* 기초구역도
http://www.juso.go.kr/addrlink/addressBuildDevNew.do?menu=bsin
"기초구역 관리번호" 라고 나오는데 기존의 우편번호를 생각하면 된다. 우편번호 시스템이 6자리에서 5자리로 변경되었다. 
현재 행정동이 없어지면서 구역별로 위치를 잡는 기준이 없어졌다. 새롭게 사용할 수 있는 기준이 기초구역도이다. 
- 좌표계 :  UTM-K (GRS80)    +proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=GRS80 +units=m +no_defs


* 도로명주소 배경지도
http://www.juso.go.kr/addrlink/addressBuildDevNew.do?menu=layer
구글맵이나 다음맵 같은 경우, 지하철 노선이라든지, 강이라든지, 산 부분을 그려주는 정보이다.
- 좌표계 :  UTM-K (GRS80)    +proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=GRS80 +units=m +no_defs

* 민원행정기관 전자지도
http://www.juso.go.kr/addrlink/addressBuildDevNew.do?menu=fclts
학교, 경찰서 같은 공공기간 좌표 정보이다. 구역정보가 아니라 점 정보이므로 포인트를 표시한다.
- 좌표계 :  UTM-K (GRS80)    +proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=GRS80 +units=m +no_defs

* GIS건물정보
http://openapi.nsdi.go.kr/nsdi/index.do 에서 "GIS건물일반집합정보"를 찾는다. 
건물을 표시한다. 너무 많은 정보가 있어 지역별로 따로 파일이 있다. 
- EPSG:5174 라고 되어 있어 보정된 중부원점(Bessel)를 사용한다.
  +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
기존 중부원점을 사용하니 daump map 이나 openstreetmap 에 비해 조금 좌표가 어긋나는 문제가 있었다. 

* openstreetmap
https://download.geofabrik.de/asia/south-korea.html 
오픈스트리트맵에 있는 정보를 다운 받을 수 있다.  2018년 10월 기준으로 174MB 정도된다. 
EPSG:4326, WGS 84

 좌표계에 대한 기본 설명은 

http://www.biz-gis.com/index.php?mid=pds&document_srl=65754  를 참조하기 바랍니다.

약간의 설명을 추가하면 
지오이드 :  지구상에서 높이(해발고도)를 측정하는 기준이 되는 가상면이다. 이 면은 중력 퍼텐셜이 같은 등퍼텐셜면이고,[1] 중력 가속도를 측정할 때 기준면이 된다. (https://ko.wikipedia.org/wiki/%EC%A7%80%EC%98%A4%EC%9D%B4%EB%93%9C)
당연하지만 중력은 질량에 따라서 달라진다. 지구 내부의 밀도가 다르기 때문에 중력 가속도를 기준으로 하는 선은 엄청 찌그러진 모양일 것이다. 따라서  이것을 모델로 하면 땅의 자표를 계산하기 어려울 것이다. 

지구 타원체 :  지구의 모양에 가까운 회전타원체
(https://ko.wikipedia.org/wiki/%EC%A7%80%EA%B5%AC_%ED%83%80%EC%9B%90%EC%B2%B4) 

준거타원체 : 지역에 따라 지오이드면과의 차이가 가장 작도록 지오이드면에 근접하는 타원체를 구축하고자 하는데, 현실의 특정 지점에 지구타원체를 접합시킨 수학적 모델을 준거타원체라고 함
(http://gis.seoul.go.kr/GisWebDataStore/Gis_Edu/html/S0602/SGIS-HTML.jsp?sgis=0602&pgis=0201)
지구 타원체의 경우 중심이 x,y 축의 중심이나 준거타원체의 중심이 원점이 아니다.

데이텀 : 지구의 크기와 형태, 좌표체계의 기준점과 방향 등을 정의해 놓은 것
(http://gis.seoul.go.kr/GisWebDataStore/Gis_Edu/html/S0602/SGIS-HTML.jsp?sgis=0602&pgis=0203)

이런 좌표계들은 변환이 가능합니다. 이를 쉽게 해주는 툴이 proj4(https://proj4.org/index.html) 입니다. 여기에 대한 자세한 설정값은 https://proj4.org/usage/quickstart.html 에서 확인 할 수 있습니다.

자주 사용하는 변환은 http://www.osgeo.kr/17  를 참조하기 바랍니다.

간혹 자료를 받다보면 확장자 .prj 파일이 있는데, 이 파일이 좌표계를 표시하는 파일입니다. QGIS가 자동으로 이 파일을 읽기도 하는데, 읽지 못하는 경우 수동 설정이 필요합니다. 

 내가 회사업무 외적으로 관심을 가지는 분야가 좀 많다. TTF 폰트, 이미지 프로세싱, 한국어 형태소 분석기, 딥러닝 이런  것들을 좋아한다. 뭔가 컴퓨터 소프트웨어를 이용하기는 하지만 약간 기술적으로 다른 분야의 접목이 일어나는 분야를 좀 좋아하는 경향이 있다. 이런 것중에 요근래에 관심 가졌던 분야가 있다. 지리정보시스템이다. 위에서 열거한 분야들은 여러번 시도 했으나 뭔가 눈에 보이는 성과를 내지 못해 포기했다가 또 재미를 느껴서 다시 시도 하는 것을 반복하게 된다. 반복할 때마다 기술이 발전해서 더 쉬운 방법이 생기는 경우도 있고, 내 이해력이 좋아져서 발전하는 부분도 있다. 어째든 계속 발전하는 재미가 있다. 지리정보 분야의 경우, 정부에서 공개되는 데이터가 생김에 따라서 뭔가 이용해 볼 수 있는 가능성이 높아져서 요즘 재미를 느끼고 있다. 

 나 같은 사람은 눈에 보이는 게 없으면 재미가 없다. 우선 QGIS 부터 설치해 보자. 
https://qgis.org/ko/site/forusers/download.html
에서 QGIS 독립 설치관리자 2.X 버전을 설치하도록 하자. 3.X 버전도 있지만 2.X 를 설치하는 이유는 다음맵이나 네이버맵 플러그인  설치가 간편해서 이다. 다행히 이 프로그램은 무료이다. 그리고 윈도우즈, 맥 모두 지원한다. 

설치후 QGIS Desktop 프로그램을 실행하면 플러그인 부터 설치하도록 하자. 플러그인 설치는
http://www.qgistutorials.com/ko/docs/using_plugins.html  를 참조해서 설치하면 된다. 설치할 플로그인은 TMS for Korea 이다. 

그리고 메뉴에서 프로젝트(J)-> 새로만들기  로 프로젝트를 새로 만든다. 
메뉴에서 웹(W)->TMS for Korea ->Daum Map->아무거나 를 선택하면 화면에 아래와 같이 맵이 나온다.

왼쪽 화면 아래에 Layers Panel 이 있다. 이 Layer 패널은 포토샵 레이어 처럼 체크박스로 show, hide 시킬 수 있다. 원하는 Layer 를 추가할 수 있다.

이제, 인터넷에 공개된 정보를 다운받아서 새로운 Layer 를 추가해보자. 
국가공간정보포털(http://openapi.nsdi.go.kr/nsdi/index.do) 사이트에 가면 공개된 데이터가 많다. 이 중 SHP 파일들을 QGIS 에서 볼 수 있다. 
부동산 개방데이터에 있는 법정구역정보를 선택해 보자. 그리고 가장 최신 파일중 전체데이터 다운 받는다. 그리고 맨 아래

배포 데이터 좌표계 : Bessel/TM, EPSG:5174 이런 말이 들어 있다. 이 말을 기억해 두자. 좌표계에 대해서는 
( http://www.osgeo.kr/17 ) 를 읽어보자. 

어째든 압축을 해제하고 (SIG) (시구의 약자인듯) 파일을 다시 압축해제 한다. 그리고 확장자 SHP 파일을 QGIS 에 drag and drop 한다. 

이렇게 안 나왔다면 좌표체계가 잘못된 것이다. 아래 글을 계속 읽도록 하자. 

이런식으로 도형이 그려진다. 그리고 마우스 휠로 확대를 할 경우 두 layer 가 확대 되어 나온다. 

Layers Panel 에서 AL_00_D001_20180901(SIG) 레이어를 선택하고 마우스 오른쪽 버튼을 누른 후, 속성을 선택하면 레이어 속서을 변경할 수 있다. 여기서 좌표체계도 변경할 수 있다. 아마 자동으로 잘 선택 될 거라고 생각되는데, 잘 선택이 안되어 있다면 콤보박스 옆 아이콘을 눌러서 EPSG:5174 를 검색해서 추가한 후 수정하도록 하자.

 구역 정보만 있어서 뭔지 알 수 없으니 라벨을 추가하자. 레이어 속성창에서 라벨탭을 선택하고 Show Labels for this layer 를 선택하고 Label with 를 A2를 선택하자. 

위와 같이 설정하면 Label 값을 아래와 같이 출력할 수 있다.

더 자세한 사용법은 https://docs.qgis.org/2.18/ko/docs/user_manual/  메뉴얼을 확인하자. 고맙게도 한국어로 번역되어 있다.