i.landsat.toargrass - क्लाउड में ऑनलाइन

कार्यक्रम:

नाम

i.landsat.toar - वायुमंडल की ऊपरी चमक या परावर्तन और तापमान की गणना करता है
लैंडसैट एमएसएस/टीएम/ईटीएम+/ओएलआई

कीवर्ड

इमेजरी, रेडियोमेट्रिक रूपांतरण, चमक, परावर्तन, चमक तापमान, लैंडसैट,
वायुमंडलीय सुधार

SYNOPSIS

i.landsat.toar
i.landsat.toar --मदद
i.landsat.toar [-आरएनपी] निवेश=बेसनाम उत्पादन=बेसनाम [मेटफ़ाइल=नाम] [सेंसर=स्ट्रिंग]
[तरीका=स्ट्रिंग] [डेटा=yyyy-mm-dd] [सूर्य_ऊंचाई=नाव] [उत्पाद दिनांक=yyyy-mm-dd]
[लाभ=स्ट्रिंग] [प्रतिशत=नाव] [पिक्सेल=पूर्णांक] [रेले=नाव]
[lsatmet=स्ट्रिंग[,स्ट्रिंग,...]] [स्केल=नाव] [--अधिलेखित] [--मदद] [--वाचाल]
[--शांत] [--ui]

झंडे:
-r
सभी बैंड के लिए परावर्तन के बजाय सेंसर चमक पर आउटपुट

-n
इनपुट रैस्टर मानचित्र एक्सटेंशन के रूप में कोड के बजाय बैंड की संख्या का उपयोग करते हैं

-p
आउटपुट मेटाडेटा जानकारी प्रिंट करें

--ओवरराइट
आउटपुट फ़ाइलों को मौजूदा फ़ाइलों को अधिलेखित करने की अनुमति दें

--मदद
प्रिंट उपयोग सारांश

--शब्दशः
वर्बोज़ मॉड्यूल आउटपुट

--शांत
शांत मॉड्यूल आउटपुट

--यूआई
बलपूर्वक GUI संवाद लॉन्च करना

पैरामीटर:
निवेश=बेसनाम [आवश्यक]
इनपुट रैस्टर बैंड का आधार नाम
उदाहरण: 'बी.' बी.1, बी.2, के लिए...

उत्पादन=बेसनाम [आवश्यक]
आउटपुट रैस्टर मानचित्रों के लिए उपसर्ग
उदाहरण: 'बी.टोअर.' B.toar.1, B.toar.2, उत्पन्न करता है...

मेटफ़ाइल=नाम
लैंडसैट मेटाडेटा फ़ाइल का नाम (.met या MTL.txt)

सेंसर=स्ट्रिंग
अंतरिक्ष यान सेंसर
केवल तभी आवश्यक है जब 'मेटफ़ाइल' नहीं दिया गया हो (बुद्धि के लिए अनुशंसित)
विकल्प: एमएसएस1, एमएसएस2, एमएसएस3, एमएसएस4, एमएसएस5, टीएम4, टीएम5, टीएम7, oli8
mss1: लैंडसैट-1 एमएसएस
mss2: लैंडसैट-2 एमएसएस
mss3: लैंडसैट-3 एमएसएस
mss4: लैंडसैट-4 एमएसएस
mss5: लैंडसैट-5 एमएसएस
tm4: लैंडसैट-4 टीएम
tm5: लैंडसैट-5 टीएम
tm7: लैंडसैट-7 ईटीएम+
oli8: लैंडसैट_8 ओएलआई/टीआईआरएस

तरीका=स्ट्रिंग
वायुमंडलीय सुधार विधि
वायुमंडलीय सुधार विधि
विकल्प: अशोधित, dos1, dos2, dos2b, dos3, डॉस4
चूक: uncorrected

डेटा=yyyy-mm-dd
छवि अधिग्रहण तिथि (yyyy-mm-dd)
केवल तभी आवश्यक है जब 'मेटफ़ाइल' नहीं दिया गया हो

सूर्य_ऊंचाई=नाव
सूर्य की ऊंचाई डिग्री में
केवल तभी आवश्यक है जब 'मेटफ़ाइल' नहीं दिया गया हो

उत्पाद दिनांक=yyyy-mm-dd
छवि निर्माण दिनांक (yyyy-mm-dd)
केवल तभी आवश्यक है जब 'मेटफ़ाइल' नहीं दिया गया हो

