क्या एक नए प्रक्षेपण में बदलना, फिर वापस, डेटा सटीकता को प्रभावित करता है?

Question

मेरे पास NAD83 SC स्टेट प्लेन में एक फीचर क्लास (साउथ कैरोलिना काउंटियां, इतना बड़ा भौगोलिक क्षेत्र) है। इसे दूसरे प्रक्षेपण (NAD83 UTM 17) में बदलना होगा, फिर मूल में बदल दिया जाएगा। मैं इसे पूरा करने के लिए एस्री के प्रोजेक्ट टूल का उपयोग करूंगा ।

क्या यह दोहरी परिवर्तन पॉलीगोन निर्देशांक के स्थान में बदलाव का कारण बन सकता है, और कितना - सेंटीमीटर, मीटर, किलोमीटर से?

Comments

क्योंकि: परिवर्तन रिज़ॉल्यूशन, सिस्टम रिज़ॉल्यूशन अंतर और ज्यामिति संग्रहण रिज़ॉल्यूशन और सहिष्णुता का समन्वय करते हैं। इन "वैरिएबल्स" में से हर एक अलग है। इसलिए आपको प्रत्येक के लिए प्रलेखन पढ़ने की आवश्यकता है।

— GISI

... और यदि आप आर्कगिस का उपयोग कर रहे हैं तो संभवतः आपके डेटा के स्थानिक डोमेन के लिए उच्चतम रिज़ॉल्यूशन परिवर्तनों के क्रम में सूचीबद्ध सैकड़ों परिवर्तन समीकरण हैं।

— GISI

A -> B -> A 'का सामान्य परिणाम A ~ = A' है, लेकिन अगर आप इसे गलत करते हैं तो मिश्रण में डेटम परिवर्तन को वास्तव में बेईमानी से जोड़ा जा सकता है। एक महान सौदा इस बात पर निर्भर करता है कि समन्वय के संदर्भ कैसे परिभाषित किए जाते हैं (और इसलिए प्रत्येक समन्वय प्रणाली के मानचित्र इकाइयों में छंटनी)।

— विंस

Answer 1

मुझे नहीं पता कि कौन सा प्रोजेक्शन इंजन आर्कगिस का उपयोग करता है, लेकिन एक बहुत ही दिलचस्प सवाल भी proj.4 के लिए। तो मैं इसे GNU-R वातावरण में proj.4 प्रोजेक्शन इंजन का परीक्षण करने का प्रयास करता हूं। मैं एनएडी 83 - यूटीएम 17 कोनों और ईपीएसजी 26917 का उपयोग करता हूं और इसे 10000 और 1000000 बार पुनरावृत्ति करता हूं और शुरुआती मूल्यों के अंतर की गणना करता हूं।

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

ऐसा लगता है कि "अस्वीकृति" त्रुटि 10000 छोरों के लिए एक सेंटीमीटर सीमा के भीतर है।

"LON/LAT differences after  10000  loops"
           DLON          DLAT
1 -2.441464e-07 -1.341807e-07
2  2.441129e-07 -1.341807e-07
3  1.852679e-07 -1.691737e-08
4 -1.853157e-07 -1.691819e-08

"X/Y differences after  10000  loops"
            DX           DY
1 -0.025169783 -0.014338141
2  0.025166375 -0.014338208
3  0.002419045 -0.002016762
4 -0.002419690 -0.002016889

और अगर आप 1000000 बार लूप चलाते हैं, तो मीटर रेंज में एक त्रुटि हो जाती है।

"LON/LAT differences after  1000000  loops"
           DLON          DLAT
1 -2.441464e-05 -1.341845e-05
2  2.441128e-05 -1.341846e-05
3  1.852621e-05 -1.691837e-06
4 -1.853105e-05 -1.691828e-06

"X/Y differences after  1000000  loops"
          DX         DY
1 -2.5172288 -1.4339977
2  2.5168869 -1.4340064
3  0.2419201 -0.2017070
4 -0.2419859 -0.2017094

यहाँ स्क्रिप्ट है।

# load the package
require('proj4')

# the LON/LAT frame of NAD83 UTM 17 
lon = c(-84.00, -78.00, -84.00, -78.00 ) 
lat = c( 24.00,  24.00,  83.00,  83.00)

# build the projection conform object
ll0 = matrix(c(lon,lat),nrow=4,ncol=2)
xy0 = project(ll0,"+init=epsg:26917",ellps.default='GRS80')

# make a copy
ll1 = ll0
xy1 = xy0

# number of iterations
num = 1000000

# reproject the stuff num times
for(i in 1:num) {
 # project forward  
 xy1 = project(ll1,"+init=epsg:26917", ellps.default='GRS80')
 # project backward
 ll1 = project(xy1,"+init=epsg:26917", inverse=T, ellps.default='GRS80')
}

# build difference table ll
dll = as.data.frame(ll1-ll0)
names(dll) = c('DLON','DLAT')
# print results LON/LAT
print(paste("LON/LAT differences after ", num," loops"))
print(dll)

# build difference table xy
dxy = as.data.frame(xy1-xy0)
names(dxy) = c('DX','DY')
# print results X/Y
print(paste("X/Y differences after ", num," loops"))
print(dxy)

सांख्यिकीय वातावरण में आगे के परीक्षण आसान होने चाहिए। लिनक्स वातावरण के लिए स्क्रिप्ट और कोड स्पष्टीकरण github.com/bigopensky पर उपलब्ध हैं ।

Comments

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

— एरिका

