1

Tutorial de QGIS 2.0.1

Profa. Adj. Arlete Meneguette Universidade Estadual Paulista Campus de Presidente Prudente

arletemeneguette@fct.unesp.br https://unesp.academia.edu/ArleteMeneguette/Papers

1. Fazer download e instalar o QGIS

Para obter gratuitamente o QGIS entre em http://qgis.org/pt_PT/site/.

O aplicativo est disponvel para Windows, MacOS X, Linux, BSD e Android. Depois de baixar o executvel, instale o aplicativo em seu computador. A verso utilizada neste Tutorial a 2.0.1, a qual foi lanada em setembro de 2013. O sistema operacional utilizado o Windows 7 em um notebook 64 bits. Na dvida sobre a configurao do seu computador, baixe a verso 32 bits.

Saiba como instalar o QGIS 2.0.1 no Windows 7 acessando o excelente tutorial disponvel em http://www.processamentodigital.com.br/2013/09/18/quantum-gis-2-0-1-tutorial-para-instalacao-no-windows-7/. Aprenda tambm a incluir dois importantes complementos: Table Manager e Open Layers.

O Manual do Usurio pode ser acessado em http://qgis.org/pt_PT/docs/index.html

2

2. Criar e Organizar Pastas no seu computador

Ao longo deste Tutorial de QGIS vamos precisar fazer download e salvar arquivos, bem como criar novos arquivos. Crie em seu computador uma pasta denominada TutorialQGIS para organizar seus dados. Nessa pasta crie as seguintes subpastas: Arquivos-ZIP Arquivos-SHP Arquivos-TIF Arquivos-QSG Arquivos-QML Novos-TIF Novos-SHP

3. Fazer download de dados vetoriais do site do IBGE

Entre em http://downloads.ibge.gov.br/downloads_geociencias.htm e clique na opo Geocincias, na qual h vrias pastas, clique em malhas_digitais. Na pasta malhas_digitais h subpastas, uma delas denominada municipio_2010 e outra subpasta chamada municpio_2013 (mas nela no h dados para o Estado de So Paulo, somente para os estados onde houve emancipao e criao de novos municpios, quais sejam: Par, Santa Catarina, Rio Grande do Sul e Mato Grosso do Sul). Portanto, clique em municipio_2010, na qual h muitos arquivos ZIP, um para cada Estado. Clique com o boto direito do mouse no arquivo SP.zip. Escolha a opo Salvar link como.

Na pasta TutorialQGIS do seu computador voc criou uma subpasta denominada Arquivos-ZIP. Salve nela o arquivo SP.zip, em seguida, descompacte-o na subpasta Arquivos-SHP, o que far com que voc tenha uma pasta SP contendo vrios arquivos e aqueles que mostram a Unidade da Federao (ou seja, o Estado) apresentam a sigla UF, enquanto que os que mostram a Mesorregio tm a sigla ME, a Microrregio aparece como MI e os Municpios como MU. importante ressaltar que o Estado de So Paulo tem Geocdigo 35 e os dados esto em Coordenadas Geogrficas (GC) em Sirgas2000 (SIR), originalmente em 1:250.000 (sem supresso de pontos).

3

4. Fazer download de dados matriciais do site da Embrapa

Neste Tutorial tambm vamos fazer uso de dados SRTM disponibilizados no site da Embrapa, portanto, entre em http://www.relevobr.cnpm.embrapa.br/. Voc ter acesso a imagens SRTM (Shuttle Radar Topography Mission). As imagens do relevo realado foram geradas pela Embrapa em composies coloridas RGB (Red, Green e Blue), conforme adaptao do padro de tabela de cores da USGS (United States Geological Survey). Assim, as cotas mais baixas do terreno apresentam-se em verde escuro e as mais elevadas em tons de rosa. Os valores intermedirios distribuem-se em tons de verde, amarelo, marrom e rosa, seguindo o aumento da elevao. No mapa do Brasil clique sobre o Estado de So Paulo. Para apoiar o entendimento do relevo realado, a EMBRAPA gerou uma escala altimtrica com vinte intervalos. Como cada Estado apresenta valores especficos das altitudes mnima e mxima, cada escala de cores tambm representa valores e variaes especficas. Cabe ressaltar ainda que esta escala fixa para todos os nveis de ampliao quando navegamos em um mesmo Estado, indicando sempre o valor mnimo e mximo da totalidade do Estado.

Caso queira saber mais sobre o assunto entre posteriormente em http://www.relevobr.cnpm.embrapa.br/conteudo/infotec.htm. Para fazer download das imagens SRTM clique em Dados para Download. Escolha o Estado de So Paulo.

Escolha um dos arquivos para fazer download, por exemplo, aquele que corresponde ao mesmo recorte da carta topogrfica SF-22-Y-A em 1:250.000. Observe que os dados SRTM esto em Coordenadas Geogrficas e no Datum WGS84. Clique com o boto direito do mouse sobre e escolha Salvar link como.

4

Salve o arquivo sf-22-y-a.zip na subpasta Arquivos-ZIP no seu computador. Retorne tela anterior que mostra o Estado de So Paulo. Repita o procedimento para fazer download de todos os arquivos ZIP, salve-os na subpasta Arquivos-ZIP. Depois descompacte os arquivos ZIP na subpasta Arquivos-TIF da pasta TutorialQGIS no seu computador.

5. Criar e configurar o Projeto