लाभ=स्ट्रिंग
सभी लैंडसैट ईटीएम+ बैंड का लाभ (एच/एल) (1-5,61,62,7,8)
केवल तभी आवश्यक है जब 'मेटफ़ाइल' नहीं दिया गया हो

प्रतिशत=नाव
पथ चमक में सौर चमक का प्रतिशत
केवल तभी आवश्यक है जब 'विधि' कोई डॉस हो
चूक: 0.01

पिक्सेल=पूर्णांक
डिजिटल नंबर को डार्क ऑब्जेक्ट मानने के लिए न्यूनतम पिक्सेल
केवल तभी आवश्यक है जब 'विधि' कोई डॉस हो
चूक: 1000

रेले=नाव
रेले वातावरण (फैला हुआ आकाश विकिरण)
केवल तभी आवश्यक है जब 'विधि' DOS3 हो
चूक: 0.0

lsatmet=स्ट्रिंग [, स्ट्रिंग,...]
किसी दिए गए मेटाडेटा के लिए संग्रहीत रिटर्न मान
केवल तभी आवश्यक है जब 'मेटफाइल' और -पी दिया गया हो
विकल्प: संख्या, सृजन के, दिनांक, सूर्य_ऊंचाई, सेंसर, बैंड, सुनाज़, पहर
संख्या: लैंडसैट नंबर
निर्माण: निर्माण टाइमस्टैम्प
डेटा: दिनांक
sun_elev: सूर्य उन्नयन
सेंसरसेंसर
बैंड: बैंड गिनती
सुनाज़: सूर्य अज़ीमुथ कोण
पहर: समय

स्केल=नाव
आउटपुट के लिए स्केल फैक्टर
चूक: 1.0

वर्णन

i.landsat.toar लैंडसैट इमेजरी के कैलिब्रेटेड डिजिटल नंबर को बदलने के लिए उपयोग किया जाता है
उत्पादों को वायुमंडल के शीर्ष की चमक या वायुमंडल के शीर्ष परावर्तन और तापमान के अनुरूप
(सेंसर टीएम और ईटीएम+ का बैंड 6)। वैकल्पिक रूप से, इसका उपयोग गणना के लिए किया जा सकता है
वायुमंडलीय सुधार (DOS विधि) के साथ सतह पर चमक या परावर्तन।

आमतौर पर, ऐसा करने के लिए उत्पादन तिथि, अधिग्रहण तिथि और सौर उन्नयन शामिल हैं
आवश्यकता है। इसके अलावा, लैंडसैट-7 ईटीएम+ के लिए नौ का लाभ (उच्च या निम्न) भी आवश्यक है
संबंधित बैंड.

वैकल्पिक रूप से (अनुशंसित), डेटा को मेटाडेटा फ़ाइल (.met या MTL.txt) से पढ़ा जा सकता है
सभी लैंडसैट एमएसएस, टीएम, ईटीएम+ और ओएलआई/टीआईआरएस। हालाँकि, यदि सौर उन्नयन मान दिया गया है
मेटाडेटा फ़ाइल का ओवरराइट किया गया है. यह तब आवश्यक है जब डेटा .met फ़ाइल में हो
ग़लत या सटीक नहीं. इसके अलावा, यदि अधिग्रहण या उत्पादन की तारीखें नहीं मिलती हैं
मेटाडेटा फ़ाइल तब कमांड लाइन मानों का उपयोग किया जाता है।

ध्यान दें: इनपुट रैस्टर में कोई भी शून्य मान या QCALmin से छोटा मान शून्य पर सेट किया गया है
आउटपुट रैस्टर और यह समीकरणों में शामिल नहीं है।

uncorrected at-सेंसर मानों (विधि=अशोधित, चूक)

मानक ज्यामितीय और रेडियोमेट्रिक सुधारों के परिणामस्वरूप एक कैलिब्रेटेड डिजिटल नंबर प्राप्त होता है
(क्यूसीएएल = डीएन) छवियां। रोशनी ज्यामिति के प्रभाव को और अधिक मानकीकृत करने के लिए, QCAL
छवियों को पहले एट-सेंसर चमक में और फिर एट-सेंसर परावर्तन में परिवर्तित किया जाता है।
थर्मल बैंड को पहले QCAL से एट-सेंसर रेडिएंस में और फिर प्रभावी में परिवर्तित किया जाता है
केल्विन डिग्री में सेंसर तापमान।

