Geógrafo, mestre em Ciências Ambientais e doutorando em Ecologia Aplicada (UFLA). Trabalho com ecologia de paisagem voltada para avaliação das consequências da fragmentação insular, proporcionada por barragens hidrelétricas na Amazônia, na biodiversidade.

cezar-fonseca@hotmail.com

exec

introducao_ao_r.txt script_exercicio_4.r script_exercicio_5.r

Uma das principais consequências da fragmentação de habitat é a promoção e/ou aumento do isolamento de manchas de habitat (fragmentos), ocasionando a perda e/ou redução da conectividade pela paisagem (UEZU & METZGER, 2016; PARDINI et al., 2010; MARTENSEN et al., 2008). Tal aspecto é importante por representar a possibilidade de fluxo gênico, manutenção de (meta)populações, dispersão de sementes, redução da probabilidade de extinções locais em paisagens fragmentadas (TAYLOR et al., 1993). No entanto, a conectividade pode ser representada por duas abordagens, estrutural que compreende a ligações físicas entre manchas de habitat, por exemplo através de corredores e stepping stones (METZGER & DECAMPS, 1997 ), e conectividade funcional que consiste no grau com que a paisagem facilita ou impede o movimento das espécies e recursos entre as manchas de habitat, ou seja, a capacidade de promover o fluxo biológico pela paisagem (TAYLOR et al.,1993).

No contexto de conectividade funcional, se compreende que a movimentação de organismos em paisagens fragmentadas é comum, sendo determinada por dois principais aspectos de percepção espécie-especifico, a habilidade de deslocamento por áreas de não habitat (matriz) e a resistência que a matriz inter-habitat proporciona. Assim, medidas de conectividade funcional são fundamentais para compreender dinâmicas de populações e comunidades que ocorrem em manchas de habitat de paisagens alteradas.

A quantificação da conectividade funcional baseada na capacidade de deslocamento das espécies pela matriz, se configura como uma variável preditora importante na determinação da riqueza e abundancia de espécies, e processos ecológicos, tais como polinização e dispersão de sementes, sendo mais um aspecto a ser avaliado nos planejamentos de paisagens e pesquisas de biodiversidade.

Para isso, abordagem da teoria dos grafos sustenta o método de medida da conectividade funcional, caracterizada pelo agrupamento de fragmentos em função de regras de ligação, que representam diferentes capacidades de deslocamento pela matriz (RIBEIRO et al., 2009; METZGER et al., 2009). Atualmente tal medida pode ser realizada em diferentes softwares de SIG, através de rotinas espaciais que produzem índices de conectividade, porém para paisagens muito grandes o processamento desta métrica fica custosa, demandando muito tempo ou não conseguindo executar a tarefa completamente. Assim, é interessante a criação de uma rotina automatizada que consiga gerar essa medida de forma mais rápida e por completa.

Neste contexto, minha proposta é elaborar uma função na linguagem R que quantifique a conectividade funcional, sendo está representada pelo agrupamento de fragmentos em razão de uma determinada distancia, no qual identificará fragmentos conectados, que sejam considerados como manchas funcionais conectadas. Para a execução da função será necessário como dado de entrada um arquivo de formato shapefile (.shp) composto por polígonos que representem as manchas de habitat do estudo em questão, além disso será necessário que o arquivo esteja projetado em UTM, devido aos cálculos de área e distancia que serão realizada pela função.

A leitura deste tipo de arquivo e as medidas espaciais a serem realizadas no R serão executadas por pacotes específicos de análises espaciais que deverão ser carregados antes da execução da mesma. As operações a serem realizadas serão criação de buffers, realização de dissolve e clip dos polígonos, além da soma da área dos fragmentos in situ e conectados. Como saída a função gerará um shapefile com a identificação dos fragmentos conectados e a soma da área conectada funcionalmente para cada distancia demandada.

Oi Vitor, é isso mesmo, você compreendeu correto. Porém, eu que fiquei na dúvida no uso do igraph, pois utilizar a menor distancias entre bordas é algo obrigatório, pois minha ideia é utilizar diversas distancias que o usuário defina 'aleatóriamente' e não que seja de característica da paisagem em análise. Segundo, essa lista de distancia seria a distancia de todos os fragmentos entre eles?? O resultado do igrap me daria os fragmentos conectados com o mesmo ID, para assim eu poder identificar as manchas conectadas funcionalmente??

A rotina de processamento que eu tinha elaborado era a seguinte: entrava com o shape de floresta, gerava um ID para cada fragmento, depois calculava a área de cada fragmento, criava buffer(s) nos polígonos, utilizando a metade da distancia desejada, por exemplo, para uma distancia de 50m, definiria um buffer de 25m, considerando que para os fragmentos se conectarem basta 25m de raio de cada fragmento. Destaco que o buffer criado deve considerar a área interna dos fragmentos, para posteriormente estas serem somadas e gerarem um valor único de área para a mancha funcional (conectada). Deste buffer, realizaria um clip com o shape de floresta para que as áreas do buffer externa aos fragmentos fossem desconsideradas, ficando apenas os fragmentos com sua forma real, porém agrupados em um único polígono, os que estiverem conectados. Após disso, individualizaria os polígonos, através de alguma rotina semelhante a função Explode da barra de edição do ArcGIS,e adicionaria um ID para cada polígono. Por fim, calcularia a área destes fragmentos para obter a área conectada funcionalmente.

Seria isso, espero ter sido claro.

Abraços

O objetivo desta função será de selecionar em uma paisagem com dimensão regional e/ou de bioma, paisagens menores que apresentem diferentes quantidades de floresta e que sejam escolhidas de maneira randômica. O produto desta operação seria um conjunto de paisagens que formassem um gradiente de cobertura florestal, onde possivelmente seriam realizados levantamentos, já considerando este aspecto espacial importante em ampla escala. Concomitante, poderia ser criadas paisagens de diferentes tamanhos, considerando o efeito da escala para a seleção destas. Neste sentido, os dados de entrada seriam o dado espacial de classe de floresta projetado em UTM e a definição dos tamanhos de paisagens a serem criadas.

Oi Vitor, Obrigado pelas sugestões, modifiquei o texto , veja se está mais claro agora, valeu. Abraços!!

Olá Professor, grato pelas considerações, detalhei mais o que a função se propõe a fazer, veja se está mais compreensível agora. Algumas das etapas já consegui fazer, como carregar o shape, calcular área, fazer buffer e dissolve, no momento estou pesquisando como executar as demais. Abraços!!