실습 7. 하계망과 분수령

[1] 하계망 추출

흐름도(flow direction map) 만들기: ArcMap 구동하고 chapter7.mxd 파일 불러오기 -> Customize ~ Extensions ~ 3D Analyst와 Spatial Analyst 체크 -> Geoprocessing ~ Environments에서 Current Workpace는 chapter7 폴더 아래에 result 폴더를 생성하여 입력하고, Raster Analysis의 Cell Size는 Same as layer YDEM 선택 후 OK -> ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Flow direction 더블클릭하여 Input surface raster는 YDEM 선택, Output flow direction raster는 result 폴더 아래에 FlowDir_YDEM 입력 후 Save ~ OK

함몰지형(sink) 채우고 흐름도 다시 만들기: ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Sink 더블클릭하여 Input flow direction raster는 FlowDir_YDEM 선택, Output raster는 result 폴더에 sink로 입력 후 Save ~ OK (sink 레이어의 속성테이블 열어서 찾아낸 sink 확인) -> ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Fill 더블클릭하여 Input surface raster는 YDEM 선택, Output surface layer는 Fill_YDEM 으로 입력 후 OK -> ArcToolbox ~ Spatial Analyst Tools ~ Flow Direction 더블클릭하여 Input surface raster는 fill_ydem 선택, Output flow direction raster는 FlowDir_FYDEM으로 입력하고 OK

누적흐름도(flow accumulation map) 생성하기: ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Flow Accumulation 더블클릭하여 Input flow direction raster는 FlowDir_FYDEM 선택, Output accumulation raster는 Acc_FYDEM 입력 후 Save ~ OK

임계값 설정으로 하천 추출: ArcToolbox ~ Spatial Analyst Tools ~ Map Algebra ~ Raster Calculator 더블클릭 -> Con(“Acc_FYDEM" >= 3000, 1) 로 입력하고 Output raster는 result폴더에 stream으로 입력 후 OK -> 누적흐름량이 3,000 이상인 지역은 1의 값을 갖는 래스터로 변환됨을 확인

하계망 차수 계산: ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Stream order 더블클릭 후 Input stream raster는 stream 선택, Input flow direction raster는 FlowDir_룡뜨 선택, Output raster는 result 폴더에 StreamOr_ST 입력 후 OK -> 지류인 1차에서 본류인 4차까지 위계적으로 구분되어 나타남을 확인

래스터 하계망을 벡터로 변환: ArcToolbox ~ Conversion Tools ~ From Raster ~ Raster to Polyline 더블클릭 후 Input raster는 StreamOr_ST 선택, Output polyline features는 result 폴더에 stream.shp로 저장한 뒤 OK -> 벡터로 변환된 stream 레이어 우클릭 ~ Properties ~ Symboloby ~ Quantities ~ Graduated symbols 선택한 후 Fields Value에 GRID_CODE 선택하고 확인

유역분지(watershed)와 분수령 분석: ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Stream Link 더블클릭하여 Input stream raster는 stream 선택, Input flow direction raster는 FlowDir_FYDEM 선택, Output raster는 result 폴더에 st_link로 입력 후 OK -> ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Watershed 더블클릭하여 Input flow direction raster는 FlowDir_FYDEM 선택, Input raster or feature pour point data는 st_link 선택, Output raster는 result 폴더에 watershed로 입력 후 OK (생성된 유역분지 확인) -> Chapter7_성명_1.mxd로 저장

[2] 지형 특성 분석을 위한 고도-면적 곡선

Hypsometric Tools를 설치: ArcMap에서 새로운 빈 진도 생성 (File ~ New) -> www.esri.com 홈페이지에서 Support ~ Downloads ~ Arcscripts 메뉴에서 Hypsometry 검색한 뒤 Jerry Davis가 2010년 개발한 Hypsometric Tools를 다운로드 (AS16830.zip) -> 압축해제 후 _Hypsometry.tbx 파일을 ArcGIS 인스톨 경로(예: C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcToolbox\Toolboxes 폴더에 넣고, 압축해제 후 Scripts 폴더에 있는 모든 파일을 마찬가지로 ArcToolbox의 Scripts 폴더로 복사 -> ArcMap ~ Catalog ~ Toolboxes ~ System Toolboxes ~ _Hypsometry.tbx ~ Hypsometric 우클릭하고 Source 탭에서 Script 파일 경로를 수정 (내 경우는 경로 중간에서 Toolboxes\를 지워야 했음) -> 똑같이 Area Gradient, Area Gradient Highest Order, Get Max Watershed의 경로도 수정.

Hypsometric Tools를 이용한 고도-면적 곡선 산출: Add Data 눌러 Chapter7\map\ydem_clip.img 파일 추가 (Pyramid 구조 생성 No) -> ArcToolbox ~ Spatial Analyst Tools ~ Hydrology ~ Fill 더블클릭하여 Input surface raster는 ydem_clip 선택, Output surface layer는 clip_fill 으로 입력 후 OK -> ArcMap 오른쪽 Catalog ~ Toolboxes ~ System Toolboxes ~ Hypsometry.tbx ~ Hypsometric 더블클릭하여 watershed raster는 파일아이콘을 눌러서 result 폴더에 있는 clip_fill 선택 (아래화살표 아이콘 눌러서 선택하면 프로그램 버그로 디렉토리를 잘못 읽어서 실패하게 됨), Output raster는 result 폴더에 hypsometric 이라고 입력 후 OK -> ArcMap 메뉴 중 View ~ Graphs ~ Create Graph 눌러서 Graph type은 Scatter Plot 선택, Layer/Table은 hypsometric 선택, Y field는 PROPELEV 선택, X field는 PROPAREA 선택, Brush  심볼크기와 색상은 원하는 대로 변경하고 Next 버튼 누름 (용인시는 노년기 지형임을 알 수 있음) -> 그래프 Title 등을 적절히 적고 Finish 버튼 누름

최대 유역분지 추출: 오른쪽 Catalog ~ Toolboxes ~ System Toolboxes ~ _Hypsometry.tbx ~ Get Max Watershed 더블클릭 -> Input elevation raster는 clip_fill 선택, Output watershed elevation raster는 max로 입력 후 OK -> Chapter7_성명_2.mxd로 저장

[3] 결과물 제출

위 결과물 중에서 의미있는 내용들을 레이아웃으로 만들고 pdf로 출력하여 3장 이상을 게시판에 올릴 것. 끝.