img2grdgmt - מקוון בענן

תָכְנִית:

שֵׁם

img2grd - חלץ קבוצת משנה של קובץ img בפורמט Mercator או גיאוגרפי

תַקצִיר

img2grd imgfile grdfile באזור סוג [ ] [ [minlat/maxlat] ] [ ] [ דקות ] [ ] [ ניווט
] [ [סולם] ] [ [רָמָה] ] [ מקסלון ] [ -n]

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

תיאור

img2grd קורא קובץ בפורמט img, מחלץ תת-קבוצה וכותב אותו לקובץ רשת. ה -M
האפשרות מכתיבה אם הקרנת ה-Spherical Mercator של קובץ ה-img היא או לא
נשמר או אם יש לכתוב רשת גיאוגרפית על ידי ביטול הקרנת מרקטור. אם
רשת גיאוגרפית נבחרה, תוכל גם לבקש דגימה מחדש על המדויק -R נָתוּן.

נדרש טיעונים

imgfile
קובץ בפורמט של Mercator img כגון שדות הכבידה הימית או טופוגרפיה של קרקעית הים
מוערך מנתוני מד גובה לוויינים על ידי Sandwell and Smith. אם המשתמש הגדיר
משתנה סביבה $GMT_DATADIR, לאחר מכן img2grd ינסה למצוא imgfile in
$GMT_DATADIR; אחרת הוא ינסה להיפתח imgfile באופן ישיר.

-Ggrdfile
grdfile הוא השם של קובץ רשת הפלט.

-ר[יחידה]מערב/מזרח/דרום/צפון[/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 אוֹפְּצִיָה.
במקרה האחרון משרטטת מראה פרספקטיבי של המטוס, ללא שליש
מֵמַד.

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

-C הגדר את קואורדינטות ה-x ו-y Mercator ביחס למרכז ההקרנה [ברירת המחדל היא
ביחס לפינה השמאלית התחתונה של הרשת]. דורש -M.

-ד[minlat/maxlat]
השתמש בטווח הרוחב המורחב -80.738/+80.738. לחילופין, צרף
minlat/maxlat כהיקף הרוחב של קובץ ה-img הקלט. [ברירת המחדל היא
-72.006/72.006]. לא נדרש בדרך כלל מכיוון שאנו יכולים לקבוע את ההיקף מ
בדיקה של גודל הקובץ.

-E ניתן להשתמש מתי -M לא מוגדר לאלץ את הרשת הסופית לקבל את אותו הדבר בדיוק
אזור כמבוקש עם -R. כברירת מחדל, האזור הסופי הוא השלכה ישירה של
אזור Mercator המקורי ובדרך כלל ימשך מעט מעבר ל
טווח קו הרוחב המבוקש, ויתרה מכך, עליית הרשת בקו הרוחב לא
להתאים את תוספת קו האורך. עם זאת, הדגימה הנוספת מציגה קטן
שגיאות אינטרפולציה ויש להשתמש בהן רק אם רשת הפלט חייבת להתאים ל
האזור המבוקש ויש להם x_inc = y_inc. במקרה זה האזור שנקבע על ידי -R חייב להיות
ניתן בכפולות של התוספת (לדוגמה, -R0/45/45/72).

-I מצביע דקות כרוחב של פיקסל img קלט בדקות קו אורך.
[ברירת המחדל היא 2.0]. בדרך כלל לא נדרש מכיוון שאנו יכולים לקבוע את גודל הפיקסלים לפי
בדיקת הגודל.

-M פלט רשת מרקטור כדורית [ברירת המחדל היא רשת אורך/רוחב גיאוגרפית]. ה
הקרנת Mercator כדורית של קובץ ה-img נשמרת, כך שהאזור -R
שהוגדר על ידי המשתמש שונה מעט; האזור המשתנה מתאים לקצוות
של פיקסלים [או קבוצות של ניווט פיקסלים]. כותרת קובץ הרשת מוגדרת כך שה-x ו-
אורכי ציר y מייצגים מרחק מהקצוות המערביים והדרומיים של התמונה,
נמדד ביחידות ברירת המחדל של המשתמש, עם -Jm1 והמותאם -R. על ידי הגדרת ה
ברירת מחדל PROJ_ אליפסואיד = כדור, המשתמש יכול לעשות שכבות על עם המותאם -R
כדי שיתאימו. לִרְאוֹת דוגמאות לְהַלָן. המותאם -R כתוב גם ברשת
הערת כותרת, כך שניתן למצוא אותה מאוחר יותר. לִרְאוֹת -C להגדיר קואורדינטות ביחס אליהן
מרכז הקרנה.

-Nניווט ממוצע את הערכים בפיקסלים הקלט img לתוך ניווט by ניווט ריבועים, וליצור
פיקסל פלט אחד לכל ריבוע כזה. אם משתמשים בו עם -T3 זה ידווח על ממוצע
אילוץ בין 0 ל-1. אם משתמשים בו עם -T2 הפלט יהיה ערך נתונים ממוצע
או NaN לפי אם האילוץ הממוצע הוא > 0.5. ניווט חייב לחלק באופן שווה
לתוך מידות ה-imgfile בפיקסלים. [בְּרִירַת מֶחדָל 1 לא עושה ממוצע].

-S[סולם]
הכפל את ערכי קובץ ה-IMG ב- סולם לפני אחסון בקובץ רשת. [ברירת המחדל היא
1.0]. עבור קבצי img אחרונים: קבצי img topo מאוחסנים במטרים (מתוקנים) [-S1];
קבצי כבידה באוויר חופשי ב-mGal*10 [-S0.1 כדי לקבל mGal]; קבצי סטייה אנכית פנימה
מיקרו רדיאנים*10 [-S0.1 כדי לקבל מיקרו-רדיאנים], קבצי שיפוע כבידה אנכיים פנימה
Eotvos*50 [-S0.02 כדי לקבל Eotvos, או -S0.002 כדי לקבל mGal/km]). אם לא סולם נתון
אנו מנסים לקבוע את קנה המידה על ידי בחינת שם הקובץ לאיתור רמזים.

