r.watershedgrass - מקוון בענן

תָכְנִית:

שֵׁם

r.פרשת מים - מחשבת פרמטרים הידרולוגיים וגורמי RUSLE.

מילות מפתח

רסטר, הידרולוגיה, פרשת מים

תַקצִיר

r.פרשת מים
r.פרשת מים - עזרה
r.פרשת מים [-s4mab] גוֹבַה=שם [דכאון=שם] [תזרים=שם]
[ארץ_מופרעת=שם] [חסימה=שם] [סף=מספר שלם] [max_slope_length=לצוף]
[הצטברות=שם] [Tci=שם] [תעלת ניקוז=שם] [אַגָן=שם] [זרם=שם]
[חצי_אגן=שם] [אורך_שיפוע=שם] [שיפוע_תלול=שם] [התכנסות=מספר שלם]
[זיכרון=מספר שלם] [--להחליף] [--לעזור] [--מִלוּלִי] [--שֶׁקֶט] [--ui]

דגלים:
-s
זרימת SFD (D8) (ברירת המחדל היא MFD)
SFD: כיוון זרימה בודדת, MFD: כיוון זרימה מרובה

-4
אפשר רק זרימה אופקית ואנכית של מים

-m
אפשר אפשרות זיכרון החלפת דיסק: הפעולה איטית
נדרש רק אם דרישות הזיכרון חורגות מזיכרון RAM זמין; ראה מדריך כיצד לעשות זאת
לחשב את דרישות הזיכרון

-a
השתמש בהצטברות זרימה חיובית אפילו להערכות חסרות סבירות
ראה מדריך לתיאור מפורט של תפוקת צבירת הזרימה

-b
לייפות אזורים שטוחים
כיוון הזרימה באזורים שטוחים שונה כדי להיראות יפה יותר

--לכתוב
אפשר לקבצי פלט לדרוס קבצים קיימים

- עזרה
הדפס סיכום שימוש

--מִלוּלִי
פלט מודול מפורט

--שֶׁקֶט
פלט מודול שקט

--ui
הפעלת דו-שיח של ממשק משתמש לאלץ

פרמטרים:
גוֹבַה=שם [נדרש]
שם מפת רסטר גובה הקלט

דכאון=שם
שם מפת רסטר של שקעי קלט
כל התאים שאינם NULL ולא אפס נחשבים לשקעים אמיתיים

תזרים=שם
שם רסטר הקלט המייצג את כמות הזרימה היבשתית לתא

ארץ_מופרעת=שם
שם אחוז מפת רסטר הקלט מהקרקע המופרעת
עבור USLE

חסימה=שם
שם מפת רסטר קלט החוסמת זרימה יבשתית
עבור USLE. כל התאים שאינם NULL ולא אפס נחשבים כשטח חוסם.

סף=מספר שלם
גודל מינימלי של אגן פרשת מים חיצוני

max_slope_length=לצוף
אורך מקסימלי של זרימת פני השטח ביחידות מפה
עבור USLE

הצטברות=שם
שם למפת רסטר של צבירת פלט
מספר התאים שמתנקזים בכל תא

Tci=שם
שם עבור פלט אינדקס טופוגרפי מפת ln(a / tan(b)).

תעלת ניקוז=שם
שם עבור מפת רסטר של כיוון ניקוז הפלט

אַגָן=שם
שם עבור מפת רסטר של אגני פלט

זרם=שם
שם עבור מפת רסטר של קטעי זרם פלט

חצי_אגן=שם
שם עבור מפת רסטר של חצאי אגנים פלט
לכל חצי אגן ניתן ערך ייחודי

אורך_שיפוע=שם
שם עבור מפת רסטר באורך מדרון פלט
גורם אורך מדרון ותלול (LS) עבור USLE

שיפוע_תלול=שם
שם עבור מפת רסטר תלילות שיפוע פלט
גורם תלילות שיפוע (S) עבור USLE

התכנסות=מספר שלם
גורם התכנסות עבור MFD (1-10)
1 = הזרימה המתפצלת ביותר, 10 = הזרימה המתכנסת ביותר. מומלץ: 5
ברירת מחדל: 5

