Lab9



การวิเคราะพื้นที่เหมาะสมสำหรับการสร้างที่ฝังกลบขยะ

กำหนดประเด็นปัญหา





1.    ความเหมาะสมของพื้นที่ฝังกลบขยะในแต่ละระดับมีพื้นที่เท่าใด
2.   ในแต่ละพื้นที่ (อำเภอ) มีความเหมาะสมอยู่ในระดับใดบ้าง
กำหนดปัจจัยหรือตัวแปรวิเคราะห์
1. ความลาดชัน
2. ลักษณะชั้นหินฐาน
3. ระยะห่างจากถนนสายหลัก
4. ระยะห่างจากชุมชน
การเตรียมข้อมูล
1. การนำเข้าข้อมูลจุดความสูง การสร้างแบบจำลองความสูงเชิงเลข (DEM) การคำนวณความลาดชัน
2. การจัดการและแก้ไขข้อมูลคุณลักษณะ





 เปิดข้อมูลในโฟลเดอร์ Landfill ที่ชื่อว่า Chon_elv   และ Chon_pro ให้เลากข้อมูลเข้ามาเปิดในพื้นที่แสดงข้อมูลได้เลย



ต่อไปเราจะประมานค่าช่วงโดยใช้ Interpolation โดยใช้ค่า Z จะกำหนดกริด เป็น 40 เมตร โดยใช้วิธี IDW  ให้เราไปหาเครื่องมือใน ArcTool box ที่ชื่อว่า  Spatial Analyst Tools >> Inter polation >> IDW
   ดับเบิ้ลคลิกเปิดขึ้นมา จะมีหน้าต่างขึ้นมาให้ ในช่อง Input point features  เลือกเป็น Chon_elv ส่วนในช่อง Z value Field ให้เปลี่ยนเป็น ele  ในช่อง Output raster เลือกที่เก็บงาน ช่อง Output cell size (optional)กำหนดค่าเป็น 40  และให้ไปเซตค่าต่อที่ Environments  ให้เราคลิกเข้าไป


จะขึ้นหน้าต่างมาให้ ให้เราไปที่หัวข้อ  Processing Extent ในช่อง Extent เลือกเป็น Same as layer Chon_pro


จากนั้นไปอีกหัวข้อ คือ Raster Analysis  ในช่อง Mask เปลี่ยนเป็น Chon_pro จากนั้นกด OK



ข้อมูลที่ทำการ IDW ก็จะออกมาเป็นดังนี้


ต่อไปเราจะไปหาชั้นความสูง ให้เราไปหาเครื่องมือ ArcToolbox ที่ชื่อว่า  Spatial Analyst Tools>>Surface  ดับเบิ้ลคลิกขึ้นมา จะขึ้นหน้าต่างมาให้ ในช่อง Input raster  ใส่ข้อมูล idw  ช่อง Output raster  เลือกที่เก็บงาน ส่วนช่อง  Output measurement เลือกเป็น DEGREE และช่อง Z Factor (optional) ใส่เป็น 1 จากนั้นกด OK


ภาพที่ได้จากการทำ Slope ออกมาดังนี้

ให้เราทำการ Reclassify ต่อ ให้ไปเอาเครื่องมือใน ArcToolbox  Spatial Analyst Tools>> Reclass>>Reclassify ในช่อง Input raster  เลือกข้อมูล slope ช่อง Reclass feid  เลือกเป็น Value     ช่อง Reclassification จะมีค่ามาให้  Output raster เลือกที่เก็บงาน และให้เราไปกดที่  Classify

 จะขึ้นหน้าต่าง Classification มาให้ ให้เราเซตค่าในช่อง Method ให้เราเปลี่ยนเป็น Equal Interval จากนั้นในช่อง Classes เปลี่ยนเป็น 4 จากนั้นให้ไปเปลี่ยนเป็น Manual เหมือนเดิม ต่อไปให้เราไปเซตค่าในช่อง Break Values เป็น 5 , 15 , 30 และ 39.249115 จากนั้นกด OK  และ OK




จะได้ผลจากการ Reclassify ออกมาดังนี้



ต่อไปเราจะเตรียมชั้นหิน ให้เราไปเปิดข้อมูลเพิ่มมาอีก 1 ตัว คือ  Chon_geo  จากนั้นให้เปิดดูตารางของตัวข้อมูล โดยการคลิกขวา เลือก Open Attribute Table



