QGIS oblicza promień łuku do atrybutów

Mam zainstalowany QGIS 2.18.16. Mój problem polega na tym, jak obliczyć promień łuku / cechy łuku (linii)? Linie pochodzą z pliku .dgn, który przeczytałem w FME i napisałem do PostGis DB. Utworzyłem klucze podstawowe itp., A tabela jest w pełni edytowalna. Pomyślnie obliczyłem długości łuków, ale nie mogę dowiedzieć się, jak obliczyć promień dla tych krzywych. To ( https://www.mathopenref.com/arcradius.html ) może dać jakiś pomysł po stronie matematyki, chociaż nie mogę wymyślić, jak go używać w kalkulatorze QGIS. Więc muszę obliczyć „R” dla wszystkich linii, które mam w DB.

Poniżej znajduje się przykład. Moje dane mają łuki / krzywe na innej tabeli niż „normalne” polilinie.

qgis field-calculator calculate-values Sisuaski
źródło

Kiedy klikam łuk za pomocą narzędzia Identyfikacja, jednym z wymienionych atrybutów pochodnych jest „najbliższy promień wierzchołka”. Dlatego nie ma potrzeby obliczania promienia, wystarczy dowiedzieć się, jak uzyskać dostęp do tej właściwości w kalkulatorze polowym. Lub po prostu skopiuj go z narzędzia do identyfikacji.

csk

Można również wyświetlić właściwości wierzchołka, gdy warstwa znajduje się w trybie edycji, wybierając wierzchołki za pomocą narzędzia Węzeł. Następnie współrzędne wierzchołków (x, y, r) są wyświetlane w panelu Edytora wierzchołków.

csk

Miło wiedzieć, że można przynajmniej zobaczyć promień za pomocą narzędzia do identyfikacji. Przy niewielu łukach kopiowanie i wklejanie byłoby opcją, ale planuję użyć tego do wszystkich danych, które mam, które mają zbyt wiele funkcji łuku, aby można je było wykonać ręcznie. : /

Sisuaski,

Przeszukiwałem trochę i wydaje się, że niewiele osób używa funkcji ciągów kołowych w QGIS, więc nie ma gotowego narzędzia ani funkcji w kalkulatorze polowym. Może być konieczne zdefiniowanie niestandardowej funkcji języka Python. Jeśli chcesz pójść tą drogą, dodaj znacznik Pygis do swojego pytania.

csk

Ile masz wierzchołków dla elementu? Co widzisz Vertex Editorpo kliknięciu Node Tool?

Marco

Odpowiedzi:

Następują lekko przetestowane formuły, więc zachowaj ostrożność. Ale podążając wraz z przykładem tutaj: https://www.mathopenref.com/arcradius.html

Jeśli twoje okrągłe łuki mają wierzchołek w środkowym punkcie wzdłuż łuku (co mówię to x1, y1 na rysunku), możesz użyć go wraz z punktami początkowym i końcowym do obliczenia długości cięciwy „ W ” i wysokości „ h ”, aby uzyskać promień „R” przy użyciu następującego wzoru:

i mówiąc

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
)

w oknie dialogowym Wyrażenie kalkulatora pola będziesz mieć takie długie równanie do obliczenia promienia:

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
źródło

W 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}$$ i

$$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

Bardzo dobre! Nie wyobrażałem sobie, jak uzyskać dostęp do wierzchołków. Czy jest to możliwe w PostGIS?

Marco

Marco, uważam, że możliwe jest wydobycie dowolnego wierzchołka, którego szukasz. Zobacz link tutaj: postgis.net/docs/ST_PointN.html . Dziękujemy za najlepsze uzupełnienia i zmiany formuł.

cm1

Jest to najbardziej pomocna odpowiedź, na którą czekam, aż programiści dodadzą przycisk „łatwego obliczenia promienia”. Teraz ja (i prawdopodobnie także inni) mam przynajmniej jakiś sposób na wyodrębnienie promienia. :) Jeszcze tego nie testowałem, ale zrobię to za kilka dni, kiedy będę miał wystarczająco dużo czasu. Bardzo dziękuję za pomoc!

Sisuaski

Co z wyodrębnieniem promienia podczas wyodrębniania danych za pomocą FME? Tam masz transformator ArcPropertyExtractor, który powinien dać ci wartości promienia w atrybutach gotowych do ustawienia w DB PostGIS

firon
źródło

Byłby to jeden z możliwych sposobów postępowania, dobry sposób na obsługę wszystkich danych. Mój problem polega na tym, że w mojej organizacji nie ma wystarczającej liczby użytkowników FME i muszę maksymalnie ułatwić ten proces. Oznacza to, że cały proces, w tym łuki, obszary, linie itp., Powinien być możliwy do wykonania za pomocą Qgis. Muszę zaktualizować promień łuku do nowych obiektów, co muszę zrobić z Qgis. Dziękuję za odpowiedź, na razie będę aktualizować atrybuty promienia za pomocą FME, ale mam nadzieję, że będzie sposób, aby to zrobić za pomocą QGIS.

Sisuaski