dimfiltergmt - آنلاین در ابر

برنامه:

نام

dimfilter - فیلتر جهت فایل های شبکه ای دوبعدی در حوزه فضا (یا زمان).

خلاصه

دیمفیلتر input_file.nc فاصله_پرچم [حالت] output_file.nc
[ یقه ] [ افزایش ] [ منطقه ] [ ] [سطح] ] [ -f]

توجه داشته باشید: هیچ فاصله ای بین پرچم گزینه و آرگومان های مرتبط مجاز نیست.

شرح

دیمفیلتر a را فیلتر خواهد کرد .nc فایل در دامنه فضا (یا زمان) با تقسیم داده شده
دایره را فیلتر کنید n_بخش، با استفاده از یکی از پیچیدگی اولیه انتخاب شده یا
فیلترهای غیرپیچیده برای هر بخش، و انتخاب نتیجه نهایی با توجه به
فیلتر ثانویه انتخاب شده فواصل را با استفاده از هندسه های دکارتی یا کروی محاسبه می کند.
خروجی .nc فایل به صورت اختیاری می تواند به عنوان یک زیرمنطقه از ورودی و/یا با a تولید شود
جدید -Iافزایش می یابد. به این ترتیب، ممکن است یک "فضای اضافی" در داده های ورودی داشته باشد تا آنجا
هیچ اثر لبه ای برای شبکه خروجی نخواهد داشت. اگر فیلتر پایین گذر است، پس خروجی
ممکن است کمتر از ورودی نمونه برداری شود. -Q فقط برای حالت تحلیل خطا است
به تعداد کل ستون ها در فایل ورودی نیاز دارد که حاوی فیلتر شده است
اعماق بالاخره آدم باید این را بداند دیمفیلتر خروجی صاف مانند سایرین تولید نخواهد کرد
فیلترهای فضایی این کار را انجام می دهند زیرا حداقل میانه از را برمی گرداند N میانه های N بخش ها در
خروجی می تواند خشن باشد مگر اینکه داده های ورودی بدون نویز باشند. بنابراین، یک فیلتر اضافی
(به عنوان مثال، گاوسی از طریق فیلتر grd) از داده های فیلتر شده با DiM به طور کلی توصیه می شود.

مورد نیاز ادله

input_file.nc
شبکه داده ای که باید فیلتر شود.

-Dفاصله_پرچم
فاصله پرچم نحوه ارتباط شبکه (x,y) با فیلتر را نشان می دهد عرض، به شرح زیر است:

پرچم = 0: شبکه (x,y) واحدهای مشابه عرض، فواصل دکارتی. پرچم = 1: شبکه
(x,y) بر حسب درجه عرض در کیلومتر، فواصل دکارتی. پرچم = 2: شبکه (x,y)
در درجه، عرض در کیلومتر، dx با cos (وسط y)، فواصل دکارتی مقیاس شده است.

گزینه های بالا سریع ترین هستند زیرا آنها اجازه می دهند ماتریس وزن فقط محاسبه شود
یک بار. سه گزینه بعدی کندتر هستند زیرا وزن‌ها را برای هر کدام دوباره محاسبه می‌کنند
عرض جغرافیایی.

پرچم = 3: شبکه (x,y) بر حسب درجه، عرض در کیلومتر، dx مقیاس بندی شده توسط کسینوس (y)، دکارتی
محاسبه فاصله

پرچم = 4: شبکه (x,y) بر حسب درجه، عرض در کیلومتر، محاسبه فاصله کروی.

-F[حالت]
نوع فیلتر اولیه را تنظیم می کند. از بین فیلترهای کانولوشن و غیر پیچیدگی انتخاب کنید.
کد فیلتر و سپس قطر کامل را اضافه کنید عرض. پیچیدگی موجود
فیلترها عبارتند از:

(b) باکس گاری: همه وزن ها برابرند.

