OpenLayers के साथ डेटलाइन पर बंटवारे Linestrings

9

कुछ साल पहले मैंने द इंटरनेशनल डेट लाइन रैप को पोस्ट किया था और @jdeolive ने सुझाव दिया था कि मैं डेटलाइन पर सुविधाओं को विभाजित करता हूं। इसलिए मैंने कोशिश की।

जब मैं के साथ मेरी उपग्रह ट्रैक विभाजित करने की कोशिश splitWith डेटलाइन मैं वापस पाने पर null। मुझे पता है कि मैं सही ढंग से विभाजित कर रहा हूं क्योंकि जब मैं ग्रीनविच लाइन पर विभाजित होता हूं तो मुझे अपेक्षित परिणाम मिलते हैं।

किसी को पता है कि मैं कैसे OpenLayers के साथ डेट लाइन के साथ प्रोग्राम लाइन लाइन को ठीक से विभाजित कर सकता हूं ? अगर आपके पास है तो मैं उदाहरण कोड चाहता हूं।

मैंने कोशिश की है, wrapDateLineलेकिन यह मेरी वेक्टर परत इस तरह होने के बावजूद वेक्टर परतों पर काम नहीं करता है:

vectorLayer = new OpenLayers.Layer.Vector("GroundTracks", {
    renderers: ['Canvas', 'VML'],
    wrapDateLine: true}); // <-- shoud be wraping.

यहां छवि विवरण दर्ज करें

यहाँ मेरा कोड है:

var features = [];
var format = new OpenLayers.Format.WKT({
    'internalProjection': map.baseLayer.projection,
    'externalProjection': prjGeographic
});
var satTrack = format.read("LINESTRING (95.538611 13.286511, 94.730711 16.908947, 93.901095 20.528750, 93.043594 24.145177, 92.150978 27.757436, 91.214579 31.364666, 90.223791 34.965899, 89.165364 38.560019, 88.022401 42.145679, 86.772901 45.721205, 85.387568 49.284424, 83.826433 52.832413, 82.033480 56.361087, 79.927797 59.864504, 77.388419 63.333664, 74.227306 66.754285, 70.139140 70.102478, 64.605267 73.335774, 56.712904 76.373458, 44.881134 79.052803, 26.939886 81.047314, 02.704174 81.839241, -21.686285 81.101751, -39.887660 79.141947, -51.906937 76.480894, -59.912477 73.452897, -65.514482 70.225089, -69.645366 66.880243, -72.834535 63.461797, -75.393132 59.994131, -77.512464 56.491789, -79.315407 52.963919, -80.884039 49.416549, -82.275114 45.853820, -83.529088 42.278691, -84.675583 38.693355, -85.736827 35.099503, -86.729876 31.498490, -87.668095 27.891443, -88.562176 24.279331, -89.420849 20.663020, -90.251389 17.043303, -91.059999 13.420926, -91.852092 09.796602, -92.632515 06.171020, -93.405728 02.544857, -94.175960 -01.081217, -94.947343 -04.706542, -95.724045 -08.330456, -96.510402 -11.952298, -97.311065 -15.571400, -98.131162 -19.187081, -98.976502 -22.798638, -99.853829 -26.405335, -100.771148 -30.006378, -101.738172 -33.600889, -102.766925 -37.187866, -103.872602 -40.766117, -105.074803 -44.334175, -106.399366 -47.890158, -107.881153 -51.431559, -109.568417 -54.954914, -111.529886 -58.455253, -113.866668 -61.925160, -116.733085 -65.353081, -120.374635 -68.720132, -125.199754 -71.993642, -131.916790 -75.113368, -141.772276 -77.960803, -156.750096 -80.294831, -178.475596 -81.673196, 156.248392 -81.611421, 135.042323 -80.136505, 120.556535 -77.748172, 111.014840 -74.872356, 104.485504 -71.737081, 99.775637 -68.454400, 96.208126 -65.081545, 93.391438 -61.649716, 91.089380 -58.177038, 89.152970 -54.674643, 87.484294 -51.149703, 86.016609 -47.607042, 84.702947 -44.050030, 83.509299 -40.481112, 82.410411 -36.902133, 81.387093 -33.314533, 80.424442 -29.719485, 79.510644 -26.117981, 78.636145 -22.510889, 77.793053 -18.898997, 76.974710 -15.283040, 76.175371 -11.663718, 75.389950 -08.041709, 74.613831 -04.417680, 73.842693 -00.792294, 73.072378 02.833789, 72.298749 06.459907, 71.517566 10.085391, 70.724342 13.709564, 69.914194 17.331733, 69.081655 20.951185, 68.220447 24.567170, 67.323194 28.178891, 66.381031 31.785476, 65.383084 35.385943, 64.315735 38.979152, 63.161579 42.563725, 61.897893 46.137940, 60.494337 49.699551, 58.909396 53.245525, 57.084691 56.771602, 54.935577 60.271560, 52.334964 63.735923, 49.084320 67.149569, 44.859585 70.487030, 39.107498 73.702694, 30.852243 76.709182, 18.420695 79.329532, -00.339911 81.212453, -25.028018 81.831766)");

