(Luận văn thạc sĩ) ứng dụng công nghệ viễn thám và gis giám sát độ che phủ thực vật tại khu vực mỏ than na dương, huyện lộc bình, tỉnh lạng sơn

BỘ TÀI NGUYÊN VÀ MÔI TRƢỜNG TRƢỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƢỜNG HÀ NỘI LUẬN VĂN THẠC SĨ ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM VÀ GIS GIÁM SÁT ĐỘ CHE PHỦ THỰC VẬT TẠI KHU VỰC MỎ THAN NA DƢƠNG, HUYỆN LỘC BÌNH, TỈNH LẠNG SƠN CHUYÊN NGÀNH: KHOA HỌC MÔI TRƢỜNG NGUYỄN DUY ANH HÀ NỘI, NĂM 2019 Luan van BỘ TÀI NGUYÊN VÀ MÔI TRƢỜNG TRƢỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƢỜNG HÀ NỘI LUẬN VĂN THẠC SĨ ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM VÀ GIS GIÁM SÁT ĐỘ CHE PHỦ THỰC VẬT TẠI KHU VỰC MỎ THAN NA DƢƠNG, HUYỆN LỘC BÌNH, TỈNH LẠNG SƠN NGUYỄN DUY ANH CHUYÊN NGÀNH: KHOA HỌC MÔI TRƢỜNG MÃ SỐ: 8440301 NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS.TS HOÀNG ANH HUY TS NGUYỄN TIẾN THÀNH HÀ NỘI, NĂM 2019 Hà Nội - Năm 20 Luan van CƠNG TRÌNH ĐƢỢC HỒN THÀNH TẠI TRƢỜNG ĐẠI HỌC TÀI NGUN VÀ MƠI TRƢỜNG HÀ NỘI Cán hƣớng dẫn chính: PGS.TS Hồng Anh Huy TS Nguyễn Tiến Thành Cán chấm phản biện 1: TS Phạm Thị Mai Thảo Cán chấm phản biện 2: TS Phạm Thị Việt Anh Luận văn thạc sĩ đƣợc bảo vệ tại: HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ TRƢỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƢỜNG HÀ NỘI Ngày tháng năm 2019 Luan van THÔNG TIN LUẬN VĂN Họ tên học viên: NGUYỄN DUY ANH Lớp: CH2B.MT Khóa: 2016-2018 Cán hƣớng dẫn: PGS.TS Hoàng Anh Huy TS Nguyễn Tiến Thành Tên đề tài: Ứng dụng công nghệ viễn thám GIS giám sát độ che phủ thực vật khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn Tóm tắt luận văn: Độ che phủ thực vật nhiều yếu tố quan trọng đƣợc sử dụng để đánh giá q trình tự nhiên nhƣ xói lở, trƣợt lở, lũ lụt nhƣ tốc độ phá hủy môi trƣờng tự nhiên hoạt động nhân sinh Mục tiêu nghiên cứu giám sát đƣợc biến động độ che phủ thực vật khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn từ tƣ liệu ảnh vệ tinh LANDSAT TM VÀ OLI Trên sở số NDVI, phƣơng pháp phân loại ảnh, ứng dụng mơ hình phân giải pixel hỗn hợp tuyến tính xác định độ che phủ thực vật (FVC), sau tiến hành đánh giá biến động độ che phủ tực vật Kết nghiên cứu cho thấy: Tại khu vực mỏ than Na Dƣơng, diện tích có độ che phủ thực vật (FVC) thƣa thớt từ dƣới 10% giảm mạnh từ 26,78 xuống 13,37 mạnh từ 19,49 , diện tích có độ che phủ thực vật từ dƣới 20% giảm xuống 13,38 ảnh hƣởng hoạt động khai thác mỏ than giai đoạn 1986 – 2015 Từ kết nghiên cứu kết luận: Ứng dụng phƣơng pháp viễn thám, sử dụng tƣ liệu ảnh vệ tinh LANDSAT giúp việc đánh giá, giám sát biến động độ che phủ thực vật cách hiệu quả, nhanh chóng tiết kiệm chi phí Luan van i LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi Kết nghiên cứu luận văn đóng góp tiêng dựa số liệu thu thập, kết nghiên cứu kế thừa cơng trình khoa học khác đƣợc trích dẫn theo quy định Nếu luận văn có chép từ cơng trình khoa học khác, tác giả xin hoàn toàn chịu trách nhiệm Học viên Nguyễn Duy Anh Luan van ii LỜI CẢM ƠN Trong trình thực luận văn, tơi nhận đƣợc quan tâm, giúp đỡ nhiều cá nhân quan đơn vị Nay luận văn hoàn thành, tơi xin bày tỏ lịng biết ơn chân thành, sâu sắc tới: Tôi xin đƣợc bày tỏ sử cảm ơn trân trọng tới PGS.TS Hồng Anh Huy tận tình hƣớng dẫn, giúp đỡ tơi suốt q trình hồn thành luận văn tốt nghiệp TS Nguyễn Tiến Thành, ngƣời tận tình hƣớng dẫn tạo điều kiện thuận lợi, ý kiến đóng góp sâu sắc để tơi hồn thành luận văn tốt nghiệp Tôi xin gửi lời cảm ơn đến thầy, cô giáo Khoa Môi trƣờng, Trƣờng Đại học Tài nguyên Mơi trƣờng Hà Nội nhiệt tình giúp đỡ tơi q trình hồn thành luận văn Cuối cùng, xin chân thành cảm ơn đến gia đình, cán bộ, đồng nghiệp bạn bè tạo điều kiện tốt mặt cho suốt trình thực đề tài Mặc dù có nhiều cố gắng để thực luận văn cách hoàn chỉnh nhất, song thời gian kinh nghiệm cịn hạn chế nên khơng tránh khỏi thiếu sót, mong nhận đƣợc đóng ghóp ý kiến Thầy, Cô để luận văn đƣợc hồn thiện Một lần tơi xin trân trọng cảm ơn! Hà Nội, ngày tháng năm 2019 Học viên Nguyễn Duy Anh Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an iii MỤC LỤC LỜI CAM ĐOAN .1 LỜI CẢM ƠN ii DANH MỤC VIẾT TẮT v DANH MỤC BẢNG vi DANH MỤC HÌNH vii MỞ ĐẦU .1 Tính cấp thiết đề tài nghiên cứu .1 Mục tiêu nghiên cứu Nội dung nghiên cứu .2 CHƢƠNG 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Tổng quan viễn thám .3 1.1.1 Khái niệm viễn thám 1.1.2 Các thành phần hệ thống viễn thám 1.1.4 Ứng dụng công nghệ viễn thám 1.1.5 Tìm hiểu vệ tinh viễn thám Landsat 10 1.2 Tổng quan hệ thống thông tin địa lý (GIS) 18 1.2.1 Định nghĩa hệ thống thông tin địa lý 18 1.2.2 Chức GIS 18 1.3 Tổng quan độ che phủ thực vật .19 1.3.1 Khái niệm độ che phủ thực vật .19 1.3.2 Các phƣơng pháp xác định độ che phủ thực vật 20 1.4 Tổng quan khu vực nghiên cứu .20 1.4.1 Điều kiện tự nhiên kinh tế - xã hội huyện Lộc Bình .20 1.4.2 Khái quát chung Mỏ than Na Dƣơng 33 CHƢƠNG 2: ĐỐI TƢỢNG, PHẠM VI VÀ PHƢƠNG PHÁP NGHIÊN CỨU .39 2.1 Đối tƣợng, phạm vi nghiên cứu 39 2.2 Phƣơng pháp nghiên cứu 39 2.2.1 Tiền xử lý ảnh vệ tinh 39 2.2.2 Hiệu chỉnh xạ 39 2.2.3 Hiệu chỉnh hình học 41 2.2.4 Phân loại ảnh 43 2.2.5 Xác định độ che phủ thực vật 47 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an iv 2.2.6 Đánh giá biến động độ che phủ thực vật 49 CHƢƠNG 3: KẾT QUẢ VÀ THẢO LUẬN 51 3.1 Tƣ liệu, tài liệu phần mềm sử dụng .51 3.1.1 Tƣ liệu, tài liệu sử dụng 51 3.1.2 Phần mềm sử dụng 52 3.2 Xử lý ảnh vệ tinh phần mềm .53 3.2.1 Sử dụng phần mềm Envi xử lý ảnh vệ tinh Landsat .53 3.2.2 Sử dụng phần mềm ArcMap để đánh giá biến động biên tập đồ 69 3.3 Kết xác định độ che phủ thực vật qua hai thời kỳ .73 3.3.1 Độ che phủ phủ thực vật huyện Lộc Bình, tỉnh Lạng Sơn năm 1986 .73 3.3.2 Độ che phủ thực vật khu vực mỏ than Na Dƣơng năm 1986 74 3.3.3 Độ che phủ thực vật huyện Lộc Bình, tỉnh Lạng Sơn năm 2015 76 3.3.4 Độ che phủ thực vật khu vực mỏ than Na Dƣơng năm 2015 77 3.4 Đánh giá biến động độ che phủ thực vật qua hai thời kỳ 79 3.4.1 Đánh giá biến động độ che phủ thực vật huyện Lộc Bình qua hai thời kỳ 79 3.4.1 Đánh giá biến động độ che phủ thực vật khu vực mỏ than Na Dƣơng qua hai thời kỳ 81 3.5 Đề xuất quy trình giám sát độ che phủ thực vật khu vực nghiên cứu 83 KẾT LUẬN, KIẾN NGHỊ 85 KẾT LUẬN .85 KIẾN NGHỊ 85 TÀI LIỆU THAM KHẢO 86 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an v DANH MỤC VIẾT TẮT FVC Độ che phủ thực vật (Fractional Vegetation Cover) GIS Hệ thống thông tin địa lý (Geographic Information System) LSMM Mơ hình hỗn hợp tuyến tính (Linear spectral mixture model) NDVI Chỉ số thực vật (Normalized Difference Vegetation Index) kỹ thuật vệ tinh thăm dò Trái đất (Earth Resource Technology ERTS Sattellite) OLI Bộ thu nhận ảnh mặt đất (Operational Land Imager) LDCM Vệ tinh Landsat (Landsat Data Continuity Mission) TIRS Bộ cảm biến hồng ngoại nhiệt (Thermal Infrared Sensor) ETM+ Bộ cảm biến Landsat (Enhanced Thematic Mapper Plus) MSS Bộ quét đa phổ (Multispectral Scanner) Chỉ số thực vật nhỏ Chỉ số thực vật lớn Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an vi DANH MỤC BẢNG Bảng 1.1: Các hệ vệ tinh Landsat 10 Bảng 1.2 : Một số thông tin kênh phổ ảnh Landsat TM .11 Bảng 1.3 : Thông tin kênh phổ ảnh Landsat ETM+ 12 Bảng 1.4 : Một số thông tin kênh phổ ảnh Landsat 13 Bảng 1.5 : Khả ứng dụng kênh phổ ảnh Landsat TM 14 Bảng 1.6 : Khả ứng dụng kênh phổ ảnh Landsat ETM+ .15 Bảng 1.7 : Khả ứng dụng kênh phổ ảnh Landsat .15 Bảng 2.1: Hệ thống phân loại thực phủ khu vực nghiên cứu .44 Bảng 3.1: Bảng thành lập khóa giải đốn ảnh vệ tinh năm 1986 năm 2015 61 Bảng 3.2 : Bảng thống kê số thực vật NDVI lớp phủ 65 Bảng 3.3 : Bảng thống kê độ che phủ thực vật khu vực huyện Lộc Bình năm 1986 73 Bảng 3.4 : Bảng thống kê độ che phủ thực vật khu vực mỏ than Na Dƣơng năm 1986 75 Bảng 3.5 : Bảng thống kê độ che phủ thực vật khu vực huyện Lộc Bình năm 201576 Bảng 3.6 : Bảng thống kê độ che phủ thực vật khu vực mỏ than Na Dƣơng năm 2015 78 Bảng 3.7 : Diện tích biến động độ che phủ thực vật huyện Lộc Bình 79 năm 1986 2015 79 Bảng 3.8 : Diện tích biến động độ che phủ thực vật khu vực mỏ than Na Dƣơng năm 1986 2015 82 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 75 Ta rút bảng sau: Bảng 3.4 : Bảng thống kê độ che phủ thực vật khu vực mỏ than Na Dương năm 1986 Độ che phủ TV_FVC(%) Diện tích ( Diện tích tích lũy ( ) ) Tỷ lệ (%) Tỷ lệ % tích lũy 00 – 10 26,78 26,78 36,43 36,43 10 – 20 19,09 45,87 25,97 62,40 20 – 30 12,44 58,31 16,92 79,32 30 – 40 7,07 65,38 9,62 88,94 40 – 50 3,87 69,25 5,27 94,21 50 - 60 2,36 71,61 3,21 97,42 60 - 70 1,13 72,75 1,54 98,96 70 - 80 0,49 73,24 0,67 99,63 80 -90 0,18 73,42 0,25 99,88 90 -100 0,09 73,51 0,12 100 Tổng 73,51 100 Từ Bảng, nhận xét thấy: - Diện tích có độ che phủ thực vật (FVC) thƣa thớt dƣới 20% chiếm đến 45,85 tổng diện tích khu vực nghiên cứu 73,51 - FVC dƣới 40% chiếm 65,38 (đạt 62,43%) (đạt 88,94%) - Một số khu vực có mật độ che phủ cao 60% đến 80% chiếm diện tích tƣơng đối nhỏ 1,62 (đạt 2,21%) - Khu vực có độ che phủ dày đặc từ 80% đến 100% chiếm diện tích nhỏ 0,28 (chỉ đạt 0,37%) Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 76 3.3.3 Độ che phủ thực vật huyện Lộc Bình, tỉnh Lạng Sơn năm 2015 Hình 3.32 : Bản đồ độ che phủ thực vật huyện Lộc Bình năm 2015 Ta rút bảng sau: Bảng 3.5 : Bảng thống kê độ che phủ thực vật khu vực huyện Lộc Bình năm 2015 Độ che phủ TV_FVC(%) Diện tích ( ) Diện tích tích lũy ( ) 00 – 10 7,78 107,72 144,96 169,32 160,33 141,36 99,53 69,48 23,26 5,59 998,34 7,78 179,50 324,46 493,78 654,11 795,48 895,01 964,50 987,76 993,35 10 – 20 20 – 30 30 – 40 40 – 50 50 - 60 60 - 70 70 - 80 80 -90 90 -100 Tổng Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van Tỷ lệ (%) 7,19 10,79 14,52 16,96 16,06 14,16 9,97 6,96 2,33 0,56 100 Tỷ lệ % tích lũy 7,69 18,48 33,00 49,96 66,02 80,18 90,15 97,11 99,44 100 C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 77 Từ Bảng, nhận xét thấy: - Diện tích có độ che phủ thực vật (FVC) thƣa thớt dƣới 20% chiếm đến 179,50 tổng diện tích khu vực nghiên cứu 998.34 - FVC dƣới 40% chiếm 493,78 (đạt 18,48%) (đạt 49,96%) - Một số khu vực có mật độ che phủ cao 60% đến 80% chiếm diện tích tƣơng đối nhỏ 169,01 (đạt 9,29%) - Khu vực có độ che phủ dày đặc từ 80% đến 100% chiếm diện tích nhỏ 28,85 (chỉ đạt 2,89%) 3.3.4 Độ che phủ thực vật khu vực mỏ than Na Dương năm 2015 Hình 3.33: Bản đồ độ che phủ thực vật khu vực mỏ than Na Dương năm 2015 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 78 Ta rút bảng sau: Bảng 3.6 : Bảng thống kê độ che phủ thực vật khu vực mỏ than Na Dương năm 2015 Độ che phủ TV_FVC(%) Diện tích ( ) Diện tích tích lũy ( ) Tỷ lệ (%) Tỷ lệ % tích lũy 00 – 10 13,37 13,37 18,19 18,19 10 – 20 13,38 26,75 18,20 36,39 20 – 30 12,95 39,69 17,61 54,00 30 – 40 12,39 52,08 16,85 70,85 40 – 50 7,89 59,98 10,74 81,59 50 - 60 5,33 65,30 7,25 88,84 60 - 70 4,18 69,48 5,68 94,52 70 - 80 2,96 72,44 4,02 98,54 80 -90 0,91 73,35 1,24 99,78 90 -100 0,16 73,51 0,22 100 Tổng 73,51 100 Từ Bảng, nhận xét thấy: - Diện tích có độ che phủ thực vật (FVC) thƣa thớt dƣới 20% chiếm đến 26,75 tổng diện tích khu vực nghiên cứu 73,51 - FVC dƣới 40% chiếm 52,08 (đạt 36,39%) (đạt 70,85%) - Một số khu vực có mật độ che phủ cao 60% đến 80% chiếm diện tích tƣơng đối nhỏ 7,14 (đạt 9,70%) - Khu vực có độ che phủ dày đặc từ 80% đến 100% chiếm diện tích nhỏ 1,07 (chỉ đạt 1,46%) Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 79 3.4 Đánh giá biến động độ che phủ thực vật qua hai thời kỳ 3.4.1 Đánh giá biến động độ che phủ thực vật huyện Lộc Bình qua hai thời kỳ Hình 3.34 : Bản đồ biến động độ che phủ thực vật huyện Lộc Bình giai đoạn 1986 – 2015 Ta có đƣợc kết sau: Bảng 3.7 : Diện tích biến động độ che phủ thực vật huyện Lộc Bình năm 1986 2015 FVC Diện tích năm 2015 Diện tích năm 1986 % Biến Động % % 00 – 10 7,78 7,19 163,3 16,34 -155,5 -9,15 10 – 20 107,72 10,79 188,29 18,86 -80,57 -8,07 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 80 20 – 30 144,96 14,52 180,9 18,12 -35,94 -3,6 30 – 40 169,32 16,96 142,46 14,27 +26,86 +2,69 40 – 50 160,33 16,06 122 12,22 +38,33 +3,84 50 - 60 141,36 14,16 88,55 8,87 +52,81 +5,29 60 - 70 99,53 9,97 60,7 6,08 +38,83 +3,89 70 - 80 69,48 6,96 33,74 3,38 +35,74 +3,58 80 -90 23,26 2,33 14,08 1,41 +9,18 +0,92 90 -100 5,59 0,56 4,49 0,45 +1,1 +0,11 Trong : + Giá trị tăng đối tƣợng, - Giá trị giảm đối tƣợng Từ kết bảng 3.7 phân bố biến động hình 3.34 thấy rõ biến động độ che phủ thực vật Năm 2015, diện tích có độ che phủ thực vật (FVC) thƣa thớt dƣới 10% 20% có giảm mạnh lần lƣợt giảm từ 163,3 xuống tới 7,78 từ 188,92 xuống tới 107,72 , tập trung xã, thị trấn thuộc khu vực trung tâm huyện nhƣ: Na Dƣơng, Lộc Bình, n Khối, Sàn Viên, Nhƣợng Ban, Mẫu Sơn Có thể thấy FVC từ 30% có dấu hiệu tăng dần kể từ sau năm 1986; đặc biệt khu vực có độ che phủ trung bình từ 50% đến 60% có dấu hiệu tăng mạnh, tập trung xã thuộc khu mực miền núi nhƣ Lợi Bảo, Tĩnh Bắc, Tú Mịch Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 81 3.4.1 Đánh giá biến động độ che phủ thực vật khu vực mỏ than Na Dƣơng qua hai thời kỳ Hình 3.35 : Bản đồ biến động độ che phủ thực vật khu vực mỏ than Na Dương giai đoạn 1986 – 2015 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 82 Ta có đƣợc kết sau: Bảng 3.8 : Diện tích biến động độ che phủ thực vật khu vực mỏ than Na Dương năm 1986 2015 FVC Diện tích năm 2015 Diện tích năm 1986 % Biến Động % % 00 – 10 13,37 18,19 26,78 36,43 -13,41 -18,24 10 – 20 13,38 18,2 19,09 25,97 -5,71 -7,77 20 – 30 12,95 17,61 12,44 16,92 +0,51 +0,69 30 – 40 12,39 16,85 7,07 9,62 +5,32 +7,23 40 – 50 7,89 10,74 3,87 5,27 +4,02 +5,47 50 - 60 5,33 7,25 2,36 3,21 +2,97 +4,04 60 - 70 4,18 5,68 1,13 1,54 +3,05 +4,14 70 - 80 2,96 4,02 0,49 0,67 +2,47 +3,35 80 -90 0,91 1,24 0,18 0,25 +0,73 +0,99 90 -100 0,16 0,22 0,09 0,12 +0,07 +0,1 Trong : + Giá trị tăng đối tƣợng, - Giá trị giảm đối tƣợng Từ kết bảng 3.8 thấy rõ biến động độ che phủ thực vật khu vực mỏ than Na Dƣơng Năm 2015, diện tích có độ che phủ thực vật (FVC) thƣa thớt từ dƣới 10% giảm mạnh từ 26,78 xuống 13,37 độ che phủ thực vật từ dƣới 20% giảm mạnh từ 19,49 ảnh hƣởng hoạt động khai thác mỏ than Na Dƣơng Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van , diện tích có xuống 13,38 C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 83 3.5 Đề xuất quy trình giám sát độ che phủ thực vật khu vực nghiên cứu Tài liệu, báo cáo Thông tin thực địa Ảnh Landast Xác định điểm khảo sát Thu thập thông tin thực địa B1 Tiền xử lý ảnh B2.1 Cắt ảnh Xây dựng mẫu B2.1 Phân loại B2.3 NDVI NDVIv NDVIs Độ che phủ thực vật B2.4 B2.5 Chuyển sang Vector B2.6 Bản đồ 1986 Bản đồ 2015 Chồng ghép đồ B3.1 Đánh giá biến động Thống kê số liệu Hình 3.36 : Quy trình giám sát độ che phủ thực vật Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van B3.2 C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 84 B1: Chuẩn bị tài liệu nghiên cứu B2.1: Tiền sử lý ảnh ( hiệu chỉnh xạ, hiệu chỉnh hình học, tăng cƣờng chất lƣợng ảnh), cắt ảnh theo khu vực nghiên cứu B2.2: Kết hợp tài liệu, báo cáo, số liệu khảo sát thực địa để xác định dạng lớp phủ thực vật, chọn mẫu ảnh tƣơng ứng với loại lớp phủ Từ kết phân tích đánh giá đội xác nhóm mẫu, loại bỏ mẫu có độ tin cậy thấp B2.3: Phân loại tự động, có kiểm định ảnh Landsat TM Landsat theo phƣơng pháp Maximumlikelihood B2.4: Sử dụng cơng cụ Masking kết hợp với kết tính số thực vật NDVI (tính tốn theo giá trị phản xạ kênh phổ RED, NIR) để chiết tách số thực vật NDVIs, NDVIv loại lớp phủ B2.5: Áp dụng phƣơng pháp mơ hình phân giải pixel hỗn hợp tuyến tính (LSMM) tính tốn độ che phủ thực vật B2.6: Chuyển kết tính tốn vào phần mềm ArcGIS, thành lập đồ độ che phủ thực vật năm 1986 2015 B3.1: Chồng ghép đồ độ che phủ thực vật năm 1986-2015 B3.2: Đánh giá biến động, thống kê số liệu phục vụ công tác giám sát Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 85 KẾT LUẬN, KIẾN NGHỊ KẾT LUẬN Sau hoàn thành đề tài “Ứng dụng công nghệ viễn thám GIS giám sát độ che phủ thực vật khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn”, em xin rút số kết luận sau:  Sử dụng theo phƣơng pháp viễn thám, dùng ảnh vệ tinh có độ phân giải cao, đa thời gian để tính tốn đánh giá biến động độ che phủ thực vật cho khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn phƣơng pháp hiệu Tƣ liệu ảnh có tính cập nhập, chiết tách thơng tin lớp thảm thực vật khơng có điều kiện thực tế thƣờng xuyên  Dựa vào ảnh viễn thám, kết hợp tài liệu tham khảo sử dụng phần mềm ENVI ta tính đƣợc số thực vật NDVI cho khu vực mỏ than Na Dƣơng, huyện Lộc Binh, tỉnh Lạng Sơn dựa vào để thành lập đồ biến động độ che phủ thực vật giai đoạn 1986 – 2015  Qua đồ biến động, ta thấy độ che phủ thực vật khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn có dấu giảm mạnh hoạt động khai thác mỏ than Do đó, cần tăng cƣờng, thực có hiệu biện pháp cải tạo, phục hồi môi trƣờng sau kết thúc hoạt động khai thác KIẾN NGHỊ Tƣ liệu ảnh vệ tinh Landsat công cụ hiệu chiết tách thông tin độ che phủ thực vật Để đạt đƣợc độ xác cao hơn, nên sử dụng loại ảnh khác có độ phân giải cao Do hạn chế mặt thời gian tƣ liệu, đồ án sử dụng cảnh ảnh Landsat giai đoạn 1986 – 2015 Hƣớng phát triển đề tài sử dụng số lƣợng lớn tƣ liệu ảnh vệ tinh Landsat để đánh giá cách toàn diện khách quan không khu vực mỏ than Na Dƣơng, huyện Lộc Bình mà cịn nhiều khu vực khác Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 86 TÀI LIỆU THAM KHẢO Nguyễn Khắc Thời (2012), Giáo trình Viễn thám, Nhà xuất Đại học Nông nghiệp Hà Nội Nguyễn Văn Đài (2002), Cơ sở viễn thám, Giáo trình bậc đại học, ĐHQG Hà Nội Trịnh Lê Hùng, Vũ Danh Tuyên, Phạm Thị Thƣơng Huyền (2013), Cơ sở viễn thám, Giáo trình đại học, Đại học Tài nguyên Môi trƣờng Hà Nội Trịnh Lê Hùng, Vũ Danh Tuyên (2014), Ứng dụng công nghệ viễn thám nghiên cứu, giám sát tài ngun mơi trường, Giáo trình bậc đại học, Đại học Tài nguyên Môi trƣờng Hà Nội Hoàng Anh Huy (2016), Ứng dụng ảnh vệ tinh Landsat OLI xác định độ che phủ thực vật khu vực nội thành Hà Nội, Tạp chí Khoa học Đại học Quốc gia Hà Nội Hoàng Anh Huy (2017), Đánh giá biến động độ che phủ thực vật số khu vực đô thị ven đô Hà Nội từ tư liệu ảnh vệ tinh LANDSAT đa phổ đa thời gian, Đại học Tài nguyên Môi trƣờng Hà Nội Trung tâm quan trắc môi trƣờng tỉnh Lạng Sơn, Báo cáo giám sát môi trường định kỳ Quý IV mỏ than Na Dương xã Sàn Viên, huyện Lộc Bình, tỉnh Lạng Sơn Lu, D and Weng, Q (1986) Spectral mixture analysis of the urban landscape in Indianapolis city with Landsat ETM+ imagery Photogrammetric Engineering and Remote Sensing, 70, 1053-1062 SobrinoJ.A.,Jiménez-MozJ.C.,SịriaG.,etal(2008) Land Surface Emissivity Retrieval From Different VNIRand TIR sensors, IEEE Transactions on Geoscience and Remote Sensing, 46, 2, 316-326 10 Van der Meer, F 1999 Image classification through spectral unmixing In: Spatial Statistics for Remote Sensing, Stein, A., Van der Meer, F & Gorte, B (Eds.) Kluwer Academic Publishers, Dordrecht, pp 185-193 11 Z Zhang, M Ji (2008), Surface urban heat island in Shanghai, China: Examining the relationship between land surface temperature and impervious surface fraction derived from LANDSAT ETM+ imagery, The International Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an 87 archives of the Photogrammetric, Remote sensing and spatial information sciences, Vol 37, Beijing, pp 601 - 607 12 Hoffmann, W A., & Jackson, R (2000) Vegetation-climate feedbacks in the conversion of tropical savanna to grassland Journal of Climate, 13, 1593– 1602 13 Ward, R C., & Robinson, M (2000) Principles of Hydrology (4th edition) McGraw hill pp 450 14 Gutman, G.; Ignatov, A The derivation of the green vegetation fraction from NOAA/AVHRR data for use in numerical weather prediction models, International Journal of Remote Sensing 1998, 19 (8), 1533-1543 15 Zeng, X., Dickinson, R E., Walker, A., & Shaikh, M (2000) Derivation and evaluation of global 1-km fractional vegetation cover data for land modeling Journal of Applied Meteorology, 39, 826–839 16 Avissar, R; Pielke, R A A parameterization of heterogeneous land surfaces for atmospheric numerical models and its impact on regional meteorology Monthly Weather Review 1989, 117, 2113-2136 17 Trimble, S W Geomorphic effects of vegetation cover and management: some time and space considerations in prediction of erosion and sediment yield, in Vegetation and Erosion, edited by J B Thornes, London, John Wiley & Sons, 1990, pp 55-66 18 http://earthexplorer.usgs.gov/ 19 Juan C Jiménez-Muñoz, José A Sobrino, Antonio Plaza, Luis Guanter, José Moreno and Pablo Martínez Comparison Between Fractional Vegetation Cover Retrievals from Vegetation Indices and Spectral Mixture Analysis: Case Study of PROBA/CHRIS Data Over an Agricultural Area Sensors, 2009, 9, 768–793 20 Ying Li, Hong Wang and Xiao Bing Li Fractional Vegetation Cover Estimation Based on an Improved Selective Endmember Spectral Mixture Model PLoS One 2015 Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an LÝ LỊCH TRÍCH NGANG Họ tên: NGUYỄN DUY ANH Ngày tháng năm sinh: 27/09/1994 Nơi sinh: Hà Nội Địa liên lạc: Tổ 29, phƣờng Bồ Đề, quận Long Biên, Hà Nội Quá trình đào tạo: Đại học - Hệ đào tạo: Chính quy - Thời gian đào tạo: từ 08/2012 đến 06/2016 - Trƣờng đào tạo: Đại học Tài nguyên Môi trƣờng Hà Nội - Ngành học: Quản lý Tài nguyên Môi trƣờng - Bằng tốt nghiệp đạt loại: TB Khá Thạc sĩ - Hệ đào tạo: Chính quy - Thời gian đào tạo: từ năm 2016 đến năm 2018 - Chuyên ngành học: Khoa học Môi trƣờng - Tên luận văn: Ứng dụng công nghệ viễn thám GIS giám sát độ che phủ thực vật khu vực mỏ than Na Dƣơng, huyện Lộc Bình, tỉnh Lạng Sơn - Ngƣời hƣớng dẫn Khoa học: PGS.TS Hoàng Anh Huy TS Nguyễn Tiến Thành Q trình cơng tác: Thời gian Nơi cơng tác Cơng việc đảm nhận 11/2016 – 02/2019 Phịng Tài nguyên Môi trƣờng Ủy ban nhân dân quận Long Biên Nhân viên Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn Luan van C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn