मैंने एक मानचित्र को दोबारा तालिका में कैसे खींचा

फिर से - क्योंकि किसी भी तरह से 2 साल पहले मैं पहले से ही यह अभ्यास कर रहा था। यह एक लंबी और कठिन प्रक्रिया थी:
blogs.technet.com/b/isv_team/archive/2010/01/18/3306462.aspx
blogs.technet.com/b/isv_team/archive/2010/01/21/3307201.aspx
blogs.technet.com/b/isv_team/archive/2010/01/23/3307719.aspx
blogs.technet.com/b/isv_team/archive/2010/01/24/3307804.aspx
तब से, विज्ञान ने बहुत आगे बढ़ा दिया है। इस पोस्ट में, हम फिर से SQL Server को हमारे विशाल देश के नक्शे पर अपलोड करते हैं, इस बार मानचित्र लेखकों, नए SQL सर्वर सुविधाओं और स्वतंत्र डेवलपर्स के लिए बहुत आसान और अधिक धन्यवाद, जो आप सभी को बहुत बहुत धन्यवाद देते हैं। हमें आवश्यकता होगी:

1. उपयुक्त कार्ड
2. SQL सर्वर 2012 (एक्सप्रेस हो सकता है)
3. एक नक्शा को अपने मूल प्रारूप से SQL सर्वर में बदलने का उपकरण।

अधिकांश कार्य डेनिश MVP मोर्टन gadegaard Nielsen द्वारा लिखित अद्भुत मुफ्त Shape2SQL उपयोगिता द्वारा लिया जाएगा, ईएसआरआई में जीआईएस लीड सॉफ्टवेयर इंजीनियर के साथ काम करते हुए, भगवान उसे स्वास्थ्य प्रदान कर सकते हैं। ESRI (एनवायरनमेंटल सिस्टम रिसर्च इंस्टीट्यूट) को सॉफ्टवेयर उत्पादों के अपने आर्कजीस परिवार के लिए जाना जाता है जो रूस में भी व्यापक रूप से उपयोग किए जाते हैं। वेक्टर प्रारूप फ़ाइल प्रारूप को 90 के दशक की शुरुआत में ArcView GIS संस्करण 2 के लिए पेश किया गया था , और आज, इसकी लोकप्रियता के कारण, यह भौगोलिक सूचना प्रणालियों के बीच डेटा विनिमय के लिए वास्तविक मानक बन गया है। जैसा कि नाम से ही स्पष्ट है कि तुला एसक्यूएल सर्वर 2008 और उससे अधिक आकार की फाइलों को लोड कर सकता है। 2012 का उपयोग इस पोस्ट में किया जाता है क्योंकि इसमें दिखाई देने वाले ज्यामितीय / भौगोलिक समुच्चय का समर्थन - यह बाद में काम आएगा। 2008 में, मुझे एक कर्सर लिखना होगा।
नोट 1. SQL सर्वर में .shp आयात उपयोगिता के अलावा, SharpGIS वेबसाइट पर SqlSpatialTools v1.3.0 (348 kb) (बिल्ड 3759) में WPF और सिल्वरलाइट के लिए SqlSpatial क्वेरी टूल मैप क्लाइंट शामिल है, जो भौगोलिक पैनल की तुलना में क्षमताओं में बहुत समृद्ध है। SSMS परिणाम
नोट 2. लेखक द्वारा चेतावनी के रूप में दोनों उपकरण, व्यावसायिक उपयोग के लिए अभिप्रेत नहीं हैं, जैसा कि उसकी अपनी जरूरतों या आनंद के लिए उसके द्वारा लिखित।
नोट 3. एक खुला स्रोत .NET लाइब्रेरी जो आकार फ़ाइलों को पार्स कर सकती है, यहाँ स्थित है ।

यह रूसी संघ के प्रशासनिक-क्षेत्रीय विभाजन के साथ आकार-फ़ाइल को खोजने के लिए बना हुआ है। जीआईएस-लैब वेबसाइट पर .shp सहित विभिन्न स्वरूपों में रूसी क्षेत्रों के मानचित्र पाए गए: इस परत को CC-BY के तहत लाइसेंस प्राप्त है, जो किसी भी उपयोग की अनुमति देता है, जिसमें स्थानांतरण, संशोधन, वाणिज्यिक परियोजनाओं में उपयोग, प्राधिकरण और उपलब्धता के अधीन है। हाइपरलिंक। नक्शे के लिए डेटा स्रोत संघीय राज्य पंजीकरण, कैडस्ट्रे और कार्टोग्राफी सेवा ( रोसरेस्टर ) - संस्थाओं की सीमाएँ, ओपनस्ट्रेपोर रूस ( OSM ) - रूसी भूमि सीमा, कलिनिनग्राद क्षेत्र + मास्को की सीमा, VMap0 को छोड़कर - कलिनिनग्राद क्षेत्र की तटवर्ती + सीमाएँ और पहल समूह हैं। GIS-Lab.info परियोजना, जिसके लिए वह गहरी प्रशंसा भी व्यक्त करता है।
जीआईएस-लैब कार्ड 3 अनुमानों में प्रदान किए जाते हैं: डब्ल्यूजीएस 1984 (जीपीएस), पुलकोवो -1942 और अलबर्स-साइबेरिया। रूस के क्षेत्र को प्रदर्शित करने के लिए, तीसरा सबसे उपयुक्त है:


