341 lines
25 KiB
Plaintext
341 lines
25 KiB
Plaintext
---
|
|
title: "Arquivos com informações do assentamento"
|
|
author: "Adenor Vicente Wendling"
|
|
date: "24/02/2021"
|
|
output:
|
|
word_document:
|
|
toc: yes
|
|
html_document:
|
|
toc: yes
|
|
toc_float: yes
|
|
theme: united
|
|
highlight: zenburn
|
|
df_print: paged
|
|
---
|
|
|
|
```{r setup, include=FALSE}
|
|
knitr::opts_chunk$set(echo = TRUE)
|
|
```
|
|
|
|
# Uso de Bancos de dados para confecção de mapas de interesse.
|
|
|
|
Existem inúmeros sites que disponibilizam arquivos em formatos vetorial ou Raster gratuitamente.
|
|
|
|
O site <https://forest-gis.com/download-gis-base-de-dados/> paresenta uma coletâne dos principais mapas disponibilizados, com seus respectidos links. Pela abrangência de detalhamento dos dados, não há necessidade de repetí-los aqui.
|
|
|
|
## Divisão política e administrativa
|
|
|
|
Para usar os vetores da divisão política administrativa, onde se encontram os mapas dos estados e munícípios, recomendo utilizar os mapas do IBGE. Sempre inicio os projetos visualizando o mapa do brasil e a divisão político administrativo, pois acho que facilita a localização.
|
|
|
|
|
|
## Delimitação da área dos assentamentos.
|
|
Para a delimitação dos assentamentos que serão trabalhados no nosso exemplo, podemos utilizar os mapas do INCRA, no link <https://certificacao.incra.gov.br/csv_shp/export_shp.py>. O site oferece a opção de selecionar o estado e tipo de dados.
|
|
|
|
|
|
## CAR - Cadastro ambiental rural
|
|
Para os dados que atendem as questões ambientais legais, o governo disponibiliza o CAR, cadastro ambiental rural.
|
|
Os dados estão disponíveis no endereço <https://www.car.gov.br/publico/municipios/downloads> e estão divididos por estados e municípios.
|
|
|
|
## Características do terreno
|
|
O acesso aos dados - imagens rasters - são através do link abaixo: <http://www.webmapit.com.br/inpe/topodata/>. Neste banco de dados podemos obter dados de: altitude, declividade, orientação, relevo sombreado, com resolução de 30 x 30 m.
|
|
|
|
## Curvas de nível - elevação
|
|
O exerxito brasileiro oferece inúmeros dados através de mapas, sendo o mais importante para elaborar o estudo da área dos assentamentos é o mapa de elvação (curvas de nível do terreno), no link <https://bdgex.eb.mil.br/bdgexapp>. É preciso preencher um cadastro, mas os dados são públlicos e após o cadastro, facilmente obtidos no site.
|
|
Outra forma de opter as curvas de nível, é através do arquivo raster de **altitude**. Neste caso, segue-se os seguintes passos:
|
|
Selecione **Raster** no menu, em seguida **Extrair** e depois **contornos**
|
|
Na janela aberta, complete ou selecione os dados necessários, principalmente *equidistância entre contornos*.
|
|
No caso dos assentamentos, as curvas de nível foram definidas com uma distância de 10 m verticais, conforme figura ()
|
|
|
|
|
|
# Inserir as camadas desejadas no QGIS
|
|
Neste capítulo será apresentando um passo a passo para inserir e visualizar as camadas citadas acimea.
|
|
|
|
PAra todos as camadas, vale:
|
|
>>Depois de baixado, os arquivos devem ser descompactados, e armazenados em local específico, para seu uso.
|
|
|
|
## Divisão política e administrativa
|
|
Quando se trabalha com bases georeferencias, a visualizaçao da localização a nível municipal facilita as atividades. Por isso, uma das primeiras camadas a inserir no projeto, é a divisão política e administrativa do estado, em níveil municipal.
|
|
|
|
A base de dados que pode ser utilizadas é a do IBGE, e pode ser localizada no link <https://www.ibge.gov.br/geociencias/organizacao-do-territorio/malhas-territoriais/15774-malhas.html?=&t=downloads>. Neste link podemos selecionar o estado, e o tipo de nível de divisão que nos interessa e baixar o arquivos em formato .zip.
|
|
|
|
Depois de inserido a camada, podemos ter a seguinte visão na tela figura 4).
|
|
|
|
|
|
Para recortar o município de interesse podemos seguir os seguintes passos:
|
|
1. Clique com o botão direito sobre a camada, na aba *camadas*. Em seguida, clique em *abrir tabela de atributos*. Na janela que se abre, selecione a linha que contém os dados do município de interesse. No exemplo, foi selecionado o município de Palmas.
|
|
2. Feche a janela dos atributos.
|
|
3. Na aba *caixa de ferramentas de processamento*, busque por **recortar vetor**. Na lista que aparece, selecione com dois cliques rápidos, *recortar*. Nesta janela, selecione a camada com o mapa do seu estado (Ex.41MUS00G), marque a opção **Apenas feições selecionadas**, tanto na camada de entrada quanto na camada de sobreposição. 4. Finalmente clique em **Executar** e depois em **close**. Se tudo correu bem, foi inserida uma nova camada na aba **camadas**. Desmarque a camada de seu estado, e aparecerá apenas o município recortado.
|
|
5. Para salvar esta camada clique com o botão direito do mouse sobre a nova camada (ex. Recortada) e em seguida em *Tornar permanente*. Na nova janela, selecione o nome e caminho onde armazenar a camada (figura 5)
|
|
|
|
|
|
|
|
Agora tanto a camada com todos os municípios, quanto a camada com o mapa do município de Palmas podem ser visualizados na janela principal.
|
|
|
|
## Delimitação da área dos assentamentos.
|
|
O INCRA fornece um conjunto de arquivos shapefile (vetoriais) com a delimitação dos assentamentos do Brasil, com diversas informações úteis <https://certificacao.incra.gov.br/csv_shp/export_shp.py>
|
|
Selecione os dados e o estado do seu interesse (Ex. Projetos Assentamentos todos; Paraná). Clinque em enviar, e quando finalizado oprocessamento dos mpas, clique sobre o arquivo para donwload. Depois de baixados, e nécessário extrair os arquivos e salva-los em local de sua preferência.
|
|
|
|
Insira a camada, clicando em **camada**, depois ** adicionar camada** e depois em **adicionar camada vetorial**. Selecione o arquivo "Assentamentos.....shp", e clique em **adicionar**. No exemplo, será adicionada a camada com o mapa de todos os assentamentos do Paraná (Figura 6).
|
|
|
|
|
|
Para criar uma camada com os assentamentos de interesse, siga os passos 1 a 5 do capítulo anterior. No exemplo, selecionamos os assentamentos Sâo Lourenço, Margem do Iratim e Paraíso do Sul (Figura 7)
|
|
|
|
|
|
|
|
|
|
## CAR - Cadastro ambiental rural
|
|
|
|
Os dados disponibilizados estão divididos em vários grupos, conforme se visualiza na imagem 1.
|
|
.
|
|
|
|
### Inserir dados no projeto.
|
|
Para inserir os dados dos arquivos no projeto, deverão ser seguidos os seguintes passos
|
|
|
|
Primeiro, selecione `Camada` no menu, e depois, `AdicionarCamada` e finalmente, `adicionar camada vetorial` (Figura 2).
|
|
|
|
.
|
|
|
|
Na janela que se abre (Figura 3), faze-se necessário escolher o caminho em `base(s) de vetores`, e clicar em `Adicionar` seguido do clique em `close`. A camada será mostrada na janela do Qgis, conforme figura 4.
|
|
|
|
.
|
|
|
|
.
|
|
|
|
### Extrair área de interesse
|
|
Ao inserir a camada com algum vetor da base de dados do `CAR`, ela será de todo o município. Quando a área de interesse não for de todo o município, poderão ser extraídos apenas os vetores da área de interesse. Para isso é necessário ter o vetor da área de interesse já defenido, e inserido nas camadas da janela de visualizaao do Qgis.
|
|
No nosso caso, a área de interesse é a área toral do assentamento São Lourenço, cujo vetor tem o nome de `Assentamento_SL`.
|
|
|
|
### Alterar SRC
|
|
Para a correta funcionalidade da apresentação dos mapas, é fundamental que todos esteja com a mesma referencia geográfica, representada no `SRC`. Caso o SRC do mapa inserido seja diferente do mapa já existente, ou do SRC do projeto, é necessário alterar o mesmo.
|
|
Para alterar o SRC da camada raster, clique em **Raster** no menu principal, depois em **Projeções** e em seguida em **Reprojetar coordenadas**, conforme mostrado na Figura 3.
|
|
|
|
|
|
|
|
Na janela da reprojeção, selecione a camada cuja SRC você queira alterar, e os demais campos disponíveis: SRC de origem, SRC de destino, e *sem compactação* conforme mostra a Figura 4.
|
|
|
|
|
|
|
|
|
|
#### Salvar camada com SRC reprojetado
|
|
Após essa Reprojeção, aparecerá uma nova camada no Qgis. é importante que você salve esta nova camada em um local conhecido.
|
|
|
|
|
|
|
|
#### Remover arquivos desnecessários
|
|
Salvo o arquivo, você pode remover as camadas rasters anteriores para evitar confusão, conforme figura 6.
|
|
|
|
|
|
|
|
|
|
A extração dos vetores da área de interesse pode ser feita com o comando `extrair vetor com máscara de entrada`, selecionada diretamente na caixa de ferramentas, ou seguindo os passos a seguir.
|
|
|
|
Através dos dados do CAR, foi elaborado o mapa com indicação das àreas de preservação permanente (APP), área de reserva legal, e nascentes (figura XX)
|
|
|
|
|
|
|
|
## Características do terreno
|
|
|
|
### Os dados do **Topodata**
|
|
|
|
A seguir tem uma pequena transcrição sobre os dados do topodata. "O projeto Topodata oferece o Modelo Digital de Elevação (MDE) e suas derivações locais básicas em cobertura nacional, ora elaborados a partir dos dados SRTM disponibilizados pelo USGS na rede mundial de computadores.
|
|
|
|
Desde que o Topodata foi lançado pela primeira vez, em agosto de 2008, o processamento dos dados foi sucessivamente inspecionado e revisado, com vistas a aprimoramentos e correções. Os dados inicialmente disponibilizados seguiram fielmente as opções e especificações constantes no “Guia de utilização” associado ao Topodata. Porém, problemas na articulação entre folhas e a demanda por mais formatos levaram a um novo tratamento dos dados desde sua preparação, e detalhes do processamento de derivação geomorfométrica foram oportunamente melhorados, e estes novos produtos estiveram disponíveis desde o dia 6 de maio de 2009.
|
|
|
|
Para possibilitar uma futura expansão do Topodata, foi feita uma nova revisão dos produtos e processos, que culminou numa metodologia passível de aplicação onde quer que existam dados SRTM. Os dados atualmente disponíveis, desde novembro de 2011, foram elaborados em fiel correspondência a estes procedimentos."Fonte <http://www.dsr.inpe.br/topodata/index.php>
|
|
|
|
O acesso aos dados - imagens rasters - são através do link abaixo: <http://www.webmapit.com.br/inpe/topodata/>. Neste banco de dados podemos obter dados de: altitude, declividade, orientação, relevo sombreado, com resolução de 30 x 30 m.
|
|
|
|
## Formato e descrição doss arquivos
|
|
Os arquivos das imagens de declividade, altitude,orientação, relevo e curvaturas, são em formato disponibilizados em .gif.
|
|
A nomemnclatura segue um padrão, sendo:
|
|
Os primeiros dois números representam a latitude da imagem;
|
|
A letra S significa latude sul, e a letra N significa latitude Norte;
|
|
Os três últimos números são da longitude.
|
|
As letras finais indicam a finalidade da imagem, sendo:
|
|
`**ZN** para altitude`
|
|
`**SN** para Declividade`
|
|
`**ON** para Orientação`
|
|
`**OC** para Orientação octante`
|
|
`**RS** para Relevo sombreado`
|
|
`**VN** para Curvatura Vertical de 3 classes`
|
|
`**H3** representa imagens de curvatura Horizontal de 3 classes.`
|
|
|
|
## Passos para Baixar o arquivo de interesse
|
|
|
|
Os arquivos são acessados com facilidade, sem necessidade de conta ou cadastro. Clique na área de interesse, e a seguir, no arquivo de interesse. **feito**
|
|
|
|
Descompacte o arquivo, e o salve em local conhecido. É importante ter uma boa organização nos arquivos, com pastas separadas para arquivos originais, em edição e finalizados.
|
|
|
|
## Agora no qgis
|
|
|
|
### Inserir o arquivo
|
|
|
|
Através do caminho apontado na Figura 1.
|
|
|
|
|
|
|
|
|
|
```{r}
|
|
|
|
```
|
|
|
|
Na tela que se abre, busque o arquivo baixado. Clique em **Adicionar** e em seguida em **close**.
|
|
A imagem estará adicionada como uma nova camada, e provavelmente estará visivel (como primeira camada). Caso não seteja visível, mova a camada para o topo da lista das camadas.
|
|
|
|
|
|
|
|
|
|
|
|
Os arquvos raster do topodata estão com o SRC EPSG 4674, que é o formato SIRGAS 2000, oficial para a Amárica Latina. Esta informação pode ser obtida, clicando em **propriedades** e depois em **informações**. Quando nosso projeto está sendo elaborado em outro SRC, todos as camadas devem ser reprojetadas para este SRC. No nosso caso, estamos trabalhando com o SRC 31982.
|
|
|
|
## Alterar SRC
|
|
|
|
Para alterar o SRC da camada raster, clique em **Raster** no menu principal, depois em **Projeções** e em seguida em **Reprojetar coordenadas**, conforme mostrado na Figura 3.
|
|
|
|
|
|
|
|
|
|
Na janela da reprojeção, selecione a camada cuja SRC você queira alterar, e os demais campos disponíveis: SRC de origem, SRC de destino, e *sem compactação* conforme mostra a Figura 4.
|
|
|
|
|
|
|
|
|
|
|
|
### Salvar camada com SRC reprojetado
|
|
Após essa Reprojeção, aparecerá uma nova camada no Qgis. é importante que você salve esta nova camada em um local conhecido.
|
|
|
|
|
|
|
|
### Remover arquivos desnecessários
|
|
Salvo o arquivo, você pode remover as camadas rasters anteriores para evitar confusão, conforme figura 6.
|
|
|
|
|
|
|
|
|
|
Agora podemos iniciar os recorte da camada, já que possuem a mesma SRC.
|
|
|
|
|
|
### Recortar a área do vetor
|
|
|
|
A estas alturas, você já deve ter delimitado sua área de interesse, conforme já mostramos aqui neste repositório [Veja o documento - Como criar shapefile](CriarShapefilePRV.Rmd).
|
|
Para facilitar os trabalhos, e não sobrecarregar o computador, é recomendável recortar o raster de acordo com a área de interesse, que geralmente é muito menor do que todo o arquivo raster .
|
|
|
|
Na janela para definir os parâmetros do recorte, fique atendo para selecionar os parâmetros apontados na figura 7.
|
|
|
|
|
|
|
|
Executado o comando do recorte, aparecerá mais uma camada na tela do QGIS, agora delimitada de acordo com os parâmetros utilizados. A camada recém criada, provavelmente apresentará pontos com altitute igual a zero (em preto). Esses pontos serão removidos durante o processo de salvamento, seguindo os passos conforme figura 8. Veja, que na aba de salvar o arquivo, selecionamos o atributo **nenhum valor de dados** e inserimos os valores a serem desconsiderados na camada salva. No nosso caso, é a altitude zero que deverá se desconsiderada.
|
|
|
|
|
|
|
|
A camada terá aspectos identicos aos da figura 9.
|
|
|
|
|
|
|
|
### Configurar a camada para visualização.
|
|
A configuração consiste em atribuir cores para as diferentes camadas, de acordo com as minhas necessidades. Existe um "padrão" de paleta de corres, elaboradas pelo INPE, que podem (e devem) ser utilizados para manter o padrão e o profisionalismo do projeto.
|
|
Abaixo cópia de trechos do site <http://blog.webmapit.com.br/2013/02/topodata-paletas-qgis-para-altitude.html> para obtenção e aplicação de palheta de cores para altitudes.
|
|
|
|
`Para comemorar os 18 meses em operação do Mapa Índice TOPODATA é com satisfação que disponibilizamos um conjunto de arquivos QML para serem usados no software QGIS com dados do projeto TOPODATA.`
|
|
|
|
`Estes arquivos definem rampas de cores para o tema Altitude e foram elaborados a partir da adaptação dos valores de elevação e regras de cores contidos nos esquemas de cores SRTM, Terrain, Atlas Shader e ETOPO2 provenientes do software GRASS ( função r.colors ).`
|
|
|
|
1. Faça o download e extração dos arquivos QML <http://www.webmapit.com.br/downloads/qgis4TopodataBrasil.zip>
|
|
2. Na lista de camadas(layers), clique com o botão direito no layer correspondente ao .tif e selecione a opção "Propriedades".
|
|
3. Depois clique no botão "Carregar estilo" e selecione o arquivo .QML correspondente ao estilo que você deseja aplicar e depois clique no botão OK. Veja a Figura 10.
|
|
|
|
|
|
|
|
|
|
Depois de carregado o estilo adequado - no caso de altitude, o estilo **qgis4TopodataBrasil_Etopo2.qml** - é necessário fazer algumas adequações para o local, como por exemplo: Definir o modo de distribuição das cores - *o modo intervalos iguais* parece mais adequados; e o número de classes, conforme demostrado na Figura 11.
|
|
|
|
|
|
|
|
### A imagem Final
|
|
Após esses passos, você tera finalizado a importação e a configuração da camada de altitude do projeto. Para visualizar individualmente a altitude de cada ponto, basta clicar no ícone **Identificar feições**, e depois no ponto de deseja verificar. As definições deste ponto aparecerão no lado direito, ao alto na janela, conforme demonstrado na figura 12.
|
|
|
|
|
|
|
|
## Declividade
|
|
Os arquivos para a declividade do terreno são acrescidos das letras SN, e disponibilizado pelo site topodata.
|
|
A palheta de cores utilizada é foi do catálogo cpt_city, do tema "topografy", palleta "wiki-knutux"
|
|
|
|
|
|
## Altitude- elevação
|
|
Os arquivos para a declividade do terreno são acrescidos das letras ZN, e disponibilizado pelo site topodata.
|
|
A palheta de cores utilizada é foi do catálogo cpt_city, do tema "topografy", palleta "elevation"
|
|
|
|
# Fertilidade do solo
|
|
## V%
|
|
|
|
|
|
## pH
|
|
|
|
|
|
## Matéria Orgânica
|
|
|
|
|
|
|
|
|
|
# Impressão dos mapas
|
|
No processo de elaboração dos mapas do Assentamento São Lourenço, foram inseridos os rasters de altitude e orientação do terreno, conforme roteiro apresentado acima.
|
|
A inserção desses dois rasters gera dois mapas distintos, já que não é possível mesclar os dois mapas numa só visialização (ou pelo menos não tem funcionalidades).
|
|
|
|
### Adicionando itens ao layout
|
|
A impressão dos mapas, em formato vários formatos, pode ser feito através do `layout de impressão`, que pode ser acessado através do menu `Projeto`.
|
|
> O primeiro passa é atribuir um nome ao nome do layout (Figura 13). Lembre-se que para o mesmo projeto podem haver vários layouts.
|
|
|
|
|
|
|
|
>>** Não nos preocupamos neste momento com as configurações de cada adição. Isso será tratado no próximo capítulo.**
|
|
|
|
> O segundo passo, é desenhar (demarcar) as mimensões do mapa dentro da "folha", clicando no botão `Adicionar mapa`, conforme marcado com o ńumero 1 na figura 14.
|
|
|
|
|
|
|
|
> O terceiro passo, é o ajuste do mapa dentro da área demarcada. Para ajustar o zoom ou deslocar a área interna do mapa, clica-se no botão `mover o conteúdo do item`. Para ajustar o item (a demarcação), clica-se no botão `selecionar/mover item`, conforme demonstrado na figura 6.
|
|
|
|
> O quarto passo é a inserção de Barra de escala e seta norte, se for o caso. Para inserir a barra de escala, clica-se no botão `Inserir Barra Escala`, e para inserir a seta norte, clica-se no botão `Inserir Seta Norte`, conforme figura 14.
|
|
> O quinto passo é a inserção do Título do mapa, que será feito clicando no botão `Adicionar Rótulo`.
|
|
> O sexto passo será inclusão da Legenda, clicando no botão `Adicionar Legenda`.
|
|
|
|
### Configuração dos itens adicionados
|
|
Depois de adicionados os itens, deve-se dedicar um tempo para sua edição, conforme será visto a seguir.
|
|
Basicamente, as edições são reallizadas nas ferramentas dispostas no lado direito da tela do layout (Figura 15). Na parte superior vê-se a janela dos itens adicionados, e do meio para a parte inferior, a parte dos modelos e das propriedades de cada item. Ao selecionar um item, as propriedades do mesmo aprecem nesta parte, e é aí que se fazem as configurações.
|
|
Na aba do item, pode-se selecionar os que aparecerão na impressão, e bloquear a edição do item. Para bloquear a edição do conteúdo do interior do item, é preciso acessar a aba abaixo, de propriedades do item.
|
|
|
|
|
|
|
|
#### Configurações do Mapa
|
|
O primeiro item adicionado, também poderá ser o primeiro a ser configurado. Ao clicar no item `mapa 1` e selecionar a aba `propriedades do item`, são mostradas as configurações possíveis, como: Escala, SRC, rotação, travar camadas, travar estilo, grades, enquadramento, posição e tamanho, moldura, etc.
|
|
Neste exemplo, iremos incluir uma moldura, e travar as camadas. **Depois de travar as camadas, posso incluir outro mapa, com outro tema, sem prejuízos ao mapa atual**.
|
|
|
|
#### Configurações da legenda
|
|
Ao inclui a legenda, aparece a legenda de todas as camadas existentes no arquivo do projeto, o que geralmente não é o desejável.
|
|
Pode-se atribuir um Título ao mapa, com as configurações necessárias do título do mapa. No exemplo foi adicionado o título *Legenda*.
|
|
Em seguida, é necessário selecionar a qual mapa a legenda se relaciona (neste caso, ao mapa 1).
|
|
Em relação ao aranjo, selecionamos *simbolos à esquerda* (Figura 16).
|
|
|
|
|
|
|
|
Em relação aos itens da legenda, temos a modificação mais importante e substancial a ser feita. Aqui vamos selecionar e renomear (se necessário) aquilo que irá aparecer na impressão final.
|
|
Se não queremos que apareça tudo o que estiver aparecendo até agora na legenda, deslececiona-se o "combo" `Atualização automática` (Figura 17).
|
|
Em seguida, seleciona-se o item que será modificado, ou excluído. No exemplo, serão mantidos apenas os itens `app_Assenta_Car, Hrido_SL e Assentamento_SL`. que serão ainda editados. Para excluir os demais itens, os mesmos foram selecionados, e em seguida foi clicado no botão **-** na parte inferior da janela. Para editar (renomear) os itens que permanecem na legenda, basta clicar no mesmo, e em seguida no botão *editar* na parte inferior da janela.
|
|
|
|
|
|
|
|
#### Configuração da barra de escala
|
|
A barra de escala geralmente também necessita de alterações, pois suas configurações depende do tamanho da área do mapa. O primeiro passo e certificar-se de que está relacionada ao mapa correto, neste caso, com o mapa 1. A escolha do estilo é muito próprio de cada um, e no exemplo foi selecionado o estilo de `linhas tracejadas do meio`.
|
|
Em relação às unidades, podemos aqui trabalhar com metros, ou km. Foi escolhido a opção em Quilômetros, por ser mais adeuqada, mas a escolha depende de cada caso em específico. O mesmo é recomendado para as demais especificações. No exemplo, definiu-se por 2 segmentos à direita, a divisão em 0,5 unidades (0,5 km) (Figura 18).
|
|
|
|
.
|
|
#### Configurações do Título
|
|
A última configiração necessária para este mapa, é a fonte do título. Atribuimos fonte 12, além de fazer o alinhamento necessário para corresponder ao tamanho do mapa.
|
|
|
|
## Impressão do Mapa
|
|
Para imprimir o mapa, temos a opção de imprimir em impressora, ou em formato pdf, SVG ou png. A opção, para inserir o impresso neste script, foi `APP impressão.png` (Figura 19) com resolução de 100 dpi (para não pesar muito).
|
|
|
|
|
|
|
|
|