require(raster)
## Carregando pacotes exigidos: raster
## Warning: package 'raster' was built under R version 4.3.3
## Carregando pacotes exigidos: sp
## Warning: package 'sp' was built under R version 4.3.3
require(landscapemetrics)
## Carregando pacotes exigidos: landscapemetrics
## Warning: package 'landscapemetrics' was built under R version 4.3.3
require(magrittr)
## Carregando pacotes exigidos: magrittr
## Warning: package 'magrittr' was built under R version 4.3.3
## 
## Attaching package: 'magrittr'
## The following object is masked from 'package:raster':
## 
##     extract
library("raster")
library ("landscapemetrics")
library("magrittr")
setwd("C:/Users/813749/Downloads/dados_07_01/dados_07_01")
getwd()
## [1] "C:/Users/813749/Downloads/dados_07_01/dados_07_01"
r<-raster("C:/Users/813749/Downloads/dados_07_01/dados_07_01/AmostraBuri_17.tif")
plot(r)

check_landscape(r, verbose = TRUE)
## Warning: Caution: Coordinate reference system not metric - Units of results
## based on cellsizes and/or distances may be incorrect.
##   layer        crs   units   class n_classes OK
## 1     1 geographic degrees integer         2  ✖
nova.proj<-"+proj=utm +zone=22 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
r.m<-projectRaster(r, res=30, crs = nova.proj, method="ngb")
par(mfrow=c(1,2))
plot(r)
plot(r.m)

check_landscape(r.m, verbose = TRUE)
##   layer       crs units   class n_classes OK
## 1     1 projected     m integer         2  ✔
?landscapemetrics
## starting httpd help server ... done

P.01: Descreva brevemente 2 métricas de cada nível? R: CA (class level) CA é uma ‘métrica de área e borda’ e uma medida de composição. A área total (classe) soma a área de todas as manchas pertencentes à classe i. Mostra se a paisagem é, por ex. dominada por uma classe ou se todas as classes estiverem igualmente presentes. A CA é uma medida absoluta, dificultando comparações entre paisagens com diferentes áreas totais.

AREA_MN (class level) AREA_MN é uma ‘métrica de área e borda’. A métrica resume cada classe como a média de todas as áreas de manchas pertencentes à classe i. A métrica é uma forma simples de descrever a composição da paisagem. Especialmente junto com a área total da classe (lsm_c_ca), também pode dar uma ideia da estrutura do patch (por exemplo, muitos patches pequenos versus poucos patches grandes).

DCORE_MN (landscape level) DCORE_MN é uma ‘métrica de área central’. Ele resume a paisagem como a média de todas as manchas da paisagem. Uma célula é definida como núcleo se a célula não tiver nenhum vizinho com valor diferente dela (caso da torre). O NCORE conta as áreas centrais disjuntas, sendo que uma área central é uma ‘mancha dentro da mancha’ contendo apenas células centrais.

ENN_MN (landscape level) ENN_CV é uma ‘métrica de agregação’. Ele resume a paisagem como a média de todas as manchas da paisagem. ENN mede a distância até o patch vizinho mais próximo da mesma classe i. A distância é medida de ponta a ponta. O alcance é limitado pela resolução da célula no limite inferior e pela extensão da paisagem no limite superior.

ENN (patch level) ENN é uma ‘métrica de agregação’. A distância até o patch vizinho mais próximo da mesma classe i. A distância é medida de ponta a ponta. O alcance é limitado pela resolução da célula no limite inferior e pela extensão da paisagem no limite superior. A métrica é uma maneira simples de descrever o isolamento do patch.

PERIM (patch level) PERIM é uma ‘métrica de área e borda’. É igual ao perímetro do patch incluindo também a borda do limite da paisagem. A métrica descreve a área do patch (perímetro maior para patches maiores), mas também a forma do patch (perímetro grande para formas irregulares).

area.total<-lsm_l_ta(r.m, directions=8)

P.02: Qual é a área total da paisagem? 1.0003.86

te<-lsm_l_te(r.m, count_boundary=FALSE)
ed8<-lsm_l_ed(r.m, count_boundary = FALSE, directions = 8)
cls.area<-lsm_c_ca(r.m, directions = 8)
area.media.frag<-lsm_c_area_mn(r.m, directions = 8)
cv_tamanho<-lsm_c_area_cv(r.m, directions = 8)
sd_tamanho<-lsm_c_area_sd(r.m, directions = 8)
frag.sizes.8<-lsm_p_area(r.m, directions = 8)
f.s.8<-frag.sizes.8 %>% dplyr::filter(class == 1)

P.03: Qual é o tamanho do maior remanescente de vegetação nativa. R: 400.77m/ha

frag.sizes.4<-lsm_p_area(r.m, directions = 4)
f.s.4<-frag.sizes.4 %>% dplyr::filter(class > 0)