Execute o QGIS e verifique que ser exibida a interface do aplicativo, composta por barra de menus (no topo da tela), barra de ferramentas (abaixo da barra de menus), barra de status (na parte de baixo da tela), painel de camadas ( esquerda) e rea de visualizao de mapas ( direita).

A interface do QGIS pode estar sendo exibida em outra lngua, por exemplo, em Ingls. Vamos agora demonstrar como alterar a Lngua. Clique em Settings->Options->Locale. Habilite a caixa para sobrescrever o local do sistema. Em seguida escolha Portugus do Brasil a partir da lista.

5

Clique em OK. Feche o QGIS. Depois execute o QGIS novamente para que a alterao tenha efeito. Note que agora a interface mostrada em Portugus do Brasil. No canto inferior direito

da barra de status h um cdigo alfanumrico e um cone: .

Clique nesse cone e mostrada uma caixa de dilogo (denominada Propriedades do Projeto), na qual podemos constatar que por padro o QGIS assume Latitude e Longitude (como sistema de coordenadas) e o WGS84 (como sistema de referncia).

O European Petroleum Survey Group (EPSG) definiu os cdigos EPSG com o objetivo de sistematizar os Sistemas de Referncia de Coordenadas (SRC). Os cdigos EPSG mais utilizados no Brasil so mostrados a seguir. A lista completa e detalhada pode ser encontrada em http://www.processamentodigital.com.br/2013/07/27/lista-dos-codigos-epsg-mais-utilizados-no-brasil/

CDIGO EPSG COORDENADAS / DATUM 4225 GCS Crrego Alegre 4618 GCS SAD69 4674 GCS SIRGAS 2000 4326 GCS WGS84

6

Cabe ressaltar que GCS significa Geographic Coordinate System, ou Sistema de Coordenadas Geogrficas, assim como UTM se refere ao Sistema de Projeo Universal Transversa de Mercator. Por exemplo, para a Zona 22 Sul os cdigos EPSG correspondentes so:

CDIGO EPSG COORDENADAS / DATUM 22522 UTM zona 22 S / Crrego Alegre 29192 UTM zona 22 S / SAD69 31982 UTM zona 22 S / SIRGAS 2000 32722 UTM zona 22 S / WGS 84

Os dados que utilizaremos neste Tutorial esto definidos em diferentes sistemas de coordenadas e em diversos sistemas de referncia. Portanto, devemos fazer uso dos cdigos EPSG adequados. Tendo em vista que no Brasil devemos adotar o SIRGAS 2000 como sistema de referncia, vamos fazer essa alterao no QGIS. Habilite a caixa que permite fazer a transformao SRC on the fly, ou seja, em tempo real. Em Filtro digite 4674 que o cdigo EPSG para SIRGAS 2000. Note que agora essa opo listada, em Sistema de Referncia do World clique sobre ela e observe a sintaxe: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs. Isso significa que as coordenadas so geogrficas (latitude e longitude), o elipsoide o GRS80, alm disso, no aplicada nenhuma translao entre os sistemas, pois segundo o IBGE: No existem parmetros de transformao entre SIRGAS2000 e WGS84 porque eles so praticamente iguais, ou seja, DX = 0, DY = 0 e DZ = 0 , conforme pode ser visto em http://www.ibge.gov.br/home/geociencias/geodesia/pmrg/faq.shtm#15.

Clique em Aplicar e observe que agora o cdigo EPSG 4674 exibido no canto inferior direito na barra de status: . Clique em OK.

7

6. Salvar, Fechar e Abrir o Projeto

A partir da barra de menus, clique em Projeto->Salvar Como. exibida uma caixa de dilogo solicitando que voc indique o caminho onde ser salvo o seu Projeto. Neste Tutorial foi indicado o drive C, a pasta TutorialQGIS e a subpasta Arquivos-QSG. Digite EstadoSP como nome do Projeto, pois vamos adotar como rea de estudo o Estado de So Paulo.

Clique em Salvar. Depois clique em Projeto->Sair do QGIS. Execute o QGIS novamente e verifique que nenhum Projeto est ativo. A partir da barra de menus, clique em Projeto->Abrir Recente, escolha C:/TutorialQGIS/Arquivos-QSG/EstadoSP.qgs.

Observe que novamente, na barra de ttulo do QGIS exibido o nome do Projeto ativo.

7. Criar Grupos no Painel de Camadas

Clique com o boto direito do mouse no Painel de Camadas. Escolha a opo Adicionar novo grupo, pois vamos nos preparar para organizar os dados que sero importados a seguir.

O novo grupo criado automaticamente com o nome de grupo 1. Clique com o boto direito do mouse sobre o nome do grupo1 e escolha a opo Renomear.

8

Renomeie o grupo como CG_SIR. A escolha desse nome para o grupo se deve ao fato que vamos associar a ele os dados que baixamos do site do IBGE e que esto em Coordenadas Geogrficas em SIRGAS2000. Clique novamente com o boto direito do mouse e escolha a opo Adicionar Novo Grupo. Renomeie o novo grupo como CG_WGS. Nesse grupo estaro associados dados que esto em Coordenadas Geogrficas em WGS84, os quais foram baixados do site da Embrapa.

Outros grupos podero ser criados posteriormente, seguindo o mesmo procedimento demonstrado neste item do Tutorial de QGIS.

8. Importar arquivo SHP e alterar as propriedades