चित्र 1

पहले दो की तुलना में:


Fig.2

दुर्भाग्य से, मैं केंद्रीय मध्याह्न = 1050 पूर्व के साथ अल्बर्ट्स के बराबर शंक्वाकार प्रक्षेपण sys.spatial_reference_systems में नहीं मिला। Krasovsky गोलाकार पर आधारित है, इसलिए इसे एसक्यूएल सर्वर में ज्यामिति के रूप में खींचा जा सकता है और उदाहरण के लिए, रिपोर्टों को देखने के साधन के रूप में उपयोग किया जाता है। पहले दो को भूगोल प्रकार (SRID = 4326, 4284) में आयात किया जाता है, और उनका उपयोग कार्टिकल गणना के लिए किया जा सकता है।

नक्शा आयात करने के लिए, हमें तीन फ़ाइलों की आवश्यकता है:


3 चित्र

.Shp / .shx फ़ाइलों को ऊपर उल्लिखित Shape2SQL रनटाइम टूल को खिलाया जाता है। आकृति फ़ाइल का पूरा नाम, SQL सर्वर से कनेक्शन, तालिका का नाम जिसमें आकृति फ़ाइल से डेटा आयात किया जाएगा, तालिका में भू-स्थानिक स्तंभ का नाम, स्तंभ का नाम जो प्राथमिक कुंजी के रूप में कार्य करेगा, इंगित किया गया है। यदि तालिका मौजूद नहीं थी, तो यह बनाई गई है। भूगोल के रूप में Albers साइबेरिया अनुमानों में डेटा आयात करना विफल रहता है SQL सर्वर ऐसे SRID को नहीं जानता है और गलत तरीके से निर्देशांक का अनुवाद करता है - डेटा प्रक्षेपण या सीमा SqlGeography प्रकार द्वारा समर्थित की सीमा के बाहर हैं - इसलिए हम उन्हें ज्यामिति के रूप में आयात करते हैं। यदि वांछित है, तो आप तुरंत एक भू-स्थानिक सूचकांक का निर्माण कर सकते हैं। डेटाबेस बटन पर अपलोड पर क्लिक करें, और dbo.regions2010_albsib तालिका सफलतापूर्वक बनाई गई है और डेटा के साथ आबादी है:


चित्र 4

आप select * from rfmap.dbo.regions2010_albsib कर सकते हैं और यह सुनिश्चित कर सकते हैं कि SSMS में स्थानिक परिणाम टैब पर सब कुछ ठीक है (चित्र 1 देखें)। हालाँकि, तालिका में प्रविष्टियों की संख्या रूसी क्षेत्रों की संख्या से अधिक है। यह समझने के लिए कि मामले में क्या है, हमें तीसरी फ़ाइल की आवश्यकता है - dbf, जिसमें मानचित्र के लिए मेटाडेटा है। SQL सर्वर में dbf कैसे लोड करें, हम पिछली पोस्ट में गए थे । एक तालिका बनाएं

CREATE TABLE dbo.regions2010_captions (
ID int IDENTITY(1,1) NOT NULL,
AREA float NULL,
PERIMETER float NULL,
region varchar(50) NULL
)

और आयात / निर्यात विज़ार्ड का उपयोग करें। सिद्धांत रूप में, यह स्वचालित रूप से स्रोत और गंतव्य फ़ील्ड को समान नाम से जोड़ता है, लेकिन एक बार फिर आप इसे क्लिक करके सत्यापित कर सकते हैं


चित्रा 5

यदि क्षेत्र फ़ील्ड की लंबाई पहले से ज्ञात नहीं है, अर्थात। हाथ में कोई डीबीएफ रीडर नहीं है, और हेडर को बाइनरी संपादक के साथ आलस्य के साथ पार्स करते हुए, आप अधिकतम लंबाई निर्धारित कर सकते हैं, और आयात के बाद इसे काट सकते हैं। Dbf का आयात सफल है, और SQL सर्वर डेटाबेस में कार्टोग्राफिक ऑब्जेक्ट्स की तालिका के साथ, उनके लिए हस्ताक्षर की एक तालिका दिखाई देती है:


चित्रा 6

मानचित्र वस्तुओं के क्षेत्रों के अभिलेखों में2010_albsib हस्ताक्षर तालिका के साथ जुड़े हुए हैं ।2010_captions विशेष रूप से भौतिक क्रम से हैं, जो एक चर्चा नहीं है, क्योंकि किसी ने भी वादा नहीं किया कि आयात / निर्यात विज़ार्ड उसी क्रम में अभिलेखों को स्थानांतरित करता है जैसे वे dbf में थे। इस स्थिति से बाहर का सबसे अच्छा तरीका होगा dbf और shp में PK की उपस्थिति, जिससे उनकी तुलना करना आसान हो जाएगा और, वैसे, Shape2SQL के लेखक को सरोगेट आईडी के साथ आने से बचाएगा। लेकिन ऐसे कार्ड निपटाए जाते हैं। यह उल्लेखनीय है कि हस्ताक्षर रिकॉर्ड 1506 हैं, और कार्टोग्राफिक वस्तुओं के रिकॉर्ड - 1505 (छवि 1), हालांकि चित्रा 4 में लोडिंग प्रक्रिया के दौरान, शेप 2 एसक्यूएल ने यह भी लिखा है कि 1506 थे। एक और कहां गया? यह पता चला है कि एक भ्रूण उपकरण का एक गंभीर चूक है, जिसमें इस तथ्य में शामिल है कि यह खाली वस्तुओं को लोड नहीं करता है। वह उन्हें छोड़ देती है, और हस्ताक्षर की तुलना में कार्टोग्राफिक ऑब्जेक्ट्स की तालिका में कम प्रविष्टियां हो सकती हैं, जिसके परिणामस्वरूप भौतिक पत्राचार गोड्डम दादी के पास जाता है, क्योंकि आपको कैसे पता चलता है कि संख्या द्वारा कौन सी प्रविष्टियां बीच से फेंक दी गई थीं? यह अच्छा है कि यहां क्षेत्र और परिधि कॉलम हैं, जिनसे आप हस्ताक्षर में खाली मानचित्र ऑब्जेक्ट के अनुरूप रिकॉर्ड की गणना कर सकते हैं:


चित्र 7

हस्ताक्षर तालिका से इसे हटाने के बाद, delete from regions2010_captions where area = 0 दिखाई देता है। यह अच्छा है कि हमारे मामले में यह अंतिम लेकिन एक रिकॉर्ड था, और अंतिम की संख्या को केवल 1505 में बदला जा सकता है। फिर भी, हम ईमानदारी से बढ़ते आईडी में एक नए आदेश में हस्ताक्षर की संख्या को देखते हैं, जैसे कि कई हटाए गए रिकॉर्ड थे और वे कहीं बिखरे हुए होंगे। बीच में।

alter table regions2010_captions add n int
;with cte as (select *, row_number() over (order by ID) n1 from regions2010_captions) update cte set n = n1
alter table regions2010_captions drop column ID

मैं वास्तव में आशा करता हूं कि क्षेत्रों की भौतिक स्कैनिंग के क्रम में, क्षेत्रों2010_albsib तालिका भी नीलसन प्राप्त करती है। फिर इसे क्षेत्रों के साथ ID2010_captions तालिका (फ़ील्ड n) के साथ जोड़ा जा सकता है। कार्टोग्राफिक ऑब्जेक्ट की संख्या वास्तविक क्षेत्रों से अधिक है, क्योंकि जैसा कि हम हस्ताक्षर से देखते हैं, एक क्षेत्र में एक से अधिक बहुभुज शामिल हो सकते हैं: मास्को और ज़ेलेनोग्राड, तटीय क्षेत्र और द्वीप, आदि। उदाहरण के लिए, select s.geom from regions2010_captions c join regions2010_albsib s on s.ID = cn where c.region = ' ' :


चित्र 8

इसलिए, उन्हें उसी क्षेत्र के आधार पर जोड़कर जोड़ा जाना चाहिए:

if exists (select 1 from sys.sequences where name = 'Seq') drop sequence Seq
create sequence Seq as int start with 1 increment by 1

if object_id('regions2010', 'U') is not null drop table regions2010

select (next value for Seq) ID, c.region, geometry::UnionAggregate(s.geom) geom
into regions2010
from regions2010_captions c join regions2010_albsib s on s.ID = cn
group by c.region

select * from regions2010

