ANÁLISE TEMPORAL DA COBERTURA DA TERRA NA ESEC RASO DA CATARINA

ATRAVÉS DO USO DE IMAGENS DE SATÉLITE

Uldérico Rios Oliveira1, Patrícia Lustosa Brito2, Mauro José Alixandrini Júnior2, Júlio César
Pedrassoli2 e Ricardo Lustosa Brito3
1
Universidade do Estado da Bahia - UNEB, Barreiras – BA, Brasil, [email protected]; 2Universidade Federal da
Bahia - UFBA, Programa de Pós Graduação em Engenharia Civil – PPEC, Departamento de Engenharia de Transporte e
Geodésia – DETG, Salvador-BA, Brasil. [email protected], [email protected], [email protected];
3
UFBA, Programa de Pós Graduação Ciência Animal nos Trópicos - PPGCAT, Salvador-BA, Brasil.
[email protected].

RESUMO 1. INTRODUÇÃO

O objetivo deste trabalho é avaliar mudanças ocorridas na A degradação do meio ambiente e o uso não sustentável dos
cobertura da terra na região da Estação Ecológica - Esec recursos naturais têm sido algumas das preocupações
Raso da Catarina nos últimos anos, a partir de imagens do recorrentes de vários cientistas no mundo inteiro [1]. Estes
satélite. A Esec Raso da Catarina é uma Unidade de temas, de grande abrangência e complexidade, necessitam
Conservação Federal de proteção integral, localizada entre cada vez mais da adoção de metodologias de monitoramento
os municípios de Jeremoabo, Paulo Afonso e Rodelas, sistemático e sinóptico.
Bahia, Brasil. Foram aplicadas técnicas de processamento As áreas da região da Estação Ecológica - Esec Raso da
digital de imagens em cenas Landsat 5, 7 e 8, que distavam Catarina vêm sofrendo atividades conflitantes no seu uso e
temporalmente em até 28 anos. Foram realizadas correções ocupação do solo de forma inadequada.
geométricas e radiométricas, classificação não O estudo visa contribuir com a manutenção e conservação
supervisionada por pixel com o classificador KMedias. das áreas protegidas da Esec Raso da Catarina, consideradas
Durante 28 anos, houve crescimento bastante considerável Áreas de Preservação Permanente. A região escolhida é
nas áreas antropizadas (1.512,87 km 2 ou 18,22 da taxa de considerada de “extrema importância biológica e prioritária
crescimento) e nas áreas gramíneo lenhosa/área em para conservação”, sua vegetação natural predominante é de
processo de antropização (2.233,20 km2 ou 7,51 da taxa de Caatinga com vegetação predominantemente arbustiva
crescimento). As utilizações das imagens de satélite densa, onde existem várias espécies de ocorrência endêmica
possibilitaram a análise da evolução dos cenários da e ameaçadas de extinção [2].
cobertura da terra em distintas épocas. Dentre as espécies endêmicas e de grande importância na
Palavras-chave — Sensoriamento Remoto, detecção de região da Esec Raso da Catarina, destaca-se a arara-azul-de-
mudança, antropização, semiárido. lear (Anodorhynchus leari), espécie alvo de contrabando e
conflito com agricultores pela sua adaptação ao consumo do
ABSTRACT milho e desmatamento da sua única fonte alimentar, o licuri
(Syagus coronata).
The present study was aimed at evaluating land cover Segundo o Plano de Ação Nacional para a Conservação da
changes on Ecological Station Raso da Catarina - Esec arara-azul-de-lear, nenhum levantamento foi realizado até o
region in recent years, based on satellite earth images. Esec momento no sentido de verificar a quantidade de áreas de S.
Raso da Catarina is a Federal Conservation Unit of coronata, como também não há nenhum levantamento dos
complete protection, located among three municipalities: órgãos ambientais sobre o estado de conservação na região
Jeremoabo, Paulo Afonso and Rodelas in the State of Bahia, da Esec Raso da Catarina [3].
Brazil. Digital image processing techniques using Landsat O tema apresentado, de grande abrangência, necessita cada
images 5, 7 and 8 were applied. The images temporal vez mais da adoção de metodologias de identificação,
distance was until 28 years. Geometric and radiometric avaliação e monitoramento ambiental. Nesse caso, pode
corrections and non-supervised pixel classification, using destacar o SR como uma importante ferramenta à análise e
KMedias classifier, were performed. During 28 years, there controle das questões ambientais, sendo usado para medir e
was a huge increase the anthropic areas (1,512.87 km 2 or monitorar importantes características biofísicas e atividades
18,22 of growth rate) and the areas woody grass/areas with humanas na Terra [4].
the anthropization process (2,233.20 km2 or 7,51 of growth Assim, o presente trabalho se justifica pelo fato de fomentar
rate).The uses of the satellite images enable an analysis of a importância ecológica desta área e de buscar entender a
the evolution of land cover scenarios at different times. sua heterogeneidade espacial de tal área para o
Key words — Remote Sensing, change detection, fortalecimento de práticas sustentáveis, ajudando no suporte
anthropisation, semi-arid. e monitoramento para conservação de espécies no bioma