Na barra de ttulo do QGIS exibido o nome do Projeto ativo: . No Painel de Camadas constam somente os dois grupos que criamos anteriormente, mas no listado nenhum item ainda. A partir da barra de menus clique em Camada->Adicionar camada vetorial ou clique no cone . Uma caixa de dilogo exibida na qual possvel escolher o Tipo de Fonte (mantenha a opo Arquivo) e indicar o caminho onde o arquivo foi salvo. Altere a codificao de System para CP1252, pois assim no teremos problemas no reconhecimento de caracteres e acentuao na tabela de atributos. Agora clique em Buscar. Vamos importar dados vetoriais no formato SHP que baixamos do site do IBGE. Para tanto, na caixa de dilogo (denominada Abrir uma camada vetorial OGR suportada) indique a subpasta Arquivos-SHP que voc criou anteriormente na pasta TutorialQGIS no seu computador. Observe que o formato SHP j est selecionado na lista suspensa localizada no canto inferior direito da caixa de dilogo. Alm de SHP, muitas extenses so aceitas pelo QGIS. H vrios arquivos de extenso shp e aquele que mostra a Unidade da Federao (ou seja, o Estado) a que contm a sigla UF. Tendo selecionado o arquivo 35UFE250GC_SIR.shp clique em Abrir. De volta caixa de dilogo anterior note que o caminho escolhido e o arquivo selecionado agora aparecem como Fonte-Conjunto de dados.

Clique em Abrir. O QGIS exibe agora o polgono do Estado de So Paulo na rea de visualizao de mapas e no painel de camadas consta 35UFE250GC_SIR. Clique com o boto direito do mouse sobre o nome da camada e escolha a opo Renomear. Substitua o nome da camada por Estado SP Limites. Arraste a camada para dentro do grupo CG-SIR. Clique novamente com o boto direito do mouse sobre o nome da camada e desta vez escolha Propriedades (ou ento clique duas vezes sobre o nome da camada para acessar a caixa de dilogo de Propriedades). Note que na caixa de dilogo na aba Geral so mostrados alguns dados, dentre os quais o nome da camada e o caminho onde o arquivo original foi salvo anteriormente. Isso se deve ao fato que o QGIS no salva o arquivo, mas sim o caminho onde o arquivo est disponvel para ser exibido e editado. O Sistema de Referncia de Coordenadas (SRC) consta como sendo um datum desconhecido baseado no elipside GRS 1980. Mas ns sabemos que o sistema de referncia o SIRGAS 2000, pois essa informao foi disponibilizada pelo IBGE. Note que ao deslocar o mouse sobre a rea de visualizao de mapas as coordenadas geogrficas so exibidas na barra de status.

9

Sendo assim, em Sistema de Referncia de Coordenadas clique em Especifique. Ser aberta outra caixa de dilogo (Seletor de Sistema de Coordenadas de Referncia). Em Filtro digite 4674, depois em Sistema de Referncia do World, em Sistema de Coordenadas Geogrficas, clique em SIRGAS 2000 (ou seja, EPSG:4674), que o mesmo sistema de referncia do Projeto. Observe a sintaxe: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs.

Clique em OK. De volta caixa de Propriedades da Camada note as alteraes efetuadas.

10

Clique agora na aba Estilo. Observe que a cor mostrada neste tutorial pode estar diferente da cor no seu computador, pois o QGIS aloca cores aleatoriamente at que o usurio escolha aquela da sua preferncia. Portanto, altere a cor de preenchimento do polgono. Clique em OK na caixa Selecionar Cor, depois clique em Aplicar na caixa de Propriedades da Camada. Veja que o Estado de So Paulo agora mostrado como um polgono preenchido na cor escolhida. Entretanto, ao invs de mostrarmos o Estado de So Paulo como um polgono, podemos mostrar seus limites como linhas. Para tanto, clique sobre Preenchimento Simples e note que agora outros recursos so exibidos, por exemplo, em tipo de camada escolha Contorno Linha Simples. Note que a cor j est selecionada como preto, mas voc pode alterar se desejar, assim como pode mudar a espessura e estilo da caneta, dentre outros.

Clique em Aplicar e veja que agora o Estado de So Paulo mostrado com seus limites em preto e no mais como um polgono preenchido com uma cor slida. Clique em OK para fechar a caixa de propriedades da camada. Clique com o boto direito do mouse sobre o nome da camada Estado SP Limites e escolha a opo Abrir tabela de atributos. Observe que na tabela temos somente uma feio. Note que o nome do Estado (NM_ESTADO) So Paulo e o nome da regio (NM_REGIAO) Sudeste.

11

Clique sobre o registro cujo ID 6, CD-GEOCODU 35 e note que agora o limite do Estado de So Paulo mostrado como linha amarela (cor padro do QGIS para feies selecionadas).

Clique em para Desselecionar tudo e note que o limite volta a ser mostrado na cor definida anteriormente (preto). Feche a janela de atributos.

9. Alterar as configuraes

Se desejar alterar a cor da seleo para que no seja mostrada em amarelo, a partir da barra de menus clique em Configuraes->Opes e na aba Tela & Legenda voc poder fazer as alteraes necessrias.