(c) قوس کسینوس: وزن ها از منحنی قوس کسینوس پیروی می کنند.

(g) گاوسی: وزن ها با تابع گاوسی داده می شوند.

فیلترهای غیر کانولوشن عبارتند از:

(m) Median: مقدار میانه را برمی گرداند.

(p) حداکثر احتمال احتمال (یک برآوردگر حالت): مقدار معینی برگردانده می شود. اگر بیشتر
بیش از یک حالت یافت می شود که مقدار متوسط ​​آنها را برمی گردانیم. - یا + را به فیلتر اضافه کنید
عرض اگر می خواهید کوچکترین یا بزرگ ترین مقادیر مدال را برگردانید.

-N
نوع فیلتر ثانویه و تعداد بخش های پاپیون را تنظیم می کند. n_بخش باید
عدد صحیح و بزرگتر از 0 باشد n_بخش روی 1 تنظیم شده است، فیلتر ثانویه است
غیر تاثیرگذار. فیلترهای ثانویه موجود عبارتند از:

(l) Lower: حداقل تمام مقادیر فیلتر شده را برمی گرداند.

(u) Upper: حداکثر تمام مقادیر فیلتر شده را برمی گرداند.

(a) میانگین: میانگین تمام مقادیر فیلتر شده را برگردانید.

(m) میانه: میانه تمام مقادیر فیلتر شده را برمی گرداند.

(p) Mode: حالت تمام مقادیر فیلتر شده را برمی گرداند.

-Goutput_file.nc
output_file.nc خروجی فیلتر است.

اختیاری ادله

-I x_inc [و به صورت اختیاری y_inc] افزایش خروجی است. ضمیمه m برای نشان دادن دقیقه،
or c برای نشان دادن ثانیه ها اگر جدید x_inc, y_inc مضرب اعداد صحیح نیستند
موارد قدیمی (در داده های ورودی)، فیلتر کردن به طور قابل توجهی کندتر خواهد بود. [پیش‌فرض: همان
به عنوان ورودی.]

-R غرب, شرق, جنوبو شمال منطقه نقاط خروجی را تعریف می کند. [پیش‌فرض:
مشابه ورودی.]

-T ثبت گره را برای شبکه خروجی تغییر دهید تا برعکس شود
شبکه ورودی [پیش‌فرض مشابه شبکه ورودی را ثبت می‌کند].

-Qیقه یقه تعداد کل ستون ها در فایل جدول متنی ورودی است. برای این حالت،
انتظار دارد اعماق متشکل از چندین ستون را بخواند. هر ستون نشان دهنده یک است
شبکه فیلتر شده با عرض فیلتر، که می تواند توسط grd2xyz -Z. نتیجه
میانه، MAD و متوسط ​​خواهد بود. بنابراین، ستون با میانه ها برای تولید استفاده می شود
جزء منطقه ای و ستون با MAD ها برای انجام خطا استفاده می شود
تجزیه و تحلیل.

-V[سطح] (بیشتر ...)
سطح پرحرفی [c] را انتخاب کنید.

-f[i|o]colinfo (بیشتر ...)
انواع داده های ستون های ورودی و/یا خروجی را مشخص کنید.

-^ or تنها -
یک پیام کوتاه در مورد نحو دستور چاپ کنید، سپس از آن خارج شوید (توجه: در ویندوز
فقط استفاده کنید -).

-+ or تنها +
یک پیام کاربردی (راهنما) گسترده، از جمله توضیح هر یک را چاپ کنید
گزینه مخصوص ماژول (اما نه گزینه های رایج GMT)، سپس خارج می شود.

-? or نه استدلال
سپس یک پیام استفاده کامل (راهنما)، از جمله توضیح گزینه ها را چاپ کنید
خارج می شود.

- نسخه
نسخه GMT را چاپ کنید و خارج شوید.

--show-datadir
مسیر کامل را به دایرکتوری اشتراک گذاری GMT ​​چاپ کنید و از آن خارج شوید.

