(Luận văn) nghiên cứu, ứng dụng mô hình tích hợp markok ca và gis để dự báo biến động sử dụng đất tại tỉnh phú thọ

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 an n va p ie gh tn to NGUYỄN KHÁNH TOÀN oa nl w NGHIÊN CỨU, ỨNG DỤNG MƠ HÌNH TÍCH HỢP d MARKOV-CA VÀ GIS ĐỂ DỰ BÁO BIẾN ĐỘNG SỬ DỤNG ll u nf va an lu ĐẤT TỈNH PHÚ THỌ oi m z at nh LUẬN VĂN THẠC SĨ QUẢN LÝ ĐẤT ĐAI z m co l gm @ an Lu Hà Nội, 2018 n va ac th si 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 an n va p ie gh tn to NGUYỄN KHÁNH TOÀN w NGHIÊN CỨU, ỨNG DỤNG MƠ HÌNH TÍCH HỢP oa nl MARKOV-CA VÀ GIS ĐỂ DỰ BÁO BIẾN ĐỘNG SỬ DỤNG d ĐẤT TỈNH PHÚ THỌ ll u nf va an lu oi m LUẬN VĂN THẠC SĨ QUẢN LÝ ĐẤT ĐAI z at nh Chuyên ngành: Quản lý đất đai Mã số: 8850103 z @ m co l gm NGƯỜI HƯỚNG DẪN : PGS.TS NGUYỄN AN THỊNH an Lu Hà Nội, 2018 n va ac th si i CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI Cán hướng dẫn chính: PGS.TS Nguyễn An Thịnh Cán chấm phản biện 1: TS Nguyễn Thị Khuy Cán chấm phản biện 2: TS Nguyễn Tiến Cường Luận văn thạc sĩ bảo vệ tại: lu HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ an va TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯỜNG HÀ NỘI n Ngày 16 tháng năm 2018 p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si ii LỜI CAM ĐOAN Tôi xin cam đoan: Những kết nghiên cứu trình bày luận văn hồn tồn trung thực, tơi, khơng vi phạm điều luật sở hữu trí tuệ pháp luật Việt Nam Nếu sai, tơi hồn tồn chịu trách nhiệm trước pháp luật Hà Nội, ngày tháng năm 2018 lu Tác giả Luận văn an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si iii LỜI CẢM ƠN Được trí trường Đại học Tài nguyên Môi trường Hà Nội, Khoa Quản lý đất đai, tiến hành làm luận văn “Nghiên cứu, ứng dụng mơ hình tích hợp Markok- CA GIS để dự báo biến động sử dụng đất tỉnh Phú Thọ” Trong q trình học tập hồn thành luận văn cố gắng nỗ lực thân, nhận giúp đỡ thầy giáo, gia đình bạn bè lu Tôi xin chân thành cảm ơn Ban Giám hiệu, Ban Chủ nhiệm Khoa Quản an n va lý đất đai trường Đại học Tài nguyên Môi trường Hà Nội tồn thể trình học tập rèn luyện trường thời gian qua Đặc biệt để hoàn thành luận văn tốt nghiệp này, nhận quan p ie gh tn to thầy cô giáo dạy dỗ, quan tâm tạo điều kiện thuận lợi cho w tâm giáo viên hướng dẫn PGS.TS Nguyễn An Thịnh Thầy tận tâm oa nl hướng dẫn, giúp đỡ suốt trình thực luận văn d Tơi xin chân thành cảm ơn! lu nf va an Học viên z at nh oi lm ul Nguyễn Khánh Toàn z m co l gm @ an Lu n va ac th si iv MỤC LỤC LỜI CAM ĐOAN LỜI CẢM ƠN MỤC LỤC THÔNG TIN LUẬN VĂN DANH MỤC TỪ VIẾT TẮT DANH MỤC BẢNG lu DANH MỤC HÌNH an n va MỞ ĐẦU 1.1 Cơ sở lý luận sử dụng đất gh tn to CHƯƠNG 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU p ie 1.1.1 Sử dụng đất w 1.1.2 Cơ sở khoa học biến động sử dụng đất lớp phủ oa nl 1.1.3 Sử dụng đất, biến động sử dụng đất giới Việt Nam 12 d 1.2 Mô hình hóa Markov-CA dự báo biến động sử dụng đất 20 lu an 1.2.1 Khái niệm mô hình, mơ hình hóa mơ hình hóa khơng gian 20 nf va 1.2.2 Mơ hình hóa biến động sử dụng đất 21 lm ul 1.2.3 Khái niệm ứng dụng chuỗi Markov 22 z at nh oi 1.2.4 Mạng tự động (Cellular Automata - CA) 24 1.3 Tổng quan cơng trình nghiên cứu .26 1.3.1 Trên giới .26 z gm @ 1.3.2 Tại Việt Nam 31 CHƯƠNG ĐỐI TƯỢNG, NỘI DUNG, PHƯƠNG PHÁP NGHIÊN l co CỨU 34 m 2.1 Đối tượng, phạm vi nghiên cứu .34 an Lu 2.1.1 Đối tượng nghiên cứu .34 n va ac th si v 2.1.2 Phạm vi nghiên cứu 34 2.2 Nội dung nghiên cứu 34 2.3 Phương pháp nghiên cứu 34 2.3.1 Phương pháp phân tích tổng quan tài liệu 34 2.3.2 Phương pháp khảo sát thực địa 35 2.3.4 Phương pháp phân tích, mơ hình hóa khơng gian Markov - CA 36 CHƯƠNG KẾT QUẢ NGHIÊN CỨU 41 3.1 Các yếu tố ảnh hưởng đến sử dụng đất tỉnh Phú Thọ .41 lu 3.1.1 Nguồn lực tự nhiên tài nguyên thiên nhiên 41 an 3.1.2 Nguồn lực kinh tế - xã hội 50 va n 3.1.3 Đánh giá điều kiện tự nhiên kinh tế xã hội 57 gh tn to 3.2 Đánh giá trạng sử dụng đất giai đoạn 2005 – 2015 sở tư liệu ảnh ie viễn thám 60 p 3.2.1 Xây dựng đồ lớp phủ mặt đất tỉnh Phú Thọ năm 2005, 2010, 2015 60 nl w 3.2.2 Đánh giá trạng sử dụng đất giai đoạn 2005 - 2015 65 d oa 3.3 Dự báo biến động sử dụng đất tỉnh Phú Thọ đến năm 2020 - 2025 72 an lu 3.3.1 Quy trình nghiên cứu .72 nf va 3.3.2 Phân cấp thích hợp 74 lm ul 3.3.3 Xây dựng ma trận chuyển dịch dựa vào chuỗi Markov .85 z at nh oi KẾT LUẬN VÀ KIẾN NGHỊ .96 TÀI LIỆU THAM KHẢO .98 PHỤ LỤC .Error! Bookmark not defined z m co l gm @ an Lu n va ac th si vi THÔNG TIN LUẬN VĂN Họ tên học viên: Nguyễn Khánh Tồn Lớp: CH2B.QĐ Khố: 2B (2016 - 2018) Cán hướng dẫn: PGS.TS Nguyễn An thịnh Tên đề tài: Nghiên cứu, ứng dụng mơ hình tích hợp Markov – CA GIS để dự báo biến động sử dụng đất tỉnh Phú Thọ lu Thông tin luận văn: an n va - Tổng quan tài liệu xây dựng sở lý luận ứng dụng mơ hình - Đánh giá nhân tố ảnh hưởng tới biến động sử dụng đất tỉnh Phú Thọ - Đánh giá trạng sử dụng đất tỉnh Phú Thọ giai đoạn 2005 - 2015 p ie gh tn to Markov - CA dự báo xu hướng biến động sử dụng đất w sở tư liệu ảnh viễn thám oa nl - Dự báo biến động sử dụng đất tỉnh Phú Thọ đến năm 2025 sở d mơ hình Markov - CA GIS nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si vii DANH MỤC TỪ VIẾT TẮT lu FAO Tổ chức lương thực nông nghiệp Liên Hợp Quốc LUCC Biến động sử dụng đất lớp phủ TDB Thời điểm dự báo TCT Mốc thời gian cận trình đánh giá TCD Mốc thời gian cận trình đánh giá CA Cellular Automata - Mạng tự động MCE Multi Criteria Evaluation - Đánh giá đa tiêu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si viii DANH MỤC BẢNG Bảng 1.1 Tỉ lệ % đất tự nhiên đất nơng nghiệp tồn giới 13 Bảng 1.2 Biến động diện tích đất nơng nghiệp tồn quốc 17 Bảng 1.3 Biến động diện tích đất 2015 18 Bảng 3.1 Đặc trưng hình thái sơng địa bàn tỉnh Phú Thọ 45 Bảng 3.2 Hiện trạng sủ dụng đất tỉnh Phú Thọ ( tính đến 31/12/2016) 46 lu Bảng 3.3 Hiện trạng rừng theo chức sử dụng 48 an n va Bảng 3.4 Sản lượng khai thác số loại tài nguyên địa bàn tỉnh Phú Bảng 3.5 Diện tích, dân số mật độ dân số huyện năm 2016 51 gh tn to Thọ 50 p ie Bảng 3.6 Đóng góp khu vực kinh tế vào GDP toàn tỉnh qua w năm 55 oa nl Bảng 3.7 Hệ thống giải ảnh lớp phủ mặt đất 61 d Bảng 3.8 Hệ thống mẫu giải đoán ảnh vệ tinh khu vực nghiên cứu 63 an lu Bảng 3.9 Bảng biến động diện tích đất đai tỉnh Phú Thọ 72 nf va Bảng 3.10 Diện tích loại hình sử dụng đất tỉnh Phú Thọ năm 2005, z at nh oi lm ul 2010, 2015, 2020, 2025 (đơn vị: ha) 93 z m co l gm @ an Lu n va ac th si 86 First (earlier) land cover image: nhập liệu đồ TCD Second (later) land cover image: nhập liệu đồ TCT Prefix for output conditional probability images: tên ma trận chuyển dịch đầu Number of time periods between the first and second land cover images: khoảng thời gian thời điểm TCD TCT Number of time periods to project foward from the second image: khoảng thời gian dự báo tương lai lu Kết bước xác định ma trận chuyển dịch lớp an thông tin giai đoạn cần nghiên cứu Ma trận chuyển dịch sở để mơ va n hình dự báo thay đổi đối tượng nghiên cứu tương lai Bao gh tn to gồm: ma trận chuyển dịch xác suất, ma trận chuyển dịch diện tích lớp, ie tập hợp hình ảnh xác suất có điều kiện p Các ma trận chuyển dịch xác suất file văn ghi lại xác suất mà nl w lớp phủ thay đổi Các ma trận chuyển dịch diện tích file văn d oa ghi lại diện tích dự kiến thay đổi từ lớp đối tượng tương lai Tập nf va đổi tương lai an lu hình ảnh xác suất có điều kiện báo cáo xác suất mà loại lớp phủ biến z at nh oi lm ul z m co l gm @ an Lu Hình 3.29 Ma trận chuyển dịch diện tích lớp đối tượng n va ac th si 87 lu an va Hình 3.30 Ma trận chuyển dịch xác suất lớp n p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu Hình 3.31 Tập hình ảnh xác suất có điều kiện n va ac th si 88 3.3.4 Mơ hình hóa biến động sử dụng đất lớp phủ mặt đất dựa vào toán CA_Markov Với ưu điểm đưa yếu tố khơng gian vào tốn thống kê Markov, mơ hình CA_Markov cho phép dự báo biến đổi sử dụng đất khoảng thời gian xác định * Dự báo biến động sử dụng đất, lớp phủ mặt đất tỉnh Phú Thọ đến năm 2015 Dựa nguồn tư liệu đầu vào đồ trạng sử dụng đất năm 2005 năm 2010, mơ hình CA_markov cho phép dự báo biến động lu sử dụng đất tỉnh Phú Thọ đến năm 2015 an Trên Idrisi, chọn Modeling / Environmental-Simulation Models / va n CA_MARKOV để khởi động module Markov_CA p ie gh tn to d oa nl w nf va an lu lm ul z at nh oi Hình 3.32 Mơ hình hóa biến đổi lớp phủ mặt đất dựa vào CA_Markov Basic land cover image: nhập liệu đồ lớp phủ TCT z Markov trasition areas file: ma trận chuyển dịch lớp đối tượng @ gm Transition suitability image collection: tập ảnh phân cấp thích hợp m co tương lai) l Output land cover projection: tệp liệu đầu (bản đồ lớp phủ an Lu n va ac th si 89 Với việc sử dụng đồ phân cấp thích hợp làm ngưỡng giới hạn biến động lớp phủ mặt đất tỉnh Phú Thọ, kết mơ hình hóa thể hình sau: lu an n va p ie gh tn to d oa nl w nf va an lu lm ul z at nh oi Hình 3.33 Kết mơ hình hóa biến đổi sử dụng đất tỉnh Phú Thọ đến năm 2015 z Mục đích cơng đoạn dựa kết mơ hình hóa biến l mơ hình hóa cho giai đoạn gm @ đổi đất tỉnh Phú Thọ đến 2015 để đánh giá mức độ xác q trình m co Sử dụng chức kiểm chứng (Validate) phần mềm Idrisi: GIS an Lu Analysis / Change – Time Series / Validate, so sánh kết mơ hình hóa sử dụng đất đến năm 2015 đồ lớp phủ bề mặt năm 2015: n va ac th si 90 lu an n va Hình 3.34 Kiểm chứng kết mơ hình hóa đồ lớp phủ năm 2015 p ie gh tn to Kết thể hình 3.35 sau: d oa nl w nf va an lu z at nh oi lm ul Hình 3.35 Kết kiểm chứng kết mơ hình hóa đồ lớp phủ z gm @ năm 2015 l Hình cho thấy kết mơ hình hóa đạt tỷ lệ xác cao m co (70%) so với đồ trạng năm 2015 Kết cho phép đề tài mô an Lu hình hóa biến đổi đất tỉnh Phú Thọ đến năm 2020 2025 * Dự báo biến đổi lớp phủ mặt đất tỉnh Phú Thọ đến năm 2020 n va ac th si 91 Ứng dụng mô hình phân tích chuỗi Markov kết hợp với thuật tốn Mạng tự động để dự báo biến đổi lớp phủ mặt đất khu vực nghiên cứu tới năm 2020 2025 cho kết cụ thể sau: lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z co l gm @ m Hình 3.36 Kết dự báo lớp phủ mặt đất tỉnh Phú Thọ năm 2020 an Lu n va ac th si 92 lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ Hình 3.37 Kết dự báo lớp phủ mặt đất tỉnh Phú Thọ năm 2025 an Lu n va ac th si 93 Bảng 3.10 Diện tích loại hình sử dụng đất tỉnh Phú Thọ năm 2005, 2010, 2015, 2020, 2025 (đơn vị: ha) So sánh diện tích năm Diện tích(ha) Rừng kín Rừng t.bình Rừng thưa lu an n va Năm 2010 Năm 2015 Năm 2020 Năm 2025 2005-2020 2005-2025 45632,88 41586,52 43599,57 37600,83 37459,73 -8032,05 -8173,15 34461,09 36508,14 36495,09 27485,19 28578,71 -6975,90 -5882,38 172277,55 170745,62 178744,27 183741,66 185374,98 11464,11 13097,43 10194,66 16336,26 14446,96 15146,96 15577,37 4952,30 5382,71 18827,46 15219,27 14218,91 14327,92 12448,86 -4499,54 -6378,6 9034,02 27728,46 41475,58 46649,53 48352,06 37615,51 39318,04 62584,02 44891,01 24048,27 28085,67 25247,98 -34498,35 -37336,04 p ie gh tn to Cây bụi Mặt nước Đất x.dựng Đất trống Năm 2005 d oa nl w nf va an lu z at nh oi lm ul z gm @ Hình 3.38 Biểu đồ thể diện tích loại hình sử dụng đất tỉnh Phú Thọ l năm 2005, 2010, 2015, 2020, 2025 m co Dựa vào kết dự báo đến năm 2025 ta thấy biến động loại hình sử an Lu dụng đất tỉnh Phú Thọ tương lai sau: so với năm 2005, diện tích rừng thưa tăng 13097,43 ha, diện tích rừng kín rừng trung bình giảm lần n va ac th si 94 lượt 8173,15 5882,38 Diện tích loại bụi tăng 5382,71 , diện tích mặt nước giảm 6378,6 Diện tích đất xây dựng tăng mạnh 39318,04 Diện tích đất trống giảm 37336,04 Nhận xét khả ứng dụng mơ hình phân tích chuỗi Markov Mạng tự động dự báo biến đổi lớp phủ mặt đất: Căn kết dự báo biến đổi lớp phủ mặt đất hai thời điểm xác định, rút số nhận xét thuận lợi khó khăn việc áp dụng mơ hình phân tích chuỗi Markov Mạng tự động sau: lu - Thuận lợi: an + Mơ hình hóa khơng gian nói chung thuật tốn phân tích chuỗi va n Markov Mạng tự động nói riêng phương pháp định lượng có tính + Kết mơ hình hóa cung cấp thơng tin hữu ích xu hướng ie gh tn to logic chặt chẽ, đảm bảo kết mơ hình hóa đạt độ xác cao p biến đổi loại hình lớp phủ mặt đất cho địa phương Trên sở nl w hoạch định sách phù hợp thời gian tới d oa + Mơ hình hóa dự báo biến động sử dụng đất dựa việc tích hợp an lu chuỗi Markov với Cellular Automata GIS tận dụng ưu điểm nf va khắc phục hạn chế mơ hình riêng lẻ Mơ hình hóa tính tốn định lm ul lượng cho thấy diễn biến biến động sử dụng đất z at nh oi không gian tương lai cách logic chặt chẽ - Khó khăn hạn chế mơ hình + Mức độ chi tiết liệu khơng gian đầu vào mơ hình có vai z trò quan trọng Trong điều kiện sở hạ tầng không gian (Spatial Data gm @ Infrastructure) Việt Nam cịn nhiều hạn chế, việc tìm kiếm, thu l m hưởng lớn tới kết mơ hình hóa co thập liệu khơng gian đủ mức độ chi tiết cơng việc khó khăn, ảnh an Lu n va ac th si 95 + Mơ hình hóa khơng gian dựa vào phân tích chuỗi Markov Mạng tự động trình khép kín, gần khơng chịu tác động nhân tố bên ngồi hệ thống Trong đó, biến đổi lớp phủ mặt đất trạng sử dụng đất đối tượng phụ thuộc nhiều vào tác động yếu tố bên ngoài, đặc biệt yếu tố thể chế, sách Do vậy, việc khơng thể tích hợp yếu tố vào mơ hình hạn chế lớn mơ hình + Mức độ chi tiết lớp phủ mặt đất (số lượng đơn vị lớp phủ) lu ảnh hưởng tới kết q trình mơ hình hóa Nếu số lượng đơn vị lớp an phủ mặt đất nhiều (trên 10 đơn vị) làm cho kết mơ hình hóa thiếu va n độ tin cậy p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si 96 KẾT LUẬN VÀ KIẾN NGHỊ KẾT LUẬN Việc ứng dụng GIS kết hợp với chuỗi Markov nghiên cứu trạng sử dụng đất địa bàn tỉnh Phú Thọ giai đoạn 2005 - 2015 dự báo xu hướng biến động sử dụng đất đến năm 2025 Mơ hình hóa khơng gian đối tượng địa lý (cả tự nhiên nhân văn) mạnh công nghệ viễn thám GIS Bằng phần mềm Idrisi, kiểm lu chứng so sánh kết mơ hình hóa đến năm 2015 ảnh phân loại năm 2015 an n va cho thấy kết mơ hình hóa đạt tỷ lệ xác cao (gần 70%) Kết tn to cho phép mơ hình hóa biến đổi lớp phủ mặt đất tỉnh Phú Thọ đến năm 2020 2025 gh p ie Nhược điểm thuật tốn Markov nội suy tuyến tính để dự báo w thay đổi trạng thái pixel theo bước thời gian khác mà chưa oa nl xác định ngưỡng đánh giá (các yếu tố tự nhiên, sách phát triển d yếu tố kinh tế - xã hội) an lu Trong nghiên cứu biến đổi lớp phủ mặt đất, Idrisi Selva 17 phần nf va mềm tích hợp viễn thám GIS cho phép xử lý tư liệu viễn thám (phân loại, lm ul hiệu chỉnh phổ, ) hợp phần GIS (thành lập, biên tập đồ, ), z at nh oi hợp phần mơ hình hóa khơng gian, kết hợp với phân tích chuỗi Markov để đưa ngưỡng xác định phương pháp đánh giá đa tiêu nhằm gia tăng độ xác kết mơ hình hóa Tuy nhiên, nhược điểm z gm @ Idrisi tính tốn mơ hình số lượng đơn vị lớp phủ mặt đất nhiều (trên 10 đơn vị) làm cho kết mơ hình hóa thiếu độ tin cậy yêu cầu l xác m co liệu khơng gian đầu vào mơ hình phải đảm bảo chi tiết an Lu n va ac th si 97 KIẾN NGHỊ Với khả dự báo phát triển loại hình sử dụng đất với trợ giúp hệ thông tin địa lý, kết dự báo dựa vào phân tích chuỗi Markov mạng tự động cho phép dự báo biến động khu vực tương lai dựa sở kết đánh giá biến động qua thời kỳ.Tuy nhiên, q trình triển khai thực khóa luận, với kết khó khăn gặp phải, tác giả xin có số kiến nghị sau: lu - Cần xây dựng sở liệu đủ mạnh có chế chia sẻ an cách phù hợp để bên liên quan tiếp cận sử dụng cách hiệu va n to gh tn - Cần có đủ thời gian điều kiện liên quan để tích hợp yếu tố ie thể chế, sách vào mơ hình để mơ hình có ý nghĩa thực tiễn cao p - Cần có nghiên cứu sâu việc ứng dụng mơ hình vào vững d oa nl w công tác quản lý sử dụng đất tỉnh Phú Thọ hướng tới phát triển bền an lu Việc tích hợp mơ hình chuỗi Markov Cellular Automata cho phép nf va luận văn đưa dự báo sử dụng đất có độ tin cậy cao, nhiên có lm ul thể khơng gian hóa số biến sách phát triển kinh tế - xã hội có liên z at nh oi quan đến sử dụng đất khu vực đưa thêm biến đầu vào mô hình kết dự báo xác z m co l gm @ an Lu n va ac th si 98 TÀI LIỆU THAM KHẢO Tài liệu Tiếng Việt Bộ Tài nguyên Môi trường, (2017), Quyết định số 2311/QĐBTNMT, ngày 28 tháng năm 2017, Quyết định phê duyệt công bố kết thống kê diện tích đất đai năm 2016 Bùi Huy Hiền, Nguyễn Quang Bộ Bùi Quang Xuân (2001) Kết nghiên cứu cải tạo, sử dụng bảo vệ đất dốc sản xuất nông lâm nghiệp, Trong sách: Khoa học Công nghệ Bảo vệ sử dụng bền vững đất lu an dốc, NXB Nông nghiệp n va Cục thống kê tỉnh Phú Thọ, ( 2017) Niên giám thống kê tỉnh Phú Thọ gh tn to năm 2016 Đào Châu Thu Lê Quốc Doanh (2012) Đất dốc vùng đồi núi Việt p ie Nam - Tiềm thách thức, Trong sách: Quản lý bền vững đất nơng nghiệp - Hạn chế thối hóa phịng chống sa mạc hóa, NXB Nơng nghiệp w d nghiệp oa nl Đào Châu Thu Nguyễn Khang (1998) Đánh giá đất, NXB Nông lu nf va an Đinh Thị Bảo Hoa (2007), Nghiên cứu sử dụng hợp lý đất vùng ven đô huyện Thanh Trì, Hà Nội với hỗ trợ cơng nghệ viễn thám hệ lm ul thong tin địa lý z at nh oi Đoàn Đức Lâm Phạm Anh Tuấn, (2010) Ứng dụng GIS thành lập đồ biến động trạng sử dụng đất huyện Sông Mã, tỉnh Sơn La (giai z đoạn 1995-2005) Kỷ yếu hội thảo Ứng dụng GIS toàn quốc 2010 Thành phố gm @ Hồ Chí Minh: Nhà xuất Nơng Nghiệp l Hà Đình Tuấn Lê Quốc Doanh (2007) Nghiên cứu áp dụng m co biện pháp che phủ đất phục vụ phát triển bền vững nông nghiệp vùng cao, Việt Nam an Lu Báo cáo đề tài nghiên cứu khoa học Viện Khoa học Kỹ thuật Nông nghiệp n va ac th si 99 Hoàng Xuân Thành (2006), Thành lập đồ thảm thực vật sở phân tích, xử lý ảnh viễn thám khu vực Tủa Chùa – Lai Châu 10 Nguyễn Kim Lợi, (2005) Ứng dụng chuổi Markov GIS việc đánh giá diễn biến sử dụng đất Trong Kỷ yếu hội thảo Ứng dụng GIS toàn quốc 2011 Thành phố Hồ Chí Minh: Nhà xuất Nơng Nghiệp 11 Nguyễn Thế Đặng, Đào Châu Thu Đặng Văn Minh (2003) Đất đồi núi Việt Nam, NXB Nông nghiệp 12 Phạm Huy Hồng, (2013), Ứng dụng mơ hình Markov – Cellular lu Automata dự báo biến đổi lớp phủ mặt đất huyện Kiến Thuỳ, thành phố Hải an Phòng Đại học Quốc Gia Hà Nội va n 13 Thủ tướng Chính Phủ, (2017), Quyết định số 272/QĐ-TTg ngày 27 gh tn to tháng 02 năm 2007 V/v phê duyệt kết kiểm kê đất đai năm 2005 ie 14 Trần Đức Viên Phạm Chí Thành (1996) Nơng nghiệp đất dốc p - Thách thức tiềm năng, NXB Nông nghiệp nl w 15 Trần Văn Chính, Đỗ Nguyên Hải, Cao Việt Hà, Hoàng Văn Mùa, d oa Nguyễn Hữu Thành Nguyễn Xuân Thành (2006) Giáo trình Thổ nhưỡng an lu học, NXB Nông nghiệp nf va 16 UBND tỉnh Phú Thọ, (2012), Báo cáo thuyết minh tổng hợp quy Tài liệu Tiếng Anh z at nh oi 2015) tỉnh Phú Thọ lm ul hoạch sử dụng đất đến năm 2020 kế hoạch sử dụng đất năm kỳ đầu (2011 - 17 A.A Markov (1906), "Rasprostranenie zakona bol'shih chisel na z obschestva pri Kazanskom universitete l gm @ velichiny, zavisyaschie drug ot druga" Izvestiya Fiziko-matematicheskogo m 2013 co 18 Ellis, E (2010) Land use and land cover change, retrived April an Lu n va ac th si 100 19 FAO (1999) Land use classification for Agri - Enviromental statistics/indicators, Rome, Itatly 20 FAO (2007) Fao statistics, Food and Agriculture Organization of the United Nation, Rome, Italy 21 Meyer, W.B and Turner, B.L (1994) Changes in land use and land cover: A Global Perspective, Cambridge University Press, Cambridge 22 Muller, D and Munroe, D (2007) Issues in spatially explicit statistical land use/cover change (LUCC) models: Examples from western lu Honduras and the Central Highlands of Vietnam, Land use Policy, 24: 521- an 530 va n 23 Turner, B.L and Lambin, E (2001) The causes of land-use and gh tn to land-cover change: moving beyond the myths, Global Environmental Change, p ie 11: 261–269 24 World Bank (2012) Agriculture & Rural Development, retrieved 25 d oa nl w October 2012, from http://data.worldbank.org nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th si