क्या आप शामिल कर सकते हैं कि NAD83 UTM कोनों से आपका क्या मतलब है? यदि वे क्षेत्र के चरम पर हैं (उदाहरण के लिए उच्च अक्षांश), तो यूएस के भीतर बिंदुओं का उपयोग करने से संभवतः बेहतर परिणाम मिलेंगे।

— mkennedy

मुझे लगता है, कि WGS84 की सीमाएँ EPSG 26917 के साथ spatialreference.org/ref/epsg/nad83-utm-zone-17n पर भेज दी गई हैं । WGS84 Bounds: -84.0000, 24.0000, -78.0000, 83.0000यह ब्याज का सही क्षेत्र है। क्या मैंने गलती की है?

— हक्फिन 17

@ हुक्फिन दुह, मुझे कोड में मूल्यों को देखना चाहिए था! बेहूदे सवाल के लिए माफ कीजिएगा। UTM के बारे में सामान्यीकृत उत्तर के लिए महान मूल्य।

— mkennedy

Answer 2

Esri का अपना प्रक्षेपण इंजन है।

ब्याज के उपयुक्त क्षेत्र में उपयोग किए जाने पर अधिकांश अनुमानों और भौगोलिक / डेटा रूपांतरण विधियों का अच्छी तरह से व्यवहार किया जाता है। यदि आप किसी UTM ज़ोन से बहुत दूर हैं, तो अनुप्रस्थ मर्केटर हमेशा 'उलटा' नहीं होता है (अक्षांश-देशांतर में परिवर्तित)। पूरी दुनिया के लिए इस्तेमाल किए जाने वाले अनुमानों में डंडे या +/- 180 मेरिडियन या 'एंटी-मेरिडियन' (मेरिडियन जो अनुमानित समन्वित संदर्भ प्रणाली के केंद्र के विपरीत है) में कुछ मुद्दे हो सकते हैं।

मैंने 4 पॉइंट्स चलाए जो एसरी प्रोजेक्शन इंजन के माध्यम से दक्षिण कैरोलिना के बाहर आते हैं। 1k या 10k या 1M अंक के तनाव परीक्षण के लिए, मुझे कुछ कोड करने होंगे, जैसा कि मेरा मौजूदा समान परीक्षण सिर्फ एक 'राउंड-ट्रिप' करता है - जिसे भौगोलिक रूप से अनुमानित किया गया है। 32133 एनएडी 1983 राज्य विमान दक्षिण कैरोलिना (मीटर) है। 26917 एनएडी 1983 यूटीएम क्षेत्र 17 उत्तर है।

C:\Users\melita>inverse 32133
382000 20000
      -83.40806392522212        31.98974518135408
382000 383000
      -83.50098893136905        35.26180827475587
839100 20000
      -78.57184097446545        31.98934439195045
839100 383000
      -78.47814111839074        35.26139222680582

C:\Users\melita>forward 26917
  -83.40806392522212        31.98974518135408
       272490.5730967618        3541832.738731374
  -83.50098893136905        35.26180827475587
       272485.6257057797         3904944.98998655
  -78.57184097446545        31.98934439195045
       729409.4734382738        3541830.781689366
  -78.47814111839074        35.26139222680582
       729414.4926270114        3904946.919009762

C:\Users\melita>inverse 26917
 272490.5730967618        3541832.738731374
      -83.40806392522212        31.98974518135408
  272485.6257057797         3904944.98998655
      -83.50098893136905        35.26180827475587
  729409.4734382738        3541830.781689366
      -78.57184097446545        31.98934439195045
  729414.4926270114        3904946.919009762
      -78.47814111839074        35.26139222680582
^Z

C:\Users\melita>forward 32133
  -83.40806392522212        31.98974518135408
                382000.0                  20000.0
  -83.50098893136905        35.26180827475587
                382000.0                 383000.0
  -78.57184097446545        31.98934439195045
                839100.0        19999.99999999814
  -78.47814111839074        35.26139222680582
                839100.0        382999.9999999981

तो आप देख सकते हैं कि हमारे पास दो बिंदु थे जो 10e-09 पर वापस आए।

आर्कजीआईएस में हैंडलिंग इस तथ्य से जटिल है कि एक स्थानिक संदर्भ है। स्थानिक संदर्भ में समन्वय प्रणाली और कुछ भंडारण और विश्लेषण मूल्य शामिल हैं। डिफ़ॉल्ट रूप से, मीटर का उपयोग करने वाली समन्वय प्रणालियां 0.00 मिली के दसवें हिस्से की सटीकता के साथ संग्रहीत की जाती हैं।

प्रकटीकरण: मैं Esri के लिए काम करते हैं।

Answer 3

मुझे लगता है कि यह एक ऐसा मामला है जहां आपको कुछ परीक्षण बिंदु सुविधाओं के खिलाफ अपने प्रस्तावित वर्कफ़्लो का परीक्षण करने की आवश्यकता होती है, जो कि XY समन्वय वाले फ़ील्ड को जोड़ना आसान है।

अपने प्रारंभिक बिंदुओं के XY मानों की तुलना उन लोगों के साथ करें जिन्हें आपने प्रोजेक्ट किया है / बदल दिया है (हालांकि कई बार), और आपने अंतर को निर्धारित किया होगा।

Comments

इस बात से सहमत। इसके अलावा, ध्यान रखें कि आर्कगिस तालिका दृश्य में डिफ़ॉल्ट रूप से एक डबल डेटाटाइप के 6 दशमलव स्थानों को प्रदर्शित करता है। आप 12 दशमलव स्थानों को तालिका दृश्य में प्रदर्शित करने के लिए फ़ील्ड के गुणों को बदल सकते हैं। भौगोलिक xy मान आम तौर पर 9 दशमलव स्थान या दोगुना सटीक होते हैं।

— क्लेविस