GRID فایل فرمها

به‌طور پیش‌فرض GMT شبکه را به‌عنوان شناور دقیق تکی در یک netCDF شکایت COARDS می‌نویسد.
فرمت فایل. با این حال، GMT می‌تواند فایل‌های شبکه‌ای را در بسیاری از شبکه‌های رایج دیگر تولید کند
فرمت های فایل و همچنین به اصطلاح "بسته بندی" شبکه ها، نوشتن نقطه شناور را تسهیل می کند
داده ها به صورت اعداد صحیح 1 یا 2 بایتی. برای تعیین دقت، مقیاس و افست، کاربر باید
پسوند را اضافه کنید =id[/مقیاس/چاپ افست[/نان]]، جایی که id یک شناسه دو حرفی از شبکه است
نوع و دقت و مقیاس و چاپ افست ضریب مقیاس اختیاری و افست هستند
به همه مقادیر شبکه اعمال می شود، و نان مقداری است که برای نشان دادن داده های از دست رفته استفاده می شود. در صورت
دو شخصیت id ارائه نشده است، همانطور که در =/مقیاس از یک id=nf مفروض است. چه زمانی
شبکه های خواندن، قالب به طور کلی به طور خودکار تشخیص داده می شود. اگر نه، همان پسوند
را می توان به نام فایل های شبکه ورودی اضافه کرد. دیدن grdconvert و بخش grid-file-format of the
مرجع فنی GMT ​​و کتاب آشپزی برای اطلاعات بیشتر.

هنگام خواندن یک فایل netCDF که حاوی چندین شبکه است، GMT به طور پیش‌فرض، آن را می‌خواند
اولین شبکه 2 بعدی که می توانید در آن فایل پیدا کنید. برای ترغیب GMT به خواندن دیگری
متغیر چند بعدی در فایل شبکه، پیوست ?نام متغیر به نام فایل، جایی که
نام متغیر نام متغیر است. توجه داشته باشید که ممکن است لازم باشد از معنای خاص فرار کنید
of ? در برنامه پوسته خود با قرار دادن یک بک اسلش در مقابل آن یا با قرار دادن علامت
نام فایل و پسوند بین نقل قول یا دو نقل قول. در ?نام متغیر پسوند نیز می تواند استفاده شود
برای شبکه های خروجی برای تعیین نام متغیر متفاوت از پیش فرض: "z". دیدن
grdconvert و بخش‌های اصلاح‌کننده-برای-CF و قالب-فایل شبکه ای GMT ​​فنی
مرجع و کتاب آشپزی برای اطلاعات بیشتر، به ویژه در مورد نحوه خواندن پیوندهای 3-،
شبکه های 4 یا 5 بعدی.

جغرافیایی و TIME مختصات

هنگامی که نوع شبکه خروجی netCDF باشد، مختصات برچسب "طول جغرافیایی" می شوند.
"طول جغرافیایی" یا "زمان" بر اساس ویژگی های داده ورودی یا شبکه (در صورت وجود) یا بر اساس
-f or -R گزینه ها. مثلا هر دو -f0x -f1t و -R90w/90e/0t/3t منجر به a
شبکه طول/زمان زمانی که مختصات x، y یا z زمان باشد، در شبکه ذخیره می شود
به عنوان زمان نسبی از دوران که توسط TIME_UNIT و TIME_EPOCH در تعیین شده است gmt.conf پرونده
یا در خط فرمان علاوه بر این واحد ویژگی متغیر زمان نشان خواهد داد
هم این واحد و هم دوره

مثال ها