זיכרון=מספר שלם
זיכרון מקסימלי לשימוש עם דגל -m (ב-MB)
ברירת מחדל: 300

תיאור

r.פרשת מים יוצר קבוצה של מפות המציינת: 1) הצטברות זרימה, כיוון ניקוז,
מיקום הנחלים ואגני פרשת המים, ו-2) גורמי ה-LS וה-S של ה-Revised
משוואת אובדן קרקע אוניברסלית (RUSLE).

אורים

בלי דגל -m מוגדר, כל הניתוח מתנהל בזיכרון המתוחזק על ידי ההפעלה
מערכת. זה יכול להיות מגביל, אבל הוא מהיר מאוד. הגדרת דגל זה גורמת לתוכנית
ניהול זיכרון בדיסק המאפשר עיבוד מפות גדולות מאוד אך איטי יותר.

דֶגֶל -s לאלץ את המודול להשתמש בכיוון זרימה בודדת (SFD, D8) במקום זרימה מרובה
כיוון (MFD). MFD מופעל כברירת מחדל.

By -4 סמן את המשתמש לאפשר רק זרימה אופקית ואנכית של מים. נחל ושיפוע
האורכים זהים בערך לתפוקות מזרימת השטח המוגדרת כברירת מחדל (מאפשר
זרימת מים אופקית, אנכית ואלכסונית). הדגל הזה גם יעשה את הניקוז
אגנים נראים הומוגניים יותר.

מתי -a דגל צוין שהמודול ישתמש בהצטברות זרימה חיובית אפילו בסבירות גבוהה
מזלזל. כאשר דגל זה אינו מוגדר, תאים עם ערך צבירת זרימה כלומר
שסביר שיהיו הערכת חסר מומרים לשלילה. ראה למטה למידע מפורט
תיאור של תפוקת צבירת זרימה.

אפשרות התכנסות מציין גורם התכנסות עבור MFD. ערכים נמוכים יותר מביאים לעלייה
סטייה, הזרימה מופצת יותר. ערכים גבוהים יותר מביאים להתכנסות גבוהה יותר,
הזרימה מפוזרת פחות, והופכת דומה יותר ל-SFD.

אפשרות גוֹבַה מפרט את נתוני הגובה שעליהם מבוסס הניתוח כולו. ריק
מתעלמים מתאי (nodata), ערכים אפס וערכים שליליים הם נתוני גובה חוקיים. פערים פנימה
יש למלא לפני כן את מפת הגבהים הנמצאים בתחום העניין,
למשל עם R.fillnulls, כדי למנוע עיוותים. אין צורך למלא את מפת הגובה בכיור
מכיוון שהמודול משתמש באלגוריתם בעל עלות נמוכה.

אפשרות דכאון מציין את המפה האופציונלית של שקעים או בולענים בפועל
נוף גדול מספיק כדי להאט ולאגור נגר עילי מאירוע סערה. את כל
תאים שאינם NULL ולא אפס מציינים שקעים. מים יזרמו לתוך אבל לא
מתוך דיכאונות.

סריקה תזרים המפה מציינת את כמות הזרימה היבשתית לתא. מפה זו מציינת את הכמות
של יחידות זרימה יבשתית שכל תא יתרום למודל אגן פרשת המים.
יחידות זרימה יבשתית מייצגות את כמות הזרימה היבשתית שכל תא תורם לפני השטח
זְרִימָה. אם מושמט, מניחים ערך של אחד (1).

קלט רסטר מפה או ערך המכילים את אחוז הקרקע המופרעת (כלומר, אדמות גידול ו
אתרי בנייה) שבהם הרסטר או ערך הקלט של 17 שווה ל-17%. אם אין מפה או ערך
נתון, r.פרשת מים מניח שאין אדמה מופרעת. קלט זה משמש עבור RUSLE
חישובים.

אפשרות חסימה מפרט את השטח שיחסום את הזרימה היבשתית. לבטל את זה
חסום את זרימת השטח היבשתי והפעל מחדש את אורך המדרון עבור ה-RUSLE. כל התאים ש
אינם NULL ולא אפס מעידים על חסימת שטח.

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