-Tסוג סוג מטפל בקידוד של מידע אילוץ. סוג = 0 מציין שלא
מידע כזה מקודד בקובץ img (בשימוש עבור גרסאות שלפני 1995 של
נתוני כבידה) ומקבל את כל הנתונים. סוג > 0 מציין שמידע אילוץ הוא
מקודדים (גרסאות 1995 ואילך (נוכחיות) של קבצי ה-img), כך שאפשר
הפק קובץ רשת באופן הבא: -T1 מקבל ערכי נתונים בכל הנקודות, -T2 מקבל נתונים
ערכים בנקודות מוגבלות ו-NaN בנקודות אינטרפולציה; -T3 מקבל 1 ב
נקודות מוגבלות ו-0 בנקודות אינטרפולציה [ברירת המחדל היא 1].

-V[רָמָה] (יותר ...)
בחר רמת מילוליות [c]. מומלץ במיוחד כאן, כי זה מועיל לראות
כיצד מותאמות הקואורדינטות.

-Wמקסלון
מצביע מקסלון כמידת קו האורך המקסימלית של קובץ ה-IMG הקלט. גרסאות
מאז 1995 יש מקסלון = 360.0, בעוד לכמה קבצים קודמים היו מקסלון = 390.0.
[ברירת המחדל היא 360.0].

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

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

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

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

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

גֵאוֹגרָפִי דוגמאות

אל האני -M יש לא לכלול את האפשרות אם אתה צריך שרשת הפלט תהיה בגיאוגרפית
קואורדינטות. לחילוץ נתונים באזור -R-40/40/-70/-30 מ world_grav.img.7.2 ו
פרוייקט מחדש כדי להניב קואורדינטות גיאוגרפיות, אתה יכול לנסות

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

מכיוון שמרווח קווי הרוחב בקובץ img נמצא במרחק שווה ביחידות מרקטור, ה-
הרשת המתקבלת לא תתאים לצוין -R בדיוק, ומרווח קווי הרוחב לא
שווה למרווח קו האורך. אם אתה צריך התאמה מדויקת עם שלך -R ואותו מרווח
בקווי אורך ובקו רוחב, השתמש ב- -E אפשרות:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

מרקטור דוגמאות

מכיוון שקובצי ה-img נמצאים בהקרנת מרקטור, אין לחלץ רשת גיאוגרפית
אם התוכנית שלך היא ליצור מפת מרקטור. אם עשית את זה אתה מסיים להקרין ו
השלכה מחדש של הרשת, מאבדת פרטים באורך גל קצר. עדיף להשתמש -M ולתכנן את הרשת
באמצעות הקרנה ליניארית באותו קנה מידה כמו הקרנת Mercator הרצויה (ראה GMT
דוגמה 29). לחילוץ נתונים באזור -R-40/40/-70/-30 מ world_grav.img.7.2, לרוץ

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

