PostGIS में ज्यामितीय सफाई?

मैं कुछ बहुत बड़ी बहुभुज परतों पर कुछ प्रसंस्करण करने की कोशिश कर रहा हूं। हालाँकि मैं विभिन्न ज्यामिति त्रुटियों में भाग रहा हूँ जैसे:

NOTICE:  Ring Self-intersection at or near point 470396.52017068537 141300.52235257279
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 504154.61769969884 140782.04115761846
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 505255.50242871145 140803.34860398644
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 510312.46970004693 141215.29256710084
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 510312.46970004693 141215.29256710084
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 511839.50335641927 141115.85781738357
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 515064.03024010791 140895.68087158105
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 519233.18724611058 140881.47590733573
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 521072.73011588014 141044.83299615697
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 523331.31943088671 141144.26774587421
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 523331.31943088671 141144.26774587424
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 523395.24176999065 140725.22130063715
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1
NOTICE:  Ring Self-intersection at or near point 524531.63890961662 140810.45108610913
CONTEXT:  PL/pgSQL function st_intersection(geometry,raster,integer) line 10 at RETURN QUERY
SQL function "st_intersection" statement 1

मैंने यहां सुझाए गए फ़ंक्शन की कोशिश की है: https://trac.osgeo.org/postgis/wiki/UsersWikiCleanPolygons

जियोमेट्री की सफाई के लिए, मेरे द्वारा उपयोग किए जा रहे कोड:

UPDATE public.mytable
SET geom=cleangeometry(geom);

नतीजे के साथ:

ERROR:  GEOSisSimple: IllegalArgumentException: This method does not support GeometryCollection arguments

और भी

UPDATE public.valid_mytable
SET geom=ST_MakeValid(geom);

यह एक काम करता है, लेकिन केवल अगर मैं पहली बार अपने ज्यामिति कॉलम को ज्यामिति में बदलता हूं

ALTER TABLE public.mytable  ALTER COLUMN geom SET DATA TYPE geometry;

जो तब मुझे एक मेज के साथ छोड़ देता है जो अब मेरे अन्य कार्यों के साथ काम नहीं करती है!

ERROR:  Relate Operation called with a LWGEOMCOLLECTION type.  This is unsupported.

मैंने कॉलम को वापस ज्यामिति में बदलने की कोशिश की है (मल्टीपोलियन)

ALTER TABLE public.my_table ALTER COLUMN geom सेट डेटा प्रकार ज्यामिति (मल्टीपोलियन);

लेकिन यह विफल हो जाता है

ERROR:  Geometry type (GeometryCollection) does not match column type (MultiPolygon)

मैंने PostGIS इन एक्शन (सेकेंड एड) http://www.manning.com/obe/ के माध्यम से जाने की कोशिश की है, लेकिन मैं केवल अवैध ज्यामिति खोजने के लिए फ़ंक्शंस पा सकता हूं, लेकिन मेरा डेटा-सेट मैन्युअल रूप से इसे ठीक करने के लिए इतना बड़ा है, मैं वास्तव में कुछ ऐसा चाहिए जो उन्हें अपने आप ठीक कर दे।

जब मैं कोशिश करता हूं और ST_MakeValid () मैं परिणाम प्राप्त करता हूं, तो मैं समस्या बहुभुज को अलग करने में सक्षम हूं:

ERROR:  Geometry type (GeometryCollection) does not match column type      (MultiPolygon)
 ********** Error **********

 ERROR: Geometry type (GeometryCollection) does not match column type      (MultiPolygon)
SQL state: 22023

मैंने अपने ज्यामिति स्तंभ पर एक प्रकार की जाँच की, और यह कहा कि प्रकार "MULTIPOLYGON" था

postgis postgresql postgis-2.0

— बाजार
स्रोत

ST_MakeValid जितना हो सके उतना सही करे।

— user30184

मैं देखता हूं, धन्यवाद, मैंने वास्तव में मेरे प्रश्न में एक गलती की थी जहां मैं यह उल्लेख करना भूल गया था कि यह ST_Make_Valid था जो मेरे कॉलम के साथ समस्याओं का कारण बनता है। मैंने ST_MakeValid का उपयोग किया है, लेकिन मुझे काम करने के लिए अपने जियोम कॉलम को ज्यामिति डेटा प्रकार में बदलना होगा, और एक बार जब मैं ऐसा करूं कि मैं इसे ज्यामिति (मल्टीपॉलगन) में वापस ला सकूं