ערך שניתן על ידי max_slope_length אפשרות מציינת את האורך המרבי של פני השטח היבשתיים
זרימה במטרים. אם הזרימה היבשתית עוברת יותר מהאורך המרבי, התוכנית
מניח את האורך המקסימלי (הוא מניח שמאפייני נוף לא ניתן להבחין ב
קיים מודל הגובה הדיגיטלי שממקסם את אורך השיפוע). קלט זה משמש עבור
חישובי RUSLE והוא פרמטר רגיש.

תְפוּקָה הצטברות מפה מכילה את הערך המוחלט של כל תא במפת הפלט הזו היא
כמות זרימה יבשתית שחוצה את התא. ערך זה יהיה המספר של למעלה
תאים פלוס אחד אם לא ניתנת מפת זרימה יבשתית. אם ניתנת מפת הזרימה היבשתית, ה
הערך יהיה ביחידות זרימה יבשתית. מספרים שליליים מצביעים על כך שאותם תאים אולי
יש נגר עילי מחוץ לאזור הגיאוגרפי הנוכחי. לפיכך, כל התאים עם
ערכים שליליים אינם יכולים לחשב את תשואות הנגר והשיקוע שלהם
במדויק.

תְפוּקָה Tci מפת רסטר מכילה אינדקס טופוגרפי. TCI מחושב כ-ln(α /
tan(β)) כאשר α הוא שטח המדרון המצטבר המתנקז דרך נקודה ליחידה
אורך קו המתאר ו-tan(β) היא זווית השיפוע המקומית. ה-TCI משקף את הנטייה של
מים להצטבר בכל נקודה בתפיסה והנטייה לכוחות כבידה
להזיז את מדרון המים הזה (Quinn et al. 1991). ערך זה יהיה שלילי אם α /
tan(β) < 1.

תְפוּקָה תעלת ניקוז מפת רסטר מכילה כיוון ניקוז. מספק את ה"היבט" לכל אחד
תא נמדד CCW ממזרח. הכפלת ערכים חיוביים ב-45 תיתן את הכיוון פנימה
מעלות שהנגר העילי יעבור מאותו תא. הערך 0 (אפס) מציין
שהתא הוא אזור שקע (מוגדר על ידי מפת הקלט של שקע). ערכים שליליים
מצביעים על כך שנגר עילי עוזב את גבולות האזור הגיאוגרפי הנוכחי.
הערך המוחלט של תאים שליליים אלה מציין את כיוון הזרימה.

הפלט אַגָן המפה מכילה תווית ייחודית לכל אגן פרשת מים. כל אגן יהיה
נתון מספר שלם זוגי חיובי ייחודי. ייתכן ששטחים לאורך הקצוות לא יהיו גדולים מספיק כדי ליצור
אגן פרשת מים חיצוני. ערכי 0 מציינים שהתא אינו חלק מהשלם
אגן פרשת המים באזור הגיאוגרפי הנוכחי.

הפלט זרם מכיל קטעי זרם. הערכים תואמים לאגן פרשת המים
ערכים. ניתן לווקטור לאחר דילול (ר.דק) עם r.to.vect.

הפלט חצי_אגן מפות רסטר מאחסנות כל חצי אגן מקבל ערך ייחודי. פָּרָשַׁת הַמַיִם
אגנים מחולקים לצד שמאל וימין. התא הימני של קו פרשת המים
לאגן (במבט במעלה הזרם) ניתנים ערכים זוגיים התואמים לערכים באגן. ה
תאים בצד שמאל של אגן פרשת המים מקבלים ערכים אי-זוגיים שהם אחד פחות מ
ערכו של אגן פרשת המים.

הפלט אורך_שיפוע מפת רסטר מאחסנת גורם אורך שיפוע ותלול (LS) עבור
משוואת אובדן קרקע אוניברסלית מתוקנת (RUSLE). משוואות שנלקחו מ מתוקן אוניברסלי קרקע
פסד משוואה ל מערבי ריינג'לס (Weltz et al. 1987). מכיוון שגורם LS הוא קטן
מספר (בדרך כלל פחות מאחד), מפת הפלט של GRASS היא מסוג DCELL.

