Primeiro problema:
Você está olhando uma mistura de mínimos. Uma árvore gigantesca com uma coroa do tamanho de um hectare parece bastante , interpretada com base na densidade de pontos / grãos, como um campo sem árvores. Você terá altos valores apenas onde houver árvores pequenas e de rápido crescimento, nas bordas e nas brechas da floresta. A parte complicada é que essas árvores menores e densas têm muito mais chances de serem obscurecidas por sombras ou oclusões ou serem resolvidas a uma resolução de 1 metro, ou aglomeradas juntas porque são um grupo da mesma espécie.
A resposta de Jen está correta nesta primeira parte: jogar fora as informações do polígono é um desperdício. Há uma complicação aqui, no entanto. As árvores de cultivo aberto têm uma copa muito menos vertical e mais espalhada, todas as outras coisas iguais, do que um estande de idade par ou uma árvore em uma floresta madura. Para mais, consulte o item 3.
Segundo problema:
Idealmente, você deve trabalhar com uma comparação de maçãs com maçãs. Confiar no NDVI para um e P&B para outro introduz um viés desconhecido nos seus resultados. Se você não conseguir obter dados adequados para 1989, poderá usar dados degradados em preto e branco para 2009, ou mesmo tentar medir o viés nos dados de 2009 em relação ao preto e branco e extrapolar os resultados do NDVI para 1989.
Pode ou não ser plausível abordar esse ponto do ponto de vista trabalhista, mas há uma chance decente de que isso seja apresentado em uma revisão por pares.
Terceiro problema:
O que exatamente você está tentando medir? A densidade do kernel não é um valor sem valormétrica, fornece uma maneira de encontrar áreas de árvores jovens que crescem rapidamente e que se matam rapidamente (sujeitas às limitações de sombreamento / oclusão acima); Somente aqueles com o melhor acesso à água / luz do sol, se houver, sobreviverão em alguns anos. A cobertura do dossel seria uma melhoria na densidade do kernel para a maioria das tarefas, mas isso também tem problemas: trata um grande conjunto de árvores de 20 anos de idade pareada que mal fecharam o dossel da mesma forma que uma planta estabelecida. floresta de anos. As florestas são difíceis de quantificar de maneira a preservar as informações; Um modelo de altura da copa é ideal para muitas tarefas, mas é impossível obter historicamente. A métrica que você usa é melhor escolhida com base na elaboração de seus objetivos. O que eles são?
Editar:
O objetivo é detectar a expansão do mato nas pastagens nativas. Os métodos estatísticos ainda são perfeitamente válidos aqui, eles exigem apenas algumas elaborações e escolhas subjetivas para serem aplicadas.
- Calcule uma medida básica da cobertura do dossel. Isso pode envolver uma abordagem em grade diretamente nos polígonos da coroa ou transformar os polígonos da coroa em um raster + desfocando-os se você precisar de uma versão mais contínua.
- Tente separar classes de paisagem nas quais você deve fazer sua análise, com base na porcentagem de cobertura do dossel. As técnicas estatísticas com as quais você trabalha na floresta fechada de dossel podem ser diferentes daquelas usadas em uma pastagem quase nua, ou podem até ser excluídas defensivamente da análise. Alguma pequena área de suas paisagens incluirá "expansão de matagal" e escolher como subdividir esse efeito e ignorar dados que não são relevantes depende de você como estatístico.
- Não sei se isso funcionará durante um período de tempo de 20 anos (e funcionará melhor com épocas intermediárias adicionais), mas tente prestar atenção ao diâmetro da coroa como um proxy para a idade das árvores. Há uma pergunta definitiva que você deve fazer, se o dobro do tamanho de uma coroa existente representa "expansão" ou se requer novas árvores. Se for o último, você tem alguma idéia de se eles são novos (pelo menos para algumas classes de paisagem que você selecionou acima, onde é possível verificar um certo grau de acesso à luz solar).
- Dependendo dos seus objetivos ecológicos, pode valer a pena não apenas explorar a densidade das árvores diretamente, mas explorar a fragmentação da paisagem usando pacotes como o Fragstats .
- Plano geral: verifique se não há um conjunto de dados LIDAR do condado aguardando para ser usado como validação e avaliação de precisão para sua capacidade de distinguir coroas no conjunto de dados de 2009.