var featGreenwichLine = format.read("LINESTRING(0 -89, 0 89)");
var featDateLine = format.read("LINESTRING(180 -89, 180 89)");

features.push(featGreenwichLine);
features.push(featDateLine);
features.push(satTrack);

var resultsGreenwich = satTrack.geometry.splitWith(featGreenwichLine.geometry);
var resultsDateLine = satTrack.geometry.splitWith(featDateLine.geometry);

console.log(resultsGreenwich); //<--RETURNS EXPECTED RESULTS.
console.log(resultsDateLine);//<--RETURNS NULL.

vectorLayer.addFeatures(features);

मेरा प्रश्न इस प्रश्न का डुप्लिकेट नहीं है क्योंकि वे जानना चाहते हैं कि इसे ogr2ogr में कैसे करना है

अपडेट करें:

यह वह विशिष्ट डेटासेट है जिसके साथ मैं काम करता हूं (24 घंटे का सैटेलाइट ट्रैक): लिनेस्ट्रिंग wkt यहां पाया जा सकता है ।

यहां छवि विवरण दर्ज करें

openlayers-2  algorithm  splitting  antimeridian 
— CaptDragon
स्रोत
आप किस ओपनर का उपयोग कर रहे हैं?
— प्लक्स
@Plux 2.13.1 (नवीनतम)
— Captrragon
1
उनके एपीआई के अनुसार, रैपडेटलाइन का उपयोग केवल आधार परत पर किया जाना चाहिए, इसलिए यह कोई आश्चर्य नहीं है कि यह वेक्टर परत पर काम नहीं करता है। हालांकि, मुझे नहीं पता कि इसे वेक्टर लेयर पर कैसे काम करना है। मैं खुद को इसी तरह के मुद्दे के साथ मल्टीप्लगॉन कर रहा हूं जो डेटलाइन को पार करता है।
— प्लक्स
1
@ प्लीक्स मेरे समाधान की
— कैप्ट्रैगन
अच्छा समाधान है। यह दुर्भाग्य से मेरी समस्या पर लागू नहीं है, मेरा मानना ​​है कि मेरी समस्या जियोसर्वर-साइड पर अधिक है, और मेरे बहुभुज में निर्देशांक हैं जो डेटलाइन के "दूसरी तरफ" हैं, जैसे -180.00000000000003 90.0000010000734869, जो जियोसर्वर में कुछ मुद्दों को बनाता है, जो मुझे विश्वास है। मुझे अपने मानचित्र पर इसे प्रदर्शित करने में समस्या नहीं है (मेरे पास एक wms है जो ऐसा करता है), लेकिन मैं इन्हें एक wfs-query पर geoserver पर फ़िल्टरबॉक्स के रूप में उपयोग करना चाहता हूं :)
— Plux

जवाबों:

2

समस्या यह है कि आपकी सुविधा OpenLayers परिप्रेक्ष्य से दिनांक रेखा को पार नहीं करती है, इसलिए आपकी विभाजन रेखा आपकी सुविधा को बाधित नहीं करती है। अपने डेटा से उदाहरण:

..., -178.475596 -81.673196, 156.248392 -81.611421,...

आप -178 से 156 तक जाते हैं, और यह OpenLayers परिप्रेक्ष्य से दिनांक रेखा को पार नहीं करता है। तिथि रेखा पर विभाजन के बजाय, आपको अपने न्यूनतम एक्स मान पर विभाजित होना चाहिए।

// Build the splitting line based on the min and max coordinates of the vector to split
var minX = 999999999;
var minY = -20037508.34 // minimum value of the spherical mercator projection
var maxY = 20037508.34  // maximum value of the spherical mercator projection
//Extract the minimum X from the data as bounds seems to be rounded.
for(var i=0; i<satTrack.geometry.components.length; i++) {
    if(satTrack.geometry.components[i].x < minX)
        minX = satTrack.geometry.components[i].x;
}
var pointList = [
    new OpenLayers.Geometry.Point(minX, minY),
    new OpenLayers.Geometry.Point(minX, maxY)
];
var featDateLine = new OpenLayers.Feature.Vector(
    new OpenLayers.Geometry.LineString(pointList)
);