שים לב כי -V option אומר לנו שהטווח הותאם ל
-R-40/40/-70.0004681551/-29.9945810754. אנחנו יכולים גם להשתמש grdinfo כדי למצוא את קובץ הרשת
הכותרת מציגה את האזור שלו -R0/80/0/67.9666667 זה הטווח של x,y שנקבל
מהקרנת מרקטור כדורית באמצעות -R-40/40/-70.0004681551/-29.9945810754 ו-
-Jm1. לפיכך, כדי לקחת את ship.lonlatgrav ולהשתמש בו כדי לדגום את ה-merc_grav.nc, נוכל לעשות זאת:

gmt set PROJ_ELLIPSOID Sphere

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | gmt mapproject \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

מומלץ להשתמש בשיטה שלעיל להקרנה וביטול ההקרנה של הנתונים בכאלה
יישום, כי אז יש רק שלב אינטרפולציה אחד (ב grdtrack). אם אחד
תחילה מנסה להמיר את קובץ הרשת ל-lon,lat ואז לדגום אותו, יש שניים
שלבי אינטרפולציה (בהמרה ובדגימה).

כדי ליצור רשת lon,lat מהרשת שלמעלה נוכל להשתמש

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

במקרים מסוימים זה לא יהיה קל כמו ה -R בשתי מערכות הקואורדינטות עשוי שלא להתיישר
נו. כשזה קורה, אנחנו יכולים גם להשתמש (למעשה, אולי תמיד עדיף להשתמש)

gmt grd2xyz merc_grav.nc | gmt mapproject \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
משטח gmt -R-40/40/-70/70 -I2m -Ggrav.nc

כדי ליצור מפת Mercator של האזור הנ"ל, נניח שהערך gmt.conf שלנו עבור
PROJ_LENGTH_UNIT הוא אינץ'. ואז מאז הקובץ merc_grav.nc לעיל מוקרן עם -Jm1i
רוחבו 80 אינץ'. אנו יכולים ליצור מפה ברוחב 8 אינץ' על ידי שימוש -Jx0.1i בכל תוכניות מפות
מוחל על רשת זו (למשל, grdcontour, grdimage, grdview), ולאחר מכן עבור שכבות על אשר
עבודה ב-lon,lat (למשל, psxy, pscoast) נוכל להשתמש במותאם לעיל -R ו -Jm0.1 לקבל
שתי המערכות יתאימו.

עם זאת, אנחנו יכולים להיות חכמים יותר מזה. מבין שלקובץ ה-img הקלט יש פיקסלים 2.0
דקות ברוחב (או בדיקת nx ו-ny עם grdinfo merc_grav.nc) אנחנו מבינים את זה
merc_grav.nc השתמש ברזולוציה המלאה של קובץ ה-img ויש לו 2400 על 2039 פיקסלים, ו
ברוחב של 8 אינץ' זה 300 פיקסלים לאינץ'. אנחנו מחליטים שאנחנו לא צריכים כל כך הרבה ואנחנו
נסתפק ב-100 פיקסלים לאינץ', אז אנחנו רוצים לחשב את הנתונים בממוצע ל-3 על 3
ריבועים. (אם נרצה עלילת קו מתאר כנראה שנבחר לעשות ממוצע נתונים הרבה יותר
(למשל, 6 על 6) כדי לקבל קווי מתאר חלקים.) מכיוון ש-2039 אינו מתחלק ב-3, נקבל
שונה מותאם -R הפעם:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

הפעם אנו מוצאים את האזור המותאם הוא -R-40/40/-70.023256525/-29.9368261101 וה-
הפלט הוא 800 על 601 פיקסלים, גודל טוב יותר עבורנו. עכשיו אנחנו יכולים ליצור מלאכותי
קובץ תאורה לשימוש זה גרדיינט:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

ואם יש לנו גם קובץ CPT שנקרא "grav.cpt" נוכל ליצור מפת תבליט בגווני צבע
ככה:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

נניח שאתה רוצה להשיג רק את ערכי הנתונים המוגבלים מקובץ img, ב-lat/lon
קואורדינטות. ואז לרוץ img2grd עם -Tאפשרות 2, השתמש grd2xyz כדי לזרוק את הערכים, צינור
דרך grep -v NaN כדי לחסל NaNs, ומעבר דרך פרויקט map עם היפוך
הקרנה כאמור לעיל.

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