[QGIS]기상관측지점 위치 지도 시각화/ASOS/값 필터/라벨 수식

종관기상관측지점/ASOS

방재기상관측지점/AWS

기상관측지점 / 기상관측

QGIS 필터 / QGIS 심볼

객체 필터/ 라벨 수식/ 라벨 표현식

QGIS 라벨 붙이기/라벨 추가/ 라벨 중복

-----

   2022.03.05 메모   

- 관측지점 데이터에서 종료된 지점도 포함해서 분석을 진행하여

  최종 지도 내 갯수가 위치 중복이 있을 수 있습니다..

- 추후 수정사항 반영해서 재업로드 하겠습니다..

- 꼼꼼히 확인하지 못해 번거롭게 해드려 죄송하며, 

  작성글 보실 때 이 점 경계하셔서 보시길 바랍니다.

-----

 

기상자료개방포털에서 제공하는 데이터의 관측지점 정보를

실제 지도 위에서 시각화해보겠습니다.

 

현재 관측을 수행하고 있는 종관기상관측지점은 전국에 총 95개입니다.

 

먼저 데이터는 기상자료개방포털에서 전국 관측지점 정보를 다운받았습니다.

https://data.kma.go.kr/tmeta/stn/selectStnList.do?pgmNo=123 

기상자료개방포털[데이터:메타데이터:지점정보:관측지점정보:자료]

이후, 이전 포스팅인 '텍스트 나누기'에서 시도 정보를 추출하였습니다.

https://garden-aibigdata.tistory.com/15

[엑셀]텍스트 원하는 대로 나누기/분리하기/공백단위/띄어쓰기/기상관측 지점 시도 추출(ASOS)

 

 

---

1. 위/경도 정보가 있는 기상관측지점 데이터를 QGIS에서 불러옵니다.

   QGIS에서 '레이어' - '레이어 추가' - '구분자로 분리된 텍스트 레이어 추가'

 

  꼭 .csv 파일형태로 저장하셔야 불러올 수 있습니다!!!

---

2. 파일이름에 관측지점정보 데이터 파일을 넣고, 좌표를 설정하고 추가합니다.

 

⭐⭐⭐ 중요 ⭐⭐⭐

   X필드 ▶ 경도   

   Y필드 ▶ 위도 

⭐⭐⭐⭐⭐⭐⭐⭐

---

3. 종관기상관측지점 위치 데이터를 추가하였습니다. 이제 좀 더 잘 볼 수 있도록 시각화하겠습니다.

---

4. 읍면동 별 관측지점의 갯수를 파악하기 위해 읍면동 행정구역 shp파일을 추가하였습니다.

읍면동 행정구역 파일은 아래의 URL들에서 다운받을 수 있습니다.

http://www.gisdeveloper.co.kr/?p=2332 

대한민국 최신 행정구역(SHP) 다운로드 – GIS Developer

http://data.nsdi.go.kr/dataset/15145

행정구역_읍면동(법정동) - 오픈마켓

 

---

5. 폴리곤 별 지점 개수를 알고 싶으므로

 '벡터' - '분석 도구' - '폴리곤에 포함하는 포인트 개수 계산' 을 실행합니다.

---

6. 폴리곤 ▶ 읍면동 행정구역 레이어  /  포인트 ▶ 종관기상관측(ASOS)지점정보 레이어

 

 단순히 공간 내 포인트 개수를 알고 싶을 때는 가중치 필드를 비워둡니다!!! 

 

---

7. 계산된 레이어의 속성으로 들어갑니다.

 

---

8. '심볼'로 들어가서 '단계 구분'을 선택, 

   폴리곤 별로 포인트 개수가 계산된 필드인 'NUMPOINTS'를 값으로 선택합니다.

---

9. 라벨을 달아보니 전국에 지점이 95개라 읍면동 단위에서는 0개인 곳이 많아서 보기가 불편했습니다.

그래서 객체 필터링을 통해 걸러내려합니다.

 

---

10. 마우스 우클릭을 통해 속성 테이블을 열어봅니다.

 

---

 

11. 위에 '서식을 이용해서 객체 선택/필터링'을 클릭합니다.

---

12. NUMPOINTS 필드에서 조건에 '초과'를 선택하고 0을 입력합니다.

이러면  NUMPOINTS > 0 인 결과만 선택됩니다. 

---

13. 우측 하단에 있는 '객체 선택'을 클릭하면

아래와 같이 113개의 일치하는 객체를 찾았다고 나옵니다! 

 

0개인 읍면동을 제외한 읍면동을 선택하였습니다.

(노란색 부분입니다)

 

14. 선택한 레이어만 저장하기 위해, '내보내기' - '선택한 객체를 다른 이름으로 저장'을 수행합니다.

---

15. 이제  한 눈에 지점 정보를 파악할 수 있도록 설정합니다.

 

* 심볼에서 '단계 구분'을 해서 읍면동 별 포함하고 있는 지점 개수를 구분합니다.

 

* 라벨에서 '단일 라벨'을 통해 포함된 관측 지점 개수를 표시합니다.

 

---

시도, 읍면동, 관측지점 개수를 한 번에 보고 싶다면?

라벨에서 수식을 적용해서 할 수 있습니다.

해당 관련 포스팅은 다음에 자세히 올릴 예정이므로, 지금은 간략하게 설명하겠습니다.

---

16. 값 오른쪽에 있는 해당 버튼을 클릭하여 '표현식 대화창'으로 들어갑니다.

 

---

17. 여기서 중복으로 라벨표시를 원하는 필드를 선택합니다. 

필드를 선택하고 오른쪽에 있는 '10개 샘플'을 클릭하면 값을 일부 보여주므로 이를 통해 선택하면됩니다.

 

저는 '읍면동' + '관측지점 수'를 붙이고 싶은데

줄 바꿈으로 하고 싶어서 아래와 같이 표현식을 입력하였습니다.

근데 결과가 생각보다 별로네요 ㅎ,,,

---

시도 구분을 해서 시각화 하면?

> 시도 행정구역 shp를 추가하여 아래와 같이 설정한 뒤, 

시도 폴리곤 기준으로 포인트 개수를 계산하여 추가하였습니다.

(강원도가 잘렸어요,,,ㅜ 데이터 안에 강원도 정보 없는거 아닙니다...ㅜ)

 

 

 

 

---

사실,, 시도별로도 폴리곤 포인트 개수 구해서 시각화 하였는데,,,,

캡쳐하려는데 다 날라가서,,,

여기까지만 작성하겠습니다...ㅜ

 

다음에 다른 좋은 정보로 포스팅하겠습니다-!

 

오류 또는 보완, 수정사항 있으시면 댓글 부탁드립니다.

감사합니다😉