नेविगेटर में OSM और स्पीड बम्प मैप

चूंकि शुरुआत से शुरू करना दिलचस्प नहीं है, मैं अंत से शुरू करूंगा।

हमने कर दिखाया। हमें OpenStreetMap डेटा से स्पीड बंप मिला, उन्हें भयानक वाणिज्यिक कार्यक्रम नविटेल के साथ पार किया, इन गति धक्कों का एक वेब दर्शक बनाया, और http://latlon.org/tc/ पर शुरुआती के लिए उन्हें जोड़ने के लिए इंटरफ़ेस। और उन्होंने एक छोटी प्रेस विज्ञप्ति भी लिखी, जिसका एक लिंक मित्रों और परिचितों, मोटर चालकों को भेजा जा सकता है।

लेकिन हब्र के लिए, कोई व्यक्ति कुछ विशेष बता सकता है: यह सब अंदर कैसे काम करता है, और यह कैसे किया गया।



OpenStreetMap का मुख्य नियम मज़े करना है। वर्णित सब कुछ पैसे या कुछ अन्य अच्छाइयों के लिए नहीं किया गया था, और अपनी सुविधा के लिए भी नहीं किया गया था।



मिन्स्क मिनी-ओएसमोवकों में से एक (मिनी, क्योंकि कहीं भी, आईआरसी चैनल को छोड़कर , उन्हें घोषित नहीं किया जाता है) मुझे कार्य से सम्मानित किया गया था: “क्या आप एक उपयोगी काम करना चाहते हैं? ऑस्मोसिस से नाविक में लॉजर्स प्राप्त करें। " तब सब कुछ किसी तरह भूल गया था, और कुछ महीनों के बाद, मुझे अभी भी याद है, मंच w3bsit3-dns.com से कॉमरेड stas_symba के वीर कार्य के लिए धन्यवाद। अब कई महीनों के लिए, वह अकेले बेलारूस में स्पीड बम्प्स पर अपडेट इकट्ठा और अपलोड कर रहा था। "लेकिन आखिरकार, ऐसे डेटा को ओएसएम में तुरंत इकट्ठा करना आसान है," मैंने सोचा, और मैं निर्यात के लिए बैठ गया।

प्रारूप

चूँकि मेरे पास नाविक, नेवीटेल या कार नहीं है, इसलिए मुझे आँख बंद करके, मंचों को पढ़ने और अपने दिमाग के साथ सोचने के लिए सब कुछ करना था। आपदा के पैमाने को शुरू करने और समझने के लिए - मंचों से डेटा के प्रारूप के बारे में हम क्या हासिल करने में कामयाब रहे, इसके बारे में संक्षेप में।



Cops.txt फ़ाइल (और / या speedcam.txt) एक नियमित सीएसवी फ़ाइल है, पहली पंक्ति में हेडर होता है, और सभी बाद वाले डेटा होते हैं। शीर्षक:



idx,x,y,type,speed,dirtype,direction





जहां:

• idx - अद्वितीय वस्तु संख्या;
• x, y - निर्देशांक, प्रक्षेपण देशांतर-अक्षांश EPSG में: 4326;
• गति - वस्तु के लिए अनुमत गति;
• गंदगी का प्रकार - कार्रवाई की दिशा का प्रकार - 1 या 2 दिशाओं में;
• दिशा - कार्रवाई की दिशा;

उपकरण चयन

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



Osm2pgsql का उपयोग करते हुए होम सर्वर पर शुरू करने के लिए, मैपनिक के लिए सामान्य आधार आयात किया गया था। पहला अनुरोध जल्दी और सरल रूप से लिखा गया था:

select * from planet_osm_point where traffic_calming is not NULL;





इसके जवाब में, WKB प्रारूप में ज्यामिति सहित कई क्षेत्रों को मुझे वापस कर दिया गया, जो हेक्स में एन्कोडेड था। अच्छा बिल्कुल नहीं।



मुझे PostGIS माहोल में जाना पड़ा, जिसमें ST_X और ST_Y फ़ंक्शन पाए गए। ऐसा लगता था कि यह आवश्यक था। पुनर्लेखन:

gis=> select st_x(way), st_y(way) from planet_osm_point where traffic_calming is not NULL limit 1;

st_x | st_y

------------------+------------------

3085590.21426068 | 7159526.18035388

(1 )









अप्रत्याशित रूप से? उम्मीद है, लेकिन हमें एक स्पष्टीकरण बनाने की आवश्यकता है।



वेब मैपिंग में, दो अनुमानों का व्यापक रूप से उपयोग किया जाता है: उपयोगकर्ता के लिए ट्रांसमिशन और डिस्प्ले के लिए - अक्षांश-देशांतर, यह EPSG: 4326 और "Google" भी है, यह "स्पेक्टर पर व्यापारी" भी है, यह EPSG: 3857 है, यह EPSG: 900913, भी है। वह ईपीएसजी: 3785 है (इतने कोड क्यों हैं? बड़े निगमों, शौकियों और रजिस्टरों के बीच विवादों का एक लंबा इतिहास, एक अलग पद के योग्य)। यह प्रक्षेपण अच्छा है कि किसी भी भाषा में इससे 4326 तक संक्रमण शक्ति से गणित की दो लाइनें लेता है। यह है, मीट्रिक, और कोणीय निर्देशांक में नहीं, जिसका उपयोग स्क्रीन पर मानचित्र प्रदर्शित करने के लिए किया जाता है। और गणना की सुविधा के लिए, यह उस में है कि मेपनिक डेटाबेस संग्रहीत है।



ठीक है, ठीक है। Pereproetsiruem।

gis=> select ST_X(transform(p.way,94326)) as X, ST_Y(transform(p.way,94326)) as Y from planet_osm_point p where traffic_calming is not NULL limit 1;

x | y

------------+------------

27.7183285 | 53.9438334

(1 )







आगे क्या है? स्पीड।



OSM में कोई भी स्पीड बंप पर गति नहीं डालता है। हमें यह पता लगाना चाहिए कि इसे कहाँ प्राप्त करना है। हम पहले व्यक्ति से पूछते हैं जो सामने आया (यह वाइल्डमैन निकला ), हमें जवाब मिलता है - एसडीए के अनुसार, शहर में अनुमत गति 60 है, शहर के बाहर 90 है। और सड़क पर भी एक गति सीमा हो सकती है, जहां इसे संबंधित टैग में दर्ज किया जाएगा - मैक्सस्पीड। ठीक है, उसी समय, आपको वनवे टैग से एकतरफा होने का संकेत मिल सकता है - व्यर्थ में, शायद, प्रारूप में, इसके लिए क्षेत्र आरक्षित है?



बहुभुज की तालिका के साथ संयोजन करें - यह प्रशासनिक सीमाओं को संग्रहीत करता है, और लाइनों की तालिका - यह सड़कों की लाइनों को संग्रहीत करता है।



select ST_X(transform(p.way,94326)) as X,

ST_Y(transform(p.way,94326)) as Y, '102' as TYPE,

case when l.maxspeed is not NULL then l.maxspeed else case when t.admin_level = '8' then '60' else '90' end end as SPEED,

case when l.oneway = 'yes' then '1' else '2' end as DirType



from (planet_osm_point p

join planet_osm_line l on (l.highway is not NULL and ST_DWithin(p.way,l.way,.1)))

LEFT OUTER JOIN planet_osm_polygon t on (t.admin_level = '8' and ST_Within(p.way, t.way))

where p.traffic_calming is not NULL;









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



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

डर और प्रलेखन की कमी

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

"ठीक है, चूंकि सभी कार्ड मर्केटर प्रोजेक्शन में हैं, और यह कोण और दिशाओं को संरक्षित करता है , तो वे इसमें गिने जाते हैं!"



अजीमथ की गणना करने के लिए, PostGIS को दो बिंदुओं की आवश्यकता होती है। जाहिर है, उन्हें उसी सड़क पर झूठ बोलना चाहिए जिस पर स्पीड बंप स्थापित किया गया है, उसके दोनों तरफ कुछ दूरी पर।

