លំហាត់ទី១០៖ គណនា និងវាយតម្លៃសន្ទស្សន៍ស្ពិចត្រាល់¶
ព័ត៌មានលំហាត់
មេរៀនពាក់ព័ន្ធ៖ មេរៀនទី១០៖ សន្ទស្សន៍ស្ពិចត្រាល់ ការអនុវត្តដោយខ្លួនឯង · QGIS 3.34 LTR · ប្រហែល ១១០ នាទី
ស្ថានភាព¶
អ្នកនឹងគណនាសន្ទស្សន៍ចំនួនប្រាំ លើរូបភាព Sentinel-2 របស់អ្នក ប្រៀបធៀបនឹងគ្នា ហើយវាយតម្លៃថាសន្ទស្សន៍ណាបំពេញភារកិច្ចនីមួយៗបានល្អបំផុត។
គោលបំណង¶
- គណនា NDVI NDWI MNDWI NDBI និង BSI ក្នុង Raster Calculator។
- បង្កើតផែនទីពណ៌ (colour ramp) សម្រាប់សន្ទស្សន៍នីមួយៗ។
- យកគំរូតម្លៃលើថ្នាក់គម្របដីខុសៗគ្នា ហើយប្រៀបធៀបភាពបែងចែក។
- កំណត់សន្ទស្សន៍ណាដែលសមស្របនឹងភារកិច្ចនីមួយៗ។
ឯកសារដែលប្រើ¶
- ទិន្នន័យវគ្គ (ណែនាំ)៖ Landsat 8 ភ្នំពេញ (TOA ពីលំហាត់ទី៧) ·
RS_Data.zip→L8_PhnomPenh_2019(ឧបសម្ព័ន្ធ ខ) - ឬរូបភាព Sentinel-2 L2A ពីលំហាត់មុន (B2 B3 B4 B8 B11)
ឈ្មោះក្រុមរលកតាមឧបករណ៍
| តួនាទី | Sentinel-2 | Landsat 8/9 |
|---|---|---|
| ខៀវ · បៃតង · ក្រហម | B2 · B3 · B4 | B2 · B3 · B4 |
| NIR | B8 | B5 |
| SWIR1 · SWIR2 | B11 · B12 | B6 · B7 |
បើប្រើទិន្នន័យ Landsat 8 ភ្នំពេញ សូមប្ដូរលេខក្រុមរលកក្នុងរូបមន្ត និងបន្សំពណ៌តាមតារាងនេះ (ឧ. NDVI = (B5 − B4) / (B5 + B4))។
រូបមន្ត Landsat៖ NDVI (B5−B4)/(B5+B4) · NDWI (B3−B5)/(B3+B5) · MNDWI (B3−B6)/(B3+B6) · NDBI (B6−B5)/(B6+B5) · BSI ((B6+B4)−(B5+B2))/((B6+B4)+(B5+B2))
ពិសោធន៍មុនចាប់ផ្ដើម · ១៥ នាទី¶
សកម្មភាពទី១៖ គណនាសន្ទស្សន៍ប្រាំ · ៣៥ នាទី¶
ប្រើ Raster Calculator សម្រាប់នីមួយៗ (ជំនួស ឈ្មោះស្រទាប់ តាមតម្រូវការ)៖
| សន្ទស្សន៍ | កន្សោម |
|---|---|
| NDVI | ("B8@1" - "B4@1") / ("B8@1" + "B4@1") |
| NDWI | ("B3@1" - "B8@1") / ("B3@1" + "B8@1") |
| MNDWI | ("B3@1" - "B11@1") / ("B3@1" + "B11@1") |
| NDBI | ("B11@1" - "B8@1") / ("B11@1" + "B8@1") |
| BSI | (("B11@1"+"B4@1")-("B8@1"+"B2@1")) / (("B11@1"+"B4@1")+("B8@1"+"B2@1")) |
ប្រភេទទិន្នន័យលទ្ធផល
ក្នុងប្រអប់ Raster Calculator ត្រូវប្រាកដថាកំណត់ Output raster data type = Float32 មិនមែន Int ទេ បើមិនដូច្នោះទេ តម្លៃទសភាគនឹងត្រូវកាត់ចោល។
សកម្មភាពទី២៖ ដាក់ពណ៌ · ២០ នាទី¶
- សម្រាប់ NDVI៖ Singleband pseudocolor · Color ramp «RdYlGn» · Min
-1Max1។ - សម្រាប់ NDWI/MNDWI៖ Color ramp «Blues» · Min
-1Max1។ - សម្រាប់ NDBI/BSI៖ Color ramp «OrRd»។
- ថតរូបអេក្រង់សន្ទស្សន៍ទាំងប្រាំ។
សកម្មភាពទី៣៖ យកគំរូ និងប្រៀបធៀប · ២៥ នាទី¶
- ប្រើចំណុចគំរូពីលំហាត់ទី៣ (ឬបង្កើតថ្មី) សម្រាប់ថ្នាក់ប្រាំ។
- Sample raster values ទាញតម្លៃទាំងប្រាំសន្ទស្សន៍សម្រាប់ចំណុចនីមួយៗ។
-
គណនាមធ្យមតាមថ្នាក់ ហើយបំពេញតារាង៖
ថ្នាក់ NDVI NDWI MNDWI NDBI BSI ទឹក ព្រៃឈើ ស្រែ តំបន់សាងសង់ ដីទទេ -
សម្រាប់ថ្នាក់នីមួយៗ គូសរង្វង់លើតម្លៃខ្ពស់បំផុត (ដែលបញ្ជាក់ថ្នាក់នោះបានច្បាស់បំផុត)។
សកម្មភាពទី៤៖ វាយតម្លៃភាពបែងចែក · ១៥ នាទី¶
សំណួរ៖ ចម្លើយពីតារាងខាងលើ៖
- សន្ទស្សន៍ណាបែងចែកទឹកបានច្បាស់បំផុតពីថ្នាក់ដទៃ?
- សន្ទស្សន៍ណាបែងចែកតំបន់សាងសង់ពីដីទទេបានល្អបំផុត (ឬអត់អាចបែងចែកទាល់តែសោះ)?
- តើ NDVI តែម្នាក់ឯង គ្រប់គ្រាន់ដើម្បីបែងចែកព្រៃឈើ និងស្រែទេ?
សកម្មភាពទី៥៖ របាយការណ៍ខ្លី · ១៥ នាទី¶
សរសេរ ១ ទំព័រ៖ សន្ទស្សន៍ណាដែលអ្នកនឹងប្រើសម្រាប់ (ក) ផែនទីព្រៃឈើ (ខ) ផែនទីទឹកជំនន់ (គ) ការតាមដានការពង្រីកទីក្រុង ព្រមហេតុផលពីតារាងដែលអ្នកបានបំពេញ។
លទ្ធផលដែលរំពឹងទុក¶
ប្រៀបធៀបលទ្ធផលរបស់អ្នកជាមួយរូបខាងក្រោម ដែលបង្កើតពីទិន្នន័យវគ្គ (Landsat 8 ភ្នំពេញ ឬ Sentinel-2 ក្រុងព្រះសីហនុ)។ លទ្ធផលរបស់អ្នកមិនចាំបាច់ដូចបេះបិទទេ ប៉ុន្តែលំនាំគួរតែស្រដៀងគ្នា។
អ្វីដែលត្រូវពិនិត្យ
- NDVI បៃតងលើដើមឈើ និងដំណាំ ក្រហមលើទឹក។
- MNDWI បំបែកទឹកពីទីក្រុងបានល្អជាង NDWI។
- NDBI និង BSI ស្រដៀងគ្នា៖ សំណង់ និងដីទទេពិបាកបំបែក។
អ្វីដែលត្រូវពិនិត្យ
- ទឹក (អវិជ្ជមាន) និងដើមឈើ (ខ្ពស់) ដាច់ពីគ្នាល្អ។
- ដំណាំ/ស្មៅ ត្រួតនឹងដើមឈើមួយផ្នែក៖ NDVI តែមួយមិនគ្រប់គ្រាន់។
- សំណង់ និងដីទទេ មាន NDVI ជិតគ្នា (០–០,២)។
កំហុសទូទៅ និងដំណោះស្រាយ¶
| បញ្ហាដែលឃើញ | មូលហេតុទូទៅ | ដំណោះស្រាយ |
|---|---|---|
| NDVI ជា ០ ឬ ១ ទាំងអស់ | Output type Int | Float32 |
| NDVI > ១ | ប្រើ DN មាន offset | គណនាពីការចាំងផ្លាត |
| ទឹកមិនអវិជ្ជមាន | ទឹកល្អក់ ឬ sun glint | ប្រើ MNDWI |
| NDBI ច្រឡំដីទទេជាសំណង់ | សញ្ញាណស្រដៀងគ្នា | ប្រើរួមជាមួយ NDVI និងការបកស្រាយ |
លទ្ធផលត្រូវប្រគល់¶
- ស្រទាប់សន្ទស្សន៍ប្រាំ (.tif) ជាមួយពណ៌។
- តារាងតម្លៃមធ្យមតាមថ្នាក់ និងសន្ទស្សន៍។
- របាយការណ៍ខ្លី ១ ទំព័រ។
ពិនិត្យលទ្ធផលដោយខ្លួនឯង¶
អាចវាយលេខខ្មែរ ឬលេខអារ៉ាប់។
១. ប្រភេទទិន្នន័យលទ្ធផលដែលត្រូវកំណត់ក្នុង Raster Calculator?
២. ដែនតម្លៃទូទៅនៃសន្ទស្សន៍ normalized difference?
៣. ក្រុមរលកណាដែលប្រើក្នុងទាំង NDVI ទាំង NDBI?
៤. Color ramp ណាដែលសមស្របសម្រាប់ NDVI?
៥. សន្ទស្សន៍ណាប្រើ Green និង SWIR1?