— Mart

आप हैक ST_Buffer (जियोम, 0) का उपयोग कर सकते हैं जो कई अवैध ज्यामिति से निपटेगा। आप ST_MakeValid का भी उपयोग कर सकते हैं। अंत में, आप एक नई तालिका में चयन करने और जहाँ खंड में ST_IsValid (जियोम) डालने की कोशिश कर सकते हैं।

— जॉन पॉवेल

धन्यवाद, मैंने पहले ही बफर हैक की कोशिश की, लेकिन यह काम नहीं किया, यह एक ज्यामिति (मल्टीपोलियन) के बजाय एक ज्यामिति इनपुट चाहता था। मैं केवल मान्य बहुभुज का चयन करने का प्रयास करूंगा और देखूंगा कि कितने फ़िल्टर किए गए हैं।

— Mart

ठीक है। यह पॉलीगॉन के साथ st_makevalid उत्पादक बिंदुओं और LineStrings से आता है जो एक GeometryCollection का उत्पादन करेगा। इसके लिए एक फिक्स है जो मैं कुछ घंटों में लिखूंगा। मैं सर्फिंग करने जा रहा हूं :-)

— जॉन पॉवेल

आप केवल से बहुभुज या Multipolygons चाहते हैं ST_MakeValid आप उपयोग कर सकते ST_Dump ज्यामिति प्रकार के लिए घटक ज्यामिति और उसके बाद परीक्षण को निकालने के लिए। ST_MakeValid कभी-कभी पॉइंट्स या लाइनस्ट्रीम का उत्पादन करेगा, जहां से ज्योमेट्रीकोलेक्शन होता है। कुछ इस तरह का प्रयास करें:

SELECT 
  g.geom, 
  row_number() over() AS gid,
FROM 
  (SELECT 
     (ST_DUMP(ST_MakeValid (geom))).geom FROM your_table
  ) AS g
WHERE ST_GeometryType(g.geom) = 'ST_MultiPolygon' 
   OR ST_GeometryType(g.geom) = 'ST_Polygon';

आप OR शर्त के बजाय IN खंड का उपयोग कर सकते हैं, हालांकि परिणाम और क्वेरी योजना समान होगी। यदि आप केवल मल्टीपोलीगॉन चाहते हैं तो आप ST_Dump को ST_Multi फ़ंक्शन में लपेट सकते हैं ।

Row_number () ओवर () बस आपको एक अनूठी आईडी लौटाएगा, जो एक से शुरू होकर, ST_Dump से लौटी प्रत्येक ज्यामिति के लिए होगी। आप उसी परिणाम के साथ ST_Dump द्वारा दिए गए पथ तत्व का उपयोग भी कर सकते हैं।

आप निश्चित रूप से इसे एक CREATE TABLE cleaned_geoms AS SELECT .... टाइप स्टेटमेंट के साथ जोड़ना चाहते हैं, क्योंकि एक सीधा अपडेट ST_MakeValid wil Notl के रूप में काम करने की संभावना नहीं है (या हमेशा) I से आउटपुट में एक से एक टैपिंग का उत्पादन होता है।

यह अनुपलब्ध है क्योंकि मेरे पास वर्तमान में कोई साधन नहीं है, इसलिए गलत तरीके से किया गया कोष्ठक हो सकता है, लेकिन सामान्य सिद्धांत ध्वनि है। उम्मीद है की यह मदद करेगा।

— जॉन पॉवेल
स्रोत

आप ST_CollectionExtract को मल्टीमीटर से [मल्टी] पॉलीगॉन निकालने के लिए आज़मा सकते हैं । उन्हें MuliPolygons के रूप में मजबूर करने के लिए ST_Multi का उपयोग करें।

UPDATE public.valid_lcmsouthshapefile
  SET geom=ST_Multi(ST_CollectionExtract(ST_MakeValid(geom), 3))
  WHERE NOT ST_IsValid(geom);

अपने किए जाने के बाद, यह सुनिश्चित करने के लिए कि वे वैध रहें, एक CHECK बाधा का उपयोग करें। विवरण देखें यहाँ ।

— माइक टी
स्रोत