In darrennorris/gisdata: GIS data for processing and analysing rivers

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

Apresentação

---

Aqui vamos pesquisar desmatamento entre 2000 - 2018 ao longo de 276 quilômetros de rios na Amazônia Brasileira. As pesquisas fazem parte de atividades desenvolvidos no projeto Where is my Turtle : https://myturtlebrazil.wixsite.com/whereismyturtle .

• Objetivo não é de apresentar detalhes sobre os cálculos/métodos estatísticas ou os funções no R.

• Mas, sim, o objetivo é de apresentar um exemplo mostrado os capacidades e opções para desenvolver e integrar pesquisas cientificas no ambiente estatística de R

Com menos de 100 linhas de código vamos: 1) carregar dados; 2) olhar mapas (SIG); 3) gerar tabelas com resumos; 4) visulizar resultados em gráficos; e 5) rodar anlises estatisticas.

Pacotes

Para fazer todo isso em menos de 100 linhas, precisamos os seguintes pacotes, que deve esta instalado antes: plyr, tidyverse, sf, mapview, sjPlot, sjmisc, sjlabelled, interactions.

Portanto, deve Instalar os pacotes necessários antes de começar:

install.packages(c("plyr","tidyverse", "raster","sf", "mapview", 
                   "sjPlot", "sjmisc", "sjlabelled", "interactions"))

Carregar pacotes:

library(plyr)
library(tidyverse)
library(sf)
library(mapview)
library(sjPlot)
library(sjmisc)
library(sjlabelled)
library(interactions)

Agora podemos fazer um gráfico com função "ggplot" (pacote ggplot2), que faz parte do "tidyverse". Mais exemplos no R cookbook : http://www.cookbook-r.com/Graphs/ .

Primeiro gráfico com os dados "mtcars" no pacote ggplot2:

ggplot2::ggplot(mtcars, aes(x = wt, y =  mpg)) + geom_point()

Para entender melhor pode verificar ajudar no R:

?mtcars
?geom_point

Mais exemplos com Scatterplots aqui: http://www.cookbook-r.com/Graphs/Scatterplots_(ggplot2)/

## 1) Carregar arquivos com dados de SIG. wzxhzdk:5 Baixar arquivo com os dados (formato "GPKG", tamanho `r round(fs,1)` MB). Link: [https://github.com/darrennorris/gisdata/blob/master/inst/vector/rivers.GPKG](https://github.com/darrennorris/gisdata/blob/master/inst/vector/rivers.GPKG) . Lembrando-se de salvar o arquivo ("rivers.GPKG") em um local conhecido no seu computador. Agora avisar R sobre onde ficar o arquivo. O código abaixo vai abrir uma nova janela, e você deve buscar e selecionar o arquivo "rivers.GPKG": wzxhzdk:6 Agora vamos olhar o que tem no arquivo. wzxhzdk:7 Existem camadas diferentes com pontos e linhas: wzxhzdk:8 ### 1.1) Carregar dados (vectores) O código abaixo vai carregar os dados e criar 3 objetos "rsm", "rsl" e "fl". Agora temos dados com: pontos cada 5 km ao longo os rios (rsm) ; linha central de rios ("rsl") e pontos cada metro ao longos os rios ("fl"). wzxhzdk:9 wzxhzdk:10

## 2) Mapas ### 2.1) Mapa com pontos cada metro de rio Objeto "fl" com 276086 pontos - muitos pontos cerca de 5 minutos para concluir.... wzxhzdk:11 ### 2.2) Mapa com linha central e pontos de rios em trechos de 5km wzxhzdk:12

## 2.3) Interativo (funcione somente com internet) Mostrando agora com fundo de mapas "base" (OpenStreetMap/ESRI etc) wzxhzdk:13 wzxhzdk:14 Mais exemplos com mapas e dados espaciais no R: [sf e ggplot2](https://www.r-spatial.org/r/2018/10/25/ggplot2-sf.html) : https://www.r-spatial.org/r/2018/10/25/ggplot2-sf.html [Capitulo 8](https://geocompr.robinlovelace.net/adv-map.html) no livro [Geocomputation with R](https://geocompr.robinlovelace.net/) : https://geocompr.robinlovelace.net/adv-map.html **************