הפלט שיפוע_תלול מפת רסטר מאחסנת את גורם תלילות השיפוע (S) עבור ה-Universal
משוואת אובדן קרקע (RUSLE). משוואות שנלקחו מהמאמר בשם מתוקן מדרון תְלִילוּת
גורם ל מה היא אוניברסלי קרקע פסד משוואה (McCool et al. 1987). מכיוון שגורם S הוא a
מספר קטן (בדרך כלל פחות מאחד), מפת הפלט של GRASS היא מסוג DCELL.

AT העלות הנמוכה ביותר search אַלגוֹרִיתְם
r.פרשת מים משתמש באלגוריתם חיפוש בעלות לפחות (ראה סעיף הפניות) המיועד לכך
למזער את ההשפעה של שגיאות בנתוני DEM. לְעוּמַת r.terraflow, מספק האלגוריתם הזה
תוצאות מדויקות יותר באזורים עם שיפוע נמוך כמו גם DEMs שנבנו בטכניקות
שגגות חופה זו כמו גובה הקרקע. כינר וחב'. (2005), למשל, השתמשו
SRTM ו-IFSAR DEMs להשוואה r.פרשת מים נגד r.terraflow תוצאות בפנמה.
r.terraflow לא הצליח לשחזר מיקומי נחלים בעמקים הגדולים יותר בזמן
r.פרשת מים פעל הרבה יותר טוב. לפיכך, אם קיימת חופת יער בעמקים, SRTM, IFSAR,
ומוצרי נתונים דומים יגרמו לשגיאות גדולות ב r.terraflow פלט זרם. תַחַת
תנאים דומים, r.פרשת מים יפיק טוב יותר זרם ו חצי_אגן תוצאות. אם
מחלקי פרשת מים מכילים שיפוע שטוח עד נמוך, r.פרשת מים יפיק אגן טוב יותר
תוצאות מאשר r.terraflow. (r.terraflow משתמש באותו סוג של אלגוריתם כמו ArcGIS של ESRI
תוכנת פרשת מים שנכשלת בתנאים אלה.) כמו כן, אם פרשת מים מחלקים מכילים
חופת יער מעורבת עם אזורים לא חטופים באמצעות SRTM, IFSAR ומוצרי נתונים דומים,
r.פרשת מים יפיק תוצאות אגן טובות יותר מאשר r.terraflow. האלגוריתם מייצר
תוצאות דומות לאלו המתקבלות בריצה r.cost ו r.drain על כל תא ב-
מפת רסטר.

מְרוּבֶּה תזרים כיוון (MFD)
r.פרשת מים מציעה שתי שיטות לחישוב זרימת פני השטח: כיוון זרימה יחיד (SFD, D8)
וכיוון זרימה מרובה (MFD). עם MFD, זרימת המים מופצת לכל השכנים
תאים בעלי גובה נמוך יותר, תוך שימוש בשיפוע לכיוון תאים שכנים כגורם שקילה
לחלוקה פרופורציונלית. נתיב העלות המינימלית כלול תמיד. כתוצאה,
שקעים ומכשולים עוברים בהתכנסות זרימה חיננית לפני ה
הצפה. גורם ההתכנסות גורם להצטברות זרימה להתכנס חזק יותר עם
ערכים גבוהים יותר. הטווח הנתמך הוא 1 עד 10, מומלץ הוא גורם התכנסות של 5
(Holmgren, 1994). אם נוצרים אגני רסיס קטנים רבים עם MFD, הגדר את
גורם התכנסות לערך גבוה יותר יכול להפחית את כמות אגני הרסיס הקטנים.

לזכר מצב ו דיסק להחליף מצב
ישנן שתי גרסאות של תוכנית זו: RAM ו seg. RAM משמש כברירת מחדל, seg יכול להיות
בשימוש על ידי הגדרת ה -m דגל.

השמיים RAM הגרסה דורשת מקסימום 31 MB של זיכרון RAM עבור מיליון תאים. ביחד עם ה
כמות זיכרון המערכת (RAM) זמין, ניתן להשתמש בערך זה כדי להעריך אם
ניתן לעבד את האזור הנוכחי עם RAM הגירסה.

