포인트 주변 버퍼 만들기
콩고 민주공화국 노르키부(North Kivu) 주의 주도인 고마(Goma)는 르완다와의 국경 근처에 있습니다. 좌표는 1.66°S 29.22°E이며, Point는 이미 UTM 좌표로 goma 변수에 제공되어 있어요.
고마에서 50km 이내에 몇 개의 채굴지가 있을까요? 그리고 국립공원 면적은 얼마나 포함될까요? 이를 버퍼 연산으로 확인해 봅시다. 거리는 CRS의 단위(이 경우 미터)로 표현해야 한다는 점을 기억하세요.
참고: 불리언 시리즈(예: 공간 관계 메서드의 결과)가 있다면, 그 불리언들의 합을 구해 True 값의 개수(즉, 조건을 만족한 지오메트리 수)를 계산할 수 있어요. 이때 True는 1, False는 0으로 처리됩니다.
이 연습은 강의의 일부입니다
Python으로 지리공간 데이터 다루기
연습 안내
goma를 중심으로 50km 버퍼를 만들고, 변수goma_buffer에 할당하세요.mining_sites의 포인트 중 이 버퍼 내부에 있는 것들을 확인하고(이를mask라고 하세요), 불리언mask의 합을 구해 채굴지 개수를 세고, 그 결과를 출력하세요.- 국립공원과 버퍼의 교차 영역을 계산하세요. 교차 영역들의 면적을 합산하고, 결과를 km² 단위로 출력하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# goma is a Point
print(type(goma))
# Create a buffer of 50km around Goma
goma_buffer = ____
# The buffer is a polygon
print(type(goma_buffer))
# Check how many sites are located within the buffer
mask = ____
print(____)
# Calculate the area of national park within the buffer
print(national_parks.____ / (1000**2))