A secretaria tem dinheiro para reformar poucas escolas. Em quais municípios começar?

Montar as pastas e carregar os dois arquivos

Antes de abrir o QGIS, organize o disco. Projeto de GIS guarda caminhos, não arquivos: se você mover ou renomear um dado depois, o projeto abre com a camada quebrada.

1. Crie a pasta InfoGeo/aps01 e, dentro dela, dados_brutos e dados_tratados.
2. Coloque em dados_brutos os três arquivos do seu estado: o .csv, o .csvt de mesmo nome e o .gpkg. O .csvt não abre no QGIS e não aparece em lugar nenhum da interface; ele só precisa estar ali.
3. Abra o QGIS. Em Project ▸ New, e depois Project ▸ Save As, salve como InfoGeo/aps01/projeto.qgz. Salve de novo, com Ctrl+S, a cada etapa concluída.
4. Arraste municipios_<UF>_2022.gpkg do gerenciador de arquivos para dentro da janela do QGIS. A camada aparece no painel Layers, à esquerda, e o mapa desenha o contorno do estado.
5. Olhe o canto inferior direito da janela: deve estar escrito EPSG:4674. Esse é o sistema de coordenadas que o projeto assumiu da primeira camada carregada.
6. Agora importe a tabela de escolas: Layer ▸ Add Layer ▸ Add Delimited Text Layer…. Abre a janela Data Source Manager, já na aba certa.

Para renomear qualquer camada, clique nela no painel Layers e aperte F2. Nomes curtos e descritivos economizam confusão nos passos seguintes, porque você vai terminar a APS com seis ou sete camadas na tela.

Conferir a chave dos dois lados antes de qualquer conta

O cruzamento do passo 5 vai casar o código de município da tabela com o da malha. Cinco minutos aqui evitam refazer tudo depois.

1. Clique na camada escolas no painel Layers e aperte F6. Abre a tabela de atributos. O título da janela mostra o total de feições.
2. Localize a coluna CO_MUNICIPIO e confira que os valores têm sete dígitos.
3. Feche a tabela, clique na camada de municípios e aperte F6 de novo. Faça a mesma conferência em CD_MUN.
4. Agora o tipo. Clique duas vezes na camada escolas para abrir Layer Properties e vá na aba Fields, na lista à esquerda. A coluna Type mostra o tipo de cada campo.
5. CO_MUNICIPIO precisa aparecer como String ou text. Repita na camada de municípios: CD_MUN também deve ser String.

Agregar as escolas por município

A tabela tem uma linha por escola e você precisa de uma linha por município. Quem faz isso é o Statistics by categories: ele agrupa as linhas pelas categorias de um campo e devolve contagem, soma, média, mediana e outras estatísticas do campo de valores.

1. Abra a caixa de ferramentas com Ctrl+Alt+T, ou por Processing ▸ Toolbox. O painel Processing Toolbox abre à direita.
2. Na caixa de busca no topo do painel, digite statistics by categories. A árvore filtra enquanto você digita.
3. Clique duas vezes no resultado para abrir a janela da ferramenta. O nome completo, qgis:statisticsbycategories, aparece no rodapé da janela.
4. Preencha os campos conforme o quadro abaixo.
5. Clique em Run. A aba Log mostra o andamento e, no fim, o tempo de execução.
6. Feche a janela com Close. A camada nova está no painel Layers.

Você vai rodar essa ferramenta mais de uma vez: uma sobre todas as escolas, para ter o denominador, e outra sobre o recorte que interessa à pergunta, para ter o numerador. O passo 4 mostra como fazer o recorte.

Recortar a tabela por uma condição

Para agregar só parte das escolas, filtre a camada antes de rodar a ferramenta. O filtro é uma expressão que fica ativa na camada: o QGIS passa a enxergar apenas as linhas que atendem à condição, e as ferramentas de processamento respeitam esse recorte.

1. Clique na camada escolas e aperte Ctrl+F, ou clique com o botão direito e escolha Filter…. Abre a janela Query Builder.
2. A lista Fields, à esquerda, tem os nomes das colunas. Clique duas vezes num nome para inserir na expressão, já com as aspas corretas.
3. Para ver que valores um campo assume, selecione-o e clique em All, embaixo da lista Values. Útil para conferir se LOCALIZACAO usa 1 e 2, por exemplo.
4. Escreva a condição na caixa Provider specific filter expression, embaixo.
5. Clique em Test. Aparece uma mensagem com quantas linhas atendem à condição. Se der zero, a expressão está errada, e ainda dá tempo de corrigir.
6. Confirme com OK.
7. Com o filtro ativo, rode o Statistics by categories de novo, gravando com outro nome de arquivo.
8. Depois, volte em Ctrl+F, clique em Clear e confirme com OK. O funil desaparece e a camada volta inteira.

A sintaxe das expressões é a mesma do tema 3:

"ESGOTO_ADEQ" = 0                        -- número sem aspas
"LOCALIZACAO" = 2                        -- código: ver o LEIA-ME
"DEPENDENCIA" IN (2, 3)                  -- lista de valores
"ESGOTO_ADEQ" = 0 AND "QT_MAT_BAS" > 0   -- duas condições
"NO_MUNICIPIO" LIKE 'Santa%'             -- texto entre aspas simples

Cruzar as tabelas com a malha

Agora as tabelas agregadas ganham geometria. O Join attributes by field value copia colunas de uma tabela para outra usando um campo em comum. Você vai rodar duas vezes, uma para cada agregação, encadeando: a saída do primeiro cruzamento é a entrada do segundo.

