7-परम डेटम ट्रांसफ़ॉर्म के परिणामों के बीच भारी अंतर


12

मैं एक WGS84 lat / lon बिंदु को बदलने की कोशिश कर रहा हूं

11d42'32.10629" E 5d12'56.75371" S

तीन अलग-अलग सॉफ्टवेयर पैकेजों (प्रोज 4, जियोट्रांस और लेईको जियोऑफिस) में लेकिन मुझे परिणामों के बीच बहुत बड़े अंतर मिलते हैं - एक्स और वाई में कुछ सौ मीटर के आसपास! ध्यान दें कि ये भिन्नताएँ केवल 7 मापदंडों के साथ दिखाई देती हैं 3 के साथ नहीं।

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + से + proj = utm + ज़ोन = 32 + ellps = clrk80 + Towgs84 - 178.3,316.7,131.5, -5.278, -6.077, -6.977, -10.9792, -19.966.119 दक्षिण + इकाइयाँ = m + no_defs


GeoTrans

Delta X 178.3 
Delta Y 316.7 
Delta Z 131.5 
Rot X -5.278 
Rot Y -6.077 
Rot Z -10.9792 
SF = -19.166 / 0.999980834 (0.000019166)

लेईको जियोऑफिस

स्क्रीनशॉट


परिणाम क्रमशः हैं:

  1. 800392 9422525
  2. 800306 9422840 है
  3. 800941 9422891

मेरा मानना ​​है कि सभी 3 पैकेज 7-परम ट्रांसफॉर्म (बर्सा-वुल्फ विधि) के लिए समान गणितीय तरीकों का उपयोग करते हैं। इस भारी भिन्नता का कारण क्या हो सकता है?


जियो ऑफिस में, क्लार्क 1880 IGN है जो एक कस्टम परिभाषा है? क्या आप इसके पैरामीटर पोस्ट कर सकते हैं कि यह है या नहीं?
mkennedy

मैंने अन्य लोगों से लीका स्क्रीनशॉट प्राप्त किया है और वास्तव में अभी तक नहीं जानता।
विंडराइड सेप

जवाबों:


9

सबसे पहले, प्रोज 4 का उपयोग करता है जो ईपीएसजी 7 पैरामीटर विधि के "स्थिति वेक्टर" संस्करण को कहते हैं। यह संभव है कि जियोट्रांस और लेईको जियोऑफिस दूसरे संस्करण का उपयोग करें जिसे ईपीएसजी ने "कोऑर्डिनेट फ्रेम" कहा है। दोनों विधियां समतुल्य हैं, लेकिन रोटेशन मेट्रिक्स अलग हैं और कोणीय मापदंडों के संकेतों को बदलना होगा।

दूसरा, लेईको जियोऑफिस में परिवर्तन परिभाषा के स्क्रीन शॉट प्रदान करने के लिए धन्यवाद। या तो Proj4 परिभाषा या यह परिभाषा गलत तरीके से परिभाषित की गई है। परिभाषित मापदंडों एलीपोसिड ए और एलीपोसिड बी को स्विच किया जाना चाहिए। वर्तमान में, यह परिवर्तन WGS84 से कांगो 60 में परिवर्तित हो गया। Proj4 में, ध्यान दें कि विकल्प + Towgs84 है, इसलिए इसे FROM Congo60 TO WGS84 के रूप में परिभाषित किया गया है। परिभाषा को दिशा में बदलने के लिए, सभी मापदंडों के संकेत बदलें। जियोऑफ़िस की मदद भी देखें और देखें कि क्या एसएफ पैरामीटर प्रति मिलियन संस्करण या पहले से परिवर्तित-टू-स्केल-फैक्टर मान चाहता है या नहीं।

मुझे जियोट्रांस के बारे में पता नहीं है - क्या आपका मतलब एनजीए सॉफ्टवेयर से है? वैसे भी, उम्मीद है कि आपको मैच के लिए जियोऑफिस और प्रोज 4 मिल सकते हैं।