ST_Line_Locate_Point विधि ने रेखा को लंबाई को उसके आरंभ से बिंदु तक खोजने की अनुमति दी, और ST_Line_Interpolate_Point - इसके विपरीत, बिंदु को रेखा के साथ लंबाई का पता लगाने।



जब सीएसवी की जरूरत होती है तो मेरे अनुरोध के जवाब में साइकल एक छद्म विज्ञान की गोली खींचता है! क्या रैपर लिखना संभव है?

यह पता चला कि सब कुछ बहुत सरल है। COPY फ़ंक्शन आपको CSV को तुरंत क्वेरी का परिणाम देने की अनुमति देता है - बस आपको जो चाहिए। यहाँ इस तरह के एक छोटे बिल्ली खुशी है।



मैं पहला परिणाम निकालता हूं, मुझे नवितेल का मालिक मिलता है, मैं उसे जांचने के लिए मजबूर करता हूं ... और उफ। दिशाएं सड़कों के समानांतर नहीं हैं। सोचा "वास्तव में ..." उठता है, और पंजे खुद को अज़ीमथ की गणना करते समय 4326 में पुन: प्रक्षेपण लिखते हैं।



हां, बहादुर नविटेल डेवलपर्स एक प्रक्षेपण में दिशा पर विचार करते हैं जो दिशा को संरक्षित नहीं करता है। उनका सम्मान, प्रशंसा और सम्मान।

Kvaziitog

नतीजतन, रिक्वेस्ट सिर्फ एक ऐसा राक्षस था, जो OpenStreetMap डेटा से निकालता है कि वे शुरू में ऐसा नहीं करते थे - नविटेल प्रारूप में गति धक्कों। आप एक अनुरोध लिख सकते हैं जो होमवर्क के रूप में उनसे स्पीड कैमरा प्राप्त करता है। :)

CREATE TEMP SEQUENCE idx ;

COPY (



select nextval('idx') as idx,

ST_X(transform(p.way,94326)) as X,

ST_Y(transform(p.way,94326)) as Y, '102' as TYPE,

case when l.maxspeed is not NULL then l.maxspeed else case when t.admin_level = '8' then '60' else '90' end end as SPEED,

case when l.oneway = 'yes' then '1' else '2' end as DirType,

floor(ST_Azimuth(

transform(ST_Line_Interpolate_Point(l.way,

case when (ST_Line_Locate_Point(l.way,p.way)-0.01 < 0) then 0 else ST_Line_Locate_Point(l.way,p.way)-0.0000001 end),94326),

transform(ST_Line_Interpolate_Point(l.way,

case when (ST_Line_Locate_Point(l.way,p.way)+0.01 > 1) then 1 else ST_Line_Locate_Point(l.way,p.way)+0.0000001 end)

,94326))/(2*pi())*360) as Direction

from (planet_osm_point p

join planet_osm_line l on (l.highway is not NULL and ST_DWithin(p.way,l.way,.1)))

LEFT OUTER JOIN planet_osm_polygon t on (t.admin_level = '8' and ST_Within(p.way, t.way))

where p.traffic_calming is not NULL

) TO STDOUT WITH CSV HEADER;









लेकिन अगर उनकी बिल्ली के कान हर जगह से नहीं निकलते तो कोटियारा कोटियारा नहीं होता। इसलिए, सीआईएस के लिए अद्यतन गति धक्कों और कैमरों के पांच मिनट उतारने के अलावा, टाइल्स में इस सभी अपमान का एक प्रतिपादन आयोजित किया गया था, ताकि जिन लोगों के पास नेवीटेल के साथ नेविगेटर नहीं है, वे वार्म-अप का आयोजन कर सकें।



सप्ताहांत के करीब, कॉमरेड एंड्रीश (मुझे लगता है कि उनके पास यह भी बताने के लिए कुछ है कि क्या आप उनसे पूछते हैं :) ने पंजीकरण के बिना नई गति के धक्कों की रिपोर्टिंग के लिए एक सुविधाजनक इंटरफ़ेस लिखा था - latlon.org/tc , और यह सभी सेवाओं को जनता पर डालने का निर्णय लिया गया था। मुझे आशा है कि आप इसे पसंद करेंगे।



मजा आ गया! ;)