Upload
vitor-vieira-vasconcelos
View
228
Download
13
Embed Size (px)
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