1
वाई (नथिंग / लेटिट्यूड) मूल्य हो सकता है क्योंकि क्लार्क 1880 मूल्य काफी समान नहीं हैं, लेकिन 60 मीटर बहुत बड़ा लगता है।
mkennedy

1
@mkennedy: क्या आप वही गणना कर सकते हैं जो मैंने आर्कगिस में की थी कि क्या हम लीका के बिना संरेखित कर सकते हैं?
आंद्रेजे सेप

1
@AndreJoost, ज़रूर। ईपीएसजी 1802 और 28232 के एसरी संस्करणों का उपयोग करते हुए, मुझे पूर्व मिला: 800230.139 नॉर्थिंग: 942313.413।
mkennedy

1
मैंने यह देखने के लिए कई वेरिएंट आज़माए कि क्या मैं लाइका परिणाम को पुन: पेश कर सकता हूं और करने में सक्षम नहीं था। वही समस्या जो आपको मिली: Y 60+ m से दूर।
mkennedy

1
मुझे विश्वास है कि Proj4 (और ArcGIS!) समाधान सही हैं। प्रोज 4 स्रोत खुला है। हालांकि ArcGIS कोड खुला नहीं है, मैंने लगभग 15 वर्षों तक इससे निपटा है। मैंने Geotrans 3.0 में एक कस्टम परिवर्तन भी जोड़ा और Proj4 / ArcGIS के समान परिणाम मिला।
mkennedy

4

प्रोज 4 के लिए, मापदंडों के संकेतों को उलट करना होगा।

यह परिभाषा पृष्ठ देखें:

http://www.spatialreference.org/ref/epsg/62826405/prettywkt/

GEOGCS["Pointe Noire (deg)",
    DATUM["Congo 1960 Pointe Noire",
        SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936269,
            AUTHORITY["EPSG","7011"]],
        TOWGS84[-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849],
        AUTHORITY["EPSG","6282"]],
    PRIMEM["Greenwich",0.0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.017453292519943295],
    AXIS["Geodetic latitude",NORTH],
    AXIS["Geodetic longitude",EAST],
    AUTHORITY["EPSG","62826405"]]

मुझे यकीन नहीं है कि आपको अपने मूल्य कहां से मिले।

क्यूजीआईएस ईपीएसजी को परिभाषित करता है: 28232 प्वाइंट नायर यूटीएम 32 एस:

+proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs

और 7 मापदंडों का परिणाम उससे दूर नहीं होना चाहिए।


संपादित करें

GDAL 1.10 के साथ, मुझे निम्नलिखित परिणाम मिले:

4326-proj-3 parameters:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs
800232.21   9423131.96 -1.76
4326-proj-7 parameters from spatialrefrence.org:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849 +units=m +no_defs
800230.13   9423133.46 91.31
4326-proj-7 parameters from proj4 datum_shift.csv (EPSG:1802):
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,19.166 +units=m +no_defs
800230.13   9423133.40 -5.72

इसलिए 3 मीटर से भी कम। ध्यान दें कि ई और एस को रिक्त स्थान के बिना डिग्री का पालन करना चाहिए।

और, सिर्फ रिकॉर्ड के लिए: टॉवेल 84 को शून्य परिणामों में सेट करना:

800310.94   9422829.37 -109.32

पॉइंट नायर डिग्री से परिवर्तित करते समय EPSG: 4262 से प्वाइंट नायर UTM में परिणाम:

800311.21   9422892.49 0.00

जैसा कि आप देख सकते हैं, समान मापदंडों के साथ लीका से एक बड़ा अंतर (> 600 मीटर) है। मुझे यह जानने की आवश्यकता है कि ऐसा क्यों होता है और समान परिणाम कैसे प्राप्त करें। क्या यह प्रोज 4 में उपयोग किए जाने वाले / से भूस्थैतिक निर्देशांक में रूपांतरण की विधि से हो सकता है? इसके बहुत जटिल सूत्र हैं। हो सकता है कि लीका अलग विधि का उपयोग करें?
विंडराइड सेप

1
क्या Leica 3-पैरामीटर-रूपांतरण प्रदान करता है? और वे इसके लिए क्या गणना करते हैं?
आंद्रे जे सेप
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.