این دستور grdtrackgmt است که می تواند در ارائه دهنده هاست رایگان OnWorks با استفاده از یکی از چندین ایستگاه کاری آنلاین رایگان ما مانند Ubuntu Online، Fedora Online، شبیه ساز آنلاین ویندوز یا شبیه ساز آنلاین MAC OS اجرا شود.
برنامه:
نام
grdtrack - شبکه های نمونه در مکان های مشخص شده (x,y).
خلاصه
grdtrack [ xyfile ] grd1 grd2 ... [ f|p|m|r|R[+l] ] [ طول[u]/ds[فاصله][+a] ] [dfile
] [ خط[,خط،...][+aaz][+d][+iوارز[u]][+lطول[u]][+nnp][+oaz][+rشعاع[u] ] [ ]
منطقه ] [ روش/اصلاح کننده ها ] [ [شعاع[u]][+e|p]] [سطح] ] [ ] -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -n] [
-o] [ -s] [ -:[i|o] ]
توجه داشته باشید: هیچ فاصله ای بین پرچم گزینه و آرگومان های مرتبط مجاز نیست.
شرح
grdtrack یک یا چند فایل شبکه (یا یک فایل Sandwell/Smith IMG) و یک جدول (از
فایل یا ورودی استاندارد؛ اما ببین -E برای استثنا) با موقعیت های (x,y) [یا (lon,lat)] در
دو ستون اول (ممکن است ستون های بیشتری وجود داشته باشد). این شبکه(های) را در قسمت درون یابی می کند
در جدول قرار می گیرد و جدول را با مقادیر درون یابی اضافه شده به عنوان (یک
یا بیشتر) ستون های جدید. متناوبا، از سوی دیگر (-C، ورودی به صورت خط-قطعات و
ما پروفایل های متقابل متعامد را در هر نقطه داده یا با جدایی مساوی ایجاد می کنیم
و از شبکه(های) در امتداد این پروفایل ها نمونه برداری کنید. یک دو مکعبی [پیشفرض]، دوخطی، B-spline یا
نزدیکترین همسایه (نگاه کنید به -n) درون یابی استفاده می شود که به شرایط مرزی نیاز دارد
محدودیت های منطقه (نگاه کنید به -n; پیش فرض از شرایط "طبیعی" استفاده می کند (مشتق جزئی دوم
نرمال به لبه صفر است) مگر اینکه شبکه به طور خودکار به عنوان دوره ای شناخته شود.)
مورد نیاز ادله
-Gفایل شبکه
grdfile یک فایل شبکه باینری دو بعدی با تابع f(x,y) است. اگر شبکه مشخص شده
در قالب Sandwell/Smith Mercator است، باید فهرستی از کاما جدا شده اضافه کنید
آرگومان هایی که شامل یک مقیاس برای ضرب داده ها (معمولاً 1 یا 0.1)، حالت است
که مخفف عبارت زیر است: (0) فایل های Img بدون کد محدودیت، داده ها را برمی گرداند
در همه نقاط، (1) فایل Img با محدودیت ها کدگذاری شده، داده ها را در همه نقاط برمی گرداند، (2)
فایل img با محدودیتهای کدگذاری شده، دادهها را فقط در نقاط محدود و NaN برمیگرداند
در جاهای دیگر، و (3) فایل Img با محدودیتهای کدگذاری شده، 1 در محدودیت و 0 برمیگرداند
در جای دیگر، و به صورت اختیاری حداکثر عرض جغرافیایی در فایل IMG [80.738]. ممکن است تکرار کنید
-G هر چند بار که شبکه هایی دارید که می خواهید نمونه برداری کنید. روش دیگر، استفاده کنید -G+lفهرست
برای ارسال لیستی از نام فایل ها شبکه ها نمونه برداری می شوند و نتایج در خروجی قرار می گیرند
دستور داده شده (فرمت فایل GRID را در زیر ببینید.)
اختیاری ادله
xyfile این یک ASCII است (یا باینری، ببینید -بی) فایلی که در آن 2 ستون اول (x,y) را نگه می دارند.
موقعیت هایی که کاربر می خواهد از مجموعه داده های 2 بعدی نمونه برداری کند.
-Af|pm|r|R
برای نمونه برداری مجدد آهنگ (اگر -C تنظیم شده است) می توانیم نحوه انجام این کار را انتخاب کنیم.
ضمیمه کردن f برای حفظ نقاط اصلی، اما در صورت نیاز نقاط میانی اضافه کنید [پیشفرض]،
m as f، اما ابتدا نصف النهار (در امتداد y) سپس موازی (در امتداد x) را دنبال کنید. p as f، اما
ابتدا موازی (در امتداد y) و سپس نصف النهار (در امتداد x) را دنبال کنید. r نمونه برداری مجدد در
مکان های مساوی نقاط ورودی لزوماً در خروجی گنجانده نمی شوند و
R as r، اما فاصله داده شده را طوری تنظیم کنید که دقیقاً با طول مسیر مطابقت داشته باشد. در نهایت، ضمیمه کنید
+l اگر فاصله ها باید در امتداد خطوط رومب (لوکسودروم) اندازه گیری شود. نادیده گرفته می شود مگر اینکه
-C استفاده می شود.
-Cطول[u]/ds[فاصله][+a]
از قطعات خط ورودی برای ایجاد فاصله مساوی و (اختیاری) با فاصله یکسان استفاده کنید
مجموعه ای از پروفایل های متقاطع که در امتداد آنها از شبکه(های) نمونه برداری می کنیم [پیش فرض به سادگی نمونه است
شبکه(های) در مکان های ورودی]. دو مقیاس طولی را مشخص کنید که چگونه را کنترل می کند
نمونه برداری انجام می شود: طول طول کامل هر نمایه متقاطع را تنظیم می کند، در حالی که ds
فاصله نمونه برداری در امتداد هر نمایه متقابل است. به صورت اختیاری، ضمیمه کنید /فاصله برای
فاصله مساوی بین نمایه های متقابل [پیش فرض نمایه های متقاطع را در
مختصات ورودی]. به طور پیش فرض، تمام پروفایل های متقابل جهت یکسانی دارند. ضمیمه
+a برای تغییر جهت پروفایل های متقابل. ضمیمه واحدهای مناسب به طول;
واحد مورد استفاده را تنظیم می کند ds [و فاصله] (به واحدهای زیر مراجعه کنید). واحد پیش فرض برای
شبکه های جغرافیایی متر است در حالی که شبکه های دکارتی به معنای واحد کاربر است. خروجی
ستون ها خواهد بود خاله, LAT, ناحیه, قوس افقی در جهت گردش عقربه ساعت واقع بین نقطه ثابتی, z1, z2، ... ، zn ( zi نمونه برداری شده اند
مقادیر برای هر یک از n شبکه ها)
-Ddfile
در کنسرت با -C ما می توانیم خطوط اصلی (احتمالاً نمونه برداری مجدد) را در فایل ذخیره کنیم
dfile [پیشفرض فقط نمایههای متقابل را ذخیره میکند]. ستون ها خواهند بود خاله, LAT, ناحیه,
قوس افقی در جهت گردش عقربه ساعت واقع بین نقطه ثابتی, z1, z2، ... (مقدار نمونه برای هر شبکه)
-Eخط[,خط،...][+aaz][+d][+iوارز[u]][+lطول[u]][+nnp][+oaz][+rشعاع[تو]
به جای خواندن مختصات آهنگ ورودی، پروفایل ها را از طریق مختصات و مشخص کنید
اصلاح کننده ها فرمت هر کدام خط is شروع/متوقف کردن، که در آن شروع or متوقف کردن یا یکی هستند
خاله/LAT (x/y برای داده های دکارتی) یا یک کلید XY 2 کاراکتری که از آن استفاده می کند
فرمت فرمت توجیهی به سبک "pstext" برای تعیین یک نقطه روی نقشه به عنوان
[LCR][BMT]. علاوه بر این، می توانید از Z-، Z+ به معنای حداقل و حداکثر جهانی استفاده کنید
مکانها در شبکه (فقط در صورتی موجود است که فقط یک شبکه داده شود). به جای دوتا
مختصات می توانید یک مبدا و یکی از آنها را مشخص کنید +a, +o، یا +r. می توانید اضافه کنید
+iوارز[u] برای تنظیم فاصله نمونه برداری؛ اگر داده نشد، به طور پیش فرض به نصف می رسیم
حداقل فاصله شبکه در +a آزیموت پروفیل با طول مشخص را تعیین می کند
در حالی که از مبدا داده شده شروع می شود +o نمایه را روی مبدا متمرکز می کند. هر دو
نیاز +l. برای نمونه برداری دایره ای مشخص کنید +r برای تعریف دایره ای با شعاع داده شده
با محوریت مبدا؛ این گزینه به هر دو نیاز دارد +n or +i. +nnp تنظیم می کند
تعداد مورد نظر امتیاز، در حالی که +lطول طول کل پروفیل را نشان می دهد.
ضمیمه کردن +d برای خروجی فواصل در طول مسیر پس از مختصات. توجه: خیر
فایل آهنگ خوانده خواهد شد همچنین توجه داشته باشید که فقط یک واحد فاصله قابل انتخاب است.
دادن واحدهای مختلف منجر به خطا می شود. اگر هیچ واحدی مشخص نشده باشد ما
پیشفرض فاصله دایرهای بزرگ بر حسب کیلومتر (در صورت جغرافیایی). اگر با
داده های جغرافیایی که می توانید به آنها اضافه کنید - (زمین مسطح) یا + (ژئودزیک). وارز, طول، یا
شعاع برای تغییر حالت محاسبه فاصله [دایره بزرگ].
-N Do نه رد شدن از نقاطی که خارج از دامنه شبکه(های) قرار دارند [فقط خروجی پیش فرض
نقاط در دامنه شبکه].
-R[واحد]xmin/xmax/ymin/ymax[r] (بیشتر ...)
منطقه مورد نظر را مشخص کنید.
-Sروش/اصلاح کننده ها
در رابطه با -C، یک نمایه انباشته شده را از تمام نمایه ها در سراسر محاسبه کنید
هر بخش نحوه محاسبه انباشتگی را اضافه کنید: a = میانگین (متوسط) m =
متوسط ، p = حالت (حداکثر احتمال)، l = پایین تر، L = پایین تر اما فقط در نظر بگیرید
ارزش های مثبت، u = بالا، U = بالا اما فقط مقادیر منفی را در نظر بگیرید [a]
اصلاح کننده ها کنترل خروجی؛ یکی یا چند مورد را از میان این گزینه ها انتخاب کنید: +a : ضمیمه
مقادیر انباشته شده در تمام نمایه های متقابل. +d : انحرافات پشته را به همه اضافه کنید
پروفایل های متقابل +d : باقیمانده های داده (داده - پشته) را به همه پروفایل های متقابل اضافه کنید.
+s[پرونده] : نمایه انباشته شده را در آن ذخیره کنید پرونده [grdtrack_stacked_profile.txt]. +cواقعیت :
پاکت نامه را در نمایه انباشته به صورت +/- محاسبه کنید واقعیت *انحراف [2]. یادداشت ها: (1)
انحرافات بستگی دارد روش و st.dev هستند (aمقیاس L1 (m و p) یا نیم برد
(بالا-پایین)/2. (2) فایل نمایه انباشته شامل 1 گروه به اضافه 4 تا 6 است
ستون ها، یک گروه برای هر شبکه نمونه برداری شده. ستون اول فاصله متقاطع را نگه می دارد،
در حالی که 4 مورد اول در یک گروه مقدار انباشته، انحراف، حداقل مقدار و حداکثر را نگه می دارند
مقدار. اگر روش یکی از a|m|p سپس اطمینان پایین و بالایی را نیز می نویسیم
مرزها (نگاه کنید به +c). زمانی که یک یا چند مورد از +a, +dو +r استفاده می شود سپس ما ضمیمه
برای تمام نمایه های متقابل به انتهای هر سطر می رسد. سفارش همیشه انباشته است
مقدار (+a، به دنبال آن انحرافات (+d) و باقیمانده ها (+r). وقتی بیش از یک
از شبکه نمونه برداری می شود این دنباله از 1-3 ستون برای هر شبکه تکرار می شود.
-T[شعاع[u]][+e|p]
برای استفاده در نمونهبرداری از شبکه معمولی، و محدود به یک شبکه منفرد غیر IMG. اگر
نزدیکترین گره به نقطه ورودی NaN است، به سمت بیرون جستجو کنید تا زمانی که آن را پیدا کنیم
نزدیکترین گره غیر NaN و در عوض آن مقدار را گزارش کنید. به صورت اختیاری یک جستجو را مشخص کنید
شعاع که در نظر گرفتن را به نقاطی در این فاصله از ورودی محدود می کند
نقطه. برای گزارش موقعیت نزدیکترین گره و فاصله آن از ورودی
نقطه، الحاق +e. به جای آن نقطه ورودی را با مختصات جایگزین کنید
نزدیکترین گره، ضمیمه +p.
-V[سطح] (بیشتر ...)
سطح پرحرفی [c] را انتخاب کنید.
-Z فقط مقادیر z نمونه برداری شده را بنویسید [پیش فرض همه ستون ها را می نویسد].
-: بین ورودی/خروجی (طول و عرض جغرافیایی) و (طول و عرض جغرافیایی، طول جغرافیایی) تغییر می کند.
[پیشفرض (طول و عرض جغرافیایی)] است.
-bi[Ncols][t] (بیشتر ...)
ورودی باینری بومی را انتخاب کنید. [پیشفرض 2 ستون ورودی است].
-bo[Ncols][نوع] (بیشتر ...)
خروجی باینری بومی را انتخاب کنید. [پیشفرض یکی بیشتر از ورودی است].
-d[i|o]بدون اطلاعات (بیشتر ...)
ستون های ورودی که برابر هستند را جایگزین کنید بدون اطلاعات با NaN و برعکس را در خروجی انجام دهید.
-f[i|o]colinfo (بیشتر ...)
انواع داده های ستون های ورودی و/یا خروجی را مشخص کنید.
-g[a]x|y|d|X|Y|D|[گردنه]z[+|-]شکاف[تو] (بیشتر ...)
شکاف داده ها و شکاف های خط را تعیین کنید.
-h[i|o][n][+c][+d][+rاظهار نظر][+rعنوان] (بیشتر ...)
رد شدن یا تولید رکورد(های) سرصفحه.
-iیقه[l][sمقیاس][oچاپ افست][،...] (بیشتر ...)
ستون های ورودی را انتخاب کنید (0 ستون اول است).
-n[b|c|l|n][+a][+bBC][+c][+tآستانه] (بیشتر ...)
حالت درون یابی را برای شبکه ها انتخاب کنید.
-oیقه[،...] (بیشتر ...)
ستون های خروجی را انتخاب کنید (0 ستون اول است).
-s[یقه][a|r] (بیشتر ...)
مدیریت رکوردهای NaN را تنظیم کنید.
-^ or تنها -
یک پیام کوتاه در مورد نحو دستور چاپ کنید، سپس از آن خارج شوید (توجه: در ویندوز
فقط استفاده کنید -).
-+ or تنها +
یک پیام کاربردی (راهنما) گسترده، از جمله توضیح هر یک را چاپ کنید
گزینه مخصوص ماژول (اما نه گزینه های رایج GMT)، سپس خارج می شود.
-? or نه استدلال
سپس یک پیام استفاده کامل (راهنما)، از جمله توضیح گزینه ها را چاپ کنید
خارج می شود.
- نسخه
نسخه GMT را چاپ کنید و خارج شوید.
--show-datadir
مسیر کامل را به دایرکتوری اشتراک گذاری GMT چاپ کنید و از آن خارج شوید.
یاتاقان ها
برای واحد فاصله نقشه، ضمیمه کنید واحد d برای درجه قوس، m برای دقیقه قوس، و s برای قوس
دوم، یا e برای متر [پیش فرض]، f برای پا، k برای کیلومتر، M برای مایل قانون، n برای دریایی
مایل، و u برای پای بررسی ایالات متحده به طور پیش فرض ما چنین فاصله هایی را با استفاده از یک کروی محاسبه می کنیم
تقریب با دایره های بزرگ آماده کردن - به یک فاصله (یا واحد بدون فاصله است
داده شده) برای انجام محاسبات "زمین مسطح" (سریعتر اما کمتر دقیق) یا پیش فرض + به
انجام محاسبات ژئودزیکی دقیق (آهسته تر اما دقیق تر).
ASCII FORMAT دقیق
فرمت های خروجی ASCII داده های عددی توسط پارامترهای موجود در شما کنترل می شود gmt.conf
فایل. طول و عرض جغرافیایی بر اساس FORMAT_GEO_OUT قالب بندی می شوند، در حالی که سایر موارد
مقادیر بر اساس FORMAT_FLOAT_OUT قالب بندی می شوند. توجه داشته باشید که قالب موجود می تواند
منجر به از دست دادن دقت در خروجی می شود که می تواند منجر به مشکلات مختلفی در پایین دست شود. اگر
متوجه می شوید که خروجی با دقت کافی نوشته نشده است، تغییر به باینری را در نظر بگیرید
خروجی (-بو در صورت وجود) یا اعشار بیشتری را با استفاده از تنظیم FORMAT_FLOAT_OUT مشخص کنید.
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 بعدی.
آثار OF GRID نمونه گیری مجدد
نمونهبرداری مجدد یا نمونهبرداری از شبکهها از الگوریتمهای مختلفی استفاده میکند (نگاه کنید به -n) که ممکن است منجر شود
تحریف های احتمالی یا نتایج غیرمنتظره در مقادیر نمونه گیری مجدد. یک اثر مورد انتظار
نمونهبرداری مجدد با splines تمایل به مقادیر کمی بیشتر از مقادیر نمونهگیری مجدد است
محدودیت جهانی حداقل/حداکثر شبکه اصلی. اگر این غیر قابل قبول است، می توانید تحمیل کنید
برش دادن مقادیر مقادیر مجدد نمونه گیری شده به طوری که آنها از مقادیر حداقل/حداکثر ورودی تجاوز نکنند.
اضافه کردن +c خود را به -n گزینه.
نکات
اگر یک نقطه درون یابی روی یک گره از شبکه ورودی نباشد، یک NaN در هر گره در
همسایگی اطراف نقطه یک NaN درون یابی تولید می کند. دو مکعبی
درون یابی [پیشفرض] مشتقات اول پیوسته را به دست میدهد اما به همسایگی آن نیاز دارد
4 گره در 4 گره. درونیابی دو خطی [-n] فقط از یک محله 2 در 2 استفاده می کند، اما
فقط تداوم مرتبه صفر را ایجاد می کند. هنگامی که صاف بودن مهم است از دو مکعبی استفاده کنید. استفاده کنید
دوخطی برای به حداقل رساندن انتشار NaNs یا کمتر آستانه.
مثال ها
برای نمونهبرداری از فایل hawaii_topo.nc در امتداد مسیر SEASAT track_4.xyg (یک جدول ASCII
شامل طول جغرافیایی، عرض جغرافیایی، و گرانش ناشی از SEASAT، قبل از یک هدر
رکورد):
grdtrack track_4.xyg -Ghawaii_topo.nc -h > track_4.xygt
برای نمونه برداری از فایل Sandwell/Smith IMG topo.8.2.img (2 دقیقه عمق سنجی پیش بینی شده
در یک شبکه Mercator) و شبکه سنی Muller و همکاران.3.2.nc در امتداد مختصات lon،lat
در فایل cruise_track.xy داده شده، امتحان کنید
grdtrack cruise_track.xy -Gtopo.8.2.img,1,1 -Gage.3.2.nc > depths-age.d
برای نمونه برداری از فایل فرمت Sandwell/Smith IMG grav.18.1.img (1 دقیقه ناهنجاری های هوای آزاد در
یک شبکه Mercator) در امتداد نمایه های متقاطع ۱۰۰ کیلومتری که متعامد به بخش خط هستند
در فایل track.xy داده شده است، هر 25 کیلومتر نمایه های متقاطع ایجاد می کند و از شبکه نمونه برداری می کند.
هر 3 کیلومتر امتحان کنید
grdtrack track.xy -Ggrav.18.1.img,0.1,1 -C100k/3/25 -Ar > xprofiles.txt
برای نمونه گیری از grid data.nc در امتداد یک خط از پایین سمت چپ به گوشه سمت راست بالا،
با استفاده از فاصله شبکه 1 کیلومتری و همچنین فواصل خروجی، سعی کنید
grdtrack -ELB/RT+i1k+d -Gdata.nc > profiles.txt
با استفاده از خدمات onworks.net از grdtrackgmt به صورت آنلاین استفاده کنید