Buffer ao redor de um ponto
Considere a cidade de Goma, capital da província de Kivu do Norte, no Congo, perto da fronteira com Ruanda. Suas coordenadas são 1,66°S 29,22°E (o Point já está fornecido em coordenadas UTM como a variável goma).
Quantos locais de mineração ficam a até 50 km de Goma? E qual é a área de parque nacional nessa região? Vamos descobrir usando a operação de buffer. Lembre-se de que as distâncias devem ser expressas na unidade do CRS (ou seja, em metros, neste caso).
Observação: se você tiver uma Series booleana (por exemplo, como resultado de um método de relacionamento espacial), pode calcular quantos valores True existem (isto é, quantas geometrias passaram na verificação) somando esses booleanos, porque, nesse caso, os valores True e False são interpretados como 1 e 0.
Este exercício faz parte do curso
Trabalhando com Dados Geoespaciais em Python
Instruções do exercício
- Crie um buffer de 50 km ao redor de
gomae atribua-o à variávelgoma_buffer. - Verifique quais dos pontos de
mining_sitesestão localizados dentro desse buffer (chame isso demask). Conte esses locais somando a máscara booleanamaske imprima esse resultado. - Calcule as interseções dos parques nacionais com o buffer. Some a área dessas interseções e imprima o resultado em km².
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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))