Clicando na aba Geral voc ver outros recursos, por exemplo, possvel alterar o tamanho do cone e o tamanho da fonte exibida. Neste Tutorial o tamanho do cone foi alterado para 32 e o tamanho da fonte para 11, para permitir uma melhor leitura. Agora clique na aba Renderizao; em Rasters Melhor contraste mantenha a opo Estender para MinMax em todas as alternativas associadas com as Bandas. Note que em Limites (mnimo/mximo) est selecionada a opo Corte de Contagem de Pixel Cumulativo e que os Limites de Corte de Contagem Cumulativa de Pixel esto configurados para 2% e 98%. Altere a opo para Mnimo/Mximo.

12

A partir de agora, todo arquivo raster que for carregado no QGIS exibir seus valores mnimo e mximo. Feche a caixa de opes, clicando em OK.

Clique em Projeto->Salvar, sempre que for sair do QGIS e depois ao retornar, abra o projeto recente.

10. Importar arquivos TIF e exibir as propriedades

Agora vamos importar dados matriciais (raster) que baixamos do site da Embrapa, os quais esto na extenso TIF. Trata-se de imagens SRTM que esto salvas na sub-pasta Arquivos-

TIF da pasta TutorialQGIS. Clique no cone para adicionar uma camada raster. Uma caixa de dilogo (denominada Abrir uma fonte de dados raster GDAL suportada) exibida na qual possvel indicar o caminho onde o arquivo foi salvo, portanto, clique duas vezes em Arquivos-TIF. Observe que o formato GDAL j est selecionado na lista suspensa localizada no canto inferior direito da caixa de dilogo. Podemos verificar que muitas extenses so aceitas pelo QGIS. Selecione somente os arquivos SF-22-Y-A.TIF e SF-22-Y-B.TIF e depois clique em Abrir. Note que agora no Painel de Controle mais duas camadas so listadas, correspondentes s folhas topogrficas SF-22-Y-A e SF-22-Y-B, localizadas no extremo oeste do Estado de So Paulo, na regio conhecida como Pontal do Paranapanema. Arraste as duas camadas raster para dentro do grupo CG-WGS e expanda as duas pastas para exibir os valores mnimo e mximo de cada uma delas. A camada que representa o Estado de So Paulo serve como referncia visual para vermos onde esto localizadas as duas imagens.

13

Exiba somente a camada SF-22-Y-A, depois clique duas vezes sobre o nome da camada para exibir as Propriedades. Constate que o Sistema de Referncia de Coordenadas j est atribudo como sendo EPSG:4326 (WGS84), conforme consta no site da Embrapa.

Embora a extenso dos dois arquivos matriciais seja TIF, na realidade, trata-se de arquivos GeoTIF, pois esto georreferenciados, sendo que as coordenadas so geogrficas e o datum WGS84. Isso pode ser comprovado clicando na opo Metadata da caixa de dilogo Propriedades da camada, onde outras relevantes informaes podem ser acessadas. Por exemplo, possvel verificar que o valor mximo no histograma 629 e o mnimo 185, que a imagem tem 1800 pixels em X (colunas) e 1200 pixel em Y (linhas), que tem somente uma banda, que a origem est localizada na longitude -54.0004 e na latitude -22.0004 (ou seja, o canto inferior esquerdo da imagem). Consta tambm que o pixel tem o mesmo tamanho em X e em Y, qual seja 0.000833333 graus decimais (que correspondem a cerca de 90 metros). .

14

Clique agora na aba Estilo e observe que em Renderizao da banda o tipo de renderizao consta como sendo Banda simples cinza, sendo que realmente o valor mnimo 185 e o valor mximo 629.

Entretanto, como voc pode notar em Carregar Valores de min/max est habilitada a opo Corte de Contagem Cumulativa, na qual por default consta um intervalo de 2% a 98%. H outras opes, tambm, uma delas Mn/Mx. A outra Mdia+/-Desvio Padro que pode ser multiplicada por um fator selecionado a partir de uma lista suspensa, sendo que por default o valor 1. Clique na aba Histograma e constate que originalmente no h um grfico traado ainda.

15

Se voc clicar em Calcular Histograma ser exibido o grfico com o Valor do Pixel em X e a Frequncia em Y. Certifique-se que esteja habilitada a opo Sempre exibir marcadores de mn/mx, que pode ser selecionada em Prefs/Aes. Observe que o valor mnimo consta como 185 e o valor mximo como 629, os quais so os limiares do grfico sendo representados por linhas tracejadas nas duas extremidades dos valores de pixel.

Feche a caixa de Propriedades da Camada SF-22-Y-A. Desselecione a camada SF-22-Y-A e agora exiba somente a camada SF-22-Y-B. Clique duas vezes sobre o nome dessa camada para exibir suas Propriedades. Observe na aba Estilo quais so os valores mnimo (223) e mximo (637) dos pixels da banda simples cinza. Note que so diferentes dos valores da camada SF-22-Y-A. Feche a caixa de Propriedades da Camada SF-22-Y-B.

11. Criar um mosaico de duas imagens SRTM

At o momento importamos duas imagens SRTM que podem ser exibidas uma de cada vez ou as duas ao mesmo tempo, lado a lado, entretanto, so duas imagens separadas. Vamos agora demonstrar como produzir um mosaico das duas imagens.

16

Exiba as camadas SF-22-Y-A e SF-22-Y-B. Aplique zoom na rea para visualizar as imagens com melhor resoluo na tela e observe a diferena nos tons de cinza.

