Quero calcular as porcentagens de área coberta (e.g. floresta, agropecuária, etc) para um shapefile, e baixar um csv já só com (1) o código identificador de polígono e (2) as colunas de % de cobertura (somando 1). A ideia é como um ‘zonal histogram’.
O tutorial nesse link ensina a recortar imagens para polígonos arbitrários, mas não para calcular as porcentagens de área em si. Existe algum script para GEE ou Python fácil de se adaptar?
Obrigado.
Oi, Ricardo. Obrigado pela pergunta.
Existe uma redutor no GEE para este tipo de análise: ee.Reducer.frequencyHistogram(). Você pode inserir este redutor como argumento na função image.reduceRegion(), usando a geometria de seus shapefiles no argumento ‘geometry’.
Atenciosamente,
Felipe
Oi Felipe, obrigado pela resposta! Talvez seja pedir demais, mas você tem algum exemplo de código completo usando essas funções?
Estou ainda começando com código de GEE, e não está claro como implementar as funções.
Oi, Ricardo.
Veja se este trecho pode ajudar:
https://code.earthengine.google.com/63aa35eee0ce72aacbde4c792c806398
É uma aproximação que deve funcionar para shapefiles pequenos (veja que usei quatro feições no exemplo).
O ideal é evitar ao máximo o uso de dados vetoriais em operações no GEE.
‘Zonal histogram’ pode ser executado no GEE usando layers referentes ao seu shapefile rasterizado (usando funções como .reduce() e .group() ).
Atenciosamente,
Felipe