পোস্টজিআইএস - কীভাবে দক্ষতার সাথে একটি একক টেবিলের সমস্ত ওভারল্যাপিং বহুভুজগুলি ST_Union করতে হবে


13

আমার লক্ষ্যটি হ'ল একক টেবিলটি নিয়ে এবং st_union সমস্ত বহুভুজগুলি যা স্পর্শ করে বা একে অপরের সাথে একক বহুভুতে পরিণত হয়

আমি একজন সি # বিকাশকারী যিনি পোস্টজিআইএস সম্পর্কে শিখতে শুরু করেছেন। নীচের কোডটি ব্যবহার করে, আমি এটি সম্পাদন করতে সক্ষম হয়েছি, তবে এটি অদৃশ্য বলে মনে হচ্ছে এবং পোস্টজিআইএস-এর আরও অনেক কিছুই আমার কাছে নতুন।

আমার প্রাথমিক প্রচেষ্টা থেকে (এখনও মন্তব্যে রয়েছে), আমি একবারে কেবল পলিকে একত্রিত না করে ST_UNION এর সাথে অ্যারে_্যাগ ব্যবহার করে পুনরাবৃত্তিগুলি হ্রাস করতে সক্ষম হয়েছি।

আমি আমার উত্স 173 থেকে 133 পোল নিয়ে বাতাসে চলেছি।

sql = "DROP TABLE IF Exists tmpTable; create table tmpTable ( ID varchar(50), Geom  geometry(Geometry,4326), Touchin varchar(50) ); create index idx_tmp on tmpTable using GIST(Geom); ";
                CommandText = sql;
                ExecuteNonQuery();

                sql = "";
                for (int i = 0; i < infos.Count(); i++)
                {
                    sql += "INSERT INTO tmpTable SELECT '" + infos[i].ID + "', ST_GeomFromText('" + infos[i].wkt + "', 4326), '0';";
                }
                CommandText = sql;
                ExecuteNonQuery();

                CommandText = "update tmpTable set touchin = (select id from tmpTable as t where st_intersects(st_buffer(geom, 0.0001), (select geom from tmpTable as t2 where t2.ID = tmpTable.ID ) ) and t.ID <> tmpTable.ID limit 1)";
                ExecuteNonQuery();

                CommandText = "select count(*) from tmpTable where touchin is not null";
                long touching = (long)ExecuteScalar();
                string thisId = "";
                // string otherId = "";
                while (touching > 0)
                {
                    CommandText = "select touchin, count(*)  from tmpTable where touchin is not null group by touchin order by 2 desc limit 1";
                    //CommandText = "select id, touchin from tmpTable where touchin is not null";
                    using (var prdr = ExecuteReader())
                    {
                        CommandText = "";
                        if (prdr.Read())
                        {
                            thisId = prdr.GetString(0);
                             // otherID = prdr.GetString(1);
                            CommandText = @"update tmpTable set geom = st_union(unioned) 
                                from (select array_agg(geom) as unioned from tmpTable where touchin = '" + thisId + "' or id = '" + thisId + @"') as data
                                where id = '" + thisId + "'";
                             // CommandText = "update tmpTable set geom = st_union(geom, (select geom from tmpTable where ID = '" + otherId + "')) where id = '" + thisId + "'";
                        }
                    }

                    if (!string.IsNullOrEmpty(CommandText))
                    {
                        ExecuteNonQuery();
                        //CommandText = "update tmpTable set geom = null, touchin = null where ID = '" + otherId + "'";
                        CommandText = "update tmpTable set geom = null, touchin = null where touchin = '" + thisId + "'";
                        ExecuteNonQuery();                            
                    }

                    CommandText = "update tmpTable set touchin = (select id from tmpTable as t where st_intersects(st_buffer(geom, 0.0001), (select geom from tmpTable as t2 where t2.ID = tmpTable.ID ) ) and t.ID <> tmpTable.ID limit 1)";
                    ExecuteNonQuery();

                    CommandText = "select count(*) from tmpTable where touchin is not null";
                    touching = (long)ExecuteScalar();
                }

আমি যে ডেটাসেটটি ব্যবহার করছি তার নমুনা এখানে:

