[QGIS를 이용한 분석]_시작하기, 표고분석

안녕하세요.

QGIS를 이용한 분석을 정리해보려 합니다.

 

분석내용들에 앞서

QGIS 공식 버전을 설치합니다.

http://www.qgis.org/

QGIS 프로젝트에 오신 것을 환영합니다!

윈도우 사용자의 경우 컴퓨터 운영체체(64bit 32bit)에 맞는 설치 프로그램을 다운로드 받으시면 됩니다.

현재까지는 3.12버전까지 다운이 가능하네요!

저희는 3.8.0 버전을 사용했습니다.

 

그럼 QGIS를 이용한 표고분석에 대한 포스팅을 시작해보겠습니다.

① 레이어 정리

먼저 QGIS에서 cad파일을 불러올때는 dwg파일이 아닌 dxf파일이어야 합니다.

표고분석에 필요한 지형과 구역을 dxf파일로 각각 저장합니다.

QGIS에서 지형, 구역을 불러올 때 사용할 좌표계를 지정해 줍니다.

(Bessel 동부좌표를 사용했습니다. korean1985/modified east belt)

 

불러온 지형레이어와 구역레이어를 shp파일로 각각 내보냅니다.

(내보내기 - 다른이름으로 저장)

생성된 구역레이어는 폴리곤으로 변환합니다.

메뉴 - 벡터 - 도형도구 - "라인을 폴리곤으로"

폴리곤이 생성되며 이름은 "구역계"로 지정해줍니다.

 

② tin 생성(보간법)

다음은 지형레이어를 이용하여 tin 생성을 합니다.

공간처리toolbox - 보간법 - "tin보간"

보간에 z좌표값 이용 체크, 벡터레이어에 지형레이어로 지정, 유형은 구조선,

보간방법은 클러프토처(3차), 범위는 지형레이어 범위를 사용하고 픽셀크기는 10으로 설정합니다.

레어어에 보간출력물이라는 tin이 생성됩니다.

 

③ 래스터를 벡터로

래스터자료인 보간출력물 레이어를 벡터자료로 변환합니다.

메뉴 - 래스터 - 변환 - "폴리곤화"

하면 벡터화한 출력물이라는 레이어가 생성됩니다.

④ 구역추출

변환한 벡터레이어에서 구역계에 해당하는 자료만 추출해야합니다.

메뉴 - 벡터 - 공간처리도구 - "자르기"

입력레이어에 벡더화한출력물, 중첩레이어에 구역계 레이어를 지정합니다.

출력물이라는 레이어가 생성됩니다.

 

⑤ 속성값 분석

출력물의 속성테이블을 엽니다.

객체 수와 객체별 면적을 확인할 수 있습니다.

(면적을 확인 할 수 있을 뿐, 정확한 최소 최고값을 알기는 어렵습니다.)

속성테이블에서 편집모드로 전환합니다.

새 필드(area)를 만들어서 구역의 면적을 계산합니다.

(유형을 double로 지정하여 소수점 값이 나오도록)

필드계산기에서 기존필드를 갱신을 선택하고 , 계산할 필드를 지정해줍니다.

면적계산을 위해 $area를 더블클릭하면 표현식에 $area라고 자동 생성됩니다.

소수점 둘째자리까지 추출하기 위해서 표현식을 round($area,2)로 지정해 주었습니다.

면적이 계산되면 구역계의 구적면적과 비교해 봅니다.

(계산된 값의 합계 1,189.96㎡ , 구역계의 구적은 1,189.7㎡)

제대로 추출을 했다는 것을 알 수 있습니다.

 

⑥ 마스크 레이어로 래스터 자르기

최소 최대 값을 추출하기 위해 래스터자료인 보간 출력물을 구역계 레이어로 자릅니다.

메뉴 – 래스터 – 추출 – “마스크 레이어로 래스터 자르기”

전체 표고값은 최소값 0 최대값 212.995m

구역계 표고값은 최소값 109.978m 최대값 113.231m인 것을 알 수 있습니다.

 

⑦ 심볼조정

보간 출력물(tin) 레이어의 심볼을 조정해줍니다.

유형은 단일밴드 유사색상으로 설정해주고 원하는 색으로 지정해줍니다.

(심볼설정의 보간 유형을 불연속으로 지정했습니다.)

등간격, 등개수 등 원하는 간격으로 지정해주고 값을 설정해 줍니다.

구역 레이어의 심볼을 조정해줍니다.

채우기, 라인, 굵기 선스타일 등 원하는 스타일로 지정해줍니다.

(빨강-단순라인, 굵기0.5)

(지형레이어는 회색선으로 조정했습니다.)

 

⑧ 표고분석도 및 추출값 분석정리

(+심볼설정의 보간 유형을 선형으로 했을 때의 결과물입니다.)

전체면적은 1,189.96㎡이며

100m이상 110m미만에 속하는 면적 339.06㎡, 28.5%

110m이상 120m미만에 속하는 면적 850.9㎡, 71.5%

표고 최소값 109.98m

표고 최대값 113.23m

 

이상 표고분석을 마치겠습니다.

▼(최신) 좀 더 간편하게 표고분석 할 수 있는 방법을 정리하였으니 아래를 참고하시기 바랍니다.▼

Qgis 표고 분석방법(feat. 데이터 시각화, 입체적 보정 방법) (tistory.com)

Qgis 표고 분석방법(feat. 데이터 시각화, 입체적 보정 방법)