A partir da barra de menus clique em Raster, escolha Miscelnea->Mosaico. Na caixa Mesclar possvel escolher os Arquivos de Entrada, portanto, clique em Selecione. Em seguida selecione somente os arquivos SF-22-Y-A e SF-22-Y-B na sub-pasta Arquivos-TIF da pasta TutorialQGIS (use a tecla CTRL para facilitar a seleo de ambos os arquivos). Note que os nomes dos arquivos selecionados so listados em Nome do Ficheiro. Clique em Abrir. Agora clique em Selecionar na frente de onde estiver escrito Arquivos de Sada na caixa de dilogo Mesclar. Indique a subpasta Novos-TIF como o local onde ser salvo o mosaico e digite o nome do arquivo como SF-22-Y-A_SF-22-Y-B. Clique em Gravar. Note que agora na caixa de dilogo Mesclar so listados tanto os nomes dos arquivos de entrada como o nome do arquivo de sada e seus respectivos caminhos. Alm disso, j est selecionada a opo Carregar na tela ao concluir. Habilite o campo Opes de Criao, em Perfil mantenha a opo Padro, depois clique em . Em seguida digite TFW em Nome e digite YES em Valor.

17

Cabe ressaltar que TFW significa World File for TIF Image e nesse arquivo so armazenados dados extremamente importantes, que devem ser salvos na mesma subpasta onde o arquivo TIF estiver. No altere os demais itens, simplesmente clique em OK. Clique em OK nos dois avisos e depois clique em Fechar na caixa Mesclar. Note que no Painel de Camadas agora listada a nova camada que foi criada na mosaicagem. Desselecione as duas imagens independentes. Arraste a camada SF-22-Y-A_SF-22-Y-B para o grupo CG_WGS, exiba a camada e expanda a pasta para saber quais so os valores mnimo (208) e mximo (627). Note que a imagem mosaicada no apresenta diferena da tonalidade de cinza que notamos anteriormente, ao exibir as imagens independentes. Para se certificar, exiba as Propriedades da Camada e observe na aba Estilo que o valor mnimo realmente 208 e o valor mximo 627.

12. Criar um mosaico a partir de imagens SRTM salvas em uma pasta

No item anterior produzimos um mosaico usando duas imagens SRTM que cobrem a rea do Pontal do Paranapanema. Agora vamos criar um novo mosaico de uma rea maior, usando os arquivos que salvamos anteriormente na subpasta Arquivos-TIF da pasta TutorialQGIS. A partir do menu principal do QGIS clique em Raster, escolha Miscelnea->Mosaico. Na caixa Mesclar possvel escolher a Pasta de Entrada ao invs de arquivos, portanto, habilite essa opo:

. Em seguida clique em Selecione. Clique uma s vez no nome da subpasta Arquivos-TIF da pasta TutorialQGIS. Note que o nome da Pasta selecionada listado em Nome da Pasta. Clique em Choose. Agora clique em Selecionar na frente de onde estiver escrito Arquivos de Sada na caixa de dilogo Mesclar. Indique a subpasta Novos-TIF como o local onde ser salvo o mosaico e digite o nome do arquivo como Mosaico_EstadoSP. Clique em Gravar. Note que agora na caixa de dilogo Mesclar so listados tanto o nome da pasta contendo os arquivos de entrada como o nome do arquivo de sada e seus respectivos caminhos. Alm disso, j est selecionada a opo Carregar na tela ao concluir. Habilite o campo Opes de Criao, em Perfil mantenha a opo Padro, depois clique em . Em seguida digite TFW em Nome e digite YES em Valor.

18

No altere os demais itens, simplesmente clique em OK. Clique em OK nos dois avisos que sero mostrados ao final do processamento. Clique em Fechar na caixa Mesclar e note que no Painel de Camadas agora listada a nova camada que foi criada na mosaicagem. Arraste a camada Mosaico_EstadoSP para o grupo CG_WGS. Clique com o boto direito do mouse sobre o nome da camada e escolha a opo Aproximar a Camada sua Extenso. Note que agora todo o Estado de So Paulo exibido.

13. Trocar projeo do mosaico de imagens SRTM

Como podemos ver na aba Geral das Propriedades da Camada Mosaico_EstadoSP temos como SRC o EPSG: 4326 WGS84.

19

Para darmos continuidade ao nosso tutorial vamos trocar o sistema de referncia WGS84 para SIRGAS 2000. Sendo assim, a partir da barra de menus clique em Raster->Projees->Trocar Projeo. Na caixa trocar projeo, em Arquivo de Entrada selecione a camada Mosaico_EstadoSP. Na frente de Arquivo de Sada clique em Selecione e na caixa de dilogo que se abre escolha a pasta TutorialQGIS e a subpasta Novos-TIF, entre nessa nova subpasta, digite o nome da nova camada como Mosaico_Reprojetado e clique em Gravar. Em

clique em Selecione. Ser aberta nova caixa de dilogo onde em Filtro voc dever digitar 4326 e depois selecionar WGS84. Clique em OK na caixa onde foi selecionado o SRC. De volta caixa Trocar Projeo, habilite a opo e clique em Selecione. A caixa de seleo de SRC ser exibida novamente, na qual em Filtro voc dever digitar 4674 e depois clicar em SIRGAS2000. Clique em OK na caixa onde foi selecionado o SRC. De volta caixa anterior note que j est habilitada a opo . Habilite Mtodo de Reamostragem e escolha Bilinear.

Clique em OK na caixa Trocar Projeo, depois clique em OK nos dois avisos que sero mostrados, feche a caixa. No painel de camadas arraste Mosaico_Reprojetado para dentro do grupo CG_SIR, abaixo da camada Estado SP Limites. Expanda a pasta da nova camada contendo do mosaico reprojetado e verifique que o valor mnimo 0 e o mximo 2521.

