Xem mẫu
- Laboratory Exercise No. 1 GIS & Makov Chain Model
PHẦN CHUYÊN ĐỀ
ỨNG DỤNG GIS TRONG QUẢN LÝ
TÀI NGUYÊN THIÊN NHIÊN
Designed by Dr. Nguyen Kim Loi, July 2005 106
- Laboratory Exercise No. 1 GIS & Makov Chain Model
Chuyên đề 1
ỨNG DỤNG GIS & CHUỔI MAKOV TRONG VIỆC
ĐÁNH GIÁ DIỄN BIẾN SỬ DỤNG ĐẤT
Mục tiêu:
1) Ứng dụng GIS trong việc xây dựng bản đồ sử dụng đất, tài nguyên rừng;
2) Ứng dụng chuổi Makov & GIS trong việc xác định Ma trận chuyển đ ổi sử
dụng đất, tài nguyên rừng;
3) Dự báo thay đổi sử dụng đất và tài nguyên rừng thông qua các mô hình tóan
học.
Tài liệu tham khảo
(1) Nguyen Kim Loi, and N.Tangtham. . 2001. Effects of Land Cover and Large
Reservoir Operation on Water Balance of the Chao Pharaya River Basin, Bangkok,
Thailand. Kasetsart J. Nat. Sci. 32. p.234-244.
(2) Nguyen Kim Loi. 2002. Effect of Land Use/ Land Cover Changes and Practices on
Sediment Contribution to the Tri An Reservoir of Dong Nai Watershed, Vietnam.
M.Sc. Thesis, Graduate School, Kasetsart University, Bangkok, Thailand. 113p.
(3) Nguyễn Kim Lợi .2002. Tiếp cận mô hình hóa trong nghiên cứu thay đổi sư dụng
đất tại lưu vực sông Đồng Nai. Tập san Kỹ thuật Nông Lâm nghiệp, số 1/2002.
Đại học Nông Lâm TPHCM. p.34-40.
(4) Nguyễn Kim Lợi .2004. Anh hưởng của sự thay đổi các kiểu sử dụng đất đến
bồi lắng lòng hồ Trị An của lưu vực sông Đồng Nai.Tạp chí Khoa học đất, Số
đặc biệt 20/2004, Hà Nội.
(5) Nguyen Kim Loi. 2002. Integration of Linear Programming and GIS Land
Allocation for Quality of Life: A Case Study in Dong Nai Watershed, Vietnam.
Seminar, Kasetsart J. (Nat.Sci.) 35. Bangkok, Thailand. P.511-519.
(6) Nguyen Kim Loi., and N. Tangtham. 2004. Decision support system for
sustainable watershed management in Dong Nai watershed – Vietnam: Conceptual
framework a proposed research techniques. Paper presented in Forest and Water
in Warm Humid Asia, IUFRO Workshop, July 10-12, 2004 Kota Kinabalu,
Malaysia.
(7) Nguyen Kim Loi. and N. Tangtham. 2004. Decision support system for sustainable
watershed management in Dong Nai watershed – Vietnam: Applying Linear
Programming Technique for Land Allocation. Paper presented in International
Environmental Modelling and Software Society iEMSs 2004 International
Conference . Complexity and Integrated Resources Management Session. 14-17
June 2004 University of Osnabrck, Germany.
Designed by Dr. Nguyen Kim Loi, July 2005 107
- Laboratory Exercise No. 1 GIS & Makov Chain Model
(8) Nguyen Kim Loi., and N. Tangtham. 2005. Decision support system for sustainable
watershed management in Dong Nai watershed – Vietnam. Paper presented in
International Seminar on “Synergistic Approach to Appropriate Forestry
Technology for Sustaining Rainforest Ecosystem”, March 7 - 9, 2005, Bintulu
Kinabalu, Malaysia.
Vật liệu cần chuẩn bị:
1) Bản đồ sử dụng đất
2) Máy tính với các phần mềm hỗ trợ: Arc View, Arc GIS,…
Các kiến thức căn bản
1. Chuổi Makov
1. 1 Xác định sự thay đổi các kiểu sử dụng đất dựa trên cơ sở của mô hình
Markov Chain.
Mô hình Markov Chain đã được ứng dụng để xác định khả năng thay đổi các kiểu
sử dụng đất dựa trên sự tiến triển các kiểu sử dụng đất và các nhân tố ảnh hưởng
đến sự thay đổi. Tổng quát hóa của mô hình được minh họa như sau:
Designed by Dr. Nguyen Kim Loi, July 2005 108
- Laboratory Exercise No. 1 GIS & Makov Chain Model
Các kiểu sự dụng đất ở thời điểm to Các kiểu sự dụng đất ở thời
điểm t1
(1)Lâm γ 11
(1)Lâm nghiệp
nghiệp γ 12
γ 51
γ 13 (2)Nông nghiệp
(2)Nông nghiệp γ 52
(3)Đô thị
(3)Đô thị γ 14
(4)Đất trống
γ 53
(4)Đất trống γ 54
(5)Đất chuyên
γ 55
(5)Đất chuyên dùng
dùng
Với γ ij : là xác suất thay đổi được xác định từ việc “Overlay” bản đồ
sử dụng đất tại 2 thời điểm khác nhau. Để dự đoán phân bố các kiểu sử dụng đ ất
khác nhau vào các thời điểm tiếp theo có thể ứng mô hình Markov Chain như sau:
Ma trận về Tỉ lệ các kiểu
Tỉ lệ các kiểu sử dụng xác suất của sử dụng đất
đất ở thời điểm thứ * sự thay đổi = ở thời điểm
nhất các kiểu sử thứ
dụng đất hai
Đây có thể được viết lại dưới dạng tổng quát hóa của ma trận như sau:
γ 11, γ 12, γ 13, . . . γ 15
γ γ 22, γ 23, . . . γ 25 = [V , V , . . . , V ]
[V1, V2, . . . , V5] * 21, 1 2 5
2
1
.
.
1.2. Ứng dụng mô hình toán học để dự báo sự thay đổi các kiểu sử dụng đ ất
.
trong tương lai γ 51, γ 52, γ 53, . . . γ 55
Designed by Dr. Nguyen Kim Loi, July 2005 109
- Laboratory Exercise No. 1 GIS & Makov Chain Model
Dự báo về sự thay đổi các kiểu sử sự đất theo thời gian theo phương trình toán
học sau:
dA/dt = A* r
dA/A = r * dt
At = Ao * (1+r)t (1)
r : Tỉ lệ thay đổi của các kiểu sử dụng đất trong khoảng thời gian thu
thập số liệu.
Ao: Diện tích của kiểu sử dụng đất tại thời điểm thứ nhất
At: Diện tích của kiểu sử dụng đất tại thời điểm t.
1.3. Mô hình hóa về sự thay đổi các kiểu sử dụng đất tại lưu vực sông Đồng Nai
Như đã đề cập ở trên, mục đích của nghiên cứu nhằm phát triển mô hình về sự
thay đổi các kiểu sử dụng đất (LUCC) tại lưu vực sông Đồng Nai nhằm ứng dụng
mô hình này để dự báo hướng LUCC trong tương lai.
Sự thay đổi các kiểu sử dụng đất được xem như là một hàm số do bởi những
hoạt động của con người trong lưu vực. Hàm số đó được xác định bởi sự thay đổi
các kiểu sử dụng đất và có thể được mô hình hóa bởi một phương trình tóan h ọc
như sau:
LUCC = f(sự biến động về dân số, sự tiến bộ của khoa học kỹ thuật, kinh tế xã hội
trong lưu vực, sự nhận thức của người dân, chính sách xã hội) (2)
Thuật ngữ “Patch” (P) được sử dụng nhằm đại diện cho sự xuất hiện của các
kiểu sử dụng đất, được mô tả như sau:
P1 = Đất lâm nghiệp
P2 = Đất nông nghiệp
P3 = Khu vực dân cư
P4 = Đất trống
P5 = Đất chuyên dùng.
Sự thay đổi về sử dụng đất của mỗi “Patch” theo thời gian dựa trên cơ sở sự
tác động giữa sự biến động về dân số, sự tiến bộ của khoa học kỹ thuật, kinh tế xã
hội trong lưu vực, sự nhận thức của người dân trong lưu vực và chính sách xã hội.
Trong nghiên cứu này, tại thời điểm t1, diện tích của mỗi một “Patch” là một
hàm số với biến là diện tích của mỗi “Pacth” với hệ số ci, có thể được viết như sau:
(AP1) t1 = c1 AP1(to) (3)
(AP2) t1 = c2 AP2(to) (4)
(AP3) t1 = c3 AP3(to) (5)
( AP4) t1 = c4 AP4(to) (6)
(AP5) t1 = c5 AP5(to) (7)
Trong đó:
c1, c2, . . . c5 : hệ số LUCC;
t = thời gian;
AP1, … AP5: diện tích của P1, . . . P5
Designed by Dr. Nguyen Kim Loi, July 2005 110
- Laboratory Exercise No. 1 GIS & Makov Chain Model
Phương trình (3) đến (7), có thể được tổng quát hóa như sau:
(APn) (t+1) = cn APn (t) (8)
Tại thời điểm các năm 1990, 1995 và 2000, diện tích của lưu vực có thể
được viết như sau:
1990: A (t1) = AP1 (t1) + AP2 (t1) + AP3 (t1) + AP4 (t1) + AP5 (t1) (9)
1995: A (t2) = AP1 (t2) + AP2 (t2) + AP3 (t2) + AP4 (t2) + AP5 (t2) (10)
2000: A (t3) = AP1 (t3) + AP2 (t3) + AP3 (t3) + AP4 (t3) + AP5 (t3) (11)
A (t1) = A (t2) = A (t3) tổng diện tích của khu vực nghiên cứu.
Trong chuyên đề này, khoảng giữa thời gian to và t1, sự thay đổi giữa các
“Patch” là một ma trận được trình bày trong Bảng 1.
Bảng 1. Ma trận về thay đổi các kiểu sử dụng đất trong khoảng thời gian to đến t1
t
o P1 P2 P3 P4 P5
t1
P1 γ 11 γ 12 γ 13 γ 14 γ 15
P2 γ 21 γ 22 γ 23 γ 24 γ 25
P3 γ 31 γ 32 γ 33 γ 34 γ 35
P4 γ 41 γ 42 γ 43 γ 44 γ 45
P5 γ 51 γ 52 γ 53 γ 54 γ 55
Sự thay đổi của “Patch” P1 khoảng giữa thời gian to đến t1, có thể được viết như sau:
(AP1) t1 = c1AP1(to)
= AP1( to) – γ 12AP1(to) – γ 13AP1 (to) – γ 14AP1 (to) – γ 15AP1(to) + γ 21AP2 (to) +
γ 31AP3 (to) + γ 41AP4 (to) + γ 51AP5 (to) (12)
Tương tự như trên,sự thay đổi của “Patch” P2, P3, P4, P5, khoảng giữa thời gian
to- t1, có thể được viết như sau:
D’;/
(AP2) t1 = c2 AP2(to)
= AP2( to) ) – γ 21AP2(to) – γ 23AP2 (to) – γ 24AP2 (to) – γ 25AP2(to) + γ 12AP1 (to) + γ 32AP3
(to) + γ 42AP4 (to) + γ 52AP5 (to) (13)
(AP3) t1 = c3AP3(to)
= AP3( to) – γ 31AP3(to) – γ 32AP3 (to) – γ 34AP3 (to) – γ 35AP3(to) + γ 13AP1 (to) +
γ 23AP2 (to) + γ 43AP4 (to) + γ 53AP5 (to) (14)
(AP4) t1 = c4AP4(to)
= AP4( to) – γ 41AP4(to) – γ 42AP4 (to) – γ 43AP4 (to) – γ 45AP4(to) + γ 14AP1 (to) +
γ 24AP2 (to) + γ 34AP3 (to) + γ 54AP5 (to) (15)
Designed by Dr. Nguyen Kim Loi, July 2005 111
- Laboratory Exercise No. 1 GIS & Makov Chain Model
(AP5) t1 = c5AP5(to)
= AP5( to) – γ 51AP5(to) – γ 52AP5 (to) – γ 53AP5 (to) – γ 54AP5(to) + γ 15AP1 (to) +
γ 25AP2 (to)+ γ 35AP3 (to) + γ 45AP4 (to) (16)
Với
(APi)t1 = diện tích của “Patch” Pi tại thời điểm t1.
(AP1)to = diện tích của “Patch” Pi tại thời điểm to.
ci = hệ số thay đổi của “Patch” Pi
γ ij = hệ số thay đổi các kiểu sử dụng đất, chỉ ra khả năng thay đổi từ
“Patch” Pi đến “Patch” Pj.
Phương trình (12) “dấu (+)” chỉ sự chuyển đổi từ Patch “P 2, P3, P4, P5” đến
Patch P1, và “dấu (-)” chỉ sự chuyển đổi từ Patch P 1 đến Patch “P2, P3, P4, P5”. Những
phương trình khác tương tự “dấu (+)” và “dấu (-) giải thích sự chuyển đổi theo γ ij và
APi.
Nội dung thực hiện:
- Cơ sở dữ liệu GIS về bản đồ sử dụng đất;
- Từ Bản đồ sử dụng đất của các giai đọan khác nhau, Hãy xác định Ma trận
chuyển đổi sử dụng đất.
- Dựa vào ma trận chuyển đổi, đánh giá diễn biến sử dụng đất, tài nguyên rừng
(Đánh giá kể cả phần chính sách…)
- Trên cơ sở ma trận chuyển đổi, hãy dự báo thay đổi sử dụng đất trong tương
lai.
- Nhưng thuận lợi & khó khăn của phương pháp này.
- Các kiến nghị và đề xuất
Designed by Dr. Nguyen Kim Loi, July 2005 112
nguon tai.lieu . vn
