Nghiên cứu xác định các đối tượng rừng dựa vào ảnh vệ tinh landsat 8 và sentinel 2 tại tỉnh lào cai

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƢỜNG ĐẠI HỌC LÂM NGHIỆP - BÙI THANH VIỆT NGHIÊN CỨU XÁC ĐỊNH ĐỐI TƢỢNG RỪNG DỰA VÀO ẢNH VỆ TINH LANDSAT VÀ SENTINEL TẠI TỈNH LÀO CAI CHUYÊN NGÀNH: LÂM HỌC MÃ NGÀNH: 8620201 LUẬN VĂN THẠC SĨ LÂM HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: TS BÙI MẠNH HƢNG PGS.TS PHÙNG VĂN KHOA Hà Nội, 2019 i L Ờ I C A M Đ O A N Tôi cam đoan, cơng trình nghiên cứu riêng tơi Các số liệu, kết nêu luận văn trung thực chưa công bố công trình nghiên cứu khác Nếu nội dung nghiên cứu tơi trùng lặp với cơng trình nghiên cứu cơng bố, tơi xin hồn tồn chịu trách nhiệm tuân thủ kết luận đánh giá luận văn Hội đồng khoa học Hà Nội, ngày tháng năm 2019 Ngƣời cam đoan Bùi Thanh Việt ii L Ờ I C Ả M Ơ N Sau năm học tập nghiên cứu Trường Đại học Lâm Nghiệp Bằng kiến thức thân giúp đỡ bảo tận tình thầy (cơ) giáo Đến tơi hồn thành luận văn thạc sỹ, xin chân thành cảm ơn giúp đỡ q báu Tơi xin chân thành cảm ơn PGS.TS Phùng Văn Khoa - Trưởng khoa Quản lý tài nguyên rừng môi trường, TS Bùi Mạnh Hưng hướng dẫn tận tình tạo điều kiện, giúp đỡ suốt thời gian thực hoàn thành luận văn Tôi xin chân thành cảm ơn giúp đỡ quý thầy cô giáo Trường Đại học Lâm nghiệp Xin gửi tới gia đình, bạn bè, đồng nghiệp lịng biết ơn chân thành tạo điều kiện trình học tập hồn thành luận văn Mặc dù thân cố gắng trình thực nghiên cứu thời gian trình độ có hạn nên chắn luận văn nhiều thiếu sót hạn chế Kính mong nhận góp ý thầy cô bạn./ Hà Nội, ngày tháng năm 2019 Tác giả Bùi Thanh Việt iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN II MỤC LỤC iii DANH MỤC HÌNH ẢNH v DANH MỤC BẢNG BIỂU vi DANH MỤC CÁC TỪ VIẾT TẮT viii ĐẶT VẤN ĐỀ Chƣơng TỔNG QUAN NGHIÊN CỨU 1.1 Các khái niệm rừng 1.2 Phân loại rừng 1.3 Lịch sử phát triển viễn thám 1.4 Khái niệm viễn thám 1.4.1 Định nghĩa 1.4.2 Nguyên lý viễn thám 1.5 Đặc điểm thông số kỹ thuật ảnh vệ tinh Landsat8 Sentinel2 11 1.5.1 Ảnh Landsat8 11 1.5.2 Ảnh Sentinel2 11 1.6 Một số số dùng phân loại ảnh viễn thám 13 1.6.1 Chỉ số sai khác thực vật NDVI 13 1.6.2 Chỉ số tăng cường lớp thực vật EVI 13 1.6.3 Chỉ số thực vật có hiệu chỉnh ảnh hưởng đất SAVI 13 1.6.4 Chỉ số khác biệt nước NDWI 13 1.6.5 Chỉ số NBR (Normalized Burn Ratio) 13 1.7 Tình hình nghiên cứu ứng dụng viễn thám GIS lâm nghiệp 14 1.7.1 Trên giới 14 1.7.2 Tại Việt Nam 18 Chƣơng MỤC TIÊU, PHẠM VI, NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 24 2.1 Mục tiêu nghiên cứu 24 2.1.1 Mục tiêu tổng quát 24 2.1.2 Mục tiêu cụ thể 24 2.2 Đối tượng phạm vi nghiên cứu nghiên cứu 24 2.2.1 Đối tượng nghiên cứu 24 iv 2.2.2 Phạm vi nghiên cứu: 24 2.3 Nội dung nghiên cứu 25 2.4 Phương pháp nghiên cứu 25 2.4.1 Phương pháp kế thừa số liệu 25 2.4.2 Phương pháp điều tra thực địa 25 2.4.3 Phương pháp xử lý ảnh 28 2.4.4 Giải đoán ảnh 31 2.4.5 Đánh giá độ xác đồ 34 Chƣơng KHÁI QUÁT ĐIỀU KIỆN TỰ NHIÊN - KINH TẾ - XÃ HỘI TỈNH LÀO CAI 36 3.1 Điều kiện tự nhiên 36 3.1.1 Vị trí địa lý 36 3.1.2 Địa hình, địa mạo 37 3.1.3 Khí hậu, thủy văn 38 3.1.4 Đất đai, thổ nhưỡng 42 3.2 Đặc điểm Kinh tế - Xã hội 44 3.2.1 Nguồn nhân lực 44 3.2.2 Thực trạng Kinh tế -Xã hội 45 3.3 Thực trạng rừng đất lâm nghiệp 49 Chƣơng KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 50 4.1 Xây dựng mẫu khóa ảnh giải đốn trạng thái rừng khu vực nghiên cứu 50 4.2 Kết giải đốn ảnh theo phương pháp phân loại có kiểm định phương pháp dùng số 51 4.2.1 Kết giải đoán ảnh theo phương pháp phân loại có kiểm định 51 4.2.2 Kết phương pháp phân loại ảnh dùng số 53 4.3 Đánh giá độ xác phương pháp phân loại ảnh 72 KẾT LUẬN, TỒN TẠI VÀ KIẾN NGHỊ 77 TÀI LIỆU THAM KHẢO 79 PHỤ LỤC v D A N H M Ụ C HÌ N H Ả N H Hình 1.1: Ngun lý thu nhận ảnh viễn thám 10 Hình 2.1: Hình dạng khoảng cách chùm ô mẫu 26 Hình 2.2: Vị trí tên mẫu chùm 27 Hình 2.3 Hình dạng đo đếm rừng gỗ tự nhiên đất chưa có rừng 27 Hình 2.4 Hình dạng, vị trí đo đếm rừng gỗ trồng 28 Hình 2.5: Hình dạng, vị trí đo đếm rừng tre nứa 28 Hình 2.6: Lọc ảnh vệ tinh theo khu vực nghiên cứu 29 Hình 2.7: Tiền xử lý ảnh vệ tinh ( Lọc mây, bóng núi) 29 Hình 2.8: Tính tốn số ảnh vệ tinh 29 Hình 2.9: Xuất ảnh vệ tinh sau tiền xử lý 30 Hình 2.10: Ảnh vệ tinh Landsat khu vực nghiên cứu tiền xử lý 30 Hình 2.11: Ảnh vệ tinh Sentinel khu vực nghiên cứu tiền xử lý 31 Hình 2.12: Cây phân loại tự động phần mềm eCoginion 32 Hình 2.13: Dùng giá trị xám độ để tách trạng thái rừng 34 Hình 3.1: Bản đồ hành tỉnh Lào Cai 37 Hình 4.1: Sự phân bố mẫu khóa ảnh giải đốn 50 Hình 4.2: Kết phân loại tự động có kiểm định ảnh vệ tinh Landsat 852 Hình 4.3: Kết phân loại tự động có kiểm định ảnh vệ tinh Sentinel2 53 Hình 4.4: Kết phân loại số ảnh vệ tinh Landsat 70 Hình 4.5: Kết phân loại số ảnh ảnh vệ tinh Sentinel 71 Hình 4.6: Ma trận sai số phân loại tự động có kiểm định ảnh vệ tinh Sentinel2 73 Hình 4.7: Ma trận sai số phân loại tự động có kiểm định ảnh vệ tinh Landsat 73 Hình 4.8: Ma trận sai số phân loại theo ngưỡng số ảnh vệ tinh Sentinel2.74 Hình 4.9: Ma trận sai số phân loại theo ngưỡng số ảnh vệ tinh Landsat8 74 vi DANH MỤC BẢNG BIỂU Bảng 1.1: Tóm tắt phát triển viễn thám qua kiện Bảng 1.2: Đặc điểm ảnh vệ tinh Landsat Sentinel 12 Bảng 2.1: Mẫu khóa ảnh giải đốn ảnh… 24 Bảng 2.2: Ma trận sai số phân loại 35 Bảng 3.1: Diện tích địa hình độ cao tương đối so với mực nước biển tỉnh Lào Cai chia theo huyện………… 38 Bảng 3.2: Tổng hợp nhóm, loại đất tỉnh Lào Cai 42 Bảng 4.1: Tổng hợp số lượng mẫu khóa ảnh theo trạng thái rừng 51 Bảng 4.2:Kết phân loại tự động ảnh vệ tinh Landsat Sentinel2 51 Bảng 4.3: Thống kê giá trị số cho đối tượng nghiên cứu ảnh vệ tinh Landsat 54 Bảng 4.4: Thống kê giá trị số cho đối tượng nghiên cứu ảnh vệ tinh Sentinel 54 Bảng 4.5: Kết so sánh độc lập số NDVI trạng thái đất có rừng đất khơng có rừng ảnh vệ tinh Sentinel 55 Bảng 4.6: Kết so sánh độc lập số NDVI của trạng thái đất nông nghiệp đất trống ảnh vệ tinh Sentinel 56 Bảng 4.7: Kết so sánh độc lập số NDVI của trạng thái rừng trồng rừng tự nhiên ảnh vệ tinh Sentinel 57 Bảng 4.8: Kết so sánh độc lập số NDVI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Sentinel 58 Bảng 4.9: Kết so sánh độc lập số SAVI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Sentinel 59 Bảng 4.10: Kết so sánh độc lập số NDWI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Sentinel 60 vii Bảng 4.11: Kết so sánh độc lập số NBR trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Sentinel 61 Bảng 4.12: Tóm tắt kết so sánh khác biệt số đối tượng rừng ảnh vệ tinh Sentinel 61 Bảng 4.13: Kết so sánh độc lập số NDVI trạng thái đất có rừng đất khơng có rừng ảnh vệ tinh Landsat 62 Bảng 4.14: Kết so sánh độc lập số NDVI trạng thái đất nông nghiệp đất trống ảnh vệ tinh Landsat 63 Bảng 4.15: Kết so sánh độc lập số NDVI trạng thái rừng trồng rừng tự nhiên ảnh vệ tinh Landsat 64 Bảng 4.16: Kết so sánh độc lập số NDVI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Landsat 65 Bảng 4.17: Kết so sánh độc lập số EVI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Landsat 66 Bảng 4.18: Kết so sánh độc lập số SAVI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Landsat 67 Bảng 4.19: Kết so sánh độc lập số NDWI trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Landsat 68 Bảng 4.20: Kết so sánh độc lập số NBR trạng thái rừng hỗn giao rừng tự nhiên thường xanh ảnh vệ tinh Landsat 69 Bảng 4.21: Tóm tắt kết so sánh khác biệt số đối tượng rừng ảnh vệ tinh Landsat 69 viii D A N H STT Viết tắt M Ụ C CÁ C T Ừ V I Ế T T Ắ T Viết đầy đủ BNNPTNT Bộ Nông nghiệp phát triển nông thôn DT Đất trống EVI Enhancement vegetation index (Chỉ số tăng cường lớp thực vật) FAO (Food and Agriculture Organization of the United Nations): Tổ chức Lương thực Nông nghiệp Liên Hiệp Quốc HG Rừng hỗn giao MKA Mẫu khóa ảnh NBR Normalized Burn Ratio NDVI normalized difference vegetation index (Chỉ số sai khác thực vật) NDWI Normalized difference vegetation index (Chỉ số khác biệt nước ) 10 NIR Kênh cận hồng ngoại ảnh vệ tinh 11 NN Đất nông nghiệp 12 ÔTC Ô Tiêu chuẩn 13 RED Kênh đỏ ảnh 14 RTG Rừng trồng gỗ 15 RTCG 16 RTN 17 RS 18 SAVI 19 SWIR1 Kênh hồng ngoại song ngắn 20 SWIR2 Kênh hồng ngoại song ngắn Rừng Tự nhiên gỗ Rừng tự nhiên (Remote Sensing): Viễn thám Soild Ajusted Vegetation Index (Chỉ số thực vật có hiệu chỉnh ảnh hưởng đất) ĐẶT VẤN ĐỀ Rừng phận tổ thành quan trọng sinh hành tinh Tài nguyên rừng không mang lại giá trị kinh tế to lớn mà cịn có ý nghĩa xã hội ngày tăng, rừng có chức cân sinh thái, bảo vệ mơi trường sống, chức phịng hộ, văn hoá cảnh quan Ngày nhiệm vụ bảo vệ, phát triển kinh doanh rừng không vấn đề riêng ngành Lâm nghiệp, chuyên gia rừng mà trở thành vấn đề kinh tế có quy mơ tồn xã hội Trong thời đại cách mạng khoa học kỹ thuật ngày nay, vai trò rừng trở nên quý giá Sự phát triển văn minh nhân loại dẫn tới thay đổi môi trường sống làm tăng thêm tính phức tạp mối quan hệ phụ thuộc người vào môi trường Ở nước phát triển, q trình cơng nghiệp hố, hố học hố thị hố làm nhiễm bẩn nặng nề môi trường sống Ở nước phát triển, vùng nhiệt đới, với phương thức trồng trọt du canh khai thác rừng không hợp lý, người huỷ diệt diện tích rừng rộng lớn, phá huỷ thành phần sinh học cần thiết cho tồn thân Ngày nay, rừng trở thành bối cảnh độc đáo bên sản xuất công nghiệp bên môi trường sống người Cuộc khủng hoảng lớn môi sinh phạm vi tồn cầu chứng tỏ vai trị khơng thể thay rừng việc trì, cân sinh thái, bảo vệ môi trường sống Những chức phòng hộ rừng việc nâng cao suất mùa màng, chống xói mịn, điều hồ nguồn nước, chống cát bay, Để đáp ứng nhu cầu đời sống văn hoá, tinh thần ngày cao, ngày rừng trở thành nơi nghỉ mát, an dưỡng, du lịch…với 16,2 triệu rừng đất rừng quy hoạch cho mục đích lâm nghiệp, chiếm gần 50% diện tích tự nhiên (Theo số liệu kiểm kê rừng 2017) Nước ta nằm khu vực nhiệt đới ẩm, lượng mưa 81 núi Việt Nam, Nhà xuất trị quốc gia, Hà Nội 22 Mai Đình Yên (1998), “Các vấn đề vê môi trường Việt Nam hành động cần ưu tiên bảo vệ môi trường”, Tạp chí khoa học mơi trường, (số 12), tr 15-16 23 Nguyễn Thị Yến (2001), Giáo trình phương pháp viễn thám, Đại học Mỏ địa chất, Hà Nội 24 Ủy ban nhân dân huyện Thường Xuân (2004), Báo cáo tóm tắt quy hoạch phát triển Lâm nghiệp huyện Thường Xuân thời kỳ 2001-2010 25 Ủy ban nhân dân huyện Thường Xuân (2004), Báo cáo tóm tắt quy hoạch sử dụng đất huyện Thường Xuân thời kỳ 2003-2010 26 Ủy ban nhân dân huyện Thường Xuân (2004), Niên giám thống kê 2003 27 Ủy ban nhân dân tỉnh Thanh Hoá (2004), Báo cáo tóm tắt quy hoạch phát triển Lâm nghiệp tỉnh Thanh Hoá thời kỳ 2001-2010 28 Nguyễn văn chiến (2018) “Nghiên cứu điều tra thống kê tài nguyên rừng ảnh vệ tinh VNRedSat-1 xã Xuân Sơn huyện Tân Sơn tỉnh Phú Thọ” 29 Hoàng Hưng (2012) “Nghiên cứu lựa chọn phương pháp kiểm kê rừng phục vụ cho việc chi trả dịch vụ môi trường rừng xã Miền Đồi – Huyện Lạc Sơn – Tỉnh Hịa Bình” 30.Đinh Văn Tùng (2018) “ Nghiên cứu thay đổi diện tích rừng phịng hộ làm sở đề xuất giải pháp quản lý lưu vực bền vững địa bàn tỉnh Vĩnh Phúc” 31 Báo cáo kết rà soát , điều chỉnh quy hoạch ba loại rừng tỉnh Lào Cai, giai đoạn 2016-2025 Tiếng Anh 32 Burrough (1986), Principle of geographical information system for land resources assessment, Clarendon Press-Oxford 82 33 Charlie Navanugraha (1996), Land use/Land cover change, A case study in ThaiLand, Falculty of environmental and resource, Mahidol University, Thai land 34 Ding Yuan et al (1998), “Survey of multispectral methods for land cover change analysis” Remote sensing change detection: environmental monitoring methods and applications, Ann arbor press 35 Gautam A.p et al (2002), GIS assessment of land use/land cover changes associated with community forestry implimentation in the middle hills of Nepal, Mountain research and development 36 Luu Van Nang et al (1999) “ Potential of 1RS -1 Panchromatic Sattelite image Data for village level Land use Planning: An Example from the Sector project in VietNam”, Application of Resource Information Technologies (GIS/GPS/RS) in Forest land &Resource Management, Trl24-134, Hanoi 37 Muh dimyati, Kei Mizuno Shintaro Kobayashi and Teitaro Kitamura (1996), An analyst of land use!land cover change using the combination of MSS Landsat and land use map_A case study in YogYakata Indonesia, Kyoto University 38 Robin S.Reid, Russell Lkruska, Nyawira Muthui (2002), land use and land cover dynamics in response to changes in climate, biological and socio-political forces, the case of Southwestern Ethiopia 39 Shunji Murai (1991), Applications of Remote sensing in Asia and Oceania, Asian Association on Remote sensing 40 W.G.Rees (1993), Physical principle of remote sensing, Cambridge university press 41 Yushang Zhou and David L.Skoke (2001), Cultivated land use change analysis and modeling, A case study in the Earst region of China, 83 Center for Global Change and Earth Observations, Michigan State University 42 Yushang Zhou and David L.skole (2001), Cultivated land use change analysis and modeling, Center for Global Change and Earth Observations, Michigan State University Website 43 Bài học Hệ số Cohen’s kappa (Cohen's kappa coefficient) xem ngày 18/4/2019 Có link sau: https://vimentor.com/vi/lesson/he-so-cohen-s-kappa-cohen-s-kappacoefficient PHỤ LỤC Phụ lục Phiếu đo đếm Phiếu 1: Mốc tâm ô mẫu, điểm tham chiếu FIPI-NFI, CK5 Phiếu 1: Mốc điểm tham chiếu Ngày Tháng Năm Tên chùm ô Tên ô mẫu Tâm ô mẫu TĐX Loại mốc TĐY Khoảng cách đến tâm ô mẫu 0,01 m Hướng đến tâm ô mẫu độ Mô tả (đối tượng gần mốc) Điểm tham chiếu TĐX Loại đối tượng TĐY Khoảng cách đến tâm ô mẫu 0,01 m Hướng đến tâm ô mẫu độ Mô tả Điểm tham chiếu TĐX Loại đối tượng TĐY Khoảng cách đến tâm ô mẫu 0,01 m Hướng đến tâm ô mẫu độ Mô tả Điểm tham chiếu Loại đối tượng Mô tả Ghi TĐX TĐY Khoảng cách đến tâm ô mẫu 0,01 m Hướng đến tâm ô mẫu độ Phiếu 2: Mô tả ô mẫu Trang / FIPI-NFI, CK5 Mô tả ô mẫu Tên ô mẫu Tên chùm ô Ngày Năm Tháng Mã KN tiếp cận PP đánh giá Tên tỉnh Tên huyện Tên xã Tên thôn, Độ dốc Hướng phơi độ Số ID Tên file máy tính (Ctenchumo_Pp_nhom.jpg ) Ảnh nhóm điều tra (Ơ mẫu số chùm ơ) hướng khoảng cách m Ghi Số ID Tên file máy tính (Ctenchumo_Pp_tam.jpg) Ảnh tâm mẫu (Từ điểm đứng chụp đến tâm ô mẫu) hướng khoảng cách m Ghi Tỷ lệ che phủ tán rừng (%) Bắc Tâm Đông Tỷ lệ che phủ bụi (%) Nam Tây Tỷ lệ che phủ thảm tươi (%) Chiều cao TB bụi (0,1m) Chiều cao TB thảm tươi (0,01m) Trạng thái Tác nhân gây hại Trạng thái trước (1) Thời gian thay đổi (2) Chức rừng Mức độ thiệt hại Tác động người Cây tán rừng Số năm TĐ(1-5, 5+) (1) (2) Chủ quản lý (3) Đề xuất quản lý Đất trồng rừng chưa thành rừng rừng trồng chưa có trữ lượng (D1,35,0m H C Phiếu 4: Điều tra gỗ Trang / Tên chùm ô Tên ô mẫu gỗ SH Tên loài D1,3 (0,1 cm) Chu vi H-D Chất Hvn(0,1 m) Hdc(0,1m) (0,1 cm) (0,01 m) lượng Ghi Phụ lục 02 Phiếu 01 MKA: PHIẾU M TẢ MẪU KHĨA ẢNH VỆ TINH Mẫu khố ảnh số: Ngày điều tra: Vị trí: Người ĐT: Hướng phơi: Toạ độ: Tỉnh: Toạ độ X: Huyện: Toạ độ Y: Xã: Độ cao: Tiểu khu Hệ toạ độ: Khoảnh Mô tả thực địa Trạng thái Tiết diện ngang G1 G2 G3 G4 G5 GTB điểm bitterlich Mô tả ảnh vệ tinh Số hiệu cảnh ảnh: Chiều cao H1 H2 H3 H4 H5 HTB Thời gian thu nhận ảnh: TB điểm bitterlich Trữ lượng bình quân Độ tàn che TB: Loài ưu Ảnh chụp thực địa vị trí tâm mẫu Tên ảnh thực địa Ảnh vệ tinh Phụ lục 03: Danh sách tọa độ MKA giải đoán Tọa độ Y Tỉnh 432381 2493734 Lào Cai Bắc Hà 2 432381 2493584 Lào Cai Bắc Hà 3 432381 2493434 Lào Cai Bắc Hà 4 432531 2493434 Lào Cai Bắc Hà 5 432681 2493434 Lào Cai Bắc Hà 424381 2485734 Bắc Hà 424381 2485584 424381 2485434 424531 2485434 10 424681 2485434 400381 2469734 400381 2469584 Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai STT 11 12 Số hiệu chùm ô Số hiệu ODD 187 231 317 Tọa độ X Huyện Bắc Hà Bắc Hà Bắc Hà Bắc Hà Bảo Thắng Bảo Thắng Xã Thải Giàng Phố Thải Giàng Phố Thải Giàng Phố Thải Giàng Phố Thải Giàng Phố Nậm Mòn Nậm Mòn Nậm Mòn Nậm Mòn Nậm Mòn Trạng thái Ngoại nghiệp RTNCG DT RTNCG NN RTNCG RTG RTNCG RTG RTNCG RTG Gia Phú NN Gia Phú NN 13 400381 2469434 14 400531 2469434 15 400681 2469434 416381 2469734 17 416381 2469584 18 416381 2469434 19 416531 2469434 20 416681 2469434 432381 2469734 22 432381 2469584 23 432381 2469434 24 432531 2469434 25 432681 2469434 448381 2461734 27 448381 2461584 28 448381 2461434 29 448531 2461434 30 448681 2461434 16 21 26 318 319 368 Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Bảo Thắng Bảo Thắng Bảo Thắng Bảo Thắng Bảo Thắng Bảo Thắng Bảo Thắng Bảo Thắng Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Gia Phú DT Gia Phú DT Gia Phú DT TT Phố Lu TT Phố Lu TT Phố Lu TT Phố Lu TT Phố Lu Điện Quan Điện Quan Điện Quan Điện Quan Điện Quan Xuân Thượng Xuân Thượng Xuân Thượng Xuân Thượng Xuân Thượng RTG RTG RTG RTG RTG NN NN NN NN DT NN NN NN RTG RTG 31 456381 2445734 32 456381 2445584 33 456381 2445434 34 456531 2445434 35 456681 2445434 360381 2509734 37 360381 2509584 38 360381 2509434 39 360531 2509434 40 360681 2509434 408381 2501734 42 408381 2501584 43 408381 2501434 44 408531 2501434 45 408681 2501434 384381 2469734 47 384381 2469584 48 384381 2469434 36 41 46 451 111 147 315 Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bảo Yên Bát Xát Bát Xát Bát Xát Bát Xát Bát Xát Mường Khương Mường Khương Mường Khương Mường Khương Mường Khương Sa Pa Sa Pa Sa Pa Long Khánh Long Khánh Long Khánh Long Khánh Long Khánh A Mú Sung A Mú Sung A Mú Sung A Mú Sung A Mú Sung Lùng Vai Lùng Vai Lùng Vai Lùng Vai Lùng Vai Hầu Thào Lao Chải Lao Chải HG RTNCG DT RTNCG HG RTNCG RTNCG RTNCG RTNCG RTNCG NN RTNCG DT DT RTNCG RTNCG DT RTNCG 49 384531 2469434 50 384681 2469434 400381 2453734 52 400381 2453584 53 400381 2453434 54 400531 2453434 55 400681 2453434 400381 2485734 57 400381 2485584 58 400381 2485434 59 400531 2485434 60 400681 2485434 424381 2453734 62 424381 2453584 63 424381 2453434 64 424531 2453434 65 424681 2453434 432381 2445734 51 56 61 66 407 229 409 449 Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Sa Pa Sa Pa Sa Pa Sa Pa Sa Pa Sa Pa Sa Pa TP Lào Cai TP Lào Cai TP Lào Cai TP Lào Cai TP Lào Cai Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Hầu Thào Hầu Thào Nậm Cang Nậm Cang Nậm Cang Nậm Cang Nậm Cang Vạn Hòa Vạn Hòa Vạn Hòa Vạn Hòa Vạn Hòa Sơn Thuỷ Sơn Thuỷ Sơn Thuỷ Sơn Thuỷ Sơn Thuỷ Tân Thượng NN NN RTNCG RTNCG RTNCG RTNCG RTNCG NN RTNCG RTG HG DT RTG HG HG HG HG RTG Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn 424381 2437734 Lào Cai Văn Bàn 424381 2437584 Lào Cai Văn Bàn 424381 2437434 Lào Cai Văn Bàn 67 432381 2445584 68 432381 2445434 69 432531 2445434 70 432681 2445434 392381 2437734 72 392381 2437584 73 392381 2437434 74 392531 2437434 75 392681 2437434 408381 2437734 77 408381 2437584 78 408381 2437434 79 408531 2437434 80 408681 2437434 82 83 71 76 81 491 493 495 Tân Thượng Tân Thượng Tân Thượng Tân Thượng HG HG HG HG Nậm Xé HG Nậm Xé HG Nậm Xé HG Nậm Xé HG Nậm Xé HG Thẩm Dương Thẩm Dương Thẩm Dương Thẩm Dương Thẩm Dương Khánh Yên Trung Khánh Yên Trung Khánh Yên Trung NN RTNCG HG HG HG HG HG HG Văn Bàn 84 424531 2437434 Lào Cai 85 424681 2437434 Lào Cai Văn Bàn 400381 2429734 87 400381 2429584 88 400381 2429434 89 400531 2429434 90 400681 2429434 416381 2429734 92 416381 2429584 93 416381 2429434 94 416531 2429434 95 416681 2429434 424381 2445734 97 424381 2445584 98 424381 2445434 99 424531 2445434 100 424681 2445434 Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Lào Cai Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn Văn Bàn 86 91 96 535 536 448 Khánh Yên Trung Khánh Yên Trung HG HG Nậm Xé RTNCG Nậm Xé HG Nậm Xé HG Nậm Xé HG Nậm Xé HG Dương Quỳ Dương Quỳ Dương Quỳ Dương Quỳ Dương Quỳ Sơn Thủy Sơn Thủy Sơn Thủy Sơn Thủy Sơn Thủy HG HG HG HG HG NN DT DT DT HG ... nghiên cứu 24 2. 1.1 Mục tiêu tổng quát 24 2. 1 .2 Mục tiêu cụ thể 24 2. 2 Đối tượng phạm vi nghiên cứu nghiên cứu 24 2. 2.1 Đối tượng nghiên cứu 24 ... chất lượng giải đoán ảnh viễn thám 2. 2 Đối tƣợng phạm vi nghiên cứu nghiên cứu 2. 2.1 Đối tượng nghiên cứu - Ảnh vệ tinh Landsat8 , Sentinel2 khu vực nghiên cứu - Lớp phủ trạng rừng theo bảng khóa... 29 Hình 2 .8: Tính tốn số ảnh vệ tinh 29 Hình 2. 9: Xuất ảnh vệ tinh sau tiền xử lý 30 Hình 2. 10: Ảnh vệ tinh Landsat khu vực nghiên cứu tiền xử lý 30 Hình 2. 11: Ảnh vệ tinh Sentinel