7. ترسيم شبكة المجاري المائية

7.3. حساب شبكة المجاري المائية

الخطوة التالية هي حساب شبكة المجاري المائية (Calculate the channel network).
 
سنقوم أولاً بإنشاء طبقة بـ رتب ستراهلر (Strahler orders) للأنهار عبر تحديد الخلايا (Cells) التي تساوي رتبتها أو تزيد عن العتبة (Threshold) المحددة في الخطوة السابقة. سنستخدم هنا الرتبة 8 كعتبة.
 
باستخدام ملحق PCRaster Tools، لا يمكننا استخدام سوى طبقات الراستر (Raster) في جبر الخرائط (Map algebra). وإذا أردنا حساب طبقة منطقية (Boolean layer) جديدة تكون فيها القيم (strahler >= 8) هي الصواب (True) (أي قنواتنا المائية) والقيم (< 8) هي الخطأ (False)، فإننا نحتاج أولاً إلى إنشاء راستر ترتيبي (Ordinal raster) يتكون من خلايا ذات القيمة 8 فقط، ليمكن استخدامه لاحقاً في جبر الخرائط.
 
1.  في صندوق أدوات المعالجة (Processing Toolbox)، انتقل إلى المسار التالي: PCRaster | Data management | spatial.
 
2.  في مربع حوار Spatial، اكتب 8 في خانة Input nonspatial. اختر Ordinal لنوع بيانات المخرجات (Output data type) وملف strahler كـ طبقة القناع (Mask layer). قم بتسمية طبقة الراستر الناتجة باسم ordinal8.map.
 

Spatial Dialog to Create an Ordinal Raster with Value 8

3. انقر على تشغيل (Run) وأغلق مربع الحوار بعد اكتمال العملية.
 
سنقوم الآن بإنشاء طبقة منطقية (Boolean layer) تُعطي القيمة 1 (صواب/True) لكافة خلايا ستراهلر (Strahler cells) التي تكون قيمتها >= 8.
 
4. في صندوق أدوات المعالجة (Processing Toolbox)، انتقل إلى المسار التالي: PCRaster | Conditional and boolean operators | comparison operators.
 
5. في مربع حوار Comparison Operators، اختر strahler كراستر للمدخلات (Input raster)، والرمز >= كمعامِل مقارنة (Comparison operator)، وordinal8 كراستر المدخلات الثاني. قم بتسمية ملف الراستر المنطقي الناتج باسم channels.map.
 
 
 
6. انقر على تشغيل (Run) وأغلق مربع الحوار بعد اكتمال العملية.
 
أصبح لدينا الآن طبقة منطقية (Boolean layer) تضم القنوات المائية (Channels)، والتي يمكننا استخدامها الآن لـ تخصيص (Assign) رتب ستراهلر للأنهار (River Strahler orders).
 
7. في صندوق أدوات المعالجة (Processing Toolbox)، انتقل إلى المسار التالي: PCRaster | Conditional and boolean operators | ifthen.
 
8. في مربع حوار If Then، اختر channels كراستر لحالة الشرط المنطقي للمدخلات (Input Boolean Condition Raster)، وstrahler كراستر للقيمة "صواب" للمدخلات (Input True Raster). قم بتسمية راستر المخرجات باسم channelsstrahler.map.
 

وهذا يعني: إذا كانت طبقة القنوات تحتوي على خلايا قيمتها صواب (1)، فقم بإعطاء تلك الخلايا قيمة طبقة ستراهلر (strahler). أما كافة الخلايا الأخرى فستأخذ القيمة "لا توجد بيانات" (nodata).
 
9. انقر على  تشغيل (Run)  وأغلق مربع الحوار بعد اكتمال المعالجة (Processing).
 
شاهد هذا الفيديو للتحقق من الخطوات:


أصبح لدينا الآن طبقة راستر (Raster) تحتوي على رتب ستراهلر (Strahler orders) للقنوات المائية. يمكن تمثيل هذه البيانات بصرياً بشكل أفضل في صيغة فكتور (Vector)، وهو ما سنقوم بشرحه في القسم التالي.