grdgradientgmt - מקוון בענן

תָכְנִית:

שֵׁם

grdgradient - חישוב נגזרת כיוונית או שיפוע מרשת

תַקצִיר

גרדיינט in_grdfile out_grdfile [ Azim[/azim2] ] [ [a][c][o][n] ] [
[s|עמ']אזים/עיל[/סביבה/לפזר/ספקולרי/לזרוח] ] [ דגל ] [ [e][t][amp][/סיגמא[/לקזז]]
] [ באזור ] [ slopefile ] [ [רָמָה] ] [ -fg ] [ -n]

הערה: אין רווח בין דגל האפשרות לארגומנטים המשויכים.

תיאור

גרדיינט ניתן להשתמש כדי לחשב את הנגזרת הכיוונית בכיוון נתון (-A),
או הכיוון (-S) [והגודל (-D)] של שיפוע הווקטור של הנתונים.

ערכים משוערים בשורה/עמודה ראשונה/אחרונה של הפלט תלויים בתנאי גבול (ראה
-L).

נדרש טיעונים

in_grdfile
קובץ רשת דו-ממדי שממנו ניתן לחשב נגזרת כיוונית. (ראה פורמטים של קבצי GRID
להלן).

-Gout_grdfile
שם קובץ רשת הפלט עבור הנגזרת הכיוונית. (ראה פורמטים של קבצי GRID
להלן).

אופציונאלי טיעונים

-AAzim[/azim2]
כיוון אזימוטלי עבור נגזרת כיוונית; Azim היא הזווית ב-x,y
מישור נמדד במעלות חיוביות בכיוון השעון מצפון (כיוון +y) לכיוון
מזרח (כיוון +x). השלילי של נגזרת הכיוונית,
-[dz/dx*sin(Azim) + dz/dy*cos(Azim)], נמצא; שלילה מניבה ערכים חיוביים
כאשר השיפוע של z(x,y) הוא במורד ב- Azim כיוון, החוש הנכון עבור
הצללת הארה של תמונה (ראה grdimage ו grdview) על ידי מקור אור
מעל המישור x,y הבוהק מה- Azim כיוון. אופציונלי, ספק שניים
אזימוטים, -AAzim/azim2, ובמקרה זה ההדרגות בכל אחד מהכיוונים הללו הם
מחושב והגדול בגודלו נשמר; זה שימושי עבור
מאיר נתונים עם שני כיוונים של מבנים מסודרים, למשל, -A0/270
מאיר מצפון (למעלה) וממערב (משמאל).

-D[a][c][o][n]
מצא את כיוון השיפוע החיובי (במעלה) של הנתונים. במקום זאת
מצא את ההיבט (הכיוון של המדרון), השתמש -דא. כברירת מחדל, הנחיות הן
נמדד בכיוון השעון מצפון, כמו Azim in -A מֵעַל. לְצַרֵף c להשתמש קונבנציונאלי
זוויות קרטזיות נמדדות נגד כיוון השעון מכיוון ה-x החיובי (מזרח).
צרף o לדווח על כיוונים (0-180) במקום כיוונים (0-360). לְצַרֵף n
כדי להוסיף 90 מעלות לכל הזוויות (למשל, כדי לתת פגיעות מקומיות של פני השטח).

-E[s|p]אזים/עיל[/סביבה/לפזר/ספקולרי/לזרוח]
חישוב זוהר למברטיאני מתאים לשימוש איתו grdimage ו grdview. ה
השתקפות למברטיאן מניחה משטח אידיאלי שמחזיר את כל האור ש
פוגע בו והמשטח נראה בהיר באותה מידה מכל כיווני הצפייה. Azim
ו elev הם האזימוט והגובה של וקטור האור. אופציונלי, אספקה סביבה
לפזר ספקולרי לזרוח שהם פרמטרים השולטים במאפייני ההחזר
של פני השטח. ערכי ברירת המחדל הם: 0.55/0.6/0.4/10 להשאיר חלק מהערכים
untouched, ציין = כערך החדש. לדוגמה -E60/30/=/0.5 מגדיר את Azim elev
ו לפזר ל-60, 30 ו-0.5 ומשאיר את שאר פרמטרי ההחזר
ללא נגיעה. לְצַרֵף s להשתמש באלגוריתם למברטיאני פשוט יותר. שימו לב שעם טופס זה
אתה רק צריך לספק את פרמטרי האזימוט והגובה. לְצַרֵף p להשתמש
קירוב ליניארי של Peucker חלקית (אלגוריתם פשוט יותר אך מהיר יותר; במקרה זה
מה היא Azim ו elev מחוברים ל-315 ו-45 מעלות. זה אומר שגם אם אתה
לספק ערכים אחרים שהם יתעלמו מהם.)

-Lדגל מצב גבול דגל יכול להיות x or y or xy המציין שהנתונים הם תקופתיים בטווח של
x או y או שניהם, או דגל יכול להיות g המציין תנאים גיאוגרפיים (x ו-y הם
lon and lat). [ברירת המחדל משתמשת בתנאים "טבעיים" (נגזרת חלקית שנייה נורמלית
לקצה הוא אפס).]

-נֶטוֹ][amp][/סיגמא[/לקזז]]
נוֹרמָלִיזָצִיָה. [ברירת מחדל: אין נורמליזציה.] ההדרגות בפועל g מקוזזים ו
קנה מידה כדי לייצר שיפועים מנורמלים gn עם גודל פלט מקסימלי של amp.
If amp לא ניתן, ברירת מחדל amp = 1. אם לקזז לא ניתן, הוא מוגדר ל-
ממוצע של g. -N תשואות gn = amp * (g - לקזז)/max(abs(g - לקזז)). -לא
מנרמל באמצעות חלוקת Laplace מצטברת gn = amp * (1.0 -
exp (מ"ר(2) * (g - לקזז)/ סיגמא)) איפה סיגמא מוערך באמצעות נורמת L1 של
(g - לקזז) אם לא ניתן. -Nt מנרמל באמצעות Cauchy מצטבר
תפוקת הפצה gn = (2* amp / PI) * atan( (g - לקזז)/ סיגמא) איפה סיגמא
מוערך באמצעות נורמת L2 של (g - לקזז) אם לא ניתן.

-ר[יחידה]xmin/xmax/ymin/ymax[ר] (יותר ...)
ציין את אזור העניין. משתמש ב -R האפשרות תבחר תת-סעיף של
in_grdfile רֶשֶׁת. אם סעיף קטן זה חורג מגבולות הרשת, רק את
אזור משותף ייחלץ.

-Sslopefile
שם של קובץ רשת פלט עם גדלים סקלאריים של וקטורים גרדיאנטים. דורש -D
אלא עושה -G אופציונאלי.

-V[רָמָה] (יותר ...)
בחר רמת מילוליות [c].

-fg רשתות גיאוגרפיות (מידות קו אורך, קו רוחב) יומרו למטרים
באמצעות קירוב "אדמה שטוחה" תוך שימוש בפרמטרים האליפסואידים הנוכחיים.

-n[b|c|l|n][+a][+bBC][+c][+tסף] (יותר ...)
בחר מצב אינטרפולציה עבור רשתות.

-^ or רק -
הדפס הודעה קצרה על התחביר של הפקודה, ואז צא (הערה: ב-Windows
להשתמש רק -).

-+ or רק +
הדפס הודעת שימוש נרחבת (עזרה), כולל הסבר על כל
אפשרות ספציפית למודול (אך לא האפשרויות הנפוצות של GMT), ואז יוצאת.

-? or לא טיעונים
לאחר מכן הדפס הודעת שימוש מלאה (עזרה), כולל הסבר על האפשרויות
יציאות.

--גִרְסָה
הדפס גרסת GMT וצא.

--show-datadir
הדפס את הנתיב המלא לספריית השיתוף של GMT וצא.

GRID מרחק יחידות

אם לרשת אין מטר כיחידה האופקית, הוסף +uיחידה לקובץ הקלט
שם להמרה מהיחידה שצוינה למטר. אם הרשת שלך היא גיאוגרפית, המר
מרחקים למטרים על ידי אספקה -fg במקום.

הצבעים

אם אתה לא יודע מה -N אפשרויות לשימוש כדי ליצור עבורו קובץ אינטנסיביות grdimage or
grdview, ניסיון ראשון טוב הוא -לא0.6.

בדרך כלל 255 גוונים הם די והותר למטרות הדמיה. אתה יכול לחסוך 75% בדיסק
רווח על ידי הוספה של =nb/a לשם קובץ הפלט out_grdfile.

אם אתה רוצה לעשות כמה מפות מוארות של תת-אזורים של מערך נתונים גדול, ואתה
צריך שאפקטי התאורה יהיו עקביים בכל המפות, השתמש ב- -N אפשרות ו
לספק את אותו הערך של סיגמא ו לקזז ל גרדיינט עבור כל מפה. ניחוש טוב הוא
לקזז = 0 ו- סיגמא נמצא על ידי grdinfo -L2 or -L1 מיושם על גרדיאנט לא מנורמל.

אם אתה פשוט צריך את x- או y-נגזרות של הרשת, שימוש גרדמאת.

GRID קובץ פורמטים

כברירת מחדל, GMT כותב את הרשת כצף דיוק יחיד ב-NETCDF של תלונת COARDS
פורמט קובץ. עם זאת, GMT מסוגל לייצר קבצי רשת בהרבה רשתות נפוצות אחרות
פורמטים של קבצים וגם מקל על מה שנקרא "אריזה" של רשתות, כתיבת נקודה צפה
נתונים כמספרים שלמים של 1 או 2 בתים. כדי לציין את הדיוק, קנה המידה וההיסט, המשתמש צריך
להוסיף את הסיומת =id[/סולם/לקזז[/נאן]], איפה id הוא מזהה בן שתי אותיות של הרשת
סוג ודיוק, ו סולם ו לקזז הם גורם קנה מידה אופציונלי וקיזוז להיות
מיושם על כל ערכי הרשת, ו נאן הוא הערך המשמש לציון נתונים חסרים. במקרה
שתי הדמויות id לא מסופק, כמו ב =/סולם מ id=nf מונחת. מתי
קריאת רשתות, בדרך כלל הפורמט מזוהה אוטומטית. אם לא, אותה סיומת
ניתן להוסיף לשמות קבצי רשת קלט. לִרְאוֹת grdconvert ו-Section grid-fil-פורמט של
עיון טכני של GMT וספר בישול למידע נוסף.

בעת קריאת קובץ netCDF המכיל רשתות מרובות, GMT יקרא, כברירת מחדל, את
רשת דו-ממדית ראשונה שיכולה למצוא בקובץ הזה. לשדל את GMT לקרוא אחר
משתנה רב ממדי בקובץ הגריד, צרף ?varname לשם הקובץ, איפה
varname הוא שם המשתנה. שימו לב שייתכן שתצטרכו לברוח מהמשמעות המיוחדת
of ? בתוכנית המעטפת שלך על ידי הצבת קו נטוי אחורי לפניה, או על ידי הצבת
שם קובץ וסיומת בין מרכאות או מרכאות כפולות. ה ?varname ניתן להשתמש גם בסיומת
עבור רשתות פלט כדי לציין שם משתנה שונה מברירת המחדל: "z". לִרְאוֹת
grdconvert ו-Sections modifiers-for-CF ו-grid-fil-format של GMT Technical
עיון וספר בישול למידע נוסף, במיוחד כיצד לקרוא חיבורים של 3-,
רשתות 4 או 5 מימדיות.

דוגמאות

כדי ליצור קובץ להארת הנתונים ב-geoid.nc באמצעות מעברי מעבר מנורמלים ב
הטווח [-0.6,0.6] המחקה מקורות אור בכיוון צפון ומערב:

gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V

כדי למצוא את כיווני האזימוט של בד קרקעית הים בקובץ topo.nc:

gmt grdgradient topo.nc -Dno -Gazimuths.nc -V

ביבליוגרפיה

Horn, BKP, Hill-Shading and the Reflectance Map, Proceedings of the IEEE, Vol. 69, מס'
1, ינואר 1981, עמ' 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

השתמש ב-grdgradientgmt באינטרנט באמצעות שירותי onworks.net