[QGIS] 경사도 tif 와 행정동 shp 데이터 활용하기

 

활용 데이터

 

여기서 경사도 데이터를 수집했는데, tif 경남 지역 전체말고 부산의 경사도만 잘라내서 사용하고 싶다.

https://www.vworld.kr/dtmk/dtmk_ntads_s002.do?dsId=30017

브이월드

 

 

 

 

 

 

행정동 SHP에서 부산광역시만 잘라내기

 

어떻게 해야 부산의 경사도만 잘라서 사용할 수 있을까~ 

 

 

 

1. 전체 행정동 SHP에서 부산만 선택 - Edit → Select → Select Features by Expression 

2. 부산만 잘라내기 - Raster → Extraction → Clip Raster by Mask Layer

• Input layer: slope_gn
• Mask layer: BND_ADM_DONG_PG
• "Selected features only" 체크 켜기 ← 필수 (이게 206개만 마스크로 씀)
• Run

3. 점으로 변환 Processing Toolbox → Raster pixels to points → Input은 방금 잘라낸 부산 래스터 → Run

4. 좌표 컬럼 추가 Add geometry attributes → Input은 3번 점 레이어 → Run (xcoord, ycoord 생김)

5. CSV로 내보내기 4번 결과 레이어 우클릭 → Export → Save Features As → Format을 CSV → 저장

 

 

 

 

부산의 시도 코드는 21!

BND_ADM_DONG_PG CSV에서 ADM_CD 21로 시작한다는 것을 알 수 있다.

 

 

 

 

 

 

 

 

메뉴에서 Edit > Select > Select Features by Expression

 

 

 

아래 구문으로 했더니 먹히지 않는다...ADM_CD가 숫자가 아닌 문자열인걸까?

 

left(trim(to_string("ADM_CD")), 2) = '21' 을 입력하고 Select Features를 누르면

206개가 나온다!

 

 

 

 

 

 

 

 

 

Raster > Extraction > Clip raster by mask layer

 

 

 

 

 

 

 

3번 점으로 변환하고 생긴 Vector Points를

Add geometry attributes을 Input layer로 넣을것! Run

 

Encoding; UTF-8 선택 (이렇게 해야 한글이 있으면 깨지지 않음)

Geometry type; No Geometry 선택

 

만약 컴퓨터 경로를 설정 안 하고 냅다 파일명만 작성하면 아래와 같은 에러가 뜨니 유의하자! 

 

 

 

 

 

 

결과를 확인해보니, 경사도(value)와 좌표만 나온 것을 알 수 있다. 기존 행정동 데이터에는 시군구 주소가 있었기에 같이 출력하면 보기 편할 것 같다는 생각을 했다.

 

점 레이어와 행정동 경계를 결합하여 CSV 추출하기