Perimetros Urbanos Nespereira

Fonte: QGIS
Saltar para a navegação Saltar para a pesquisa

Perímetros urbanos de Nespereira

Objetivos

FIXME

Dados

Guião

Parte 1: preparar o projeto incluindo criar uma camada que represente Nespereira

  1. Criar a pasta base "Cinfaes_Nespereira", e as subpastas "dados_input", "dados_output" e "resultados"
  2. Descarregar os dados de input (cartografia de Cinfães e BD.gpkg)
  3. Descomprimir a cartografia de Cinfães e copiar para a pasta "dados_input" estas 7 shapefiles que representam 22 ficheiros:
    1. edificio (contém 46779 itens)
    2. Rodovias (contém 56715 itens)
    3. campo_jogos (contém 46 itens, neste caso não são 3 mas sim 4 ficheiros já que inclui o ficheiro .prj que especifica o sistema de coordenadas dessa camada)
    4. deposito (contém 68 itens)
    5. equipame (contém 104 itens)
    6. piscina (contém 5 itens)
    7. zverde (contém 100 itens)

    Verificar que na pasta "dados_input" ficaram 22 ficheiros que ocupam 43 MB

  4. Copiar a geopackage BD.gpkg para "dados_input" e mudar-lhe o nome para "BD_Cinfaes.gpkg", a seguir usar o navegador do QGIS para ver o conteúdo da geopackage e clicando com o botão direito do rato sobre o nome da camada (estamos a falar da camada armazenada no disco na geopackage), apagar ("Delete Layer...") as camadas "COS" (é uma camada matricial conforme mostrado pelo ícone parecido com um tabuleiro de xadrez), "lagoasAA", "precipAA", e "soloAA" já que não vamos necessitar delas
  5. Criar um novo projeto QGIS, verificar que o sistema de referência de coordenadas é "EPSG: 3763 ETRS89 / Portugal TM06" (menu "Projeto / Propriedades / Sistema de referência de coordenadas (SRC) do projeto"), a seguir gravar dentro da pasta "Cinfaes_Nespereira" usando o nome "Cinfaes_Nespereira.qgz"
  6. Adicionar ao projeto as camadas edificio, Rodovias, campo_jogos, deposito, equipame, piscina, zverde e verificar que a quantidade de itens de cada camada está correta (ver quantidades dadas acima)
  7. Adicionar o mapa de contexto do Open Street Maps (menu "Web / QuickMapServices / OSM / OSM Standard", se esse menu não existir é necessário instalar o módulo "QuickMapServices" indo ao menu "Módulos / Gerir e instalar módulos...")
  8. Verificar que todas as camadas parecem estar corretamente georreferenciadas em Cinfães
  9. Para todas as camadas em que aparece um ponto de interrogação à direita do nome da camada (todas exceto campo_jogos) significa que essa camada não tem definido o sistema de referência de coordenadas, embora por omissão seja assumido o sistema de referência de coordenadas do projeto (pelo que as camadas neste momento ficam bem georreferenciadas), se mais tarde criarmos uma nova camada a partir de uma dessas camadas (por exemplo usando "selecionar elementos por expressão") então a nova camada vai ficar com o SRC errado (ficará com o SRC standard WGS84) por isso é importante corrigir antes de se criar camadas derivadas. Clicar em cada um dos pontos de interrogação e escolher "EPSG: 3763 ETRS89 / Portugal TM06" e premir OK
  10. Criar a camada "ED01" que foi criada a partir da camada edificio excluindo os edifícios do tipo Anexo, Cobertura/Telheiro, Precário. Para se conseguir isso usar "selecionar elementos por expressão" do painel "Selection toolbar" podendo-se criar o que se quer usando por exemplo a expressão NOT ("tipo" IN ( 'Anexo' , 'Cobertura/Telheiro' , 'Precário' )), deve dar 39763 itens, a seguir usar o menu "Editar / Copiar Elementos" seguido de "Editar / Colar elementos como / Camada temporária" e colocar no "Nome da camada" o valor "ED01"
  11. Da mesma forma criar a camada "VIA01" a partir da camada "Rodovias" de modo a só ter as estradas principais, isto é, as que são do tipo "Estrada/rua" (descartando as do tipo estrada de terra batida, via de apoio, via para peões). Pode-se usar por exemplo a expressão "tipo" = 'Estrada/rua', deve dar 18361 itens. Note-se que embora seja óbvio que estrada de terra batida e via para peões não representa uma estrada principal o mesmo já não se aplica a via de apoio. Para ter a certeza que via de apoio não é uma estrada relevante deve-se selecionar esse tipo de vias (usando uma expressão) e a seguir fazer zoom no mapa de modo a ver o contexto, só aí se conclui que vias de apoio não representam estradas que devam ser consideradas (muitas vezes são por exemplo estradas que ligam a estrada principal a uma propriedade).
  12. Tornar as 2 camadas temporárias permanentes, guardando-as na subpasta "dados_output" numa geopackage chamada "Cinfaes_camadas.gpkg": clicar no ícone da memória que identifica a camada como sendo temporária, escolher o formato "Geopackage", no nome do ficheiro clicar nas reticências, definir a localização e o nome do ficheiro, a seguir escolher para nome da camada o mesmo nome que está no projeto (isto é, respetivamente "ED01" e "VIA01"), e premir OK
  13. No caso de por lapso se ter criado a camada no ficheiro geopackage correto mas se ter usado o nome da camada errado (fácil de acontecer porque por omissão o nome da camada é inicializado com o nome do ficheiro do geopackage) eis a forma de corrigir o nome da camada:
    1. Selecionar a camada, clicar com o botão direito do rato, "Exportar / Guardar elementos como..." guardar no mesmo geopackage mas agora com o nome correto e deixar ativa a opção "Adicionar ficheiro guardado ao mapa", irá aparecer uma nova camada com o nome do ficheiro geopackage, um espaço e o nome da camada
    2. Usando o navegador navegar até à geopackage, selecionar a camada com o nome errado, clicar com o botão da direita do rato, escolher "Eliminar Camada..." (estamos a eliminar a camada dentro da geopackage armazenada no disco), deve aparecer uma mensagem a dizer que a camada existe no projeto atual e a perguntar se quer removê-la do projeto e eliminá-la, selecionar "Yes" (dessa forma não só a apaga no disco como também a remove do projeto)
    3. No projeto alterar o nome da camada nova de modo a remover o prefixo com "o nome do geopackage e um espaço", dessa forma o nome no projeto fica igual ao nome no geopackage
  14. Guardar o projeto
  15. Adicionar ao projeto a camada "COS_vetorial" que se encontra no ficheiro "BD_Cinfaes.gpkg" que se encontra na pasta "dados_input"
  16. A partir da camada "COS_vetorial" criar a camada temporária "Cinfaes" que contém somente o concelho de Cinfães, usar por exemplo "selecionar elementos por expressão" usando a expressão "Concelho" = 'Cinfães', deve dar 14 itens
  17. A partir da camada "Cinfaes" criar a camada temporária "Nespereira" que contenha somente a freguesia de Nespereira (irá dar 1 item)
  18. Gravar as camadas temporárias "Cinfaes" e "Nespereira" na geopackage "dados_output/Cinfaes_camadas.gpkg"
  19. Gravar o projeto, sair do QGIS, e fazer uma cópia de segurança do resultado obtido até aqui.