ให้เราสร้างฟิลด์ขึ้นมาใหม่ ชื่อว่า Score  ให้ไปคลิกที่   Add Field…


จะขึ้นหน้าต่าง Add Field มาให้ ให้เราตั้งชื่อ จากนั้นกด OK ก็จะได้ฟิลด์ใหม่มา


เมื่อหน้าต่างตารางขึ้นมา ให้เราเพิ่มฟิลด์ให้ตัวตาราง โดยไป เลือก  
 Select By Attributes…


จะขึ้นหน้าต่าง Select by Attribute มาให้ และในช่อง SELECT'FROM Chon_geo WHERE  ให้เราใส่ข้อมูลจาก Method คือ DESC_T1 และ = / OR 

สรุปคือ  "DESC_T1" = 'หินกรวดมน หินฟิลไลต์ หินทรายเนื้อภูเขาไฟกึ่งแปรสภาพ และหินชนวน' OR "DESC_T1" = 'หินทราย หินทรายเนื้อภูเขาไฟ หินปูนเนื้อดิน และหินปูนเนื้อไข่ปลา หินไรโอลิติกทัฟฟ์และแอนดิซิติกทัฟฟ์กึ่งแปลสภาพ หินดินดานเลนซ์ หินปูนหินเชิร์ต' OR "DESC_T1" = 'หินทราย หินปูนเนื้อดิน หินดินดาน และหินเชิร์ต' OR "DESC_T1" = 'หินทรายอาร์ดคส สีขาว มีหินกรวดมนและหินดินดานแทรกสลับ' OR "DESC_T1" = 'หินฟิลไลต์ หินฟิลไลต์เนื้อปนถ่าน และหินฟิลไลต์เนื้อซิลิกา'


ใส่ข้อมูลดังภาพ จากนั้นกด Apply


จะแสดงข้อมูลในตารางเป็นไฮไลต์มาให้



จากนั้นไปคลิกขวาที่ฟิลด์ Score ที่เราสร้าง   Field Calculator…


จะมีหน้าต่าง Field Calculator ขึ้นมาให้ จากนั้นให้เราใส่เลข 1 ลไปในช่อง Score จากนั้นกด OK


ตารางก็จะเปลี่ยนเป็นเลข 1 ตามที่เราคำนวณค่าเอาไว้



ต่อไปให้เราไปเลือก    Select By Attribute จะขึ้นหน้าต่างมาให้ ตรงช่อง SELECT'FROM Chon_geo WERE ให้เราใส่ข้อมูลจาก Method คือ DESC_T1 และ = / OR 

สรุปคือ "DESC_T1" = 'ตะกอนชายฝั่งทะเลโดยอิทธิพล น้ำขึ้น-ลง โคลน ตม ทรายละเอียด ป่าชายเลน พรุ พื้นที่ชื้นแฉะ ชวากทะเล' OR "DESC_T1" = 'ตะกอนชายฝั่งทะเลโดยอิทธิพลคลื่น ทรายหาด สันดอนทราย เนินทราย' OR "DESC_T1" = 'ตะกอนที่ราบสะสมตัวโดยทางน้ำ ตะกอนกรวด ทราย ทรายแป้ง ดิน สะสมตามร่องน้ำ และที่ราบน้ำท่วมถึง' OR "DESC_T1" = 'ตะกอนสะสมตัวตามเชิงเขา ตะกอนผุพังอยู่กับที่ เศษหิน ทราย ดิน ชั้นศิลาแลง' OR "DESC_T1" = 'ตะกอนเนินตะพัก ชั้นกรวด ทราย ทรายแป้ง ดิน และหินลูกรัง'


เมื่อใส่เสร็จให้กด Apply ก็จะมีข้อมูลแสดงขึ้นมาให้


ต่อไปให้เราใส่ค่าคะแนน โดยการไปคลิกขวาที่ Score  Field Calculator… 


จะมีหน้าต่าง Field Calculator มาให้ ในช่อง Score ใส่เลข 2 จากนั้นกด OK



ต่อไปให้เราไปเลือก    Select By Attribute จะขึ้นหน้าต่างมาให้ ตรงช่อง SELECT'FROM Chon_geo WERE  ให้เราใส่ข้อมูลจาก Method คือ DESC_T1 และ = / OR 