रेडियोमेट्रिक अंशांकन QCAL को परिवर्तित करता है at-सेंसर चमक, एक रेडियोमेट्रिक मात्रा
समीकरणों का उपयोग करके W/(m² * sr * µm) में मापा जाता है:

· लाभ = (Lmax - Lmin) / (QCALmax - QCALmin)

· पूर्वाग्रह = एलमिन - लाभ * क्यूसीएएलमिन

· चमक = लाभ * QCAL + पूर्वाग्रह
कहा पे, Lmax और lmin अंशांकन स्थिरांक हैं, और QCALmax और QCALmin हैं
QCAL में पुनः मापी गई चमक की सीमा के उच्चतम और निम्नतम बिंदु।

फिर, गणना करने के लिए at-सेंसर परावर्तन समीकरण हैं:

· sun_radiance = [Esun * syn(e)] / (PI * d^2)

· परावर्तन = दीप्ति / सूर्य_प्रभा
कहा पे, d पृथ्वी-सूर्य की दूरी खगोलीय इकाइयों में है, e सौर उन्नयन कोण है,
और Esun औसत सौर बाह्य वायुमंडलीय विकिरण W/(m² * µm) में है।

सरलीकृत सतह पर मानों (विधि=क्या करें[1-4])

वायुमंडलीय सुधार और परावर्तन अंशांकन पथ की चमक को हटा देते हैं, अर्थात
वायुमंडल से भटकती रोशनी, और सौर रोशनी का वर्णक्रमीय प्रभाव। आउटपुट के लिए
ये सरल सतह पर चमक और सतह पर परावर्तन, समीकरण हैं (के लिए नहीं)।
थर्मल बैंड):

· sun_radiance = TAUv * [Esun * syn(e) * TAUz + Esky] / (PI * d^2)

· दीप्ति_पथ = दीप्ति_अंधेरा - प्रतिशत * सूर्य_चमक

· चमक = (at-sensor_radiance - radiance_path)

· परावर्तन = दीप्ति / सूर्य_प्रभा
कहा पे, प्रतिशत 0.0 और 1.0 के बीच का मान है (आमतौर पर 0.01), एस्की फैला हुआ आकाश है
विकिरण, ताउज़ सूर्य से सूर्य तक के मार्ग में वायुमंडलीय संचरण है
ज़मीन की सतह, और ताऊव जमीन से पथ के साथ वायुमंडलीय संचरण है
सेंसर की सतह। दीप्ति_अंधेरा सबसे अंधेरे से गणना की गई एट-सेंसर चमक है
ऑब्जेक्ट, यानी पूरी छवि के लिए कम से कम 'डार्क_पैरामीटर' (आमतौर पर 1000) पिक्सल वाला डीएन।
मूल्य हैं,

· DOS1: TAUv = 1.0, TAUz = 1.0 और Esky = 0.0

· DOS2: अधिकतम तरंग वाले सभी बैंड के लिए TAUv = 1.0, Esky = 0.0, और TAUz = syn(e)
लंबाई 1 से कम। (यानी बैंड 4-6 एमएसएस, 1-4 टीएम, और 1-4 ईटीएम+) अन्य बैंड TAUz =
1.0

· DOS3: TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = rayleigh

· DOS4: TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = PI *
दीप्ति_अंधेरा
ध्यान दें: तब आउटपुट चमक अछूती रहती है (यानी नकारात्मक होने पर 0.0 पर कोई सेट नहीं)।
वे संभावित नकारात्मक मूल्य हैं। हालाँकि, आउटपुट परावर्तन 0.0 पर सेट है
ऋणात्मक मान प्राप्त हुआ।

टिप्पणियाँ

आउटपुट रैस्टर सेल मानों को इसके साथ पुनः स्केल किया जा सकता है स्केल पैरामीटर (उदाहरण के लिए, 100 इंच के साथ
परावर्तन आउटपुट का उपयोग करने का मामला i.gensigset).

On लैंडसैट-8 मेटाडेटा पट्टिका
NASA ने लैंडसैट डेटा के लिए L1G मेटाडेटा फ़ाइल (LDCM-DFCB-004.pdf) की संरचना की रिपोर्ट दी
निरंतरता मिशन (यानी लैंडसैट-8)।