Resultado da Parte 1

Arquivo com o resultado da Parte 1, o projeto e os dados: Cinfaes_Nespereira-20211025.zip (108 MB)

Tempo de execução (quando se tenta repetir pela segunda vez): cerca de FIXME minutos

Vídeo com a resolução do exercício: [FIXME FIXME.wmv] (FIXME minutos e FIXME segundos, FIXME MB)

Parte 2: Áreas edificadas consolidadas e em consolidação de Nespereira

  1. Reabrir o projeto ou usar o arquivo ZIP Cinfaes_Nespereira-20211025.zip com o resultado da Fase 1 (NOTA: se se usar esse arquivo ZIP ver-se-á que no projeto existe uma camada chamada "Cont_AAD_CAOP2020" que contém somente 5 itens em vez dos 3223 itens que se encontram na camada do mesmo nome na geopackage "dados_input/BD_Cinfaes.gpkg", trata-se de um lapso de uma camada que foi criada e não foi apagada, remover essa camada do projeto. Para além disso se aparecer uma pasta chamada "__MACOSX" removê-la, ela aparece quando se cria arquivos ZIP num computador da Apple)
  2. Adicionar à camada "ED01" um novo atributo chamada area_m2 do tipo número decimal (real) com a $area
  3. Criar a camada "ED02" com os elementos de "ED01" que possuem uma área maior que 30 m2 (selecionar pela expressão com area_m2 > 30), deve dar 23295 itens
  4. Criar a camada "EDN01" com os elementos de "ED02" que pertencem a Nespereira, para isso no painel "Selection Toolbar" selecionar o ícone mais à esquerda ("Selecionar por localização", atenção que este ícone não aparece se o módulo de "Processing" estiver desativado, esse módulo é desativado automaticamente quando se arranca o QGIS a seguir a uma operação que originou um crash do QGIS enquanto um algoritmo do módulo "Processing" estava a ser executado, isso é um problema que é frequente acontecer). Na janela de "Selecionar por localização" escolher:
    1. Selecionar elementos de: ED02
    2. Where the features (geometric predicate)": intersecta
    3. Ao comparar com os elementos de: Nespereira

    O resultado é que esses elementos ficam selecionados, a seguir copiar e colar como uma nova camada temporária chamada "EDN01" que deve ter 2715 itens, guardar essa camada no geopackage "cinfaes-camadas.gpkg" numa camada com o mesmo nome

  5. Criar uma camada tampão sobre EDN01 com 25 metros e agregada num único item: menu "Vetor / Ferramentas de geoprocessamento / Buffer...", escolhendo:
    1. Camada de entrada: EDN01
    2. Distância: 25 metros
    3. Ativar a checkbox "Resultado da agregação"

    Deve resultar num único elemento, mudar o nome da camada resultante de "Buffered" para "BUFF01"

  6. Só para testar, fazer o mesmo mas desta vez sem ativar a checkbox "Resultado da agregação", comparar o resultado, agora temos 2715 itens em vez de 1, apagar esta camada acabada de criar com 2715 itens
  7. Unir áreas com uso urbano estabilizado:
    1. Juntar as camadas vetoriais "campo_jogos, deposito, equipame, piscina, zverde" usando o menu "Vetor / Ferramentas de gestão de dados / Juntar camadas vetoriais...", aceitando o nome da camada dado por omissão "Juntado", deve dar 323 itens
    2. Agregar os itens da camada "Juntado" usando o menu "Vetor / Ferramentas de geoprocessamento / Agregar...", aceitar o nome por omissão de "Agregado", tem de obviamente dar 1 único elemento (um multi polígono)
    3. Na camada "Agregado" converter esse multi polígono usando o menu "Vetor / Ferramentas de geometria / Múltiplas partes para partes únicas..." criando a camada "PartesSimplesCinfaes".
      NOTA: este algoritmo pega numa camada vetorial com geometrias multi-partes e gera uma nova camada vetorial na qual todas as geometrias contêm uma parte única. Os elementos com geometrias multi-parte são divididos em tantos elementos diferentes quantas as partes que a geometria contém, e são utilizados os mesmos atributos para cada um deles. Neste caso o resultado são 254 itens (polígonos simples).
    4. Intersetar "PartesSimplesCinfaes" com Nespereira (menu "Vetor / Ferramentas de geoprocessamento / Interseção..."), criando a camada "USOURB01", deve dar 13 itens
      NOTA: este algoritmo interseta 2 camadas e é simétrico, tanto faz intersetar a camada "PartesSimplesCinfaes" com a camada "Nespereira" como fazê-lo na ordem inversa
  8. Unir áreas de uso urbano com buffer do edificado maior que 30 m2 separado de um máximo de 50 m:
    1. Juntar as camadas vetoriais "USOURB01" e "BUFF01" (menu "Vetor / ferramentas de gestão de dados / Juntar camadas vetoriais...") de modo a criar a camada "USOURB02"
    2. FIXME

Resultado

Ficheiros com o resultado:

  • [FIXME FIXME] (FIXME, FIXME MB)

Arquivo com o resultado, o projeto e os dados: [FIXME FIXME.zip] (FIXME MB)

Estrutura do arquivo FIXME.zip:

  • FIXME

Tempo de execução (quando se tenta repetir pela segunda vez): cerca de FIXME minutos

Vídeo com a resolução do exercício: [FIXME FIXME.wmv] (FIXME minutos e FIXME segundos, FIXME MB)

Obtida de "https://ccti.ismai.pt/qgis/index.php?title=Perimetros_Urbanos_Nespereira&oldid=376"