השמיים RAM הגרסה משתמשת בזיכרון וירטואלי המנוהל על ידי מערכת ההפעלה כדי לאחסן את כל הנתונים
מבנים ומהיר יותר מה seg גִרְסָה; seg משתמש בספריית הפילוח GRASS
המנהלת נתונים בקבצי דיסק. seg משתמש רק בזיכרון המערכת (RAM) ככל שצוין
עם זיכרון אפשרות, המאפשרת לתהליכים אחרים לפעול על אותה מערכת, גם כאשר
האזור הגיאוגרפי הנוכחי הוא עצום.

בשל דרישות הזיכרון של שתי התוכניות, די קל להיגמר הזיכרון כאשר
עבודה עם אזורי מפה ענקיים. אם ה RAM בגרסה נגמר הזיכרון והרזולוציה
לא ניתן להגדיל את גודל האזור הגיאוגרפי הנוכחי, גם אם צריך יותר זיכרון
נוסף למחשב, או שיש להגדיל את גודל שטח ההחלפה. אם seg נגמר
זיכרון, יש לפנות שטח דיסק נוסף כדי שהתוכנית תפעל. ה
r.terraflow המודול תוכנן במיוחד עם מחשבה על אזורים ענקיים ועשוי להיות שימושי
כאן כחלופה, אם כי דרישות שטח דיסק של r.terraflow הם מספר פעמים
גבוה יותר מ seg.

גדול אזורים עם רב תאים
הגבול העליון של RAM הגרסה היא 2 מיליארד (231 - 1) תאים, ואילו הגבול העליון
עבור seg הגרסה היא 9 מיליארד מיליארד (263 - 1 = 9.223372e+18) תאים.
במצבים מסוימים, גודל האזור (מספר התאים) עשוי להיות גדול מדי עבור הכמות של
זמן או זיכרון זמינים. רץ r.פרשת מים לאחר מכן עשוי לדרוש שימוש בחומר גס יותר
פתרון הבעיה. כדי שהתוצאות יהיו דומות יותר לנתוני השטח העדינים יותר, צור מפה
שכבה המכילה את ערכי הגובה הנמוכים ביותר ברזולוציה הגסה יותר. זה נעשה על ידי:
1) הגדרת האזור הגיאוגרפי הנוכחי שווה לשכבת מפת הגובה עם G.Region,
ו-2) השתמש ב- r.שכנים or r.resamp.stats פקודה כדי למצוא את הערך הנמוך ביותר עבור אזור
שווה בגודל לרזולוציה הרצויה. לדוגמה, אם הרזולוציה של הגובה
הנתונים הם 30 מטר והרזולוציה של האזור הגיאוגרפי עבור r.פרשת מים יהיה בן 90
מטר: השתמש בפונקציית המינימום עבור שכונה של 3 על 3. לאחר השינוי ל-
רזולוציה שבה r.פרשת מים יופעל, r.פרשת מים יש להפעיל באמצעות הערכים
מ שכונה שכבת מפת פלט המייצגת את הגובה המינימלי בתוך
אזור של התא הגס יותר.

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

מאסק ו לא נתונים
השמיים r.פרשת מים התוכנית אינה מחייבת את המשתמש לקבל את האזור הגיאוגרפי הנוכחי
מלא בערכי גובה. אזורים ללא נתוני גובה (מסיכות או תאים NULL) הם
התעלמו. אין צורך ליצור מפת רסטר (או סיווג רסטר מחדש) בשם
MASK עבור תאי NULL. אזורים ללא נתוני גובה יטופלו כאילו הם מחוץ ל-
קצה האזור. אזורים כאלה יצמצמו את הזיכרון הדרוש להפעלת התוכנית.
מיסוך אזורים לא חשובים יכול להפחית משמעותית את זמן העיבוד אם פרשת המים
העניין תופס אחוז קטן מהשטח הכולל.

פערים (תאי NULL) במפת הגובה הממוקמים בתוך אזור העניין יעשו זאת
משפיעים מאוד על הניתוח: מים יזרמו לתוך הפערים הללו אך לא יצאו מהם. הפערים האלה
יש למלא מראש, למשל עם R.fillnulls.

