10. QGIS 3.4 강좌 16 스캔한 종이지도를 디지털지도로 제작하기
상단 래스터 메뉴에서 GDAL 지오레퍼런 보이지 않을 경우
상단 메뉴 플러그인 에서 GDAL 검색 후 체크
1. 스캔 종이지도에 구글 맵이나 현장 GPS 로 측정한 위치자료 표기
- 본 강좌는 구글에서 확인된 경위도 좌표 이용.
- 종이지도 불러오기
- 종이지도에서 사용되는 좌표계선택 : WGS84 =EPSG 4326
2. 종이지도에서 확인되는 지점의 위치좌표 등록
- 위치정확도가 좋은 최소 6개이상의 지점 등록
* 여기서 6개 이상 지점을 GCP 라고 함
GCP : Ground Control Point
3. 종이지도를 디지털지도로 변환
1) 상단메뉴 설정 > Transformation Settings > 변환유형 : 선형 , 대상 공간 참조 시스템 : EPSG:4326 - WGS 84
2) 정확도 확인 및 변환식 설정
데이터 확인후 상단메뉴 ▶ 클릭
4. 디지털지도를 불러오기
래스터 클릭 > "무등산기하보정_변화됨.tif" 으로 저장한 파일 불러오기
QGIS 는 오픈소스 GIS 프로그램입니다.
전에 소개했던 OSGeo4W 를 설치했다면 이미 깔려있습니다.
qgis 플러그인 중 GDAL 지오레퍼런서에 대해 예제를 통해서 살펴보겠습니다.
시작 > 모든프로그램 > OSGeo4W > qgis desktop 을 실행시키고
메뉴> 플러그인 > 플러그인 설치 및 관리... 를 누릅니다.
설치됨 항목에서 GDAL 지오레퍼런서를 체크해서 플러그인을 활성화시키면
메뉴 > 래스터 > 지오레퍼런서 메뉴가 생깁니다.
지오레퍼런서 메뉴 > 설정 > 지오레퍼런스 설정 에서 "ID 표시" 와 "좌표보기"를 체크.
1. 래스터 열기 아이콘을 눌러서 매칭할 이미지를 열면 좌표계 선택 창이 뜹니다.
2. 필터창에 4326을 입력하면
3. 세계의 좌표계 란에 WGS84(EPSG:4326)이 나옵니다. 선택하고 확인을 누릅니다.
1.포인트 추가 아이콘을 누르고
2. 좌표를 알고 있는 지점을 마우스 좌클릭하면
3. 지도 좌표 입력창이 뜨는데 x, y 값(경도, 위도)을 써줍니다.
아래가 좌표점을 하나 입력한 후 모습입니다.
지도 아래쪽 GCP 테이블에 포인트에 대한 값이 나옵니다.
id, srcX, srcY, defX, defY,..(포인트ID, 이미지X 픽셀값, 이미지Y픽셀값, 경도값, 위도값)
계속해서 총 16개의 포인트를 추가한 화면입니다.
이미지 픽셀값의 원점은 좌상단입니다.
꽤 번거로운 작업입니다.
그래서 GCP 테이블을 .points 파일로 저장하고 불러오는 기능이 있습니다.
아래는 예제에 쓰인 1:50,000 축척 지형도(33602_추자)에 해당하는 points 파일입니다.
지형도 도엽코드 규칙을 이용하면 각 지형도에 해당하는 points 파일을 따로 만들어서
좀더 쉽게 지오레퍼런싱을 할 수 있습니다.
mapX,mapY,pixelX,pixelY,enable
126.25,33.75,543.68037159763355248,-8381.70678816568033653,1
126.33333333299999879,33.75,2671.43259171597628665,-8396.54119289940899762,1
126.41666666700000121,33.75,4799.14044497041413706,-8409.49250177514841198,1
126.5,33.75,6926.38575147928986553,-8420.59362366863933858,1
126.25,33.83333333329999704,562.22660828402388233,-5834.97752426035458484,1
126.33333333299999879,33.83333333329999704,2687.62172781065100935,-5847.46628639053233201,1
126.41666666700000121,33.83333333329999704,4813.47939408284037199,-5860.41759526627174637,1
126.5,33.83333333329999704,6938.8745136094676127,-5873.36890414201207022,1
126.25,33.91666666670000296,581.19102485207122299,-3287.75280473372777124,1
126.33333333299999879,33.91666666670000296,2704.27341065088785399,-3300.24156686390506366,1
126.41666666700000121,33.91666666670000296,4828.28088994082827412,-3313.65542248520705471,1
126.5,33.91666666670000296,6951.36327573964445037,-3326.60673136094646907,1
126.25,34,599.23034792899420609,-740.06553846153838094,1
126.33333333299999879,34,2720.92509349112424388,-754.86703431952651044,1
126.41666666700000121,34,4842.15729230769193236,-768.28088994082827412,1
126.5,34,6963.85203786982219754,-778.45691834319518421,1
출처: https://gisngps.tistory.com/8 [GIS&GPS]
* QGIS 사용자 지침서 한글 문서 : https://docs.qgis.org/2.18/ko/docs/user_manual/introduction/getting_started.html
* 국가공간정보포털 : http://openapi.nsdi.go.kr/nsdi/index.do
댓글 ( 4)
댓글 남기기