Manipulação de dados no SPSS e QGIS

Preview:

Citation preview

Prática de Manipulação de Dados no SPSS e Qgis

Vitor Vieira Vasconcelos

BH1350 – Métodos e Técnicas de Análise da Informação para o Planejamento

Junho de 2016

O que nós vimos na aula passada

Inferência estatística

Normalidade de dados

Gráficos de caixas e bigodes (boxplot)

Valores discrepantes (outliers)

Gráficos Q-Q

Transformações de dados

Conteúdo SPSS

Prática com dados da Amostra do Censo Editor de Sintaxe

Ponderação e Agregação

Qgis

Comparação entre Censo 2000 e 2010 Intercessão e ponderação por área

Dados de chuva das estações para os setores censitários Interpolação de pontos para raster

Estatísticas Zonais (raster para poligono)

Declividade de ocorrências de deslizamento Amostra de raster para pontos

Contagem e média de atributos de pontos em polígonos

Proximidade de ocorrências

Dados da AmostraResultados Gerais da Amostra:

ftp://geoftp.ibge.gov.br/recortes_para_fins_estatisticos/malha_de_areas_de_ponderacao/censo_demografico_2010/

Dados da AmostraMalha em shapefile:

ftp://geoftp.ibge.gov.br/recortes_para_fins_estatisticos/malha_de_areas_de_ponderacao/censo_demografico_2010/

Formato em ASCII

Script para formatação dos microdadoshttp://metodologiapolitica.com/leitura-dos-microdados-ibge-spss/#sthash.jUGqmTB7.dpbs

Janela de sintaxe No Editor de Dados

Arquivo-> Novo -> Sintaxe ->

No Editor de Sintaxe

Arquivo-> Abrir -> Sintaxe

Arquivos do tipo: “todos os arquivos”

Abrir: Leitura-SPSS-Censo-2010-domicílios.txt

Local e nome do Arquivo

Selecionando SBC e criando novo conjunto

Dados-> Selecionar Casos

Alterar V002 para tipo Numérico

Abrir arquivo “Layout_microdados_amostra.xls”

Transformação de dados Transformar-> Calcular Variável

Dividir V0010 por 10000000000 (13 casas decimais)

Obs: qualquer outra variável da amostra que tenha casas decimais precisaria ser transformada também

Ponderação dos casos Dados-> Ponderar Casos

V0201 – Condição de OcupaçãoV2011 – Aluguel

Verificar Atributos de Tipo de Medida das Variáveis

Agregar dados

Dados-> Agregar

Funções de AgregaçãoPeso V0201

V2011

Dados Agregados

Salvar como Dbase IV: sbc_agregado.dbf

União com Shapefile no Qgis Camada-> Adicionar Camada ->vetorial->

35SEE250GC_SIR_area_de_ponderacao.shp

Camada-> Adicionar Camada -> vetorial-> sbc_agregado.dbf

União com Shapefile no Qgis Clique com botão direito na camada:

35SEE250GC_SIR_area_de_ponderacao.shp

Propriedades-> Uniões

Relacionar V0010 (de sbc_agregado) com AREA_POND

Selecionar e salvar áreas de ponderação

Clique duas vezesno atributo unidode casos, para ordenar

Selecione as áreasde ponderaçãocom casos

Clique com botãodireito na camada:

Salvar como-> Salvar somente

feiçõesselecionadas

Validar Geometria

Vetor-> Geometrias-> Verificar a validade da Geometria

Limpar Geometria

Processar-> Caixa de Ferramentas-> Comandos GRASS GIS 7-> vector-> v.clean

Comparação entre Censos de 2000 e 2010

Criar novo projeto

Abrir camadas com população por setores:

sbc_setores_2000_pop.shp (Censo 2000 – urbano)

sbc_setores_2010_pop.shp (Censo 2010)

Vamos manter sempre a projeção UTM

Algoritmos necessitam das camadas na mesmaprojeção

Projeção UTM é necessáriapara cálculo de área emmetros