1934
Caatinga (A. leari e S. coronata), fortalecendo o sertanejo Foram aplicadas técnicas de processamento digital de
que encontra nessas espécies uma fonte adicional de renda. imagens, como correção radiométrica nas imagens Landsat
Portanto, neste estudo o objetivo é analisar temporalmente 5, 7 e 8, e geométrica nas imagens Landsat 5 e 7 no sistema
mudanças ocorridas na cobertura da terra na região da Esec Quantum GIS versão 2.14.1.
Raso da Catarina a partir de imagens do satélite Landsat 5, A correção radiométrica foi realizada através da correção
Landsat 7 e Landsat 8. atmosférica denominada DOS1 – Dark Object Subtraction,
presente no complemento SCP - Semi-Automatic
2. MATERIAIS E MÉTODOS Classification Plugin, para as imagens Landsat.
Para a correção geométrica utilizou-se o registro das
2.1. Área de Estudo imagens Landsat 5 e 7, tomando como base a imagem
Abrange parte dos municípios de Jeremoabo, Paulo Afonso Landsat 8 (obtida já georreferenciada).
e Rodelas, no estado da Bahia, Brasil, região do entorno da O processamento digital das imagens foi realizado
Esec Raso da Catarina (retângulo envolvente), conforme utilizando o sistema SPRING - Sistema de Processamento
Figura 1. de Informações Geográficas/INPE versão 5.3. Foram usadas
as bandas e composição R(4) G(3) B(2) para Landsat 5 e 7,
R(5) G(4) B(3) para a Landsat 8.

3.4. Classificação das Imagens

Com as imagens geradas após o pré-processamento e
composição foram realizados diferentes testes de
segmentações utilizando o método Crescimento de Regiões
para com diferentes similaridades e áreas (similaridade-área:
10-10; 15-15; 20-20; 25-25; 30-30; 35-35; 40-40; 45-45 e
50-50), as quais não foram satisfatórias nos processos de
segmentação para os diferentes sensores. Após esta etapa foi
executada a classificação não supervisionada por pixel com
o classificador KMedias, e também foram realizados
diferentes testes para os temas e interações (tema-interação:
3-3; 4-3; 5-3; 6-3; 7-3; 8-3; 9-3; 10-3; 3-4; 4-4; 5-4; 6-4; 7-
Figura 01. Localização da área de estudo (retângulo envolvente), 4; 8-4; 9-4; 10-4; 3-5; 4-5; 5-5; 6-5; 7-5; 8-5; 9-5 e 10-5).
região da Esec Raso da Catarina, Jeremoabo, Paulo Afonso e Dentre os testes realizados o que melhor representou a real
Rodelas. situação foi com o tema 5 e interação 3 para as imagens
Landsat 5, 7 e 8, com base na informação espectral das
O retângulo envolvente é a área utilizada para as análises do bandas selecionadas.
estudo. A área de interesse foi delimitada a partir do buffer O algoritmo KMedias é um classificador não supervisionado
da Estação Ecológica – Esec, que envolve toda a Zona de que utiliza uma abordagem de agrupamento. O espaço de
Amortecimento da Esec, as comunidades (indígenas, atributos da imagem é partido em K grupos. A partir dessa
assentamentos, povoados e cidade do seu entorno). situação inicial, cada pixel da imagem é alocado ao centro
A Esec Raso da Catarina é uma Unidade de Conservação mais próximo segundo a distância euclidiana [5].
Federal de proteção integral, criada pelo Decreto nº As cinco classes consideradas foram: 1. Área antropizada; 2.
89.268/84, como Reserva Ecológica Raso da Catarina, Gramíneo lenhosa/área em processo de antropização; 3.
sendo recategorizada para Esec Raso da Catarina pela Vegetação rasteira/Caatinga arborizada 4. Vegetação
Portaria n° 373/2001. Localizada, nos municípios de densa/perenifólia; e 5. Corpos d’água.
Jeremoabo, Paulo Afonso e Rodelas. Tem extensão de 8.720
Km2 [4], sendo a segunda maior unidade de conservação de 3.5. Validação e ajustes
proteção integral no estado da Bahia [3] (Figura 1). Após realizada a classificações, efetuou-se uma análise de
exatidão, através da matriz de erro onde se calculou os erros
2.2. Aquisição de dados de comissão (pixels de outras classes que foram atribuídos à
Para a realização deste trabalho foram utilizadas imagens classe de referência) e omissão (pixels pertencentes a uma
orbitais Landsat 5, 7 e 8, para um período de 28 anos classe de referência que foram atribuídos a outras classes)
(07/10/21987, 05/10/2001 e 21/11/2015 - linha 67, cenas para cada classe, bem como a precisão do produtor,
216), obtidas por meio de download gratuito do catálogo de acurácias de exatidão global e índice de Kappa [6].
imagens do Instituto Nacional de Pesquisas Espaciais - Após a validação da classificação foi feita uma edição da
INPE e United States Geological Survey - USGS. classificação para eliminar possíveis erros e omissões
visando o refinamento das regiões classificadas,
2.3. Pré-processamento considerando verificação e interpretação visual do resultado.

