क्या मैं Google से Mercator पहेली हल कर सकता हूं?

हाल ही में, Google से एक शांत मानचित्रण पहेली के बारे में एक लेख आया था। लगभग 20 मिनट इसे हल करने के बाद, मैं यह समझना चाहता था कि क्या यह स्वचालित रूप से करना संभव है?



पहली बात यह है कि सभी बहुभुज प्राप्त करें। पृष्ठ का स्रोत कोड पढ़ने और डेवलपर कंसोल खोलने के बाद, निम्न कोड का जन्म हुआ:

for (polyIdx=0;i<polys.length;polyIdx++){ var paths = polys[polyIdx].getPaths().getArray(); var arMultiPolygons=[]; for (j=0;j<paths.length;j++){ var arPolygons=[]; for (i=0;i<paths[j].length;i++) { arPolygons.push([paths[j].getArray()[i].lng(),paths[j].getArray()[i].lat()]); } arMultiPolygons.push(p) } } 

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

 from selenium import webdriver from shapely.geometry import polygon,multipolygon from shapely.wkt import loads import psycopg2 con = psycopg2.connect(database="WB",user='postgres',password='pass', host='127.0.0.1', port='5432') cur = con.cursor() #  ur.execute("create table googlePolygons ( idx integer, googlegeom geometry );") browser = webdriver.Firefox() browser.get("http://gmaps-samples.googlecode.com/svn/trunk/poly/puzzledrag.html") polysCount = browser.execute_script('return polys.length;') for iPoly in xrange(polysCount): element = browser.execute_script('var paths = polys[%s].getPaths().getArray();var arMultiPolygons=[];' 'for (j=0;j<paths.length;j++){ var arPolygons=[];' 'for (i=0;i<paths[j].length;i++) ' '{arPolygons.push([paths[j].getArray()[i].lng(),paths[j].getArray()[i].lat()]);} arMultiPolygons.push(p)}' 'return arMultiPolygons;'% str(iPoly)) polygons = [] #    for i in element: polygons.append(polygon.Polygon(i)) #   cur.execute("insert into googlePolygons (googlegeom,idx) values (geomfromewkt('"+multipolygon.MultiPolygon(polygons).wkt+"')," + str(iPoly)+");") con.commit() 

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

• एक बहुभुज लिया गया था, एक केन्द्रक को उसमें से हटा दिया गया था, फिर उस देश का केंद्रक जिसके साथ सामंजस्य किया गया था, जोड़ा जाएगा (दूसरे शब्दों में, वह बहुभुज को स्थानांतरित कर दिया गया)
• इसके अलावा, यह उस अनुपात पर विचार किया जाएगा जिसके द्वारा आपको बहुभुज की सीमा पर प्रत्येक बिंदु को सेंट्रो से ज़ूम इन / आउट करना होगा (प्रक्षेपण के कारण पैमाने में वृद्धि / कमी)
• क्षेत्र माना जाता था

और एल्गोरिथ्म सिर्फ इस तरह का है क्योंकि सुडौल नहीं जानता है कि भौगोलिक वस्तुओं के क्षेत्र की गणना कैसे करें (कम से कम मुझे नहीं मिला है कि कैसे)। इसलिए, सभी समान, मैंने पोस्टगिस को बंधन छोड़ने का फैसला किया, क्योंकि वहां आप भूगोल को भूगोल में बदल सकते हैं।
ठीक है, फिर सब कुछ सरल है: आपको इंटरनेट पर देशों के बहुभुज के साथ एक आकृति फ़ाइल ढूंढनी होगी, इसे पोस्टगिस पर अपलोड करना होगा और तुलना करना होगा।
आकृतियाँ मैंने यहाँ से ली हैं । उन्हें डेटाबेस में कैसे अपलोड किया जाए, इसे बहुत लिखा गया है, समस्या केवल एन्कोडिंग में हो सकती है, जो कि यूटीएफ -8 नहीं है।

सबसे सरल हिस्सा बना हुआ है: क्षेत्रों की तुलना करें

 cur.execute("select idx,st_area(geography(st_setsrid(googlegeom,4326)))/1000000, " "st_astext(st_centroid(st_setsrid(googlegeom,4326))) " "from googlePolygons order by idx;") unRealC = cur.fetchall() cur.execute("select name, st_area(geography(st_setsrid(geom,4326)))/1000000, " "st_astext(st_centroid(st_setsrid(geom,4326))) from tm_world_borders;") realC = cur.fetchall() for urc in unRealC: for rc in realC: if float(abs(rc[1]-urc[1])) /rc[1]<0.02: browser.execute_script('pl = new google.maps.Polyline(' '{path: [new google.maps.LatLng('+str(loads(rc[2]).y)+','+ str(loads(rc[2]).x)+'),' ' new google.maps.LatLng('+str(loads(urc[2]).y)+','+ str(loads(urc[2]).x)+')],' 'map: map});') 

मैंने float(abs(rc[1]-urc[1])) /rc[1]<0.02 - float(abs(rc[1]-urc[1])) /rc[1]<0.02 त्रुटि] सेट किया। एक पंक्ति को जोड़ा गया जो मूल बहुभुज के प्रतिशत को उसके संभावित वास्तविक स्थान के केंद्र से जोड़ता है। loads() - एक वस्तु के लिए WKT कन्वर्ट।
यहाँ क्या हुआ:

क्या काम नहीं किया: ग्रीनलैंड, दक्षिण अफ्रीका, आइसलैंड, थाईलैंड। अपनी डिस्क पर पाए गए एक और आकार का उपयोग करते हुए, मैं उसी 2% के साथ हासिल करने में कामयाब रहा जो केवल दक्षिण अफ्रीका और थाईलैंड की पहचान नहीं थी।

उपरोक्त सभी से, यह इस प्रकार है कि मैं सफल नहीं हुआ, इस तथ्य के कारण कि मुझे नहीं पता है कि Google किस डेटा द्वारा राज्यों की सीमाओं का निर्माण करता है, और यह भी ज्ञात नहीं है कि उन्होंने इन सीमाओं को कितना सरल बनाया।

जैसा कि मुझे लगता है कि परिणाम प्राप्त करना संभव होगा: बहुभुज के आकार की समानता के लिए एक चेक जोड़ें। लेकिन एक ही Postgis के बुनियादी तरीकों को नहीं किया जा सकता है। एल्गोरिथ्म मुझे लगता है कि यह है: ज्यामिति का एक मजबूत सरलीकरण, और फिर स्वीकार्य सीमा में कोणों की तुलना करना।