فرض کنید north_pacific_dbdb5.nc یک فایل 5 دقیقه ای عمق سنجی از 140E تا 260E است و
0 تا 50 شمالی، و می‌خواهید میانه‌های مقادیر را در شعاع 300 کیلومتری (600 کیلومتر کامل) پیدا کنید
عرض) از نقاط خروجی، که انتخاب می کنید از 150E تا 250E و 10N تا 40N باشند، و
شما مقادیر خروجی را هر 0.5 درجه می خواهید. برای جلوگیری از سوگیری رسانه ها توسط
در صفحه شیب دار، می خواهید دایره فیلتر را به 6 بخش تقسیم کرده و آن را انتخاب کنید
کمترین مقدار در بین 6 میانه با استفاده از محاسبات فاصله کروی، شما نیاز دارید:

دیم فیلتر gmt north_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150/250/10/40 -I0.5 -V

فرض کنید cape_verde.nc یک فایل 0.5 دقیقه عمق سنجی از 32W تا 15W و 8N به
25N، و می خواهید ویژگی های در مقیاس کوچک را حذف کنید تا یک تورم در یک
منطقه از 27.5 وات تا 20.5 وات و 12.5 نیوتن تا 19.5 نیوتن، و شما مقدار خروجی را هر یک می خواهید
2 دقیقه با استفاده از محاسبات فاصله دکارتی، شما نیاز دارید:

دیم فیلتر gmt cape_verde.nc -Gt.nc -Fm220 -Nl8 -D2 -R-27.5/-20.5/12.5/19.5 -I2m -V
gmt grdfilter t.nc -Gcape_swell.nc -Fg50 -D2 -V

فرض کنید که طیف وسیعی از عرض فیلتر را برای یک منطقه مشخص پیدا کرده اید و آن را فیلتر کرده اید
داده های عمق سنجی با استفاده از محدوده عرض فیلتر (مثلاً f100.nc f110.nc f120.nc
f130.nc، و می خواهید یک روند منطقه ای را با استفاده از محدوده عرض فیلتر و
می‌خواهید در هر نقطه داده، تخمین‌های میانگین انحراف مطلق (MAD) را بدست آورید. سپس، شما
باید انجام شود:

gmt grd2xyz f100.nc -Z > f100.d
gmt grd2xyz f110.nc -Z > f110.d
gmt grd2xyz f120.nc -Z > f120.d
gmt grd2xyz f130.nc -Z > f130.d
چسباندن f100.d f110.d f120.d f130.d > depths.d
gmt dimfilter depths.d -Q4 > output.z

محدودیت ها

هنگام کار با شبکه های جغرافیایی (lat، lon)، هر سه فیلتر پیچشی (boxcar،
قوس کسینوس و گاوسی) به درستی وزن فیلتر را برای تغییر در نرمال می کنند.
اندازه gridbox با عرض جغرافیایی، و به درستی تعیین کنید که کدام گره ها برای آن مورد نیاز هستند
پیچیدگی زمانی که فیلتر "دایره" از مرز تناوبی (0-360) عبور کند یا حاوی یک
قطب جغرافیایی با این حال، فیلترهای فضایی، مانند فیلترهای مدیان و مد، استفاده نمی کنند
وزن ها و بنابراین فقط باید در شبکه های دکارتی (یا در عرض های جغرافیایی بسیار کم) استفاده شود.
اگر می‌خواهید چنین فیلترهای فضایی را اعمال کنید، باید داده‌های خود را به یک ناحیه مساوی بفرستید
طرح ریزی و دیم فیلتر را روی شبکه دکارتی حاصل اجرا کنید.

اسکریپت قالب

dim.template.sh یک اسکریپت پوسته اسکلت است که می تواند برای راه اندازی یک DiM کامل استفاده شود.
تجزیه و تحلیل، از جمله تجزیه و تحلیل MAD.

منابع

Kim، S.-S.، و Wessel، P. (2008)، فیلتر میانه جهت دار برای منطقه ای-باقیمانده
جداسازی حمام سنجی، ژئوشیمی. ژئوفیز. Geosyst., 9، Q03005،
doi: 10.1029/2007GC001850.

از dimfiltergmt به صورت آنلاین با استفاده از خدمات onworks.net استفاده کنید