ערכים אפס (0) ושליליים יטופלו כנתוני גובה (לא no_data).

עוד תהליך of תפוקה שכבות
אזורים בעייתיים, כלומר אותם חלקים של אגן עם סבירות נמוכה של הערכת זרימה
הצטברות, ניתן לזהות בקלות עם למשל
r.mapcalc "problems = if(flow_acc < 0, basin, null())"
אם אזור העניין מכיל אזורים בעייתיים כאלה, וזה לא רצוי, ה
יש להרחיב את האזור החישובי עד לאזור התפיסה של אזור העניין
כלול לחלוטין.

כדי לבודד רשת נהר בודדת באמצעות הפלט של מודול זה, מספר
ניתן לשקול גישות.

1 השתמש בדגימה מחודשת של מפת הרסטר של אגנים כמסכה.
שיטת המפה הווקטורית המקבילה דומה לשימוש v.select or v.overlay.

2 השתמש r.cost מודול עם נקודה בנהר כנקודת התחלה.

3 השתמש v.net.iso מודול עם צומת בנהר כנקודת התחלה.

ניתן לזהות דרכן את כל רשתות הנהרות הבודדות בפלט מקטעי הנחל
נקודות היציאה האולטימטיביות שלהם. נקודות אלו הן כל התאים בפלט מקטעי הזרם
עם כיוון ניקוז שלילי. נקודות אלו יכולות לשמש כנקודות התחלה עבור
r.outlet.water or v.net.iso.

ליצור נהר מיל פילוח ממפת זרמים וקטורית, נסה את v.net.iso or
v.lrs.segment מודולים.

ניתן לשנות בקלות את פלט מקטעי הזרם לאחר דילול ר.דק. כל
קטע הנחל במפה הווקטורית יקבל את הערך של האגן המשויך. לבודד
תת אגנים וזרמים עבור אגן גדול יותר, ניתן ליצור MASK עבור האגן הגדול יותר
r.outlet.water. פלט מקטעי הנחל משמש כמדריך היכן למקם את השקע
נקודה המשמשת כקלט ל r.outlet.water. סף האגן חייב להיות מספיק
קטן כדי לבודד רשת נחל ותתי אגנים בתוך האגן הגדול יותר.

דוגמאות

דוגמאות אלה משתמשות במערך הנתונים לדוגמה של Spearfish.

המרת r.פרשת מים מזרים פלט מפה למפה וקטורית
אם אתה רוצה רשת זרם מפורטת, הגדר את אפשרות הסף קטנה כדי ליצור הרבה
אגני היבול, שכן רק נחל אחד מוצג לכל נחל. דגל r.to.vect -v
שומר על מזהה התפיסה כמספר הקטגוריה הווקטורית.
r.watershed elev=elevation.dem stream=rwater.stream
r.to.vect -v in=rwater.stream out=rwater_stream

הגדר טבלת צבעים שונה עבור מפת הצבירה:
MAP=rwater.accum
r.watershed elev=elevation.dem accum=$MAP
eval `r.univar -g "$MAP"`
stddev_x_2=`echo $stddev | awk '{print $1 * 2}'`
stddev_div_2=`echo $stddev | awk '{print $1 / 2}'`
r.colors $MAP col=rules << EOF
0% אדום
-$stddev_x_2 אדום
-$stddev צהוב
-$stddev_div_2 ציאן
-$mean_of_abs כחול
0 לבן
$mean_of_abs כחול
$stddev_div_2 ציאן
$stddev צהוב
$stddev_x_2 אדום
100% אדום
EOF

