זוהי הפקודה grdrastergmt שניתן להריץ בספק האירוח החינמי של OnWorks באמצעות אחת מתחנות העבודה המקוונות החינמיות שלנו, כגון Ubuntu Online, Fedora Online, אמולטור מקוון של Windows או אמולטור מקוון של MAC OS.
תָכְנִית:
שֵׁם
grdraster - חילוץ תת-אזור מרסטר בינארי ושמירה כרשת GMT
תַקצִיר
גרדרסטר [ קובץ מספר | "טֶקסט תַבְנִית" ] באזור [ grdfile ] [ תוֹסֶפֶת ] [ פרמטרים ]
[[רָמָה] ] [ -בו] [ תעשה] [ -o]
הערה: אין רווח בין דגל האפשרות לארגומנטים המשויכים.
תיאור
גרדרסטר קורא קובץ בשם grdraster.info מתוך ספריית העבודה הנוכחית, ה-
ספריות שאליהן מצביעים משתני הסביבה $GMT_USERDIR ו $GMT_DATADIR, או ב
$GMT_SHAREDIR/dbase (בסדר הזה). הקובץ grdraster.info מגדיר מערכים בינאריים של נתונים
מאוחסן בפורמט סריקה בקבצי נתונים. לכל קובץ ניתנת קובץ מספר במידע
קובץ. גרדרסטר מגלה כיצד לטעון את נתוני הרסטר לקובץ רשת המשתרע על פני אזור
מוגדר על ידי -Rכברירת מחדל, מרווח הרשת שווה למרווח הרסטר. -I אפשרות עשויה להיות
משמש לדגימה משנה של נתוני הרסטר. עם זאת, לא מתבצעים סינון או אינטרפולציה;
x_inc ו y_inc של הרשת חייב להיות כפולות של המרווחים של קובץ הרסטר ו
גרדרסטר פשוט לוקח כל נקודה n. הפלט של גרדרסטר הוא או רשת או פיקסל
רישום בהתאם לרישום הרסטר בו נעשה שימוש. זה תלוי במערכת GMT
אדם לתחזק את grdraster.info קובץ בהתאם לרסטרים הזמינים ב
כל אתר. מערכי נתוני רסטר אינם מסופקים עם GMT אך ניתן להשיגם באמצעות ftp אנונימי.
ועל גבי תקליטור (ראה דף README בספריית dbase). גרדרסטר יפרט את הזמינים
קבצים אם לא ניתנים ארגומנטים. לבסוף, גרדרסטר יכתוב שלישיות xyz ל-stdout אם
לא ניתן שם קובץ רשת פלט
נדרש טיעונים
קובץ מספר
אם מספר שלם התואם לאחד מהקבצים המפורטים ב- grdraster.info הקובץ ניתן
נשתמש במערך הנתונים הזה, אחרת נתאים את תבנית הטקסט הנתונה לנתונים
תיאור הקבוצה על מנת לקבוע את מערך הנתונים.
-ר[יחידה]מערב/מזרח/דרום/צפון[/zmin/zmax][r]
מערב, מזרח, דרום, ו צפון ציין את אזור העניין, ותוכל לציין
אותם במעלות עשרוניות או בפורמט [+-]dd:mm[:ss.xxx][W|E|S|N]. לְצַרֵף r אם נמוך יותר
קואורדינטות המפה השמאלית והימנית העליונה ניתנות במקום w/e/s/n. השניים
קצרנות -רג ו -מחקר ופיתוח מייצג תחום גלובלי (0/360 ו-180/+180 בקו אורך
בהתאמה, עם -90/+90 בקו הרוחב). לחלופין ליצירת רשת, תן
Rקודלון/lat/nx/ny, שם קוד הוא שילוב של 2 תווים של L, C, R (עבור שמאל,
במרכז, או ימינה) ו-T, M, B עבור למעלה, באמצע או למטה. למשל, BL עבור שמאל למטה.
זה מציין איזו נקודה באזור מלבני את לון/lat קואורדינטה מתייחס
ל, וממדי הרשת nx ו ny עם מרווחי רשת דרך -I משמש ליצירה
האזור המקביל. לחלופין, ציין את השם של קובץ רשת קיים
ו -R הגדרות (ומרווח רשת, אם רלוונטי) מועתקות מהרשת.
שימוש -Rיחידה מצפה לקואורדינטות מוערכות (קרטזיות) תואמות לבחירה -J
ואנחנו מקרינים באופן הפוך כדי לקבוע אזור גיאוגרפי מלבני בפועל. ל
נקודת מבט (-p), אופציונלי להוסיף /zmin/zmax. במקרה של נקודת מבט
(-p), טווח z (zmin, zmax) ניתן להוסיף לציון הממד השלישי. זֶה
צריך להיעשות רק בעת שימוש ב -Jz אפשרות, לא כאשר משתמשים רק ב- -p אוֹפְּצִיָה.
במקרה האחרון משרטטת מראה פרספקטיבי של המטוס, ללא שליש
ממד. אם r מצורף, ניתן גם לציין היטל מפה כדי להגדיר את
צורת האזור שלך. אזור הפלט יעוגל למספר השלם הקרוב ביותר
צעד רשת בשני הממדים.
אופציונאלי טיעונים
-Ggrdfile
שם קובץ רשת הפלט. אם לא מוגדר, הרשת תיכתב כ-ASCII (או בינארי;
לִרְאוֹת -בו) שלישיות xyz ל- stdout במקום זאת.
-Iשינק[יחידה][=|+][/yinc[יחידה[+]]
x_inc [ובאופציה y_inc] הוא מרווח הרשת. לחלופין, הוסף סיומת
שינוי. גיאוגרפי (מעלות) קואורדינטות: צרף m כדי לציין דקות קשת או s
כדי לציין שניות קשת. אם אחת מהיחידות e, f, k, M, n or u מצורף
במקום זאת, ההנחה היא שהתוספת ניתנת במטר, רגל, ק"מ, מייל, ימי
מייל או רגל סקר בארה"ב, בהתאמה, ויומרו למקבילה
מעלות קו אורך בקו הרוחב האמצעי של האזור (ההמרה תלויה ב
PROJ_ELLIPSOID). אם /y_inc נתון אך מוגדר ל-0 הוא יאופס שווה ל x_inc;
אחרת הוא יומר לקו רוחב מעלות. הכל קואורדינטות: אם = is
מצורף ואז המקסימום המתאים x (מזרח) או y (צפון) עשוי להיות מותאם מעט
כדי להתאים בדיוק לתוספת הנתונה [כברירת מחדל, התוספת עשויה להיות מותאמת
מעט כדי להתאים לתחום הנתון]. לבסוף, במקום לתת תוספת אפשר
ציין את מספר of צמתים רצוי על ידי הוספה + למספר השלם שסופק
טַעֲנָה; לאחר מכן, התוספת מחושבת מחדש ממספר הצמתים וה-
תְחוּם. ערך התוספת המתקבל תלוי אם בחרת ב-a
רשת רשומה קו רשת או פיקסל; ראה פורמטים של קבצי אפליקציה לפרטים.
הערה: אם -Rgrdfile נמצא בשימוש אז מרווח הרשת כבר אותחל; להשתמש
-I לעקוף את הערכים.
-Jפרמטרים (יותר ...)
בחר הקרנת מפה.
-V[רָמָה] (יותר ...)
בחר רמת מילוליות [c].
-בו[ncols][סוג] (יותר ...)
בחר פלט בינארי מקורי.
תעשהאין מידע (יותר ...)
החלף עמודות פלט השוות ל-NaN עם אין מידע.
-oצווארונים[,...] (יותר ...)
בחר עמודות פלט (0 הוא העמודה הראשונה).
אפשרות זו חלה רק אם לא -G הוגדרה אפשרות.
-^ or רק -
הדפס הודעה קצרה על התחביר של הפקודה, ואז צא (הערה: ב-Windows
להשתמש רק -).
-+ or רק +
הדפס הודעת שימוש נרחבת (עזרה), כולל הסבר על כל
אפשרות ספציפית למודול (אך לא האפשרויות הנפוצות של GMT), ואז יוצאת.
-? or לא טיעונים
לאחר מכן הדפס הודעת שימוש מלאה (עזרה), כולל הסבר על האפשרויות
יציאות.
--גִרְסָה
הדפס גרסת GMT וצא.
--show-datadir
הדפס את הנתיב המלא לספריית השיתוף של GMT וצא.
דוגמאות
כדי לחלץ נתונים מרסטר 1, תוך לקיחת נקודה אחת כל 30 דקות, באזור מורחב
מעבר ל-360 מעלות כדי לאפשר סינון מאוחר יותר, הפעל
gmt grdraster 1 -R-4/364/-62/62 -I30m -Gdata.nc
כדי להשיג נתונים עבור היטל מרקטור אלכסוני עלינו לחלץ נתונים נוספים שהם
בשימוש בפועל. זה הכרחי מכיוון שהפלט של גרדרסטר יש לו קצוות המוגדרים על ידי
מקבילים ומרידיאנים, בעוד שהמפה האלכסונית באופן כללי לא עושה זאת. לכן, כדי לקבל את כל
נתונים מנתוני ETOPO2 הדרושים ליצירת מפת קווי מתאר עבור האזור המוגדר על ידי החלק התחתון שלו
פינות שמאליות ועליונות והימניות הרצויות, השתמשו
gmt grdraster ETOPO2 -R160/20/220/30r -Joc190/25.5/292/69/1 -Gdata.nc
כדי לחלץ נתונים מתערובת התבליט של Geoware ב-2 דקות ולכתוב אותם כקובץ בינארי כפול
שלישיות xyz מדויקות לפלט סטנדרטי:
gmt grdraster "2 דקות גיאו-וור" -R20/25/-10/5 -bo > triplets.b
השתמש ב-grdrastergmt באופן מקוון באמצעות שירותי onworks.net