1935
 3. RESULTADOS certa conservação dos recursos naturais envolvendo
principalmente a vegetação.
A classificação digital não supervisionada permitiu analisar Com a geração dos mapas de cobertura da terra em 2001 e
a distribuição espacial e quantificar a cobertura da terra na 2015 (Figura 2) e na Tabela 2 que apresenta o cálculo das
Esec Raso da Catarina e no seu entorno. Os resultados são áreas de cada classe e a diferença para cada classe em
evidenciados em três mapas que correspondem aos anos de ambos os anos analisados da cobertura, observam-se as
1987, 2001 e 2015, Figura 2 e na Tabela 01, que consistem modificações consideradas nas distribuições das classes na
na quantificação das formas da cobertura da terra, ao longo área de estudo.
do período de tempo avaliado, o qual permitiu a Primeiramente, nota-se uma diminuição considerável da
individualização das classes. cobertura de vegetação da área de estudo. As áreas de
vegetação rasteira/Caatinga arborizada e vegetação
densa/perenifólia diminuíram, em relação com os anos de
1987 e 2001, deram espaço as áreas de classes de gramíneo
lenhosa/área em processo de antropização 2.513,00 Km²
(35,01%) e áreas antropizadas 1.166,40 km² (16,25%). Tal
acréscimo foi causado principalmente, pela expansão das
áreas agrícolas, áreas agrícolas irrigadas e a expansão das
comunidades, fatos esses que alteraram a dinâmica do solo.
No ano de 2015, a cobertura da terra na área de estudo é
caracterizada pelas intensificações das tendências de
mudanças já descritas em 2001.
Dentro de um determinado espaço de tempo, 28 anos, houve
um crescimento bastante considerável nas áreas
antropizadas (1.512,87 km2) e nas áreas gramíneo
Figura 2. Cobertura da terra na Estação Ecológica do Raso da lenhosa/área em processo de antropização (2.233,20 km2)
Catarina e do seu entorno nos anos de 1987, 2001 e 2015 . apresentando as maiores taxas de incrementos de aumento,
18,22 e 7,51 vezes aumentaram, respectivamente.
Tabela 01. Comparação da cobertura da terra de 1987, 2001 e Ainda neste estudo comparativo houve diminuição do
2015. perímetro referente às áreas com cobertura vegetal com
Área (km2) Diferenças (km2) redução de 2.230,39 km2 para vegetação rasteira/Caatinga
Classes
1987- 2001- 1987-
arborizada e 1.284,72 km2 para vegetação densa/perenifólia
1987 2001 2015
2001 2015 2015 diminuíram - 0,55 e – 0,20 vezes, respectivamente.
1. Área Os corpos d’água não apresentam modificações
83,02 1.166,45 1.512,87 1.083,25 346,42 1.429,85
antropizada significantes em seus valores ao longo dos anos.
A partir da análise das matrizes de erro podem-se obter os
2. Gramíneo
lenhosa/área em valores dos erros de omissão e comissão para cada classe
342,80 2.513,00 2.576,00 2.170,20 63,00 2.233,20
processo de
antropização
das imagens, precisão do produtor, acurácias de exatidão
global e índice de Kappa: Para a classificação da imagem
3. Vegetação Landsat 5, o erro de omissão foi observada para a classe 1
rasteira/Caatinga 5.007,84 3.177,85 2.777,45 -1.829,99 -400,40 -2.230,39 com 14%, classe 5 com 10% e classe 2 com 3%, os
arborizada
melhores desempenhos do classificador de precisão do
4. Vegetação produtor foram as classes 3 e 4 tiveram, ambas com 100%;
densa/ 1.604,42 320,00 319,70 -1.284,42 -0,30 -1.284.72 Para comissão o maior erro foi observado para a classe 2
perenifólia
com 22%, os melhores desempenhos do classificador de
5. Corpos d’água 1,75 1,0 1,33 - 0,75 0,33 -0,42 precisão do produtor foram as classes 1 e 5 com 100%;
acurácia de exatidão global de 99,8% e índice de Kappa
Área total das
7.039,83 7.178,30 7.187,35 igual a 0,99. Landsat 7, a omissão observou-se para a classe
classes
4 com 11%, classe 5 com 3%, classes 1 e 3 com 2%, a
classe 2 teve o melhor desempenho do classificador, com
Por meio da análise do mapa da cobertura da terra em 1987, 100% de precisão do produtor; Na comissão o maior erro foi
foi possível verificar que no ano de 1987 a cobertura vegetal observado para a classe 5 com 30%, classe 3 com 15%,
formada pelas classes de vegetação rasteira/Caatinga classe 4 com 2% e classe 2 com 1%, a classe 1 tive melhor
arborizada 5.007,84 km² (71,14%) e vegetação desempenho do classificador, com 100% de precisão do
densa/perenifólia 1.604,42 km² (22,80%), eram as maiores produtor; exatidão global de 94,6% e índice de Kappa 0,93.
contribuições na área de estudo (somatório de 93,94%), Landsat 8, o maior erro de omissão foi observado para a
mostrando assim, que neste período ainda se verificava uma classe 1 com 4%, classe 5 com 3% e classe 4 com apenas

