grdgradientgmt - מקוון בענן

תָכְנִית:

שֵׁם

grdgradient - חישוב נגזרת כיוונית או גרדיאנט מרשת

תַקצִיר

גרדיינט קובץ_in_grd קובץ_out_grd [ Azim[/אזים2] ] [ [a][c][o][n] ] [
[ש|פ]אזים/עלה[/סביבה/לפזר/ספקולרי/לזרוח] ] [ דגל ] [ [e][t][amp][/סיגמא[/לקזז]]
] [ באזור ] [ סלופפיל ] [ [רָמָה] ] [ -fg ] [ -n]

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

תיאור

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

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

נדרש טיעונים

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

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

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

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

-ד[א][ג][או][נ]
מצא את כיוון הגרדיאנט החיובי (במעלה השיפוע) של הנתונים. במקום זאת
מצא את ההיבט (כיוון השיפוע היורד), השתמש -דאכברירת מחדל, הוראות הגעה הן
נמדד עם כיוון השעון מצפון, כפי ש 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 להשאיר חלק מהערכים
ללא שינוי, ציין = כערך החדש. לדוגמה -E60/30/=/0.5 מגדיר את Azim elev
ו לפזר ל-60, 30 ו-0.5 ומשאיר את פרמטרי ההחזרה האחרים
לא נגעו. הוסף s להשתמש באלגוריתם למברטי פשוט יותר. שימו לב שעם צורה זו
אתה רק צריך לספק את פרמטרי האזימוט והגובה. הוסף p להשתמש
קירוב ליניארי חלקית לפי פוקר (אלגוריתם פשוט אך מהיר יותר; במקרה זה
מה היא Azim ו elev מחוברים בקפידה ל-315 ו-45 מעלות. משמעות הדבר היא שגם אם אתם
לספק ערכים אחרים, הם יתעלמו מהם.)

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

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

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

-Sסלופפיל
שם קובץ רשת הפלט עם גודל סקלרי של וקטורי גרדיאנט. דורש -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_grd.

אם אתם רוצים ליצור מספר מפות מוארות של תת-אזורים של מערך נתונים גדול, ואתם
צריך שאפקטי התאורה יהיו עקביים בכל המפות, השתמש ב -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 באמצעות גרדיאנטים מנורמלים exp ב-
הטווח [-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

ביבליוגרפיה

הורן, BKP, הצללת גבעות ומפת ההחזרה, הליכי ה-IEEE, כרך 69, מס'
1, ינואר 1981, עמ' 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

השתמש ב- grdgradientgmt באופן מקוון באמצעות שירותי onworks.net