#script para exercício de dimensão de paisagem

#Indique seu diretório de trabalho (pasta onde você vai organizar os arquivos da disciplina/exercícios)

setwd("~/Downloads/dados_paisagem")
 #função que indica para o R qual é sua pasta de trabalho

#Bibliotecas #### Biblioteca é um conjunto de funções que espande as funções R (ou outros programas)

library(landscapemetrics) # library(): função para carregar no R as bibliotecas que você precisa
                          #landscapemetrics é o principal pacote que vamos usar para o exercício de hoje
library(raster) #biblioteca para ler e manipular arquivos em Raster... vocês vão descobrir o que é isso no futuro
## Carregando pacotes exigidos: sp
library(sf)     #biblioteca para ler e manipular arquivos shapefile
## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
#library(here)   #pacote para facilitar a localização de arquivos
library (rgdal) # adicionei, pois pedir
## Please note that rgdal will be retired by the end of 2023,
## plan transition to sf/stars/terra functions using GDAL and PROJ
## at your earliest convenience.
## 
## rgdal: version: 1.5-28, (SVN revision 1158)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 3.0.4, released 2020/01/28
## Path to GDAL shared files: /usr/share/gdal
## GDAL binary built with GEOS: TRUE 
## Loaded PROJ runtime: Rel. 6.3.1, February 10th, 2020, [PJ_VERSION: 631]
## Path to PROJ shared files: /usr/share/proj
## Linking to sp version:1.4-5
## To mute warnings of possible GDAL/OSR exportToProj4() degradation,
## use options("rgdal_show_exportToProj4_warnings"="none") before loading sp or rgdal.

#dados #### ## importar #### Aqui vamos carregar para o R os dados que utilizaremos no exercício

Aqui removi o “here” que não estava ajudando

raster("dimensao_paisagem.tif") -> paisagem_geral               #raster(): função para carregar arquigo

st_read("pontos_dimensao_paisagem.shp", quiet = TRUE) -> ponto_focal

manipular e visualizar

check_landscape(paisagem_geral)           #Funçao para checar se sua paisagem esta ok para as analises
##   layer       crs units   class n_classes OK
## 1     1 projected     m integer         3  ✓
plot(paisagem_geral)                      #A cara da nossa paisagem :)
plot(ponto_focal, col = "red", add = T)   #Nosso ponto de amostra

#Análises #### ## paisagem quadrada ####

sample_lsm(landscape = paisagem_geral,
           y = ponto_focal,
           size = 3000,                        #tamanho da paisagem em metros. Modifique aqui para tamanho desejado. para mais de um tamanho basta usar c(1000,2000,3000) para 1km, 2km, 3km
           level = 'class',                    #Para qual nível da paisagem calcular a métrica
           metric = 'area') -> buff_quad_3km   #Qual tipo de métrica vamos escolher. Área no caso do exercício

paisagem redonda

sample_lsm(landscape = paisagem_geral,
           y = ponto_focal,
           size = 3000,
           shape = "circle",                   #Aqui indicamos no argumento 'shape' qual o formato queremos.
           level = 'class',                    
           metric = 'area',) -> buff_circ_3km

Perguntas