20

Exiba as Propriedades da Camada para certificar-se que foi mesmo feita a troca de projeo e que o SRC consta como EPSG: 4674.

Verifique tambm os valores mnimo e mximo na aba Estilo.

Clique em Projeto->Salvar.

21

14. Recortar o mosaico usando mscara

Para recortar o mosaico de imagens SRTM usando como mscara a representao vetorial do Estado de So Paulo, a partir da barra de menus clique em Raster->Extrao->Cortador. Na caixa de dilogo Cortador clique em Selecionar para indicar o arquivo de entrada que ser recortado. Selecione Mosaico_Reprojetado na lista suspensa. Clique em Abrir. De volta caixa Cortador clique em Selecionar para indicar o caminho e o nome do arquivo de sada. Escolha a subpasta Novos-TIF e digite o nome do novo arquivo como Mosaico_Recortado. Clique em Gravar. De volta caixa Cortador observe que no Modo Clipping h duas opes, clique em Camada Mscara. Indique qual camada ser utilizada como mscara para recortar o mosaico. H vrias opes disponveis no nosso projeto, escolha Estado_SP.

Clique em OK. Quando o processo tiver sido concludo ser mostrado um aviso, clique em OK. Note que no Painel de Camadas agora listada a nova camada que foi recortada. Arraste a camada Mosaico_Recortado para o grupo CG_SIR, abaixo da camada Estado SP Limites. Expanda a pasta para mostrar os valores mnimo e mximo.

22

Observe que no entorno do Estado de So Paulo todos os pixels so mostrados em preto. Use

a ferramenta para Identificar Resultados, clique no pixel em preto e verifique que o valor 0 (zero).

Feche a caixa de Identificar Resultados. Exiba as Propriedades da Camada e note que na aba Transparncia tem um campo no preenchido associado a Sem Valor de Dados Adicionais. Digite 0 (zero corresponde a preto) no campo que est sem preenchimento. Clique em Aplicar. Observe que agora os pixels no entorno do Estado de So Paulo so mostrados em branco.

Aproveite para verificar na aba Geral que o SRC consta realmente como EPSG: 4674.

23

15. Produzir mapa hipsomtrico

Nas Propriedades da Camada Mosaico_Recortado veja que na aba Estilo o valor mnimo 0 e o valor mximo 2476, indicando as variaes dos tons de cinza. Note que na aba Estilo est selecionada a opo Banda simples cinza. Altere para Banda simples falsa-cor a partir da lista suspensa. Veja que agora consta que o valor mnimo 29.325 e o valor mximo 1051.6, ou seja, valores que variam de 0 a 1100, mas para mantermos uma equidistncia de 200 metros podemos adotar 1200 como valor mximo.Em Gerar novo mapa de cores escolha a paleta YlOrBr (Yellow Orange Brown) e em Modo escolha Intervalo Igual. Em Classes escolha 5. Altere o valor mnimo do mapa de cores para 200 e o mximo para 1200 e clique em Classificar. Observe que agora mostrada uma prvia de legenda com os valores numricos associados com cada uma das 6 classes.

Clique em Rtulo e edite os textos para que indiquem os limiares das classes: 0 a 200m, 200 a 400m, 400 a 600m, 600 a 800m, 800 a 1000m, 1000 a 1200m.

24

Clique em Aplicar e depois em OK. O resultado obtido uma Representao Hipsomtrica do Estado de So Paulo, com 6 classes de altitude, com equidistncia de 200 metros, variando de 0 a 1200 metros de altitude.

16. Duplicar camada

A camada que contm a Representao Hipsomtrica do Estado de So Paulo se chama Mosaico_Recortado e no induz o usurio a saber qual seu contedo. Clique com o boto direito do mouse sobre o nome da camada e escolha a opo Duplicar. Desselecione a camada Mosaico_Reprojetado. Observe que no painel de camadas exibido o nome da nova camada como sendo Mosaico_Reprojetado copiar, portanto apague esse nome e digite Mapa_Hipsometrico.

Tal produto cartogrfico associa faixas de altitude das res emersas com cores hipsomtricas, enquanto que para as reas submersas so utilizadas cores batimtricas (tons de azul, sendo mais claro nas reas mais rasas e azul mais escuro nas reas mais profundas).

17. Salvar e aplicar Estilo

Na nova camada os pixels que eram pretos nas bordas no mosaico original agora so mostrados em amarelo! Exiba as Propriedades da Camada e na aba Estilo observe que foi mantido o mesmo estilo que tnhamos definido anteriormente.

25

Clique na aba Transparncia, digite 0 na frente de Sem Valor de Dados Adicionais e em seguida clique em Aplicar. Veja que os pixels das bordas que eram mostrados em amarelo agora no aparecem mais.

Para que possamos conservar as caractersticas definidas para este Mapa Hipsomtrico e reutilizar a mesma aparncia em outro mapa, clique em Salvar estilo. Indique o local onde deseja salvar o arquivo QML (escolha a subpasta Arquivos-QML na pasta TutorialQGIS). Em seguida digite o nome do estilo: Hipsometrica. Clique em Salvar. No futuro se voc precisar criar outro mapa com as mesmas caractersticas deste, basta clicar em Carregar Estilo e indicar o local onde salvou o arquivo Hipsometrica.qml. Clique em OK na caixa de propriedades da camada. Com esses procedimentos demonstramos que o usurio pode customizar a aparncia dos seus mapas como desejar, levando em conta os dados de entrada da sua rea de estudo.

