QGIS calcula o raio de um arco para atributos

Eu tenho o QGIS 2.18.16 instalado. Meu problema é: como calcular o raio de um recurso de arco / curva (linha)? As linhas são originadas do arquivo .dgn, que li com o FME e escrevi no PostGis DB. Criei chaves primárias, etc., e a tabela é totalmente editável. Eu calculei comprimentos para os arcos com sucesso, mas não consigo descobrir como calcular o raio para essas curvas. Este ( https://www.mathopenref.com/arcradius.html ) pode dar uma ideia do lado da matemática, embora eu não consigo descobrir como usá-lo na calculadora QGIS. Então, eu preciso calcular o "R" para todas as linhas que tenho no banco de dados.

Abaixo está um exemplo. Meus dados têm arcos / curvas em uma tabela diferente das polilinhas "normais".

qgis field-calculator calculate-values Sisuaski
fonte

Quando clico em um arco com a ferramenta Identify, um dos atributos derivados listados é o "raio do vértice mais próximo". Portanto, não há necessidade de calcular o raio, você só precisa descobrir como acessar essa propriedade na calculadora de campo. Ou simplesmente copie-o da ferramenta de identificação.

CSK

Você também pode visualizar as propriedades de um vértice quando uma camada está no modo de edição selecionando vértices com a Ferramenta Nó. Em seguida, as coordenadas do vértice (x, y, r) são exibidas no painel Editor de vértices.

csk

É bom saber que é pelo menos possível ver o raio através da ferramenta de identificação. Com apenas alguns arcos, copiar e colar seria uma opção, mas pretendo usar isso para todos os dados que tenho, que possuem muitos recursos de arco para serem feitos manualmente. : /

Sisuaski

Pesquisei um pouco no Google e parece que muitas pessoas não usam recursos de cadeia circular no QGIS; portanto, não há uma ferramenta ou função pronta na calculadora de campo. Pode ser necessário definir uma função Python personalizada. Se você quiser seguir esse caminho, adicione a tag pygis à sua pergunta.

CSK

Quantos vértices você possui para o recurso? O que você vê Vertex Editorquando clica no Node Tool?

12483 Marco Marco

Respostas:

Seguem-se fórmulas levemente testadas, portanto, prossiga com cuidado. Mas seguindo um exemplo aqui: https://www.mathopenref.com/arcradius.html

Se seus arcos circulares tiverem um vértice no ponto do meio do arco (que eu estou dizendo é x1, y1 na figura), você poderá usá-lo juntamente com os pontos inicial e final para calcular o comprimento do acorde " W " e a altura " h " para obter o raio "R" usando a seguinte fórmula:

e dizendo

W =

sqrt(
     ( $x_at(-1) - $x_at(0) )^2 +
     ( $y_at(-1) - $y_at(0) )^2
    )

H =

sqrt(
    ( $x_at(1) - ( $x_at(-1) + $x_at(0) )/2 )^2
    + ( $y_at(1) - ( $y_at(-1) + $y_at(0) )/2 )^2
)

na sua caixa de diálogo Expressão do campo, você teria essa equação longa para o cálculo do raio:

R =

sqrt( 
    ( $x_at(1) - ($x_at(-1) + $x_at(0) )/2 )^2 
    + ( $y_at(1) - ($y_at(-1) + $y_at(0))/2 )^2
     ) /2
+
(   ( $x_at(-1) - $x_at(0) )^2 )  + 
    ( $y_at(-1) - $y_at(0) )^2 )  )
/ ( 8 * sqrt(
            ( ($x_at(1) - ($x_at(-1) + $x_at(0) )/2 )^2
            + ( $y_at(1) - ($y_at(-1) + $y_at(0) )/2 )^2
             )
   )

cm1
fonte

No LaTe X $$W= \sqrt {(x_{-1} - x_0)^2 +(y_{-1}-y_0)^2}$$ $$H = \sqrt { (x_1- \frac {(x_{-1}+x_0)} 2 )^2 + (y_1- \frac {(y_{-1}+y_0)} 2 )^2}$$ e

$$R = \frac {\sqrt { (x_1- \frac {(x_{-1}+x_0)}  2 )^2 + (y_1- \frac {(y_{-1}+y_0)} 2 )^2}} {2} + \frac {(x_{-1} - x_0)^2 +(y_{-1}-y_0)^2} {8 \sqrt { (x_1- \frac {(x_{-1}+x_0)}  2 )^2 + (y_1- \frac {(y_{-1}+y_0)} 2 )^2}}$$

Marco

Muito bom! Eu não imaginava como acessar os vértices. É possível no PostGIS?

17008 Marco Marco

Marco, acredito que é possível extrair qualquer vértice que você está procurando. Veja o link aqui: postgis.net/docs/ST_PointN.html . Obrigado pelas mais excelentes adições e edições de fórmula.

Cm1

Enquanto espero que os desenvolvedores adicionem um botão para "calcular facilmente o raio", esta é a resposta mais útil que eu poderia esperar. Agora eu (e possivelmente outros também) tenho pelo menos uma maneira de extrair raio. :) Ainda não testei isso, mas em alguns dias terei tempo suficiente. Muito obrigado pela sua ajuda!

Sisuaski

Que tal extrair o raio enquanto extrai os dados com o FME? Lá você tem o transformador ArcPropertyExtractor, que deve fornecer os valores de raio nos atributos prontos para serem configurados no seu PostGIS DB

fgiron
fonte

Essa seria uma maneira possível de proceder, uma boa maneira de lidar com todos os dados. Meu problema é que não há usuários suficientes de FME em minha organização e preciso tornar o processo o mais fácil possível. Isso significa que todo o processo, incluindo arcos, áreas, linhas etc., deve poder ser realizado com o Qgis. A capacidade de atualizar o raio do recurso Arc para novos objetos é o que eu preciso fazer com o Qgis. Obrigado pela sua resposta, entretanto, atualizarei os atributos de raio com o FME por enquanto, mas espero que haja uma maneira de fazer isso com o QGIS.

Sisuaski