INSERT INTO tmpTable SELECT '872538', ST_GeomFromText('POLYGON((-101.455035985 26.8835084441,-101.455035985 26.8924915559,-101.444964015 26.8924915559,-101.444964015 26.8835084441,-101.455035985 26.8835084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872550', ST_GeomFromText('POLYGON((-93.9484752173 46.0755084441,-93.9484752173 46.0844915559,-93.9355247827 46.0844915559,-93.9355247827 46.0755084441,-93.9484752173 46.0755084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872552', ST_GeomFromText('POLYGON((-116.060688575 47.8105084441,-116.060688575 47.8194915559,-116.047311425 47.8194915559,-116.047311425 47.8105084441,-116.060688575 47.8105084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872553', ST_GeomFromText('POLYGON((-116.043688832 47.8125084441,-116.043688832 47.8214915559,-116.030311168 47.8214915559,-116.030311168 47.8125084441,-116.043688832 47.8125084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872557', ST_GeomFromText('POLYGON((-80.6380222359 26.5725084441,-80.6380222359 26.5814915559,-80.6279777641 26.5814915559,-80.6279777641 26.5725084441,-80.6380222359 26.5725084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872558', ST_GeomFromText('POLYGON((-80.6520223675 26.5755084441,-80.6520223675 26.5844915559,-80.6419776325 26.5844915559,-80.6419776325 26.5755084441,-80.6520223675 26.5755084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872559', ST_GeomFromText('POLYGON((-80.6400224991 26.5785084441,-80.6400224991 26.5874915559,-80.6299775009 26.5874915559,-80.6299775009 26.5785084441,-80.6400224991 26.5785084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872560', ST_GeomFromText('POLYGON((-80.6530226307 26.5815084441,-80.6530226307 26.5904915559,-80.6429773693 26.5904915559,-80.6429773693 26.5815084441,-80.6530226307 26.5815084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872568', ST_GeomFromText('POLYGON((-90.7892258584 30.7365084441,-90.7892258584 30.7454915559,-90.7787741416 30.7454915559,-90.7787741416 30.7365084441,-90.7892258584 30.7365084441))', 4326), '0';
INSERT INTO tmpTable SELECT '872569', ST_GeomFromText('POLYGON((-90.7832259127 30.7375084441,-90.7832259127 30.7464915559,-90.7727740873 30.7464915559,-90.7727740873 30.7375084441,-90.7832259127 30.7375084441))', 4326), '0';

আপনার প্রশ্নে আসল ডেটা নিজেই প্রয়োজন?
পল

@ পল - নিশ্চিত ছিলেন না যে এটি সহায়ক হবে কিনা।
ক্যারল অ্যান্ডারমার্টেন লাইবস্টার 18

উত্তর:


20

সহজ উপায়টি ST_Unionপুরো টেবিলের দিকে যাবে:

SELECT ST_Union(geom) FROM tmpTable;

এটি আপনাকে একটি বিশাল দেবে MultiPolygon, যা সম্ভবত আপনি চান তা নয়। আপনি স্বতন্ত্র দ্রবীভূত উপাদানগুলির সাথে পেতে পারেন ST_Dump। তাই:

SELECT (ST_Dump(geom)).geom FROM (SELECT ST_Union(geom) AS geom FROM tmpTable) sq;

এটি আপনাকে স্পর্শকারী ইনপুটগুলির প্রতিটি সেটের জন্য পৃথক বহুভুজ দেয় , তবে সামান্য দূরত্বের দ্বারা পৃথক করা ইনপুটগুলির গোষ্ঠীগুলি পৃথক জ্যামিতি হিসাবে থাকবে। আপনার অ্যাক্সেস রয়েছে এমন PostGIS 2.2.0rc1 , আপনি জ্যামিতি যে ঘনিষ্ঠ একটি একক মধ্যে একসাথে মার্জ করতে পারবেন GeometryCollectionব্যবহার ST_ClusterWithin :

SELECT unnest(ST_ClusterWithin(geom, 0.0001)) AS grp FROM tmpTable;

আপনি যদি চান Polygonsমধ্যে GeometryCollectionধ্বংস হতে যাচ্ছে তখন আপনি চালাতে পারেন ST_UnaryUnionপ্রতিটি GeometryCollection, ফলে মত:

SELECT ST_UnaryUnion(grp) FROM
(SELECT unnest(ST_ClusterWithin(geom, 0.0001)) AS grp FROM tmpTable) sq;

এটি অবশ্যই খুব দ্রুত, তবে দুটি জিনিস: (1) ফলাফলগুলিতে আমি কী আইডি ক্ষেত্রটি সংরক্ষণ করতে পারি? এটি কোনটি গুরুত্বপূর্ণ তা নয়, তবে আমার ফলাফল নেওয়া এবং এটি থেকে অন্যান্য ডেটা নেওয়া দরকার। (২) আবার এস ইনফার বাফার যুক্ত করার কোনও উপায় আছে?
ক্যারল অ্যান্ডারমার্টেন লেবস্টার

1
(1) সহজে নয়, তবে বৈশিষ্ট্যগুলি ফিরে পাওয়ার একটি সহজ উপায় স্থানগতভাবে আপনার ফলাফল বহুভুজগুলিকে আপনার ইনপুট বহুভুজগুলির অভ্যন্তর বিন্দুর বিপরীতে যোগ দেওয়া। (২) জ্যামিতিগুলি হ্যান্ডল করার জন্য কিছু ব্যাখ্যা যুক্ত করা হয়েছে যা কাছাকাছি হলেও স্পর্শ করে না।
dbaston

সহায়তার জন্য ধন্যবাদ - আমার কাছে বর্তমানে ২.২ নেই, সুতরাং আমি যখন এটি আপগ্রেড করব তখন আমাকে এটি পুনরায় দেখাতে হবে। আপাতত বাফারকে বাদ দেওয়া কোনও ডিল ব্রেকার নয়।
ক্যারল অ্যান্ডারমার্টেন লাইবস্টার 21

আমি সম্মতি জানাই এটি সহজতম। আমি অবাক হচ্ছি যে যদি কোনও পুনরাবৃত্ত কোয়েরি করার কোনও উপায় থাকে যা
স্পর্শকৃত রত্নগুলি

1
@ ক্রিসমার্কস, পুনরাবৃত্তির সমাধান সম্পর্কে কিছু ধারণার জন্য gis.stackexchange.com/a/94228/18189 দেখুন
dbaston
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.