mapprojectgmt - מקוון בענן

תָכְנִית:

שֵׁם

mapproject - בצע טרנספורמציות מפות קדימה והפוכות, המרות נתונים וגיאודזיה

תַקצִיר

פרויקט map [ שולחנות ] פרמטרים באזור [ b|B|f|F|o|O[0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[נתון] ] [ [יחידה] ] [ [x0/y0/][[+|-]יחידה][+|-] ] [ ] [ line.xy[/[+|-]יחידה][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]החל מ-[/ל] ] [ [רָמָה] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

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

תיאור

פרויקט map קורא (קו אורך, קו רוחב) עמדות מ שולחנות [או קלט סטנדרטי] ו
מחשב (x,y) קואורדינטות באמצעות הקרנת המפה והסולמות שצוינו. אופציונלי, זה
יכול לקרוא (x,y) מיקומים ולחשב ערכי (קו אורך, רוחב) תוך ביצוע הפוך
טרנספורמציה. זה יכול לשמש כדי להפוך נקודות ליניאריות (x,y) המתקבלות על ידי דיגיטציה
מפה של הקרנה ידועה לקואורדינטות גיאוגרפיות. יכול גם לחשב מרחקים לאורך
מסלול, לנקודה קבועה, או הגישה הקרובה ביותר לקו. לבסוף, ניתן להשתמש בו לביצוע
המרות נתונים שונים. שדות נתונים נוספים מותרים לאחר 2 העמודות הראשונות
אשר חייב להיות (קו אורך, קו רוחב) או (x, y). ראה אפשרות -: על איך לקרוא
קבצי (קו רוחב, קו אורך).

נדרש טיעונים

-Jפרמטרים (יותר ...)
בחר הקרנת מפה.

-ר[יחידה]xmin/xmax/ymin/ymax[ר] (יותר ...)
ציין את אזור העניין. מקרה מיוחד להקרנת UTM: אם -C משמש
ו -R לא נתון אז האזור מוגדר לחפוף לאזור UTM הנתון אז
באשר לשמירת הפתרון האליפסואידי המלא (ראה הגבלות למידע נוסף
מֵידָע).

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

שולחן ASCII אחד או יותר (או בינארי, ראה -דוּ[ncols][סוג]) קבצי טבלת נתונים המכילים א
מספר עמודות הנתונים. אם לא ניתנות טבלאות, אנו קוראים מקלט סטנדרטי.

-Ab|B|f|F|o|O[0/lat0]
-אף מחשב את האזימוט (קדימה) מנקודה קבועה lon/lat לכל נקודת נתונים.
השתמש -אב כדי לקבל אחורה אזימוט מנקודות נתונים לנקודה קבועה. להשתמש -או כדי לקבל
כיוונים (-90/90) ולא אזימוטים (0/360). אוֹתִיוֹת גְדוֹלוֹת F, B or O יצטרך
המרת קווי רוחב גיאודטיים לגיאוצנטריים והערכת אזימוט של גיאודזה
(בהנחה שהאליפסואיד הנוכחי אינו כדור). אם לא ניתנת נקודה קבועה אז
אנו מחשבים את האזימוט (או אחורי האזימוט) מהנקודה הקודמת.

-C[dx/dy]
הגדר את מרכז הקואורדינטות המוקרנות להיות במרכז הקרנת המפה [ברירת המחדל היא
בפינה השמאלית התחתונה]. לחלופין, הוסף קיזוזים ביחידות הצפויות שיתווספו (או
מופחת מתי -I מוגדר) ל (מ) הקואורדינטות הצפויות, כגון false
מזרחים וצפוניים עבור אזורי הקרנה מסוימים [0/0]. היחידה המשמשת עבור
קיזוז הוא יחידת מרחק העלילה בתוקף (ראה PROJ_LENGTH_UNIT) אלא אם כן -F is
בשימוש, ובמקרה זה הקיזוזים הם במטרים.

-Dc|i|p
עוקף זמנית את PROJ_LENGTH_UNIT והשתמש c (ס"מ), i (אינץ'), או p (נקודות)
במקום זאת. לא ניתן להשתמש עם -F.

-E[נתון]
המר מגיאודטי (לאורך, רוחב, גובה) לכדור הארץ קבוע (ECEF)
(x,y,z) קואורדינטות (הוסף -I עבור ההמרה ההפוכה). הוסף מזהה נתונים (ראה -Qd)
או לתת אליפסואיד:dx,dy,dz איפה אליפסואיד יכול להיות מזהה אליפסואיד (ראה -Qe) או
נתון כ a[,*inv_f*], איפה a הוא הציר הראשי למחצה ו inv_f הוא הפוך
רידוד (0 אם הושמט). אם נתון הוא - או לא נתון, אנו מניחים WGS-84.

-F[יחידה]
כוח קנה מידה 1:1, כלומר פלט (או קלט, ראה -I) הנתונים מוצגים בפועל
מטרים. כדי לציין יחידות אחרות, הוסף את היחידה הרצויה (ראה יחידות). לְלֹא -F,
הפלט (או הקלט, ראה -I) נמצאים ביחידות שצוינו על ידי PROJ_LENGTH_UNIT (אבל
לִרְאוֹת -D).

-G[x0/y0/][[+|-]יחידה][+|-]
חשב מרחקים לאורך המסלול or לנקודה האופציונלית שהוגדרה עם -Gx0/y0. לְצַרֵף
יחידת המרחק (ראה יחידות), כולל c (מרחק קרטזי באמצעות קלט
קואורדינטות) או C (מרחק קרטזי באמצעות קואורדינטות מוקרנות). ה C יחידה
דורש -R ו -J להיות מוכן. ללא נקודה קבועה אנו מחשבים מצטבר
מרחקים לאורך המסילה. לְצַרֵף - כדי להשיג מרחק מצטבר בין עוקבים
נקודות. לְצַרֵף + כדי לציין את הנקודה השנייה באמצעות שתי עמודות נוספות בקובץ הקלט.

-I בצע את הטרנספורמציה ההפוכה, כלומר קבל (קו אורך, קו רוחב) מנתוני (x,y).

-Lline.xy[/[+|-]יחידה][+]
קבע את המרחק הקצר ביותר מנקודות נתוני הקלט לקו/ים המופיעים ב
קובץ ASCII מרובה מקטעים line.xy. המרחק והקואורדינטות של ה
הנקודה הקרובה ביותר תצורף לפלט כשלוש עמודות חדשות. צרף את
יחידת מרחק (ראה יחידות), כולל c (מרחק קרטזי באמצעות קואורדינטות קלט)
or C (מרחק קרטזי באמצעות קואורדינטות מוקרנות). ה C יחידה דורשת -R ו
-J להיות מוכן. לבסוף, צרף + לדווח על מזהה קטע השורה והשבר
מספר נקודה במקום lon/lat של הנקודה הקרובה ביותר.

-N[a|c|g|m]
המר מקווי רוחב גאודטיים (באמצעות האליפסואיד הנוכחי; ראה PROJ_ELLIPSOID)
לאחד מארבעה קווי רוחב שונים (קוי האורך אינם מושפעים). בחר
החל מ- aאוטלי, cלא פורמלי, gאאוצנטרי, ו mקווי רוחב ארדיוניים [גיאוצנטריים]. להשתמש -I
להמיר מקווי רוחב עזר לקווי רוחב גיאודטיים.

-ש[ד|ה רשום את כל פרמטרי ההקרנה. לרשימת נתונים בלבד, השתמש -Qd. רק לרשימה
אליפסואידים, שימוש -Qe.

-S הדחיק נקודות הנופלות מחוץ לאזור.

-ת[ח]החל מ-[/ל]
תיאום המרות בין נתונים החל מ- ו ל באמצעות התקן מולודנסקי
טרנספורמציה. להשתמש -ת אם לעמודת הקלט השלישית יש גובה מעל אליפסואיד [ברירת מחדל
מניח שגובה = 0, כלומר, על האליפסואיד]. ציין נתונים באמצעות מזהה הדאטום (ראה
-Qd) או לתת אליפסואיד:dx,dy,dz איפה אליפסואיד יכול להיות מזהה אליפסואיד (ראה -Qe) או
נתון כ a[,*inv_f*], איפה a הוא הציר הראשי למחצה ו inv_f הוא הפוך
רידוד (0 אם הושמט). אם נתון הוא - או לא נתון, אנו מניחים WGS-84. -T יכול להיות
בשימוש בשילוב עם -R -J כדי לשנות את הנתון לפני הקרנת קואורדינטות
(לְהוֹסִיף -I כדי להחיל את המרת הנתונים לאחר ההקרנה ההפוכה). תוודא ש
ההגדרה PROJ_ELLIPSOID נכונה עבור המקרה שלך.

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

-W[w|h]
מדפיס רוחב וגובה המפה בפלט סטנדרטי. לא נקראים קבצי קלט. רק ל
פלט את הרוחב או הגובה, הוסף w or h, בהתאמה. היחידות של ה
מידות עשויות להשתנות באמצעות -D.

-דוּ[ncols][t] (יותר ...)
בחר קלט בינארי מקורי. [ברירת המחדל היא 2 עמודות קלט].

-בו[ncols][סוג] (יותר ...)
בחר פלט בינארי מקורי. [ברירת המחדל זהה לקלט].

-d[i|o]אין מידע (יותר ...)
החלף עמודות קלט שוות אין מידע עם NaN ועשה הפוך על הפלט.

-f[i|o]colinfo (יותר ...)
ציין סוגי נתונים של עמודות קלט ו/או פלט.

-g[a]x|y|d|X|Y|D|[col]z[+|-]פער[u] (יותר ...)
קבע פערי נתונים ומעברי שורות.

-h[i|o][n][+c][+d][+rהֶעָרָה][+rכותרת] (יותר ...)
דלג או הפק רשומות כותרות.

-iצווארונים[ל][ססולם][oלקזז][,...] (יותר ...)
בחר עמודות קלט (0 הוא העמודה הראשונה).

-oצווארונים[,...] (יותר ...)
בחר עמודות פלט (0 הוא העמודה הראשונה).

-p[x|y|z]Azim/elev[/zlevel][+w0/lat0[/z0]][+vx0/y0] (יותר ...)
בחר תצוגת פרספקטיבה.

-ס[צווארונים][a|r] (יותר ...)
הגדרת טיפול ברשומות NaN.

-:[i|o] (יותר ...)
החלף את העמודה הראשונה והשנייה בקלט ו/או פלט.

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

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

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

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

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

יחידות

ליחידת מרחק מפה, צרף יחידה d עבור דרגת קשת, m עבור דקת קשת, ו s עבור קשת
שני, או e עבור מד [ברירת מחדל], f לרגל, k עבור ק"מ, M עבור מייל חוקי, n עבור ימיים
מייל, ו u עבור רגל סקר בארה"ב. כברירת מחדל אנו מחשבים מרחקים כאלה באמצעות כדור
קירוב עם עיגולים גדולים. הכנה מראש - למרחק (או שהיחידה אינה מרחק הוא
נתון) כדי לבצע חישובי "אדמה שטוחה" (מהירים יותר אך פחות מדויקים) או להקדים + ל
לבצע חישובים גיאודזיים מדויקים (איטיים יותר אך מדויקים יותר).

ASCII פורמט דיוק

פורמטים של פלט ASCII של נתונים מספריים נשלטים על ידי פרמטרים שלך gmt.conf
קוֹבֶץ. קווי אורך וקווי רוחב מעוצבים לפי FORMAT_GEO_OUT, בעוד אחרים
הערכים מעוצבים לפי FORMAT_FLOAT_OUT. שים לב שהפורמט בפועל יכול
להוביל לאובדן דיוק בתפוקה, מה שעלול להוביל לבעיות שונות במורד הזרם. אם
אתה מוצא שהפלט לא כתוב בדיוק מספיק, שקול לעבור לבינארי
פלט (-בו אם זמין) או ציין יותר עשרוניות באמצעות ההגדרה FORMAT_FLOAT_OUT.

דוגמאות

כדי להמיר קואורדינטות UTM במטרים למיקומים גיאוגרפיים, בהינתן קובץ utm.txt ו
לדעת את אזור ה-UTM (והאזור או ההמיספרה), נסה

gmt mapproject utm.txt -Ju+11/1:1 -C -I -F

כדי להפוך קובץ עם (קו אורך, קו רוחב) למיקומי (x,y) בס"מ על מרקטור
רשת עבור קנה מידה נתון של 0.5 ס"מ לכל מעלה, הפעל

gmt mapproject lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

כדי להפוך כמה קבצים בני 2 עמודות, בינאריים, עם דיוק כפול עם (קו רוחב, קו אורך)
למיקומים (x,y) באינץ' על רשת מרקטור רוחבית (קו אורך מרכזי 75W) עבור
קנה מידה = 1:500000 ודכא את הנקודות הללו שיפלו מחוץ לאזור המפה, רוץ

gmt mapproject tracks.* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

להמרת הקואורדינטות הגיאודטיות (אורך, רוחב, גובה) בקובץ old.dat מה-NAD27
נתון CONUS (נתון מזהה 131 המשתמש באליפסואיד Clarke-1866) ל-WGS 84, הפעלה

gmt mapproject old.dat -Th131 > new.dat

כדי לחשב את המרחק הקרוב ביותר (בק"מ) בין כל נקודה בקובץ הקלט quakes.dat
ואת מקטעי הקו הניתנים בקובץ ASCII מרובה מקטעים coastline.xy, הפעל

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

הגבלות

אזור הקלט המלבני מוגדר עם -R באופן כללי ימופו לא מלבני
רֶשֶׁת. אֶלָא אִם -C מוגדר, לנקודה השמאלית ביותר ברשת זו יש xvalue = 0.0, וה-
לנקודה הנמוכה ביותר יהיה yvalue = 0.0. לכן, לפני שאתה עושה דיגיטציה של מפה, הפעל את האקסטרים
מפה קואורדינטות דרך פרויקט map באמצעות הסולם המתאים וראה אילו (x,y) ערכי
הם ממופים על. השתמש בערכים אלה בעת הגדרת דיגיטציה על מנת לקבל את
טרנספורמציה הפוכה פועלת כהלכה, או לחילופין, השתמש awk לקנה מידה ולהזיז את
(x,y) ערכים לפני ההמרה.

עבור הקרנה מסוימת, ניתן להשתמש בפתרון כדורי למרות שהמשתמש בחר ב-
אליפסואיד. זה מתרחש כאשר המשתמשים -R הגדרה מרמזת על אזור החורג מהדומיין
שבו הרחבות הסדרה אליפסואידית תקפות. אלו התנאים: (1)
למברט קונפורמל קוניק (-JL)ו- Albers Equal-Area (-ג'יי) ישתמש בתמיסה הכדורית
כאשר קנה המידה של המפה עולה על 1.0E7. (2) מרקטור רוחבי (ג'ט) ו-UTM (-JU) ירצה
השתמש בפתרון הכדורי כאשר הגבול המערבי או המזרחי נתון -R זה יותר מ
10 מעלות מהמרידיאן המרכזי, ו-(3) זהה עבור קאסיני (-JC) אבל עם מגבלה של
רק 4 מעלות.

אליפסואידים ו ספירואידים

GMT ישתמש בנוסחאות אליפסואידיות אם הן יושמו והמשתמש בחר ב-
אליפסואיד כצורת הייחוס (ראה PROJ_ELLIPSOID). המשתמש צריך להיות מודע לכמה
מלכודות פוטנציאליות: (1) עבור תחזיות מסוימות, כגון Transverse Mercator, Albers, ו
החרוט הקונפורמי של למברט אנו משתמשים בביטויים אליפסואידים כאשר האזורים הממופים הם
קטן, ועבור לביטויים הכדוריים (והחלפת העזר המתאים
קווי רוחב) עבור מפות גדולות יותר. הנוסחאות האליפסואידיות משמשות כדלקמן: (א) רוחבי
Mercator: כאשר כל הנקודות נמצאות בטווח של 10 מעלות ממרידיאן מרכזי, (ב) תחזיות חרוטיות
כאשר טווח האורך הוא פחות מ-90 מעלות, (ג) הקרנת Cassini כאשר כל הנקודות
נמצאים בטווח של 4 מעלות מהמרידיאן המרכזי. (2) כאשר אתה מנסה להתאים היסטורית כלשהי
נתונים (למשל, קואורדינטות שהושגו עם השלכה מסוימת והתייחסות מסוימת
אליפסואיד) אתה עשוי לגלות ש-GMT נותן תוצאות מעט שונות. אחד סביר
מקור אי ההתאמה הזה הוא שחישובים ישנים יותר השתמשו לעתים קרובות בספרות פחות משמעותיות. ל
לדוגמה, הדוגמאות של סניידר משתמשות לעתים קרובות באליפסואיד קלארק 1866 (שהוגדר על ידו כבעל
a השטחה f = 1/294.98). מ-f נקבל את האקסצנטריות בריבוע להיות 0.00676862818
(זה מה ש-GMT משתמש), בעוד סניידר מסיים ומשתמש ב-0.00676866. ההבדל הזה יכול
לתת פערים של כמה עשרות ס"מ. אם אתה צריך לשחזר קואורדינטות מוקרנות
עם האקסצנטריות המעט שונה הזו, עליך לציין את האליפסואיד שלך עם ה-
אותם פרמטרים כמו קלארק 1866, אבל עם f = 1/294.97861076. כמו כן, שים לב לכך שמבוגרים יותר
ניתן להפנות נתונים לדאטונים שונים, אלא אם כן אתה יודע באיזה נתון נעשה שימוש
המר את כל הנתונים לנתון משותף שאתה עלול לחוות אי התאמה של עשרות עד מאות
מטרים. (3) לבסוף, שים לב שלPROJ_SCALE_FACTOR יש ערכי ברירת מחדל מסוימים עבור חלקם
תחזיות כך שייתכן שתצטרך לעקוף את ההגדרה כדי להתאים לתוצאות שהופקו
עם הגדרות אחרות.

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