Tutorial do Projeto de Extensão de Monitoramento Ambiental no
Laboratório de Geoprocessamento de Ciências Ambientais –
LABGEO/CA
Prof. Herondino dos Santos Filho
Passo 1: Salvar Camada raster como... na Extensão da visão do
Mapa no QGIS
Passo 2: Utilizar a banda x para identificar a região de interesse
- Retirar a região de interesse
- Fazendo o sombreamento (Primeira Imagem)
-Colorindo a Segunda imagem
Passo 3: Para a edição do polígono delimitador da ressaca
Opção 1:
Para a delimitação da área de ressaca poderá ser feita a combinação
em transparência da sintética (33% de transparência) + Sombreada
(Estilo ->Renderização da cor: modo de mistura -50 brilho)
Opção 2:
Para essa perspectiva pode ser feita a troca das telas Colorida e a
sintética através da barra de espaços
Opção 3:
Pode ser feita o mesmo procedimento da opção 1 com a sintética
(33% de transparência) e colorida sem transparência
Passo 4: Recorte do polígono da imagem
Passo 5: Instalação do OTB Toolbox (Ferramentas Orfeo)
Opção 1: através do instalador osgeo4w pelo link:
https://www.youtube.com/watch?v=WBCYG6eddcA&list=LLExiPfMzVS
oHMPbT5vBZb4g&index=2 do ONF Intenacional
Opção 2: baixar do site oficial pelo endereço https://www.orfeo-
toolbox.org/download/
Nesta opção deve ser descompactado e colado no c:\OTB
Nas duas situações é necessário entrar no menu Processar em
Opções, como mostra a figura abaixo:
Obs: deve ser feito o direcionamento da pasta das aplicações OTB e
da pasta das ferramentas da linha de comando OTB.
Opção 1:
Pelo instalador o caminho da pasta de aplicações OTB é:
C:\OSGeo4W64\apps\orfeotoolbox\applications
E da pasta de ferramentas da linhas de comando OTB é:
C:\OSGeo4W64\bin
Opção 2:
Para opção em que é baixado pelo site oficial e instalado no c:\OTB
O caminho da pasta de aplicações OTB é:???
E da pasta de ferramentas da linhas de comando OTB é:???
Passo 6: Classificação pelo OTB Toolbox Orfeo
- Criar uma camada Shapefile tipo: polígono
-Prováveis classes encontradas nas amostras para treinamento
- Estatística da imagem a ser classificada
- Treinamento
-Adicionar também o Stats.xml
- Saída: a Matriz de Confusão e o Modelo
- Classificação da Imagem
7. Poligonização
Das áreas úmidas seguindo a metodologia:
- Curvas de níveis a partir do MDT da Base;
1. Encontrar as curvas de níveis envolvidas (3 ou 4);
2. Verificar a curva de nível predominante e realizar a sua
edição;
- Sombreamento e Falsa-Cor (utilizando a cor de altimetria-
elevation);
8. Densificar Geometria
9. Mesclar Geometrias (Feições)
Fechar seleção com o botão direito!
Ir em Editar e selecionar Mesclar feições selecionadas.
10. Recorte do Polígono
Opção 1:--config GDALWARP_IGNORE_BAD_CUTLINE YES
Opção 2:
from osgeo import ogr, gdal
import subprocess
# Define I/O files
InputImage = "Input.tif"
OutImage = "Out.tif"
Shapefile = "SomeName.shp"`
# Open the shapefile in read-only mode and retrieve the spatial reference
Driver = ogr.GetDriverByName("ESRI Shapefile")
VectorDataset = Driver.Open(Shapefile, 0)
layer = VectorDataset.GetLayer()
feature = layer.GetNextFeature()
geom = feature.GetGeometryRef()
spatialRef = geom.GetSpatialReference()`
# Get bounding area of the layer
minX, maxX, minY, maxY = geom.GetEnvelope()`
# Clip the InputImage to the bounding area of the input feature (using Nearest Neighbour interpolation).
OutTile = gdal.Warp(OutImage, InputImage, format="GTiff", outputBounds=[minX, minY, maxX,
maxY], xRes=30, yRes=30, dstSRS=spatialRef, resampleAlg=gdal.GRA_NearestNeighbour)`
# Build overviews for the OutImage
subprocess.call("gdaladdo --config COMPRESS_OVERVIEW DEFLATE "+OutImage+" 2 4 8 16 32
64", shell=True)
print("Done.")
11. Unido polígonos de camadas distintas
12. Composição colorida
Composição
Bandas 1(gray);2;3
Bandas 4;3;2 -> Falsa Cor 3D
Bandas 3;2;1 -> Real
Bandas 3;4;2 -> Falsa cor
13. Realce
Carregar valores máximo e mínimo
Top Related