operação "recortar" (clip) e modificações de atributos do shapefile

[gavi...@gmail.com]

      Prezados, estou com uma dificuldade na operação "clip" - recortar.

  Estou trabalhando com os dados do Cadastro ambiental rural (www.car.gov.br), e quero separar os shp de vegetação nativa que incidem sobre os shp de propriedades de tamanhos diferentes.....iniciei a operação com "intersecção", depois tentei "extrair por localização", "overlay", e foi o algoritmo "recortar/cliping" o que surtiu melhor resultado...ocorre que na feição "recortada/clipada", os atributos (área do fragmento florestal) não corresponde com a realidade - no shp o fragmento incidente numa determinada propriedade tem 0,5 hectare, mas nos atributos ele transportou o atributo 'original" que é referente a área muito maior...

  alguém sabe solucionar isto?

  Att

[Instrutor GIS]

A ferramenta que realiza o transporte de atributos é o intersect, não o clip. Aliás,  essa é a única característica que diferencia uma ferramenta de outra, a saber, a capacidade de preservar os atributos das duas camadas envolvidas.

Melhora você rever seus procedimentos,  parece que você está usando a ferramenta errada.

--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para qgisb...@googlegroups.com.
Acesse esse grupo em https://groups.google.com/group/qgisbrasil.
Para mais opções, acesse https://groups.google.com/d/optout.

[gavi...@gmail.com]

Prezado, como está/

  Obrigado pela informação e pronta resposta!  O "intersect" deu certo em uma operação, depois acusou o erro 

Erro no geoprocessamento GEOS: Uma ou mais feições de entrada têm geometria inválida

         Eliminei os fragmentos dos poligonos, mas não deu certo....o que pode estar havendo?

[Kyle Felipe]

O problema de trabalhar com o CAR é Erros na geometria...
Tem muito erro, e tem de fazer a correção antes de iniciar o trabalho...

[Felipe Gavioli]

Olá Kyle!

  Esta correção seria aplicar o "validar geometria" nos shps? Ou voce recomenda transformar o poligono em linha e corrigir eventuais aberturas, etc.?

--

Comunidade QGIS Brasil: http://qgisbrasil.org/
---

Você recebeu essa mensagem porque está inscrito em um tópico no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse tópico, acesse https://groups.google.com/d/topic/qgisbrasil/yEPX81BWpMc/unsubscribe.
Para cancelar inscrição nesse grupo e todos os seus tópicos, envie um e-mail para qgisbrasil+...@googlegroups.com.

[Kyle Felipe]

aplicar o validar geometria....

[Felipe Gavioli]

Usei o "verificador de topologia" e de fato acusou erros de sobreposição dos polígonos....agora tô batendo cabeça com o v.clean para tentar corrigir

   Att

[Jorge Santos]

a operação break do v.clean deve quebrar os polígonos nas sobreposições.

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.

Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.

[Felipe Gavioli]

Valeu Jorge, vou testar

[Felipe Gavioli]

Jorge, de fato o "break" acaba com as sobreposições, mas cria um outro problema - duplica uma série de polígonos....na ferramenta "v.clean" estou usando "rmdupl" (remove duplicated), mas não tem dado resultado satisfatório....estou tentando via "verificar geometrias", para ver se limpo os polígonos....alguém já passou por experiencia semelhante?

Em qui, 20 de set de 2018 às 08:56, Jorge Santos <jorgepsa...@gmail.com> escreveu:

[Jorge Santos]

Você precisa usar o v.clean.advanced para determinar mais de um método. Eu tenho um tutorial em algum lugar perdido.

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

[Felipe Gavioli]

Poisé, encontrei este tutorial que achei muito bom (em inglês)

https://www.geodose.com/2017/10/how-to-clean-topology-error-in-qgis.html

[Jorge Santos]

É esse mesmo. É o melhor tutorial sobre o v.clean que vi até hoje. 

[Felipe Gavioli]

Colegas, 

  De fato a base de dados do CAR tem MUITOS erros......e a ferramenta v.clean.advanced não é capaz de corrigir.....chega determinado ponto que a função "bpol" não limpa mais sobreposições......tentei também o "verificar geometria", mas a cada erro reparado são gerados outros.....tá complicado...alguém tem experiência com tratamento dos dados do CAR:?