1. Na caixa de busca do Processing Toolbox, digite join attributes by field value e abra a ferramenta com duplo clique.

O prefixo importa. As duas tabelas agregadas trazem colunas com o mesmo nome, count e sum. Sem prefixos diferentes, a segunda sobrescreve a primeira ou entra com nome duplicado, e depois você não sabe qual é qual.

Deixe Discard records which could not be joined desmarcado nos dois cruzamentos. Município sem correspondência precisa continuar na camada, com nulo na coluna, para que você possa contá-lo e decidir o que fazer com ele.

Construir o indicador com o Field Calculator

Campo calculado é uma coluna nova cujo conteúdo vem de uma expressão sobre as colunas que já existem. É onde o indicador nasce.

1. Abra a tabela de atributos da camada cruzada com F6.
2. Na barra de ícones da tabela, clique no ábaco, cujo nome ao passar o mouse é Open Field Calculator. Com a tabela em foco, o atalho é Ctrl+I.
3. A camada entra automaticamente em modo de edição, e o ícone de lápis fica pressionado. Ela vai continuar assim até você salvar.
4. Preencha o alto da janela conforme o quadro abaixo e escreva a expressão na caixa grande.
5. Confirme com OK. A coluna nova aparece no fim da tabela.
6. Salve as alterações: clique no ícone de disquete da barra da tabela, Save Edits, ou aperte Ctrl+S.
7. Saia do modo de edição clicando no lápis. Se ele continuar pressionado, o cálculo ainda não está gravado no arquivo.

A parte de baixo da janela tem três colunas. À esquerda, a caixa onde você escreve a expressão. No meio, a lista de funções, com busca no topo e um grupo Fields and Values com os campos da camada, que entram na expressão com duplo clique. À direita, a ajuda da função selecionada.

A sintaxe é a mesma dos temas 3 e 4:

"TOT_count"                              -- campo entre aspas duplas
"SEM_count" / "TOT_count"                -- razão entre colunas
100.0 * "SEM_count" / "TOT_count"        -- percentual
round(100.0 * "SEM_count" / "TOT_count", 1)   -- uma casa
if("SEM_count" IS NULL, 0, "SEM_count")  -- trata o nulo
"SEM_sum" / ("AREA_KM2" / 100)           -- parênteses mudam a ordem

Ordenar, ler o topo e levar as listas para o documento

1. Na tabela de atributos, clique no cabeçalho da coluna do indicador. Um clique ordena crescente; clique de novo para decrescente. A setinha no cabeçalho indica o sentido.
2. Clique no número da primeira linha, à esquerda, segure Shift e clique no número da décima. As dez ficam destacadas em amarelo, na tabela e no mapa.
3. Ctrl+C copia as linhas selecionadas com o cabeçalho das colunas. Cole direto no Word ou no Excel.
4. Se preferir gravar em arquivo, feche a tabela e vá em botão direito na camada ▸ Export ▸ Save Selected Features As…, escolhendo Comma Separated Value como formato.
5. Para limpar a seleção antes do próximo critério, use Ctrl+Shift+A.

Repita para cada critério que a Pergunta 1 e a Pergunta 2 pedem, guardando cada lista com um nome que identifique o critério. Você vai precisar das três lado a lado para comparar.

Somar a área dos municípios escolhidos

Somar área em coordenadas geográficas não dá resultado em quilômetro quadrado, porque grau não é unidade de comprimento. Antes de medir, o dado precisa estar num sistema projetado adequado para área, como você viu no tema 2.

1. Vá em Project ▸ Properties ▸ General e, no bloco Measurements, escolha o elipsoide e deixe Area units em Square kilometers. É essa configuração que o $area vai usar.
2. Selecione os municípios que você vai recomendar. Na tabela, com Ctrl ou Shift nos números das linhas; ou por expressão, com Ctrl+F3, que abre Select Features by Expression.
3. Vá em Vector ▸ Data Management Tools ▸ Reproject layer….
4. Em Input layer, escolha a camada. Para usar só a seleção, clique no ícone à direita da lista suspensa e marque Selected features only.
5. Em Target CRS, escolha o sistema projetado adequado para medir área. Grave a saída em arquivo e clique em Run.
6. Na camada reprojetada, abra o Field Calculator e crie um campo do tipo Decimal number (real) com a expressão $area.
7. Abra o painel de estatísticas em View ▸ Panels ▸ Statistics. No alto do painel, escolha a camada e o campo. Leia a linha Sum.

A camada traz também AREA_KM2, que é a área oficial publicada pelo IBGE. Comparar o seu número com o oficial é uma boa conferência, e a Pergunta 3 pede que você declare em que sistema mediu.

Montar o coroplético

1. Clique duas vezes na camada no painel Layers para abrir Layer Properties, e vá na aba Symbology, a segunda da lista à esquerda.
2. Na lista suspensa do topo, troque Single Symbol por Graduated. A janela muda inteira.
3. Em Value, escolha o campo do indicador. Em Color ramp, uma rampa sequencial: cor clara para valor baixo, cor escura para valor alto.
4. Escolha o Mode e o número de Classes, e clique em Classify, embaixo da lista. As faixas aparecem na tabela do meio.
5. Para definir as faixas você mesmo, troque o Mode para Manual e clique duas vezes num valor da coluna Values para editar o limite. O limite superior de uma classe é o inferior da seguinte.
6. Clique duas vezes no texto da coluna Legend para reescrevê-lo. É esse texto que sai na legenda do mapa, e ele deve dizer o que a cor significa, não repetir o intervalo numérico bruto.
7. Confirme com Apply para ver no mapa sem fechar a janela, e depois OK.