यह कोड SQL Server 2012 में दो नई विशेषताओं को दिखाता है: अनुक्रम और भू-स्थानिक समुच्चय।


चित्र 9

रंग से देखते हुए, यह सच लगता है। इसलिए, हस्ताक्षर और बहुभुज को जोड़ने के साथ, हमसे गलती नहीं हुई।
इस तथ्य के बावजूद कि एक ही क्षेत्र से संबंधित बहुभुजों को एक ही रंग से चित्रित किया जाता है, वे अलग-अलग वस्तुएं हैं, जिनमें से एक ही ज्यामितीय संग्रह क्षेत्र (बहुभुज) है। आप क्वेरी का select ID, Region, geom.STGeometryType(), geom.STNumGeometries() from regions2010 :


चित्र १०

चयनित क्षेत्र के लिए संग्रह से तालिका में उन्हें तैनात करने के लिए, मैंने एक सरल कार्य लिखा:

if object_id('dbo.ufnGeometries', 'TF') is not null drop function dbo.ufnGeometries
go
create function dbo.ufnGeometries(@g geometry) returns @t table (id int identity(1, 1), g geometry) as begin
declare @i int = 1
while @i <= @g.STNumGeometries() begin
insert @t (g) values (@g.STGeometryN(@i))
set @i += 1
end
return
end

जो मूल रूप से, मूल रूप से ,2010 की तालिका से मूल क्षेत्रों में फिर से लौटने की अनुमति देता है।

select r.ID, r.Region, f.ID, fg from regions2010 r cross apply dbo.ufnGeometries(geom) f


चित्रा 11

वास्तव में, यह वास्तव में उसका नहीं है। ध्यान दें कि पहले से ही 1503 विस्तृत बहुभुज हैं, 1505 नहीं। मैं इस तथ्य से समझाता हूं कि एक ही क्षेत्र में दो जोड़े बहुभुज या तीन प्रतिच्छेदित बहुभुज हैं, जिन्हें एक में जोड़ा गया था। आइए देखें:

with cte as (select s.ID, s.geom, c.region from regions2010_albsib s join regions2010_captions c on s.ID = cn)
select cte1.ID, cte1.geom, cte2.ID, cte2.geom, cte1.geom.STIntersection(cte2.geom).STGeometryType() from cte cte1 join cte cte2 on cte1.ID < cte2.ID and cte1.region = cte2.region
where cte1.geom.STIntersects(cte2.geom) = 1

दरअसल, कुल सात ऐसे चौराहे हैं, जिनमें से 5 बिंदु चौराहे हैं, जो संयुक्त होने पर, एक बहुभुज में नहीं गिरते हैं, और दो जोड़े आसन्न बहुभुज, अर्थात्। टूटी हुई रेखा के रूप में एक आम सीमा होना। इनमें से प्रत्येक जोड़े को एक बहुभुज में जोड़ा गया था, इसलिए उनकी संख्या 2 से कम हो गई थी:


चित्र 12

निष्कर्ष में, हम कम () विधि के व्यावहारिक प्रभाव पर विचार करते हैं। कृपया इस तथ्य पर ध्यान दें कि 2-3 मिनट की देरी के साथ क्षेत्रों2010 का नक्शा (चित्र 9) तैयार किया गया है। यह एक महत्वपूर्ण संख्या के कारण होता है - हजारों या दसियों हज़ारों के क्रम से - प्रत्येक क्षेत्र के समोच्च को परिभाषित करता है। रिपोर्ट में दिखाए जाने वाले आकार के नक्शे के लिए, इस तरह की छान-बीन अतिश्योक्तिपूर्ण है, अनुचित रूप से प्रतिपादन समय को बढ़ाती है। मोटे करने के लिए, Reduce () विधि का उपयोग किया जाता है, जो डगलस-पीकर एल्गोरिथम के अनुसार काम करता है। पैरामीटर का मूल्य जितना बड़ा होगा, उतना ही मजबूत होगा (हमने यहां सहिष्णु तरीकों की जांच की)। उदाहरण के लिए, select ID, Region, geom.STNumPoints() [Points Total], geom.Reduce(100).STNumPoints() Reduce100, geom.Reduce(500).STNumPoints() Reduce500, geom.Reduce(1000).STNumPoints() Reduce1000 from regions2010


चित्रा 13

इस प्रकार, select ID, Region, geom.Reduce(1000) from regions2010 लिए दिखाई देने वाले गुणवत्ता के नक्शे को सहेजते समय select ID, Region, geom.Reduce(1000) from regions2010 तेजी से select ID, Region, geom.Reduce(1000) from regions2010 प्रदान किया जाता है:


चित्र 14