26

18. Recompor a aparncia da camada

Exiba agora as Propriedades da camada Mosaico_Reprojetado e constate que, apesar de termos alterado a aparncia visual da camada, na aba Metadados o valor mnimo continua sendo de 0 e o valor mximo continua como 2476.

Vamos agora demonstrar como gerar uma representao com tais valores. Clique na aba Estilo e altere o Tipo de Renderizao para banda simples cinza. Digite o valor mnimo como 0 e o valor mximo como 2476, em Melhorar Contraste escolha Estender para MinMax.Clique em Aplicar e veja o resultado! O mosaico recortado agora tem novamente a aparncia de antes.

19. Produzir mapa de relevo sombreado

Vamos gerar um mapa de relevo sombreado, fazendo uso da camada Mosaico_Recortado. Clique em Raster->Anlise do Terreno->Sombreamento. Na caixa de dilogo podemos escolher tanto a camada com as elevaes quanto a camada de sada, sendo assim, escolha como camada de entrada Mosaico_Recortado, a partir da lista suspensa. Para selecionar a camada de sada escolha a subpasta Novos-TIF e digite o nome como Mosaico_Sombreado. Clique em Gravar. De volta caixa Sombreamento, altere o azimute para 315 e o ngulo vertical como 45.

27

Clique em OK. Arraste a camada Mosaico_Sombreado para dentro do grupo CG-SIR. Exiba o resultado e expanda a pasta para mostrar os limiares mnimo e mximo. Exiba as Propriedades da Camada e constate que o SRC continua sendo EPSG 4674.

Feche as Propriedades da Camada Mosaico_Sombreado. Repita o procedimento de gerao de imagem sombreada, mas desta vez salve o resultado como Relevo_Sombreado2 e digite 135 para o Azimute e 45 para o ngulo. Clique em OK.

28

Arraste a nova camada para dentro do grupo CG_SIR. Exiba o resultado obtido. Se desejar, assegure-se que o SRC tambm pe EPSG: 4674.

Aplique zoom na rea do Pontal do Paranapanema para exibir mais detalhes da imagem. Inicialmente amplie a camada Mosaico_Sombreado e depois faa o mesmo com a camada Mosaico_Sombreado2 para comparar ainda melhor as duas representaes do relevo sombreado. Verifique que a segunda d a impresso de relevo invertido, ou seja, os fundos de vale parecem ser cristas e vice versa. Isso se deve ao valor numrico definido para o Azimute (o qual varia de 0 a 360 no sentido horrio contado a partir do Norte).

29

20. Remover camada

Clique com o boto direito do mouse sobre o nome da camada Mosaico_Sombreado2 e escolha a opo Remover.

21. Produzir mapa de relevo colorido

Vamos demonstrar agora como gerar mapas de relevo colorido. Na barra de menus clique em Raster->Anlise do Terreno->Relevo. Na caixa de dilogo escolha como camada de entrada Mosaico_Recortado a partir da lista suspensa. Para selecionar a camada de sada escolha de novo a subpasta Novos-TIF e digite o nome como Mosaico_Colorido. Clique em Salvar. De volta janela Relevo, em Cores do Relevo, clique em criar automaticamente e note que as cores so definidas pelo QGIS:

Observe que o limite valor mnimo de 0 e o valor mximo de 2476 so aqueles fornecidos na aba Metadados da caixa de Propriedades da Camada. Mas voc se lembra de que quando da produo do mapa hipsomtrico foi exibido o valor mnimo de 29.325 e o valor mximo de 1051.6, ou seja, valores que variam de 0 a 1100, mas para mantermos uma equidistncia de 200 metros podemos adotar 1200 como valor mximo. Podemos editar os valores dos limites e alterar tambm as cores. Para editar o limite de uma classe clique duas vezes sobre o valor numrico, isso far com que uma caixa de dilogo seja exibida onde possvel apagar o limite antigo e digitar o limite novo, depois clicar em OK.

Repita o mesmo procedimento tanto para o limite inferior quanto para o limite exterior de cada uma das classes, com equidistncia de 100 metros, obtendo ao final 6 classes: 0 200, 200 400, 400 600, 600 800, 800 1000, 1000 1200. Entretanto, note que h mais classes do que necessitamos, sendo assim, clique na classe excedente e depois clique em para suprimir a classe selecionada.

30

Em seguida, clique na cor associada com cada classe e escolha uma nova cor a partir da caixa de seleo de cor, evitando usar verde, pois normalmente associada com vegetao.

Clique em OK na caixa de dilogo Relevo. Arraste a camada Mosaico_Colorido para dentro do grupo CG_SIR. Clique com o boto direito do mouse sobre o nome da nova camada e escolha a opo aproximar a camada sua extenso.

Aplique zoom na rea do Pontal do Paranapanema para exibir mais detalhes da imagem.

31

Desselecione a camada Mosaico_Colorido e exiba agora a camada Mapa_Hipsomtrico, mantendo a mesma rea e mesmo nvel de zoom. Compare a figura mostrada acima com a figura abaixo.

22. Sobrepor camadas raster e alterar a transparncia

