यह grdfiltergmt कमांड है जिसे हमारे कई मुफ्त ऑनलाइन वर्कस्टेशन जैसे कि Ubuntu ऑनलाइन, Fedora ऑनलाइन, Windows ऑनलाइन एमुलेटर या MAC OS ऑनलाइन एमुलेटर का उपयोग करके OnWorks मुफ्त होस्टिंग प्रदाता में चलाया जा सकता है।
कार्यक्रम:
नाम
grdfilter - अंतरिक्ष (या समय) डोमेन में ग्रिड फ़िल्टर करें
SYNOPSIS
grdfilter इंग्रिड दूरी_ध्वज xचौडाई[/चौड़ाई2][संशोधक] आउटग्रिड [ वेतन वृद्धि ] [ i|p|r ]
[ क्षेत्र ] [ ] [[स्तर] ] [ -f]
नोट: विकल्प ध्वज और संबंधित तर्कों के बीच किसी स्थान की अनुमति नहीं है।
वर्णन
grdfilter चयनित कन्वोल्यूशन में से किसी एक का उपयोग करके समय डोमेन में ग्रिड फ़ाइल को फ़िल्टर करेगा
या गैर-कन्वल्यूशन आइसोट्रोपिक या आयताकार फिल्टर और कार्टेशियन का उपयोग करके दूरी की गणना करें
या गोलाकार ज्यामितियाँ। आउटपुट ग्रिड फ़ाइल को वैकल्पिक रूप से उप-क्षेत्र के रूप में भी जनरेट किया जा सकता है
इनपुट के (के माध्यम से) -R) और/या नई वेतन वृद्धि के साथ (के माध्यम से -I) या पंजीकरण (के माध्यम से) -T)। इसमें
वैसे, इनपुट डेटा में "अतिरिक्त स्थान" हो सकता है ताकि किनारों का उपयोग न किया जा सके और
आउटपुट इनपुट किनारों की आधी चौड़ाई के भीतर हो सकता है। यदि फ़िल्टर लो-पास है,
तो आउटपुट को इनपुट की तुलना में कम बार नमूना लिया जा सकता है।
आवश्यक बहस
इंग्रिड फ़िल्टर किए जाने वाले बिंदुओं की ग्रिड फ़ाइल. (नीचे ग्रिड फ़ाइल स्वरूप देखें).
-Dदूरी_ध्वज
दूरी झंडा बताता है कि ग्रिड (x,y) फ़िल्टर से कैसे संबंधित है चौडाई के रूप में इस प्रकार है:
झंडा = p: ग्रिड (px,py) के साथ चौडाई पिक्सेल की विषम संख्या; कार्तीय दूरी.
झंडा = 0: ग्रिड (x,y) के समान इकाइयाँ चौडाई, कार्तीय दूरियाँ।
झंडा = 1: ग्रिड (x,y) डिग्री में, चौडाई किलोमीटर में, कार्तीय दूरियाँ।
झंडा = 2: ग्रिड (x,y) डिग्री में, चौडाई किमी में, dx को cos(मध्य y) द्वारा मापा गया, कार्टेशियन
दूरियां।
उपरोक्त विकल्प सबसे तेज़ हैं क्योंकि वे केवल वज़न मैट्रिक्स की गणना करने की अनुमति देते हैं
एक बार। अगले तीन विकल्प धीमे हैं क्योंकि वे प्रत्येक के लिए वजन की पुनः गणना करते हैं
अक्षांश.
झंडा = 3: ग्रिड (x,y) डिग्री में, चौडाई किमी में, dx को कोसाइन (y), कार्टेशियन द्वारा बढ़ाया गया
दूरी की गणना.
झंडा = 4: ग्रिड (x,y) डिग्री में, चौडाई किमी में, गोलाकार दूरी की गणना।
झंडा = 5: मर्केटर में ग्रिड (x,y) -जे.एम1 छवि इकाइयाँ, चौडाई किमी में, गोलाकार दूरी
गणना।
-एफएक्सचौडाई[/चौड़ाई2][संशोधक]
फ़िल्टर प्रकार सेट करता है। कन्वोल्यूशन और नॉन-कन्वोल्शन फ़िल्टर में से चुनें। कोई भी इस्तेमाल करें
फ़िल्टर कोड x (नीचे सूचीबद्ध) उसके बाद पूरा व्यास चौडाई. इससे एक
आइसोट्रोपिक फ़िल्टर; जोड़ें /चौड़ाई2 एक आयताकार फिल्टर के लिए (आवश्यक) -डीपी or -डी0).
कुछ फ़िल्टर वैकल्पिक तर्कों और संशोधकों की अनुमति देते हैं।
कन्वोल्यूशन फिल्टर (और उनके कोड) हैं:
(b) बॉक्सकार: सभी वजन बराबर हैं।
(c) कोसाइन आर्क: वज़न कोसाइन आर्क वक्र का अनुसरण करता है।
(g) गॉसियन: भार गॉसियन फ़ंक्शन द्वारा दिए जाते हैं, जहाँ चौडाई 6 गुना है
पारंपरिक गौसियन सिग्मा.
(f) कस्टम: भार फ़िल्टर भार ग्रिड में पूर्व-गणना किए गए मानों द्वारा दिए जाते हैं
पट्टिका भार, जिसके आयाम विषम होने चाहिए; इसके लिए भी आवश्यकता होती है -डी0 और आउटपुट रिक्ति
इनपुट स्पेसिंग से मेल खाना चाहिए या पूर्णांक गुणज होना चाहिए।
(o) ऑपरेटर: भार फ़िल्टर भार ग्रिड में पूर्व-गणना किए गए मानों द्वारा दिए जाते हैं
पट्टिका भार, जिसके आयाम विषम होने चाहिए; इसके लिए भी आवश्यकता होती है -डी0 और आउटपुट रिक्ति
इनपुट स्पेसिंग से मेल खाना चाहिए या पूर्णांक गुणज होना चाहिए। भारों का योग माना जाता है
शून्य इसलिए भार योग और सामान्यीकरण का कोई संचय नहीं किया जाएगा।
गैर-कन्वल्यूशन फिल्टर (और उनके कोड) हैं:
(m) माध्यिका: माध्यिका मान लौटाता है। किसी अन्य क्वांटाइल को चुनने के लिए, जोड़ें +qमात्रात्मक in
0-1 रेंज [डिफ़ॉल्ट 0.5 है, यानी, माध्यिका]।
(p) अधिकतम संभावना संभावना (एक मोड अनुमानक): रिटर्न मोडल मान। यदि अधिक
एक से ज़्यादा मोड मिलने पर हम उनका औसत मान लौटाते हैं। +l or +u यदि आप चाहें तो
मोडल मानों में से सबसे निचले या सबसे ऊपरी मान को वापस करना चाहते हैं।
(h) हिस्टोग्राम मोड (एक अन्य मोड अनुमानक): केंद्र के रूप में मोडल मान लौटाएँ
हिस्टोग्राम में प्रमुख शिखर का। / जोड़ेंबिनविड्थ बिनिंग निर्दिष्ट करने के लिए
अंतराल. संशोधक का उपयोग करें +c डिब्बों को के गुणजों पर केन्द्रित करने के लिए बिनविड्थ [डिफ़ॉल्ट है
बिन किनारे जो के गुणज हैं बिनविड्थ] यदि एक से अधिक बहुलक पाए जाते हैं तो हम
उनका औसत मान लौटाएँ। जोड़ें +l or +u यदि आप वापस लौटना चाहते हैं
मोडल मानों में सबसे निचला या सबसे ऊपरी मान।
(l) निचला: सभी मानों में से न्यूनतम लौटाएँ।
(L) निचला: केवल सभी सकारात्मक मानों में से न्यूनतम लौटाएँ।
(u) ऊपरी: सभी मानों में से अधिकतम लौटाएँ।
(U) ऊपरी: केवल अधिकतम या सभी नकारात्मक मान लौटाएँ।
की दशा में L|U यह संभव है कि कोई भी डेटा प्रारंभिक संकेत परीक्षण में उत्तीर्ण न हो; में
उस स्थिति में फ़िल्टर NaN लौटाएगा.
-Gआउटग्रिड
आउटग्रिड फ़िल्टर की आउटपुट ग्रिड फ़ाइल है। (नीचे ग्रिड फ़ाइल स्वरूप देखें)।
वैकल्पिक बहस
-Ixinc[इकाई][=|+][/यिनसी[इकाई][=|+]]
x_inc [और वैकल्पिक रूप से वाई_इंक] आउटपुट इंक्रीमेंट है। संलग्न m चाप को इंगित करने के लिए
मिनट, या s आर्क सेकंड को दर्शाने के लिए। यदि नया x_inc, वाई_इंक पूर्णांक नहीं हैं
पुराने के गुणकों (इनपुट डेटा में), फ़िल्टरिंग काफी हद तक होगी
धीमा. [डिफ़ॉल्ट: इनपुट के समान.]
-नि|पी|आर
निर्धारित करें कि इनपुट ग्रिड में NaN-मान फ़िल्टर किए गए आउटपुट को कैसे प्रभावित करते हैं: जोड़ें i सेवा मेरे
फ़िल्टर किए गए मान की गणना में सभी NaN को अनदेखा करें [डिफ़ॉल्ट], r के रूप में ही है i
सिवाय इसके कि यदि इनपुट नोड NaN था तो आउटपुट नोड NaN पर सेट किया जाएगा (केवल
यदि दोनों ग्रिड सह-पंजीकृत हैं तो लागू होता है), और p जो फ़िल्टर किए गए मान को बाध्य करेगा
यदि फ़िल्टर सर्कल के अंदर NaN-मान वाले कोई ग्रिड-नोड पाए जाते हैं तो NaN होना चाहिए।
-R पश्चिम, पूर्व, दक्षिण, तथा उत्तर आउटपुट बिंदुओं के क्षेत्र को परिभाषित करता है। [गलती करना:
इनपुट के समान।]
-T आउटपुट ग्रिड के लिए नोड पंजीकरण को टॉगल करें ताकि इसके विपरीत हो सके
इनपुट ग्रिड [डिफ़ॉल्ट इनपुट ग्रिड के समान पंजीकरण देता है]।
-वी[स्तर] (अधिक ...)
वर्बोसिटी स्तर [सी] का चयन करें।
-एफ [मैं | ओ]कॉलिनफो (अधिक ...)
डेटा प्रकार के इनपुट और/या आउटपुट कॉलम निर्दिष्ट करें।
-^ or केवल -
कमांड के सिंटैक्स के बारे में एक छोटा संदेश प्रिंट करें, फिर बाहर निकलें (नोट: विंडोज़ पर
बस का उपयोग करें -).
-+ or केवल +
व्यापक उपयोग (सहायता) संदेश प्रिंट करें, जिसमें किसी की व्याख्या भी शामिल है
मॉड्यूल-विशिष्ट विकल्प (लेकिन GMT सामान्य विकल्प नहीं), फिर बाहर निकल जाता है।
-? or नहीं तर्क
विकल्पों की व्याख्या सहित एक पूर्ण उपयोग (सहायता) संदेश प्रिंट करें, फिर
बाहर निकलता है।
--संस्करण
GMT संस्करण प्रिंट करें और बाहर निकलें।
--शो-डेटादिर
GMT शेयर निर्देशिका के लिए पूरा पथ प्रिंट करें और बाहर निकलें।
ग्रिड फ़ाइल प्रारूप
डिफ़ॉल्ट रूप से GMT ग्रिड को लिखता है क्योंकि COARDS-शिकायत netCDF में एकल परिशुद्धता तैरती है
फाइल प्रारूप। हालाँकि, GMT कई अन्य सामान्य रूप से उपयोग किए जाने वाले ग्रिड में ग्रिड फ़ाइलों का उत्पादन करने में सक्षम है
फ़ाइल स्वरूपों और ग्रिड के तथाकथित "पैकिंग" की सुविधा भी देता है, फ्लोटिंग पॉइंट को लिखता है
डेटा 1- या 2-बाइट पूर्णांक के रूप में। सटीक, स्केल और ऑफ़सेट निर्दिष्ट करने के लिए, उपयोगकर्ता को चाहिए
प्रत्यय जोड़ें =id[/स्केल/ओफ़्सेट[/नान]], कहां id ग्रिड का दो अक्षर का पहचानकर्ता है
प्रकार और सटीक, और स्केल और ओफ़्सेट वैकल्पिक पैमाने के कारक हैं और होने के लिए ऑफसेट हैं
सभी ग्रिड मानों पर लागू होता है, और नान लापता डेटा को इंगित करने के लिए उपयोग किया जाने वाला मान है। यदि
दो अक्षर id प्रदान नहीं किया गया है, जैसा कि =/स्केल एक से id=nf माना गया है। कब
ग्रिड पढ़ना, प्रारूप आमतौर पर स्वचालित रूप से पहचाना जाता है। यदि नहीं, तो वही प्रत्यय
इनपुट ग्रिड फ़ाइल नामों में जोड़ा जा सकता है। देखो जीआरडी कन्वर्ट और अनुभाग ग्रिड-फ़ाइल-प्रारूप
अधिक जानकारी के लिए GMT तकनीकी संदर्भ और कुकबुक।
एक से अधिक ग्रिड वाली netCDF फ़ाइल को पढ़ते समय, GMT डिफ़ॉल्ट रूप से पढ़ेगा:
पहला 2-आयामी ग्रिड जो उस फ़ाइल में मिल सकता है। GMT को दूसरे पढ़ने के लिए मनाना
ग्रिड फ़ाइल में बहु-आयामी चर, संलग्न करें ?VARNAME फ़ाइल नाम के लिए, जहां
VARNAME चर का नाम है। ध्यान दें कि आपको विशेष अर्थ से बचने की आवश्यकता हो सकती है
of ? अपने शेल प्रोग्राम में इसके सामने बैकस्लैश लगाकर, या रखकर
फ़ाइल नाम और प्रत्यय उद्धरण या दोहरे उद्धरण चिह्नों के बीच। NS ?VARNAME प्रत्यय का भी प्रयोग किया जा सकता है
आउटपुट ग्रिड के लिए डिफ़ॉल्ट से भिन्न चर नाम निर्दिष्ट करने के लिए: "z"। देखो
जीआरडी कन्वर्ट और अनुभाग संशोधक-के लिए-CF और GMT तकनीकी के ग्रिड-फ़ाइल-प्रारूप
अधिक जानकारी के लिए संदर्भ और कुकबुक, विशेष रूप से 3- के स्प्लिसेस को कैसे पढ़ें,
4-, या 5-आयामी ग्रिड।
भौगोलिक और TIME COORDINATES
जब आउटपुट ग्रिड प्रकार netCDF होता है, तो निर्देशांक "देशांतर" लेबल किए जाएंगे,
इनपुट डेटा या ग्रिड (यदि कोई हो) की विशेषताओं के आधार पर "अक्षांश", या "समय" या
-f or -R विकल्प। उदाहरण के लिए, दोनों -f0x -f1t और -R90w/90e/0t/3t का परिणाम a . होगा
देशांतर/समय ग्रिड। जब x, y, या z निर्देशांक समय होता है, तो इसे ग्रिड में संग्रहीत किया जाएगा
TIME_UNIT और TIME_EPOCH द्वारा निर्दिष्ट युग के बाद से सापेक्ष समय के रूप में जीएमटी.conf पट्टिका
या कमांड लाइन पर। इसके साथ में इकाई समय चर की विशेषता इंगित करेगी
इस इकाई और युग दोनों।
उदाहरण
मान लीजिए कि north_pacific_etopo5.nc 140E से 260E तक 5 मिनट की बाथिमेट्री की फ़ाइल है
और 0N से 50N तक, और आप 300 किमी त्रिज्या (600 किमी) के भीतर मानों के मध्यिका ज्ञात करना चाहते हैं
आउटपुट बिंदुओं की पूर्ण चौड़ाई), जिसे आप 150E से 250E और 10N से 40N तक चुनते हैं,
और आप हर 0.5 डिग्री पर आउटपुट मान चाहते हैं। गोलाकार दूरी की गणना का उपयोग करते हुए,
आप की जरूरत है:
gmt grdfilter north_pacific_etopo5.nc -Gfiltered_pacific.nc -Fm600 \
-डी4 -आर150/250/10/40 -आई0.5 -वी
यदि हम उच्च-पास परिणाम चाहते हैं तो हम संगत निम्न-पास परिणाम कर सकते हैं
मोटे ग्रिड अंतराल का उपयोग करके फ़िल्टर करें क्योंकि grdfilter परिणाम को उसी पर पुनः नमूना करेगा
रिज़ॉल्यूशन को इनपुट ग्रिड के रूप में उपयोग करें ताकि हम अवशिष्टों की गणना कर सकें, उदाहरण के लिए,
gmt grdfilter north_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600 \
-डी4 -आर150/250/10/40 -आई0.5 -वी
यहां, residual_pacific.nc ग्रिड का 5 मिनट का रिज़ॉल्यूशन मूल ग्रिड के समान ही होगा।
कस्टम अनिसोट्रोपिक गॉसियन फ़िल्टर exp का उपयोग करके ripples.nc में डेटासेट फ़िल्टर करने के लिए
(-0.5*r^2) जिसकी केंद्र से दूरी r, (2x^2 + y^2 -2xy)/6 द्वारा दी गई है, प्रमुख के साथ
अक्ष को क्षैतिज से 63 डिग्री के कोण पर रखें, प्रयास करें
gmt grdmath -R-10/10/-10/10 -I1 X 2 POW 2 MUL Y 2 POW ADD XY MUL 2 MUL \
सब 6 डिव नेग 2 डिव एक्सपी डुप योग डिव = gfilter.nc
gmt grdfilter ripples.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V
सीमाएं
का उपयोग करने के लिए -D5 विकल्प में इनपुट मर्केटर ग्रिड को img2mercgrd का उपयोग करके बनाया जाना चाहिए -C
विकल्प इसलिए y-मानों का मूल भूमध्य रेखा है (अर्थात, x = y = 0 lon = के अनुरूप है)
अक्षांश = 0).
onworks.net सेवाओं का उपयोग करके grdfiltergmt का ऑनलाइन उपयोग करें