สรุปคือ "DESC_T1" = 'หินฟิลไลต์ หินฟิลไลต์เนื้อปนถ่าน และหินฟิลไลต์เนื้อซิลิกา' OR "DESC_T1" = 'หินไบโอไทต์แกรนิต ทัวมารีนแกรนิต ไบโอไทต์มัสโคไวต์แกรนิต แกรโนไดโอไรต์ มัสโคไวต์ทัวมารีนแกรนิต ไบโอไทต์ทัวมารีนแกรนิต' OR "DESC_T1" = ' '




เมื่อใส่เสร็จให้กด Apply ก็จะมีข้อมูลแสดงขึ้นมาให้


ต่อไปให้เราใส่ค่าคะแนน โดยการไปคลิกขวาที่ Score  Field Calculator… และใส่เลข 3 ในช่อง Score  จากนั้นกด OK


ก็จะมีค่า Score ที่เราใส่ไปแสดงขึ้นมาในตาราง



จากนั้นให้เราไปหาเครื่องมือใน ArcToolbox ที่ชื่อว่า Conversion Tools>> To Raster >>polygon to Raster จะมีหน้าต่างขึ้นมาให้ ในช่อง InputFeatures  ใส่ข้อมูล Chon_geo  ช่อง  Value Feild เลือกใส่ Score ในช่อง Outputraster Dataset เลือกที่เก็บงาน ในช่อง Cell assignment (optional) ใส่ค่าเดิม ส่วนในช่อง Priority feild ใส่ค่าเดิม และช่อง Cellsize (optional) เปลี่ยนเลขเป็น 40 จากนั้นกด OK



จะได้ข้อมูลออกมาดังนี้


ต่อไปเราจะทำข้อมูลเส้นถนน ไปหาเครื่องมือใน ArcToolbox ที่ชื่อว่า  -- Spatial Analyst Tools --Distance >> Euclidance Distance จะขึ้นหน้าต่างมาให้ จากนั้นไปที่  Environments


จะขึ้นหน้าต่างมาให้เราเซตค่า ไปที่หัวข้อProcessing Extent  เลือกเป็น Same as layer Chon_pro



และอีกหัวข้อคือ Raster Analysis  ในช่อง Mask เลือกเป็น Chon_pro จากนั้นกด OK และ  OK



จะได้ข้อมูลแสดงออกมาดังนี้



ให้เราทำการ Reclassify ต่อ ให้ไปเอาเครื่องมือใน ArcToolbox   Spatial Analyst Tools >> Reclass>> Reclassify  ในช่อง Input raster  เลือกข้อมูล tran1 ช่อง Reclass field  เลือกเป็น Value     ช่อง Reclassification จะมีค่ามาให้ Output raster  เลือกที่เก็บงาน และให้เราไปกดที่   Classify

จะขึ้นหน้าต่าง Classification มาให้ ให้เราเซตค่าในช่อง Method ให้เราเปลี่ยนเป็น Equal Interval จากนั้นในช่อง Classes เปลี่ยนเป็น 4 จากนั้นให้ไปเปลี่ยนเป็น Manual เหมือนเดิม ต่อไปให้เราไปเซตค่าในช่อง Break Values เป็น 1000 , 2000 , 3000 และ 21549.904297 จากนั้นกด OK  และ OK



ภาพที่แสดงออกมาจะเป็นแบบนี้




ต่อไปเราจะทำข้อมูลเส้นถนน ไปหาเครื่องมือใน ArcToolbox ที่ชื่อว่า  Processing Extent จะขึ้นหน้าต่างมาให้ โดยเซตค่า ในช่อง Input ใส่ข้อมูล Chon_vill ช่อง Output ใส่ที่เก็บงาน และช่อง Output cell size ใส่เลข 40  จากนั้นไปที่    



จะขึ้นหน้าต่างมาให้เราเซตค่า ไปที่หัวข้อ  เลือกเป็น Same as layer Chon_pro



และอีกหัวข้อคือ  Raster Analysis ในช่อง Mask เลือกเป็น Chon_pro จากนั้นกด OK และ  OK



จะได้ข้อมูลออกมาดังนี้


ไม่มีความคิดเห็น:

แสดงความคิดเห็น