No item anterior demonstramos como ligar e desligar cada camada, exibindo uma de cada vez. Vamos agora exibir duas camadas raster simultaneamente para termos os benefcios que cada uma delas tem a oferecer. Verifique se a camada Mosaico_Recortado est acima da camada Relevo_Sombreado no painel de camadas. Caso no esteja clique sobre o nome da camada e arraste de acordo. Exiba apenas a camada Relevo_Sombreado.

32

Em seguida exiba tambm a camada Mosaico_Recortado. Apesar da camada Relevo_Sombreado estar habilitada seu contedo no pode ser visto porque est abaixo da outra camada. Exiba as propriedades da camada Mosaico_Recortado e na aba transparncia constate que o valor de 0%.

Altere o valor para 25% e clique em Aplicar. Observe que agora possvel ver as duas camadas simultaneamente. Esse artifcio acrescenta uma impresso de textura ao mosaico de imagens SRTM.

33

Feche a caixa de propriedades da camada Mosaico_Recortado. Repita o procedimento, mas desta vez exiba a camada Relevo_Colorido com 40% de transparncia, tendo abaixo dela a camada Relevo_Sombreado.

Feche a caixa de propriedades da camada. Repita o procedimento, usando desta vez a camada Mapa_Hipsometrico com transparncia de 10%.

34

Tendo finalizado esta etapa edite novamente as propriedades das camadas para que a transparncia seja de 0% de novo. Explore os recursos do QGIS, gere novas camadas alterando os parmetros, por exemplo, crie novos mapas de relevo sombreado com outros valores de azimute e ngulo vertical. Analise outras reas do Estado de So Paulo, alm dessa que demonstramos, localizada no Pontal do Paranapanema.

Lembre-se de vez em quando de clicar em Projeto->Salvar.

23. Produzir mapa de aspecto

O mapa de aspecto mostra a orientao das vertentes representadas por 8 setores com amplitude de 45, que indicam a direo do movimento da gua da chuva quando escorre sobre a superfcie do terreno. Na barra de menus clique em Raster e escolha a opo Anlise do Terreno, em seguida escolha Aspecto. Na caixa de dilogo podemos escolher tanto a camada com as elevaes quanto a camada de sada, sendo assim, escolha como camada de entrada Mosaico_Recortado, a partir da lista suspensa. Para selecionar a camada de sada escolha a subpasta Novos-TIF e digite o nome como Relevo_Aspecto. Clique em Salvar. De volta caixa Aspecto, mantenha os demais parmetros sugeridos pelo QGIS e clique em OK.

Arraste a nova camada para dentro do grupo CG_SIR. Exiba o resultado obtido e expanda a pasta para verificar os valores mnimo e mximo da camada, lembrando que o limiar mnimo teoricamente seria de 0 graus e o limiar mximo de 360 graus.

35

Amplie a imagem na rea do Pontal do Paranapanema.

Exiba as propriedades da camada e clique na aba Estilo, altere o tipo de renderizao de Banda simples cinza para Banda simples falsa cor, escolha a paleta Spectral e mantenha os parmetros sugeridos pelo QGIS. Clique em Classificar e depois em Aplicar. Analise o resultado e observe atentamente os valores numricos sugeridos pelo QGIS e as cores que foram associadas na imagem. Compare a figura mostrada acima com a figura mostrada abaixo neste Tutorial, tente estabelecer um vnculo com os pontos cardeais (Norte, Sul, Este e Oeste) e colaterais (Nordeste, Sudeste, Sudoeste e Noroeste).

36

Note que na imagem acima a cor vermelha est associada com os menores valores numricos (7.568010), enquanto que a cor azul est associada com os maiores valores (354.600000). A cor amarela est associada com os valores intermedirios (181.084005). Para gerar outra verso do Mapa de Aspecto, em Modo escolha Intervalo igual, defina 9 classes, para o valor mnimo digite 0 e o valor mximo digite 360, clique em Classificar e veja uma prvia da legenda. Altere a primeira e a ltima classe para que sejam representadas com cores iguais, altere as cores das demais classes para que tenham um bom contraste entre si. Altere os limiares superiores de cada classe, sendo 22.5 para a primeira classe (renomeie como Norte), 67.5 para a segunda (Nordeste), 112.5 para a terceira (Este), 157.5 para a quarta (Sudeste), 202.5 para a quinta (Sul), 247.5 para a sexta (Sudoeste), 292.5 para a stima (Oeste), 337.5 (Noroeste) para a oitava e 360 para a nona classe (Norte novamente). Clique em Aplicar. Note que a imagem resultante parecida com a anterior, mas h diferenas em funo do fatiamento que fizemos com a criao de 9 classes que se diferenciam em 45. A classe Norte aparece duas vezes, uma delas com valores de 0 a 22.5 e tambm de 337.5 a 360 graus. Clique em OK na caixa de propriedades.

Expanda a pasta da camada para exibir em detalhe as caractersticas das classes.

37

24. Inserir rosa dos ventos

Clique em Exibir->Decoraes->Rosa dos Ventos. Habilite a opo Mostrar rosa dos ventos, escolha a posio que preferir (por exemplo: inferior direito) e depois clique em OK.

Observe que agora a rosa dos ventos exibida na tela.

Clique em Projeto->Salvar. Feche o QGIS.

Novos tutoriais sero disponibilizados ao longo do tempo. Tendo dvidas sobre os procedimentos adotados e mesmo sugestes de

aprimoramento ou correes, entre em contato pelo email arletemeneguette@fct.unesp.br

38