norge0 / arvores_lx_r Goto Github PK

#Instalar e carregar biblioteca para dados espaciais

install.packages("raster")
library(raster)
install.packages("rgdal")
library(rgdal

#Definir pasta de trabalho
setwd("G:/data5")

#Importar shapefiles
concelhos = shapefile("concelhos.shp")
arvores_lisboa= shapefile("Arvoredo.shp")

#Ver os dados
plot(concelhos)
plot(arvores_lisboa)

#Inspecionar campos
summary(concelhos)
summary(arvores_lisboa)

#Selecionar apenas choupos 'Populus alba' e 'Populus nigra’
#Identificar valores únicos do campo espécie (ESPECIE_VA que é a coluna4 [,4])
unique(arvores_lisboa$ESPECIE_VA)

#São as linhas 5 e 30 na tabela
[5] "Populus nigra"
[30] "Populus alba"

#Selecionar apenas('Populus alba' e 'Populus nigra’)
choupos_lisboa = subset(arvores_lisboa, ESPECIE_VA == "Populus alba" | ESPECIE_VA == "Populus nigra")

#Para visualizar o resultado
plot(choupos_lisboa)

...