การวิเคราะห์พื้นที่เหมาะสมสำหรับการสร้างที่ฝังกลบขยะ
กำหนดประเด็นปัญหา
1. ความเหมาะสมในแต่ละระดับมีพื้นที่เท่าใด
2. ในแต่ละพื้นที่
(อำเภอ) มีความเหมาะสมอยู่ในระดับใดบ้าง
กำหนดปัจจัยหรือตัวแปรที่ใช้วิเคราะห์
1. ความลาดชัน
2. ลักษณะชั้นหินฐาน
3. ระยะห่างจากถนนสายหลัก
4. ระยะห่างจากชุมชน
การเตรียมข้อมูล
ได้แก่
การนำเข้าข้อมูลจุดความสูง การสร้างแบบจำลองความสูงเชิงเลข (DEM) การคำนวณความลาดชัน
การคำนวณ
Euclidean
distance
การจัดกลุ่มข้อมูลใหม่ (Reclassify)
ค่าคะแนนและค่าถ่วงน้ำหนักของแต่ละตัวแปร
การวิเคราะห์ข้อมูล
การวิเคราะห์ซ้อนทับ
(Overlay analysis)
Suitable
land
fill = W1*S1 + W2*S2 + …
+ Wn*Sn
จัดกลุ่มพื้นที่เหมาะสม
จัดกลุ่มพื้นที่เหมาะสมสำหรับการสร้างที่ฝังกลบขยะออกเป็น
3
กลุ่ม
ได้แก่
เหมาะสมน้อย X –
SD
เหมาะสมปานกลาง X - SD
≤ S ≤ X +
SD
เหมาะสมมาก
เปิดข้อมูล Chon_elv ที่เป็นข้อมูลแสดงค่าความสูง และข้อมูล Chon_pro ที่เป็นข้อมูลจังหวัด
ทำการประมาณค่าช่วงด้วยข้อมูลจุด หรือ IDW โดยเรียกคำสั่งมาจาก ArcToolbox > Spatial Analyst Tools > Interpolation > IDW
หน้าต่างของ IDW ช่องที่ 1 Input point features ใส่ข้อมูลที่แสคงค่าความสูง Chon_elv
ช่องที่ 2 Z value field เลือกเป็น ele
ช่องที่ 3 Output raster เลือกที่ที่จะจัดเก็บ
ช่องที่ 4 Output cell size กำหนดค่าเป็น 40 จากนั้นทำการตั้งค่า คลิก Environments...
เลือกตั้งค่าที่ Processing Extent ช่อง Extent เลือก Same as layer Chon_pro จากนั้นเลือกมาด้านล่างเพื่อทำการตั้งค่า Raster Analysis
Raster Analysis ในช่อง Mask เลือก Chon_pro จากนั้นคลิก OK
ตรวจสอบความถูกต้องของข้อมูล จากนั้นกด OK
จากนั้นจะได้ผลลัพธ์ขากการทำ IDW
ทำการหาค่า Slope โดยเลือกจากเครื่องมือ ArcToolbox > Spatial Analyst Tools > Surface > Slope
หน้าต่างของ Slope ช่องที่ 1 Input raster ให้เลือกชั้นข้อมูลที่ต้องการหาค่า
ช่องที่ 2 Output raster เลือกว่าจะจัดเก็บที่ไหน
ช่องที่ 3 Output measurrement เลือกเป็น Degree หรือแบบองศา จากนั้นคลิก OK
ผลลัพธ์ที่ได้จากการทำ Slope
ต่อมาคือการใช้คำสั่ง Reclassify โดยเลือกจาก ArcToolbox > Spatial Analyst Tools > Reclass > Reclassify
ช่อง Input เลือกข้อมูลที่ต้องการทำ Reclassify จากนั้นเข้าไปตั้งค่า โดยคลิกที่ Classify
ช่อง Classes ให้เลือก Class เป็น 4 และทำการเปลี่ยนค่าที่ Breale Values ตามลำดับคือ 5 , 15 , 30 และค่าสูงสุดในชั้นข้อมูล จากนั้นคลิก OK
เปลี่ยนค่าในช่อง New Values คือ 3 , 2 , 1 , 0 ตามลำดับ และในช่อง Output raster เลือกที่ ที่ต้องการจัดเก็บ จากนั้นคลิก OK
จะได้ข้อมูล Slope ที่ทำ Reclassify เรียบร้อยแล้ว
ต่อมาคือการ ทำชั้นข้อมูลหิน
ทำการเปิดชั้นข้อมูลหินขึ้นมา Chon_geo
ชั้นข้อมูลจะแสดงรายละเอียดข้อมูลของชั้นหิน
ทำชั้นข้อมูลให้ คลิกที่ Table Option > Add Field
ในช่อง Name ทำการตั้งชื่อชั้นข้อมูล ช่อง Type เลือกเป็น Short Integer และช่อง Field Properties เลือกจำนวนอักขระ จากนั้นคลิก OK
เมื่อได้ Field ใหม่แล้วจากนั้นทำการหาข้อมูล โดยเลือกใช้คำสั่ง Select By Attributes...
เริ่มจากการหาข้อมูลของ หินทราย หินทรายแป้ง และหินกรวด โดยการคลิกที่ "DESC_T1" คลิกแสดงข้อมูล Get Unique Values
วิธีการเลือกข้อมูล คือ คลิกเลือก "DESC_T1" = (ข้อมูลชั้นหิน1) OR "DESC_T1" = (ข้อมูลชั้นหิน2) OR ... เมื่อเลือกข้อมูลเรียบร้อยให้คลิกที่ Apply เพื่อแสดงข้อมูล
หลังจากที่ทำการเลือกชั้นข้อมูลเสร็จแล้ว ก็จะทำการใส่ค่าชั้นข้อมูล โดยการคลิกขวาที่ชั้นข้อมูล จากนั้นเลือก Field Calculator...
ใส่ค่า Class จากนั้นคลิก OK ค่าของ Class ก็จะปรากฏอยู่ในชั้นข้อมูล โดยให้เลข 1 แทนค่าของส่วนหินทราย หินทรายแป้ง และหินกรวด
ส่วนต่องมาที่จะทำการใส่ค่าข้อมูลคือ ส่วนที่ 2 หินบะซอลต์ หินปูน และหินตะกอน
โดยการคลิกที่ "DESC_T1" คลิกแสดงข้อมูล Get Unique Values
วิธีการเลือกข้อมูล คือ คลิกเลือก "DESC_T1" = (ข้อมูลชั้นหิน1) OR "DESC_T1" = (ข้อมูลชั้นหิน2) OR ... เมื่อเลือกข้อมูลเรียบร้อยให้คลิกที่ Apply เพื่อแสดงข้อมูล
หลังจากที่ทำการเลือกชั้นข้อมูลเสร็จแล้ว ก็จะทำการใส่ค่าชั้นข้อมูล โดยการคลิกขวาที่ชั้นข้อมูล จากนั้นเลือก Field Calculator...
ใส่ค่า Class จากนั้นคลิก OK ค่าของ Class ก็จะปรากฏอยู่ในชั้นข้อมูล โดยให้เลข 2 แทนหินบะซอลต์ หินปูน และหินตะกอน
ส่วนต่องมาที่จะทำการใส่ค่าข้อมูลคือ ส่วนที่ 3 คือหินโคลน หินดินดาน หินแกรนิต หินไนส์ และควอร์ตไซต์
โดยการคลิกที่ "DESC_T1" คลิกแสดงข้อมูล Get Unique Values
วิธีการเลือกข้อมูล คือ คลิกเลือก "DESC_T1" = (ข้อมูลชั้นหิน1) OR "DESC_T1" = (ข้อมูลชั้นหิน2) OR ... เมื่อเลือกข้อมูลเรียบร้อยให้คลิกที่ Apply เพื่อแสดงข้อมูล
หลังจากที่ทำการเลือกชั้นข้อมูลเสร็จแล้ว ก็จะทำการใส่ค่าชั้นข้อมูล โดยการคลิกขวาที่ชั้นข้อมูล จากนั้นเลือก Field Calculator...
ใส่ค่า Class จากนั้นคลิก OK ค่าของ Class ก็จะปรากฏอยู่ในชั้นข้อมูล โดยให้เลข 2 หินโคลน หินดินดาน หินแกรนิต หินไนส์ และควอร์ตไซต์
ตรวจสอบข้อมูลเพื่อความเรียบร้อย
ทำการแปลงข้อมูล Polygon เป็น Raster โดยการใช้คำสั่ง Polygon to Raster ที่อยู่ใน AceToolbox > Conversion Tools > To Raster > Polygon to Raster
หน้าต่างของ Polygon to Raster ช่องแรกคือช่อง Input Features เลือกข้อมูลที่ต้องการแปลง
ช่องที่ 2 Value field เลือกชั้นข้อมูล class_1 ในตารางที่ทำการใส่ค่าเรียบร้อยแล้ว
ช่องที่ 3 Output Raster Dataset เลือกที่ที่ต้องการจัดเก็บ
และช่องสุดท้าย Cellsize เลือกค่าเป็น 40 จากนั้นคลิก OK
ผลลัพธ์ที่ได้จากการแปลงค่า Polygon เป็น Raster
การจัดข้อมูลเส้น เลือกข้อมูล Polygon = Chon_pro และ ข้อมูล Line = Chon_tran
ทำการคำนวณค่าโดยใช้คำสั่ง Euclidean Distance ที่อยู่ใน AceToolbox > Spatial Analyst Tools > Distance > Euclidean Distance
หน้าต่างของ Euclidean Distance ช่องที่ 1 Input raster or feature source data ให้ใส่ข้อมูลเส้น
ช่องที่ 2 Output distance raster ใส่ที่ที่ต้องการจัดเก็บ
และช่องที่ 3 Output cell size ใส่ค่าข้อมูล จากนั้นไปตั้งค่าที่ Environments...
เลือกตั้งค่าที่ Processing Extent ช่อง Extent เลือก Same as layer Chon_pro จากนั้นเลือกมาด้านล่างเพื่อทำการตั้งค่า Raster Analysis
Raster Analysis ในช่อง Mask เลือก Chon_pro จากนั้นคลิก OK
ตรวจสอบความถูกต้องจากนั้นคลิก OK
ผลลัพธ์ที่ได้จากการทำ Euclidean Distance
ต่อมาคือการใช้คำสั่ง Reclassify โดยเลือกจาก ArcToolbox > Spatial Analyst Tools > Reclass > Reclassify
ช่อง Input เลือกข้อมูลที่ต้องการทำ Reclassify จากนั้นเข้าไปตั้งค่า โดยคลิกที่ Classify
ช่อง Classes ให้เลือก Class เป็น 4 และทำการเปลี่ยนค่าที่ Breale Values ตามลำดับคือ 1000 , 2000 , 3000 และค่าสูงสุดในชั้นข้อมูล จากนั้นคลิก OK
เปลี่ยนค่าในช่อง New Values คือ 0 , 1 , 2 , 3ตามลำดับ และในช่อง Output raster เลือกที่ ที่ต้องการจัดเก็บ จากนั้นคลิก OK
จะได้ผลลัพธ์ของการทำ Chon_tran ที่ทำการ Reclassify แล้ว
ต่อมาจะทำการหาระยะจุดห่าน จุดห่างจากหมู่บ้าน โดยการเปิดชั้นข้อมูลจุดขึ้นมา Chon_vill
ทำการคำนวณค่าโดยใช้คำสั่ง Euclidean Distance ที่อยู่ใน AceToolbox > Spatial Analyst Tools > Distance > Euclidean Distance
หน้าต่างของ Euclidean Distance ช่องที่ 1 Input raster or feature source data ให้ใส่ข้อมูลจุด
ช่องที่ 2 Output distance raster ใส่ที่ที่ต้องการจัดเก็บ
และช่องที่ 3 Output cell size ใส่ค่าข้อมูล จากนั้นไปตั้งค่าที่ Environments...
เลือกตั้งค่าที่ Processing Extent ช่อง Extent เลือก Same as layer Chon_pro จากนั้นเลือกมาด้านล่างเพื่อทำการตั้งค่า Raster Analysis
Raster Analysis ในช่อง Mask เลือก Chon_pro จากนั้นคลิก OK
ตรวจสอบความถูกต้องจากนั้นคลิก OK
ผลลัพธ์ที่ได้จากการทำ Euclidean Distance
ต่อมาคือการใช้คำสั่ง Reclassify โดยเลือกจาก ArcToolbox > Spatial Analyst Tools > Reclass > Reclassify
ช่อง Input เลือกข้อมูลที่ต้องการทำ Reclassify จากนั้นเข้าไปตั้งค่า โดยคลิกที่ Classify
ช่อง Classes ให้เลือก Class เป็น 4 และทำการเปลี่ยนค่าที่ Breale Values ตามลำดับคือ 1000 , 2000 , 3000 และค่าสูงสุดในชั้นข้อมูล จากนั้นคลิก OK
เปลี่ยนค่าในช่อง New Values คือ 0 , 1 , 2 , 3 ตามลำดับ และในช่อง Output raster เลือกที่ ที่ต้องการจัดเก็บ จากนั้นคลิก OK
จะได้ผลลัพธ์ของ eu_vill ที่ทำ Reclassify
จะใช้คำสั่ง Raster Calculator ที่อยู่ในคำสั่ง ArcTollbox > Spatial Analyst > Map Algebra > Raster Calculator
พิมพ์สูตรการคำนวณพื้นที่ จากนั้นเลือกที่จัดเก็บที่ช่อง Output raster และคลิก OK
ผลลัพธ์ที่ได้จากการทำ Raster Calculator หรือการคำนวณพื้นที่
ทำการจัดข้อมูลโดยใช้คำสั่ง Reclassify โดยเลือกจาก ArcToolbox > Spatial Analyst Tools > Reclass > Reclassify
ช่อง Input เลือกข้อมูลที่ต้องการทำ Reclassify จากนั้นเข้าไปตั้งค่า โดยคลิกที่ Classify
ช่อง Classes ให้เลือก Class เป็น 4 และทำการเปลี่ยนค่าที่ Breale Values ตามลำดับคือ 15.09 , 27.47 , 36 จากนั้นคลิก OK
เปลี่ยนค่าในช่อง New Values คือ 1 , 2 , 3 ตามลำดับ และในช่อง Output raster เลือกที่ ที่ต้องการจัดเก็บ จากนั้นคลิก OK
ผลลัพธ์ที่ได้จากการทำ Reclassify สีเขียวเข้มคือพื้นที่ที่เหมาะสมน้อย
สีเขียว(อ่อนปานกลาง) คือพื้นที่ที่เหมาะสมปานกลาง ส่วนสีเขียวอ่อนคือพื้นที่ที่เหมาะสมมากที่สุด
เปิดดูข้อมูลด้วยการคลิกขวาที่ Layer ข้อมูลจากนั้นเลือกที่ Open Attribute Table
ในช่อง COUNT คือช่องที่นับจำนวนพิกเซล ว่ากลุ่มนั้นๆ มีกี่พิกเซล หน่วยของพื้นมีหน่วยเป็นไร้
Video How to 9 : การวิเคราะห์เเบบจำลองเเรสเตอร์ (Raster model)