Ao salvar uma camada:

Clique com botão direito nacamada

Salvar como->

Alterar SRC para

EPSG 29193 SAD69 UTM 23S

Vamos manter sempre a projeção UTM

Manter a projeçãodo projeto sendo a mesma dos shapefiles

Arquivo->Propriedades do Projeto -> SRC

Calculando areas Abra a tabela de atributos de

“sbc_setores_2010_pop.shp”

Clique no icone “Calculadora de Campos”

Calcule $area emum novo campo

Calculando densidade populacional Abra a tabela de atributos de

“sbc_setores_2010_pop.shp”

Clique no icone “Calculadora de Campos”

Calcule

(“pop” / $area) * 1000000

em outro campo “dens_km2”

Crie o campo de áreae de densidadedemográfica para a camada de 2000

Interseção de Polígonos1. Vetor->Ferramentas de

Geoprocessamento -> Intersecção

2. Processar->Caixa de Ferramentas->

1. Comandos Grass-> Vector -> v.overlay

2. Geoalgoritmos-> Ferramentas de sobreposição de vetor-> Intersecção

3. Saga-> Vector Polygon Tools-> Interseção

Interseção de Polígonos1. Vetor->Ferramentas de

Geoprocessamento -> Intersecção

2. Processar->Caixa de Ferramentas->

1. Comandos Grass-> Vector -> v.overlay

2. Geoalgoritmos-> Ferramentas de sobreposição de vetor-> Intersecção

3. Saga-> Vector Polygon Tools-> Interseção

Não

funcionam se

houver erro

de geometria!

Criar ambiente GRASS no QGIS Instalar complemento GRASS 7

Complemento-> GRASS-> Novo Mapset

1

2

3

Criar ambiente GRASS no QGIS Complemento-> GRASS-> Novo Mapset

4 5

6

Importando shapefile no GRASS

• Botão Ferramentas do GRASS

• Complementos-> GRASS -> Abrir Ferramenta do GRASS

• Gerenciador de arquivo -> Importar para o GRASS -> v.in.ogr.qgis

Importando shapefileno GRASS

• Snapping threshold: 1

• Minimum size area: 5

• Fazer o mesmo para as

duas camadas

Interseção no GRASS• Vetor-> Análise espacial-> Sobrepor-> v.overlay.and

• Rodar e então “Ver Saída”

Exportar para Shapefile• Gerenciador de arquivo-> Exporta a partir do GRASS->

Exporta vetor a partir do GRASS -> v.out.ogr

Interpolação de Dados Instalar Complemento de Interpolação

Abrir camada: pluviometricas_sbc_utm.shp

Interpolação de Dados

Raster-> Interpolação

Coeficiente P = 4

Definir pela extensão atual

Interpolação de Dados

Interpolação de Dados

Processar-> Caixa de Ferramentas-> Comandos GRASS GIS 7-> Vector-> v.surf.rst

Interpolação de Dados

Extração Zonal de Atributos

Instalar complemento Estatística por Zona

Raster-> Estatística por Zonas

Amostragem e Contagem de Pontos

Instalar complemento“Point Sampling Tool”

Adicionar camadas:

Vetorial: ocorrencias_sbc_utm.shp

Raster: srtm_sbc_utm_declividade.tif

Amostragem e Contagem de Pontos

Complementos-> Analyses-> Point Sampling Tool

Amostragem e Contagem de Pontos

Vetor-> Analisar-> Pontos no Polígono

Proximidade de pontos

Criar um atributo de valor “1” na camada“ocorrencias_sbc_utm.shp”

Gravar alterações na camada e fechar edição

Proximidade de pontos Processar-> Caixa de Ferramentas-> SAGA-> Raster Creation Tools->

Rasterize

Proximidade de pontos Raster-> Análise-> Proximidade (Distância Raster)

Proximidade de pontos Raster-> Estatística por Zonas

Próxima Aula

Apresentação:

Proposta de trabalho

Base territorial

Análise exploratória