Sử dụng ảnh vệ tinh landsat 8 trong xây dựng bản đồ phủ mặt đất huyện tam đảo tỉnh vĩnh phúc

TRƯỜNG ĐẠI HỌC LÂM NGHIỆP VIỆT NAM KHOA QUẢN LÝ TÀI NGUN RỪNG & MƠI TRƯỜNG KHĨA LUẬN TỐT NGHIỆP SỬ DỤNG ẢNH VỆ TINH LANDSAT TRONG XÂY DỰNG BẢN ĐỒ LỚP PHỦ MẶT ĐẤT HUYỆN TAM ĐẢO, TỈNH VĨNH PHÚC NGÀNH: QUẢN LÝ TÀI NGUYÊN THIÊN NHIÊN MÃ NGÀNH: 7908532 Giảng viên hướng dẫn : ThS Lê Thái Sơn Sinh viên thực : Nguyễn Huy Hoàng Mã sinh viên : 1553090858 : 60 – QLTNTN (c) Lớp Khóa học : 2015 - 2019 Hà Nội, tháng năm 2020 LỜI CẢM ƠN Em xin gửi lời cảm ơn chân thành tới thầy cô khoa Quản lý tài nguyên rừng & môi trường, bạn lớp K60_ Quản lý tài nguyên thiên nhiên (Chuẩn) em sinh viên lớp K61_ Quản lý tài nguyên thiên nhiên (Chuẩn) tạo điều kiện giúp đỡ em trình học tập Đặc biệt, em xin gửi lời cảm ơn sâu sắc tới thầy ThS Lê Thái Sơn, giảng viên Bộ môn Quản lý môi trường – Khoa Quản lý tài nguyên rừng môi trường – Đại học Lâm nghiệp Việt Nam, người hướng dẫn, giúp đỡ em nhiều trình học tập q trình hồn thành khóa luận Em xin gửi lời cảm ơn tới gia đình có cơng ni dưỡng tạo điều kiện tốt để em học tập đến ngày hơm Cuối cùng, em xin gửi lời kính chúc sức khỏe đến tồn thể q thầy, trường Đại học Lâm nghiệp, chúc quý thầy, cô thành công công tác giảng dạy truyền đạt kiến thức bổ ích cho hệ sau chúng em Em xin chân thành cảm ơn! Sinh viên thực Nguyễn Huy Hồng i TĨM TẮT Khóa luận tốt nghiệp “ Sử dụng ảnh Landsat xây dựng đồ lớp phủ mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc” thực từ tháng 02/2020 đến tháng 05/2020 Phương pháp tiếp cận đề tài sử dụng công nghệ viễn thám kết hợp với GIS Nội dung đề tài nghiên cứu vấn đề sau:  Điều tra, nghiên cứu đặc điểm lớp phủ mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc  Thu thập liệu ảnh vệ tinh Landsat 8, số liệu thống kê  Tiến hành giải đoán ảnh  Thành lập đồ lớp phủ bề mặt đánh giá biến động qua năm Sau trình thực đề tài thu số kết sau:  Bản đồ lớp phủ mặt đất khu vực huyện Tam Đảo, tỉnh Vĩnh Phúc năm 2014 năm 2020 với loại: giao thông, đất rừng, mặt nước, đất nông nghiệp, đất trống khu dân cư  Bản đồ biến động lớp phủ mặt đất huyện Tam Đảo giai đoạn từ năm 2014 đến năm 2020 Với kết đạt được, thấy việc sử dụng cơng nghệ viễn thám GIS phương pháp có hiệu với độ xác cao, tiết kiệm chi phí việc phân loại phân tích biến động lớp phủ mặt đất ii MỤC LỤC LỜI CẢM ƠN i TÓM TẮT ii MỤC LỤC iii DANH MỤC TỪ VIẾT TẮT v DANH MỤC BẢNG vi DANH MỤC HÌNH vii Chương 1: ĐẶT VẤN ĐỀ 1.1 Tính cấp thiết đề tài 1.2 Cấu trúc đề tài Chương 2: TỔNG QUAN NGHIÊN CỨU 2.1 Khu vực nghiên cứu 2.1.1 Điều kiện tự nhiên 2.1.2 Điều kiện kinh tế - xã hội 2.2 Tổng quan viễn thám 2.2.1 Khái niệm 2.2.2 Nguyên lý hoạt động viễn thám 2.3 Tổng quan GIS 2.3.1 Định nghĩa 2.3.2 Chức GIS: 2.4 Các khái niệm 2.4.1 Lớp phủ mặt đất (Lớp thực phủ - Land cover) 2.4.2 Phân loại lớp phủ mặt đất 11 2.4.3 Khái quát đồ biến động lớp phủ bề mặt 13 2.5 Giới thiệu vệ tinh LANDSAT 15 Chương 3:VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 17 3.1 Dữ liệu thu thập 17 3.2 Mục tiêu nghiên cứu 17 3.2.1 Mục tiêu chung 17 iii 3.2.2 Mục tiêu cụ thể 17 3.3 Nội dung nghiên cứu 17 3.4 Phương pháp nghiên cứu 18 3.4.1 Khảo sát thực địa 18 3.4.2 Hệ thống phân loại lớp phủ khu vực nghiên cứu 21 3.4.3 Lựa chọn phương pháp phân loại ảnh 22 3.4.4 Xử lý liệu ảnh 24 a Gộp kênh ảnh 24 b Cắt ảnh 24 3.4.5 Giải đốn ảnh khóa giải đốn ảnh 24 3.4.6 Đánh giá độ xác xử lý ảnh sau phân loại 26 a Đánh giá độ xác sau phân loại ma trận sai số 26 b Xử lý ảnh sau phân loại 28 3.4.7 Thành lập đồ 28 3.4.8 Phương pháp thành lập đồ biến động lớp phủ mặt đất 28 Chương 4: KẾT QUẢ VÀ THẢO LUẬN 30 4.1 Thành lập đồ lớp phủ mặt đất 30 4.1.1 Kết phân loại ảnh viễn thám 30 4.1.2 Kết đánh giá độ xác hệ số Kappa 33 4.1.3 Bản đồ lớp phủ mặt đất 35 4.2 Bản đồ biến động lớp phủ bề mặt 38 4.3 Khả ứng dụng ảnh Landsat công nghệ GIS thực tiễn.41 4.4 Thảo luận 42 4.5 Kiến nghị 44 Chương 5: KẾT LUẬN – ĐỀ NGHỊ 45 5.1.Kết luận 45 5.2 Đề nghị 45 TÀI LIỆU THAM KHẢO PHỤ LỤC iv DANH MỤC TỪ VIẾT TẮT CORINE: Coordination of information on the environment ĐPGKG: Độ phân giải không gian ERTS: Earth Resource Technology Sattellite FAOLCC: Food and Agriculture Organization Land Cover Classification GIS: Geographic Information System GPS: Global Position System LDCM: Landsat Data Continuity Mission NDVI: Normalized Difference Vegetation Index USGS: United States Geological Survey v DANH MỤC BẢNG Bảng 1: Hệ thống phân loại lớp phủ mặt đất huyện Tam Đảo 13 Bảng 2: Đặc trưng Bộ cảm ảnh vệ tinh LANDSAT 16 Bảng 1: Dữ liệu ảnh thu thập 17 Bảng 2: Thống kê số điểm mẫu loại lớp phủ 20 Bảng 3: Hệ thống phân loại lớp phủ khu vực nghiên cứu 22 Bảng 4: Khóa giải đốn cho khu vực nghiên cứu 25 Bảng 1: Kết phân loại ảnh huyện Tam Đảo năm 2014 32 Bảng 2: Kết phân loại ảnh huyện Tam Đảo năm 2020 33 Bảng 3: Kết ma trận sai số phân loại ảnh 34 Bảng 4: Thống kê biến động lớp phủ bề mặt 40 vi DANH MỤC HÌNH Hình 1: Khu vực nghiên cứu Hình 2: Nguyên lý thu nhận liệu viễn thám Hình 3: Sơ đồ tổng quát lớp phủ bề mặt đất 10 Hình 1: Phương pháp nghiên cứu 19 Hình 2: Các điểm mẫu khảo sát thực địa 20 Hình 1: Sơ đồ vị trí lấy mẫu giải đốn ảnh năm 2014 30 Hình 2: Sơ đồ vị trí lấy mẫu giải đốn ảnh năm 2020 30 Hình 3: Kết phân loại ảnh huyện Tam Đảo năm 2014 phương pháp Maximum Likelihood 31 Hình 4: Kết phân loại ảnh huyện Tam Đảo năm 2020 phương pháp Maximum Likelihood 32 Hình 5: Bản đồ phân loại lớp phủ bề mặt huyện Tam Đảo năm 2014 36 Hình 6: Bản đồ phân loại lớp phủ bề mặt huyện Tam Đảo năm 2020 37 Hình 7: Bản đồ biến động lớp phủ bề mặt huyện giai đoạn năm 2014-2020.38 Hình 8: Biểu đồ biến động lớp phủ bề mặt giai đoạn năm 2014-2020 39 vii Chương 1: ĐẶT VẤN ĐỀ 1.1 Tính cấp thiết đề tài Đất đai tài nguyên thiên nhiên vô quý giá, tư liệu sản xuất đặc biệt Đất đai môi trường sống người sinh vật, địa bàn phân bố dân cư, xây dựng công trình kinh tế, văn hố, an ninh quốc phịng Ngày nay, tăng dân số, phát triển đô thị, tăng trưởng kinh tế xã hội số vấn đề khác tác động lớn tới đất đai, đặc biệt nơi có tốc độ phát triển nhanh huyện Tam Đảo thuộc tỉnh Vĩnh Phúc Trước áp lực đó, đất đai lớp phủ mặt đất biến động không ngừng với phát triển xã hội Đây nguồn tài nguyên đặc biệt khai thác sử dụng làm tăng thêm mặt số lượng Do việc theo dõi, nghiên cứu, quản lý sử dụng loại tài nguyên cách hiệu hợp lý vấn đề quan trọng Huyện Tam Đảo huyện trọng điểm du lịch tỉnh Vĩnh Phúc Đặc biệt năm qua với chuyển mạnh mẽ theo hướng cơng nghiệp hóa – đại hóa chung đất nước, huyện thay đổi nhanh chóng theo hướng giảm diện tích đất nơng nghiệp, tăng diện tích đất phi nơng nghiệp, chủ yếu đất sử dụng vào mục đích du lịch – dịch vụ (Theo cổng thông tin – giao tiếp điện tử tỉnh Vĩnh Phúc) Với phát triển mạnh mẽ khoa học vũ trụ ảnh viễn thám xuất ngày tỏ rõ tính ưu việt cơng tác điều tra, quản lý tài nguyên Đặc biệt xuất nhiều loại tư liệu ảnh viễn thám: SPOT, LANDSAT, có độ phân giải khơng gian phân giải phổ cao Hơn nữa, tư liệu ảnh cịn có khả chụp lập thể, đặc biệt cập nhật thơng tin nhanh chóng thơng qua việc thu nhận ảnh vệ tinh nhiều thời điểm khác tạo thành ảnh đa thời gian dạng số, sản phẩm dễ dàng sử dụng phần mềm phân tích ảnh đại có khả tích hợp thuận tiện hệ thống thơng tin địa lý GIS Đặc biệt, việc phóng vệ tinh Landsat phóng vào ngày 11 tháng 02 năm 2013 mở hướng ứng dụng ảnh viễn thám Việt Nam Nhận thức tầm quan trọng thay đổi trình sử dụng đất với thay đổi chất lượng sống, đề tài nghiên cứu “Sử dụng ảnh Landsat xây dựng đồ lớp phủ mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc” thực 1.2 Cấu trúc đề tài Tồn đề tài trình bày chương:  Chương I: ĐẶT VẤN ĐỀ  Chương II: TỔNG QUAN NGHIÊN CỨU  Chương III: VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU  Chương IV: KẾT QUẢ VÀ THẢO LUẬN  Chương V: KẾT LUẬN – KIẾN NGHỊ Chương 5: KẾT LUẬN – ĐỀ NGHỊ 5.1.Kết luận Huyện Tam Đảo huyện thuộc Vĩnh Phúc có tốc độ thị hóa nhanh Cụ thể là: đất nơng nghiệp năm 2014 3.579,93 giảm xuống cịn 2.735,46 tính đến năm 2020, giảm 844,47 ha, diện tích đất rừng năm 2014 là16.683 giảm xuống cịn 14.939,37 tính đến năm 2020 Thay vào lớp đất nông nghiệp đất rừng khu dân cư Diện tích khu dân cư năm 2014 2.135,70 tăng lên 4.277,97 tính đến năm 2020, tăng 2.142,27 Công nghệ viễn thám kết hợp với GIS cho hiệu cao khách quan đánh giá biến đổi lớp phủ mặt đất Kết thực nghiệm rõ, việc kết hợp công nghệ viễn thám GIS hữu hiệu để xác định diện tích biến động, mức độ biến động phần xu hướng biến động đối tượng Đây sở giúp nhà hoạch định, quản lý tài nguyên – môi trường, quy hoạch sử dụng đất đánh giá xác trạng khu vực Nghiên cứu ứng dụng công nghệ viễn thám kết hợp với GIS để thành lập đồ lớp phủ bề mặt đất huyện Tam Đảo năm 2014 năm 2020 đồ biến động lớp phủ bề mặt huyện Tam Đảo giai đoạn năm 2014 – 2020, Bên cạnh đó, việc sử dụng liệu ảnh có độ phân giải chưa cao gây nhiều khó khăn q trình giải đốn Vì vậy, cần phải điều tra mẫu thực địa để kết giải đoán đạt độ tin cậy Việc sử dụng tư liệu viễn thám thành lập đồ lớp phủ tương đối đơn giản nhanh chóng Nếu đầu tư ứng dụng rộng rãi tiết kiệm chi phí, cơng sức, thời gian mà kết thu lại có hiệu cao 5.2 Đề nghị Do hạn chế kinh phí nên đề tài sử dụng ảnh miễn phí với độ phân giải trung bình, chất lượng khơng cao, đó, kết nhận chưa đạt độ xác cao Để đạt độ 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ế thời gian nguồn liệu nên đề tài thành lập đồ lớp phủ năm 2014, năm 2020 đồ biến động lớp phủ bề mặt giai đoạn năm 2014 – 2020 Để đạt kết có giá trị cao làm nguồn liệu cho nhà hoạch định sách, nên sử dụng nhiều ảnh nhiều thời điểm thu hẹp biên độ thời gian đánh giá biến động Nên kết hợp nhiều phương pháp phân loại loại liệu ảnh viễn thám khác để thực giải đoán nhằm đạt kết tốt 45 TÀI LIỆU THAM KHẢO Tài liệu Tiếng Việt: Lê Quang Tuấn, Hà Quý Quỳnh, Lê Minh Hạnh, Lê Xuân Cảnh, 2011, Ứng dụng công nghệ viễn thám hệ thống thông tin địa lý để nghiên cứu lớp phủ thực vật Vườn Quốc gia Tam Đảo, Viện Sinh thái Tài nguyên sinh vật Lê Văn Trung, 2010, Viễn Thám, NXB Đại Học Quốc Gia TPHCM Nguyễn Hải Hòa, Nguyễn Văn Quốc, 2017, Sử dụng ảnh viễn thám Landsat GIS xây dựng đồ biến động diện tích rừng vùng đệm Vườn Quốc gia Xuân Sơn, Trường Đại học Lâm nghiệp Nguyễn Khắc Thời, 2011, Giáo trình viễn thám Nguyễn Kim Lợi, Lê Đình Cảnh, Trần Thống Nhất, 2009, Hệ thống thông tin địa lý nâng cao, NXB Nông nghiệp Nguyễn Trường Sơn, 2008, Nghiên cứu sử dụng ảnh vệ tinh công nghệ GIS việc giám sát trạng tài nguyên rừng, Báo cáo khoa học, Trung tâm viễn thám quốc gia, Hà Nội Nguyễn Xuân Trung Hiếu, 2013, Ứng dụng viễn thám GIS thành lập đồ biến động loại thực phủ địa bàn thành phố Huế - tỉnh Thừa Thiên Huế, KLTN, ĐH Nông Lâm, TP HCM Tài liệu Website: http://tamdao.vinhphuc.gov.vn/Pages/Default.aspx http://tnmtvinhphuc.gov.vn/ https://earthexplorer.usgs.gov/ https://vi.wikipedia.org/wiki/Tam_%C4%90%E1%BA%A3o Phụ lục 01: Một số điểm mẫu đặc trưng STT Tọa độ 21⁰24’16,3”N 105⁰38’17,5”E 21⁰24’14,56”N 105⁰38’13,38”E Loại lớp phủ Hình ảnh Ghi Mặt nước Đập Xạ Hương, Minh Quang, Tam Đảo Đất trống Đập Xạ Hương, Minh Quang, Tam Đảo 21⁰25’37,0”N 105⁰35’21,2”E Giao thông ĐT302, Làng Mạ, Tam Quan, Tam Đảo 21⁰25’1,51”N 105⁰35’59,84”E Đất nông nghiệp Núc Thượng, Hồ Sơn, Tam Đảo 21⁰23’18,99”N 105⁰36’52,93”E Giao thông QL2B Km9 Hợp Châu, Tam Đảo 21⁰27’22,86”N 105⁰38’43,92”E Khu dân cư Thị trấn Tam Đảo, Tam Đảo Đất rừng Vĩnh Ninh, Đạo Trù, Tam Đảo 21⁰25’31,17”N 105⁰36’43,56”E Mặt nước Làng Hà, Hồ Sơn, Tam Đảo 21⁰29’47,63”N 105⁰31’0,29”E Đất nông nghiệp Yên Dương, Tam Đảo 10 21⁰27’2,57”N 105⁰33’36,55”E Khu dân cư Sơn Thanh, Đại Đình, Tam Đảo 21⁰31’47,50”N 105⁰30’43,56”E Quan Nội, Tam Quan, 11 21⁰25’27,02”N 105⁰34’49,74”E 12 21⁰25’40,86”N 105⁰36’13,17”E Mặt nước Đồng Bùa, Tam Quan, Tam Đảo 13 21⁰30’24,27”N 105⁰32’27,13”E Đất nông nghiệp Thị Tứ, Đạo Trù, Tam Đảo Đất trống Tam Đảo 14 21⁰25’5,53”N 105⁰34’48,33”E Đất rừng Quan Ngoại, Tam Quan, Tam Đảo 15 21⁰26’49,71”N 105⁰35’34,29”E Khu dân cư Đồng Hội, Đại Đình, Tam Đảo 16 21⁰24’13,06”N 105⁰37’35,91”E Đất trống Khu 9, Hợp Châu, Tam Đảo 17 21⁰31’44,61”N 105⁰31’29,78”E Mặt nước Vĩnh Ninh, Đạo Trù, Tam Đảo 18 21⁰23’52,79”N 105⁰36’51,48”E Đất trống Núc Hạ, Hồ Sơn, Tam Đảo 19 21⁰28’55,39”N 105⁰31’41,22”E Giao thông TL302, Bồ Lý, Tam Đảo 20 21⁰25’9,42”N 105⁰36’23,95”E Đất nông nghiệp Làng Hà, Hồ Sơn, Tam Đảo BẢNG I: Thống kê diện tích đất huyện Tam Đảo Nguồn:Quy hoạch, kế hoạch sử dụng đất huyện Tam Đảo – Sở Tài nguyên & Môi trường Vĩnh Phúc STT 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 2.1 2.2 2.3 2.4 CHỈ TIÊU SỬ DỤNG ĐẤT TỔNG DIỆN TÍCH TỰ NHIÊN Đất nơng nghiệp Đất trồng lúa Trong đó: Đất chuyên trồng lúa nước Đất trồng hàng năm khác Đất trồng lâu năm Đất rừng phòng hộ Đất rừng đặc dụng Đất rừng sản xuất Đất nuôi trồng thủy sản Đất nông nghiệp khác Đất phi nơng nghiệp Đất quốc phịng Đất an ninh Đất khu công nghiệp Đất thương mại, dịch vụ HIỆN TRẠNG NĂM 2014 Diện tích (ha) Cơ cấu (%) 23469,90 100,00 20359,34 86,75 2700,45 11,51 2700,45 11,51 764,07 3,26 1955,40 8,33 417,26 1,78 12364,57 52,68 1943,67 8,28 93,45 0,40 120,47 0,51 3084,53 13,14 249,59 1,06 217,98 0,93 103,62 0,44 30,72 0,13 KẾ HOẠCH NĂM 2020 Diện tích (ha) Cơ cấu (%) 23469,90 100,00 17822,67 75,94 2166,08 9,23 2166,08 9,23 430,91 1,84 1547,14 6,59 322,73 1,38 12276,82 52,31 773,09 3,29 98,45 0,42 207,45 0,88 5633,12 24,00 942,24 4,01 237,08 1,01 294,24 4,01 243,26 1,04 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13 2.14 2.15 2.16 2.17 2.18 2.19 2.20 2.21 2.22 Đất sở sản xuất phi nơng nghiệp Đất sử dụng cho hoạt động khống sản Đất phát triển hạ tầng cấp quốc gia, cấp tỉnh, cấp huyện, cấp xã Đất có di tích lịch sử - văn hóa Đất bãi thải, xử lý chất thải Đất nông thôn Đất đô thị Đất xây dựng trụ sở tổ chức quan Đất xây dựng trụ sở tổ chức nghiệp Đất sở tôn giáo Đất làm nghĩa trang, nghĩa địa Đất sản xuất vật liệu xây dựng Đất sinh hoạt cộng đồng Đất vui chơi, giải trí cơng cộng Đất sở tín ngưỡng Đất sơng, ngịi, kênh, rạch, suối Đất có mặt nước chun dùng Đất phi nơng nghiệp khác Đất chưa sử dụng 23,76 0,10 108,80 21,00 0,46 0,09 1115,84 4,75 2048,94 8,73 89,03 0,78 686,03 4,77 11,85 0,38 0,00 2,92 0,02 0,05 138,69 12,51 683,50 48,71 14,99 5,59 0,05 3,68 0,21 0,06 18,45 0,08 19,46 0,08 32,79 52,40 66,14 14,49 2,10 11,30 227,70 125,70 0,03 26,03 0,14 0,22 0,28 0,06 0,01 0,05 0,97 0,53 0,00 0,11 56,98 64,47 148,98 29,91 2,10 17,55 195,49 124,16 0,03 14,11 0,24 0,27 0,63 0,13 0,01 0,07 0,83 0,53 0,00 0,06 BẢNG II: Tọa độ điểm mẫu thực địa lớp giao thông, mặt nước đất trống STT Giao thông Mặt nước Đất trống Tọa độ X Tọa độ Y Tọa độ X Tọa độ Y Tọa độ X Tọa độ Y 105°37'32,332"E 21°23'45,262"N 105°31'25,841"E 21°30'33,072"N 105°38'24,583"E 21°23'52,109"N 105°35'12,979"E 21°25'40,405"N 105°31'23,243"E 21°30'47,124"N 105°38'41,926"E 21°23'26,225"N 105°34'19,558"E 21°24'24,257"N 105°31'36,825"E 21°30'20,652"N 105°33'43,060"E 21°29'21,476"N 105°36'5,317"E 21°23'58,327"N 105°37'39,965"E 21°24'17,726"N 105°34'0,988"E 21°28'35,767"N 105°36'41,620"E 21°23'37,542"N 105°38'17,503"E 21°24'14,561"N 105°34'38,690"E 21 27'50,036"N 105°36'32,794"E 21°23'36,284"N 105°39'26,593"E 21°23'5,291"N 105°35'36,874"E 21°26'13,668"N 105°36'47,971"E 21°22'34,009"N 105°36'43,560"E 21°25'31,172"N 105°37'13,192"E 21°23'54,118"N 105°37'45,534"E 21°23'28,012"N 105°31'7,140"E 21°30'30,845"N 105°36'24,795 E 21°24'25,430"N 105°38'28,180"E 21°22'6,010"N 105°31'12,220"E 21°29'52,379"N 105° 35' 4,319" E 21°25'6,593"N 10 105°36'32,780"E 21°23'5,217"N 105°31'29,783"E 21°31'44,612"N 105°31'46,561"E 21°30'43,089"N 11 105°34'41,487"E 21°26'53,060"N 105°31'29,244"E 21°31'44,368"N 105°37'35,907"E 21°24'13,058"N 12 105°34'13,130"E 21°27'5,292"N 105°33'41,724"E 21°29'8,999"N 105°29'58,219"E 21°31'18,969"N 13 105°33'38,022"E 21°26'44,365"N 105°32'56,630"E 21°27'35,226"N 105°30'1,645"E 21°30'5,294"N 14 105°30'33,739"E 21°29'34,109"N 105°36'13,168"E 21°25'40,855"N 105°32'8,477"E 21°28'5,797"N 15 105°37'1,595"E 21°26'15,710"N 105°38'48,043"E 21°23'31,787"N 105°33'33,346"E 21°26'8,482"N 16 105°38'52,341"E 21°22'58,066"N 105°36'0,120"E 21°24'28,692"N 105°33'10,495"E 21°25'53,412"N 17 105°31'34,607"E 21°29'46,241"N 105°35'43,402"E 21°26'1,177"N 105°36'51,478"E 21°23'52,786"N 18 105°38'36,580"E 21°21'33,608"N 105°30'9,510"E 21°30'58,108"N 105°39'30,208"E 21°22'33,386"N 19 105°37'18,844"E 21°23'40,302"N 105°34'31,203"E 21°27'55,541"N 105°38'43,619"E 21°21'28,852"N 20 105°37'11,197"E 21°23'57,369"N 105°34'33,925"E 21°27'46,496"N 105°39'30,804"E 21°22'2,403"N 21 105°34'39,239"E 21°24'52,128"N 105°34'38,374"E 21°25'23,841"N 105°35'48,682"E 21°24'24,368"N 22 105°34'32,786"E 21°25'33,974"N 105°36'54,011"E 21°22'1,554"N 105°38'36,404"E 21°23'35,817"N 23 105°35'57,159"E 21°25'34,194"N 105°36'13,753"E 21°24'20,548"N 105°31'14,630"E 21°31'3,558"N 24 105°35'46,892"E 21°24'10,067"N 105°35'12,044"E 21°25'27,272"N 105°32'29,671"E 21°27'18,711"N 25 105°35'12,294"E 21°24'44,784"N 105°29'55,496"E 21°30'58,258"N 105°33'47,392"E 21°27'29,642"N 26 105°36'0,878"E 21°24'59,586"N 105°29'31,072"E 21°30'42,668"N 105°37'53,492"E 21°22'52,952"N 27 105°31'43,170"E 21°27'54,834"N 105°29'18,215"E 21°30'20,589"N 105°36'6,899"E 21°24'42,636"N 28 105°34'3,050"E 21°26'55,857"N 105°29'52,798"E 21°29'55,272"N 105°33'36,885"E 21°29'46,938"N 29 105°33'51,703"E 21°27'9,320"N 105°31'16,806"E 21°28'26,635"N 105°32'0,187"E 21°31'12,881"N 30 105°32'4,898"E 21°29'44,229"N 105°31'9,689"E 21°28'42,147"N 105°30'38,002"E 21°29'37,614"N 31 105°31'59,118"E 21°30'12,819"N 105°31'28,238"E 21°28'35,377"N 105°38'13,384"E 21°24'14,564"N 32 105°30'59,816"E 21°31'34,243"N 105°32'5,349"E 21°28'38,831"N 105°31'54,488"E 21°31'24,117"N 33 105°34'58,361"E 21°27'35,253"N 105°33'26,550"E 21°29'31,378"N 34 105°33'21,297"E 21°26'2,297"N 105°32'8,483"E 21°30'46,675"N 35 105°33'55,858"E 21°26'2,526"N 105°31'36,787"E 21°30'29,909"N 36 105°37'25,703"E 21°26'43,096"N 105°31'16,619"E 21°30'47,872"N 37 105°37'37,451"E 21°26'49,507"N 105°39'10,912"E 21°21'37,854"N 38 105°37'27,633"E 21°25'16,387"N 105°37'11,573"E 21°24'53,484"N 39 105°36'56,460"E 21°25'34,359"N 105°32'9,079"E 21°26'47,435"N 40 105°34'3,685"E 21°24'38,602"N 105°31'35,484"E 21°27'52,257"N 41 105°33'34,434"E 21°26'20,248"N 105°31'42,482"E 21°29'37,224"N 42 105°38'27,894"E 21°23'8,353"N 105°38'56,716"E 21°24'22,944"N 43 105°38' 22,445"E 21°23'23,219"N 105°31'14,176"E 21°31'37,778"N 44 105°35'45,539"E 21°25'36,497"N 105°30'51,383"E 21°30'48,968"N 105°35'15,087"E 21°24'11,096"N BẢNG III: Tọa độ điểm mẫu thực địa lớp khu dân cư, đất nông nghiệp đất rừng Khu dân cư STT Đất nông nghiệp Đất rừng Tọa độ X Tọa độ Y Tọa độ X Tọa độ Y Tọa độ X Tọa độ Y 105°34'20,482"E 21°25'34,957"N 105°33'44,179"E 21°28'52,270"N 105°35'51,644"E 21°26'41,953"N 105°36'24,451"E 21°23'32,946"N 105°33'18,263"E 21°29'10,950"N 105°32'59,456"E 21°27'52,484"N 105°36'59,767"E 21°22'47,662"N 105°32'1,320"E 21°30'0,018"N 105°33'37,996"E 21°28'30,384"N 105°38'10,374"E 21°22'43,068"N 105°30'32,548"E 21°31'17,832"N 105°32'34,022"E 21°30'43,427"N 105°38'43,923"E 21°27'22,857"N 105°32'29,731"E 21°29'44,228"N 105°30'48,748"E 21°30'19,359"N 105°36'47,855"E 21°24'55,840"N 105°31'45,581"E 21°29'0,265"N 105°32'12,596"E 21°28'59,063"N 105°36'31,086"E 21°25'24,766"N 105°33'21,678"E 21°27'4,993"N 105°34'45,280"E 21°27'16,576"N 105°35'6,857"E 21°26'46,151"N 105°34'16,295"E 21°27'19,073"N 105°32'0,331"E 21°27'50,643"N 105°34'17,468"E 21°27'50,555"N 105°34'11,708"E 21°26'42,018"N 105°30'54,002"E 21°31'32,878"N 10 105°33'20,225"E 21°26'28,133"N 105°36'33,070"E 21°25'1,380"N 105°36'58,268"E 21°24'29,107"N 11 105°31'25,586"E 21°28'58,440"N 105°35'39,188"E 21°24'16,038"N 105°39'5,580"E 21°22'53,103"N 12 105°31'57,094"E 21°29'36,769"N 105°36'13,766"E 21°23'13,571"N 105°40'2,377"E 21°22'15,875"N 13 105°33'0,979"E 21°29'24,310"N 105°38'24,785"E 21°23'16,055"N 105°38'25,468"E 21°24'5,545"N 14 105°32'47,270"E 21°29'45,427"N 105°38'17,250"E 21°22'17,548"N 105°36'52,138"E 21°25'41,417"N 15 105°32'9,272"E 21°30'15,422"N 105°31'0,294"E 21°29'47,633"N 105°37'48,848"E 21°24'37,525"N 16 105°31'7,288"E 21°31'17,296"N 105°38'25,332"E 21°23'11,112"N 105°38'45,114"E 21°21'55,524"N 17 105°29'44,210"E 21°30'14,119"N 105°37'47,233"E 21°23'48,138"N 105°34'2,898"E 21°29'6,830"N 18 105°34'23,837"E 21°26'55,824"N 105°36'26,044"E 21°24'16,426"N 105°37'51,311"E 21°27'2,951"N 19 105°33'4,813"E 21°26'59,269"N 105°35'22,037"E 21°25'36,257"N 105°37'27,275"E 21°26'31,200"N 20 105°37'35,564"E 21°25'34,918"N 105°35'52,399"E 21°23'51,160"N 105°37'35,254"E 21°25'22,325"N 21 105°35'6,220"E 21°27'32,151"N 105°34'43,312"E 21°27'1,588"N 105°30'50,870"E 21°29'47,038"N 22 105°36'10,805"E 21°26'8,060"N 105°37'36,125"E 21°23'24,964"N 105°34'48,334"E 21°25'5,548"N 23 105°34'39,899"E 21°26'30,841"N 105°36'11,408"E 21°24'32,829"N 105°31'37,135"E 21°31'6,926"N 24 105°33'36,551"E 21°27' 2,566"N 105°35'47,253"E 21°25'21,538"N 105°37'3,702"E 21°24'10,155"N 25 105°39'31,761"E 21°21'15,326"N 105°35'0,683"E 21°26'37,862"N 105°37'26,730"E 21°24'22,978"N 26 105°36'7,392"E 21°25'6,727"N 105°34'56,583"E 21°24'51,166"N 105°35'19,503"E 21°26'29,942"N 27 105°32'44,265"E 21°27'5,751"N 105°34'55,006"E 21°25'57,005"N 105°33'52,479"E 21°27'47,693"N 28 105°35'34,287"E 21°26'49,708"N 105°32'37,433"E 21°29'20,993"N 105°35'12,590"E 21°28'20,119"N 29 105°29'52,957"E 21°30'35,649"N 105°38'0,027"E 21°23'5,954"N 105°36'53,807"E 21°25'0,964"N 30 105°32'33,028"E 21°27'32,257"N 105°32'27,121"E 21°30'24,273"N 105°35'46,676"E 21°24'50,880"N 31 105°32'10,090"E 21°28'34,809"N 105°37'3,243"E 21°22'22,104"N 105°35'32,975"E 21°25'11,583"N 32 105°33'42,589"E 21°27'48,293"N 105°38'4,156"E 21°22'0,535"N 105°36'4,011"E 21°24'19,838 N 33 105°32'30,063"E 21°29'5,467"N 105°35'59,839"E 21°25'1,513"N 105° 36' 5,225" E 21°25'32,121"N 34 105°38'17,423"E 21°21'25,565"N 105°36'8,625"E 21°25'55,724"N 105°37'45,976"E 21°26'49,381"N 35 105°37'3,672"E 21°25'52,188"N 105°35'46,643"E 21°26'14,735"N 105°37'22,081"E 21°26'45,654"N 36 105°37'14,896"E 21°24'49,325"N 105°30'33,433"E 21°31'44,164"N 105°36'31,224"E 21°26'0,585"N 37 105°37'10,629"E 21°24'18,360"N 105°32'19,058"E 21°26'54,681"N 105°35'43,692"E 21°26'55,281"N 38 105°34'56,630"E 21°25'1,009"N 105°34'30,907"E 21°26'15,345"N 105°35'19,361"E 21°27'28,083"N 39 105°34'47,086"E 21°25'27,719"N 105°36'26,103"E 21°25'14,698"N 105°34'36,464"E 21°28'20,032"N 40 105°35'28,366"E 21°27'2,238"N 105°35'55,542"E 21°23'32,572"N 105°38'58,414"E 21°23'9,548"N 41 105°34'3,558"E 21°26'1,000"N 105°36'45,362"E 21°22'52,462"N 105°39'20,904"E 21°22'21,478"N 42 105°35'17,104"E 21°26'4,619"N 105°34'13,338"E 21°25'1,051"N 105°35'17,742"E 21°25'28,313"N 43 105°33'9,489"E 21°28'45,550"N 105°35'1,457"E 21°26'23,845"N 105°33'52,911"E 21°26'17,678"N 44 105°34'1,468"E 21°28'53,552"N 105°33'30,314"E 21°26'34,482"N 105°31'37,995"E 21°29'54,585"N 45 105°30'38,474"E 21°29'51,864"N 105°30'38,846"E 21°29'28,522"N 105°30'23,427"E 21°29'56,997"N 46 105°32'7,370"E 21°30'35,475"N 105°29'55,960"E 21°30'11,680"N 105°29'30,527"E 21°30'32,597"N 47 105°30'17,181"E 21°31'23,580"N 105°36'23,950"E 21°25'9,423"N 105°30'31,908"E 21°31'30,210"N 48 105°31'7,337"E 21°29'22,946"N 105°29'43,569"E 21°30'45,795"N 105°31'8,090"E 21°29'0,699"N 49 105°36'39,817"E 21°24'10,795"N 105°30'38,598"E 21°30'2,945"N 105°36'49,421"E 21°28'48,092"N 50 105°35'4,580"E 21°24'36,952"N 105°31'31,555"E 21°28'17,447"N 105°36'7,425"E 21°28'35,158"N 51 105°39'3,060"E 21°22'8,793"N 105°34'12,931"E 21°24'11,451"N 105°35'57,178"E 21°27'11,226"N 52 105°32'12,218"E 21°27'40,278"N 105°38'40,386"E 21°21'41,375"N 105°36'31,175"E 21°25'42,772"N 53 105°30'11,908"E 21°30'20,626"N 105°39'34,817"E 21°23'13,416"N 105°30'43,562"E 21°31'47,501"N 54 105°32'32,603"E 21°27'54,307"N 105°34'28,669"E 21°24'48,134"N 55 105°31'56,128"E 21°27'35,492"N 105°34'27,658"E 21°25'27,396"N 105°34'9,473"E 21°24'53,107"N 56 57 105°31'52,197"E 21°28'14,213"N 58 105°33'14,144"E 21°29'58,636"N 59 105°30'55,892"E 21°31'56,309"N 60 105°30'20,137"E 21°30'56,440"N 61 105°33'18,115"E 21°29'30,270"N 62 105°35'35,353"E 21°26'3,147"N 63 105°39'13,066"E 21°21'20,440"N 64 105°37'58,345"E 21°22'45,241"N 65 105°34'33,919"E 21°24'20,673"N 66 105°31'3,083"E 21°31'0,148"N ... lớp phủ mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc? ?? nhằm thành lập đồ lớp phủ bề mặt huyện Tam Đảo năm 2020 đánh giá tình hình biến động diện tích lớp phủ bề mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc. .. Góp phần xây dựng công cụ quản lý đất đai ứng dụng GIS viễn thám cho khu vực huyện Tam Đảo, tỉnh Vĩnh Phúc 3.2.2 Mục tiêu cụ thể  Thành lập đồ lớp phủ mặt đất huyện Tam Đảo, tỉnh Vĩnh Phúc năm... Bản đồ lớp phủ mặt đất khu vực huyện Tam Đảo, tỉnh Vĩnh Phúc năm 2014 năm 2020 với loại: giao thông, đất rừng, mặt nước, đất nông nghiệp, đất trống khu dân cư  Bản đồ biến động lớp phủ mặt đất