צור מפת זרמים מפורטת יותר באמצעות מפת הצבירה והמר אותה לוקטור
מפת פלט. חתך הצבירה, ולכן הממד הפרקטלי, הוא שרירותי; ב
בדוגמה זו אנו משתמשים במספר הממוצע של המפה של תאי התפיסה במעלה הזרם (מחושב ב-
דוגמה לעיל מאת r.univar) כערך החתך. זה עובד רק עם SFD, לא עם MFD.
r.watershed elev=elevation.dem accum=rwater.accum
r.mapcalc 'MASK = if(!isnull(elevation.dem))'
r.mapcalc "rwater.course =
if( abs(rwater.accum) > $mean_of_abs,
abs(rwater.accum),
ריק() )"
r.colors -g rwater.course col=bcyr
g.remove -f type=שם רסטר=MASK
# מדלל is נדרש לפני המרת סריקה קווים ל וקטור
r.thhin in=rwater.course out=rwater.course.Thin
r.colors -gn rwater.course.צבע דק=אפור
r.to.vect in=rwater.course.Thin out=rwater_course type=line
v.db.dropcolumn map=rwater_course column=label

צור קו פרשת המים אגנים מַפָּה ו להמיר ל a וקטור מצולע מַפָּה
r.watershed elev=elevation.dem basin=rwater.basin thresh=15000
r.to.vect -s in=rwater.basin out=rwater_basins type=area
v.db.dropcolumn map=rwater_basins column=label
v.db.renamecolumn map=rwater_basins column=value,catchment

הצג פלט בצורה יפה
r.relief map=elevation.dem
d.shade shade=elevation.dem.shade color=rwater.basin בהיר=40
d.vect rwater_course color=orange

ביבליוגרפיה

· Ehlschlaeger C. (1989). שימוש מה היא AT חיפוש אַלגוֹרִיתְם ל לפתח הידרולוגי מודלים
החל מ- דִיגִיטָלי גובה נתונים, הליכים of ברמה בינלאומית גיאוגרפי מֵידָע
מערכות (IGIS) סִימפּוֹזִיוֹן "89, עמ' 275-281 (בולטימור, MD, 18-19 במרץ 1989).
כתובת האתר: http://chuck.ehlschlaeger.info/older/IGIS/paper.html

· Holmgren P. (1994). מְרוּבֶּה תזרים כיוון אלגוריתמים ל נגר דוּגמָנוּת in
רשת מבוסס גוֹבַה הדגמה An אֶמפִּירִי הערכה. הידרולוגי תהליכים כרך
8(4), 327-334.
DOI: 10.1002/hyp.3360080405

· Kinner D., Mitasova H., Harmon R., Toma L., Stallard R. (2005). מבוסס GIS זרם
רשת אָנָלִיזָה ל השמיים צ'אגרס נהר אַגָן, רפובליקה of פנמה. השמיים ריו Chagres:
A רב תחומי פּרוֹפִיל of a טְרוֹפִּי פָּרָשַׁת הַמַיִם, ר' הרמון (עורך),
שפרינגר/קלוור, עמ' 83-95.
כתובת האתר: http://www4.ncsu.edu/~hmitaso/measwork/panama/panama.html

· McCool et al. (1987). מתוקן מדרון תְלִילוּת גורם ל מה היא אוניברסלי קרקע פסד
משוואה, עסקות of מה היא ASAE כרך 30(5).

· Metz M., Mitasova H., Harmon R. (2011). יָעִיל הוֹצָאָה of תעלת ניקוז רשתות
החל מ- מַסִיבִי, מבוסס מכ"ם גוֹבַה מודלים עם הכי פחות עלות נתיב search, הידרול.
כדור הארץ סיסט. Sci. כרך 15, 667-678.
DOI: 10.5194/hess-15-667-2011

· Quinn P., K. Beven K., Chevallier P., Planchon O. (1991). השמיים נבואה of
מדרון גבעה תזרים שבילים ל מופץ הידרולוגית דוּגמָנוּת באמצעות דִיגִיטָלי
גובה מודלים, הידרולוגי תהליכים כרך 5(1), עמ' 59-79.
DOI: 10.1002/hyp.3360050106

· Weltz MA, Renard KG, Simanton JR (1987). מתוקן אוניברסלי קרקע פסד
משוואה ל מערבי ריינג'לס, ארה"ב/מקסיקו סִימפּוֹזִיוֹן of אסטרטגיות ל
מִיוּן ו ניהול שוטף of יליד צמחייה ל מזון הפקה In צָחִיחַ
אזורים (טוסון, AZ, 12-16 באוקטובר 1987).

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