QGIS एक चाप की त्रिज्या को विशेषताओं की गणना करता है

Question

मेरे पास QGIS 2.18.16 है। मेरी समस्या यह है कि, आर्क / वक्र सुविधा (लाइन) की त्रिज्या की गणना कैसे करें? लाइनें .dgn फ़ाइल से उत्पन्न होती हैं, जिसे मैंने FME के ​​साथ पढ़ा और PostGis DB को लिखा। मैंने प्राथमिक कुंजी आदि बनाई है, और तालिका पूरी तरह से संपादन योग्य है। मैंने सफलतापूर्वक आर्क्स के लिए लंबाई की गणना की, लेकिन इन वक्रों के लिए त्रिज्या की गणना करने का तरीका पता नहीं लगा सका। यह ( https://www.mathopenref.com/arcradius.html ) गणित पक्ष पर कुछ विचार दे सकता है, हालांकि मैं यह पता नहीं लगा सकता कि इसे क्यूजीआईएस कैलकुलेटर में कैसे उपयोग किया जाए। इसलिए मुझे डीबी में मौजूद सभी लाइनों के लिए "आर" की गणना करने की आवश्यकता है।

Bellow एक उदाहरण है। मेरे डेटा में "सामान्य" पॉलीइन्स की तुलना में अलग-अलग टेबल पर आर्क / कर्व्स हैं।

Comments

जब मैं आइडेंटिफाई टूल के साथ एक आर्क पर क्लिक करता हूं, तो सूचीबद्ध विशेषताओं में से एक "निकटतम शीर्ष त्रिज्या" है। इसलिए त्रिज्या की गणना करने की कोई आवश्यकता नहीं है, आपको केवल यह पता लगाने की आवश्यकता है कि फ़ील्ड कैलकुलेटर में इस संपत्ति का उपयोग कैसे किया जाए। या बस इसे पहचान उपकरण से कॉपी करें।

— csk

जब आप नोड मोड के साथ कोने का चयन करके एक परत संपादित मोड में होते हैं, तो आप एक शीर्ष के गुणों को भी देख सकते हैं। तब वर्टेक्स संपादक पैनल में वर्टेक्स निर्देशांक (x, y, r) प्रदर्शित करता है।

— csk

यह जानकर अच्छा लगा कि पहचान उपकरण के माध्यम से त्रिज्या को देखना कम से कम संभव है। केवल कुछ आर्क कॉपी-पेस्ट के साथ एक विकल्प होगा, लेकिन मैं इसे पूरे डेटा के लिए उपयोग करने की योजना बना रहा हूं, जिसमें मैन्युअल रूप से किए जाने के लिए बहुत अधिक आर्क विशेषताएं हैं। : /

— सिसुस्की

मैं थोड़ा सा घूम गया, और ऐसा लगता है कि बहुत से लोग क्यूजीआईएस में परिपत्र स्ट्रिंग विशेषताओं का उपयोग नहीं करते हैं, इसलिए क्षेत्र कैलकुलेटर में तैयार उपकरण या फ़ंक्शन नहीं है। आपको एक कस्टम पायथन फ़ंक्शन को परिभाषित करने की आवश्यकता हो सकती है। यदि आप उस मार्ग पर जाना चाहते हैं, तो अपने प्रश्न में pygis टैग जोड़ें।

— csk

आपके पास सुविधा के लिए कितने वर्टेक्स हैं? Vertex Editorजब आप क्लिक करते हैं तो आप क्या देखते हैं Node Tool?

— मार्को

Answer 1

हल्के से परीक्षण किए गए सूत्र अनुसरण करते हैं, इसलिए सावधानी से आगे बढ़ें। लेकिन यहाँ एक उदाहरण के साथ निम्नलिखित: https://www.mathopenref.com/arcradius.html

यदि आपके परिपत्र चाप में चाप के साथ मध्य बिंदु पर एक शीर्ष है (जो मैं कह रहा हूं कि चित्र में X1, y1 है), तो आप इसे कॉर्ड लंबाई " डब्ल्यू " और ऊंचाई कीगणना करने के लिए प्रारंभ और अंत बिंदुओं के साथ उपयोग कर सकते हैं।" एच " निम्नलिखित सूत्र का उपयोग करके त्रिज्या "आर" प्राप्त करने के लिए:

और कह रहा है

W =

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

तथा

ज =

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

फ़ील्ड कैलकुलेटर के अपने अभिव्यक्ति डायलॉग में आपको त्रिज्या गणना के लिए यह लंबा समीकरण होगा:

आर =

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

Comments

लाटे एक्स $$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}$$ और $$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}}$$

— मार्को

बहुत अच्छा! मैंने कल्पना नहीं की थी कि शीर्षकों तक कैसे पहुंचा जाए। क्या यह PostGIS में संभव है?

— मार्को

मार्को, मेरा मानना ​​है कि जो भी वर्टिकल आप ढूंढ रहे हैं उसे निकालना संभव है। : लिंक यहाँ देखें postgis.net/docs/ST_PointN.html । सबसे उत्कृष्ट सूत्र परिवर्धन और संपादन के लिए धन्यवाद।

— सेमी 1

डेवलपर्स के लिए "आसान गणना त्रिज्या" के लिए एक बटन जोड़ने की प्रतीक्षा करते समय, यह सबसे उपयोगी उत्तर है जिसके लिए मैं आशा कर सकता हूं। अब मैं (और संभवतः अन्य भी) त्रिज्या निकालने के लिए कम से कम कुछ तरीका है। :) मैंने अभी तक इसका परीक्षण नहीं किया है, लेकिन मैं कुछ दिनों में जब मेरे पास पर्याप्त समय होगा। आपकी मदद के लिए बहुत बहुत धन्यवाद!

— सिसुस्की

Answer 2

FME के ​​साथ डेटा निकालने के दौरान त्रिज्या निकालने के बारे में क्या? वहां आपके पास ArcPropertyExtractor ट्रांसफॉर्मर है, जो आपको अपने GGB DB में सेट करने के लिए तैयार विशेषताओं में त्रिज्या मान देना चाहिए

Comments

यह आगे बढ़ने का एक संभव तरीका होगा, पूरे डेटा को संभालने का एक अच्छा तरीका। मेरी समस्या यह है कि मेरे संगठन में पर्याप्त FME उपयोगकर्ता नहीं हैं और मुझे इस प्रक्रिया को यथासंभव आसान बनाने की आवश्यकता है। इसका मतलब है कि पूरी प्रक्रिया, जिसमें आर्क, क्षेत्र, रेखाएं आदि शामिल हैं, को Qgis के साथ करने में सक्षम होना चाहिए। नई वस्तुओं के लिए आर्क सुविधा त्रिज्या को अद्यतन करने की क्षमता है जो मुझे Qgis के साथ करने की आवश्यकता है। हालांकि आपके उत्तर के लिए धन्यवाद, मैं अभी के लिए FME के ​​साथ त्रिज्या विशेषताओं को अपडेट कर रहा हूं, लेकिन आशा है कि इसे GGIS के साथ करने का एक तरीका होगा।

— सिसुस्की