1936
 1%, a melhor classificação foi encontrada para a classe 2 e período de tempo avaliado, o qual permitiu a
4, com 100% de precisão do produtor; para a comissão o individualização das classes em uma análise temporal.
maior erro foi observado apenas na classe 4 com 0,5%, as As utilizações de imagens de satélite possibilitaram a análise
classes 1, 2, 3 e 5 tiveram melhores desempenhos do da evolução dos cenários da cobertura da terra em distintas
classificador, com 100% de precisão do produtor; exatidão épocas, mostrando a real situação do processo de
global 99,8% e índice de Kappa 0,99. antropização na região. Fica evidente que a principal
ameaça para a espécie da arara-azul-de-lear (Anodorhynchus
4. DISCUSSÃO leari) é a redução das áreas de licurizeiros (Syagrus
coronata), pois estas áreas estão inseridas em local a onde
O grande crescimento nos últimos anos da agricultura, e se verifica a presença humana em diversas comunidades do
áreas de pastagem para a criação (de bovinos, caprinos e entorto na Esec Raso da Catarina, consequentemente o
ovinos) vem acarretando essa degradação ambiental na desmatamento é a origem de uma série de conflitos
região da Esec Raso da Catarina [3]. Esses dados também ambientais gerados na região, inseridos nas áreas
foram observados na classificação da cobertura da terra em antropizada ou gramíneo lenhosa/área em processo de
estudo realizado por Oliveira et al.[7], com imagens Landsat antropização.
8 no mesmo período para o ano de 2014, na qual a região,
os três municípios (Jeremoabo, Paulo Afonso e Rodelas), 6. REFERÊNCIAS
apresentaram maior quantidade de pastagem (4552 km²) que
a de vegetação (4050,8 km²). [1] IPCC. Painel Intergovernamental sobre Mudança do
Em análise e comparação com os dados apresentados pelo Clima., Mudança do clima 2007: Impactos, adaptação e
Instituto do Meio Ambiente e Recursos Hídricos [8], a partir vulnerabilidade. Disponível em: <http://
de arquivo formato shapefile da vegetação do estado da www.mct.gov.br/index.php/content/view/50400.html>.
Bahia, pode-se observar que as classes referidas na Acesso em: 19. jun. 2018.
classificação da cobertura da terra tiveram resultados bem [2] SÁ, I.B.; RICHÉ, G.R.; FOTIUS, G.A. Fatores
semelhantes, uma vez que as áreas da classe vegetação abióticos: áreas e ações prioritárias para conservação da
rasteira/Caatinga arborizada estão conforme as áreas caatinga. In: SILVA, J.M.C. et al. (Org.). Biodiversidade da
apresentadas. caatinga: e ações prioritárias para conservação. Brasília:
Portanto, o método da classificação não supervisionada por p.37-44, 2003.
pixel com o classificador Kmedias obteve desempenho [3] LUGARINI, C.; BARBOSA, A.E.A.; OLIVEIRA, K.G.
excelente para as três imagens analisadas e observadas. (Orgs.), Plano de Ação Nacional para a Conservação da
Assim, os resultados comprovam o excelente desempenho Arara-azul-de-Lear. 2ª Ed. Série Espécies Ameaçadas.
da classificação a partir da relação das matrizes de erros de Instituto Chico Mendes de Conservação da Biodiversidade,
omissão e comissão, precisão do produtor, acurácia de Brasília, p. 144. 2012.
exatidão global (99,8%; 94,6% e 99,8%) das classes e os [4] JENSEN, J.R., Sensoriamento remoto do ambiente: uma
valores dos índices de Kappa (0,99; 0,93 e 0,99) que foram perspectiva em recursos terrestres. Tradução: EPIPHANIO,
comparados e considerados como excelentes de acordo com J. C. N. (coordenador). São José dos Campos, p.598, 2009.
a tabela de referência de Landis e Koch [6]. [5] MATHER P.M., Computer processing of remotely
sensed images; an introduction. 2nd Ed. John Wiley & Sons,
5. CONCLUSÕES New York, 1999.
[6] Lenney, M.P.; Woodcock, C.E.; Collins, J.B.; Hamdi,
O mapeamento a partir de técnicas de SR permitiu a H., The Status of Agricultural Lands in Egypt: The Use of
visualização da dinâmica temporal da cobertura da terra na Multitemporal NDVI Features Derived from Landsat TM.
região da Esec Raso da Catarina. Com a utilização das Remote Sensing of Environment, v. 56, n. 1, p. 8-20, 1996.
imagens Landsat 5, 7 e 8, através da classificação não [7] OLIVEIRA, U.R.; BRITO, P.L.; ZALOTI, F.A.;
supervisionada por pixel com o classificador KMedias foi HADLICH, G.M. Sensoriamento Remoto na Classificação
possível reconstituir os cenários da cobertura da terra e do Uso e Ocupação da Região da Estação Ecológica Raso da
quantificá-los. Catarina, Bahia, Brasil. In: X Jornada de Educação em
As análises de suas principais classes como área antropizada Sensoriamento Remoto no âmbito do Mercosul. Educação,
e cobertura vegetal (gramíneo lenhosa/área; vegetação Sensoriamento e Remoto e Geotecnologias 2015, Lençois -
densa/perenifólia e vegetação rasteira/caatinga arborizada) BA, 2015.
permitiram uma modelagem cartográfica do cenário atual e [8] INEMA - Instituto do Meio Ambiente e Recursos
temporal, para o intervalo de 28 anos, apresentados em três Hídricos. Mapas Tematicos. (2012). Disponível em:
mapas que correspondem aos anos de 1987, 2001 e 2015, http://www.inema.ba.gov.br/servicos/mapas-
Figura 02 e Tabela 01. Tais dados consistem na tematicos/?dl_page=2>. Acesso em: 29. set. 2018.
quantificação das formas da cobertura da terra, ao longo do

Powered by TCPDF (www.tcpdf.org)

1937