[Kyle Felipe]

 Os dados do car estão em (eka) SHP??
Se estiverem, coloque eles em um geopakage/spatialite/postgis e tenta reparar por lá.... esses formatos possuem mais ferramentas para se trabalhar do que o shp (que já tem 28 anos que não recebe melhorias)

[Felipe Gavioli]

Sim, os dados do car estão em shp. ok, vou tentar esta estratégia, obrigado

[Felipe Gavioli]

Kyle, deixa eu te perguntar:

  No QGis já consegui salvar os arquivos em formato geopakage.....você me sugere tentar o reparo das sobreposições, duplicações etc com as funções v.clean.advanced do grass ou mesmo "reparar geometria" (aplicadas no aquivo .gpkg) ou há algum complemento mais específico para tratar estes dados?

  Obrigado

   

[Instrutor GIS]

Esses dados são públicos? Se sim, pode subir eles para a nuvem?

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

[Kyle Felipe]

@Felipe Gavioli Eu nunca usei o v.clean, somento o "reparar Geometria"...
O problema é que, erros como sobreposição são muito comuns no CAR (MUITO MESMO, em alguns casos chegam a 100% da propriedade menor) e se remover a feição sobreposta, perde-se o dado. e dificil definir quem tem o direito sobre a área definida lá dentro.
Já os problemas de geometria são mais fáceis de se corrigir, são os primeiros que eu tento corrigir com as ferramentas.....

A grande quantidade de poligonos torna muito custoso o processo de correção, mesmo Desktops mais parrudas levam muito tempo para corrigir, o legal seria fazer por partes.

Tive o mesmo problema com as unidades de conservação aqui de minas e com os limites municipais...

Assim que eu tiver um tempo quero testar a V.CLEAN.

@Jorge Santos, sim, os dados são públicos, tem um site onde pode-se baixá-los por municipio...

[Felipe Gavioli]

Fala Jorge, sim sim, os dados são públicos

[Felipe Gavioli]

Sim, isto é fato.....há muitos polígonos envolvidos, tentei reparar tudo de uma vez e o QGIS trava, é inviável....

    O V.Clean funcionou para limpar as sobreposições, mas dai ele duplica os polígonos.....quando eu removo, via v.clean os poligonos duplicados, as sobreposições voltam......não achei a ferramenta boa (ou não estou usando direito, tem esta possibilidade também). 

[Instrutor GIS]

Envie o link para eu verificar os erros e fazer um vídeo comas correções;

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

[Instrutor GIS]

Eu vejo que esta é uma operação impossível para o Geoprocessamento por conta dessa base de imóveis do INCRA. Há uma infinidade de sobreposições, pontas de polígonos, geometrias inválidas de toda a sorte, etc. Um trabalho impossível de ser concluído se você está realizando intersecções com esse grupo de camadas. Agora, se você admitir a análise espacial considerando processar a sobreposição de cada imóvel individual, aí são outros quinhentos (vai saber o que isso significa).

Como são limites de imóveis, eles não podem ser corrigidos porque os validadores vão modificar os limites. Por incrível que pareça, esse dado do INCRA não pode ser arrumado.

Minha ideia inicial é seguir esse fluxo de trabalho:

1 - Corrigir todos os erros topológicos que por acaso possam existir na vegetação;

2 - Tentar transformar a vegetação em geometria multiparte;

3 - Separar, um a um, todos os imóveis do INCRA;

4 - Executar o Geoprocesso Intersect em lote entre o imóvel do tiozinho e toda a vegetação multiparte. Os inputs são a vegetação e cada um desses imóveis horrorosamente "organizados".

Vou verificar se isso vai dar certo.

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

Em sex, 28 de set de 2018 às 13:25, Felipe Gavioli <gavi...@gmail.com> escreveu:

Fala Jorge!

   Seguem os shp em anexo. O "imovel_agregado" são os shps dos poligonos do imóveis rurais do CAR para um grupo de municipios paulistas.......o "vegetação nativa" agregado são os shps dos poligonos de vegetação nativa, para os imóveis destes mesmos municipios paulistas...

    minha intenção é separar os imóveis por grupos de tamanhos (minifundios, pequenas , médias e grandes propriedade), e cruzar estas categorias com o shp de vegetação nativa, para avaliar quanto de vegetação incide em cada grupo de imóveis destes.....a separação por categorias consegui fazer via "extração por atributo" tranquilo.....ocorre que na intersecção do grupo "minifundios" com a "vegetação nativa" o QGIS acusa erro, e foi dai que verifiquei erros de topologia....

   Muito obrigado pelo auxilio!  A propósito, dei uma olhada no seu site, excelente conteúdo!  vou divulgar entre conhecidos, ]

   Att

  OBS. os shp eu coletei em www.car.gov,.br 

[Manoel Guedes]

Fiz a formação do CAR pelo  Secretaria do Meio Ambiente do Ceará - SEMACE  e digo logo, esse cadastro visa qualquer coisa, menos precisão e veracidade...

[Felipe Gavioli]

Olá Jorge, 

 

  Deixa eu ver se entendi:

1. A base do CAR que te passei apresenta erros impossíveis de serem corrigidos, pq a correção implicaria em descaracterizar a informação. 

2. Assim, a solução que vc propõe é intersectar não em grupos, mas considerando cada poligono de imóvel individualmente, pq assim seria possível "isolar temporariamente" e durante a operação da intersecção, estas sobreposições existentes

     Valeu jorge, vou tentar caminhar por esta via que voce indicou também , obrigado

[Instrutor GIS]

Perai que já estou fazendo na prática e você vai entender minha estratégia doida.

Esse é o procedimento para corrigir a vegetação: https://drive.google.com/file/d/1pF_gPcc0S2WzdbKHRcS2SyrXd9via384/view?usp=sharing

Tô fazendo a separação e o intersect vem depois,

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

[Instrutor GIS]

O processo deu certo. É duro, mas parece que cada um mede o seu quintal à sua forma. O que tem na minha cabeça é que o INCRA fornece uma base agregada, mas não é para que façamos nada com ela, ou seja, nada no sentido de análises espaciais. Qualquer procedimento precisa ser executado nos imóveis individuais. O conjunto inteiro fica mesmo zuado. Abaixo estão os vídeos.

Como eu já expliquei, primeiro eu arrumo a vegetação com dissolve e fixar geometrias:

https://drive.google.com/file/d/1pF_gPcc0S2WzdbKHRcS2SyrXd9via384/view?usp=sharing

Em seguida, eu crio um campo para o QGIS 3 dividir a camada por atributos, para em seguida usar o Renomear Tudo que é free para limar o prefixo.

https://drive.google.com/file/d/1XtCtiSkazh4SXwckK9NVvclNZOzmI6jA/view?usp=sharing

Por último, eu executo o intersect e corto a vegetação nos limites de cada propriedade.

https://drive.google.com/file/d/1_bz6J5fSsdmxpJLDvXIIoK5ngOG3EpaK/view?usp=sharing

Este foi o melhor caminho que encontrei. Sugestões e melhorias são bem-vindas. O QGIS 3 está respondendo muito bem a essa complicação de comandos.

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002

Em sex, 28 de set de 2018 às 15:13, Felipe Gavioli <gavi...@gmail.com> escreveu:

[Felipe Gavioli]

Valeu Jorge!!  Vou estudar os videos e trabalhar com base nesta sua metodologia!! Muito obrigado!!

  Só p enriquece o debate - a Embrapa fez um trabalho tambem com estes mesmos dados do CAR, e disponibilizou o roteiro metodológico neste link - https://www.embrapa.br/car/material-metodos

Mais uma vez obrigado

[Felipe Gavioli]

Excelentes videos, muito bom, uma baita aula!!

[Felipe Gavioli]

Achei curioso vc ter usado o intersect pelo SAGA Gis e não o caminho convencional dentro do QGIS mesmo

[Instrutor GIS]

Os Geoprocessos do SAGA GIS são mais rápidos, funcionam de forma mais eficiente e raramente dão pau.

No Curso Avançado de QGIS eu recomendo o uso deles.

Jorge Santos
Instrutor de Geotecnologias
http://www.instrutorgis.com.br
jorgep...@outlook.com
WhatsApp: (61) 99616-5665
Skype: jorgepsantos2002
https://www.linkedin.com/in/jorgepsantos
Twitter: http://twitter.com/jpsantos2002