मैंने यहां एक उदाहरण बनाया है जो आपके उपग्रह ट्रैक को 2 विशेषताओं में सफलतापूर्वक विभाजित करता है: http://jsfiddle.net/6XJ5A/

अब अपने अपडेट में कई लाइन के साथ डब्ल्यूकेटी का उपयोग करने के लिए, एक सीधी रेखा का उपयोग करने के बजाय, आपको पूरे डेटासेट के माध्यम से जाना चाहिए और डेटलाइन के पार जाने वाले सभी निर्देशांक के साथ अपनी विभाजन रेखा का निर्माण करना होगा। एक मल्टीलाइन के अंदर छोटी लाइन का निर्माण करके, आप सभी निर्देशांक पर विभाजित कर सकते हैं जो डेटलाइन के पार जाना चाहिए। यहाँ अद्यतन उदाहरण है: http://jsfiddle.net/Jc274/

और कोड:

// Build the splitting line based on the min and max coordinates of the vector to split
var pointList = [];
var lastPoint = satTrack.geometry.components[0];
//Extract the minimum X from the data as bounds seems to be rounded.
for (var i = 1; i < satTrack.geometry.components.length; i++) {
    if (Math.abs(satTrack.geometry.components[i].x - lastPoint.x) > 10000000) {
        pointList.push(satTrack.geometry.components[i]);
    }
    lastPoint = satTrack.geometry.components[i];
}

var lineList = [];
for(var i=0; i<pointList.length; i++) {
    lineList.push(new OpenLayers.Geometry.LineString([
        new OpenLayers.Geometry.Point(pointList[i].x, pointList[i].y-0.00001), 
        new OpenLayers.Geometry.Point(pointList[i].x, pointList[i].y+0.00001)
    ]));
}

var featDateLine = new OpenLayers.Feature.Vector(
new OpenLayers.Geometry.MultiLineString(lineList), null, split_style);

यह आपको सभी बिंदुओं पर एक स्प्लिटेड लाइन लौटाएगा जो डेटलाइन को "क्रॉस" करेगा

ध्यान दें कि मैं 2 निर्देशांक को जोड़ने के लिए मानचित्र पर जाने वाली रेखा को हटाने के लिए निर्देशांक के माध्यम से भी लूप करता हूं:

for (var i = 0; i < resultsDateLine.length; i++) {
    // Remove the first (or last) point of the line, the one that cross the dateline
    if (Math.abs(resultsDateLine[i].components[0].x - resultsDateLine[i].components[1].x) > 10000000) {
        resultsDateLine[i].removeComponent(resultsDateLine[i].components[0]);
    }
    if (Math.abs(resultsDateLine[i].components[resultsDateLine[i].components.length - 1].x - resultsDateLine[i].components[resultsDateLine[i].components.length - 2].x) > 10000000) {
        resultsDateLine[i].removeComponent(resultsDateLine[i].components[resultsDateLine[i].components.length - 1]);
    }
    features.push(new OpenLayers.Feature.Vector(resultsDateLine[i], null, style_array[i]));
}

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

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

— जूलियन-सैमुअल लैक्रोस
स्रोत
प्रयास के लिए +1, लेकिन आपके नक्शे उदाहरण समस्या को हल करने के लिए प्रतीत नहीं होते हैं। आपके उदाहरण मानचित्र पर उपग्रह ट्रैक अभी भी उपग्रह ट्रैक के प्राकृतिक पथ का अनुसरण करने के बजाय पूरे विश्व को पार करते हैं। मेरा कुछ छूट रहा है?
— कैप्टनड्रैगन
आप समस्या को गलत समझ सकते हैं क्योंकि ..., -178.475596 -81.673196, 156.248392 -81.611421,...बिल्कुल डेटलाइन को पार कर जाता है। यहां देखें
— CaptDragon
मुझे कोड और मेरे स्पष्टीकरण को अपडेट करने दें। मुझे पता है कि इसे डेटा लाइन को पार करना चाहिए, लेकिन OpenLayers इसका समर्थन नहीं करता है। OL के दृष्टिकोण से, यह डेटलाइन को पार नहीं कर रहा है।
— जूलियन-सैमुअल लैक्रिक्स
ये सही है। यही समस्या है। मुझे OpenLayers को ट्रिक करने और लाइन को विभाजित करने की आवश्यकता होगी ताकि यह किनारे तक सही जाए और फिर दूसरी तरफ जारी रहे जहां यह माना जाता है।
— कैप्टनड्रैगन
2