NASA MIN_MAX_RADIANCE समूह में डिजिटल परिवर्तन के लिए आवश्यक जानकारी रखता है
चमक मूल्यों में संख्याएँ (डीएन)। तब, i.landsat.toar संभावित मानक को प्रतिस्थापित करता है
मेटाडेटा मानों के साथ मान। परिणाम मेटाडा द्वारा रिपोर्ट किए गए मानों से मेल खाते हैं
RADIOMETRIC_RESCALING समूह में फ़ाइल।

इसके अलावा, नासा अधिकतम-न्यूनतम मूल्यों और में सभी बैंडों के लिए परावर्तन के समान मूल्यों की रिपोर्ट करता है
लाभ-पूर्वाग्रह मूल्य। यह अजीब है कि सभी बैंडों की परावर्तन सीमा समान होती है। भी,
उन्होंने वेब पेज में लिखा कि पहले डीएन से सीधे परावर्तन की गणना करें
RADIOMETRIC_RESCALING मान और दूसरा पाप (सूर्य_ऊंचाई) से विभाजित।

यह एक साधारण रीस्केलिंग है

· परावर्तन = चमक / सूर्य_चमक = (DN * RADIANCE_MULT + RADIANCE_ADD) /
सूर्य_चमक

· अब परावर्तन = DN * REFLECTANCE_MULT + REFLECTANCE_ADD

· फिर REFLECTANCE_MULT = RADIANCE_MULT / sun_radiance

· और REFLECTANCE_ADD = RADIANCE_ADD / sun_radiance

समस्या तब उत्पन्न होती है जब हमें sun_radiance और की गणना करने के लिए ESUN मानों (प्रदान नहीं किए गए) की आवश्यकता होती है
करने योग्य। हम मानते हैं कि REFLECTANCE_MAXIMUM RADIANCE_MAXIMUM से मेल खाता है, तो

· REFLECTANCE_MAXIMUM / पाप(ई) = RADIANCE_MAXIMUM / sun_radiance

· एसुन = (पीआई * डी^2) * रेडिएंस_मैक्सिमम / रिफ्लेक्टेंस_मैक्सिमम
जहां d मेटाडेटा फ़ाइल द्वारा प्रदान की गई या इसके अंदर गणना की गई पृथ्वी-सूर्य की दूरी है
कार्यक्रम.

RSI i.landsat.toar Lmax, Lmin, और Esun का उपयोग जारी रखने के लिए NASA पुनर्स्केलिंग को वापस लौटाता है
डीएन को चमक और चमक को परावर्तन में बदलने के लिए स्थिरांक की गणना करने के लिए मान
"पारंपरिक" समीकरण और सरल वायुमंडलीय सुधार। ध्यान दें: जब अधिकतम
मान प्रदान नहीं किए गए हैं, i.landsat.toar Lmax, Lmin और Esun की गणना करने का प्रयास करता है
RADIOMETRIC_RESCALING (परीक्षणों में परिणाम समान थे)।

कैलिब्रेशन स्थिर
वर्बोज़ मोड में (ध्वज --शब्दशः), प्रोग्राम बुनियादी उपग्रह डेटा लिखता है और
परिवर्तनों में प्रयुक्त पैरामीटर।

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

· लैंडसैट-1 एमएसएस: कभी नहीं

· लैंडसैट-2 एमएसएस: 16 जुलाई, 1975

· लैंडसैट-3 एमएसएस: 1 जून, 1978

· लैंडसैट-4 एमएसएस: 26 अगस्त, 1982 और 1 अप्रैल, 1983

· लैंडसैट-4 टीएम: 1 अगस्त 1983 और 15 जनवरी 1984

· लैंडसैट-5 एमएसएस: 6 अप्रैल, 1984 और 9 नवंबर, 1984

· लैंडसैट-5 टीएम: 4 मई 2003 और 2 अप्रैल 2007

· लैंडसैट-7 ईटीएम+: 1 जुलाई 2000

· लैंडसैट-8 ओएलआई/टीआईआरएस: 2013 में लॉन्च किया गया

उदाहरण