## 3) Resumo de dados na tabela de atributos ### 3.1) Desmatamento por rio e por seis zonas Agora vamos fazer alguns resumos com os dados na tabela de atributos. No arquivo "fl" temos dados mostrando valores de [desmatamento](http://earthenginepartners.appspot.com/science-2013-global-forest) , "Global Forest Change" desenvolvido por Hansen e co-autores 2013 (https://science.sciencemag.org/content/342/6160/850) e obtidos de no site https://earthenginepartners.appspot.com/science-2013-global-forest/download_v1.6.html . Os valores calculados são uma proporção de área desmatado nos anos 2000-2018 em raios de 500 metros, 1 quilometro, 5 quilômetros e 10 quilômetros. Os dados em pontos cada metro ao longo de 276 quilômetros de rios (276086 pontos) . Vamos calcular comprimento de rio e valor media de desmatamento em seis zonas diferentes: wzxhzdk:15 ### 3.2) Desmatamento por rio e por 52 subzonas Agora, repetindo o mesmo processo, mas acrescentando resumos para 52 subzonas (cerca de 5 km cada) para um analise mais preciso. wzxhzdk:16 Visualizar tabela, aproveitando funções disponíveis no pacote "sjPlot": wzxhzdk:17 Exportar tabela com resumos ("dfsubzona") em formato "html" para word: wzxhzdk:18 Exportar tabela com resumos ("dfsubzona") em formato para excel: wzxhzdk:19

## 4) Gráficos ### 4.1) mostrando as diferenças em dematamento entre zonas wzxhzdk:20 Agora, através de ajustes no código, vamos fazer um gráfico mais elegante e mais informativo..... wzxhzdk:21

## 4.2) Gráfico comparando desmatamento entre raios de 500m, 1km, 5km e 10km Primeiro preciso reorganizar os dados para facilitar apresentação visual: wzxhzdk:22 Agora o gráfico: wzxhzdk:23

## 5) Analises estatistica, comparando desmatamento entre zonas e raios Veja guias mostrando analises estatisticas no R para entender melhor: [Capítulo 12](https://cdr.ibpad.com.br/modelos.html) no livro “Ciência de Dados com R–Introdução......”: https://cdr.ibpad.com.br/modelos.html Modelo adequado para dados de proporção (valores de 0-1). Aqui um modelo comparando efeito de raio (500 m, 1 km, 5 km, 10 km) e zona sobre proporção de desmatamento: wzxhzdk:24 Agora, mostrando três formas diferentes de apresentar resumo do modelo: Primeiramente com base R: wzxhzdk:25 E com funções disponíveis no pacote "sjPlot": wzxhzdk:26 Gráfico de valores do modelo, no estilo "Forest-plot" wzxhzdk:27 E com funçao "cat_plot" disponivel no pacote "interactions": wzxhzdk:28

## 6) Finalizar. Salvar uma copiar de dados para uso futuro. wzxhzdk:29 ## 7) Leitura Para quem buscar mais detalhes sobre os funções/métodos estatísticas no [R](https://cran.r-project.org/): Livro [Processamento e Análise de Dados](https://www.msperlin.com/padfeR/index.html) : https://www.msperlin.com/padfeR/index.html Livro [Ciência de Dados com R](https://cdr.ibpad.com.br/) : https://cdr.ibpad.com.br/ Além disso, existem fontes variadas de ajudar online: [cursos](https://pt.coursera.org/learn/r-programming) : https://pt.coursera.org/learn/r-programming; [R-br](http://r-br.2285057.n4.nabble.com/) a lista Brasileira oficial de discussão do programa R: http://r-br.2285057.n4.nabble.com/, Capítulos 1 a 3 no livro “Processamento e Análise de Dados.....”: https://www.msperlin.com/padfeR/index.html Capítulos 1 a 2 no livro “Ciência de Dados com R–Introdução......”: https://cdr.ibpad.com.br/cdr-intro.pdf Ingles Instruções de instalação: https://moderndive.netlify.com/1-getting-started.html Com 3 aulas básicas: https://rladiessydney.org/courses/ryouwithme/01-basicbasics-0/ Foruns: [rcomputing](https://www.facebook.com/rcomputing/) : https://www.facebook.com/rcomputing/ ; https://www.r-bloggers.com/ exemplo com PCA: https://www.r-bloggers.com/computing-and-visualizing-pca-in-r/

darrennorris/gisdata documentation built on March 7, 2024, 5:50 p.m.