मुझे @ गेनथ गीथब पर एक महान समाधान मिला। इसे Arc.js कहा जाता है और यह महान सर्कल मार्गों की गणना के लिए है। इतना ही नहीं, यह डेटलाइन पर लाइन को विभाजित भी करेगा और आपको दो लाइनस्ट्रेस प्रदान करेगा जो डेटलाइन पर मिलते हैं, जिसे ओपन लाइयर आसानी से मैप कर सकते हैं। मुझे आशा है कि वह इनाम का दावा करने के लिए आगे आएगा।

यहाँ मेरे परिणाम हैं:

यहां छवि विवरण दर्ज करें यहां छवि विवरण दर्ज करें

— CaptDragon
स्रोत
1

SplitWith समारोह पृथ्वी की 3 आयामी आकार के बारे में पता नहीं है। यह केवल 2-आयामी दुनिया में संचालित होता है। आपके मामले में, आपके सभी LINESTRINGX निर्देशांक -180 और 180 के बीच हैं। इसलिए OpenLayers के दो-आयामी परिप्रेक्ष्य से, लाइन स्ट्रिंग वास्तव में कभी भी आपकी विभाजन ज्यामिति (तिथि रेखा) को पार नहीं करती है और यह आपको वापस लौटकर बताती है null।

मेरा मानना ​​है कि बंटवारे के लिए आपको कुछ कस्टम कोड लिखने होंगे। मैं एक एल्गोरिथ्म की कल्पना कर रहा हूं, जो आपके कोने पर लूप करता है, इस तरह आउटपुट लाइन स्ट्रिंग्स का निर्माण करता है:

  • प्रत्येक आसन्न जोड़ी के लिए तय करें कि उनके बीच का खंड तिथि रेखा को पार करता है या नहीं।
  • यदि यह नहीं है, तो उस खंड को इस तरह रखें और इसे "वर्तमान" आउटपुट लाइन स्ट्रिंग में जोड़ें।
  • यदि ऐसा होता है , तो खंड को दो भागों में विभाजित करें। "करंट" लाइन स्ट्रिंग में एक भाग जोड़ें, एक नया "करंट" लाइन स्ट्रिंग शुरू करें, और दूसरा भाग इस नए में जोड़ें।

यह निर्धारित करने के लिए एक उचित अनुमानी है कि यदि जोड़ी की एक रेखा तिथि रेखा को पार करती है, तो यह देखना है कि क्या एक्स निर्देशांक का अंतर 180 डिग्री से अधिक है। (हालांकि, यह गलत हो सकता है, उदाहरण के लिए, ध्रुवीय क्षेत्रों में। हो सकता है कि आप बहुत भाग्यशाली हों, जिनके पास वास्तव में कोई उच्च अक्षांश न हो।)

एक खंड को दो भागों में विभाजित करने का कार्य रैखिक प्रक्षेप के रूप में सरल हो सकता है (यदि आप ट्रैक की सटीकता के बारे में बहुत अधिक परवाह नहीं करते हैं)। जब आप यह पता लगा लेते हैं कि सेगमेंट की तारीख लाइन को पार कर जाती है तो आप दूसरे वर्टेक्स की एक कॉपी बनाते हैं और उसके एक्स कोऑर्डिनेट (360 डिग्री जोड़कर या घटाकर) ले जाते हैं और फिर वाई कोऑर्डिनेट करते हैं।

संपादित करें : यहाँ एक JSFiddle है जो आपके डेटा पर उपरोक्त एल्गोरिथ्म को प्रदर्शित करता है: http://jsfiddle.net/85vjS/

— क्रिस्टोफ़र स्ट्रीट डे
स्रोत
0

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

var featDateLine = format.read("LINESTRING(179.99 -89, 179.99 89)");

और शायद

var featDateLine = format.read("LINESTRING(-179.99 -89, -179.99 89)");

दूसरी तरफ के लिए।

एक अन्य समाधान CRS में काम करना है जो डेटलाइन पर "आउट ऑफ बाउंड" नहीं है। फिर आपको समस्या के बिना अपने डेटा को विभाजित करने में सक्षम होना चाहिए।

— radouxju
स्रोत
मैंने पहले ही कोशिश की है LINESTRING(179.99 -89, 179.99 89)और LINESTRING(-179.99 -89, -179.99 89)कोई फायदा नहीं हुआ है। सीआरएस के बारे में, दुर्भाग्य से, यह मेरे उद्देश्य के लिए काम नहीं करेगा क्योंकि मैं कई बार दुनिया भर में जाने वाले सैटेलाइट ट्रैक्स की मैपिंग कर रहा हूं। इसलिए सभी CRS कहीं अलग हो गए हैं और मेरे पास एक ही मुद्दा होगा कहीं भी मैं इसे विभाजित कर दूंगा। आपके सहयोग के लिए धन्यवाद।
— CaptDragon
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.