मेटाडाटा पट्टिका उदाहरण
लैंडसैट-7 ईटीएम+ के डिजिटल नंबरों को बैंड रैस्टर 203_30.1, 203_30.2 [...] में परिवर्तित करें
आउटपुट फ़ाइलों 203_30.1_toar, 203_30.2_toar [...] और में असंशोधित एटी-सेंसर परावर्तन
आउटपुट फ़ाइलों में सेंसर तापमान 293_39.61_toar और 293_39.62_toar:
i.landsat.toar इनपुट=203_30। आउटपुट=_टार \
मेटफ़ाइल=p203r030_7x20010620.met
or
i.landsat.toar इनपुट=L5121060_06020060714। \
आउटपुट=L5121060_06020060714_toar \
मेटफ़ाइल=L5121060_06020060714_MTL.txt
or
i.landsat.toar इनपुट=LC80160352013134LGN03_B आउटपुट=toar \
metfile=LC80160352013134LGN03_MTL.txt sensor=oli8 date=2013-05-14

डॉस1 उदाहरण
DOS1 का उपयोग करके परावर्तन के लिए DN:
# i.landsat.toar की इनपुट योजना से मिलान करने के लिए चैनलों का नाम बदलें या एक प्रतिलिपि बनाएं:
g.copy raster=lsat7_2002_10,lsat7_2002.1
g.copy raster=lsat7_2002_20,lsat7_2002.2
g.copy raster=lsat7_2002_30,lsat7_2002.3
g.copy raster=lsat7_2002_40,lsat7_2002.4
g.copy raster=lsat7_2002_50,lsat7_2002.5
g.copy raster=lsat7_2002_61,lsat7_2002.61
g.copy raster=lsat7_2002_62,lsat7_2002.62
g.copy raster=lsat7_2002_70,lsat7_2002.7
g.copy raster=lsat7_2002_80,lsat7_2002.8
DOS1 का उपयोग करके DN से परावर्तन मानों की गणना (मेटाडेटा से प्राप्त)।
p016r035_7x20020524.met.gz):
i.landsat.toar इनपुट=lsat7_2002. आउटपुट=lsat7_2002_toar. सेंसर=tm7\
विधि=dos1 दिनांक=2002-05-24 रवि_ऊंचाई=64.7730999 \
उत्पाद दिनांक=2004-02-12 लाभ=HHHLHLHHL
परिणामी लैंडसैट चैनलों के नाम lsat7_2002_toar.1 .. lsat7_2002_toar.8 हैं।

संदर्भ

· चंदर जी., बीएल मार्खम और डीएल हेल्डर, 2009: पर्यावरण की रिमोट सेंसिंग,
उड़ान। 113

· चंदर जीएच और बी. मार्खम, 2003: जियोसाइंस और रिमोट पर आईईईई लेनदेन
सेंसिंग, वॉल्यूम. 41, नहीं. 11।

· चावेज़ पीएस, जूनियर। 1996. छवि-आधारित वायुमंडलीय सुधार - पुनरीक्षित और
सुधार हुआ. फोटोग्राममेट्रिक इंजीनियरिंग और रिमोट सेंसिंग 62(9): 1025-1036।

· हुआंग एट अल: एट-सैटेलाइट रिफ्लेक्शन, 2002: ए फर्स्ट ऑर्डर नॉर्मलाइज़ेशन ऑफ़
लैंडसैट 7 ईटीएम+ छवियाँ।

· आर. आयरिश: लैंडसैट 7. विज्ञान डेटा उपयोगकर्ता पुस्तिका। फ़रवरी 17, 2007; 15 मई 2011.

· मार्खम बीएल और जेएल बार्कर, 1986: लैंडसैट एमएसएस और टीएम पोस्ट-कैलिब्रेशन डायनेमिक
रेंज, बाहरी वायुमंडलीय परावर्तन और उपग्रह पर तापमान। ईओएसएटी लैंडसैट
तकनीकी नोट्स, नंबर 1.

· मोरन एमएस, आरडी जैक्सन, पीएन स्लेटर और पीएम टीलेट, 1992: रिमोट सेंसिंग
पर्यावरण, खंड. 41.

· सॉन्ग एट अल, 2001: लैंडसैट टीएम डेटा का उपयोग करके वर्गीकरण और परिवर्तन का पता लगाना, कब
और वायुमंडलीय प्रभावों को कैसे ठीक करें? पर्यावरण का रिमोट सेंसिंग, संस्करण। 75.

onworks.net सेवाओं का उपयोग करके i.landsat.toargrass का ऑनलाइन उपयोग करें