أكثر

الحصول على قيم ميل غير صحيحة من ASTER DEM في QGIS؟

الحصول على قيم ميل غير صحيحة من ASTER DEM في QGIS؟


لدي طبقة DEM من القمر الصناعي ASTER وأريد الحصول على معلومات المنحدر منه.

أستخدم تحليل التضاريس> المنحدر ، لكن الناتج خاطئ جدًا. لماذا هذا؟

هنا DEM:

وهنا "المنحدر" هذا خاطئ تمامًا. يجب أن تمدد القيم من 0 - 89.9097 لكنها لا تفعل ذلك.


Aster DEMs موجودة في EPSG 4326 (خط الطول WGS84). وفقًا لوثائق GDAL DEM:

بالنسبة لإسقاطات LatLong بالقرب من خط الاستواء ، حيث تتشابه وحدات خط العرض ووحدات خط الطول ، يمكن تحويل وحدات الارتفاع (z) لتكون متوافقة باستخدام المقياس = 370400 (إذا كان الارتفاع بالأقدام) أو المقياس = 111120 (إذا كان الارتفاع في أمتار). بالنسبة للمواقع غير القريبة من خط الاستواء ، سيكون من الأفضل إعادة تقسيم شبكتك باستخدامغدالواربقبل استعمالgdaldem.

في هذه الإجابة لمشكلة تتعلق بتظليل التلال (لا يؤثر عامل المقياس و Z على تحليل الظل في QGIS) ، يوفر المستخدمLuke طريقة لحساب المقياس الصحيح ("نسبة الوحدات الرأسية إلى الأفقية") بضرب جيب التمام لخط العرض ضرب 111320. يؤثر المقياس على حسابات الميل.

في الصورة أدناه ، تمكنت من إنتاج منحدرات في درجات من DEM في EPSG 4326 عند خطوط عرض تتراوح من 19 درجة 45 شمالاً و 20 درجة شمالاً. تم ضبط المقياس (-s) على104550.67.

إذا كنت ترغب في العمل مع DEM في CRS الذي يحدد وحدات قياس الإحداثيات والارتفاعات بالأمتار ، يجب عليك إعادة إسقاط DEM باستخدام البيانات النقطية / الإسقاطات / الالتواء (إعادة الإسقاط). على سبيل المثال ، إعادة الإسقاط من EPSG 4326 (خط الطول wgs84) إلى EPSG 32164 (UTM WGS84 14N). انظر أدناه:

في هذه الحالة ، عند حساب الميل ، سيتم تعيين "نسبة الوحدات الرأسية إلى الأفقية" (المقياس) على 1.


ربما لأن الإسقاط لم يكن بالأمتار.

المنحدر مع الإسقاط DEM بالأمتار (UTM / 31N: EPSG 25831):

المنحدر بنفس المارك الألماني المسقط بالدرجات (WGS 84: EPSG 4326) و z = 1:

المنحدر بنفس DEM المسقط بالدرجات (WGS 84: EPSG 4326) و z = 111320 (عند خط الاستواء ، وحدة أفقية واحدة ، 1 ° ، تتوافق مع حوالي 111320 م):


شاهد الفيديو: How to clip Dem Easily By QGIS