Tutorial Grass6 UFAM

Geoprocessamento com Software Livre:
GRASS-GIS e QGIS
Carlos Henrique Grohmann

[email protected]
Instituto de Geociencias - USP
Sao Paulo, 12 de outubro de 2012
Este documento n
ao e parte da documentacao oficial do software descrito. Mais informacoes podem ser encontradas nos web sites do GRASS-GIS (http://grass.osgeo.org)e do
QGIS (http://qgis.org).
Os detalhes, dados, resultados, etc. descritos neste documento foram escritos e verificados
pelo autor; porem, erros quanto ao conte
udo sao possveis. O autor nao assume qualquer responsabilidade por falhas e suas consequencias. Indicacoes de possveis erros sao bem-vindas.
Este documento foi escrito em LATEX, e esta disponvel em formato PDF em
http://carlosgrohmann.com/ ou como fonte LATEX por solicitacao ao autor.
c 2012 Carlos Henrique Grohmann.

Copyright
garantida a permiss
E
ao para copiar, distribuir e/ou modificar este documento sob os
termos da Licenca de Documentac
ao Livre GNU (GNU Free Documentation License), Versao 1.2 ou qualquer vers
ao posterior publicada pela Free Software Foundation; sem Secoes
Invariantes, Textos de Capa Frontal, e sem Textos de Quarta Capa. Uma copia da licenca e
includa na sec
ao intitulada GNU Free Documentation License.

SUMARIO
Sum
ario
Sum
ario
Lista de Figuras
iii
Lista de Tabelas
vi
1 Introdu
c
ao
1.1 Fontes de informac
ao . . . . . . . . . . . . .
1.2 Design e estrutura de SIGs . . . . . . . . .
1.3 Dimens
ao de dados em SIG . . . . . . . . .
1.4 Estrutura dos comandos no GRASS . . . .
1.5 Obtendo ajuda com relac
ao aos modulos . .
1.6 Organizac
ao dos projetos no GRASS-GIS .
1.7 O mapset PERMANENT . . . . . . . . . .
1.7.1 Arquivos do mapset PERMANENT
1.8 Organizac
ao dos projetos no QGIS . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
2 Mapas Raster - vis

ao geral
1
1
2
4
5
6
6
7
8
8
9
3 Mapas Vetoriais - vis

ao geral
3.1 Manipulac
ao de atributos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10
12
4 Elementos de cartografia
4.1 Escala . . . . . . . . . . . . . . . . . . . . . .
4.2 Ge
oide, elips
oide e datum . . . . . . . . . . .
4.3 Sistemas de coordenadas . . . . . . . . . . . .
4.4 Coordenadas planas . . . . . . . . . . . . . .
4.5 Coordenadas geogr
aficas (latitude-longitude)
4.6 Projec
oes cartogr
aficas . . . . . . . . . . . . .
4.7 Projec
oes cilndricas . . . . . . . . . . . . . .
4.7.1 Cilndrica conformal (Mercator) . . .
4.7.2 Cilndrica equidistante . . . . . . . . .
4.8 Projec
oes c
onicas . . . . . . . . . . . . . . . .
4.8.1 C
onica conformal . . . . . . . . . . . .
4.8.2 C
onica igual-
area . . . . . . . . . . . .
4.8.3 C
onica equidistante . . . . . . . . . .
4.8.4 Polic
onica . . . . . . . . . . . . . . . .
4.9 Projec
oes azimutais . . . . . . . . . . . . . . .
4.9.1 Azimutal igual-
area . . . . . . . . . .
4.9.2 Azimutal equidistante . . . . . . . . .
4.9.3 Azimutal ortogr
afica . . . . . . . . . .
4.9.4 Azimutal gnom
onica . . . . . . . . . .
4.10 Projec
oes compromisso . . . . . . . . . . . . .
14
14
15
18
18
18
22
24
24
24
25
25
26
26
27
27
27
28
28
29
29
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

SUMARIO
ii
.
.
.
.
.
.
29
30
31
32
32
33
5 Bem-vindo(a) ao GRASS-GIS!
5.1 Criando uma Location e um mapset . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2 Apresentando o GIS Manager e o Map Display . . . . . . . . . . . . . . . . . . . . . .
36
36
43
6 Dados Raster
6.1 Modelos de Elevac
ao SRTM . . . .
6.1.1 Modelos SRTM30 PLUS . .
6.2 Download e importac
ao dos dados
6.3 Raster patching . . . . . . . . . . .
6.4 Trabalhando com regions . . . . .
6.4.1 Relevo sombreado . . . . .
.
.
.
.
.
.
45
45
46
47
51
56
56
7 Dados vetoriais
7.1 Arquivos ESRI shape . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
7.2 Dados vetoriais no QGIS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
63
63
63
8 Importa
c
ao de dados vetoriais no GRASS
8.1 Arquivos ASCII . . . . . . . . . . . . . . . .
8.1.1 Gerac
ao de pontos aleat
orios . . . .
8.1.2 Amostragem da altitude . . . . . . .
8.1.3 Exportac
ao como ASCII . . . . . . .
8.1.4 Importac
ao do arquivo ASCII . . . .
.
.
.
.
.
73
78
78
80
83
83
.
.
.
.
.
85
85
88
89
89
93
.
.
.
.
.
.
.
95
95
96
97
97
97
98
99
4.11
4.12
4.13
4.14
4.15
4.10.1 Projec
ao de Robinson . .
Coordenadas UTM . . . . . . . .
Carta Internacional do Mundo ao
Sistema cartogr
afico nacional . .
Padr
ao de exatid
ao cartogr
afica .
Mapeamento sistem
atico . . . . .
9 An
alise Digital de Terreno
9.1 Hipsometria . . . . . . . . . .
9.2 Relevo sombreado . . . . . .
9.3 Indices de forma . . . . . . .
9.3.1 Declividade e Aspecto
9.3.2 Curvaturas . . . . . .
.
.
.
.
.
.
.
.
.
.
. . . . . . . .
. . . . . . . .
Milionesimo .
. . . . . . . .
. . . . . . . .
. . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
10 Interpola
c
ao de superfcies
10.1 Reamostragem . . . . . . . . . . . . . .
10.2 Inverso da potencia da dist
ancia - IDW
10.3 Splines regularizadas com tens
ao - RST
10.4 Interpolac
ao bilinear e bic
ubica . . . . .
10.5 Outros m
odulos de interpolac
ao . . . . .
10.6 Reamostragem dos dados SRTM . . . .
10.7 Extrac
ao de curvas de nvel . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
LISTA DE FIGURAS
10.8 Interpolac
ao de MDE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10.8.1 Recorte de dados vetoriais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10.8.2 Interpolac
ao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iii
99
99
101
11 Convers
ao entre proje
c
oes no GRASS
105
12 Visualiza
c
ao n-dimensional
107
13 O QGIS como interface gr

afica alternativa
110
14 Produzindo mapas para impress

ao
113
15 Comandos do GRASS
115
16 Refer
encias Bibliogr
aficas
130
17 Defini
c
ao de Software Livre
132
18 GNU Free Documentation License
135
Lista de Figuras
1.1
1.2
1.3
1.4
2.1
4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
4.9
4.10
4.11
4.12
4.13
4.14
4.15
4.16
4.17
Geometia e atributos de dados em SIG . . . . . . . . . . . . . . . . . . . . . . . . . . .

Representac
ao do relevo de uma area por estrutura discreta (curvas de nvel vetoriais)
e contnua (superfcie raster gerada por interpolacao). . . . . . . . . . . . . . . . . . .
Dimens
ao de dados em SIG . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Exemplo da estrutura projetos no GRASS-GIS . . . . . . . . . . . . . . . . . . . . . .
Resoluc
ao espacial de um raster . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Exemplos de escalas gr
aficas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Representac
ao das ondulac
oes do geoide, construda a partir de dados de anomalia
gravimetrica (NASA/JPL/University of Texas Center for Space Research). . . . . . .
Relac
oes entre a superfcie da Terra, o geoide e o elipsoide. . . . . . . . . . . . . . . .
Elementos do elips
oide de referencia. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Datum topocentrico (NAD27) e geocentrico (WGS84). . . . . . . . . . . . . . . . . . .
Sistema de coordenadas planas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Sistema de coordenadas geogr
aficas (latitude-longitude). . . . . . . . . . . . . . . . . .
Longitude, latitude geodetica e latitude geocentrica. . . . . . . . . . . . . . . . . . . .
Variac
ao do comprimento do grau de longitude com a latitude. . . . . . . . . . . . . .
Classificac
ao de projec
oes cartograficas . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
oes Secantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao cilndrica conformal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao cilndrica equidistante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao c
onica conforme de Lambert . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao c
onica igual-
area de Albers . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao c
onica equidistante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao polic
onica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
4
5
7
9
14
15
15
16
17
18
19
19
20
23
23
24
25
25
26
26
27
iv
LISTA DE FIGURAS
4.18
4.19
4.20
4.21
4.22
4.23
4.24
4.25
4.26
4.27
5.1
5.2
5.3
5.4
5.5
5.6
5.7
5.8
5.9
5.10
5.11
5.12
5.13
5.14
5.15
5.16
5.17
6.1
6.2
6.3
6.4
6.5
6.6
6.7
6.8
6.9
6.10
6.11
6.12
6.13
6.14
6.15
6.16
Projec
ao azimutal igual-
area de Lambert . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao polar azimutal equidistante . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao azimutal ortogr
afica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao azimutal gnom
onica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Projec
ao de Robinson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Divis
ao das zonas UTM. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Origem das coordenadas nas zonas UTM. . . . . . . . . . . . . . . . . . . . . . . .
Articulac
ao da Carta Internacional ao Milionesimo em territorio brasileiro. . . . .
Articulac
ao do sistema cartogr
afico nacional (visao condensada). . . . . . . . . . .
Articulac
ao do sistema cartogr
afico nacional (visao expandida). . . . . . . . . . . .
Tela de terminal com a mensagem de boas-vindas ao GRASS-GIS. . . . . . . . . .
Janela inicial do GRASS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Definic
ao do diret
orio GISBASE e do nome da Location. . . . . . . . . . . . . . .
Defina como ser
a o metodo de criacao da Location. . . . . . . . . . . . . . . . . .
Defina a projec
ao e o sistema de coordenadas. . . . . . . . . . . . . . . . . . . . . .
Escolha dos par
ametros da projecao. . . . . . . . . . . . . . . . . . . . . . . . . . .
Escolha do datum WGS-84. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Confira se est
a tudo certo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Quase l
a. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Defina os limites e a resoluc
ao da region padrao. . . . . . . . . . . . . . . . . . . .
Location pronta, vamos ao mapset. . . . . . . . . . . . . . . . . . . . . . . . . . .
Defina o nome de nosso mapset. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Interface gr
afica do GRASS-GIS 6.4 . . . . . . . . . . . . . . . . . . . . . . . . . .
Janela do GIS Manager. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Busca de um m
odulo no gism. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Barra de ferramentas do GIS Manager. . . . . . . . . . . . . . . . . . . . . . . . . .
Barra de ferramentas do Map Display. . . . . . . . . . . . . . . . . . . . . . . . . .
Metodos de reamostragem dos dados SRTM utilizados pela NASA e pelo USGS . .
Modelo de elevac
ao global SRTM30 PLUS . . . . . . . . . . . . . . . . . . . . . . .
Selecione a importac
ao atraves da biblioteca GDAL . . . . . . . . . . . . . . . . .
Opc
oes do comando r.in.gdal. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Insira um novo layer na pilha. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
As imagens SRTM, exibidas no Map Display . . . . . . . . . . . . . . . . . . . . .
Localize o comando r.patch. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Abas Required e Optional do comando r.patch. . . . . . . . . . . . . . . . . . .
Abas Command output e Manual do comando r.patch. . . . . . . . . . . . . .
Mapa raster resultante do comando r.patch. . . . . . . . . . . . . . . . . . . . . .
Localize a opc
ao de exibir os metadados do arquivo que criamos. . . . . . . . . . .
Metadados do arquivo criado com r.patch. . . . . . . . . . . . . . . . . . . . . . .
Localize o comando r.shaded.relief. . . . . . . . . . . . . . . . . . . . . . . . . .
Opc
oes do comando r.shaded.relief. . . . . . . . . . . . . . . . . . . . . . . . .
Resultado do comando r.shaded.relief, com resolucao original de 0:0:03 graus.
Resultado do comando r.shaded.relief, com resolucao de 0:0:30 graus. . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
27
28
28
29
29
30
31
32
34
35
36
37
37
38
38
39
39
40
40
41
41
42
42
43
44
44
44
46
47
48
49
50
50
51
52
52
53
54
54
56
57
57
59
LISTA DE FIGURAS
6.17
6.18
6.19
6.20
6.21
7.1
7.2
7.3
7.4
7.5
7.6
7.7
7.8
7.9
7.10
8.1
8.2
8.3
8.4
8.5
8.6
8.7
8.8
8.9
8.10
8.11
8.12
8.13
9.1
9.2
9.3
9.4
9.5
9.6
9.7
9.8
9.9
9.10
9.11
10.1
10.2
Selecione uma sub-

area da imagem original. . . . . . . . . . . . . . . . . . . . . . . . .
Display com resoluc
ao de 30 segundos (g.region res=0:0:30). . . . . . . . . . . . . .
Display com resoluc
ao de 05 minutos (g.region res=0:05). . . . . . . . . . . . . . . .
Display com resoluc
Display com resoluc
Interface gr
afica do QGIS com botao para adicionar dados vetoriais selecionado. . . .
Selecione o layer LM_UF. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
QGIS com layers LM_UF e VETOR_250000 carregados. . . . . . . . . . . . . . . . . . . .
Informac
oes sobre um dos polgonos do layer VETOR_250000. . . . . . . . . . . . . . . .
Janela de di
alogo das propriedades do layer. . . . . . . . . . . . . . . . . . . . . . . . .
Layer VETOR_250000 com transparencia no preenchimento. . . . . . . . . . . . . . . . .
Janela de informac
oes com dois resultados. . . . . . . . . . . . . . . . . . . . . . . . .
Recobrimento da base vetorial contnua em escala 1:250.000. . . . . . . . . . . . . . . .
Dados das
areas 1-2008 e 4-2009 exibidos no QGIS. Os dados da area 4-2009 estao no
alto, pr
oximos da seta do mouse. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dados das
areas 1-2008 e 4-2009 exibidos no QGIS com a projecao corrigida. . . . . .
Selecione a importac
ao atraves da biblioteca OGR . . . . . . . . . . . . . . . . . . . .
Importac
ao de shapefile pela interface simplificada do comando v.in.ogr. . . . . . . .
Opc
oes do comando v.in.ogr. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Vetores de hidrografia importados para o GRASS. . . . . . . . . . . . . . . . . . . . .
Detalhe do layer contorno_hd_permanente_LL, com areas preenchidas e centroides. .
Opc
oes do comando d.vect. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Convers
ao de tipos de geometrias vetoriais durante a importacao. . . . . . . . . . . . .
Opc
oes do comando v.random. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Pontos aleat
orios gerados com v.random. . . . . . . . . . . . . . . . . . . . . . . . . .
Mensagem de erro do gerenciador de tabelas de atributos do GRASS. . . . . . . . . .
Conex
ao de tabela com layer vetorial pelo gerenciador de tabelas do GRASS. . . . . .
Coluna elev atualizada com os valores do SRTM. . . . . . . . . . . . . . . . . . . . .
Mapa de pontos importado com v.in.asciisobreposto ao mapa gerado com v.random.
Reclassificac
ao dos dados SRTM em classes de 50 metros. . . . . . . . . . . . . . . . .
Mapa hipsometrico com classes de 50 metros e tabela de cores elevation. . . . . . . .
Escolha interativa das cores com r.colors.rules. . . . . . . . . . . . . . . . . . . . .
Adicione um layer de mapa sobreposto a relevo sombreado. . . . . . . . . . . . . . . .
Mapa hipsometrico sobreposto a relevo sombreado. . . . . . . . . . . . . . . . . . . . .
Mapa de declividade. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Mapa de aspecto (orientac
ao de vertentes). . . . . . . . . . . . . . . . . . . . . . . . .
Janela do r.mapcalculator. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Mapas aspecto com escala de cores aspect (esquerda) e definida manualmente (direita).
Relac
oes entre as curvaturas Vertical e Horizontal . . . . . . . . . . . . . . . . . . . . .
Mapas de curvatura vertical (profile curvature, esquerda) e horizontal (tangential curvature, direita). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Reamostragem (mudanca na resolucao espacial) por vizinhos mais proximos. . . . . .
Reamostragem (transformac
ao geometrica) por vizinhos mais proximos. . . . . . . . .
59
61
61
62
62
64
65
65
66
67
68
68
69
70
72
73
74
75
75
76
76
77
79
79
80
81
82
84
86
87
88
88
89
90
90
92
93
94
94
95
96
vi
LISTA DE TABELAS
10.3 Esquema da interpolac

ao por IDW. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10.4 Histogramas dos mapas de aspecto, a partir dos dados SRTM originais (esquerda) e
reamostrados com r.resamp.interp (direita). . . . . . . . . . . . . . . . . . . . . . .
10.5 Curvas de nvel para cotas definidas entre 0 m e 540 m (esquerda) e comparacao com
curvas de escala 1:250.000 do IBGE (direita, curvas do IBGE em vermelho). . . . . . .
10.6 Operac
oes de overlay. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10.7 Resultado da interpolac
ao por RST (esquerda) e modelo SRTM (direita). . . . . . . .
10.8 Resultado da subtrac
ao entre o modelo SRTM e a interpolacao por RST. . . . . . . .
11.1 Icone para ajuste do sistema de coordenadas (CRS - Coordinate Reference System) no
QGIS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
12.1 Janela do NVIZ. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
12.2 Controles de iluminac
ao do NVIZ. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
12.3 Selecione para alterar a cor da superfcie. . . . . . . . . . . . . . . . . . . . . . . . . .
12.4 Sobreposic
ao do mapa de diferencas à superfcies topografica. . . . . . . . . . . . . . .
13.1 Altere a cor da linha do mapa vetorial. . . . . . . . . . . . . . . . . . . . . . . . . . . .
13.2 Mapas raster e vetorial exibidos no QGIS. . . . . . . . . . . . . . . . . . . . . . . . . .
13.3 Janela do GRASS tools. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
13.4 Mapas rasters disponveis e propriedades do layer srtm_diff_rst. . . . . . . . . . . .
14.1 Janela do comando ps.map . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
14.2 Mapa visualizado no Documento Viewer . . . . . . . . . . . . . . . . . . . . . . . . . .
96
99
100
101
103
104
105
107
108
108
109
110
111
111
112
114
114
Lista de Tabelas
1
2
3
4
5
6
7
Estrutura dos comandos do GRASS . . . . . . . . . . . . .

Comandos do GRASS para importacao de mapas raster . .
Comandos do GRASS para importacao de mapas vetoriais .
Dimens
oes de alguns elips
oides comumente utilizados . . . .
Alguns datums comumente utilizados . . . . . . . . . . . . .
Variac
ao dos comprimentos de graus de latitude e longitude
Nomenclaturas de mapas no Sistema Cartografico Nacional
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
5
10
12
16
17
20
34
1 - Introduc
ao
1.
Introduc
ao
Hoje em dia, programas Livres e de codigo aberto se tornaram sinonimos de inovacao e
progresso1 . A liberdade de uso, distribuicao e modificacao dos programas e de seus codigosfonte asseguram o interc
ambio de ideias entre usuarios e desenvolvedores, o que reflete em um
sistema est
avel, com atualizac
oes constantes e orientadas à demanda dos usuarios.
O GRASS-GIS (Geographic Resources Analysis Support System), e um Sistema de Informacoes Geogr
aficas (SIG), com m
odulos para processamento e georreferenciamento de imagens multiespectrais, manipulac
ao de dados em formato raster e vetorial, com interfaces para
bancos de dados PostgreSQL, MySQL, SQLite, DBF, e ODBC, alem de poder ser conectado
a UMN/Mapserver, R-stats, gstat, Matlab, Octave, Povray, Paraview e outros programas.
O QGIS (Quantum GIS ) e um SIG com interface amigavel e multiplataforma, com suporte a diversos formatos de arquivos raster e vetorial (atraves das bibliotecas GDAL/OGR) e
de bancos de dados. Suporta leitura de arquivos diretamente do GRASS e possui uma toolbox
que permite rodar diversos comandos do GRASS, funcionando tambem como uma interface
grafica alternativa `
a interface nativa do GRASS.
Este tutorial tem como objetivo introduzir o usuario ao uso do GRASS-GIS e do QGIS,
com foco principal no GRASS-GIS.
Os exemplos foram produzidos com dados da regiao de Belo Horizonte e Catas Altas
(MG) e podem ser obtidos em http://carlosgrohmann.com/. Neste tutorial foi utilizado o
GRASS-GIS vers
ao 6.4 e QGIS vers
ao 1.7.4 em ambiente Linux; diferencas na interface grafica
podem ocorrer de acordo com a versao do programa e o sistema operacional utilizados. O
ambiente utilizado foi o OSGeo-Live 5.5 2 , um sistema baseado no Xubuntu Linux3 , que pode
ser executado a partir de um DVD, drive USB ou via maquina virtual. O sistema traz diversos
aplicativos pre-configurados para an
alises geoespaciais, como visualizacao, armazenagem em
bancos de dados, manipulac
ao de dados e producao de mapas.
1.1.
Fontes de informac
ao
Algumas fontes de informac

ao, manuais e tutoriais que podem ajudar a aprender e usar
o GRASS-GIS e o QGIS.
GRASS-GIS web site oficial: http://grass.osgeo.org
GRASS Wiki http://grass.osgeo.org/wiki/
Lista de discuss
ao eletr
onica: http://grass.osgeo.org/community/support.php
1
Leia mais sobre Software Livre em http://www.gnu.org

http://live.osgeo.org/en/index.html
3
http://xubuntu.org/
2
1.2
Design e estrutura de SIGs
QGIS web site oficial: http://qgis.org

QGIS Wiki http://www.qgis.org/wiki/
Lista de discuss
ao eletr
onica: http://qgis.org/en/community/mailing-lists.html
Dassau et al. (2005).
An introduction to the practical use of the Free Geographi-
cal Information System GRASS 6.0, GDF Hannover bR, disponvel em http://www.
gdf-hannover.de/literature.
Neteler & Mitasova (2008). Open Source GIS: A GRASS GIS Approach. Boston, Kluwer
Academic Publishers.
Grohmann (2008). Introduc
ao ao geoprocessamento e `
a an
alise digital de terreno com
software livre. Instituto de Geociencias, Universidade de Sao Paulo, disponvel em
http://carlosgrohmann.com/.
Nanni et al. (2012). Quantum GIS - Guia do Usu
ario, Vers
ao 1.7.4 Wroclaw, disponvel
em http://qgisbrasil.files.wordpress.com/2012/05/guia_do_usuario_174_pt_
br1.pdf
1.2.
Design e estrutura de Sistemas de Informac

ao Geogr
afica
De maneira geral, um Sistema de Informacoes Geograficas pode ser caracterizado por

quatro componentes principais:
Entrada de dados
Administrac
ao
Analise
Apresentac
ao
Os tipos de dados nesse sistema de quatro componentes podem ser classificados em geom
etricos, atributos e gr
aficos:
Dados geom
etricos descrevem a situacao espacial dos objetos, sua forma e sua posicao
no espaco. As relac
oes entre pontos individuais, linhas e areas sao dadas por um sistema de
coordenadas de referencia, e sua relacoes com o mundo real. Dados geometricos podem ser do
tipo raster (matricial, pixel ) ou vetorial (pontos, linhas, polgonos, areas, solidos - ver figura
1.1):
1.2
Design e estrutura de SIGs
Dados raster s
ao dados regularmente espacados no espaco, em uma estrutura de matriz com
celulas quadradas (normalmente) e de mesmo tamanho. Este tipo de dado e utilizado
em imagens de satelite, modelos de relevo e mapas escaneados, por exemplo.
Dados vetoriais s
ao usados para informacoes de objetos lineares (linhas e polilinhas) ou areas
definidas por linhas fechadas (polgonos).
Dados pontuais podem ser considerados uma forma especial de dados vetoriais (pontos vetoriais), e s
ao usados para armazenar informacoes pontuais espalhadas em uma area,
como pontos de coletas de amostras, localizacao de stios arqueologicos, etc.
Atributos (categorias) s
ao, como o nome sugere, atributos relacionados aos tipos de
dados mencionados anteriormente, e normalmente sao armazenados em um banco de dados
interno ao sistema SIG ou em um banco de dados externo atraves de uma interface tipo DBMI
(Database Management Interface).
Gr
aficos descrevem os metodos usados para desenhar os objetos espaciais nos varios
dispositivos de sada existentes (tela do monitor, impressora, etc).
Figura 1.1. Geometia e atributos de dados em SIG
1.3
Dimens
ao de dados em SIG
Em um SIG, as relac
oes entre esses tipos de dados (geometricos, atributos e graficos)
definem duas estruturas b
asicas de dados:
Estruturas contnuas objetos laminares e ilimitados no espaco (superfcies)
Estruturas discretas objetos definidos como linhas e areas
Em um SIG, m
odulos especficos permitem a conversao entre estruturas diferentes. Curvas
de nvel, por exemplo, podem ser armazenadas como vetores (estruturas discretas) e podem ser
convertidos para uma representac
ao contnua do terreno (raster) por interpolacao. Nesse caso
a qualidade do modelo resultante ir
a depender da qualidade e resolucao dos dados originais
(Fig. 1.2).
Figura 1.2. Representac

ao do relevo de uma a
rea por estrutura discreta (curvas de nvel vetoriais) e
contnua (superfcie raster gerada por interpolaca
o).
1.3.
Dimens
ao de dados em SIG
Dados espaciais podem ser representados por pontos, linhas,

areas/superfcies e volumes
(Fig. 1.3).
Pontos podem representar a localizacao de objetos reais (como uma casa) ou virtuais
(como o epicentro de um terremoto), e nao tem dimensao fsica. Linhas tem apenas uma
dimensao fsica, o comprimento, enquanto que as areas tem duas dimensoes (2D: comprimento
e largura). No caso de superfcies contnuas, como em um modelo da superfcie topografica,
onde ha um terceiro par
ametro alem das coordenadas X e Y (Z=elevacao, por exemplo),
dizemos que os dados tem 2,5D. Isto pode parecer um pouco estranho, ja que os programas
1.4
Estrutura dos comandos no GRASS
de SIG nos permitem visualizar as superfcies em perspectivas tri-dimensionais, porem apesar

do parametro Z ser utilizado para gerar um imagem 3D, a superfcie nao possui um valor de
espessura (a terceira dimens
ao fsica). Os u
nicos dados realmente tri-dimensionais sao aqueles
que possuem quatro par
ametros (XYZ + atributo) como no caso de solidos representando
edifcios ou em modelos que utilizam voxels, que sao o equivalente dos pixels em tres dimensoes,
representados por cubos no espaco contendo um valor de atributo .
Figura 1.3. Dimens

ao de dados em SIG
1.4.
Estrutura dos comandos no GRASS
O GRASS-GIS e um Sistema de Informacoes Geograficas com funcoes orientadas a dados

do tipo raster ou vetoriais. Cada funcao e conduzida por um modulo proprio, o que faz
do GRASS um SIG bem estruturado e transparente. Alem disso, como apenas os modulos
necessarios est
ao em execuc
ao, os recursos do computador sao preservados, e garantem um
melhor desempenho do sistema.
A estrutura dos comandos do GRASS e bastante simples e clara. Os comandos (chamados
de m
odulos no GRASS) tem nomes auto-explicativos e seu tipo pode ser reconhecido pelo
prefixo. Por exemplo, o m
odulo para obter informacoes sobre mapas raster e r.info e o
modulo para converter de formato vetorial para raster e o v.to.rast. Na tabela 1 vemos a
estrutura geral dos nomes dos m
odulos em mais detalhe.
Alem dos m
odulos existentes no GRASS, todos os programas do ambiente Unix/Linux
estao disponveis com uso de linhas de comando na janela de terminal (ou de console). Apesar
de muitos usu
arios iniciantes considerarem difcil usar estas funcionalidades, a programacao de
shell scripts e simples de se aprender, e pode ser usada para automatizar tarefas ou extender
as funcionalidades do ambiente SIG de acordo com as necessidades de cada um.
1.5.
Obtendo ajuda com relac

ao aos m
odulos
Arquivos de ajuda est

ao disponveis para quase todos os 400 modulos do GRASS, onde sao
descritos a funcionalidade do m
odulo e a sintaxe de uso. Uma ajuda breve pode ser obtido com
o uso do par
ametro -help, na linha de comando, como em d.rast -help. Descricoes mais
1.6
Organizac
ao dos projetos no GRASS-GIS
Prefixo
Func
ao
Significado dos comandos
d.*
r.*
i.*
v.*
g.*
ps.*
db.*
r3.*
display
raster
imagem
vetor
geral
postscript
database
voxel
Exibicao e consulta grafica no Map Display

Processamento de dados raster
Processamento de imagens
Processamento de dados vetoriais
Comandos gerais para manutencao de arquivos
Criacao de mapas em formato postscript
Gerenciamento de banco de dados
Processamento de dados raster 3D
Tabela 1. Estrutura dos comandos do GRASS
detalhadas encontram-se nas p

aginas de manual para cada modulo. Essas paginas podem ser
acessadas on-line no web site do GRASS ou com o comando g.manual, como em g.manual
d.rast, ou clicando no bot
ao Help na janela de cada comando. A pagina de manual ser
a
aberta em um navegador de internet.
1.6.
Organizac
ao dos projetos no GRASS-GIS
A organizac
ao dos projetos no GRASS-GIS segue uma hieraquia baseada em Locations
e mapsets. A Location compreende toda a area de trabalho do projeto, enquanto que o
mapset e a porc
ao ativa e utilizada para analise, que pode ser do mesmo tamanho ou menor
que a location. Ainda, v
arios mapsets podem ser definidos para a mesma location.
Em ambientes multiusu
ario, v
arias pessoas podem trabalhar na mesma Location ao
mesmo tempo, mas n
ao no mesmo mapset. Para contornar essa limitacao, dados de interesse comum (tais como modelos de relevo, imagens de satelite etc) podem ser armazenados
em um mapset especial ao qual todos os usuarios tem acesso, chamado de PERMANENT, que e
criado automaticamente ao se criar uma nova Location.
Um conceito importante dentro do GRASS e o de region, que define, dentro do mapset,
a area de interesse e a resoluc
ao espacial dos mapas raster. Tanto a resolucao espacial quanto
as coordenadas do ret
angulo envolvente da region podem ser facilmente alteradas sem a
necessidade de reinicializac
ao do sistema ou a criacao de novos projetos; e possivel salvar as
configurac
oes da region para acess
a-la facilmente quando necessario.
preciso frisar que todas as an
E
alises envolvendo mapas raster (analise de terreno, algebra
de mapas, interpolac
ao de superfcies etc) sao efetuadas de acordo com as configuracoes da
region ativa, e que esta n
ao necessariamente corresponde com as configuracoes do Display.
Os dados referentes aos projetos sao armazenados em um diretorio (pasta) chamado
pelo GRASS de GISBASE. As Locations serao subdiretorios de GISBASE, os mapsets serao
subdiretorios dentro de cada Location e assim por diante.
1.7
O mapset PERMANENT
Figura 1.4. Exemplo da estrutura de projetos no GRASS-GIS, extrado de Dassau et al. (2005).
Normalmente o diret
orio GISBASE (p.ex., grassdata) e criado antes de se rodar o GRASS
pela primeira vez:
> cd <Enter>
> mkdir grassdata <Enter>
No Linux, o comando cd (change directory) leva ao diretorio inicial do usuario. Com mkdir
(make directory), o diret
orio grassdata e criado.
Uma vez que cada parte dos layers (p.ex., geometria, tabela de atributos, objetos graficos) e armazenada em subdiret
orios diferentes, a manutencao dos dados dos projetos (p.ex.,
copiar, apagar, renomear) deve ser feita apenas com os comandos apropriados (g.copy,
g.remove, g.rename).
1.7.
O mapset PERMANENT
Todas as informac
oes sobre projecao cartografica, resolucao espacial e extensao da area
do projeto s
ao armazenadas no mapset PERMANENT, que e gerado automaticamente pelo
GRASS ao se criar uma Location. Se necessario dados de interesse comum podem ser armazenados neste mapset. Como apenas o usuario que criou a Location tem permissao de
escrita ao mapset PERMANENT, tem-se a garantia de que os dados nao serao apagados ou
1.8
Organizac
ao dos projetos no QGIS
modificados acidentalmente. Outros usuario (e tambem a pessoa com permissao de escrita

ao mapset PERMANENT) devem criar mapsets adicionais para criar, salvar e alterar seus
proprios arquivos e an
alises baseadas nos dados do mapset PERMANENT.
1.7.1
Arquivos do mapset PERMANENT
Se nenhum mapa for armazenado em PERMANENT, o diretorio ira conter apenas informacoes sobre o projeto:
DEFAULT WIND Especificacoes dos limites e resolucao espacial da region padrao
(default).;
MYNAME Nome do projeto. p.ex., doutorado, relevoAmazonia;
PROJ INFO Especificac
oes sobre a projecao cartografica, datum e elipsoide;
PROJ UNITS Unidade de medidas (metros, graus);
WIND Especificac
oes dos limites e resolucao espacial da region atual.
1.8.
Organizac
ao dos projetos no QGIS
No QGIS n
ao e necess
ario importar arquivos. Os dados sao acessados diretamente, em
arquivos locais ou via rede, pela biblioteca GDAL/OGR. A [des]organizacao dos dados fica a
cargo do usu
ario, tal como no ArcGIS ou gvSIG.
O QGIS tambem e capaz de realizar projecoes entre datuns e sistemas de coordenadas
em tempo real (projec
ao on-the-fly), nao sendo necessario converter todos os dados para
um datum e sistema de coordenadas u
nico, como no caso do GRASS. Ao mesmo tempo em
que a projec
ao on-the-fly e muito u
til, a obrigatoriedade da conversao dos dados para um
referencial cartogr
afico u
nico e vista por muitos como benefica, pois minimiza a possibilidade
de problemas nas an
alises.
2 - Mapas Raster - vis

ao geral
2.
Mapas Raster - vis

ao geral
Mapas raster s
ao compostos por celulas, sendo que cada celula (ou pixel: picture element)
recebe o valor de um atributo, que representa um fenomeno (por exemplo temperatura ou

altitude). A resoluc
ao espacial de um raster equivale à dimensao da celula no mundo real.
Quanto maior a resoluc
ao espacial, maior e o nvel de detalhe que o raster pode exibir, e
portanto menor o tamanho do pixel no mundo real (Fig. 2.1).
Figura 2.1. Resoluca

o espacial de um raster
As celulas s
ao organizadas em linhas e colunas, e seu valor pode ser acessado pelas coordenadas absolutas da matriz (linha/coluna) ou pelas coordenadas geograficas. Mapas raster
sao normalmente utilizados para aplicacoes de sensoriamento remoto como analise de imagens orbitais ou fotografias aereas e para interpolacao de dados irregularmente distribudos
no espaco, tal como a gerac
ao de Modelos Digitais de Terreno (MDT) a partir de pontos
coletados por GPS ou por levantamentos a laser aeroportados (LiDAR).
Uma das desvantagens de se trabalhar com dados raster e a quantidade de memoria
necessaria e o tempo de processamento das analises, que cresce exponencialmente com o
aumento da resoluc
ao espacial utilizada (particularmente no caso de modelos de elevacao de
alta resoluc
ao como os levantamentos a laser), porem com o contnuo aumento na capacidade
de processamento e de armazenamento dos compuatdores atuais, essa questao ja nao e mais
considerada como essencial. Outro ponto a ser lembrado e que dados raster nao possuem
relacoes de vizinhanca (Topologia: o polgono A esta à direita do polgono B, e assim por
diante), uma vez que cada pixel e definido de maneira independente dos outros.
No GRASS e preciso importar arquivos raster, que serao convertidos para um formato
proprio. No QGIS, os arquivos s
ao acessados diretamente, sem necessidade de importacao. O
GRASS pode importar diversos tipos de mapas raster. De maneira geral, podemos agrupa-los
em tres categorias:
3 - Mapas Vetoriais - vis

ao geral
10
Imagem: Rasters em formato de imagem sempre possuem valores positivos e inteiros, tais
como os formatos PPM, PNG, JPEG, e GIF.
ASCII: O formato ASCII podem conter valores positivos, negativos, inteiros ou decimais.
Os arquivos ASCII-GRID do ArcInfo sao exemplos desse formato.
Bin
ario: Rasters bin
arios podem ter valores inteiros ou decimais, positivos ou negativos, que
podem estar em diferentes canais (bandas), em diferentes resolucoes espaciais. Exemplos
sao os arquivos [Geo]TIFF e ERDAS/IMG.
Ao se importar um mapa raster no GRASS, deve-se sempre lembrar que os mapas sao
importados com sua extens
ao e resolucao originais, mas exportados de acordo com a extensao
e resolucao definidas pela region ativa.
A tabela 2 lista os comandos usados para importar diferentes formatos raster:
Comando do GRASS
Formato Raster
r.in.ascii
r.in.bin
r.in.gdal
GRASS ASCII
BIL, arquivos binarios do GMT, LANDSAT TM5/ETM+
ArcInfo ASCII/Binary GRID, BIL, ERDAS (LAN, IMG),
USGS DOQ, JPEG, SAR CEOS, EOSAT, GeoTIFF,
PPM/PNM, SDTS DEM, GIF, PNG
arquivos SRTM em formato HGT (usa GDAL)
ESRI ARC/INFO ascii raster
Imagens Terra-ASTER e ASTER DEMs
r.in.srtm
r.in.arc
r.in.aster
Tabela 2. Comandos do GRASS para importaca

o de mapas raster
O modulo mais utilizado para importacao (e exportacao) de rasters e certamente o

r.in.gdal, que usa a biblioteca GDAL (Geospatial Data Abstraction Library 4 ).
3.
Mapas Vetoriais - vis

ao geral
Mapas vetoriais representam linhas, polilinhas (polylines) ou polgonos.
Uma linha conecta dois pontos extremos (end points nos ou vertices), cada um com suas
respectivas coordenadas XY ou XYZ. Cada objeto vetorial pode possuir varios (ou nenhum)
atributos, armazenados em um banco de dados. Para uso interno ao GRASS, os bancos de
dados dBase e sqlite est
ao disponveis. Existem interfaces para acesso à DBMS (Database
Management System) externos, como PostgreSQL, Oracle, mySQL, Access, etc. Para o armazenamento externo ao GRASS de dados geometricos e atributos, existem interfaces para
PostGIS ou Oracle Spatial.
Em comparac
ao com dados raster, os dados vetoriais utilizam menos memoria e tem menor
tempo de processamento na maioria das analises. Alem disso, no GRASS, os vetores possuem
4
http://www.gdal.org/
3 - Mapas Vetoriais - vis

ao geral
11
relacoes de topologia, o que significa que cada objeto vetorial sabe quais vertices possui e/ou
com quais
areas faz vizinhanca.
Mais precisamente, isso significa que elementos de mapas vetoriais compartilham componentes. Por exemplo, em um sistema nao topologico, se duas areas (polgonos) possuem uma
borda (aresta) em comum, essa borda teria que ser digitalizada e armazenada duas vezes, o
que nao ocorrem em um sistema topologico, onde a borda em comum e compartilhada pelas duas areas. Esse tipo de representacao ajuda a manter mapas com uma geometria mais
limpa e permite an
alises que n
ao podem ser realizadas com dados nao topologicos (tambem
chamados de vetores tipo espaguete, uma vez que cada polgono tem que ser digitalizado
inteiro). No GRASS-GIS, vetores com topologia sao referidos como de nvel 2, e dados tipo
espaguete s
ao de nvel 1.
Entretanto, a informac
ao de topologia nem sempre e necessaria, e os requisitos de memoria
e espaco em disco podem se tornar impraticaveis com determinados tipos de dados (como os
resultantes de levantamentos topograficos por laser aeroportados - LiDAR). Sendo assim,
existem dois m
odulos do GRASS que permitem trabalhar com dados vetoriais de nvel 1 (sem
topologia). O m
odulo v.in.ascii e usado para importacao de dados pontuais e pode criar
um layer vetorial sem construir a tabela de topologia. Normalmente e usado em conjunto
com v.surf.rst, que realiza a interpolacao de valores a partir de pontos ou isolinhas.
Os seguinte objetos vetoriais est
ao definidos no GRASS-GIS:
ponto: um ponto no espaco;
linha: uma seq
uencia de pontos (vertices) conectados, com dois pontos extremos (endpoints)
chamados de n
os (nodes);
borda (boundary): a linha que define uma area;
centr
oide: um ponto dentro de uma borda fechada;
area:
a composic
ao topol
ogica de borda + centr
oide;
face:
uma
area tridimensional;
n
ucleo (kernel): um centr
oide tridimensional dentro de um volume (ainda nao implementado);
volume: um corpo tridimensional, a composicao topologica das faces + n
ucleo (ainda nao
implementado);
O modulo v.type pode ser usado para conversao entre tipos vetoriais (se possvel). O
modulo v.build e usado pra gerar topologia, e opcionalmente permite a extracao de objetos
esp
urios em um novo mapa. Erros de topologia podem ser corrigidos manualmente dentro do
3.1
12
Manipulac
ao de atributos
ambiente do m
odulo v.digit ou, ate certo ponto, automaticamente com v.clean. Polgonos
adjacentes pode ser detectados com v.to.db (opcao sides).
Ao se importar um mapa vetorial, e preciso considerar que existem diversos formatos e
padroes diferentes de arquivos. A tabela 3 lista os comandos usados para importacao de dados
vetoriais. Provavelmente o formato mais utilizado e o ESRI SHAPE (o famoso shapefile). O
modulo v.in.ogr e o mais utilizado, por oferecer uma interface com a maioria dos formatos
vetoriais. Adicionalmente, esse m
odulo tem opcoes para criar novas Locations ou estender a
region padr
ao de acordo com os limites do mapa a ser importado. Outros modulos existem
para casos mais especficos, como o v.in.ascii, usado para importar dados de arquivos texto
com as coordenadas e os atributos de pontos no espaco (como no caso de dados LiDAR), e o
v.in.db, para importar a partir de um banco de dados contendo as coordenadas e atributos
dos dados.
O modulo v.external permite que mapas externos sejam ligados virtualmente ao mapset,
nao necessitando de importac
ao (e consequente conversao para o formato nativo do GRASS).
Apenas uma pseudo-topologia e gerada e a geometria do vetor nao e importada. Os comandos
v.out.* s
ao usados para exportar os mapas vetoriais em varios formatos.
Comando do GRASS
Formato Vetorial
v.in.ogr
SHAPE file, UK.NTF, SDTS, TIGER, S57, MapInfo-File,

DGN, VRT, AVCBin, REC, Memory, GML, ODBC
(veja tambem: http://www.gdal.org/ogr/ogr_formats.html)
Arquivos texto em formato x|y[|z] e
vetores ASCII exportados por versoes antigas do GRASS
Arquivos DXF (CAD Desktop eXchange Format)
Importa Waypoints, Tracks e Routes de um GPS Garmin
ArcInfo-E00
Cria vetores a partir de bancos de dados com coordenadas x|y[|z]
v.in.ascii
v.in.dxf
v.in.garmin
v.in.e00
v.in.db
Tabela 3. Comandos do GRASS para importaca

o de mapas vetoriais
3.1.
Manipulac
ao de atributos
Apos a importac
ao, o mapa estara disponvel no formato vetorial nativo do GRASS.
A geometria, topologia e os atributos sao armazenados em um banco de dados interno do
GRASS (em formato DBF por padr
ao). O GRASS pode ser conectado a um ou varios bancos
de dados. Os comandos bd.* fornecem funcionalidades basicas de SQL (Structured Query
Language) para manutenc
ao de atributos, enquanto que os comandos v.db.* operam nos
mapas vetoriais,
Categorias: o n
umero da categoria e o ID (IDentificador) do vetor e e usado para
conectar (link) atributo(s) `
a cada objeto vetorial (cada objeto vetorial pode ter zero,
uma, duas ou mais categorias). Os n
umeros de categoria sao armazenados tanto no
3.1
Manipulac
ao de atributos
13
arquivo de geometria quanto na(s) tabela(s) de atributos para cada objeto vetorial
(normalmente e a coluna cat). A visualizacao e manutencao dos n
umeros de categoria
pode ser feita com v.category. Para poder ligar um vetor à mais de uma tabela de
atributos, mais de um n
umero de categoria e necessario.
Layers (camadas): e possvel ligar (link ) os objetos geograficos em um mapa vetorial a
uma ou mais tabelas. Cada link com uma tabela diferente e chamado de layer. Um
link define qual o driver de acesso ao banco de dados (MySQL, PostgreSQL, etc), qual
o banco de dados e qual tabela desse banco de dados sera utilizada. Cada n
umero de
categoria no arquivo de geometria corresponde à um linha na tabela de atributos (a
coluna de ligac
ao e normalmente a coluna cat). A visualizacao e manutencao dos
layers pode ser feita com v.db.connect.
Os layers do GRASS n
ao contem nenhum objeto geografico, mas na verdade consistem
de links para tabelas de atributos nas quais os objetos podem ter zero, uma ou mais categorias. Se um objeto vetorial tem zero categorias em um layer, entao ele nao e exibido
nesse layer. Os objetos vetoriais podem, portanto, aparecer em um determinado layer,
mas n
ao em outro(s). Esse sistema permite que objetos topologicamente conectados,
mas tematicamente distintos (como florestas e lagos, por exemplo) sejam colocados em
um mesmo mapa. Outra possibilidade e a ligacao de atributos temporais à localidades
que n
ao mudam com o passar do tempo.
Por padr
ao, o primeiro layer e ativo, ou seja, a primeira tabela corresponde ao primeiro
layer. Outras tabelas s
ao ligadas aos layers seguintes.
Suporte `
a SQL: o driver ao banco de dados DBF (padrao interno do GRASS), prove
apenas funcionalidades limitadas de SQL (uma vez que o DBF nao e um banco de dados SQL), enquanto que outros DBMS externos ao GRASS como PostgreSQL, MySQL,
SQLite etc, proveem funcionalidades de SQL completas, ja que os comandos SQL sao
enviados diretamente `
a interface do banco de dados. Comandos SQL podem ser executados diretamente com db.execute, db.select e outros modulos db.*.
14
4 - Elementos de cartografia
4.
Elementos de cartografia
Antes de criar um projeto e preciso definir qual projecao cartografica sera utilizada. Alem
disso, no GRASS n
ao e possvel fazer transformacoes de coordenadas on-the-fly, ou seja,
utilizar dados em diferentes sistemas de coordenadas ao mesmo tempo. Nas secoes seguintes
veremos brevemente alguns conceitos de geodesia, cartografia e os parametros das projecoes
mais utilizadas.
4.1.
Escala
A escala de um mapa define a relacao entre a medida de um objeto ou lugar representado

no papel e sua medida real:
E=
d
D
onde d e a dist
ancia medida no mapa e D a distancia real. A escala pode ser representada
em sua forma numerica, como em 1/10.000 e 1:10.000 ou de forma gr
afica. Uma escala
de 1:10.000 significa que os objetos e distancias reais sao 10.000 vezes maiores do que suas
representac
oes em mapa (ou seja, 1 cm no mapa equivale a 100 m ou 10.000 cm no
terreno).
A escala gr
afica (fig. 4.1) e uma representacao de varias distancias do terreno sobre uma
linha reta graduada, de acordo com a unidade de medida utilizada (quilometros, milhas, etc).
O segmento `
a direita do valor zero e chamado de escala prim
aria, e o segmento à esquerda
do zero recebe o nome de tal
ao ou escala de fracionamento, e e divido em sub-m
ultiplos da
unidade de medida utilizada, graduado da direita para a esquerda.
0
10
20
0
5
talo
30
40
50 km
10
20
30
40
50 milhas
10
15
20 km
escala primria
Figura 4.1. Exemplos de escalas gr

aficas.
A escala e uma divis

ao matem
atica entre os valores real e em mapa, portanto uma escala
de 1:1.000 e maior que 1:100.000:
1
= 0.001
1000
>
1
= 0.00001
100000
Em outras palavras, quanto maior a escala do mapa, maior o nvel de detalhe apresentado.
4.2
15
Ge
oide, elips
oide e datum
4.2.
Ge
oide, elips
oide e datum
O geoide e considerado como uma representacao matematica bastante precisa da forma da

a superfcie equipotencial do campo gravitacional da Terra, que coincide
Terra (fig. 4.2). E
com o nvel medio dos oceanos (fig. 4.3). A determinacao dos parametros do geoide e alvo
de estudo da geodesia, a ciencia que trata do levantamento e da representacao da forma e da
superfcies da Terra, incluindo o campo gravitacional. Devido à complexidade matematica do
geoide, a forma da Terra e geralmente representada por elipsoides.
Figura 4.2. Representac

ao das ondulaco
es do ge
oide, construda a partir de dados de anomalia gravimetrica (NASA/JPL/University of Texas Center for Space Research).
Geide
Elipside
Continente
Oceano
Figura 4.3. Relaco

es entre a superfcie da Terra, o ge
oide e o elips
oide.
Como a simplificac
ao da forma da Terra para uma esfera nao e precisa a suficiente para
gerar mapas em escala maior (com maior detalhe) do que 1:2.000.000, e necessario utilizar
elips
oides de revoluc
ao, que possuem achatamento nos polos, para tentar definir uma superfcie
mais simples de se trabalhar do que o geoide, e ainda assim precisa o suficiente para fins
cartograficos.
Matematicamente, o elips
oide de referencia e normalmente um esferoide oblato (achatado)
com dois eixos diferentes, o raio equatorial ou semi-eixo maior (a) e o raio polar ou semi-
4.2
16
Ge
oide, elips
oide e datum
eixo menor (b). As relac

oes entre os semi-eixos definem as propriedades de achatamento
(flattening, f ) e ecentricidade (e):
f=
(a b)
(b)
=1
a
a
e=
a2 b2
a
Usualmente, os elips
oides s
ao definidos em termos do semi-eixo maior (a) e do inverso do
achatamento (1/f ). A figura 4.4 ilustra os elementos do elipsoide de referencia e a tabela 4
lista os par
ametros de alguns elips
oides comumente utilizados.
Eixo menor
Esfera
Elipside
Semi-eixo
menor (b)
Eixo maior
Semi-eixo
maior (a)
Figura 4.4. Elementos do elips

oide de referencia.
Elips
oide
Semi-eixo maior a (m)
Semi-eixo menor b (m)
Achatamento (1/f )
Clarke 1866
Internacional 1924
GRS 1980
WGS 1984
SGR 1967
Esfera (6371 km)
6
6
6
6
6
6
6
6
6
6
6
6
1/294,9786982
1/297,0
1/298,257222101
1/298,257223563
1/298,25
378
378
378
378
378
371
206,4
388,0
137,0
137,0
160,0
000,0
356
356
356
356
356
371
583,8
911,9
752,3141
752,3142
776,0
000,0
Tabela 4. Dimens
oes de alguns elips
oides comumente utilizados
Um datum (no plural, datums ou data), e uma referencia a partir da qual sao realizadas
medidas. Em cartografia datums s
ao utilizados para definir sistemas de cordenadas, sejam
eles globais ou locais. Existem datums topocentricos, cuja origem e um ponto na superfcie
terrestre, e datums geocentricos, com origem no centro de massa da Terra.
4.2
17
Ge
oide, elips
oide e datum
Na figura 4.5 o datum NAD27 (North American Datum 1927 ) e do tipo topocentrico, e
define um sistema de coordenadas local. Veja que o elipsoide se ajusta bem à area da America
do Norte, mas n
ao se ajusta `
a regi
ao da Australia. O datum WGS84 (World Geodetic System
1984 ) e do tipo geocentrico, e busca um ajuste medio ao planeta. Por definir um sistema de
coordenadas global, e utilizado no mundo todo.
superfcie da
Terra
Amrica do
Norte
NAD27
WGS84
Austrlia
NAD27 - Datum topocntrico

Sistema de coordenadas local
WGS84 - Datum geocntrico

Sistema de coordenadas global
Figura 4.5. Datum topocentrico (NAD27) e geocentrico (WGS84).
O datum SAD-69 (South American Datum 1969 ), amplamente utilizado no Brasil, tem
sua origem no vertice do Sistema Geodesico Brasileiro localizado no municpio de Chua (MG).
Recentemente completamos uma fase de transicao do datum SAD-69 (topocentrico) para o
novo datum oficial do Brasil, o SIRGAS5 (Sistema de Referencia Geocentrico para a America
do Sul ), que como o nome indica, trata-se de um datum geocentrico.
O sistema de referencia do SIRGAS coincide com o IERS6 (International Earth rotation
and Reference systems Service) e com o ITRF7 (International Terrestrial Reference Frame),
e tem os par
ametros do elips
oide Geodetic Reference System (GRS-80), sendo considerado
identico ao WGS-84 para fins cartograficos. A tabela 4 lista as caractersticas de alguns
datuns comumente utilizados.
Datum
Regi
ao de uso
Origem
Elips
oide
WGS 84
NAD 83
ED 50
SAD 69
SIRGAS
Global
America do Norte, Caribe
Europa, Africa
(norte)
America do Sul
America da Sul
Centro de massa da Terra

Potsdam
Chu
a
WGS 84
GRS 80
Internacional 1924
SGR 67
GRS 80
Tabela 5. Alguns datums comumente utilizados
http://www.ibge.gov.br/sirgas/
http://www.iers.org/
7
http://itrf.ensg.ign.fr/
6
4.3
18
Sistemas de coordenadas
4.3.
Sistemas de coordenadas
De modo geral, podemos agrup

a-los em sistemas globais e locais. Sistemas globais requerem coordenadas esfericas, enquanto que sistemas locais podem trabalhar com coordenadas
planas em diferentes unidades de medida (metros, quilometros, pes, milhas etc).
4.4.
Coordenadas planas
O sistema de coordenadas planas e o mais simples de utilizar. Assim como em um grafico

cartesiano tipo XY, a origem dos dados e determinada pelo cruzamento de dois eixos perpendiculares e os valores (coordenadas) aumentam para a direita (em X) e para cima (em Y). O
eixo horizontal e chamado de abcissa, e o eixo vertical, de ordenada (fig. 4.6).
xy = (3,2)
75.000 E, 50.000 N
100.000
cordenada norte (N)
150.000 E, 140.000 N
50.000
xy = (6,4)
ordenada
Coordenadas planas (metros, milhas etc)

150.000
Grfico cartesiano (XY)
4
abcissa
50.000
100.000
150.000
coordenada leste (E)
Figura 4.6. Sistema de coordenadas planas.
4.5.
Coordenadas geogr
Neste sistema os planos de referencia sao o Equador e o Meridiano de Greenwich. A Terra

e dividida em 180 de longitude a partir de Greenwich para Leste e para Oeste, e em 90 de
latitude para Sul e Norte do Equador (fig. 4.7). A altitude e medida a partir do centro de
massa da Terra. Em geodesia, a latitude e representada pela letra grega phi () e a longitude
por lambda ().
Um paralelo e o crculo imagin
ario que une pontos de mesma latitude, formado pela
interseccao de um plano paralelo ao Equador com a superfcie terrestre. Um meridiano e um
crculo imagin
ario perpendicular ao Equador e que passa pelos polos da Terra. Por convencao
internacional, determinou-se que o meridiano que passa pelo observatorio de Greenwich, na
Inglaterra, seja usado como referencia origem das coordenadas geograficas, sendo chamado de
meridiano origem, meridiano fundamental ou primeiro meridiano (prime meridian).
4.5
19
Coordenadas geogr
W
Eq
ad
or
Mer. G
r een
w ic
Latitude ()
60 E, 55 N
60 N
40 N
55
60
0
20 E 40 E
20 N
0
60 E
Longitude ()
S
Figura 4.7. Sistema de coordenadas geogr
aficas (latitude-longitude).
Quando falamos de latitude, normalmente estamos nos referindo à latitude geodetica (ou
latitude geogr
afica), que e definida como o angulo entre o plano equatorial e uma linha normal
ao elipsoide de referencia, que passe pelo ponto de interesse. O angulo entre o plano equatorial
e uma linha que ligue o ponto de interesse ao centro da Terra e chamado de latitude geocentrica,
representada pela letra grega psi () (fig. 4.8). A longitude geocentrica nao difere da longitude
geodetica, uma vez que os planos de mesma latitude (paralelos) podem ser considerados como
crculos perfeitos, ao contr
ario dos planos formados pelos meridianos.
Plo
Tangente
Elipside
P
Normal ao elipside
no ponto P
Equador
Latitude Geodtica
no ponto P ()
Longitude Geodtica
no ponto P ()
Latitude Geocntrica
no ponto P ()
Figura 4.8. Longitude, latitude geodetica e latitude geocentrica.
4.5
20
Coordenadas geogr
Considerando os par
ametros do elipsoide de referencia, podemos calcular o comprimento
do arco equivalente a 1 de latitude ou longitude. Se a Terra fosse uma esfera perfeita, a
distancia entre dois paralelos (o arco de meridiano) seria sempre a mesma, mas como aproximamos sua forma de um elips
oide, a curvatura da superfcie varia mais rapidamente proximo
do Equador do que nos p
olos, o que faz com que o grau de latitude proximo ao Equador seja
ligeiramente menor do que pr
oximo dos polos. O comprimento do grau de longitude varia de
acordo com o cosseno da latitude (tabela 6)8 , de aproximadamente 111,32 km no Equador
ate zero nos p
olos (fig. 4.9).
Figura 4.9. Variaca

o do comprimento do grau de longitude com a latitude.
latitude
arco de 1 de latitude (km)
arco de 1 de longitude (km)
0
10
20
30
40
50
60
70
80
90
110.574
110.608
110.704
110.852
111.035
111.229
111.412
111.562
111.660
111.694
111.319 km
109.639 km
104.647 km
96.486 km
85.394 km
71.696 km
55.800 km
30.187 km
19.393 km
0.000 km
km
km
km
km
km
km
km
km
km
km
Tabela 6. Variaca
o dos comprimentos de graus de latitude e longitude
A distancia mais curta entre dois pontos na superfcie de uma esfera (ou elipsoide) e dada
por um crculo m
aximo, formado pela interseccao entre a superfcie da esfera e um plano
que passa pelo centro da esfera. Todos os meridianos sao crculos maximos, assim como o
Equador. Um crculo mnimo e definido por um plano que intersecta a esfera mas nao passa
pelo seu centro. Todos os paralelos sao crculos mnimos, exceto quando a latitude for igual
a 90 , situac
ao em que um plano perpendicular ao Equador nao ira cortar a esfera, mas ser
a
8
Os valores da tabela 6 foram calculados no stio Length of degree calculator - National GeospatialIntelligence Agency http://www.nga.mil/MSISiteContent/StaticFiles/Calculators/degree.html
4.5
Coordenadas geogr
21
tangente à ela. Uma linha que une pontos de mesma direcao, ou seja, que cruza os meridianos
com o mesmo
angulo e chamada de linha de rumo ou loxodr
omia.
Convers
ao entre graus-minutos-segundos e graus decimais
As coordenadas esfericas podem ser especificadas no sistema sexagesimal, em graus, minutos e segundos (cada grau e dividido em 60 minutos e cada minuto em 60 segundos) mais uma
letra para indicar o quadrante como em 45 3157W, ou em graus decimais (valores positivos
para Norte e Leste e negativos para Sul e Oeste) como em -45,5325 .
A convers
ao entre valores no sistema sexagesimal (graus-minutos-segundos) para o sistema
decimal (graus decimais) e simples. Basta lembrar que em um grau ha 60 minutos e que em
cada minuto h
a 60 segundos (ou seja, em um grau ha 3600 segundos). Assim, dividimos os
minutos e os segundos por suas fracoes de grau e somamos os valores.
Como exemplo, vamos converter 45 3157W em graus decimais (lembre-se que coordenadas Oeste e Sul recebem valores negativos):
45 +
31
57
+
= 45 + 0.516666 + 0.015833 = 45.5325
60
3600
Para converter valores em graus decimais para a notacao graus-minutos-segundos, sao necessarios alguns passos a mais. Vamos converter -45,5325 de volta para o sistema sexagesimal:
1. A parte inteira do n
umero passa a ser o valor em graus.
45.5325 45
2. Multiplica-se a parte decimal por 60. A parte inteira do resultado passa a ser o valor
dos minutos.
0.5325 60 = 31.95 31
3. Multiplica-se a nova parte decimal por 60. O resultado sera o valor dos segundos (este
valor pode ter uma parte decimal ou ser arredondado).
0.95 60 = 57 57
4. Junta-se os resultados em um valor u
nico.
45.5325 = 45 3157W
4.6
Projec
oes cartogr
aficas
4.6.
22
Projec
oes cartogr
aficas
A projec
ao e necess
aria para transferir a forma tridimensional da Terra para o plano,
bidimensional (mapas). A fim de evitar, tanto quanto possvel, a distorcao (inevitavel),
existem diferentes modelos, que podem ser aplicados dependendo da situacao, area geografica,
e ate mesmo interesse poltico-econ
omico.
Podemos classificar as projec
oes de acordo com a superfcie de projecao e de acordo com
as propriedades que s
ao preservadas nos mapas (area, forma, distancias, etc).
De acordo com a superfcie de projecao (fig. 4.10):
Azimutal ou Plana: o plano de projecao e posicionado como tangente à superfcie da
Terra (do elips
oide). Pode-se imaginar uma fonte de luz no infinito, iluminando atraves
do globo, e projetando as sombras dos meridianos e paralelos no plano do mapa.
Cilndrica: o plano do mapa e tangente ao Equador, e enrolado no globo como um
cilindro. Os meridianos e paralelos sao projetados no plano e criam uma grade retangular.
C
onica: neste caso, o plano do mapa e enrolado como um cone, com uma linha de
tangencia com a Terra (paralelo padrao - standard parallel ). Freq
uentemente o vertice
do cone e alinhado com os dois polos, os meridianos se dispoes radialmente (com o
mesmo
angulo) e os paralelos formam crculos concentricos.
A superfcie de projec
ao pode ser orientada perpendicularmente ao eixo da Terra (Normal), paralela (Transversa) ou diagonalmente (Oblqua).
A superfcie de projec
ao pode tangenciar a Terra, caso em que ha um paralelo padrao
nas projec
oes c
onicas, ou pode cortar o globo (ser secante ao elipsoide), caso em que
haver
ao dois paralelos padr
ao (fig. 4.11).
De acordo com as propriedades preservadas:
Direca
o: azimutal
Forma (localmente): conformal ou ortomorfica
Area:
igual-
area ou equivalente ou autalica
Dist
ancia: equidistante
4.6
23
Projec
oes cartogr
aficas
Cilndricas
linha de
tangncia
linha de
tangncia
linha de
tangncia
Equatorial ou Normal:
eixo do cilindro paralelo
ao eixo da Terra
Transversa:
eixo do cilindro perpendicular
ao eixo da Terra
Cnicas
linha de
tangncia
Oblqua:
eixo do cilindro inclinado em
relao ao eixo da Terra
linha de tangncia
linha de
tangncia
Transversa:
eixo do cone perpendicular
ao eixo da Terra
Normal:
eixo do cone paralelo ao
eixo da Terra
Oblqua:
eixo do cone inclinado em
relao ao eixo da Terra
Planas (azimutais)
ponto de
tangncia
ponto de tangncia
ponto de
tangncia
Polar:
plano tangente ao plo
Equatorial:
plano tangente ao equador
Oblqua:
plano tangente em um
ponto qualquer
Figura 4.10. Classificaca

o de projec
oes cartogr
aficas de acordo com a superfcie de projeca
o (IBGE, 1998).
Cilndricas
Cnicas
linhas
secantes
linhas
secantes
Figura 4.11. Projeco

es secantes.
4.7
Projec
oes cilndricas
24
A seguir, veremos algumas das principais projecoes. Os crculos verdes sobre os mapas
sao indicatrizes de Tissot (cada crculo e uma indicatriz), e representam areas circulares com
o mesmo tamanho dispostos sobre o globo de forma a mostrar as deformacoes de area e forma
(ou ambas) que as projec
oes acarretam9 .
4.7.
Projec
oes cilndricas
Nas projec
oes cilndricas, os meridianos sao linhas verticais igualmente espacadas e os
paralelos s
ao linhas horizontais. Essas projecoes diferenciam-se pelo tipo de deformacao ao
longo da latitude ( direc
ao N-S).
4.7.1
Cilndrica conformal (Mercator)
Esta e um projec
ao conformal, portanto preserva-se a forma dos elementos do mapa.
Na projec
ao de Mercator, a escala N-S e igual à escala E-W, o que leva a uma deformacao
excessiva em altas latitudes.
Figura 4.12. Projeca

o cilndrica conforme de Mercator.
4.7.2
Cilndrica equidistante
Tambem chamada de equiretangular ou plate carre (literalmente, placa quadrada), nao e

nem conformal nem de igual-
area. Os meridianos e paralelos sao igualmente espacados e a
escala N-S e constante, portanto longitude e latitude sao mapeadas diretamente para x e y e
o mapa-mundi e um ret
angulo com proporcao 2:1.
9
As figuras das projec

oes foram retiradas de http://matplotlib.github.com/basemap/users/mapsetup.
html
4.8
25
Projec
oes c
onicas

o cilndrica equidistante.
4.8.
Projec
oes c
onicas
As projec
oes c
onicas s
ao mais usadas em areas com extensao maior E-W, fora do Equador.
4.8.1
C
onica conformal
Uma projec
ao secante com dois paralelos padrao (projecao conica conforme de Lambert).

o c
onica conforme de Lambert.
4.8
26
Projec
oes c
onicas
4.8.2
C
onica igual-
area
Projec
ao secante com dois paralelos padrao. Na projecao de Albers, a distorcao entre os
paralelos padr
ao e mnima.

o c
onica igual-
area de Albers.
4.8.3
C
onica equidistante
Nesta projec
ao, os paralelos s
ao igualmente espacados. Utilizada para areas proximas do
paralelo padr
ao.

o c
onica equidistante.
4.9
27
Projec
oes azimutais
4.8.4
Polic
onica
Indicada para
areas com grande extensao N-S. Nao e conformal ou igual-area. Baseada
em m
ultiplos cones, cada um com um paralelo padrao. Utilizada em mapas da America do
Sul e do Brasil.

o polic
onica.
4.9.
Projec
oes azimutais
4.9.1
Azimutal igual-
area
Nesta projec
ao, as direc
oes v
alidas a partir do ponto central. Toda linha reta que passe
pelo centro e um crculo m
aximo. Em geologia estrutural, e chamada de projecao de SchmidtLambert.

o azimutal igual-
area de Lambert.
4.9
28
Projec
oes azimutais
4.9.2
Azimutal equidistante
As dist
ancias e direc
oes s
ao v
alidas para qualquer ponto a partir do centro do mapa. Toda
linha reta que passe pelo centro e um crculo maximo.

o polar azimutal equidistante.
4.9.3
Azimutal ortogr
afica
Ponto de perspectiva no infinito. Mostra o planeta como um fotografia a partir da Lua.

o azimutal ortogr
afica.
4.10
29
Projec
oes compromisso
4.9.4
Azimutal gnom
onica
Ponto de perspectiva no centro do globo. Nao e conformal, nem igual-area nem equidistante. Crculos m
aximos s
ao linhas retas.

o azimutal gnom
onica.
4.10.
Projec
oes compromisso
4.10.1
Projec
ao de Robinson
Nao e conformal, nem igual-

area, foi usada pela National Geographic Society para seus
mapas-mundi.

o azimutal de Robinson.
4.11
4.11.
30
Coordenadas UTM
Coordenadas UTM
No sistema UTM (Universal Transverse Mercator ), utiliza-se a projecao cilndrica transversa de Mercator. A Terra e dividida em 60 zonas (fusos) com 6 de longitude cada, entre
84 de latitude norte e 80 de latitude sul, com excecao de algumas areas nas regioes proximas
aos polos (fig. 4.23). Cada zona UTM e mapeada por um cilindro transverso, com linha de
tangencia no meridiano central da zona.
As coordenadas s
ao indicadas por E (Este) e N (Norte) e tem sua origem no Meridiano
central de cada zona e no Equador. Para evitar valores negativos, o meridiano central recebe
o valor de 10.000 km e o Equador, 500 km (fig. 4.24).
Os meridianos centrais est
ao localizados em 3 , 9 , 15 , e assim por diante, para Leste
e Oeste. As zonas s
ao ainda divididas em faixas de sul para norte, indicadas por letras
(fig. 4.23). As faixas marcadas com letras de A ate M estao no hemisferio Sul, e de N ate Z
no hemisferio Norte.
Note que a notac
ao completa das coordenadas UTM inclui a letra da faixa, como em
23K 500.000E / 7.500.000N. Na figura 4.24 fica claro que o mesmo valor numerico de
coordenadas ocorre nos dois hemisferios10 , portanto indicar a letra da faixa e fundamental.
Um ponto importante a lembrar e que os programas de SIG usualmente diferenciam as zonas
UTM apenas por Norte ou Sul, como em 23N e 23S. Isso funciona para os programas, mas
devemos ficar atentos para n
ao cometer o erro de indica uma coordenada 23S para uma
localidade no hemisferio Sul, j
a que a zona 23S esta no hemisferio Norte.
31X 33X
35X 37X
X
W
32V
01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
U
U
U
T
23
K
K
J
01
26
60
02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59
E
E
E
D
C
A
D
C
B
Figura 4.23. Divis

ao das zonas UTM.
10
Na verdade, o mesmo valor numerico ocorre nas 60 zonas, nos dois hemisferios, ou seja, 120 vezes.
4.12
31
Carta Internacional do Mundo ao Milionesimo
Zona UTM 23 (48 W - 42 W)

Hemisfrio Norte
45 W
Y Z
Hemisfrio Sul
meridiano central
45 W
Equador
10.000.000 N
(metros)
X
W
V
L
23 W (500.000 E, 7.500.000 N)
hemisfrio Norte
K
J
23 K (500.000 E, 7.500.000 N)
hemisfrio Sul
U
T
5.000.000 N
(metros)
P
N
Equador
500.000 E
(metros)
Origem das coordenadas
para o hemisfrio norte
(0 E, 0 N)
500.000 E
(metros)
Origem das coordenadas
para o hemisfrio sul
(0 E, 0 N)
Figura 4.24. Origem das coordenadas nas zonas UTM.
4.12.
Carta Internacional do Mundo ao Milion

esimo
Em 1891, o alem
ao Albrecht Penck propos uma divisao do mundo em cartas na escala
1:1.000.000. Essa divis
ao e feita em fusos de 6 de longitude, numerados de 1 a 60 a partir
do antimeridiano de Greenwich para Leste. Cada fuso e dividido em faixas de 4 de latitude,
marcadas, a partir do Equador, pelas letras A, B, C, D, E, F,... mais a notacao do hemisferio
(N ou S). A figura 4.25 mostra a articulacao da Carta Internacional ao Milionesimo em
territorio brasileiro.
4.13
Sistema cartogr
afico nacional
32
Figura 4.25. Articulaca

o da Carta Internacional ao Milionesimo em territ
orio brasileiro.
4.13.
Sistema cartogr
afico nacional
4.14.
Padr
ao de exatid
ao cartogr
afica
O decreto federal n 89.817 de 20 de junho de 1984 (Anexo ??) estabelece as Normas

Tecnicas da Cartografia Nacional. Entre outras deliberacoes, determina-se que os produtos
cartograficos devem obedecer ao Padrao de Exatidao Cartografica (PEC), onde 90% dos
pontos no terreno n
ao podem ter erro de posicionamento (planimetrico e altimetrico) maior
do que o estabelecido pelo PEC, e que os mapas sao classificados de acordo com os seguinte
criterios:
Classe A
Erro planimetrico: 0,5 mm na escala da carta
Erro altimetrico: metade da equidistancia entre as curvas de nvel
4.15
Mapeamento sistem
atico
33
Classe B
Erro altimetrico: tres quintos da equidistancia entre as curvas de nvel
Classe C
Erro altimetrico: tres quartos da equidistancia entre as curvas de nvel
4.15.
Mapeamento sistem
atico
O mapeamento sistem
atico topografico do Brasil e de responsabilidade do IBGE e compreende as seguintes escalas: 1:1.000.000, 1:500.000, 1:250.000, 1:100.000, 1:50.000 e 1:25.000.
Mapeamentos em escala maior s
ao considerados cadastrais e de responsabilidade dos estados
e municpios. Suas escalas normais variam de 1:10.000 ate 1:2.000.
Os mapas com escala de 1:250.000 a 1:25.000 utilizam coordenadas UTM, enquanto que
os mapas em 1:1.000.000 e 1:500.000 estao em coordenadas geograficas, na projecao Conica
Conforme de Lambert.
A articulac
ao dos mapas segue uma nomenclatura que indica tanto a escala do mapa
quanto sua localizac
ao dentro da carta ao milionesimo. O ndice e baseado em divisoes sucessivas dos mapas, onde letras ou n
umeros representam cada etapa. O mapa inicial na escala
1:1.000.00 e dividido em 4 mapas 1:500.000 (V,X,Y,Z). Cada mapa 1:500.00 e dividido em
4 mapas 1:250.000 (A,B,C,D) e cada um desses em 6 mapas 1:100.000 (I,II,III,IV,V,VI).
Os mapas 1:100.000 s
ao divididos em 4 mapas 1:50.000 (1,2,3,4) e finalmente cada mapa
1:50.000 e dividido em 4 mapas 1:25.000 (NO,NE,SE,SO).
As figuras 4.26 e 4.27 ilustram o processo. Veja que a partir do ndice de nomenclatura
de um mapa e possvel determinar sua escala:
1 : 1.000.000
SF-23 Rio de Janeiro
1 : 500.000
SF-23-X
1 : 250.000
SF-23-X-B Ponte Nova
1 : 100.000
SF-23-X-B-I
1 : 50.000
SF-23-X-B-I-1 Catas Altas
1 : 25.000
SF-23-X-B-I-1-SO
Outro modelo utilizado para indexar os mapas no Brasil e o Mapa Indice (MI), que
atribui uma numerac
ao sequencial aos mapas de acordo com a escala. Na escala 1:50.000,
acrescenta-se ao c
odigo uma numeracao (1,2,3,4) e para 1:25.000 acrescenta-se um indicador
(NO,NE,SE,SO), tal como no esquema anterior:
4.15
34
Mapeamento sistem
atico
1 : 1.000.000
numerac
ao de 1 a 46
1 : 250.000
numerac
ao de 1 a 550
1 : 100.000
numerac
ao de 1 a 3036 (p.ex., MI 2574)
1 : 50.000
c
odigo MI + numeracao 1, 2, 3, 4
1 : 25.000
c
odigo MI 1:50.000 + indicador NO, NE, SE, SO
A tabela 7 apresenta ambos os sistemas de indexacao:

Escala
1 : 1.000.000
1 : 500.000
1 : 250.000
1 : 100.000
1 : 50.000
1 : 25.000
Indice de nomenclatura
Indice MI
SF-23 Rio de Janeiro

SF-23-X
SF-23-X-B Ponte Nova
SF-23-X-B-I
SF-23-X-B-I-1 Catas Altas
SF-23-X-B-I-1-SO
MI-39
MI-0466
MI-2574
MI-2574-1
MI-2574-1-SO
Tabela 7. Nomenclaturas de mapas no Sistema Cartogr

afico Nacional

o do sistema cartogr
afico nacional (vis
ao condensada).
4.15
Mapeamento sistem
atico

o do sistema cartogr
afico nacional (vis
ao expandida).
35
36
5 - Bem-vindo(a) ao GRASS-GIS!
5.
Bem-vindo(a) ao GRASS-GIS!
5.1.
Criando uma Location e um mapset
O primeiro passo ao trabalhar com GRASS e a criacao de uma Location e de ao menos

um mapset. Ao rodar o programa pela primeira vez, uma janela de terminal (tambem
chamada de console) como a da figura 5.1 e exibida. Essa mensagem somente e mostrada na
primeira vez em que o GRASS e iniciado e como o OSGeo-Live ja vem com alguns dados de
exemplo, n
ao veremos essa tela e sim a janela verde-clara da figura 5.2.
Figura 5.1. Tela de terminal com a mensagem de boas-vindas ao GRASS-GIS.
Na figura 5.2 temos o campo superior central (GIS Data Directory), com o caminho
para o diret
orio (pasta) GISBASE, onde sao armazenados os dados dos projetos. No campo
inferior da esquerda est
ao listadas as Locations; cada uma e um subdiretorio de GISBASE.
No campo inferior central, s
ao listados os mapsets; novamente, cada um e um subdiretorio
da sua Location.
No campo da direita podemos criar novas Locations e mapsets ou
renomear/deletar mapsets.
Para criar uma Location, clique no botao Location wizard.
A primeira etapa e
definir o caminho para o diret

orio GISBASE (GIS Data Directory, ja predefinido como
/home/usr/grassdata), e o nome da Location. O ttulo e opcional.
Na tela da figura 5.4 decidimos como criar a Location. Temos varias opcoes: selecionar
parametros de listagens; utilizar c
odigos pre-definidos pela European Petroleum Survey Group
(EPSG) para diversas combinac
oes de datums e projecoes cartograficas; ler os parametros a
partir de arquivos georreferenciados (como os GeoTIFFs), de arquivos WKT (Well-Known
Text) ou PRJ (arquivos com informacoes de projecao que acompanham dados vetoriais em
5.1
37
Figura 5.2. Janela inicial do GRASS.
Figura 5.3. Definica

o do diret
orio GISBASE e do nome da Location.
formato shapefile); definir projec

oes e datuns manualmente ou ainda criar sistemas arbitrarios
XY. Selecione a primeira opc
ao.
O primeiro par
ametro a ser definido e a projecao e o sistema de coordenadas. Na figura 5.5
temos uma listagem das projec
oes suportadas. Podemos fazer uma busca por palavras-chave
ou escolher da lista. O GRASS oferece suporte a diversas projecoes cartograficas e datums,
atraves da biblioteca PROJ.4
11
http://proj.maptools.org
11 .
Neste exemplo, usaremos Latitude/Longitude.
5.1
38
Figura 5.4. Defina como ser

a o metodo de criaca
o da Location.
Figura 5.5. Defina a projeca

o e o sistema de coordenadas.
Agora definimos se a projec

ao sera associada a um datum e um elipsoide ou se apenas a
um elipsoide (fig. 5.6). As opc
oes disponveis variam de acordo com a escolha passo anterior.
Selecione Datum with associated ellipsoid.
5.1
Figura 5.6. Escolha dos par

ametros da projeca
o.
A seguir voce ter

a que escolher o datum. Na tela da figura 5.7, selecione wgs84.
Figura 5.7. Escolha do datum WGS-84.
39
5.1
40
Agora uma tela com um resumo das suas escolhas e exibida (fig. 5.8). Se tudo estiver
certo, clique em Finish.
Figura 5.8. Confira se est

a tudo certo.
Definida a Location, voce ser

a perguntado se quer definir os limites e a resolucao espacial
da region padr
ao (fig. 5.9). Clique em Yes.
Figura 5.9. Quase l

a.
Na janela da figura 5.10, os campos ao redor do mapa indicam os limites geograficos da

region padr
ao. no nosso exemplo, usaremos North:
0, South:
-5, West:
-65, East:
-55. A resoluc
ao N-S e E-W ser
a de 0.0008333, que corresponde a 3 segundos de arco, que
e a resoluc
ao dos dados SRTM.
Esses valores s
ao em graus decimais. Quando ja estivermos trabalhando no GRASS, dentro
do terminal, os valores podem estar no formato <graus:minutos:segundos N/S/E/W> como
em <19:36:00S> ou em graus decimais (com valores negativos para os hemisferios Sul e
Oeste), como em <-19.6>, e a resolucao como em <0:00:03> ou <0.0008333>.
5.1
41
Figura 5.10. Defina os limites e a resoluca

o da region padr
ao.
Muito bem, agora a Location esta pronta. Clique em Create mapset (fig. 5.11) e defina
nosso mapset como manaus (fig. 5.12).
Figura 5.11. Location pronta, vamos ao mapset.
5.1
42
Figura 5.12. Defina o nome de nosso mapset.
Agora sim.
Criado o mapset, clique em Start GRASS! Temos agora tres jane-
las(fig. 5.13): uma de terminal (em baixo, à esquerda), onde podemos digitar comandos;
uma janela chamada de GIS Manager (em cima, à esquerda), onde teremos os arquivos
abertos para visualizac
ao e os menus de comandos; e uma janela de Map Display (em cima,
à direita), onde os dados s
ao exibidos. Diversos Map Displays podem ser abertos ao mesmo
tempo. Uma descric
ao mais detalhada das funcionalidades do GIS Manager e do Map Display
serao discutidas no pr
oximo captulo.
Figura 5.13. Interface gr

afica do GRASS-GIS 6.4
5.2
Apresentando o GIS Manager e o Map Display
5.2.
43
Neste captulo veremos brevemente as funcionalidades do GIS Manager e do Map Display.

No GIS Manager (ou, simplesmente, gism, fig. 5.14) encontramos todos os comandos separados em menus, e alguns dos mais utilizados estao agrupados em uma barra de ferramentas
(fig. 5.16). Abaixo da barra de ferramentas temos a area onde sao organizados os diversos
layers, vetoriais e raster, alem de composicoes coloridas (RBG ou IHS), graficos e elementos
cartograficos, como barras de escala e seta de norte. Note que o empilhamento vertical dos
layers no gism reflete a ordem em que serao exibidos no Map Display.
Figura 5.14. Janela do GIS Manager.
Na porc
ao inferior da janela do gism, temos tres abas: Map layers, Command console e
Search module. Na aba Map layers ficarao os layers carregados e que podem ser exibidos no
Map Display. A aba Command console e utilizada para exibir informacoes de sada (output) de
alguns comandos e tambem possui um console, ou terminal, onde podemos digitar comandos
(mais utilizado em ambiente Windows, onde a janela de terminal que vimos na figura 5.13
nao e muito utilizada). Na aba Search module, podemos buscar modulos (comandos) a partir
de parte do nome ou da descric
ao de sua funcionalidade (fig. 5.15).
Na figura 5.16, temos barra de ferramentas do gism, com a indicacao das funcionalidades
e na figura 5.17, a barra de ferramentas do Display Manager.
5.2
Figura 5.15. Busca de um m

odulo no gism.
Figura 5.16. Barra de ferramentas do GIS Manager.
Figura 5.17. Barra de ferramentas do Map Display.
44
45
6 - Dados Raster
6.
Dados Raster
Vamos comecar a utilizar o GRASS com dados tipo raster. Isso nos permitira aprender
como baixar dados de elevac

ao SRTM, como trabalhar com esses dados, e ainda, a trabalhar
com o conceito de region, algho exclusivo do GRASS e extremamente poderoso.
6.1.
Modelos de Elevac
ao SRTM
A miss
ao Shuttle Radar Topography Mission (SRTM Farr & Kobrick, 2000; van Zyl, 2001;
Rabus et al., 2003; Farr et al., 2007) foi realizada pela agencia espacial norte-americana (National Aeronautics and Space Administration NASA), a National Geospatial-Intellingence
Agency (NGA), o Departamento de Defesa dos Estados Unidos (DoD) e as agencias espaciais alema (Deustches Zentrum f
ur Luft- und Raumfahrt DLG) e italiana (Agenzia Spaziale
Italiana ASI) em fevereiro de 2000, para mapear o relevo da area continental da Terra com
interferometria de radar de abertura sintetica (InSAR), entre 60 de latitude norte e 54 de
latitude sul, o que corresponde `
a aproximadamente 80% das areas emersas do planeta.
Os princpios da InSAR foram desenvolvidos no incio dos anos 1980 (Zebker & Goldstein,
1986). Assim como a estereoscopia optica se baseia na paralaxe causada pela diferenca de
posicao do observador com relac
ao ao ponto de interesse (Lillesand & Kiefer, 1987), a altitude
de um ponto na superfcie pode ser determinada a partir da diferenca de fase entre reflexoes
de um mesmo sinal de radar captado por dois receptores distantes um do outro (Rabus et al.,
2003).
No caso da miss
ao SRTM, os canais principais de transmissao e recepcao para as bandas
C e X foram posicionados na
area de carga do onibus espacial Endeavour, enquanto que
as antenas secund
arias (somente recepcao) foram colocadas na extremidade de um mastro
retratil de 60 m de comprimento. Para uma descricao detalhada de interferometria de radar
de visada lateral, consultar, por exemplo, Bamler & Hartl (1998).
Os dados brutos foram processados pela NASA e, apesar de obtidos com resolucao espacial de 1 segundo de arco (aproximadamente 30 m no equador), estao disponveis com este
nvel de detalhe apenas para a
area dos Estados Unidos; para os outros pases, houve uma
reamostragem dos dados para 3 segundos de arco ( 90 m) de resolucao espacial. O produto
final possui precis
ao vertical global de 16 m e horizontal de 20 m (Rabus et al., 2003).
Para a America do Sul, a precis
ao vertical e de 6,2 m e a horizontal de 9,0 m (Rodriguez
et al., 2006).
O produto e georreferenciado ao datum WGS84 em coordenadas geograficas decimais e
pode ser obtido via Internet
12
13
12,13 .
http://seamless.usgs.gov
ftp://e0srp01u.ecs.nasa.gov/srtm
Os dados para o territorio brasileiro processados pela
6.1
Modelos de Elevac
ao SRTM
46
Embrapa (Miranda, 2005) est

ao disponveis segundo a articulacao das folhas topograficas de
escala 1:250.000
14 .
Deve-se levar em conta que o resultado das tecnicas empregadas na missao SRTM e a
geracao de Modelos Digitais de Elevac
ao (MDEs), pois os sinais de radar sao refletidos,
por exemplo, pelo dossel das
arvores em areas densamente florestadas e nao pelo terreno
subjacente. Os termos Modelo Digital de Terreno (MDT) e Modelo Numerico de Terreno
(MNT) devem ser reservados para casos onde o modelo e produzido a partir de valores de
altitude do nvel do terreno, obtidos, por exemplo, em mapas topograficos (curvas de nvel),
levantamentos por GPS ou por altimetria a laser (LiDAR).
Os modelos SRTM apresentam
areas sem valores de elevacao, chamadas de vazios (voids),
decorrentes de problemas na reflex
ao dos sinais de radar em corpos dagua ou em vertentes
com mergulho elevado (Farr & Kobrick, 2000; Farr et al., 2007). A interpolacao dos dados com
raio de busca pequeno (ou seja, utilizando poucos valores) significa que vazios maiores que a
janela de busca n
ao ser
ao completamente preenchidos. Por outro lado, utilizar um n
umero
maior de dados aumenta drasticamente o tempo de interpolacao e resulta em superfcies
muito suaves, sem o detalhamento desejado. No caso da presenca de vazios muito grandes,
uma alternativa vi
avel e preencher os vazios antes da reamostragem.
Uma caracterstica importante a ser observada e a metodologia empregada na reamostragem dos dados originais para 3 segundos de arco. Os dados distribudos pelo Servico Geologico
dos Estado Unidos (USGS)12 foram gerados a partir da media dos valores originais, em uma
janela de 3x3 pixels, o que resulta em um modelo suavizado com relacao ao original. Os dados
disponibilizados pela NASA13 foram gerados a partir do valor do pixel central da janela 3x3
e portanto possuem informac
oes dos dados originais (Fig. 6.1).
Figura 6.1. Metodos de reamostragem dos dados SRTM utilizados pela NASA e pelo USGS, para geraca
o
dos modelos com resoluca
o de 3 segundos de arco (adaptado de Steiner, 2007)
14
http://www.relevobr.cnpm.embrapa.br
6.2
47
Download e importac
ao dos dados
6.1.1
Modelos SRTM30 PLUS
Os modelos SRTM30 PLUS V3 (Becker & Sandwell, 2007) consistem em um modelo global
com resoluc
ao espacial de 0 0 30 (aproximadamente 1 km), criado a partir da fusao de dados
de diversos projetos (Fig. s6.2).
Nas areas continentais, o modelo utiliza principalmente os dados do USGS SRTM30,
uma reamostragem para 1km dos dados SRTM originais. Para as areas em altas latitudes,
onde nao h
a dados SRTM, foi utilizado o modelo GTOPO30 15 . A area oceanica e baseada no
modelo global de Smith & Sandwell (1997), com resolucao espacial de 0 01, entre as latitudes
72 . Modelos de mais alta resolucao foram incorporados a partir dos projetos LDEO Ridge
Multibeam Synthesis Project
16
e NGDC Coastal Relief Model
17 .
A batimetria do artico
deriva da International Bathymetric Chart of the Oceans (Jakobsson et al., 2000).
Figura 6.2. Modelo de elevaca

o global SRTM30 PLUS (Poppe et al., 2005)
6.2.
Download e importac
ao dos dados
Agora que j
a vimos um pouco do basico sobre o GRASS, podemos prosseguir com a
importacao da imagem SRTM, usando o comando r.in.gdal18 . Os dados utilizados neste
tutorial podem ser obtidos no site do Consortium for Spatial Information (CGIAR-CSI), que
disponibiliza dados pre-processados (com remocao de rudos etc) em arquivos para areas de
5 x5 . O link direto para os arquivos e:
http://goo.gl/YEI1e19
15
http://edcdaac.usgs.gov/gtopo30/gtopo30.asp
http://ocean-ridge.ldeo.columbia.edu/general/html/home.html
17
http://www.ngdc.noaa.gov/mgg/coastal/coastal.html
18
se voce tiver baixado a imagem SRTM como um arquivo com extens
ao .HGT (ou .HGT.ZIP), pode usar o
comando r.in.srtm para a importaca
o.
19
http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_24_13.zip
16
6.2
Download e importac
ao dos dados
48
http://goo.gl/eYdy420
Apos baixar os arquivos, continue com a importacao. Va em File Import raster
data Common import formats (Fig. 6.3). Note que à medida em que voce movimenta
o ponteiro do mouse sobre os tens dos menus, o nome do comando aparece na parte inferior
esquerda da janela do gism (nesse caso, r.in.gdal). Alem de importar o arquivo para
a Location ativa, pode-se ampliar os limites da region ativa ou mesmo criar uma nova
Location, de acordo com a imagem que esta sendo importada. Na figura 6.4 vemos as
opcoes do m
odulo. Selecione o arquivo com extensao .TIF (GeoTIFF) para importacao. Ser
a
necessario marcar a opc
ao Overrride projection (use locations projection) para que
a importac
ao seja realizada, uma vez que a imagem TIF nao possui informacao de elipsoide.
Importe os dois arquivos (srtm_24_13.tif e srtm_25_13/tif).
Figura 6.3. Selecione a importaca

o atraves da biblioteca GDAL
20
http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_25_13.zip
6.2
Download e importac
ao dos dados
Figura 6.4. Opco

es do comando r.in.gdal.
49
6.2
Download e importac
ao dos dados
50
Com as imagens importadas, vamos exibi-las no Map Display. Nas versoes mais atuais
do GRASS, v
arios comandos j
a se encarregam de incluir o layer resultante na pilha de layers
do gism. Caso isso n
ao ocorra, adicione um layer tipo raster clicando sobre o quinto cone
da barra de ferramentas superior do gism (fig. 5.16). Na aba textbfRequired, um menu tipo
drop-down mostra uma lista com os rasters disponveis (fig. 6.5). Selecione um dos arquivos
que importamos e clique em OK. Faca o mesmo para o segundo arquivo.
Figura 6.5. Insira um novo layer na pilha.
No Map Display, use as ferramentas de Zoom in, Zoom out e Pan (veja a figura 5.17)
para enquadrar as duas imagens na tela. Se estiver usando os mesmos arquivos de exemplo,
voce dever
a ver o mesmo que a figura 6.6. As cores representam a elevacao do modelo, e
rapidamente pode-se notar algumas areas onde a variacao de cores nao e suave como seria de
se esperar, mas e brusca.
Figura 6.6. As imagens SRTM, exibidas no Map Display
6.3
51
Raster patching
6.3.
Raster patching
Nosso primeiro passa ser

a juntar os dois arquivos SRTM em um raster u
nico. Para isso
iremos usar o comando r.patch, que voce encontra em Raster Overlay rasters
Patch raster maps (fig. 6.7).
Figura 6.7. Localize o comando r.patch.
Vamos aproveitar e ver algumas das opcoes que comumente temos nos comandos do
GRASS. Quando a janela de di
alogo de um comando e aberta, vemos algumas abas na parte
superior (a quantidade de abas varia de acordo com o comando). Na figura 6.8 temos a primeira aba, Required, que como o nome indica, sao os parametros obrigatorios para que o
modulo possa ser executado. Na aba Optional (fig. 6.8) temos algumas opcoes, como deixar
que o GRASS sobrescreva arquivos j
a existentes (o padrao e nao sobrescrever). A aba Command output (fig. 6.9) ir
a mostrar as informacoes de sada do modulo e a aba Manual
(fig. 6.9) traz explicac
oes sobre todas as opcoes e funcionamento do comando.
6.3
52
Raster patching
Figura 6.8. Abas Required e Optional do comando r.patch.
Figura 6.9. Abas Command output e Manual do comando r.patch.
Na aba Required, o primeiro campo que temos a preencher e com os mapas que queremos
juntar. Selecione os dois rasters que importamos (um por vez). No segundo campo indique o
nome do arquivo de sada, no nosso exemplo, usamos srtm_patch.
Aqui vale uma dica: computadores sao burros, e ter que trabalhar em plataformas diferentes (Windows, Mac, Linux) s
o complica as coisas. Entao tente minimizar as chances de
problemas futuros e use nomes inteligentes em seus arquivos. Nao use espacos entre as palavras, use um _ (underscore/underline) ou use o que programadores chama de CamelCase,
onde cada palavra comeca com letra mai
uscula, sem espaco entre as palavras. Como os no-
6.3
53
Raster patching
mes de arquivos ser

ao usados internamente pelos programas de SIG (ou seja, serao variaveis
em termos de programac
ao), os nomes nao devem (nao podem) comecar por n
umeros, mas
podem ter n
umeros no meio ou no final do nome. Finalmente, tente usar nomes que indiquem
que arquivo e esse e o que foi feito com ele, afinal nossa memoria nao e tao confiavel a longo
prazo. Um raster que representa um modelo de terreno pode muito bem ser chamado de
DEM (Digital Elevation Model ) ou MDT (Modelo Digital de Terreno) como em DEM_manaus ou
manausDEM. Se esse raster passar por uma filtragem, por exemplo, o arquivo de sada pode
ser manausDEM_filter_3x3 (aqui 3x3 representa o tamanho do filtro aplicado). Bem melhor
do que nomear os arquivos de mapa1, mapa2, mapa_teste, mapa_agora_vai...
Apos rodar o comando, o novo raster srtm_patch ja deve ser carregado e exibido no Map
display (fig. 6.10). Veja que agora a transicao de cores e suave, pois antes uma tabela de cores
era aplicada a cada layer, e agora temos apenas um layer.
Figura 6.10. Mapa raster resultante do comando r.patch.
Va agora ao gism e clique sobre o cone que esta logo à direita do nome do layer do
arquivo que criamos (clicar com o botao direito do mouse sobre o nome do arquivo tambem
funciona) (fig. 6.11). Veja que aqui temos varias opcao à disposicao, como remover o arquivo
da pilha (sem deletar o arquivo do computador), renomear o nome de exibicao (o arquivo
nao e renomeado), mudar a opacidade (deixar o mapa transparente), mudar as propriedades
graficas do que e exibido (cores, espessura de linhas etc) e varias outras. No momento, clique
na u
ltima opc
ao, Metadata. Os metadados sao dados sobre os dados, e indicam propriedades
importantes, como limites geogr
aficos, hora de criacao do arquivo e historico de acoes (como
indicar que o arquivo foi criado a partir da juncao de dois outros arquivos). Ao selecionar
essa opcao, o gism deve ativar a aba Command console, com a sada do comando r.info
(fig. 6.12).
6.3
54
Raster patching
Figura 6.11. Localize a opca

o de exibir os metadados do arquivo que criamos.
Figura 6.12. Metadados do arquivo criado com r.patch.
6.3
Raster patching
Esta e a informac
ao que foi exibida na janela do gism:
+----------------------------------------------------------------------------+
| Layer:
srtm_patch@manaus
Date: Sun Sep 9 07:56:40 2012
|
| Mapset:
manaus
Login of Creator: user
|
| Location: cursoUFAM
|
| DataBase: /home/user/grassdata
|
| Title:
( srtm_patch )
|
| Timestamp: none
|
|----------------------------------------------------------------------------|
|
|
|
Type of Map: raster
Number of Categories: 543
|
|
Data Type:
CELL
|
|
Rows:
6001
|
|
Columns:
12002
|
|
Total Cells: 72024002
|
|
Projection: Latitude-Longitude
|
|
N:
0
S: 5:00:02.27988S
Res: 0:00:02.99988
|
|
E: 54:59:58.07976W
W: 65:00:02.63952W
Res: 0:00:02.99988
|
|
Range of data:
min = -38 max = 543
|
|
|
|
Data Description:
|
|
generated by r.patch
|
|
|
|
Comments:
|
|
r.patch input="srtm_25_13,srtm_24_13" output="srtm_patch"
|
|
|
+----------------------------------------------------------------------------+
55
6.4
Trabalhando com regions
6.4.
56
O conceito de region e exclusivo do GRASS, e muito poderoso. A region define quais

os limites da
area geogr
afica sendo utilizada para analises e qual a resolucao espacial que ser
a
adotada.
Ao alterar a resoluc
ao da region, o GRASS faz uma reamostragem dos dados pelo metodo
do vizinho mais pr
oximo, ou seja, o valor a ser usado para um pixel sera o valor do pixel original
mais proximo de seu centro.
Alterar os limites da region permite, por exemplo, executar analises em sub-areas de
imagens que s
ao muito grandes e cujo processamento seria muito demorado. Assim podemos
testar nossas ideias rapidamente, antes de aplica-las para uma area maior.
6.4.1
Relevo sombreado
Vamos ver isso na pr

atica. Vamos fazer uma imagem de relevo sombreado (shaded relief )
dos nossos dados SRTM. V
a em Raster Terrain analysis Shaded relief (fig. 6.13).
Figura 6.13. Localize o comando r.shaded.relief.
Na aba Required, entre com o mapa srtm_patch. As opcoes importantes estao na

aba Optional (fig. 6.14). Este comando vai simular uma iluminacao solar sobre o modelo
de elevacao. Para isso precisamos informar a altitude do sol com relacao ao horizonte (no
exemplo, 25 ) e tambem o seu azimute (45 ). Adicionalmente, podemos utilizar um exagero
vertical (5), muito u
til em
areas sem desnvel acentuado. Por fim, esse comando espera que
as unidades dos dados na horizontal e na vertical sejam as mesmas (isso vale para outros
SIGs), mas no nosso caso elas n
ao s
ao. As unidades horizontais estao em graus e as verticais
em metros. Para corrigir esse problema, devemos selecionar meters na u
ltima opcao, que faz
o ajuste entre as escalas horizontal e vertical.
6.4
57
Finalmente, de um nome ao arquivo de sada. Para nao perder o costume de usar nomes
explicativos, uma boa opc
ao e srtm_patch_shade_045_25_z5. Isso significa que o mapa base
foi o srtm_patch, que foi aplicado um sombreamento shade, com azimute de 045 graus,
sempre bom utilizar tres
inclinacao de 25 graus e exagero vertical (em z ) de 5 vezes. E
algarismos para representar o azimute, para evitar confusao com a inclinacao. O resultado
esta na figura (6.15).
Figura 6.14. Opco

es do comando r.shaded.relief.
Figura 6.15. Resultado do comando r.shaded.relief, com resoluca

o original de 0:0:03 graus.
6.4
58
Nao da para ver muita coisa. Isso porque nossa area e muito grande, e a resolucao original
dos dados (3 arcos de segundo, ou aproximadamente 90 m) e muito detalhada para observar
a area como um todo. Para poder ver as formas de relevo maiores, precisamos diminuir a
resolucao espacial.
Va na janela de terminal e digite g.region -p e pressione Enter:
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
3 (Latitude-Longitude)
0
wgs84
wgs84
0
5:00:02.27988S
65:00:02.63952W
54:59:58.07976W
0:00:02.99988
0:00:02.99988
6001
12002
72024002
Aqui o comando g.region imprime (com a opcao -p) as configuracoes da region ativa:
os limites geogr
aficos, a resoluc
ao espacial e a quantidade de linhas e colunas.
Vamos alterar a resoluc
ao espacial para cerca de 1 km (30 segundos de arco). Para isso,
informamos o par
ametro res=xxxx ao comando g.region: Ao trabalhar com projetos em
Latitude/Longitude, podemos usar valores no formato <graus:minutos:segundos N/S/E/W>
como em <19:36:00S> ou em graus decimais (com valores negativos para os hemisferios Sul
e Oeste), como em <-19.6>, e a resolucao como em <0:00:03> ou <0.0008333>.
GRASS 6.4.1
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
(cursoUFAM):~ > g.region -p res=0:0:30

0
wgs84
wgs84
0
5:00:02.27988S
65:00:02.63952W
54:59:58.07976W
0:00:30.0038
0:00:30.0038
600
1200
720000
6.4
59
Agora refaca a imagem de relevo sombreado, usando os mesmos parametros de iluminacao

e exagero vertical. Na figura 6.16 podemos ver bem melhor o relevo da regiao.
Figura 6.16. Resultado do comando r.shaded.relief, com resoluca

o de 0:0:30 graus.
Mas ainda temos uma

area muito grande para trabalhar. Vamos escolher uma sub-area interativamente. No Map Display, usa a ferramenta de Zoom e selecione uma area (fig. 6.17).
Figura 6.17. Selecione uma sub-

area da imagem original.
Depois de escolher a sub-

area, v
a em Zoom options Set computational region
from display. Volte ao terminal e veja o resultado:
6.4
GRASS 6.4.1
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
60
(cursoUFAM):~ > g.region -p res=0:0:03

0
wgs84
wgs84
1:39:00.75236S
3:51:01.755508S
61:40:58.12232W
57:16:26.112226W
0:00:03.00038
0:00:02.999813
2640
5291
13968240
Como o tamanho das celulas esta um pouco diferente nas direcoes N-S e E-W, vamos
forcar a region a se alinhar (com a opcao -a) ao tamanho das celulas, usando g.region
-pa:
GRASS 6.4.1
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
(cursoUFAM):~ > g.region -pa res=0:0:03

0
wgs84
wgs84
1:39S
3:51:03S
61:41W
57:16:24W
0:00:03
0:00:03
2641
5292
13976172
Bem melhor. Assim voce pode selecionar uma area interativamente ate que fique do seu
agrado. E pode ainda salvar essa region para acesso rapido, com g.region -pa res=0:0:03
save=manaus. Para ativar a region ja salva, use g.region region=manaus.
Podemos ver como os par
ametros da region afetam nossos dados se configurarmos o
Map Display para exibir os mapa de acordo com as configuracoes da region ativa (por
padrao, os mapas raster s
ao exibidos com suas propriedades originais). Para isso, localize na
barra inferior do Map Display um menu que originalmente indicar Coordinates. Esse menu
possui diversas opc
ao do Map Display. Quando a opcao Coordinates e a selecionada, a parte
esquerda da barra inferior do Map Display mostra as coordenadas do ponteiro do mouse.
Selecione Display mode e veja que na parte esquerda da barra agora ha uma check-box que
indica Constrain display resolution to computational settings. Se essa caixa estiver
6.4
61
marcada, a resoluc
ao dos mapas no display ira seguir as definicoes da region (os dados nao sao
alterados). Faca o teste, marque a caixa, altere a resolucao da region e force o re-desenho dos
dados usando o segundo cone da barra de ferramentas do Map Display (Re-render map).
As figuras a seguir mostram esse comportamento.
Figura 6.18. Display com resoluca

o de 30 segundos (g.region res=0:0:30).

o de 05 minutos (g.region res=0:05).
6.4
62


Antes de passarmos para os dados vetoriais, vamos ver um comando do terminal muito u
til,
o man. Esse comando mostra o manual de outros comandos. Se voce digitar, por exemplo, man
g.region, ver
a uma tela com todas as opcoes do g.region. Navegue por essa tela usando as
setas para cima e para baixo do teclado. Quando quiser sair, pressione q.
Outra maneira de usar os comandos do GRASS e digitar apenas o comando no terminal
e pressionar Enter. Isso abrir
a a janela do comando digitado.
7 - Dados vetoriais
7.
63
Dados vetoriais
Neste captulo veremos um pouco mais sobre dados em formato vetorial. Vamos aproveitar
a oportunidade para comecar a trabalhar com o QGIS.
7.1.
Arquivos ESRI shape
Dentre os diversos formatos de mapas vetoriais existentes, provavelmente o mais utilizado atualmente seja o shapefile, criado pela empresa ESRI21 , desenvolvedora do pacote SIG
c
proprietario ArcGIS
. Com uma implementacao simples e especificacao tecnica aberta22 , e
muito difcil encontrar, hoje em dia, programas SIG que nao suportem o formato SHP.
Um shapefile consiste de um arquivo principal, um arquivo de ndice e uma tabela em
formato dBASE. O arquivo principal contem os registros de todos os elementos graficos (shapes
pontos, linhas, polgonos) com uma lista dos vertices de cada elemento. No arquivo de ndice,
cada registro determina a posic
ao do registro correspondente no arquivo principal. A tabela
dBASE contem os atributos das entidades geometricas. Note que devido à limitacoes do
formato dBASE IV, os nomes das colunas na tabela de atributos nao podem exceder onze
caracteres.
Adicionalmente, pode-se ter um arquivo que especifique a projecao cartografica e ainda
arquivos auxiliares. Todos os arquivos devem ter o mesmo nome, sendo diferenciados pela
extensao. Por exemplo:
estradas.shp arquivo principal, armazena os elementos graficos
estradas.shx arquivo de ndice dos elementos graficos
estradas.dbf tabela de atributos em formato dBASE IV
estradas.prj arquivo de projec
ao cartografica
estradas.sbn arquivo auxiliar - ndice espacial dos elementos graficos
estradas.sbx arquivo auxiliar - ndice espacial dos elementos graficos
7.2.
Dados vetoriais no QGIS
Um dado muito interessante e o mapa-ndice que mostra a situacao atual do mapeamento

sistematico topogr
afico do Brasil. Sao diversos layers vetoriais, que trazem quais mapas
topograficos est
ao disponveis de acordo com a escala. Baixe o arquivo em http://goo.gl/
X4Y3I23 . A nota explicativa est
a em http://goo.gl/WKDuK24 .
21
http://www.esri.com
http://www.esri.com/library/whitepapers/pdfs/shapefile.pdf
23
ftp://geoftp.ibge.gov.br/mapeamento_sistematico/mapa_indice_digital_4ed/produto_mapa_
indice_digital/shape.zip
24
ftp://geoftp.ibge.gov.br/mapeamento_sistematico/mapa_indice_digital_4ed/leia_me_4ed.pdf
22
7.2
64
Apos baixar os dados, descompacte o arquivo shape.zip e inicie o QGIS. A interface do

QGIS possui uma
area `
a esquerda para a pilha de layers carregados e uma area à direita para
a exibicao dos mapas.
Tanto no QGIS como na maioria dos SIG (ArcGIS etc) os dados (shapefiles, geotiffs etc)
nao sao abertos usando a t
ao comum combinacao File Open. Esse caminho e usado para
abrir arquivos de projeto, que s
ao usados quando se deseja salvar todo o esquema de trabalho,
incluindo os diversos layers carregados, suas propriedades graficas (cores, simbologia) e demais
configurac
oes. Os dados vetoriais ou raster sao adicionados no projeto atual. Isso e feito
clicando no cone com um smbolo de +(Add vector layer - fig. 7.1). O QGIS tem botoes
diferentes para cada tipo de dado (vetor, raster etc), enquanto que o ArcGIS, por exemplo,
tem apenas um bot
ao para adicionar dados, independente do tipo.
Figura 7.1. Interface gr

afica do QGIS com bot
ao para adicionar dados vetoriais selecionado.
No arquivo shape.zip que baixamos, existem alguns dados de referencia alem das cartas
topograficas disponveis por escala. Sao dados em escala 1:1.000.000:
Malha Municipal (LM_MUNICIPIO_2007)
Hidrografia (HD_MID_CURSO_DAGUA e HD_MID_MASSA_DAGUA)
Rodovias (ST_MID_RODOVIA)
Capitais (LC_MID_CAPITAL)
Unidades federativas (LM_UF)
7.2
65
Adicione, por exemplo, o layer LM_UF (fig. 7.2) e voce tera um mapa do Brasil dividido
por estado (cada estado e um polgono). Adicione depois o layer VETOR_250000, que mostra
quais mapas topogr
aficos em escala 1:250.000 estao disponveis em formato vetorial (fig. 7.3).
As cores que voce ver
a certamente serao diferentes das exibidas nas figuras, pois elas sao
atribudas aleatoriamente.
Figura 7.2. Selecione o layer LM_UF.
Figura 7.3. QGIS com layers LM_UF e VETOR_250000 carregados.
Aproveite para explorar a interface do QGIS. O cone destacado na figura 7.3 e o da

ferramenta para identificar feic
oes e mostrar os atributos dessas feicoes que estao na tabela
de atributos. O layer que ser
a analizado e o que estiver selecionado na pilha de layers. A
selecao e feita ao clicar sobre o nome do layer, e uma faiza azul sobre o nome mostra qual
layer esta selecionado. No nosso exemplo, selecionamos o layer VETOR_250000, e ao clicarmos
7.2
66
sobre um dos polgonos vemos uma janela como a da figura 7.4. Todas essas informacoes
estao armazenadas na tabela de atributos (o arquivo .dbf) do shapefile.
Figura 7.4. Informaco

es sobre um dos polgonos do layer VETOR_250000.
Mas fica um pouco difcil visualizar os dois layers um sobre o outro. vamos entao alterar
as propriedades de exibic
ao do layer VETOR_250000. De um duplo-clique sobre o nome do
layer (ou clique com o bot
ao direito e selecione a opcao Properties). A janela de dialogo
das propriedades e mostrada na figura 7.5. Aqui podemos, por exmplo, alterar a cor do
preenchimento dos polgonos, a cor e espessura das linhas entre varias outras opcoes. Altere
a transparencia do preenchimento dos polgonos e clique em Apply. Assim voce pode testar
varias configurac
oes ate encontrar a que mais lhe agrada. Quando estiver satisfeito, clique
em OK.
Na figura 7.6 vemos que alguns dos polgonos estao mais escuros que outros. Qual sera a
causa? Como os preenchimentos est
ao com transparencia, uma cor mais escura pode indicar
que ha dois ou mais polgonos sobrepostos, o que nesse caso indicaria um erro nos dados, pois
7.2
67
Figura 7.5. Janela de di

alogo das propriedades do layer.
ha duplicac
ao de informac
ao. Para confirmar a suspeita, use a ferramenta de informacao e
clique sobre um desses polgonos mais escuros (pode ser necessario usar a ferramenta Zoom
antes). Na janela de resultados (7.7) vemos que realmente temos dois elementos no mesmo
lugar.
Muito bem. Agora que j
a estamos acostumados com o QGIS, vamos baixar alguns dados
para trabalhar. O IBGE disponibiliza arquivos shapefile do que e chamado de Base vetorial
contnua em escala 1:250.000. Ainda nao ha um recobrimento completo do Brasil (fig. 7.8)25 ,
mas na regi
ao Amaz
onica temos os dados necessarios.
Esses dados est
ao separados por temas (estradas, rios etc) e e preciso baixar os arquivos
individualmente. Para nossos exerccos vamos baixar os arquivos de Hidrografia e Curvas de
Nvel das
areas 1-2008 e 4-2009.
25
ftp://geoftp.ibge.gov.br/mapeamento_sistematico/base_vetorial_continua_escala_250mil/
mapabrasil250.pdf
7.2
68
Figura 7.6. Layer VETOR_250000 com transparencia no preenchimento.
Figura 7.7. Janela de informaco

es com dois resultados.
Area
1-2008
contorno_hd_permanente_final.zip (4.946 kB)
http://goo.gl/Zbihm26
rio_permanente_ou_aluviao_final.zip (34.461 kB)
http://goo.gl/ETXlF27
curvas_de_nivel_final.zip (23.762 kB)
http://goo.gl/OAcok28
26
ftp://geoftp.ibge.gov.br/mapeamento_sistematico/base_vetorial_continua_escala_250mil/area_
1_2008/contorno_hd_permanente_final.zip
27
1_2008/rio_permanente_ou_aluviao_final.zip
28
1_2008/curvas_de_nivel_final.zip
7.2
69
Figura 7.8. Recobrimento da base vetorial contnua em escala 1:250.000.
Area
4-2009
hid_limite_massa_dagua_l.zip (15.506 kB)
http://goo.gl/jYcYq29
hid_massa_dagua_a.zip (2.989 kB)
http://goo.gl/EOvTN30
rel_curva_nivel_l.zip (20.0938 kB)
http://goo.gl/iiI1831
29
4_2009/hid_limite_massa_dagua_l.zip
30
4_2009/hid_massa_dagua_a.zip
31
4_2009/rel_curva_nivel_l.zip
7.2
70
Apos descompactar os arquivos, insira-os no QGIS. Voce vera que os dados das areas
1-2008 e 4-2009 n
ao encaixam como seria de se esperar. Ao usar a ferramenta Zoom Full,
vemos que deve haver um problema, pois os dados da area 4-2009 sao exibidos bem pequenos
à norte da
area 1-2009 (fig. 7.9).
Esse tipo de problema, infelizmente, nao e tao raro de acontecer. Normalmente acontece
quando o sistema de coordenadas ou a projecao cartografica dos dados nao e compatvel.
Para verificar essas informac
oes, abra os arquivos de projecao (.prj) dos shapefiles. Use um
editor de textos simples, como o Mousepad, que vem instalado no OSGeo-Live 5.5 e esta no
menu Applications Accessories Mousepad.
Figura 7.9. Dados das a

reas 1-2008 e 4-2009 exibidos no QGIS. Os dados da a
rea 4-2009 est
ao no alto,
pr
oximos da seta do mouse.
Este e o conte
udo do arquivo hid_limite_massa_dagua_l.prj (area 4-2009):
GEOGCS["GCS_SIRGAS_2000",DATUM["D_SIRGAS_2000",
SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
E este e o conte
udo do arquivo contorno_hd_permanente_final.prj (area 1-2008):
PROJCS["South_America_Albers_Equal_Area_Conic",GEOGCS["GCS_SIRGAS_2000",
DATUM["D_SIRGAS_2000",SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],
PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",-59.0],PARAMETER["Standard_Parallel_1",1.0],
PARAMETER["Standard_Parallel_2",-14.0],
PARAMETER["Latitude_Of_Origin",5.0],UNIT["Kilometer",1000.0]]
7.2
71
Nosso problema e que os dados da area 1-2008 estao em uma projecao pouco usual (Albers
conica), com coordenadas em quil
ometros, enquanto que os dados da area 4-2009 estao em
coordenadas geogr
aficas sem uma projecao especfica. Para poder trabalhar com esses dados
no GRASS, temos que converte-los para uma mesma projecao e para isso vamos utilizar o
comando ogr2ogr, que faz parte da biblioteca GDAL/OGR.
Esse comando pode ser utilizado para converter dados vetoriais em diferentes formatos e
alterar ou atribuir projec
oes cartograficas. Na janela de terminal, sera preciso navegar ate
o diretorio onde est
ao os arquivos. Isso e feito com o comando cd change directory. Para
descobrir em qual diret
orio voce est
a, use pwd print work directory. Para listar os elementos
do diretorio, use ls.
Por exemplo:
GRASS 6.4.1 (cursoUFAM):~ > pwd
/home/user
GRASS 6.4.1 (cursoUFAM):~ > ls
data
gisvm
launchassist.sh saga
Desktop
grassdata ossim
sharedOSX
Downloads gvSIG
qgis-examples
welcome_message.txt
GRASS 6.4.1 (cursoUFAM):~ > cd sharedOSX/
GRASS 6.4.1 (cursoUFAM):~/sharedOSX > ls
contorno_hd_permanente_final hid_massa_dagua_a
shape
curvas_de_nivel_final
hid_massa_dagua_a.zip
srtm_24_13
hid_limite_massa_dagua_l
hid_trecho_drenagem_l.zip
srtm_25_13
hid_limite_massa_dagua_l.zip hid_trecho_massa_dagua_a.zip
GRASS 6.4.1 (cursoUFAM):~/sharedOSX > cd contorno_hd_permanente_final/
GRASS 6.4.1 (cursoUFAM):~/sharedOSX/contorno_hd_permanente_final > ls
contorno_hd_permanente_FINAL.dbf contorno_hd_permanente_FINAL.shp
contorno_hd_permanente_FINAL.prj contorno_hd_permanente_FINAL.shx
GRASS 6.4.1 (cursoUFAM):~/sharedOSX/contorno_hd_permanente_final >
Agora que j
a estamos no mesmo diretorio dos dados, vamos à conversao. No nosso exemplo, o par
ametro -t_srs EPSG:4674 indica para transformar o sistema de referencia atual
para o sistema EPSG:4674. Os c
odigos EPSG sao pre-definidos pela European Petroleum
Survey Group para diversas combinacoes de datums e projecoes cartograficas. Esses codigos
podem ser consultados em http://spatialreference.org/ref/epsg/. Como parametros
finais ao comando ogr2ogr, informamos primeiro o nome do arquivo de sada e depois o
nome do arquivo a ser processado.
ogr2ogr -t_srs EPSG:4674 contorno_hd_permanente_LL.shp contorno_hd_permanente_FINAL.shp
7.2
72
Experimente agora adicionar o arquivo contorno_hd_permanente_LL.shp no QGIS e veja

o resultado (fig. 7.10).
Figura 7.10. Dados das

areas 1-2008 e 4-2009 exibidos no QGIS com a projeca
o corrigida.
8 - Importac
8.
73
Importac
Agora que j
a temos os shapefiles com a projecao cartografica correta, vamos importa-los
no GRASS. V
a em File Import vector data Common import formats (figura
8.1). Veja que a janela que e aberta (Import vector data - fig. 8.2), e apenas uma interface
simplificada para o comando v.in.ogr, que ainda pode ser acessado no botao Command
dialog, na parte inferior da janela.
Figura 8.1. Selecione a importac

ao atraves da biblioteca OGR
No dialogo de importac
ao, escolhemos o tipo de arquivo de uma lista com os formatos
suportados pela biblioteca OGR (o ESRI Shapefile e o padrao) e localizamos o arquivo com
extensao .shp. Caso o shapefile n
ao possua o arquivo de projecao cartografica (.prj), mas
sabemos que os dados est
ao no mesmo sistema de coordenadas da nossa Location, devemos
marcar a opc
ao Override dataset projection (use locations projection).
Na figura 8.3 vemos a janela de dialogo completa do modulo v.in.ogr. Alem de importar
o arquivo para a Location ativa, pode-se importar apenas as entidades localizadas dentro de
determinadas coordenadas (subset), ou entidades selecionadas pelos seus atributos, por meio
de express
oes SQL. Note que esta janela de dialogo possui diversas abas de opcoes, apesar
disso ficar um pouco escondido e marcado apenas pela setinha preta na parte superior da
janela, ao lado das abas. A figura 8.3 mostra a janela com todas as opcoes.
Na primeira aba, o campo (OGR datasource name) corresponde ao diret
orio (pasta)
onde estao os arquivos. Aqui tambem indicamos o nome do mapa que sera criado. Na segunda
aba (Selection), o primeiro campo (OGR layer name) e onde escolhemos qual arquivo ser
a
8 - Importac
74
Figura 8.2. Importaca

o de shapefile pela interface simplificada do comando v.in.ogr.
importado (cuidado, se nenhum nome for marcado aqui, todos os arquivos shapefile existentes
no diretorio ser
ao importados!). O nome do shapefile deve ser informado sem a extensao
shp. Na aba Attributes, por exemplo, temos a opcao Do not create attribute table, mais
utilizada para importar arquivos com dados de elevacao 3D, que ocupariam muito espaco
caso a tabela de atributos fosse criada. A opcao Override dataset projection, esta na aba
Optional.
Importe os shapefiles usando o dialogo simplificado. A figura 8.4 mostra os dados de
hidrografia sobre o modelo de elevacao SRTM.
Voce pode notar que os dados do shapefile contorno_hd_permanente_LL aparecem um
pouco diferentes no Map Display. Use a ferramenta de Zoom para examinar em detalhe
(fig: 8.5). O que vemos e que as
areas estao preenchidas em cinza e ha diversos smbolos x
sobre o mapa. Isso aconteceu porque o shapefile que importamos contem polgonos. Enquanto
que em um shapefile n
ao e possvel ter mais de um tipo de geometria (ponto, linha, polgono)
por arquivo, o formato vetorial do GRASS nao tem essa limitacao. Quando um arquivo de
polgonos e importado, temos as linhas que representam os limites dos polgonos (boundaries),
as
areas dos polgonos e os centr
oides de cada area (algo como o centro de massa do polgono).
Se quisermos, podemos pedir ao GRASS que exiba somente as linhas ou os centroides ou
combinacoes das diversas geometrias.
8 - Importac
Figura 8.3. Opco

es do comando v.in.ogr.
Figura 8.4. Vetores de hidrografia importados para o GRASS.
75
8 - Importac
76
Figura 8.5. Detalhe do layer contorno_hd_permanente_LL, com a

reas preenchidas e centr
oides.
Para isso, abra o di

alogo de Propriedades do layer contorno_hd_permanente_LL clicando
no botao à direta do nome do layer no gism (ou um clique com o botao direito do mouse sobre o
nome) e escolha a opc
ao Properties da lista. Na janela que e aberta e a do comando d.vect,
responsavel pela exibic
ao de dados vetoriais (para dados raster, o comando e d.rast). Na
aba Selection, na opc
ao Feature type, deixe marcada apenas a caixa de boundary (fig: 8.6).
Aqui voce pode clicar em Apply para ver o resultado no Map Display sem ter que fechar a
janela do d.vect, e assim mudar os parametros ate chegar a um resultado satisfatorio.
Figura 8.6. Opco

es do comando d.vect.
77
8 - Importac
Isso resolve nosso problema mas apenas em parte. Se quisermos realmente usar esses
dados como vetores de linhas e n
ao de polgonos, temos que forcar uma conversao durante
a importac
ao. Para isso e necess
ario abrir o dialogo completo do comando v.in.ogr e na
aba Selection, marcar a opc
ao import area boundaries as lines (fig: 8.7). No exemplo,
importamos novamente o arquivo contorno_hd_permanente_LL mas desta vez com o nome
de contorno_hd_permanente_LL_lines.
Figura 8.7. Convers

ao de tipos de geometrias vetoriais durante a importaca
o.
Podemos ver a diferenca ao analisarmos os metadados de cada layer. Aqui mostramos

apenas a parte relevante da sada do comando v.info, que e utilizado quando selecionamos
a opcao Metadata. O layer contorno_hd_permanente_LL apresenta, alem de boundaries,
areas e ilhas (
areas dentro de
areas):
|
|
|
|
|
|
Type of Map:
Number
Number
Number
Number
of
of
of
of
vector (level: 2)
points:
lines:
boundaries:
centroids:
0
0
9130
6970
Number
Number
Number
Number
of
of
of
of
areas:
islands:
faces:
kernels:
8434
8045
0
0
|
|
|
|
|
|
Ja o layer contorno_hd_permanente_LL_lines possui apenas linhas:

|
|
|
|
|
|
Type of Map:
Number
Number
Number
Number
of
of
of
of
vector (level: 2)
points:
lines:
boundaries:
centroids:
0
16790
0
0
Number
Number
Number
Number
of
of
of
of
areas:
islands:
faces:
kernels:
0
0
0
0
|
|
|
|
|
|
8.1
Arquivos ASCII
8.1.
78
Arquivos ASCII
A importac
ao de pontos vetoriais a partir de arquivos texto e uma tarefa bastante comum
em SIG. S
ao arquivos de texto puro (tambem chamado de ASCII32 ), sem formatacao (negrito,
italico, etc). Podem ter extens
ao .txt (text), .dat (data), .csv (comma separated values),
.asc, entre outras. Os dados s
ao dispostos um por linha, com os atributos (latitude, longitude,
elevacao, etc) separados por espaco, vrgula, tabulacao, ou outros caracteres especiais, como
&e |.
Ao inves de simplesmente importar um arquivo XYZ qualquer, vamos aproveitar esta
oportunidade para ver mais alguns comandos do GRASS. Vamos criar um arquivo de pontos
vetoriais aleatoriamente distribudos em uma area, atribuir um valor de elevacao para cada
ponto (a partir do modelo SRTM), exportar esses pontos em formato ASCII e finalmente
importar o arquivo XYZ.
8.1.1
Gerac
ao de pontos aleat
orios
Primeiramente, vamos ajustar os limites da region ativa para concidir com os limites
do mapa SRTM. Com o mapas selecionados no gism, va no Map Display em Zoom to
selected map layer(s) e depois em Zoom to... Set computational region from
display. Se quiser conferir se a operacao funcionou, use a opcao Zoom to... Zoom to
computational region e confira se a area exibida n
ao se altera. Se preferir, voce pode fazer
o mesmo digitando, na janela de terminal g.region -p rast=srtm_patch:
GRASS 6.4.1
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
(cursoUFAM):~ > g.region -p rast=srtm_patch

0
wgs84
wgs84
0
5:00:02.27988S
65:00:02.63952W
54:59:58.07976W
0:00:02.99988
0:00:02.99988
6001
12002
72024002
Agora vamos criar os pontos aleatorios, com o comando v.random. Va em Vector

Generate points Generate random points. Crie um novo layer chamado pontos_rand, com 2500 pontos (Fig. 8.8).
32
http://pt.wikipedia.org/wiki/ASCII
8.1
79
Arquivos ASCII
Figura 8.8. Opco

es do comando v.random.
Certamente seu resultado ser

a diferente da figura 8.9, uma vez que a posicao dos pontos
e definida aleatoriamente. Se voce n
ao gostar da distribuicao dos pontos por algum motivo,
como areas de grande concentrac
ao de pontos e outras meio vazias, selecione a opcao Allow
output files to overwrite existing files, na aba Optional, e rode o comando novamente,
ate atingir um resultado satisfat
orio.
Figura 8.9. Pontos aleat

orios gerados com v.random.
8.1
Arquivos ASCII
8.1.2
80
Amostragem da altitude
Agora vamos amostrar a altitude de cada ponto, a partir do modelo SRTM. Com o mapa
pontos_rand selecionado no gism, procure a opcao Show attribute columns , e clique no
cone correspondente. Voce dever
a ver a mensagem de erro (Fig. 8.10):
Figura 8.10. Mensagem de erro do gerenciador de tabelas de atributos do GRASS.
O que significa que o mapa foi criado, mas a tabela de atributos nao. Quando mapas
vetoriais s
ao criados dessa maneira, e preciso criar a tabela de atributos e povoa-la com uma
linha por categoria. Isto pode ser feito de duas maneiras, usando o gerenciador de tabelas do
GRASS ou o comando v.db.addtable (na verdade o gerenciador de tabelas apenas faz uma
interface com o v.db.addtable). Vamos usar o gerenciador de tabelas, que ja deve ter sido
aberto automaticamente ap
os clicarmos um OK na mensagem de erro da figura 8.10.
Va na aba Manage layers. Aqui voce pode gerenciar as conexoes de bancos de dados dos
layers vetoriais. Primeiro precisamos criar uma tabela e depois conecta-la ao mapa pontos_rand. No lado direito da janela, Table description, preencha o campo Table name com
pontos_rand (para seguir a convencao da tabela ter o mesmo nome do mapa) e em Key
column deixe cat (categoria). Clique em Create table. Nao feche a janela do gerenciador
de tabelas.
Com isso criamos a tabela, agora vamos conecta-la ao mapa. No lado esquerdo da janela,
Layer description, deixe Layer em 1, selecione a tabela que criamos em Table, mantenha
Key column como cat e marque a caixa da opcao Insert record for each category
into table, para popular a tabela. Clique em Add layer. Na barra inferior da janela uma
mensagem mostrar
a Number of loaded records:2500, e na parte superior em List of layers
veremos qua a tabela foi conectada usando o driver dbf (fig. 8.11).
Se voce for agora na aba Browse data do gerenciador de tabelas vera uma interface tipo
planilha, com uma linha para cada objeto do layer e suas colunas de atributos. Por enquanto
nossa tabela s
o tem a coluna cat, com valores sequenciais para os pontos.
8.1
Arquivos ASCII
81
Figura 8.11. Conex

ao de tabela com layer vetorial pelo gerenciador de tabelas do GRASS.
Agora sim podemos extrair a altitude de cada ponto e vamos fazer isso de duas maneiras,
incluindo valores na tabela e depois criando atributos de elevacao 3D que nao sao armazenados
na tabela.
Primeiro abra o gerenciador de tabelas e na aba Manage tables, na parte inferior da
janela em Manage columns, preencha Column name com elev, deixe o campo Data
type como integer (ser
ao valores inteiros, se fossem valores decimais usaramos float)
e clique em Add. Agora abra o comando v.what.rast em Vector Upadte point
attibutes from raster Sample raster maps at point locations. No dialogo do
comando v.what.rast, escolha pontos_rand como mapa vetorial de entrada, srtm_patch
como raster a ser questionado e elev como a coluna do banco de dados que sera atualizada
com os valores do raster em cada ponto.
Abra novamente o gerenciador de tabelas e veja que agora temos duas colunas de dados,
uma com a categoria de cada ponto e outra com o valor de elevacao extrado do SRTM
(fig. 8.12).
8.1
82
Arquivos ASCII
Figura 8.12. Coluna elev atualizada com os valores do SRTM.
Agora vamos usar o comando v.drape. Va em Vector Develop Map Create 3D

vector over raster. Na aba Required, especifique como mapa de entrada pontos_rand e
pontos_rand_Z como mapa resultante. Na aba Optional, deixa apenas point selecionado
em Feature type e escolha srtm_patch como raster a ser amostrado. O mapa gerado por
este comando e um vetor 3D, ou seja, os valores amostrados do raster sao armazenados como
eleva
c
ao, e n
ao como atributos na tabela.
Apos rodar este comando, se voce usar a ferramenta Query do Map Display (Query
raster/vector map(s) (display mode)), ao clicar sobre os pontos do mapa, voce vera na
aba Command console do gism, as propriedades de cada um, como coordenadas, categoria
e elevacao:
v.what --v -a map=pontos_rand_Z@manaus east_north=-60.727280,-2.225000 distance=0.055263
East: 60:43:38.208W
North: 2:13:30S
Map: pontos_rand_Z
Mapset: manaus
Type: Point
Id: 2482
Point height: 58.000000
Layer: 1
Category: 2482
Driver: dbf
Database: /home/user/grassdata/cursoUFAM/manaus/dbf/
Table: pontos_rand_Z_1
Key column: cat
cat : 2482
elev : 58
8.1
Arquivos ASCII
8.1.3
83
Exportac
ao como ASCII
Com o mapa de pontos criado e convertido para 3D, podemos exporta-lo como ASCII, com
o comando v.out.ascii. V
a em File Export vector Map ASCII points/GRASS
ASCII vector export. Na aba Required defina pontos_rand_Z como input e na aba
Optional, pontos_rand_Z_ascii.txt como resultado (output). Deixe o tipo de formato de
sada como point.
Se voce n
ao definir um diret
orio especfico para o arquivo de sada, ele sera gravado no
seu diretorio atual (descubra qual e usando pwd no terminal).
Apesar de em sistemas Linux/Mac OSX nao ser exatamente obrigatorio definir uma extensao para o arquivo (pois o sistema sabe que se trata de um arquivo texto ASCII e vai
usar um editor de textos simples para abr-lo), no Windows e necessario indicar uma extensao,
como .txt. Como boa pr
atica, use uma extensao que descreva seus dados, como .txt ou
.asc.
Abra o arquivo pontos_rand_Z com um editor de textos (Mousepad, por exemplo). Voce
vera as colunas separadas pelo sinal |, como X | Y | Z | cat:
-60.858969|-1.55195453|72|1
-61.27686897|-2.13385107|48|2
-62.32804517|-1.09972788|16|3
-57.53513346|-3.79940948|25|4
-64.65947378|-4.04771712|81|5
...
8.1.4
Importac
ao do arquivo ASCII
Vamos agora importar o arquivo ASCII (nosso objetivo inicial), com o comando
v.in.ascii. V
a em File Import vector data ASCII points/GRASS ASCII
vector import.
Na aba Required defina pontos_rand_Z_ascii como entrada e pontos_ascii_import
como resultado. Na aba Points, use a a coluna 3 como valor Z, para que o mapa seja criado
como 3D. Note as opc
oes Do not create table in points mode e Do not build topology in points
mode, muito usadas quando o arquivo tem muitos pontos (da ordem de dezenas de milhares
a milhoes), o que torna invi
avel trabalhar com tabelas DBF.
Por fim, exiba os dois mapas de pontos sobrepostos, no gism. Faca com que pontos_rand_Z fique por baixo. Altere as cores de linha e do preenchimento do smbolo para vermelho e mude tamanho do smbolo para 10. Altere o smbolo de pontos_ascii_import para
basic/diamond, com tamanho 5 (gig. 8.13).
8.1
Arquivos ASCII
Figura 8.13. Mapa de pontos importado com v.in.asciisobreposto ao mapa gerado com v.random.
84
85
9 - An
9.
An
9.1.
Hipsometria
Vamos iniciar a an
alise de terreno com o mapa hipsometrico, onde temos a representacao
do relevo em classes de altitude (no nosso caso, classes de 50 metros). Vamos utilizar os
comandos r.reclass para criar o mapa com as classes de altitude e o comando r.colors
para definir a cor de cada classe.
Va em Raster Change category values and labels Reclassify. Na aba Required, especifique o nome do mapa a ser reclassificado (no nosso caso, srtm_pacth) e o
nome do novo mapa a ser criado, como srtm_patch_hipso50. Na aba Optional, podemos
indicar um arquivo de texto que contenha as regras de classificacao ou entrar com os valores
manualmente. A figura 9.1 mostra a regra de classificacao utilizada.
As regras de classificac
ao devem ser digitadas de acordo com o seguinte esquema:
intervalo de dados
n
umero da categoria
r
otulo da categoria
No nosso exemplo, usamos as seguintes regras:

-50 thru 50 = 1 0 - 50m
50 thru 100 = 2 50 - 100m
100 thru 150 = 3 100 - 150m
150 thru 200 = 4 150 - 200m
200 thru 250 = 5 200 - 250m
250 thru 300 = 6 250 - 300m
300 thru 350 = 7 300 - 350m
350 thru 400 = 8 350 - 400m
400 thru 450 = 9 400 - 450m
450 thru 500 = 10 450 - 500m
500 thru 550 = 11 500 - 550m
Note que um intervalo dos dados e definido por menor valor thru maior valor. Outro
ponto importante e que o comando r.reclass na verdade n
ao gera um layer novo, mas
apenas uma tabela de reclassificac
ao que e interpretada pelo GRASS cada vez que e necessario
acessar o layer novo (reclassificado). Isto e feito para economizar espaco em disco e faz com
que, enquanto o layer reclassificado existir, nao seja possvel apagar ou renomear o layer
original, pois a reclassificac
ao depende dele para ser acessada.
9.1
86
Hipsometria
Figura 9.1. Reclassificaca

o dos dados SRTM em classes de 50 metros.
Para alterar as cores do mapa reclassificado, vamos usar o comando r.colors. Podemos
escolher esquemas pre-definidos (Raster Manage colors Color tables) ou utilizar
regras proprias (r.colors.rules, Raster Manage colors Color rules).
O comando r.colors oferece v
arias tabelas de cores pre-definidas:
aspect:
aspectcolr:
bcyr:
bgyr:
byg:
byr:
celsius:
corine:
curvature:
differences:
elevation:
etopo2:
evi:
gdd:
grey:
grey.eq:
grey.log:
grey1.0:
aspect oriented grey colors

aspect oriented rainbow colors
blue through cyan through yellow to red
blue through green through yellow to red
blue through yellow to green
blue through yellow to red
blue to red for degree Celsius temperature
EU Corine land cover colors
for terrain curvatures (from v.surf.rst and r.slope.aspect)
differences oriented colors
maps relative ranges of raster values to elevation color ramp
colors for ETOPO2 worldwide bathymetry/topography
enhanced vegetative index colors
accumulated growing degree days
grey scale
histogram-equalized grey scale
histogram logarithmic transformed grey scale
grey scale for raster values between 0.0-1.0
9.1
87
Hipsometria
grey255:
gyr:
haxby:
ndvi:
population:
precipitation:
precipitation_monthly:
rainbow:
ramp:
random:
rstcurv:
rules:
ryb:
ryg:
sepia:
slope:
srtm:
terrain:
wave:
grey scale for raster values between 0-255

green through yellow to red
relative colors for bathymetry or topography
Normalized Difference Vegetation Index colors
color table covering human population classification breaks
precipitation color table (0..2000mm)
precipitation color table (0..1000mm)
rainbow color table
color ramp
random color table
terrain curvature (from r.resamp.rst)
create new color table based on user-specified rules read from stdin
red through yellow to blue
red through yellow to green
yellowish-brown through to white
r.slope.aspect-type slope colors for raster values 0-90
color palette for Shuttle Radar Topography Mission elevation
global elevation color table covering -11000 to +8850m
color wave
As regras de classificac
ao podem ser digitadas de acordo da seguinte maneira:
n
umero da categoria
valor_RED
valor_GREEN
valor_BLUE
Onde os valores para as cores estao entre 0-255. Existem outras opcoes para definir
cores com r.colors, como utilizar porcentagens de Vermelho/Verde/Azul (RGB) ao inves
de valores entre 0-255, utilizar nomes de cores ou esquemas pre-definidos. Para ver todas as
possibilidades, consulte a p
agina de manual do comando.
Na figura 9.2 temos o mapa hipsometrico com a tabela de cores elevation.
Figura 9.2. Mapa hipsometrico com classes de 50 metros e tabela de cores elevation.
9.2
Relevo sombreado
88
O comando r.colors.rules fornece uma interface interativa, onde e possvel adicionar

as classes (ou intervalos de valores) e escolher as cores, bem como obter uma pre-visualizacao
do resultado final (fig. 9.3).
Figura 9.3. Escolha interativa das cores com r.colors.rules.
9.2.
Relevo sombreado
Vimos brevemente os mapas de relevo sombreado na secao 6.4.1.

Uma opc
ao interessante e visualizar o relevo sombreado com a escala de cores sobreposta.
Para isso, no gism, clique no cone Add various raster map layers e escolha Add shaded
relief map layer (fig. 9.4). Na janela do comando d.shademap, selecione srtm_patch_shade_045_25_z5_1km como imagem de relevo e srtm_patch_reclass como imagem a ser
sobreposta ao relevo drape map. Visualize o resultado no Map Display (fig. 9.5).
Figura 9.4. Adicione um layer de mapa sobreposto a relevo sombreado.
A aparencia da imagem produzida varia bastante de acordo com o azimute do iluminante.

Feicoes lineares do relevo (que podem estar ligadas à estruturas geologicas, por exemplo) sao
destacadas quando o iluminante est
a posicionado na direcao perpendicular à feicao. Caso o
9.3 Indices de forma
89
Figura 9.5. Mapa hipsometrico sobreposto a relevo sombreado.
interesse seja uma an

alise de elementos lineares interpretados em imagens de relevo sombreado,
deve-se produzir v
arios mapas, com iluminante posicionado em diferentes orientacoes.
A posic
ao do iluminante pode ate mesmo fazer com que tenhamos a impressao derelevo
invertido. Pessoas que moram no hemisferio sul, por exemplo, estao acostumadas com iluminacao vinda de norte em aerofotografias e em imagens orbitais. Um mapa de relevo sombreado
com iluminac
ao vinda de sul pode parecer invertido para alguns.
9.3. Indices de forma

Diversos par
ametros morfometricos podem ser usados para descrever a superfcie topogr
afica, tais como declividade da vertente (slope), orientacao da vertente (azimute ou aspecto,
aspect) e curvaturas associadas `
a forma concava/convexa da forma de relevo.
Para calcular esses ndices, vamos usar o comando r.slope.aspect (Raster Terrain
analysis Slope and aspect). Nosso mapa de entrada e srtm_patch. Na aba Output, defina os nomes dos mapas de sada, de acordo com o parametro calculado: srtm_patch_slope,
srtm_patch_aspect, srtm_patch_pcurv, srtm_patch_tcurv. Nao se esqueca de conferir as
configurac
oes da region ativa antes de executar o comando.
9.3.1
Declividade e Aspecto
Por padr
ao o mapa de declividade (fig. ??) tem valores em graus, ou opcionalmente em
porcentagem. O mapa de orientac
ao de vertentes (fig. ??) tem valores iniciando em Leste,
crescendo em sentido anti-hor
ario: 90 corresponde a Norte, 180 a Oeste, 270 a Sul e 360
a Leste. O valor 0 (zero) e usado para identificar areas de orientacao indefinida quando a a
90
superfcie e plana (declividade zero). Caso a unidade dos valores de elevacao do layer original
nao esteja em metros, e preciso usar a opcao zfactor para converter os valores.
Figura 9.6. Mapa de declividade.
Figura 9.7. Mapa de aspecto (orientaca

o de vertentes).
Dependendo da aplicac
ao, o usu
ario pode necessitar de uma mapa com declividade agrupadas em classes. Isto pode ser feito com o comando r.reclass. Um exemplo de classificacao
e dado na p
agina de manual do comando r.slope.aspect:
categoria
1
2
3
4
5
6
intervalo
(em graus)
0-1
2-3
4-5
6-8
9-11
12-14
r
otulos (labels)
(em porcentagem)
0-2%
3-5%
6-10%
11-15%
16-20%
21-25%
15-90
91
26% e acima
Os calculos deste comando s

ao realizados de acordo com as formulas de Horn (1981). Para
a determinac
ao da declividade e do aspecto, o algoritmo utiliza uma matriz 3x3 ao redor de
cada pixel, portanto n
ao e possvel calcular os valores para as celulas adjacentes às bordas do
layer. Essas celulas recebem valores de declividade e aspecto zero.
De acordo com o manual do comando, caso o layer utilizado para calcular o aspecto seja do
tipo integer (apenas valores inteiros), deve ocorrer maior freq
uencia de valores nas direcoes
0 , 45 , 90 , 180 , 225 , 270 , 315 e 360 . Caso o layer seja do tipo float (valores fracionais),
esse tipo de problema n
ao deve acontecer. Veremos como resolver esse problema sa secao
10.6.
Celulas com declividade muito baixa tambem terao maior freq
uencia de valores nas
direcoes 0 , 45 , 90 , 180 , 225 , 270 , 315 e 360 . Para evitar que isso ocorra, pode-se
utilizar a opc
ao min_slp_allowed, que fara com que todas as celulas com declividade menor
que o valor estabelecido recebam valores nulos.
Como vimos, o mapa de aspecto produzido pelo GRASS tem valores entre 0-360 iniciando
no Leste e crescendo em sentido anti-horario. Essa organizacao dos valores e esperada por
alguns outros comandos (como os de modelagem hidrologica), e enquanto que a escala de cores
utilizada e muito boa para interpretacoes visuais, se quisermos fazer uma analise estatstica
ou correlac
ao com o mundo real, devemos transformar os valores para a orientacao geografica
(compass-oriented aspect values). Podemos fazer isso com a ferramenta de algebra de mapas
raster, r.mapcalc. V
a em Raster Raster map calculator. Aqui entramos com uma
expressao que ser
a avaliada pelo comando e executada, para um ou mais rasters. Pode-se
utilizar desde express
oes simples ate esquemas complexos, envolvendo testes logicos. Uma
boa leitura sobre essa poderosa ferramenta e (Shapiro & Westervelt, 1991).
No campo maior da parte inferior da janela, Expression, podemos digitar as expressoes
ou podemos utilizar a terceira opc
ao da parte direita da janela, Insert mapcalc function
para nos auxiliar. Procure por if(x,a,b) e a insira. Essa expressao ira primeiro verificar a
condicao x. Se x for verdadeiro, ent
ao a instrucao a sera executada. Se x for falso, a instrucao
b sera executada.
92
Para converter os valores de aspecto, usaremos a seguinte formula:

if(A <= 90, 90-A, 360+90-A)
Portanto se os valores do mapa A forem menores ou iguais a 90, iremos subtrai-los de 90,
caso os valores sejam maiores do que 90, eles serao subtrados de 450.
Preencha os campos da janela (fig: 9.8) de modo que o mapa de sada seja srtm_patch_aspect_compass e a express
ao seja:
if(srtm_patch_aspect@manaus <= 90, 90 - srtm_patch_aspect@manaus ,
450 - srtm_patch_aspect@manaus)
Figura 9.8. Janela do r.mapcalculator.
Apos a execuc
ao do r.mapcalc, o GRASS atribui ao novo mapa uma escala de cores
padrao, mas voce pode mudar isso usando r.colors e atribuir, por exemplo, a escala aspect
ou aspectcolr. Se fizer isso, ir
a notar que a escala de cores aspect deixa o mapa um pouco
estranho, ate com aparencia de relevo invertido (fig: 9.9). Como falamos na secao 9.2, estamos
acostumados a imagens com iluminante a partir do norte, entao podemos acertar essa escala
de cores usando o seguinte esquema em r.colors:
0 white
180 black
360 white
93
Figura 9.9. Mapas aspecto com escala de cores aspect (esquerda) e definida manualmente (direita).
9.3.2
Curvaturas
Alem da declividade e da orientacao de vertentes, podemos usar o comando

r.slope.aspect para calcular a curvatura vertical (ou curvatura de perfil, profile curvature),
a curvatura horizontal (ou curvatura tangencial, tangential curvature), e as derivadas parcias
de primeira e segunda ordem nas direcoes X e Y.
A curvatura vertical e a taxa de variacao da declividade medida na direcao perpendicular `
a
vertente (na direc
ao do aspecto), e exprime o comportamento de aceleracao/desaceleracao de
fluxos sobre a superfcie topogr
afica (fig. 9.10). Valores positivos indicam curvaturas concavas
(desacelerac
ao de fluxo) e valores negativos indicam curvaturas convexas (aceleracao de fluxo).
As curvaturas tem valores expressos na forma 1/metros, ou seja, um valor de 0.05 significa
um raio de curvatura de 20m.
A curvatura horizontal e medida na direcao perpendicular à declividade da vertente, e
exprime a relac
ao de convergencia/divergencia de fluxos sobre a superfcie (fig. 9.10). Valores
positivos indicam vertentes c
oncavas (fluxo convergente) e valores negativos indicam vertentes
convexas (fluxo divergente).
Para mais informac
oes sobre o calculo dos ndices de forma discutidos, recomenda-se a
consulta de Ruhe (1975); Evans (1980); Zevenbergen & Thorne (1987); Wood (1996); Florinsky
(1998); Florinsky et al. (2002); Shary et al. (2002).
Figura 9.10. Relaco

es entre as curvaturas Vertical e Horizontal
Figura 9.11. Mapas de curvatura vertical (profile curvature, esquerda) e horizontal (tangential curvature,
direita).
94
95
10 - Interpolac
ao de superfcies
10.
Interpolac
ao de superfcies
Uma das tarefas mais comuns em SIGs e a interpolacao de valores pontuais (distribudos
regular ou irregularmente no espaco) em uma superfcie contnua, que pode ser usada em
diversas modelagens e simulac
oes. Provavelmente o exemplo mais comum seja a geracao de
Modelos Digitais de Elevac
ao (MDEs), a partir de valores pontuais de elevacao ou curvas de
nvel digitalizadas de mapas topogr
aficos.
Entre as aplicac
oes principais para interpolacao estao:
Alterar a resoluc
ao espacial de dados raster (reamostragem)
Metodos dos vizinhos mais proximos (nearest neighbor - NN)
Interpolac
ao por Splines (Regularized Splines with Tension - RST)
Preenchimento de
areas sem dados (interpola
c
ao)
Inverso da potencia da distancia (Inverse Distance Weighted - IDW)
Interpolac
ao por Splines (Regularized Splines with Tension - RST)
10.1.
Reamostragem
A reamostragem de mapas raster pode ser feita de varias maneiras, como por vizinhos
mais proximos (m
odulo r.resample), por Splines (r.resamp.rst) ou outros metodos de
interpolac
ao (r.resamp.interp). No primeiro caso, o mapa resultante tera a resolucao da
region ativa, e o valor de cada pixel sera dado pelo valor do pixel original que estiver mais
proximo do centro do novo pixel (figs. 10.1 e 10.1). Nos demais casos, pode-se ajustar diversos
parametros da interpolac
ao para um melhor resultado.
10
10
Reamostragem
Figura 10.1. Reamostragem (mudanca na resoluca

o espacial) por vizinhos mais pr
oximos.
10.2
96
Inverso da potencia da dist

ancia - IDW
10
10 11 12 13
14 15 16 17
14
Reamostragem
14
18 19 20 21
18
22 23 24 25
11
15
19
23
12
16
20
24
17
21
25
24
Figura 10.2. Reamostragem (transformaca

o geometrica) por vizinhos mais pr
oximos.
10.2.
Inverso da pot
encia da dist
ancia - IDW
Trata-se de um metodo de interpolacao local, que se baseia na premissa de que pontos

proximos possuem valores mais similares que pontos distantes. A variacao entre os valores
e modelada segundo potencias (quadrado ou cubo) do inverso da distancia entre os pontos
(Fig. 10.3).
Modulos do GRASS para interpolacao por IDW (o resultado sera sempre um mapa raster):
v.surf.idw dados vetoriais (pontos, curvas de nvel vetoriais)
r.surf.idw dados raster
ponto a ser
estimado (Z')
raio de busca
Z1
d1
Z2
d2
d6
Z6
Z5
d3
d5
d4
Z3
Z4
2
Z'=(1/d )Z / 1/d
Figura 10.3. Esquema da interpolaca

o por IDW.
10.3
Splines regularizadas com tens

ao - RST
10.3.
97
Splines regularizadas com tens

ao - RST
O metodo de interpolac
ao/aproximacao por splines tambem e apropriado para o preenchimento de vazios nos dados (como a area de um lago, por exemplo) A ideia e gerar uma
superfcie que passe exatamente pelos pontos de dados ou proximo deles o suficiente para
gerar uma superfcie contnua e que preencha os espacos vazios (data blanks, voids). O ajuste
dos parametros de interpolac
ao muda o comportamento da superfcie calculada, desde uma
membrana el
astica passando exatamente por todos os pontos ate uma superfcie mais rgida,
que nao passa pelos pontos e define apenas uma tendencia geral dos dados. Para um melhor
resultado, recomenda-se a leitura das paginas de manual dos modulos.
Informac
oes adicionais e mais detalhadas podem ser encontradas nos seguintes trabalhos:
Smith & Wessel (1990); Mitasova & Mitas (1993); Mitasova & Hofierka (1993); Hofierka et al.
(2002).
Modulos do GRASS para interpolacao por RST:
v.surf.rst dados vetoriais (pontos, curvas de nvel vetoriais)
v.vol.rst para gerac
ao de volumes tridimensionais (grid 3D)
10.4.
Interpolac
ao bilinear e bic
ubica
Este metodo e bastante utilizado para a reamostragem de dados raster e para a geracao
de superfcies a partir de bandos de dados muito volumosos como os obtidos por altimetria
laser:
r.resamp.interp reamostragem de dados raster
v.surf.bspline pontos vetoriais
10.5.
Outros m
odulos de interpolac
ao
Outros m
odulos do GRASS para interpolacao de dados raster:
r.surf.area Estimac
ao de
areas para mapas raster
r.surf.contour C
alculo de superfcies (MDEs) a partir de curvas de nvel em formato
raster
r.surf.fractal Gera uma superfcie fractal de determinada dimensao fractal
r.surf.gauss Usa um gerador de n
umeros aleatorios gaussiano para gerar uma superfcie de valores de distribuicao gaussiana com media e desvio padrao determinados
pelo usu
ario
10.6
98
Reamostragem dos dados SRTM
r.surf.random Gera um mapa raster com valores aleatorios com intervalo (range)
definido pelo usu
ario.
10.6.
Reamostragem dos dados SRTM
Vamos aproveitar a oportunidade e corrigir alguns problemas do mapa de aspecto que

fizemos anteriormente. Como utilizamos o modelo SRTM para calcular o aspecto, e os dados
originais est
ao em formato de n
umeros inteiros, temos artefatos nas direcoes 0 , 45 , 90 ,
180 , 225 , 270 , 315 e 360 . Vamos fazer uma reamostragem dos dados originais mas sem
alterar muito a resoluc
ao espacial. Isso ja sera suficiente para obter uma superfcie um pouco
mais contnua e com dados tipo float.
No terminal, digite g.region -p rast=srtm_patch e veja que o tamanho dos pixels nao
e exatamente 0:0:03:
GRASS 6.4.1
projection:
zone:
datum:
ellipsoid:
north:
south:
west:
east:
nsres:
ewres:
rows:
cols:
cells:
(cursoUFAM):~ > g.region -p rast=srtm_patch

0
wgs84
wgs84
0
5:00:02.27988S
65:00:02.63952W
54:59:58.07976W
0:00:02.99988
0:00:02.99988
6001
12002
72024002
A mudanca da resoluc
ao espacial para 0:0:03 ja dever ser suficiente.
terminal, rode g.region -pa rast=srtm_patch res=0:0:03.
Ainda no
Agora abra o comando
r.resamp.interp em Raster Develop raster map Resample using multiple

methods. Nosso mapa de entrada sera srtm_patch e o de sada srtm_patch_bilinear.
Aqui bilinear corresponde ao metodo de interpolacao que pode ser escolhido na aba Optional.
Com a reamostragem completa, refaca o mapa de aspecto, usando o mapa reamostrado
como entrada e srtm_patch_interp_aspect como sada. Examine o histograma de cada um
desses mapas com a ferramenta Analize Create histogram of image or raster file do
Map Display. Na figura 10.4 voce ve os histogramas dos mapas de aspecto a partir dos dados
originais e reamostrados.
10.7
Extrac
ao de curvas de nvel
99
Figura 10.4. Histogramas dos mapas de aspecto, a partir dos dados SRTM originais (esquerda) e reamostrados com r.resamp.interp (direita).
10.7.
Extrac
ao de curvas de nvel
Agora que j
a temos o modelo SRTM corrigido, podemos extrair as curvas de nvel do
MDE. Esta e uma operac
ao muito u
til quando se trabalha em areas carentes de levantamentos
topograficos de maior detalhe. Usaremos o comando r.contour, disponvel no menu Raster
Generate contour lines. Como opcoes deste comando, podemos selecionar exatamente
quais curvas de nvel queremos ou indicar os nveis mnimo e maximo, e o intervalo altimetrico
entre as curvas.
No exemplo usamos o raster srtm_patch_bilinear e criamos o mapa vetorial srtm_bilinear_contours.
Com base nas informacoes obtidas anteriormente com o comando
r.info, definimos como valor mnimo a cota de 0 m, e como maximo 540 m, com intervalo entre as curvas de 20 m, o mesmo utilizado nos mapas topograficos em escala 1:50.000.
Note que n
ao e necess
ario usar valores que estejam dentro do intervalo de valores do mapa.
No nosso caso, o valor mnimo do raster e de -26 m, mas utilizamos como cota mnima 0 m,
para obter um espacamento mais usual das curvas de nvel.
Dentre os shapefiles que baixamos do IBGE, importe o arquivo rel_curva_nivel_l, da
area 4-2009, e o compare as curvas que foram geradas por r.contour. A figura 10.5 mostra
os resultados obtido.
10.8.
Interpolac
ao de MDE
10.8.1
Recorte de dados vetoriais
Antes de partir para a interpolacao do DEM, vamos aproveitar e ver outra ferramenta
importante, que usa l
ogica booleana para sobreposicao de vetores. O comando v.overlay
(Vector Overlay vector maps Overlay vector maps) usa dois mapas vetoriais
10.8
Interpolac
ao de MDE
100
Figura 10.5. Curvas de nvel para cotas definidas entre 0 m e 540 m (esquerda) e comparaca
o com curvas
de escala 1:250.000 do IBGE (direita, curvas do IBGE em vermelho).
de entrada para produzir um mapa de sada a partir das operacoes and, or, not, xor
(fig. 10.6):
and operac
ao de intersec
c
ao. Os elementos do mapa A e de mapa B.
or operac
ao de uni
ao. Os elementos de A ou de B.
not Os elementos de A mas n
ao de B.
xor Os elementos de A ou de B, mas n
ao os de A e B
10.8
101
Interpolac
ao de MDE
Figura 10.6. Operaco

es de overlay.
Vamos usar a operac

ao de intersec
c
ao para recortar as curvas de nvel apenas em uma
region menor do que toda a extens
ao dos nossos dados:
g.region -p n=2:50S s=3:30S w=60:40W e=59:50W res=0:0:03 save=iranduba
Para o recorte precisamos de um polgono com os limites da region, e o comando
v.in.region faz justamente isso. V
a em Vector Generate area for current region e
crie o mapa region_manacapuru.
Agora use v.overlay. Na aba Required, entre com o mapa srtm_bilinear_contours
como entrada A (ainput), region_manacapuru como entrada B (binput) e curvas_srtm_manacapuru. Na aba Optional defina o tipo de ainput como line e a operacao de overlay
como and. Essa operac
ao deve demorar um pouco, ja que o arquivo de entrada e bem grande.
10.8.2
Interpolac
ao
Antes de iniciar a interpolac

ao e sempre bom checar os limites da region de interesse, j
a
que alguns par
ametros do comando tem seus valores padrao determinados a partir da resolucao
espacial da region ativa. Agora sim inicie o comando v.surf.rst, pelo menu Raster
Interpolate surfaces Regularized spline tension. Este comando tem muitas opcoes.
10.8
Interpolac
ao de MDE
102
Usaremos o mapa srtm_bilinear_contours como dado de entrada, mas antes precisamos

descobrir como os dados de altitude de cada curva de nvel estao armazenados, se estao na
tabela ou como atributo 3D. Use v.info no terminal ou a opcao de Metadada no gism para
ver as propriedades desse layer. Com v.info vemos que o mapa e 3D e tambem possui uma
tabela associada. Com v.describe podemos ver as caractersticas da tabela de atributos:
v.info srtm_bilinear_contours
Map is 3D:
Yes
Number of dblinks:
1
GRASS 6.4.1 (cursoUFAM):~ > db.describe srtm_bilinear_contours
table:srtm_bilinear_contours
description:
insert:yes
delete:yes
ncols:2
nrows:28
column:cat
description:
type:INTEGER
len:11
scale:0
precision:10
default:
nullok:yes
select:yes
update:yes
column:level
description:
type:DOUBLE PRECISION
len:20
scale:6
precision:18
default:
nullok:yes
select:yes
update:yes
Entao tambem temos uma coluna level na tabela de atributos com os valores de elevacao. Podemos escolher qualquer um dos dois para a interpolacao. Na janela do comando
v.surf.rst, entre com srtm_bilinear_contours na aba Required. Na aba Selection escolhemos de onde buscar os valores de elevacao: se a opcao Use z-coordinates (3D vector
only) for selecionada, mude o layer para 0 (zero). A opcao WHERE e utilizada para selecionar uma parte dos dados, como por exemplo: income < 1000 and inhab > 10000. Se for
10.8
Interpolac
ao de MDE
103
usar os dados da tabela, deixe layer em 1 e na aba Parameters entre com o nome da coluna
em zcolumn. Na aba Output, chame o mapa resultante de topografia_RST, e mantenha os
demais valores como os padr
oes sugeridos pelo programa (veja que na aba Parameters, as
opcoes dmax e dmin s
ao definidas automaticamente a partir do valor da resolucao da region
ativa). Ainda pode-se calcular ao mesmo tempo os mapas de declividade, orientacao de vertentes e curvaturas (sec
ao 9.3). Exiba o resultado no Map Display e compare com o modelo
SRTM (fig. 10.7).
Aparentemente os mapas s
ao muito similares. As cores sao diferentes porque no mapa
interpolado, o intervalo de valores para distribuir as cores e menor. Vamos ver isso com
mais detalhe usando a ferramenta Map Calculator para algebra de mapas raster (comando
r.mapcalc). No mapcalc, faca uma subtracao dos dois mapas (SRTM menos RST) e chame
o resultado de srtm_diff_rst. Use r.colors para aplicar a escala de cores difference e
exiba o resultado (fig. 10.8).
Figura 10.7. Resultado da interpolac

ao por RST (esquerda) e modelo SRTM (direita).
10.8
104
Interpolac
ao de MDE
Figura 10.8. Resultado da subtraca

o entre o modelo SRTM e a interpolaca
o por RST.
Veja as diferencas numericas com r.info para cada mapa:

r.info map=srtm_patch_bilinear
Range of data:
min = -26.5713794650012
r.info map=topografia_RST
Range of data:
min = -17.82484
max = 542.924676987079
max = 133.0811
r.info map=srtm_diff_rst
Range of data:
min = -22.2308455429801
max = 25.8248443603516
11 - Convers
ao entre projec
oes no GRASS
11.
105
Convers
ao entre projec
oes no GRASS
Ate agora trabalhamos em uma Location com coordenadas em Latitude-Longitude, porque nossos dados est
ao assim. E se quisermos trabalhar em outro sistema de coordenadas ou
outra projec
ao cartogr
afica?
No QGIS, temos a opc
ao de projecao on-the-fly, que converte as coordenadas e projeta os
dados de maneira transparente, sem necessidade de criar novos arquivos. No GRASS isso nao
e possvel, os dados devem estar todos em um mesmo sistema de coordenadas. Precisamos
entao converte-los antes de importar no GRASS ou projeta-los dentro do GRASS, entre duas
Locations.
Vamos usar o QGIS para ver como nossos dados ficariam no sistema UTM. No QGIS, abra
o raster do GRASS srtm_patch_bilinear. No canto direito inferior da janela do QGIS h
a
um cone de um globo e a indicac
ao do sistema de coordenadas no codigo EPSG (fig. 11.1).
Figura 11.1. Icone para ajuste do sistema de coordenadas (CRS - Coordinate Reference System) no QGIS.
Esse cone e um atalho para o di

alogo de Propriedades do Projeto, onde marcamos a opcao
Enable on the fly CRS transformation e escolhemos qual o sistema de coordenadas
vamos utilizar. O sistema UTM est
a em textbfProjected Coordinate Systems Universal
Transverse Mercator (UTM), mas qual das opcoes escolher? O datum pode ser tanto WGS84
como SIRGAS 2000, pois s
ao identicos para fins cartograficos, mas qual a zona UTM?
Para resolver essa quest
ao, feche (cancele) o dialogo e abra o shapefile do mapa-ndice
do mapeamento sistem
atico (que usamos na secao 7) VETOR_1000000, que mostra os mapas
em escala 1:1.000.000 do territ
orio brasileiro. Use a ferramenta de questao do QGIS para
identificar o c
odigo dos mapas que recobrem nossa area. Voce vera que a metade Oeste do
nosso SRTM est
a na folha SA-20 e a metade Leste na folha SA-21. Lembrando que a zona
UTM corresponde ao fuso da carta internacional ao milionesimo, nossos dados estao nas zonas
UTM 20 e 21, hemisferio Sul.
Mas isso n
ao responde completamente nossa questao, pois usaremos a zona 20 ou 21?
Faca um teste r
apido: abra o di
alogo das Propriedades do Projeto, marque a projecao on-
11 - Convers
ao entre projec
oes no GRASS
106
the-fly e selecione o sistema UTM, WGS84, zona UTM 21 Sul. Veja como as coordenadas
se comportam ao percorrer a
area com o ponteiro do mouse (se o mapa sumir da tela ap
os
aplicar a nova projec
ao, clique com o botao direito sobre o raster e selecione a opcao Zoom
to layer extent). Veja que na porc
ao oeste do mapa aparecem coordenadas negativas. Isso
acontece porque as coordenadas UTM aumentam para leste e nossos dados se estendem para
oeste alem do limite da zona UTM 21.
Refaca o teste com a zona UTM 20 Sul. Veja que as coordenadas nao apresentam mais
valores negativos. Aproveite e anote os valores aproximados das coordenadas nos limites do
raster.
Note que nossos dados ainda ultrapassam os limites da zona UTM 20, ou seja, nossa
area e grande demais para trabalharmos em UTM. O correto seria utilizar uma projecao que
atenda às necessidades de uma
area grande (como uma projecao conica). Para este exerccio,
usaremos mesmo a projec
ao UTM, mas lembre-se que isso nao e exatamente correto do ponto
de vista cartogr
afico.
Vamos ent
ao criar uma nova Location (por exemplo, cursoUFAMutm) e um novo mapset
em coordenadas UTM. Feche todas as janelas do GRASS e o inicie novamente. Na janela
inicial, refaca os passos da sec
ao 5.1. Utilize o datum WGS-84, resolucao de 90 m, e use os
limites do raster srtm_patch_bilinear que anotamos no QGIS (n=10.000.000, s=9.440.000,
w=270.000, e=1.400.000). Crie um mapset e inicie a sessao do GRASS.
Para reprojetar o SRTM, v
a em Raster Develop raster map Reproject raster
map. Na janela do comando r.proj, na aba Required, informe a Location que contem o
mapa a ser reprojetado. Na aba Source indique o mapset e qual mapa sera utilizado. Na aba
Target, defina o nome do mapa de sada (por padrao o mesmo nome do mapa de entrada),
o metodo de interpolac
ao a ser aplicado (use bilinear) e a resolucao do mapa resultante (90
m).
Reprojete tambem os mapas topografia_RST e srtm_diff_rst.
Se quiser, antes de rodar o comando, voce pode examinar os limites do mapa original no
novo sistema de coordenadas com a opcao p ou g do r.proj:
r.proj -g input=srtm_patch_bilinear location=cursoUFAM mapset=manaus method=bilinear resolution=90
n=10000000 s=9446905.82560594 w=278154.50026392 e=1393219.77169675 rows=6001 cols=12002
Input Projection Parameters: +proj=longlat +no_defs
+a=6378137 +rf=298.257223563 +towgs84=0,0,0,0,0,0,0
Input Unit Factor: 1
Output Projection Parameters: +proj=utm +south +no_defs
+zone=20 +a=6378137 +rf=298.257223563 +towgs84=0,0,0,0,0,0,0
Output Unit Factor: 1
Input map <srtm_patch_bilinear@manaus> in location <cursoUFAM>:
107
12 - Visualizac
ao n-dimensional
12.
Visualizac
ao n-dimensional
Nesta sec
ao veremos um pouco sobre o visualizador n-dimensional do GRASS, o NVIZ.
Nele e possvel visualizar superfcies (pixels - 2.5D) e volumes (voxels - 3D), bem como dados
vetoriais.
A experiencia mostra que e melhor trabalhar no NVIZ em Locations definidas com sistemas de coordenadas planas, pois n
ao ha necessidade de conversao dos valores de altitude. O
NVIZ usa as definic
oes da region ativa, e pode demorar para processar os dados dependendo
do tamanho da
area e resoluc
ao espacial. Escolha uma sub-area no mapa com a ferramenta
Zoom e defina as definic
oes da region a partir do Map Display.
Inicie o NVIZ em File NVIZ (requires Tcl/Tk). Selecione o raster srtm_patch_bilinear na opc
ao elevation e rode o comando. Voce devera ver uma janela como a da
figura 12.1.
Figura 12.1. Janela do NVIZ.
A posic
ao do observador pode ser controlada arrastando a seta para que ela aponte para
a direcao que se deseja observar. A perspectiva aproxima ou afasta a cena, a altura (height)
modifica sua posic
ao na vertical e a opcao z-exag altera a escala vertical da superfcie.
Experimente um pouco os comandos para se acostumar com eles. Lembre-se que quanto
maior a resoluc
ao espacial da region ativa, mais pesado fica para o programa, entao por
vezes voce poder
a ter que clicar em Draw para forcar a imagem a ser redesenhada.
Na janela do NVIZ, no menu Appearance, existem varias opcoes para inserir/modificar
mapas raster, vetoriais, iluminac
ao, etc. Selecione a opcao Lighting e altere a posicao da
fonte de luz do modelo (fig. 12.2). Tambem e possvel controlar a cor da luz ambiente a partir
dos controles de vermelho, verde e azul.
12 - Visualizac
ao n-dimensional
108
Figura 12.2. Controles de iluminaca

o do NVIZ.
Agora vamos alterar a cor da superfcie exibida. Abra o painel de superfcies raster
(textbfVisualize Raster surfaces), clique no sub-menu Surface attributes e depois em
color:srtm_patch_bilinear@manaus (fig. 12.3). Na janela que abrir clique em New Map,
para selecionar outro raster que ser
a exibido sobreposto à superfcie da topografia (drape).
Na proxima janela, clique sobre o nome do mapset onde esta o mapa (no nosso caso so temos
uma possibilidade, manaus) e selecione o mapa srtm_diff_rst .
Figura 12.3. Selecione para alterar a cor da superfcie.
Na figura 12.4 vemos o resultado da sobreposicao do mapa de diferencas entre o modelo

SRTM e o interpolado por RST sobre a topografia, onde pode-se ver bem em quais areas estao
as maiores diferencas entre as superfcies. Veja que qualquer mapa raster pode ser sobreposto
a outro no NVIZ, seja um modelo de terreno, imagem de satelite ou fotografia aerea. Veja
tambem que onde havia a seta para definir a direcao de observacao, ha uma cruz, que foi
12 - Visualizac
ao n-dimensional
109
selecionada pela opc

ao center (acima da cruz, marcada em vermelho), usada para alterar o
posicionamento do centro da imagem observada.
Figura 12.4. Sobreposic

ao do mapa de diferencas a
` superfcies topogr
afica.
13 - O QGIS como interface gr

afica alternativa
13.
110
O QGIS como interface gr

afica alternativa
O QGIS, alem de poder trabalhar nativamente com formatos ShapeFile, Geotiff e varios
outros, possui um plugin para interagir com dados e ferramentas do GRASS, o que o torna
interessante como interface gr
afica alternativa ao gism.
Ao iniciar o QGIS, v
a ao menu Plugins Manage Plugins e selecione o plugin GRASS.
Com o plugin selecionado, teremos acesso à ferramentas para incluir mapas raster ou vetoriais
na area de trabalho do QGIS. Algumas ferramentas so se tornam ativas depois que tivermos
includo algum mapa no QGIS.
Vamos comecar incluindo um mapa raster. Use o cone da barra de ferramentas ou v
a
em Plugins GRASS Add GRASS raster layer. Na janela que se abre, selecione a
Location cursoUFAM, o Mapset manaus, e o mapa srtm_patch_bilinear. O layer sera exibido
no painel central do QGIS e o nome sera includo no painel da esquerda (Layers). Adicione
tambem o layer vetorial HID_Limite_Massa_Dagua_L. Note que no painel Layers, a ordem
vertical dos layers funciona da mesma maneira que no GRASS. Aqui tambem podemos ligar
e desligar a visualizac
ao do layer ou alterar suas propriedades, com um duplo-clique sobre
seu nome.
De um duplo-clique sobre o nome do layer HID_Limite_Massa_Dagua_L, e altere a cor de
exibicao da linha para azul (Fig. 13.1). Voce pode fazer isso no botao Change que possui
um quadrado colorido indicando a cor da linha (na area central da janela) ou no botao
Change com uma ferramenta (na porcao esquerda da janela), onde voce pode alterar outras
propriedades, como o tipo de linha (tracejada etc).
Figura 13.1. Altere a cor da linha do mapa vetorial.

afica alternativa
111
Figura 13.2. Mapas raster e vetorial exibidos no QGIS.
Por enquanto, apenas podemos visualizar os mapas do GRASS, mas a ferramentaGRASS

Tools, que d
a acesso `
a diversos comandos do GRASS a partir do QGIS, nao esta acessvel.
Para que esta ferramenta seja ativada, temos duas opcoes: iniciar o QGIS a partir do terminal
do GRASS (simplesmente digite qgis no terminal) ou usar o comando Open Mapset do plugin
e selecionar qual Mapset queremos acessar. Como ja estamos trabalhando no QGIS, vamos
usar a segund opc
ao. V
a em Plugins GRASS Open Mapset e na janela que se
abrir, selecione a Location cursoUFAM e o Mapset manaus. Com isso o GRASS Tools deve
ficar disponvel. Abra o GRASS Tools (Fig. 13.3).
Figura 13.3. Janela do GRASS tools.

afica alternativa
112
A janela do GRASS Tools tem tres abas: Modules Tree, Modules List e Browser. As
abas modules, d
ao acesso a v
arios comandos do GRASS de maneiras diferentes, e na aba
browser podemos navegar pela
arvore de layers raster e vetoriais, visualizar as informacoes
sobre cada layer (figs. 13.4), inclu-los na area de trabalho com um clique no botao direito
(opcao Add selected map to canvas), apaga-los (permanentemente, cuidado!), renomea-los ou
definir os limites da region ativa de acordo com o layer desejado.
Figura 13.4. Mapas rasters disponveis e propriedades do layer srtm_diff_rst.
14 - Produzindo mapas para impress

ao
14.
113
Produzindo mapas para impress

ao
O comando ps.map permite que criemos mapas complexos, que serao salvos em arquivos
tipo postscript 33 . As informac
oes necessarias para construir o mapa (escala, grade de coordenadas, etc) s
ao salvas em um arquivo de texto simples, o que permite que criemos varios
mapas com exatamente o mesmo layout, apenas alterando o(s) layer(s) raster ou vetoriais
desejados.
Abaixo vemos um exemplo de um script simples para ps.map. O comando possui diversas
opcoes, que podem ser consultadas na sua pagina de manual. Salve o texto abaixo em um
arquivo texto e rode o comando ps.map, encontrado em File Postscript plot. O mapa
produzido pode ser visualizado, por exemplo, com o aplicativo Documento Viewer no Linux
(fig. 14.2).
rast srtm_patch_bilinear
vlines HID_Limite_Massa_Dagua_L
color blue
width 0.1
end
grid 20000
color black
numbers 2
end
colortable y
raster srtm_patch_bilinear
range 0 540
width 0.2
height 3.5
where 7 3.0
color black
nodata y
end
scale 1:500000
paper a4
end
end
33
http://en.wikipedia.org/wiki/PostScript
14 - Produzindo mapas para impress

ao
Figura 14.1. Janela do comando ps.map
Figura 14.2. Mapa visualizado no Documento Viewer
114
15 - Comandos do GRASS
15.
115
Comandos do GRASS
A listagem abaixo apresenta todos os comandos do GRASS 6.4. Listagens atualizadas podem ser
encontradas no site oficial do GRASS34 , na secao Docs Manuals.
comandos d.* (display)
d.ask
Prompts the user to select a GRASS data base file from among files displayed
in a menu on the graphics monitor.
d.barscale
Displays a barscale on the graphics monitor.
d.colorlist
Output a list of all available display colors with a configurable separator (default
is comma).
d.colors
Allows the user to interactively change the color table of a raster map layer
displayed on the graphics monitor.
d.colortable
Displays the color table associated with a raster map layer.
d.correlate
Prints a graph of the correlation between data layers (in pairs).
d.erase
Erase the contents of the active display frame with user defined color
d.extend
Sets window region so that all currently displayed raster and vector maps can
be shown in a monitor.
d.extract
Selects and extracts vectors with mouse into new vector map.
d.font
Selects the font in which text will be displayed on the users graphics monitor.
d.frame
Manages display frames on the users graphics monitor.
d.geodesic
Displays a geodesic line, tracing the shortest distance between two geographic
points along a great circle, in a longitude/latitude data set.
d.graph
Program for generating and displaying simple graphics on the display monitor.
d.grid
Overlays a user-specified grid in the active display frame on the graphics monitor.
d.his
Displays the result obtained by combining hue, intensity, and saturation (his)
values from user-specified input raster map layers.
d.histogram
Displays a histogram in the form of a pie or bar chart for a user-specified raster
map.
d.info
Display information about the active display monitor
d.labels
Displays text labels (created with v.label) to the active frame on the graphics
monitor.
d.legend
Displays a legend for a raster map in the active frame of the graphics monitor.
d.linegraph
Generates and displays simple line graphs in the active graphics monitor display
frame.
d.mapgraph
Generates and displays simple graphics on map layers drawn in the active
graphics monitor display frame.
d.measure
Measures the lengths and areas of features drawn by the user in the active
display frame on the graphics monitor.
d.menu
Creates and displays a menu within the active frame on the graphics monitor.
34
http://grass.osgeo.org
d.m
d.mon
d.monsize
d.mvmon
d.nviz
d.out.file
d.out.gpsdrive
d.out.png
d.path
d.polar
d.profile
d.rast.arrow
d.rast.edit
d.rast
d.rast.leg
d.rast.num
d.redraw
d.resize
d.rgb
d.rhumbline
d.save
d.shadedmap
d.slide.show
d.split.frame
d.split
d.text
d.thematic.area
d.title
d.vect.chart
d.vect
d.vect.thematic
d.what.rast
116
Display manager for GRASS

To establish and control use of a graphics display monitor.
Selects/starts specified monitor at specified size
Moves displayed maps to another monitor
Creates fly-through script to run in NVIZ.
Saves the contents of the active display monitor to a graphics file.
Export display monitor to a GpsDrive compatible backdrop image
Saves active display monitor to PNG file in home directory
Finds shortest path for selected starting and ending node.
Draws polar diagram of angle map such as aspect or flow directions
Interactive profile plotting utility with optional output.
Draws arrows representing cell aspect direction for a raster map containing
aspect data.
Interactively edit cell values in a raster map.
Displays user-specified raster map in the active graphics frame.
Displays a raster map and its legend on a graphics window
Overlays cell category values on a raster map layer displayed to the graphics
monitor.
Redraws the current display frame in the GRASS monitor
Resizes active display monitor
Displays three user-specified raster maps as red, green, and blue overlays in the
active graphics frame.
Displays the rhumbline joining two user-specified points, in the active frame on
the users graphics monitor.
Creates a list of commands for recreating screen graphics.
Drapes a color raster over a shaded relief map using d.his
Slide show of GRASS raster/vector maps.
Split the display into subframes.
Divides active display into two frames & displays maps/executes commands in
each frame.
Draws text in the active display frame on the graphics monitor using the current
font.
Displays a thematic vector area map in the active frame on the graphics monitor.
Create a TITLE for a raster map in a form suitable for display with d.text.
Displays charts of vector data in the active frame on the graphics monitor.
Displays user-specified vector map in the active graphics frame.
Displays thematic vector map
Allows the user to interactively query the category contents of multiple raster
map layers at user specified locations within the current geographic region.
d.what.vect
d.where
d.zoom
117
Allows the user to interactively query a vector map layer at user-selected locations within the current geographic region.
Identifies the geographic coordinates associated with point locations in the
active frame on the graphics monitor.
Allows the user to change the current geographic region settings interactively,
with a mouse.
comandos db.* (database)

db.columns
db.connect
db.copy
db.describe
db.drivers
db.dropcol
db.droptable
db.execute
db.in.ogr
db.login
db.out.ogr
db.select
db.tables
db.test
List all columns for a given table.

Prints/sets general DB connection for current mapset and exits.
Copy a table.
Describes a table in detail.
Lists all database drivers.
Drops a column from selected attribute table
Drops an attribute table.
Executes any SQL statement.
Imports attribute tables in various formats.
Sets user/password for driver/database.
Exports attribute tables into various formats.
Selects data from attribute table (performs SQL query statement(s)).
Lists all tables for a given database.
Test database driver, database must exist and set by db.connect.
comandos g.* (general )

g.access
g.ask
g.copy
g.dirseps
g.extension
g.filename
g.findetc
g.findfile
g.gisenv
g.gui
g.list
g.manual
g.mapset
g.mapsets
g.message
g.mkfontcap
g.mlist
Controls access to the current mapset for other users on the system.
Prompts the user for the names of GRASS data base files.
Copies available data files in the users current mapset search path and location
to the appropriate element directories under the users current mapset.
Internal GRASS utility for converting directory separator characters.
Tool to maintain GRASS extensions in local GRASS installation.
Prints GRASS data base file names.
Searches for GRASS support files.
Searches for GRASS data base files and sets variables for the shell.
Outputs and modifies the users current GRASS variable settings.
Launches a GRASS graphical user interface (GUI) session.
Lists available GRASS data base files of the user-specified data type.
Display the HTML man pages of GRASS
Changes current mapset.
Modifies the users current mapset search path.
Prints a message, warning, progress info, or fatal error in the GRASS way.
Generates the font configuration file by scanning various directories for fonts
Lists available GRASS data base files of the user-specified data type optionally
using the search pattern.
g.mremove
g.parser
g.pnmcomp
g.proj
g.region
g.remove
g.rename
g.setproj
g.tempfile
g.transform
g.version
118
Removes data base element files from the users current mapset.
Provides automated parser, GUI, and help support for GRASS scipts.
Overlays multiple PPM image files
Converts co-ordinate system descriptions (i.e. projection information) between
various formats (including GRASS format).
Manages the boundary definitions for the geographic region.
Removes data base element files from the users current mapset.
Renames data base element files in the users current mapset.
Interactively reset the locations projection settings.
Creates a temporary file and prints the file name.
Computes a coordinate transformation based on the control points.
Displays version and copyright information.
comandos i.* (imagery)

i.atcorr
i.cca
i.class
i.cluster
i.fft
i.fusion.brovey
i.gensig
i.gensigset
i.group
i.his.rgb
i.ifft
i.image.mosaic
i.in.spotvgt
i.landsat.acca
i.landsat.rgb
i.landsat.toar
i.maxlik
i.modis.qc
i.oif
i.ortho.photo
i.pca
i.photo.2image
i.photo.2target
Performs atmospheric correction using the 6S algorithm.

Canonical components analysis (cca) program for image processing.
Generates spectral signatures for an image by allowing the user to outline
regions of interest.
Generates spectral signatures for land cover types in an image using a clustering
algorithm.
Fast Fourier Transform (FFT) for image processing.
Brovey transform to merge multispectral and high-res panchromatic channels
Generates statistics for i.maxlik from raster map.
Generates statistics for i.smap from raster map.
Creates, edits, and lists groups and subgroups of imagery files.
Transforms raster maps from HIS (Hue-Intensity-Saturation) color space to
RGB (Red-Green-Blue) color space.
Inverse Fast Fourier Transform (IFFT) for image processing.
Mosaics up to 4 images and extends colormap; creates map *.mosaic
Import of SPOT VGT NDVI file into a raster map
Performs Landsat TM/ETM+ Automatic Cloud Cover Assessment (ACCA).
Performs auto-balancing of colors for LANDSAT images.
Calculates top-of-atmosphere radiance or reflectance and temperature for Landsat MSS/TM/ETM+.
Classifies the cell spectral reflectances in imagery data.
Extracts quality control parameters from Modis QC layers.
Calculates Optimum-Index-Factor table for LANDSAT TM bands 1-5, & 7
Menu driver for the photo imagery programs.
Principal components analysis (PCA) for image processing.
Interactively mark fiducial or reseau points on an image.
Creates control points on an image to be ortho-rectified.
i.photo.camera
i.photo.init
i.photo.rectify
i.points
i.rectify
i.rgb.his
i.smap
i.spectral
i.target
i.tasscap
i.topo.corr
i.vpoints
i.zc
119
Interactively select and modify the imagery group camera reference file.
Interactively creates or modifies entries in a camera initial exposure station file
for imagery group referenced by a sub-block.
Orthorectifies an image by using the image to photo coordinate transformation
matrix.
Mark ground control points on image to be rectified.
Rectifies an image by computing a coordinate transformation for each pixel in
the image based on the control points.
Transforms raster maps from RGB (Red-Green-Blue) color space to HIS (HueIntensity-Saturation) color space.
Performs contextual image classification using sequential maximum a posteriori
(SMAP) estimation.
Displays spectral response at user specified locations in group or images.
Targets an imagery group to a GRASS location and mapset.
Tasseled Cap (Kauth Thomas) transformation for LANDSAT-TM data
Computes topographic correction of reflectance.
Set registration points for an imagery group from a vector map or keyboard
entry.
Zero-crossing edge detectionraster function for image processing.
comandos m.* (miscellaneous)

m.cogo
m.nviz.image
m.proj
A simple utility for converting bearing and distance measurements to coordinates and vice versa.
Creates a 3D rendering of GIS data.
Convert coordinates from one projection to another (cs2cs frontend).
comando ps.* (postscript)

ps.map
Hardcopy PostScript map output utility.
comandos r.* (raster )

r.average
r.basins.fill
r.bilinear
r.bitpattern
r.blend
r.buffer
r.carve
r.category
Finds the average of values in a cover map within areas assigned the same
category value in a user-specified base map.
Generates watershed subbasins raster map.
Bilinear interpolation utility for raster map layers.
Compares bit patterns with a raster map.
Blends color components of two raster maps by a given ratio.
Creates a raster map layer showing buffer zones surrounding cells that contain
non-NULL category values.
Takes vector stream data, transforms it to raster and subtracts depth from the
output DEM.
Manages category values and labels associated with user-specified raster map
layers.
r.circle
r.clump
r.coin
r.colors
r.colors.out
r.colors.stddev
r.composite
r.compress
r.contour
r.cost
r.covar
r.cross
r.describe
r.digit
r.distance
r.drain
r.external
r.fill.dir
r.fillnulls
r.flow
r.grow.distance
r.grow
r.gwflow
r.his
r.horizon
r.in.arc
r.in.ascii
r.in.aster
r.in.bin
120
Creates a raster map containing concentric rings around a given point.

Recategorizes data in a raster map by grouping cells that form physically discrete areas into unique categories.
Tabulates the mutual occurrence (coincidence) of categories for two raster map
layers.
Creates/modifies the color table associated with a raster map layer.
Exports the color table associated with a raster map layer.
Set color rules based on stddev from a maps mean value.
Combines red, green and blue raster maps into a single composite raster map.
Compresses and decompresses raster maps.
Produces a vector map of specified contours from a raster map.
Creates a raster map showing the cumulative cost of moving between different geographic locations on an input raster map whose cell category values
represent cost.
Outputs a covariance/correlation matrix for user-specified raster map layer(s).
Creates a cross product of the category values from multiple raster map layers.
Prints terse list of category values found in a raster map layer.
Interactive tool used to draw and save vector features on a graphics monitor
using a pointing device (mouse) and save to a raster map.
Locates the closest points between objects in two raster maps.
Traces a flow through an elevation model on a raster map.
Link GDAL supported raster file to a binary raster map layer.
Filters and generates a depressionless elevation map and a flow direction map
from a given elevation raster map.
Fills no-data areas in raster maps using spline interpolation.
Construction of slope curves (flowlines), flowpath lengths, and flowline densities
(upslope areas) from a raster digital elevation model (DEM).
Generates a raster map layer of distance to features in input layer.
Generates a raster map layer with contiguous areas grown by one cell.
Numerical calculation program for transient, confined and unconfined
groundwater flow in two dimensions.
Generates red, green and blue raster map layers combining hue, intensity and
saturation (HIS) values from user-specified input raster map layers.
Horizon angle computation from a digital elevation model.
Converts an ESRI ARC/INFO ascii raster file (GRID) into a (binary) raster
map layer.
Converts ASCII raster file to binary raster map layer.
Georeference, rectify and import Terra-ASTER imagery and relative DEMs
using gdalwarp.
Import a binary raster file into a GRASS raster map.
r.info
r.in.gdal
r.in.gridatb
r.in.mat
r.in.poly
r.in.srtm
r.in.wms
r.in.xyz
r.kappa
r.lake
r.le.patch
r.le.pixel
r.le.setup
r.le.trace
r.li.cwed
r.li.daemon
r.li.dominance
r.li.edgedensity
r.li
r.li.mpa
r.li.mps
r.li.padcv
r.li.padrange
r.li.padsd
r.li.patchdensity
r.li.patchnum
r.li.pielou
r.li.renyi
r.li.richness
r.li.setup
r.li.shannon
r.li.shape
r.li.simpson
121
Output basic information about a raster map layer.

Import GDAL supported raster file into a binary raster map layer.
Imports GRIDATB.FOR map file (TOPMODEL) into GRASS raster map
Imports a binary MAT-File(v4) to a GRASS raster.
Creates raster maps from ASCII polygon/line/point data files.
Import SRTM HGT files into GRASS
Downloads and imports data from WMS servers.
Create a raster map from an assemblage of many coordinates using univariate
statistics.
Calculate error matrix and kappa parameter for accuracy assessment of classification result.
Fills lake at given point to given level.
Calculates attribute, patch size, core (interior) size, shape, fractal dimension,
and perimeter measures for sets of patches in a landscape.
Contains a set of measures for attributes, diversity, texture, juxtaposition, and
edge.
Interactive tool used to setup the sampling and analysis framework that will
be used by the other r.le programs.
Displays the boundary of each r.le patch and shows how the boundary is traced,
displays the attribute, size, perimeter and shape indices for each patch and saves
the data in an output file.
Calculates contrast weighted edge density index on a raster map
Support module for r.li landscape index calculations.
Calculates dominances diversity index on a raster map
Calculates edge density index on a raster map, using a 4 neighbour algorithm
Landscape structure analysis package overview.
Calculates mean pixel attribute index on a raster map
Calculates mean patch size index on a raster map, using a 4 neighbour algorithm
Calculates coefficient of variation of patch area on a raster map
Calculates range of patch area size on a raster map
Calculates standard deviation of patch area a raster map
Calculates patch density index on a raster map, using a 4 neighbour algorithm
Calculates patch number index on a raster map, using a 4 neighbour algorithm.
Calculates Pielous diversity index on a raster map
Calculates Renyis diversity index on a raster map
Calculates dominances diversity index on a raster map
Configuration editor for r.li.index
Calculates Shannons diversity index on a raster map
Calculates shape index on a raster map
Calculates Simpsons diversity index on a raster map
r.los
r.mapcalc
r.mapcalculator
r.mask
r.median
r.mfilter.fp
r.mfilter
r.mode
r.neighbors
r.null
r.out.arc
r.out.ascii
r.out.bin
r.out.gdal
r.out.gdal.sh
r.out.gridatb
r.out.mat
r.out.mpeg
r.out.png
r.out.pov
r.out.ppm3
r.out.ppm
r.out.tiff
r.out.vrml
r.out.vtk
r.out.xyz
r.pack
r.param.scale
r.patch
r.plane
r.profile
r.proj
122
Line-of-sight raster analysis program.

Raster map calculator.
Calculate new raster map from a r.mapcalc expression.
Creates a MASK for limiting raster operation.
Finds the median of values in a cover map within areas assigned the same
category value in a user-specified base map.
Raster map matrix filter.
Performs raster map matrix filter.
Finds the mode of values in a cover map within areas assigned the same category
value in a user-specified base map.
Makes each cell category value a function of the category values assigned to
the cells around it, and stores new cell values in an output raster map layer.
Manages NULL-values of given raster map.
Converts a raster map layer into an ESRI ARCGRID file.
Converts a raster map layer into an ASCII text file.
Exports a GRASS raster map to a binary array.
Exports GRASS raster maps into GDAL supported formats.
Exports GRASS raster into GDAL supported formats.
Exports GRASS raster map to GRIDATB.FOR map file (TOPMODEL)
Exports a GRASS raster to a binary MAT-File.
Raster File Series to MPEG Conversion.
Export GRASS raster as non-georeferenced PNG image.
Converts a raster map layer into a height-field file for POVRAY.
Converts 3 GRASS raster layers (R,G,B) to a PPM image file at the pixel
resolution of the CURRENTLY DEFINED REGION.
Converts a GRASS raster map to a PPM image file at the pixel resolution of
the currently defined region.
Exports a GRASS raster map to a 8/24bit TIFF image file at the pixel resolution of the currently defined region.
Export a raster map to the Virtual Reality Modeling Language (VRML)
Converts raster maps into the VTK-Ascii format
Export a raster map to a text file as x,y,z values based on cell centers.
Packs up a raster map and support files for copying.
Extracts terrain parameters from a DEM.
Creates a composite raster map layer by using known category values from one
(or more) map layer(s) to fill in areas of no datain another map layer.
Creates raster plane map given dip (inclination), aspect (azimuth) and one
point.
Outputs the raster map layer values lying on user-defined line(s).
Re-projects a raster map from one location to the current location.
r.quant
r.quantile
r.random.cells
r.random
r.random.surface
r.reclass.area
r.reclass
r.recode
r.region
r.regression.line
r.report
r.resamp.interp
r.resample
r.resamp.rst
r.resamp.stats
r.rescale.eq
r.rescale
r.ros
r.series
r.shaded.relief
r.sim.sediment
r.sim.water
r.slope.aspect
r.spread
r.spreadpath
r.statistics
r.stats
r.sum
123
Produces the quantization file for a floating-point map.

Compute quantiles using two passes.
Generates random cell values with spatial dependence.
Creates a raster map layer and vector point map containing randomly located
points.
Generates random surface(s) with spatial dependence.
Reclasses a raster map greater or less than user specified area size (in hectares).
Creates a new map layer whose category values are based upon a reclassification
of the categories in an existing raster map layer.
Recodes categorical raster maps.
Sets the boundary definitions for a raster map.
Calculates linear regression from two raster maps: y = a + b*x.
Reports statistics for raster map layers.
Resamples raster map layers to a finer grid using interpolation.
GRASS raster map layer data resampling capability.
Reinterpolates and optionally computes topographic analysis from input raster
map to a new raster map (possibly with different resolution) using regularized
spline with tension and smoothing.
Resamples raster map layers to a coarser grid using aggregation.
Rescales histogram equalized the range of category values in a raster map layer.
Rescales the range of category values in a raster map layer.
Generates three, or four raster map layers showing 1) the base (perpendicular)
rate of spread (ROS), 2) the maximum (forward) ROS, 3) the direction of the
maximum ROS, and optionally 4) the maximum potential spotting distance.
Makes each output cell value a function of the values assigned to the corresponding cells in the input raster map layers.
Creates shaded relief map from an elevation map (DEM).
Sediment transport and erosion/deposition simulation using path sampling
method (SIMWE).
Overland flow hydrologic simulation using path sampling method (SIMWE).
Generates raster maps of slope, aspect, curvatures and partial derivatives from
a elevation raster map.
Simulates elliptically anisotropic spread on a graphics window and generates a
raster map of the cumulative time of spread, given raster maps containing the
rates of spread (ROS), the ROS directions and the spread origins.
Recursively traces the least cost path backwards to cells from which the cumulative cost was determined.
Calculates category or object oriented statistics.
Generates area statistics for raster map layers.
Sums up the raster cell values.
r.sun
r.sunmask
r.support
r.support.stats
r.surf.area
r.surf.contour
r.surf.fractal
r.surf.gauss
r.surf.idw2
r.surf.idw
r.surf.random
r.terraflow
r.texture
r.thin
r.tileset
r.timestamp
r.topidx
r.topmodel
r.to.rast3elev
r.to.rast3
r.to.vect
r.transect
r.univar
r.univar.sh
r.unpack
r.volume
r.walk
r.water.outlet
r.watershed
r.what.color
r.what
124
Solar irradiance and irradiation model.

Calculates cast shadow areas from sun position and elevation raster map.
Allows creation and/or modification of raster map layer support files.
Update raster map statistics
Surface area estimation for rasters.
Surface generation program from rasterized contours.
Creates a fractal surface of a given fractal dimension.
GRASS module to produce a raster map layer of gaussian deviates whose mean
and standard deviation can be expressed by the user. It uses a gaussian random
number generator.
Surface generation program.
Surface interpolation utility for raster map.
Produces a raster map of uniform random deviates whose range can be expressed by the user.
Flow computation for massive grids (float version).
Generate images with textural features from a raster map.
Thins non-zero cells that denote linear features in a raster map.
Produces tilings of the source projection for use in the destination region and
projection.
Modifies a timestamp for a raster map.
Creates topographic index map from elevation raster map.
Simulates TOPMODEL which is a physically based hydrologic model.
Creates a 3D volume map based on 2D elevation and value raster maps.
Converts 2D raster map slices to one 3D raster volume map.
Converts a raster map into a vector map layer.
Outputs raster map layer values lying along user defined transect line(s).
Calculates univariate statistics from the non-null cells of a raster map.
calculates univariate statistics from a GRASS raster map
Unpacks a raster map packed with r.pack.
Calculates the volume of data clumps, and (optionally) produces a GRASS
vector points map containing the calculated centroids of these clumps.
Outputs a raster map layer showing the anisotropic cumulative cost of moving
between different geographic locations on an input elevation raster map layer
whose cell category values represent elevation combined with an input raster
map layer whose cell values represent friction cost.
Watershed basin creation program.
Watershed basin analysis program.
Queries colors for a raster map layer.
Queries raster map layers on their category values and category labels.
comandos r3.* (raster 3D)

r3.cross.rast
Creates cross section 2D raster map from 3D raster map based on 2D elevation
map.
r3.gwflow
r3.in.ascii
r3.info
r3.in.v5d
r3.mapcalc
r3.mapcalculator
r3.mask
r3.mkdspf
r3.null
r3.out.ascii
r3.out.v5d
r3.out.vtk
r3.stats
r3.timestamp
r3.to.rast
r3.univar
125
Calculates numerically transient, confined groundwater flow in three dimensions.

Converts a 3D ASCII raster text file into a (binary) 3D raster map layer.
Outputs basic information about a user-specified 3D raster map layer.
Imports 3-dimensional Vis5D files (i.e. the V5D file with 1 variable and 1 time
step).
3D raster map calculator.
Calculates new grid3D volume from r3.mapcalc expression.
Establishes the current working 3D raster mask.
Creates a display file from an existing grid3 file according to specified threshold
levels.
Explicitly create the 3D NULL-value bitmap file.
Converts a 3D raster map layer into an ASCII text file.
Exports GRASS 3D raster map to 3-dimensional Vis5D file.
Converts 3D raster maps into the VTK-ASCII format.
Generates volume statistics for 3D raster maps.
Print/add/remove a timestamp for a 3D raster map
Converts 3D raster maps to 2D raster maps
Calculates univariate statistics from the non-null 3d cells of a raster3d map.
comandos v.* (vector )

v.buffer
v.build.all
v.build
v.build.polylines
v.category
v.centroids
v.class
v.clean
v.colors
v.convert.all
v.convert
v.db.addcol
v.db.addtable
v.db.connect
v.db.dropcol
v.db.droptable
v.db.join
v.db.reconnect.all
Creates a buffer around features of given type (areas must contain centroid).
Rebuilds topology on all vector maps in the current mapset.
Creates topology for GRASS vector map.
Builds polylines from lines or boundaries.
Attach, delete or report vector categories to map geometry.
Adds missing centroids to closed boundaries.
Classifies attribute data, e.g. for thematic mapping
Toolset for cleaning topology of vector map.
Set color rules for features in a vector using a numeric attribute column.
Converts all older versions of GRASS vector maps in current mapset to current
format.
Imports older versions of GRASS vector maps.
Adds one or more columns to the attribute table connected to a given vector
map.
Creates and connects a new attribute table to a given layer of an existing vector
map.
Prints/sets DB connection for a vector map to attribute table.
Drops a column from the attribute table connected to a given vector map.
Removes existing attribute table of a vector map.
Allows to join a table to a vector map table.
Reconnects vectors to a new database.
v.db.renamecol
v.db.select
v.db.univar
v.db.update
v.delaunay
v.digit
v.dissolve
v.distance
v.drape
v.edit
v.external
v.extract
v.extrude
v.generalize
v.hull
v.in.ascii
v.in.db
v.in.dxf
v.in.e00
v.info
v.in.garmin
v.in.geonames
v.in.gns
v.in.gpsbabel
v.in.lines
v.in.mapgen
v.in.ogr
v.in.region
v.in.sites.all
v.in.sites
v.in.wfs
v.kcv
126
Renames a column in the attribute table connected to a given vector map.

Prints vector map attributes.
Calculates univariate statistics on selected table column for a GRASS vector
map.
Allows to update a column in the attribute table connected to a vector map.
Creates a Delaunay triangulation from an input vector map containing points
or centroids.
Interactive editing and digitization of vector maps.
Dissolves boundaries between adjacent areas sharing a common category number or attribute.
Finds the nearest element in vector map to for elements in vector map from.
Converts vector map to 3D by sampling of elevation raster map.
Edits a vector map, allows adding, deleting and modifying selected vector features.
Creates a new vector as a read-only link to OGR layer.
Selects vector objects from an existing vector map and creates a new map
containing only the selected objects.
Extrudes flat vector object to 3D with defined height.
Vector based generalization.
Produces a convex hull for a given vector map.
Creates a vector map from ASCII points file or ASCII vector file.
Creates new vector (points) map from database table containing coordinates.
Converts files in DXF format to GRASS vector map format.
Import E00 file into a vector map.
Outputs basic information about a user-specified vector map.
Download waypoints, routes, and tracks from a Garmin GPS receiver into a
vector map.
Imports geonames.org country files into a GRASS vector points map.
Imports US-NGA GEOnet Names Server (GNS) country files into a GRASS
vector points map.
Import waypoints, routes, and tracks from a GPS receiver or GPS download
file into a vector map.
Import ASCII x,y[,z] coordinates as a series of lines.
Imports Mapgen or Matlab-ASCII vector maps into GRASS.
Convert OGR vector layers to GRASS vector map.
Creates a vector polygon from the current region extent.
Converts all old GRASS < Ver5.7 sites maps in current mapset to vector maps.
Converts a GRASS site lists file into a vector map.
Import GetFeature from WFS
Randomly partition points into test/train sets.
v.kernel
127
Generates a raster density map from vector point data using a moving kernel
or optionally generates a vector density map on a vector network.
v.label
Creates paint labels for a vector map from attached attributes.
v.label.sa
Create optimally placed labels for vector map(s)
v.lidar.correction
Correction of the v.lidar.growing output. It is the last of the three algorithms
for LIDAR filtering.
v.lidar.edgedetection Detects the objects edges from a LIDAR data set.
v.lidar.growing
Building contour determination and Region Growing algorithm for determining
the building inside
v.lrs.create
Create Linear Reference System
v.lrs.label
Create stationing from input lines, and linear reference system
v.lrs.segment
Creates points/segments from input lines, linear reference system and positions
read from stdin or a file.
v.lrs.where
Finds line id and real km+offset for given points in vector map using linear
reference system.
v.mkgrid
Creates a GRASS vector map of a user-defined grid.
v.neighbors
Makes each cell value a function of the attribute values assigned to the vector
points or centroids around it, and stores new cell values in an output raster
map layer.
v.net.alloc
Allocate subnets for nearest centres (direction from centre).
v.net.allpairs
Computes the shortest path between all pairs of nodes in the network.
v.net.bridge
Computes bridges and articulation points in the network.
v.net.centrality
Computes degree, centrality, betweeness, closeness and eigenvector centrality
measures in the network.
v.net.components
Computes strongly and weakly connected components in the network.
v.net.connectivity
Computes vertex connectivity between two sets of nodes in the network.
v.net.distance
Computes shortest distance via the network between the given sets of features.
v.net.flow
Computes the maximum flow between two sets of nodes in the network.
v.net
Performs network maintenance.
v.net.iso
Splits net by cost isolines.
v.net.path
Finds shortest path on vector network.
v.net.salesman
Creates a cycle connecting given nodes (Traveling salesman problem).
v.net.spanningtree
Computes minimum spanning tree for the network.
v.net.steiner
Create Steiner tree for the network and given terminals
v.net.timetable
Finds shortest path using timetables.
v.net.visibility
Visibility graph construction.
v.normal
Tests for normality for points.
v.out.ascii
Converts a GRASS binary vector map to a GRASS ASCII vector map.
v.out.dxf
Exports GRASS vector map layers to DXF file format.
v.out.gpsbabel
Exports a vector map to a GPS receiver or file format supported by GpsBabel.
v.outlier
v.out.ogr
v.out.pov
v.out.svg
v.out.vtk
v.overlay
v.parallel
v.patch
v.perturb
v.proj
v.qcount
v.random
v.rast.stats
v.reclass
v.report
v.sample
v.segment
v.select
v.split
v.support
v.surf.bspline
v.surf.idw
v.surf.rst
v.to.3d
v.to.db
v.to.points
v.to.rast3
v.to.rast
v.transform
v.type
v.univar
128
Removes outliers from vector point data.

Converts GRASS vector map to one of the supported OGR vector formats.
Converts to POV-Ray format, GRASS x,y,z -> POV-Ray x,z,y
Exports a GRASS vector map to SVG.
Converts a GRASS binary vector map to VTK ASCII output.
Overlays two vector maps.
Creates parallel line to input vector lines.
Create a new vector map layer by combining other vector map layers.
Random location perturbations of GRASS vector points
Re-projects a vector map from one location to the current location.
Indices for quadrat counts of sites lists.
Randomly generate a 2D/3D vector points map.
Calculates univariate statistics from a raster map based on vector polygons and
uploads statistics to new attribute columns.
Changes vector category values for an existing vector map according to results
of SQL queries or a value in attribute table column.
Reports geometry statistics for vectors.
Samples a raster map at vector point locations.
Creates points/segments from input vector lines and positions.
Selects features from vector map (A) by features from other vector map (B).
Split lines to shorter segments.
Updates vector map metadata.
Bicubic or bilinear spline interpolation with Tykhonov regularization.
Surface interpolation from vector point data by Inverse Distance Squared
Weighting.
Spatial approximation and topographic analysis from given point or isoline data
in vector format to floating point raster format using regularized spline with
tension.
Performs transformation of 2D vector features to 3D.
Populates database values from vector features.
Create points along input lines in new vector with 2 layers.
Converts a binary GRASS vector map (only points) layer into a 3D GRASS
raster map layer.
Converts (rasterize) a vector map into a raster map.
Performs an affine transformation (shift, scale and rotate, or GPCs) on vector
map.
Change the type of geometry elements.
Calculates univariate statistics for attribute. Variance and standard deviation
is calculated only for points if specified.
v.univar.sh
v.vol.rst
v.voronoi
v.what
v.what.rast
v.what.vect
129
Calculates univariate statistics on selected table column for a GRASS vector

map.
Interpolates point data to a 3D raster map using regularized spline with tension
(RST) algorithm.
Creates a Voronoi diagram from an input vector map containing points or
centroids.
Queries a vector map layer at given locations.
Uploads raster values at positions of vector points to the table.
Uploads vector values at positions of vector points to the table.
16 - Referencias Bibliogr
aficas
16.
130
Refer
encias Bibliogr
aficas
Bamler, R. & Hartl, P., 1998. Synthetic aperture radar interferometry. Inverse Problems, 14:R1R54.
Becker, J. J. & Sandwell, D. T., 2007. SRTM30 PLUS: Data fusion of SRTM land topography with measured
and estimated seafloor topography (Version 3.0). http://topex.ucsd.edu/WWW_html/srtm30_plus.html,
[accessed 09/04/2008].
Dassau, O., Holl, S., Neteler, M., & Redslob, M., 2005. An introduction to the practical use of the Free
Geographical Information System GRASS 6.0. GDF Hannover bR.
Evans, I. S., 1980. An integrated system of terrain analysis and slope mapping. Zeitschrift f
ur Geomorphologie,
Suppl.-Bd. 36:274295.
Farr, T. G. & Kobrick, M., 2000. Shuttle Radar Topography Mission produces a wealth of data. EOS
(Transactions, American Geophysical Union), 81:583585.
Farr, T. G., Rosen, P. A., Caro, E., Crippen, R., Duren, R., Hensley, S., Kobrick, M., Paller, M., Rodriguez,
E., Roth, L., Seal, D., Shaffer, S., Shimada, J., Umland, J., Werner, M., Oskin, M., Burbank, D., & Alsdorf,
D., 2007. The Shuttle Radar Topography Mission. Review of Geophysics, 45:RG2004.
Florinsky, I. V., 1998. Accuracy of local topographic variables derived from digital elevation models. International Journal of Geographical Information Science, 12:4761.
Florinsky, I. V., Eilers, R. G., Manning, G., & Fuller, L. G., 2002. Prediction of soil properties by digital
terrain modelling. Environmental Modelling & Software, 17:295311.
Grohmann, C. H., 2008. Introduca
o ao geoprocessamento e a
` an
alise digital de terreno com software livre.
Technical report, Instituto de Geociencias, Universidade de S
ao Paulo, S
ao Paulo.
Hofierka, J., Mitasova, H., Parajka, J., & Mitas, L., 2002. Multivariate interpolation of precipitation using
Regularized Spline with Tension. Transactions in GIS, 6:135150.
Horn, B. K. P., 1981. Hill Shading and the Reflectance Map. Proceedings of the IEEE, 69:1447.
IBGE, 1998. Noco
es b
asicas de Cartografia. Diretoria de Geociencias. Departamento de Cartografia.
Jakobsson, M., Cherkis, N. Z., Woodward, J., Macnab, R., & Coakley, B., 2000. New grid of Arctic bathymetry
aids scientists and mapmakers. EOS (Transactions, American Geophysical Union), 81:8996.
Lillesand, T. M. & Kiefer, R. W., 1987. Remote sensing and image interpretation. John Wiley & Sons (2nd
ed). New York.
Miranda, E. E., 2005. Brasil em Relevo. Embrapa Monitoramento por Satelite, Campinas. Documento eletr
onico. http://www.relevobr.cnpm.embrapa.br.
Mitasova, H. & Hofierka, J., 1993. Interpolation by regularized spline with tension:II. Application to terrain
modeling and surface geometry analysis. Mathematical Geology, 25:657669.
Mitasova, H. & Mitas, L., 1993. Interpolation by regularized spline with tension:I Theory and implementation.
Mathematical Geology, 25:641655.
Nanni, A. S., Descovi Filho, L., Virtuoso, M. A., Montenegro, D., Willrich, G., Machado, P. H., Sperb, R.,
Dantas, G. S., & Calanzans, Y., 2012. Quantum GIS - Guia do Usu
ario, Vers
ao 1.7.4 Wroclaw.
16 - Referencias Bibliogr
aficas
131
Neteler, M. & Mitasova, H., 2008. Open Source GIS : A GRASS GIS Approach, Third Edition (The International Series in Engineering and Computer Science), volume 773. Springer, New York.
Poppe, L. J., Williams, S. J., & Paskevich, V. F., 2005. USGS east-coast sediment analysis: procedures,
database, and GIS data. US Geological Survey Open-File Report, 1001:00358.
Rabus, B., Eineder, M., Roth, R., & Bamler, R., 2003. The Shuttle Radar Topography Missiona new class
of digital elevation models acquired by spaceborne radar. ISPRS Journal of Photogrammetry & Remote
Sensing, 57:241262.
Rodriguez, E., Morris, C. H., & Belz, J. E., 2006. A global assessment of the SRTM performance. Photogrammetric Engineering and Remote Sensing, 72:249260.
Ruhe, R. V., 1975. Geomorphology. Houghton Mifflin, Boston.
Shapiro, M. & Westervelt, J., 1991. r.mapcalc: An Algebra for GIS and Image Processing. Technical report,
U.S. Army Construction Engineering Research Laboratory.
Shary, P. A., Sharaya, L. S., & Mitusov, A. V., 2002. Fundamental quantitative methods of land surface
analysis. Geoderma, 107:132.
Smith, W. H. F. & Sandwell, D. T., 1997. Global seafloor topography from satellite altimetry and ship depth
soundings. Science, 277:19571962.
Smith, W. H. F. & Wessel, P., 1990. Gridding with continuous curvatures splines in tension. Geophysics,
55:293305.
Steiner, S. S., 2007. Aquisic
ao e processamento de dados morfometricos derivados do modelo digital de elevaca
o
SRTM. Masters thesis, Instituto de Geociencias, Universidade de S
ao Paulo, S
ao Paulo.
van Zyl, J. J., 2001. The Shuttle Radar Topography Mission (SRTM): a breakthrough in remote sensing of
topography. Acta Astronautica, 48:559565.
Wood, J., 1996. The Geomorphological characterization of Digital Elevation Models. PhD thesis, University of
Leicester, U. K. , Department of Geography, PhD Thesis.
Zebker, H. A. & Goldstein, R. M., 1986. Topographic mapping from interferometric synthetic aperture radar
observations. Journal of Geophysical Research, 91:49934999.
Zevenbergen, L. W. & Thorne, C. R., 1987. Quantitative analysis of land surface topography. Earth Surface
Processes and Landforms, 12:4756.
17 - Definic
17.
132
Definic
N
os mantemos esta definic
ao do Software Livre para mostrar claramente o que deve ser verdadeiro `
a
respeito de um dado programa de software para que ele seja considerado software livre.
Software Livre e uma quest
ao de liberdade, n
ao de preco. Para entender o conceito, voce deve pensar
em liberdade de express
ao, n
ao em cerveja gr
atis.
Software livrese refere a
` liberdade dos usu
arios executarem, copiarem, distriburem, estudarem, modificarem e aperfeicoarem o software. Mais precisamente, ele se refere a quatro tipos de liberdade, para os
usu
arios do software:
A liberdade de executar o programa, para qualquer prop

osito (liberdade n 0)
A liberdade de estudar como o programa funciona, e adapt
a-lo para as suas necessidades (liberdade n
1). Acesso ao c
odigo-fonte e um pre-requisito para esta liberdade.
A liberdade de redistribuir c
opias de modo que voce possa ajudar ao seu pr
oximo (liberdade n 2).
A liberdade de aperfeicoar o programa, e liberar os seus aperfeicoamentos, de modo que toda a comunidade se beneficie (liberdade n 3). Acesso ao c
odigo-fonte e um pre-requisito para esta liberdade.
Um programa e software livre se os usu
arios tem todas essas liberdades. Portanto, voce deve ser livre
para redistribuir c
opias, seja com ou sem modificaco
es, seja de graca ou cobrando uma taxa pela distribuica
o,
para qualquer um em qualquer lugar. Ser livre para fazer essas coisas significa (entre outras coisas) que voce
n
ao tem que pedir ou pagar pela permiss
ao.
Voce deve tambem ter a liberdade de fazer modifcaco
es e us
a-las privativamente no seu trabalho ou lazer,
sem nem mesmo mencionar que elas existem. Se voce publicar as modificaco
es, voce n
ao deve ser obrigado a
avisar a ninguem em particular, ou de nenhum modo em especial.
A liberdade de utilizar um programa significa a liberdade para qualquer tipo de pessoa fsica ou jurdica
utilizar o software em qualquer tipo de sistema computacional, para qualquer tipo de trabalho ou atividade,
sem que seja necess
ario comunicar ao desenvolvedor ou a qualquer outra entidade em especial.
A liberdade de redistribuir c
opias deve incluir formas bin
arias ou execut
aveis do programa, assim como o
c
odigo-fonte, tanto para as vers
oes originais quanto para as modificadas. Est
a ok se n
ao for possvel produzir
uma forma bin
aria ou execut
avel (pois algumas linguagens de programaca
o n
ao suportam este recurso), mas
deve ser concedida a liberdade de redistribuir essas formas caso seja desenvolvido um meio de cria-las.
De modo que a liberdade de fazer modificaco
es, e de publicar vers
oes aperfeicoadas, tenha algum significado, deve-se ter acesso ao c
odigo-fonte do programa. Portanto, acesso ao c
odigo-fonte e uma condica
o
necess
aria ao software livre.
Para que essas liberdades sejam reais, elas tem que ser irrevog
aveis desde que voce n
ao faca nada errado;
caso o desenvolvedor do software tenha o poder de revogar a licenca, mesmo que voce n
ao tenha dado motivo,
o software n
ao e livre.
Entretanto, certos tipos de regras sobre a maneira de distribuir software livre s
ao aceit
aveis, quando
elas n
ao entram em conflito com as liberdades principais. Por exemplo, copyleft (apresentado de forma bem
17 - Definic
133
simples) e a regra de que, quando redistribuindo um programa, voce n

ao pode adicionar restrico
es para negar
para outras pessoas as liberdades principais. Esta regra n
ao entra em conflito com as liberdades; na verdade,
ela as protege.
Portanto, voce pode ter pago para receber c
opias do software GNU, ou voce pode ter obtido c
opias sem
nenhum custo. Mas independente de como voce obteve a sua c
opia, voce sempre tem a liberdade de copiar e
modificar o software, ou mesmo de vender c
opias.
Software Livre n
ao significa n
ao-comercial. Um programa livre deve estar disponvel para uso
comercial, desenvolvimento comercial, e distribuica
o comercial. O desenvolvimento comercial de software livre
n
ao e incomum; tais softwares livres comerciais s
ao muito importantes.
Regras sobre como empacotar uma vers
ao modificada s
ao aceit
aveis, se elas n
ao acabam bloqueando a
sua liberdade de liberar vers
oes modificadas. Regras como se voce tornou o programa disponvel deste modo,
voce tambem tem que torn
a-lo disponvel deste outro modotambem podem ser aceitas, da mesma forma.
(Note que tal regra ainda deixa para voce a escolha de tornar o programa disponvel ou n
ao.) Tambem e
aceit
avel uma licenca que exija que, caso voce tenha distribudo uma vers
ao modificada e um desenvolvedor
anterior peca por uma c
opia dele, voce deva enviar uma.
No projeto GNU, n
os usamos copyleft35 para proteger estas liberdades legalmente para todos. Mas
tambem existe software livre que n
ao e copyleft. N
os acreditamos que hajam raz
oes importantes pelas quais e
36
melhor usar o copyleft, mas se o seu programa e free-software mas n
ao e copyleft, n
os ainda podemos utiliz
a-lo.
Veja Categorias de Software Livre37 para uma descrica
o de como software livre, software copylefte
outras categoria se relacionam umas com as outras.
` vezes regras de controle de exportaca
As
o e sans
oes de comercio podem limitar a sua liberdade de distribuir
c
opias de programas internacionalmente. Desenvolvedores de software n
ao tem o poder para eliminar ou
sobrepor estas restric
oes, mas o que eles podem e devem fazer e se recusar a imp
o-las como condico
es para o
uso dos seus programas. Deste modo, as restrico
es n
ao afetam as atividades e as pessoas fora da jurisdica
o
destes governos.
Quando falando sobre o software livre, e melhor evitar o uso de termos como dado ou de graca,
porque estes termos implicam que a quest
ao e de preco, n
ao de liberdade. Alguns temos comuns como
piratariaenglobam opini
oes que n
os esperamos voce n
ao ir
a endossar. Veja frases e palavras confusas38 que
e melhor evitar para uma discuss
ao desses termos. N
os tambem temos uma lista de traduco
es39 do termo
software livre para v
arias lnguas.
Finalmente, note que criterios como os estabelecidos nesta definica
o do software livre requerem cuidadosa
deliberaca
o quanto a
` sua interpretaca
o. Para decidir se uma licenca se qualifica como de software livre,
n
os a julgamos baseados nestes criterios para determinar se ela se segue o nosso esprito assim como as
palavras exatas. Se uma licenca inclui restrico
es impensadas, n
os a rejeitamos, mesmo que n
os n
ao tenhamos
` vezes um requerimento de alguma licenca levanta uma quest
antecipado a quest
ao nestes criterios. As
ao
35
http://www.gnu.org/copyleft/copyleft.pt.html
http://www.gnu.org/philosophy/pragmatic.html
37
http://www.gnu.org/philosophy/categories.html
38
http://www.gnu.org/philosophy/words-to-avoid.pt.html
39
http://www.gnu.org/philosophy/fs-translations.html
36
17 - Definic
134
que requer excessiva deliberaca

o, incluindo discuss
oes com advogados, antes que n
os possamos decidir se o
requerimento e aceit
avel. Quando n
os chegamos a uma conclus
ao sobre uma nova quest
ao, n
os frequentemente atualizamos estes criterios para tornar mais f
acil determinar porque certas licencas se qualificam ou n
ao.
Se voce est
a interessado em saber se uma licenca em especial se qualifica como uma licenca de software
livre, veja a nossa lista de licencas40 . Se a licenca com a qual voce est
a preocupado n
ao est
a listada, voce pode
nos questionar enviando e-mail para [email protected].
Outro grupo iniciou o uso do termo software aberto41 para significar algo pr
oximo (mas n
ao identico) a
software livre. N
os preferimos o termo software livre porque, uma vez que voce tenha aprendido que ele se
refere a
` liberdade e n
ao ao preco, voce se preocupar
a com a quest
ao da liberdade.
Por favor envie d
uvidas ou quest
oes sobre FSF e/ou GNU para [email protected].
c
Copyright 1996,
1997, 1998, 1999, 2000 Free Software Foundation, Inc., 51 Franklin St, Fifth Floor,
Boston, MA 02110, USA
A c
opia fiel e a distribuica
o deste artigo completo e permitida em qualquer meio, desde que esta nota seja
preservada.
Atualizado: 2006/05/20 09:22:22 Author: hicham
Traduzido por: Fernando Lozano [email protected]
40
41
http://www.gnu.org/licenses/license-list.html
http://www.gnu.org/philosophy/free-software-for-freedom.html
18 - GNU Free Documentation License
18.
135
GNU Free Documentation License
Esta e uma traduca

o n
ao oficial da Licenca de Documentaca
o Livre GNU em Portugues do Brasil. Ela
n
ao e publicada pela Free Software Foundation, e n
ao se aplica legalmente a
` distribuica
o de documentos que
usem a GNU LDL apenas o texto original em Ingles da GNU FDL se aplica. Contudo, esperamos que esta
traduca
o ajude os falantes desta lngua a entender melhor a GNU GFL.
This is an unofficial translation of the GNU Free Documentation License into Brazilian Portuguese. It
was not published by the Free Software Foundation, and does not legally state the distribution terms for
documentation that uses the GNU FDL only the original English text of the GNU FDL does that. However,
we hope that this translation will help Brazilian Portuguese speakers understand the GNU FDL better.
NOTA DO TRADUTOR
Traduca
o: Norton T. Roman ([email protected])
Revis
ao: Jo
ao S. O. Bueno Calligaris ([email protected])
Ultima Atualizac
ao: 01 de Maio de 2005.
C
opias exatas e distribuica
o deste documento s
ao permitidas em qualquer meio desde que a nota de
direitos autorais (copyright) e esta nota sejam preservadas.
136
Licen
ca de Documentac
ao Livre GNU
Vers
ao 1.2, Novembro de 2002
c 2000, 2001, 2002 Free Software Foundation, Inc. 59 Temple Place, Suite 330, Boston, MA
Copyright
permitido a qualquer um copiar e distribuir c
02111-1307 USA. E
opias exatas deste documento de licenca,
embora n
ao seja permitido alter
a-lo.
0. PREAMBULO
O prop
osito desta Licenca e fazer com que um manual, livro-texto, ou outro documento funcional e u
til
seja livre, garantindo a todos a liberdade efetiva de copi
a-lo e redistribui-lo, com ou sem modificaco
es, tanto
comercialmente como n
ao comercialmente. Em segundo lugar, esta Licenca fornece ao autor e ao editor um meio
de obter credito pelo seu trabalho, n
ao sendo, ao mesmo tempo, considerados respons
aveis por modificaco
es
feitas por outros.
Esta licenca e um tipo de esquerdo de c
opia (copyleft), o que significa que trabalhos derivados do
documentos devem, por sua vez, ser livres no mesmo sentido. Ela complementa a Licenca P
ublica Geral GNU,
a qual e uma licenca de esquerdo de c
opia criada para programas livres.
Criamos esta Licenca para que seja usada em manuais para programas livres, porque programas livres precisam de documentaca
o livre: um programa livre deveria vir com manuais que oferecam as mesmas liberdades
que o programa oferece. Mas esta Licenca n
ao est
a limitada a manuais de programas de computador; ela pode
ser usada para qualquer trabalho de texto, independentemente do assunto ou se e publicado como um livro impresso. N
os recomendamos esta Licenca principalmente para trabalhos cujo prop
osito e instruca
o ou referencia.
1. APLICABILIDADE E DEFINIC
OES
Esta licenca se aplica a qualquer manual ou outro trabalho, em qualquer meio, que contenha uma nota
introduzida pelo detentor dos direitos autorais dizendo que o documento pode ser distribudo sob os termos
desta. Tal nota garante uma licenca mundial, livre de royalties, de duraca
o ilimitada, para usar este trabalho
sob as condico
es aqui colocadas. O Documento, abaixo, se refere a qualquer tal manual ou trabalho. Qualquer
membro do p
ublico e um licenciado, e ser
a tratado por voce. Voce aceita a licenca se copiar, modificar ou
distribuir o trabalho de um modo que necessite de permiss
ao de acordo com a lei de direitos autorais.
Uma Vers
ao Modificada do Documento se refere a qualquer trabalho contendo o Documento ou uma
parte deste, quer seja copiado sem modificaco
es, quer com modificaco
es e/ou traduzido para outra lngua.
Uma Seca
o Secund
aria e um apendice com nome ou uma seca
o inicial do Documento que trata exclusivamente da relaca
o dos editores ou autores do Documento com seu assunto geral (ou temas relacionados) e
n
ao contem nada que possa estar diretamente dentro do assunto geral. Assim, se o Documento e em parte um
livro-texto de matem
atica, uma Seca
o Secund
aria n
ao pode explicar nada de matem
atica. Tal relaca
o pode
ser uma conex
ao hist
orica com o assunto ou com temas relacionados, ou tratar de quest
oes legais, comerciais,
filos
oficas, eticas ou polticas com relaca
o a eles.
Seco
es Invariantes s
ao certas Seco
es Secund
arias cujos ttulos s
ao designados como sendo de Seco
es
invariantes na nota que afirma que o Documento e publicado sob esta Licenca. Se uma seca
o n
ao se encaixa na
definica
o acima de Secund
aria, ent
ao n
ao se permite que seja designada como Invariante. O Documento pode
n
ao conter nenhuma Seca
o Invariante. Se o documento n
ao identificar quaisquer Seco
es Invariantes, ent
ao n
ao
h
a nenhuma.
Textos de Capa s
ao certas passagens de texto que s
ao listada como Textos de Capa Frontal ou Texto de
Quarta Capa, na nota que afirma que o Documento e publicado sob esta Licenca. Um Texto de Capa Frontal
pode ter no m
aximo 5 palavras, e um Texto de Quarta Capa pode ter no m
aximo 25 palavras.
Uma c
opia Transparente do Documento significa uma c
opia que pode ser lida pelo computador, representada em um formato cuja especificaca
o esteja disponvel ao p
ublico geral, que seja apropriada para a imediata
revis
ao do documento usando-se editores de texto genericos ou (para imagens compostas de pixeis) programas
137
gr
aficos genericos ou (para desenhos) algum editor de desenhos amplamente disponvel, e que seja apropriado
para inclus
ao em formatadores de texto ou para traduc
ao autom
atica para uma variedade de formatos apropriados para inclus
ao em formatadores de texto. Uma c
opia feita em outro formato de arquivo Transparente
cuja marcaca
o, ou ausencia desta, foi manipulada para impedir ou desencorajar modificaca
o subseq
uente pelos leitores n
ao e Transparente. Um formato de imagem n
ao e Transparente se usado em lugar de qualquer
quantidade substancial de texto. Uma c
opia que n
ao e Transparente e chamada Opaca.
Exemplos de formatos apropriados para c
opias Transparentes incluem ASCII puro sem marcaca
o, formato
de entrada Texinfo, LaTex, SGML ou XML usando um DTD publicamente disponvel, e HTML padr
ao simples,
PostScript ou PDF projetados para modificaca
o por humanos. Exemplos de formatos de imagem transparentes
incluem PNG, XCF e JPG. Formatos Opacos incluem formatos propriet
arios que podem ser lidos e editados
somente por processadores de texto propriet
arios, SGML ou XML para os quais o DTD e/ou ferramentas de
processamento n
ao s
ao largamente disponibilizadas, e HTML, Postscript ou PDF gerados automaticamente
com prop
osito apenas de sada por alguns processadores de texto.
P
agina de Ttulo significa, para um livro impresso, a pr
opria p
agina do ttulo, alem das p
aginas subseq
uentes necess
arias para conter, de forma legvel, o material que esta Licenca requer que apareca na p
agina
do ttulo. Para trabalhos em formatos que n
ao tem uma p
agina de ttulo assim, P
agina de Ttulo significa o
texto pr
oximo a
` ocorrencia mais proeminente do ttulo do trabalho, precedendo o incio do corpo do texto.
Uma seca
o Intitulada XYZ significa uma sub-unidade com nome do Documento cujo ttulo ou e precisamente XYZ ou contem XYZ em parenteses seguindo o texto que traduz XYZ em outra lngua. (Aqui
XYZ representa o nome de uma seca
o especfica mencionado acima, tal como Agradecimentos, Dedicat
oria,
Apoio, ou Hist
orico.) Preservar o Ttulo de uma seca
o assim quando voce modifica o Documento significa
que ela continua sendo uma seca
o Intitulada XYZ de acordo com esta definic
ao.
O Documento pode incluir Notas de Garantia em seguida a
` nota que afirma que esta Licenca se aplica ao
Documento. Estas Notas de Garantia s
ao tidas como inclusas por referencia nesta Licenca, mas somente com
relaca
o a
`s notas de garantia: qualquer outra implicaca
o que estas Notas de Garantia possam ter e anulada e
n
ao tem efeito algum no conte
udo desta Licenca.
2. COPIAS
LITERAIS
Voce pode copiar e distribuir o Documento em qualquer meio, comercialmente ou n
ao-comercialmente,
desde que esta licenca, as notas de direitos autorais (copyright), e a nota de licenca afirmando que esta Licenca
se aplica ao Documento sejam reproduzidas em todas as c
opias, e que voce n
ao inclua outras condico
es,
quaisquer que sejam, a
`s condico
es desta Licenca. Voce n
ao pode usar de medidas tecnicas para obstruir
ou controlar a leitura ou c
opia futura das c
opias que voce fizer ou distribuir. Contudo, voce pode aceitar
compensaca
o em troca das c
opias. Se voce distribuir um n
umero suficientemente grande de c
opias, voce deve
tambem respeitar as condico
es na seca
o 3.
Voce pode tambem emprestar c
opias, sob as mesmas condico
es acima mencionadas, e voce tambem as
pode mostrar publicamente.
3. COPIANDO EM QUANTIDADE
Se voce publicar c
opias impressas (ou c
opias em um meio que normalmente tem capas impressas) do
documento, em n
umero maior que 100, e a nota de licenca do Documento requer Textos de Capa, voce
deve encadernar as c
opias em capas que carreguem, de forma clara e legvel, todos estes Textos de Capa:
Textos de Capa Frontal na capa frontal, e Textos de Quarta Capa na quarta capa. Ambas as capas devem
tambem identificar, de forma clara e legvel, voce como o editor das c
opias. A capa frontal deve apresentar
o ttulo completo com todas as palavras deste igualmente proeminentes e visveis. Voce pode adicionar outro
material nas capas. C
opias com mudancas limitadas a
`s capas, desde que preservando o ttulo do Documento
e satisfazendo estas condico
es, podem ser tratadas como c
opias literais em outros aspectos.
138
Se os textos necess
arios a qualquer uma das capas s
ao demasiado volumosos para serem includos de forma
legvel, voce deve colocar os primeiros listados (quantos couberem razoavelmente) na pr
opria capa, e continuar
o resto nas p
aginas adjacentes.
Se voce publicar ou distribuir c
opias Opacas do Documento em n
umero maior que 100, voce deve ou
incluir uma c
opia Transparente legvel por computador juntamente com cada c
opia Opaca, ou dizer em, ou
juntamente com, cada c
opia Opaca um endereco de rede a partir do qual o p
ublico geral possa acessar e obter,
usando protocolos de rede p
ublicos padr
ao, uma c
opia Transparente completa do Documento, livre de material
adicionado. Se voce decidir pela segunda opca
o, voce deve seguir passos razoavelmente prudentes, quando
comecar a distribuir as c
opias Opacas em quantidade, para garantir que esta c
opia transparente permanecer
a
acessvel no local indicado por pelo menos um ano ap
os a u
ltima vez que voce distribuir uma c
opia Opaca
(diretamente ou atraves de seus agentes ou distribuidor) desta edica
o ao p
ublico.
solicitado, mas n
E
ao exigido, que voce contate os autores do Documento muito antes de redistribuir qualquer n
umero grande de c
opias, para dar a eles uma chance de lhe fornecer uma vers
ao atualizada do Documento.
4. MODIFICAC
OES
Voce pode copiar e distribuir uma Vers
ao Modificada do Documento sob as condico
es das seco
es 2 e 3
acima, desde que voce forneca a Vers
ao Modificada estritamente sob esta Licenca, com a Vers
ao Modificada
no papel de Documento, permitindo assim a distribuica
o e modificaca
o da Vers
ao Modificada a quem quer que
possua uma c
opia desta. Alem disso, voce deve executar os seguintes procedimentos na Vers
ao Modificada:
A. Use na P
agina de Ttulo (e nas capas, se alguma) um ttulo distinto do ttulo do Documento, e dos de
vers
oes anteriores (os quais devem, se houver algum, ser listados na seca
o Hist
orico do Documento). Voce
pode usar o mesmo ttulo que uma vers
ao previa se o editor original daquela vers
ao assim o permitir.
B. Liste na P
agina de Ttulo, como autores, uma ou mais pessoas ou entidades respons
aveis pela autoria
ou modificac
oes na Vers
ao Modificada, juntamente com pelo menos cinco dos autores principais do Documento
(todos seus autores principais, se houver menos que cinco), a menos que estes lhe desobriguem desta exigencia.
C. Mencione na P
agina de Ttulo o nome do editor da Vers
ao Modificada, como seu editor.
D. Preserve todas as notas de direitos autorais (copyright) do Documento.
E. Adicione uma nota apropriada de direitos autorais para suas modificaco
es, adjacente a
`s outras notas
de direitos autorais.
F. Inclua, imediatamente ap
os as notas de direitos autorais, uma nota de licenca dando ao p
ublico permiss
ao
para usar a Vers
ao Modificada sob os termos desta Licenca, na forma mostrada no Adendo abaixo.
G. Preserve naquela nota de licenca a lista completa de Seco
es Invariantes e Textos de Capa requeridos
dados na nota de licenca do Documento.
H. Inclua uma c
opia inalterada desta Licenca.
I. Preserve a seca
o intitulada Hist
orico, preserve seu ttulo, e adicione a esta um item mencionando pelo
menos o ttulo, ano, novos autores, e editor da Vers
ao Modificada conforme includo na P
agina de Ttulo. Se
n
ao houver uma seca
o intitulada Hist
orico no Documento, crie uma mencionando o ttulo, ano, autores e
editor do Documento como mostrado na P
agina de Ttulo, em seguida adicione um item descrevendo a Vers
ao
Modificada como mencionado na sentenca anterior.
J. Preserve o endereco de rede, se algum, dado no Documento para acesso p
ublico a uma c
opia Transparente
deste e, da mesma maneira, os enderecos de rede dados no Documento para vers
oes previas nas quais este se
baseia. Estes podem ser colocados na seca
o Hist
orico. Voce pode omitir um endereco de rede para um
trabalho que foi publicado pelo menos quatro anos antes do Documento em si, ou se o editor original da vers
ao
a
` qual o endereco se refere der permiss
ao.
K. Para qualquer seca
o intitulada Agradecimentos ou Dedicat
oria, preserve o ttulo da seca
o, e preserve
dentro da seca
o toda a subst
ancia e tom de cada um dos agradecimentos e/ou dedicat
orias l
a mencionados.
139
L. Preserve todas as Seco

es Invariantes do Documento, inalteradas no seu texto e ttulos. N
umeros de
seca
o ou o equivalente n
ao s
ao considerados parte dos ttulos das sec
oes.
M. Apague qualquer seca
o intitulada Apoio. Tal seca
o n
ao ser includa na Vers
ao Modificada.
N. N
ao modifique o ttulo de qualquer seca
o a ser intitulada Apoio ou que resulte em conflito com ttulo
de qualquer Sec
ao Invariante.
O. Preserve quaisquer notas de garantia.
Se a Vers
ao Modificada incluir novas seco
es iniciais ou apendices que sejam qualificados como Seco
es
Secund
arias, e n
ao contiver material copiado do Documento, voce pode, a seu criterio, tornar algumas dessas
ou todas essas seco
es em invariantes. Para fazer isso, adicione seus ttulos a
` lista de Seco
es Invariantes na nota
de licenca da Vers
ao Modificada. Estes ttulos devem ser distintos de quaisquer outros ttulos de seco
es.
Voce pode incluir uma seca
o intitulada Apoio, dado que ela contenha nada alem de apoio recebido para
sua Vers
ao Modificada por v
arias fontes por exemplo, notas do revisor ou de que o texto foi aprovado por
uma organizaca
o como a definica
o autoritativa de um padr
ao.
Voce pode adicionar uma passagem de ate cinco palavras como Texto de Capa Frontal, e uma passagem
de ate 25 palavras como Texto de Quarta Capa, ao fim da lista de Textos de Capa na Vers
ao Modificada.
Somente uma passagem de Texto de Capa Frontal e uma de Texto de Quarta Capa pode ser adicionado por
(ou atraves de arranjos feitos por) uma entidade qualquer. Se o Documento j
a incluir um texto de capa para a
mesma capa, previamente includo por voce ou por arranjo feito pela mesma entidade em cujo nome voce est
a
agindo, voce n
ao pode adicionar outro; mas voce pode substituir o antigo, com permiss
ao explcita do editor
anterior, que o incluiu.
O(s) autor(es) e editor(es) do Documento, por esta Licenca, n
ao d
ao permiss
ao para seus nomes serem
usados para publicidade ou defesa ou apoio implcito para qualquer Vers
ao Modificada.
5. COMBINANDO DOCUMENTOS
Voce pode combinar o documento com outros documentos publicados sob esta Licenca, sob os termos
definidos na seca
o 4 acima para vers
oes modificadas, desde que voce inclua na combinaca
o todas as Seco
es
Invariantes de todos os documentos originais, sem modificaco
es, e as liste como Seco
es Invariantes de seu
trabalho combinado, na sua nota de licenca, e que voce preserve todas as Notas de Garantia.
O trabalho combinado somente precisa conter uma c
opia desta Licenca, e m
ultiplas Seco
es Invariantes
identicas podem ser substitudas por uma u
nica c
opia. Se houver m
ultiplas Seco
es Invariantes com o mesmo
nome, porem com conte
udos diferentes, torne o ttulo de cada uma destas seco
es u
nico, adicionando ao fim
dele, entre parenteses, o nome do autor ou editor original desta seca
o, se conhecido, ou ent
ao um n
umero
u
nico. Faca o mesmo ajuste nos ttulos de sec
ao na lista de Seco
es Invariantes na nota de licenca do trabalho
combinado.
Na combinaca
o, voce deve combinar quaisquer sec
oes intituladas Hist
orico nos v
arios documentos
originais, formando uma seca
o intitulada Hist
orico; do mesmo modo, combine quaisquer seco
es intituladas
Agradecimentos, e quaisquer seco
es intituladas Dedicat
oria. Voce deve apagar todas as seco
es intituladas
Apoio.
6. COLEC
OES
DE DOCUMENTOS
Voce pode fazer uma coleca
o consistindo do Documento e outros documentos publicados sob esta Licenca,
e substituir as c
opias individuais desta Licenca, nos v
arios documentos, por uma u
nica c
opia a ser includa na
coleca
o, desde que voce siga as regras desta Licenca para c
opias literais de cada documento em todos os outros
aspectos.
Voce pode extrair um u
nico documento desta coleca
o, e distribu-lo individualmente sob esta Licenca,
desde que voce insira uma c
opia desta Licenca no documento extrado, e siga esta Licenca em todos os outros
140
aspectos com relaca

o a
` c
opia literal do documento.
A TRABALHOS INDEPENDENTES
7. AGREGAC
AO
Uma compilaca
o do Documento ou seus derivados com outros documentos ou trabalhos separados e independentes, dentro de ou junto a um volume de um meio de armazenagem ou distribuica
o, configura um
agregado se os direitos autorais resultantes da compilaca
o n
ao forem usados para limitar os direitos legais
dos usu
arios desta alem do que os trabalhos individuais permitem. Quando o Documento e includo em um
agregado, esta Licenca n
ao se aplica aos outros trabalhos no agregado que n
ao forem, por sua vez, derivados
do Documento.
Se o requerimento do Texto de Capa da seca
o 3 for aplic
avel a estas c
opias do documento, ent
ao, se
o Documento for menor que metade do agregado inteiro, os Textos de Capa do Documento podem ser
colocados em capas que encerrem o Documento dentro do agregado, ou o equivalente eletr
onico das capas
se o Documento estiver em formato eletr
onico. Do contr
ario, eles devem aparecer como capas impressas que
envolvam o agregado inteiro.
8. TRADUC
OES
Uma traduca
o e considerada como sendo um tipo de modificaca
o, ent
ao voce pode distribuir traduco
es do
Documento sob os termos da seca
o 4. A substituica
o de Seco
es Invariantes por traduco
es requer permiss
ao
especial dos detentores dos direitos autorais, embora voce possa incluir traduco
es de algumas ou todas as
Seco
es Invariantes juntamente a
`s vers
oes originais destas. Voce pode incluir uma traduca
o desta Licenca, e
todas as notas de licenca no Documento, e qualquer Nota de Garantia, desde que voce tambem inclua a vers
ao
original em Ingles desta Licenca e as vers
oes originais das notas de licenca e garantia. Em caso de discord
ancia
entre a traduca
o e a vers
ao original desta Licenca ou nota de licenca ou garantia, a vers
ao original prevalecer
a.
Se uma seca
o no Documento for intitulada Agradecimentos, Dedicat
oria, ou Hist
orico, o requerimento
(seca
o 4) de Preservar seu Ttulo (seca
o 1) tipicamente exigir
a a mudanca do ttulo em si.
9. TERMINO
Voce n
ao pode copiar, modifica, sub-licenciar, ou distribuir o Documento a
` exceca
o do modo expressamente provido por esta Licenca. Qualquer outra tentativa de copiar, modificar, sub-licenciar ou distribuir
o Documento e anulada, e implicar
a em termino autom
atico de seus direitos sob esta Licenca. Contudo, as
partes que receberam c
opias, ou direitos, de voce sob esta Licenca n
ao ter
ao suas licencas terminadas enquanto
tais partes permanecerem em total acordo com a Licenca.
10. REVISOES
FUTURAS DESTA LICENCA
A Free Software Foundation pode publicar novas vers

oes revisadas da Licenca de Documentaca
o Livre
GNU de tempos em tempos. Tais vers
oes ser
ao similares em esprito `
a vers
ao presente, embora possam diferir
em detalhes para abordar novos problemas ou quest
oes. Veja http://www.gnu.org/copyleft/.
A cada vers
ao da Licenca e dado um n
umero de vers
ao distinto. Se o Documento especificar que um
n
umero de vers
ao particular desta Licenca ou qualquer vers
ao posterior se aplica a ele, voce tem a opca
o
de seguir os termos e condic
oes ou da vers
ao especificada ou de qualquer vers
ao posterior que tenha sido
publicada (n
ao como rascunho) pela Free Software Foundation. Se o documento n
ao especificar um n
umero de
vers
ao desta Licenca, voce pode escolher qualquer vers
ao j
a publicada (n
ao como rascunho) pela Free Software
Foundation.
ADENDO: Como usar esta Licenca em seus documentos
141
Para usar esta Licenca em um documento que voce escreveu, inclua uma c
opia desta no documento e
adicione as seguintes notas de direitos autorais e licenca logo ap
os a p
agina de ttulo:
garantida a permiss
Copyright (c) ANO SEU NOME. E
ao para copiar, distribuir e/ou modificar
este documento sob os termos da Licenca de Documentaca
o Livre GNU (GNU Free Documentation License), Vers
ao 1.2 ou qualquer vers
ao posterior publicada pela Free Software Foundation;
sem Seco
es Invariantes, Textos de Capa Frontal, e sem Textos de Quarta Capa. Uma c
opia da
licenca e includa na seca
o intitulada GNU Free Documentation License.
Se voce tiver Seco
es Invariantes, Textos de Capa Frontal e Textos de Quarta Capa, substitua a linha sem
Seco
es...Capa. por isto:
com as Seco
es Invariantes sendo LISTE SEUS TITULOS, com os Textos de Capa Frontal sendo
LISTA, e com os Textos de Quarta Capa sendo LISTA.
Se voce tiver Seco
es Invariantes sem Textos de Capa, ou alguma outra combinaca
o dos tres, mescle estas
duas alternativas para se adequar a
` situac
ao.
Se seu documento contiver exemplos n
ao-triviais de c
odigo de programaca
o, recomendamos publicar estes
exemplos paralelamente, sob a licenca de software livre que voce escolher, como por exemplo a Licenca P
ublica
Geral GNU (GNU General Public License), para permitir seu uso em software livre.

Tutorial Grass6 UFAM

Enviado por

Direitos autorais:

Formatos disponíveis

Tutorial Grass6 UFAM

Enviado por

Dados do documento

Descrição original:

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Tutorial Grass6 UFAM

Enviado por

Direitos autorais:

Formatos disponíveis

Geoprocessamento com Software Livre:

Carlos Henrique Grohmann

Sao Paulo, 12 de outubro de 2012

c 2012 Carlos Henrique Grohmann.

2 Mapas Raster - vis

3 Mapas Vetoriais - vis

13 O QGIS como interface gr

14 Produzindo mapas para impress

18 GNU Free Documentation License

Geometia e atributos de dados em SIG . . . . . . . . . . . . . . . . . . . . . . . . . . .

Selecione uma sub-

10.3 Esquema da interpolac

Estrutura dos comandos do GRASS . . . . . . . . . . . . .

Algumas fontes de informac

Leia mais sobre Software Livre em http://www.gnu.org

Design e estrutura de SIGs

QGIS web site oficial: http://qgis.org

Dassau et al. (2005).

An introduction to the practical use of the Free Geographi-

Design e estrutura de Sistemas de Informac

De maneira geral, um Sistema de Informacoes Geograficas pode ser caracterizado por

Design e estrutura de SIGs

Figura 1.1. Geometia e atributos de dados em SIG

Figura 1.2. Representac

Dados espaciais podem ser representados por pontos, linhas,

Estrutura dos comandos no GRASS

de SIG nos permitem visualizar as superfcies em perspectivas tri-dimensionais, porem apesar

Figura 1.3. Dimens

Estrutura dos comandos no GRASS

O GRASS-GIS e um Sistema de Informacoes Geograficas com funcoes orientadas a dados

Obtendo ajuda com relac

Arquivos de ajuda est

Significado dos comandos

Exibicao e consulta grafica no Map Display

Tabela 1. Estrutura dos comandos do GRASS

detalhadas encontram-se nas p

modificados acidentalmente. Outros usuario (e tambem a pessoa com permissao de escrita

Arquivos do mapset PERMANENT

2 - Mapas Raster - vis

Mapas Raster - vis

recebe o valor de um atributo, que representa um fenomeno (por exemplo temperatura ou

Figura 2.1. Resoluca

3 - Mapas Vetoriais - vis

Tabela 2. Comandos do GRASS para importaca

O modulo mais utilizado para importacao (e exportacao) de rasters e certamente o

Mapas Vetoriais - vis

3 - Mapas Vetoriais - vis

SHAPE file, UK.NTF, SDTS, TIGER, S57, MapInfo-File,

Tabela 3. Comandos do GRASS para importaca

A escala de um mapa define a relacao entre a medida de um objeto ou lugar representado

Figura 4.1. Exemplos de escalas gr

A escala e uma divis

O geoide e considerado como uma representacao matematica bastante precisa da forma da

Figura 4.2. Representac

Figura 4.3. Relaco

eixo menor (b). As relac

Figura 4.4. Elementos do elips

Semi-eixo maior a (m)

Semi-eixo menor b (m)