Thử nghiệm phương pháp xây dựng bản đồ hiện trạng lớp phủ thực vật trong lưu vực từ ảnh vệ tinh landsat 8

Quản lý tài nguyên rừng & Môi trường 40 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 3 2015 THỬ NGHIỆM PHƯƠNG PHÁP XÂY DỰNG BẢN ĐỒ HIỆN TRẠNG LỚP PHỦ THỰC VẬT TRONG LƯU VỰC TỪ ẢNH VỆ TINH LANDSAT 8 Ngu[.]

Quản lý tài nguyên rừng & Môi trường THỬ NGHIỆM PHƯƠNG PHÁP XÂY DỰNG BẢN ĐỒ HIỆN TRẠNG LỚP PHỦ THỰC VẬT TRONG LƯU VỰC TỪ ẢNH VỆ TINH LANDSAT Nguyễn Quang Giáp ThS Trường Đại học Lâm nghiệp TĨM TẮT Cơng nghệ viễn thám sử dụng ảnh vệ tinh quang học có độ phân giải trung bình tỏ có nhiều ưu điểm triển vọng điều tra rừng, việc xây dựng đồ trạng cho khu vực có diện tích lớn (như cấp lưu vực) phục vụ theo dõi diễn biến diện tích trạng thái chất lượng rừng Sử dụng ảnh vệ tinh Landsat8, mơ hình số độ cao ASTER (DEM), đồ tài liệu phù trợ, kết hợp với phương pháp điều tra rừng truyền thống nghiên cứu xác định ranh giới lưu vực, xây dựng khóa giải đốn trạng thái rừng đất chưa có rừng lưu vực cách thiết lập quan hệ giá trị kênh phổ theo số thực vật khác biệt chuẩn hóa (NDVI) với trạng thái rừng đất chưa có rừng, thành lập đồ trạng cho lưu vực nghiên cứu với độ xác 81% đề xuất bước kỹ thuật thành lập đồ trạng cấp lưu vực từ ảnh vệ tinh quang học có độ phân giải trung bình mà cụ thể ảnh Landsat Từ khoá: DEM, trạng rừng, Landsat8, lưu vực, NDVI I ĐẶT VẤN ĐỀ nói chung tài nguyên rừng nói riêng Kỹ Bản đồ trạng lớp phủ thực vật cơng cụ quan trọng cho quản lý tài ngun nói chung cho thực sách nơng – lâm nghiệp nói riêng Trong thập kỷ gần đây, công tác quản lý chưa bền vững dẫn đến rừng nước ta có suy giảm rõ rệt số lượng lẫn chất lượng Thảm thực vật rừng bị suy thoái nguyên nhân gây hàng loạt tượng thời tiết cực đoan lũ lụt, hạn hán, mùa… kéo theo tai biến môi trường làm ảnh hưởng lớn đến trình sản xuất sinh hoạt người dân thuật viễn thám với khả quan sát đối Việc điều tra theo dõi diễn biến tài nguyên rừng năm trước chủ yếu phương pháp thủ cơng nên có nhiều nhược điểm độ xác khơng cao, thơng tin khơng cập nhật với tình hình rừng đất Lâm nghiệp ln biến động Hiện nay, với phát triển khoa học công nghệ viễn thám hệ thông tin địa lý, mở hướng nghiên cứu tiếp cận cho việc quản lý hệ thống tài nguyên thiên nhiên 40 tượng độ phân giải khác xạ, phổ, không gian, thời gian cho phép quan sát xác định nhanh chóng lớp phủ thực vật rừng, từ dễ dàng xác định biến động rừng xu hướng biến động, thuận lợi cho việc quản lý cập nhật với độ tin cậy xác định mà phương pháp điều tra xây dựng đồ trạng truyền thống khơng có Một loại ảnh vệ tinh quang học hồn tồn miễn phí, cập nhật thường xuyên ứng dụng nhiều ngành, nhiều lĩnh vực tỏ có nhiều ưu điểm triển vọng áp dụng việc giải đoán trạng lớp phủ thực vật quy mô lưu vực ảnh vệ tinh Landsat Xuất phát từ ý nghĩa thực tiễn việc “Thử nghiệm phương pháp xây dựng đồ trạng lớp phủ thực vật lưu vực từ ảnh vệ tinh Landsat 8” thực để góp phần cung cấp sở cho công tác xây dựng đồ trạng lớp phủ thực vật cấp lưu vực nói riêng cho cơng tác xây dựng đồ trạng rừng nói chung TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3-2015 Quản lý tài nguyên rừng & Môi trường II VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU kiểm kê rừng tỉnh Bắc Cạn 2.1 Đối tượng điều tra khảo sát - Nghiên cứu khoanh vẽ xác định ranh giới đặc trưng lưu vực lựa chọn 2.3 Nội dung nghiên cứu Đối tượng điều tra, khảo sát nghiên cứu toàn rừng đất chưa có rừng 20 tiểu lưu vực nằm lưu vực sông Cầu địa bàn: huyện Bạch Thông, huyện Chợ Đồn Thị xã Bắc Cạn, tỉnh Bắc Cạn - Nghiên cứu xây dựng đồ trạng lớp phủ thực vật lưu vực lựa chọn 2.4 Phương pháp nghiên cứu 2.4.1 Phương pháp khoanh vẽ xác định ranh giới đặc trưng lưu vực lựa chọn 2.2 Vật liệu nghiên cứu Vật liệu nghiên cứu chủ yếu mơ hình số hóa độ cao ASTER GDEM (được tạo Bộ Công nghiệp, Thương mại Kinh tế Nhật Bản phố hợp với NASA Mỹ), ảnh vệ tinh Landsat khu vực tỉnh Bắc Cạn chụp ngày 18 tháng 12 năm 2013, độ phân giải không gian sau hiệu chỉnh tăng cường 15m nắn chỉnh trực giao sai lệch không 30m đồ nền, đồ quy hoạch loại rừng, đồ trạng rừng chu kỳ IV, đồ kết a Phương pháp khoanh vẽ xác định ranh giới lưu vực lựa chọn Ranh giới lưu vực khoanh vẽ tự động từ mơ hình số độ cao (DEM) theo trình tự bước sau (Thông tư số 60/2012/TTBNNPTNT): + Xác định điểm đầu lưu vực: điểm đầu lưu vực nghiên cứu xác định sau Bảng Bảng tổng hợp vị trí điểm đầu 20 lưu vực nghiên cứu TT Điểm đầu Kinh độ Vĩ độ x y BK01 105.897 22.101 592333 2444377 BK02 105.879 22.202 590391 2455533 BK03 105.879 22.201 590408 2455480 BK04 105.854 22.228 587819 2458493 BK05 105.837 22.244 586002 2460176 BK06 105.873 22.249 589799 2460776 BK07 105.878 22.276 590212 2463807 BK08 105.859 22.281 588328 2464313 BK09 105.766 22.167 578739 2451680 10 BK10 105.770 22.168 579212 2451758 11 BK11 105.732 22.126 575339 2447075 12 BK12 105.674 22.105 569299 2444712 Vị trí Thác Bạc - Thơn Thách Giềng - Xã Xuất Hóa Thị Xã Bắc Kạn Thơn Nà Tu - Xã Cẩm Giàng - Huyện Bạch Thông Suối Quăn - Thôn Nà Tu - Xã Cẩm Giàng Huyện Bạch Thông Thôn Pá Yếu - Xã Hạ Vị - Huyện Bạch Thông Suối Nà Giảo - Thôn Pác Chang - Xã Lục Bình Huyện Bạch Thơng Thơn Cóc Pái - Xã Tân Tiến - Huyện Bạch Thông Thôn Chi Quảng B - Xã Phương Linh - Huyện Bạch Thông Suối Vi Hương - Thôn Nà Sang - Xã Vi Hương Huyện Bạch Thông Suối Bản Chịt - Thôn Nà Đán - Xã Đôn Phong Huyện Bạch Thông Suối Bản Đán - Thôn Nà Đán - Xã Đôn Phong Huyện Bạch Thông Suối Bản Mún - Thôn Bản Mún - Xã Dương Phong - Huyện Bạch Thông Thôn Bản Pè - Xã Dương Phong - Huyện Bạch Thông TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3-2015 41 Quản lý tài nguyên rừng & Môi trường TT Điểm đầu Kinh độ Vĩ độ x y 13 BK13 105.659 22.132 567774 2447716 14 BK14 105.660 22.167 567856 2451627 15 BK15 105.662 22.181 568096 2453064 16 BK16 105.667 22.182 568566 2453203 17 BK17 105.643 22.189 566124 2454031 18 BK18 105.629 22.196 564623 2454822 19 BK19 105.637 22.202 565491 2455466 20 BK20 105.635 22.201 565293 2455375 Vị trí Suối Cốc Lùng - Thôn Nà Kệt - Xã Đông Viên Huyện Chợ Đồn Suối Khuổi Gia - Thôn Pác Già - Xã Rã Bản Huyện Chợ Đồn Suối Khuổi Hun - Thôn Khuổi Nhang - Xã Rã Bản - Huyện Chợ Đồn Suối Khuổi Xá - Thôn Nà Thung - Xã Rã Bản Huyện Chợ Đồn Thôn Choong - Xã Phương Viên - Huyện Chợ Đồn Suối Khuổi Đải - Thôn Khuổi Đải - Xã Phương Viên - Huyện Chợ Đồn Suối Nà Đon - Thôn Nà Chúa - Xã Phương Viên Huyện Chợ Đồn Suối Nà Chúa - Thôn Nà Chúa - Xã Phương Viên - Huyện Chợ Đồn + Hiệu chỉnh DEM: DEM tạo thường có vùng trũng không mong muốn không với thực tế Vì vậy, cần phải hiệu chỉnh theo phương pháp “điền trũng” hàm “fill sinks” phần mềm ArcGIS - Chỉ số hình dạng trịn lưu vực: Kc = 0,28P/A0,5 + Xác định hướng dòng chảy: đồ hướng dòng chảy sản phẩm trung gian để tạo đồ tích lũy dịng chảy, thực hàm “flow direction” phần mềm ArcGIS - Độ cao trung bình lưu vực Hbq (m), Độ dốc trung bình lưu vực (Jtb) xác định phần mềm ArcGIS từ mơ hình số độ cao (DM) cho lưu vực + Xây dựng đồ tích lũy dịng chảy lưu vực: đồ tích lũy dịng chảy sở để vẽ lưu vực tự động cách sử dụng hàm “flow accumulation” phần mềm ArcGIS 2.4.2 Phương pháp xây dựng đồ trạng lớp phủ thực vật lưu vực lựa chọn + Xác định ranh giới lưu vực từ điểm đầu lưu vực: tiến hành vẽ lưu vực tự động từ điểm đầu lưu vực xác định trước hàm “watershed” phần mềm ArcGIS - Diện tích lưu vực (A): tính hàm CartesianArea(obj, "hectare") phần mềm Mapinfo quy km2 Sử dụng phần mềm chuyên dụng (chức năng: Multiresolution segmentation phần mềm Ecognition) để khoanh diện tích đồng ảnh nằm ranh giới quy hoạch ba loại rừng lưu vực thành lơ trạng thái có kích thước nhỏ, hay cịn gọi lơ phụ Kết xuất thành tệp đồ gồm ranh giới lô trạng thái rừng với thông tin giá trị trung bình, sai tiêu chuẩn kênh - Chu vi lưu vực (P): tính hàm Perimeter(obj, "km") phần mềm Mapinfo Tính số thực vật NDVI cho lô trạng theo công thức: b Phương pháp xác định đặc trưng lưu vực nghiên cứu 42 Trong đó: Kc số hình dạng trịn lưu vực; P chu vi lưu vực (km); A diện tích lưu vực (km2) Bước 1: Xây dựng khoá giải đoán trạng lớp phủ thực vật lưu vực TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3-2015 Quản lý tài nguyên rừng & Môi trường NDVI  Nir  Re d (2.1) Nir  Re d Trong đó: Nir giá trị phản xạ phổ kênh cận hồng ngoại (Band 5), Red giá trị phả xạ phổ kênh đỏ (Band 4) ảnh Landsat Theo phương pháp điều tra mẫu điển hình, vào danh sách trạng thái rừng lưu vực theo Thông tư số 34/2009/TTBNNPTNT, kết khoanh vi lô trạng thái phân bố kiểu trạng thái rừng (bản đồ trạng chu kỳ IV) đề tài tiến hành điều tra khu vực 100 ô tiêu chuẩn, 210 điểm điều tra phụ để xác định trạng thái trữ lượng rừng, 200 điểm điều tra phụ để kiểm tra độ xác đồ giải đốn Các ô tiêu chuẩn điểm điều tra phụ phân bố tất trạng thái rừng, tuỳ thuộc vào phân bố nhiều hay mức độ tiếp cận đến trạng thái rừng dễ hay khó mà số lượng tiêu chuẩn điểm điều tra phụ trạng thái rừng có khác Đối với tiêu chuẩn: xác định vị trí máy GPS có độ xác cao, đo đường kính ngang ngực, chiều cao vút ngọn, chiều cao cành thước đo chuyên dụng theo quy định điều tra thu thập số liệu ngoại nghiệp Các điểm điều tra phụ: vị trí xác định máy GPS, trạng thái xác định mục trắc hiệu chỉnh lại thông qua trữ lượng lâm phần, xác định tổng tiết diện ngang thước biteclich Tại điểm điều tra phụ xác định tổng tiết diện ngang lâm phần vị trí, vị trí vị trí hướng Đơng, Tây, Nam, Bắc cách vị trí 25m Xác định chiều cao lâm phần cách đo chiều cao có đường kính trung bình lâm phần Các tiêu chuẩn điểm điều tra nằm ô khoanh vi ảnh đại diện cho lơ giá trị phản xạ phổ Chồng xếp vị trí tiêu chuẩn điểm điều tra phụ lên đồ ranh giới lô phụ gán tên trạng thái rừng cho lô phụ chứa ô tiêu chuẩn mặt đất theo tên trạng thái rừng ô tiêu chuẩn, sau tính giá trị số thực vật khác biệt chuẩn hóa (NDVI) cho lơ rừng chứa ô tiêu chuẩn điểm điều tra Thống kê phân ngưỡng giá trị NDVI theo trạng thái rừng khoá xác định tên trạng thái rừng Bước Giải đoán ảnh lưu vực thành lập đồ trạng lớp phủ thực vật lưu vực Sử dụng khoá xác định tên trạng thái rừng, khoá xác định trữ lượng rừng để giải đoán ảnh thành lập đồ trạng rừng khu vực Bước Kiểm tra độ xác kết giải đốn Sử dụng kết đo đếm điểm điều tra phụ (200 điểm) so sánh kết với kết giải đốn ảnh để đánh giá độ xác giải đoán ảnh Cụ thể: - Xác định sai số giải đốn theo trạng thái cơng thức: S%(i)=Nij*100/Ni (2.2) Trong đó: S%(i) sai số giải đốn trạng thái i; Nij tổng số điểm kết giải đoán sai so với thực địa trạng thái i; Ni tổng số điểm kiểm tra trạng thái i - Xác định sai số giải đoán chung cơng thức S%= ∑S%(i)*Ni/N (2.3) Trong đó: S% sai số chung giải đoán; S%(i) sai số giải đoán trạng thái i; Ni tổng số điểm kiểm tra trạng thái i; N tổng số điểm kiểm tra tất trạng thái (N=∑Ni) TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3-2015 43 Quản lý tài nguyên rừng & Môi trường III KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN BK03, BK09 đến BK20 3.1 Kết khoanh vẽ ranh giới lưu vực lựa chọn - Có lưu vực lưu vực nhỏ nằm lưu vực lớn Cụ thể: lưu vực BK04 BK06 nằm lưu vực BK02, lưu vực BK05 nằm lưu vực BK04, hai lưu vực BK07 BK08 nằm lưu vực BK06 Kết khoanh vẽ ranh giới 20 lưu vực ký hiệu từ BK01 đến BK20 cho thấy: - Trong 20 lưu vực nghiên cứu có 14 lưu vực độc lập lưu vực có ký hiệu: BK01, Kết tính tốn đặc trưng hình thái 20 lưu vực nghiên cứu sau: Bảng Đặc trưng hình thái 20 lưu vực nghiên cứu TT 10 11 12 13 14 15 16 17 18 19 20 Ký hiệu lưu vực BK01 BK02 BK03 BK04 BK05 BK06 BK07 BK08 BK09 BK10 BK11 BK12 BK13 BK14 BK15 BK16 BK17 BK18 BK19 BK20 A (km2) 4,57 117,40 15,28 38,76 13,61 65,13 6,11 3,29 37,96 12,73 13,41 4,86 6,15 2,92 13,26 3,33 7,66 5,53 7,91 4,23 Lưu vực có diện tích lớn BK02 (117,4 km2), nhỏ BK14 (2,92 km2), diện tích trung bình lưu vực nghiên cứu khoảng 19,2 km2 Lưu vực có chu vi lớn BK02 (55,5 km), nhỏ BK14 (7,30 km), chu vi trung bình lưu vực nghiên cứu khoảng 18,3 km Lưu vực có độ cao trung bình lớn 44 P (km) 11,0 55,5 23,2 29,4 17,1 40,3 12,5 8,8 27,5 19,3 18,1 10,2 12,2 7,3 17,4 7,6 14,1 10,7 14,2 10,3 Kc H (m) 1,4 1,4 1,7 1,3 1,3 1,4 1,4 1,4 1,2 1,5 1,4 1,3 1,4 1,2 1,3 1,2 1,4 1,3 1,4 1,4 506 387 325 465 590 378 362 294 583 483 465 372 386 430 565 583 521 726 768 818 J 24,9 18,1 19,2 18,7 21,3 18,6 19,7 15,5 25,1 23,9 23,1 19,2 16,2 21,4 23,8 28,2 19,6 24,5 22 23,3 BK20 (cao 818m), nhỏ BK08 (cao 294m), trung bình có độ cao 500m Lưu vực có độ dốc trung bình lớn BK16 (dốc 28,20), nhỏ BK08 (dốc 15,50), trung bình có độ dốc 21,30 Về hình dạng: lưu vực có hình dạng trịn BK16, lưu vực có hình dạng dài BK03 TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3-2015 Quản lý tài ài nguyên rrừng & Mơi trường Hình Ranh giới 20 lưu vực nghiên cứu Kết nội dung nghiên ên cứu c này, khoanh vẽ ranh giới xác định đ đặc trưng 20 lưu vực ực phục vụ cho nội dung nghiên cứu tiếp ếp theo 3.2 Kết xây dựng ựng đồ trạng lớp phủ thực vật lưu vực ực lựa chọn 3.2.1 Xây dựng ựng khoá giải đốn ảnh Khố giải ải đốn nhóm đối tượng t ảnh tập hợp đặc trưng vềề phản xạ phổ v thông tin bổ sung dùng ùng để đ xác định tên gọi chỉỉ thuộc tính nhóm đối tượng Kết nghiên cứu cho th thấy lưu vực nghiên cứu có thểể xác định đđược trạng thái rừng chủ yếu: (1) Rừng ừng gỗ tự nhi nhiên núi đất rộng thường xanh (LRTX LRTX) trung bình, (2) Rừng gỗ tự nhiên ên núi đđất LRTX nghèo, (3) Rừng gỗ tự nhiên núi đất ất LRTX phục hồi, ((4) Rừng hỗn giao gỗỗ tre nứa tự nhiên núi đất, (5) Rừng gỗ trồng núi đấtt trạng thái khơng có rừng là: đất ất trống núi đất đất; đất nông nghiệp mặt nước Kết ết tính tốn giá trị NDVI cho trạng ạng thái loại đất, loại rừng lưu vực thể bảng ảng 002 Bảng Thống kêê giá trị tr NDVI trung bình cho loại đất, loại rừng lưu vực TT Lo đất loại rừng Loại NDVI s n Mặt nước -0,067 0,104 15 Đất trống 0,169 0,026 13 Đất ất trồng nông nghiệp 0,260 0,033 53 Rừng ừng gỗ trồng núi đất 0,358 0,020 52 TẠP ẠP CHÍ KHOA HỌC VÀ V CÔNG NGHỆ Ệ LÂM NGHIỆP SỐ 3-2015 45 Quản lý tài nguyên rừng & Môi trường TT Loại đất loại rừng NDVI s n Rừng gỗ tự nhiên núi đất rộng thường xanh nghèo 0,398 0,004 21 Rừng gỗ tự nhiên núi đất rộng thường xanh trung bình 0,418 0,010 40 Rừng hỗn giao gỗ + tre nứa 0,445 0,006 36 Rừng gỗ tự nhiên núi đất rộng thường xanh phục hồi 0,469 0,011 25 Như vậy, giá trị NDVI thấp thuộc mặt nước tăng dần qua đối tượng đất trống, đất rồng nông nghiệp, rừng trồng gỗ, rừng thường xanh nghèo, rừng thường xanh trung bình, rừng hỗn giao rừng thường xanh phục hồi Rừng thường xanh nghèo, thường xanh trung bình có giá trị NDVI cao so với rừng trồng gỗ, nguyên nhân trạng thải rừng có số lượng bụi, thảm tươi dây leo phát triển mạnh so với rừng trồng gỗ Mặt nước hấp thụ mạnh ánh sáng vùng Để xác định khả tách đối tượng thành lớp riêng biệt giựa vào số thực vật NDVI, đề tài so sánh trung bình tổng thể cặp đối tượng với theo công thức hồng ngoại nên giá trị phổ thấp Đất trống chủ yếu cỏ xen lẫn bụi nhỏ tái sinh rải rác với khả hấp thụ ánh sáng đỏ phản xạ ánh sáng hồng ngoại tăng dần = làm giá trị NDVI đối tượng cao | − + mặt nước lại thấp so với đất trồng nơng nghiệp có mật độ dầy đặc Rừng trồng khu vực nghiên cứu chủ yếu vào giai đoạn khép tán, mật độ cao, chế độ chăm sóc tốt nên lượng bụi, thảm tươi tán rừng nên giá trị NDVI tương đối cao nhỏ so với rừng tự nhiên | Trong đó: X1, S1, n1 trung bình, độ lệch chuẩn số điểm điều tra lớp đối tượng thứ nhất, X2, S2, n2 trung bình, độ lệch chuẩn số điểm điều tra lớp đối tượng thứ hai Nếu U>1,96 tách lớp đối tượng khỏi ngược lại Bảng So sánh giá trị trung bình NDVI lớp đối tượng TT Loại đất, loại rừng Mặt nước Đất trống Đất trồng nông nghiệp Rừng gỗ trồng núi đất Rừng gỗ tự nhiên núi đất rộng thường xanh nghèo Rừng gỗ tự nhiên núi đất rộng thường xanh trung bình Rừng hỗn giao gỗ + tre nứa Rừng gỗ tự nhiên núi đất rộng thường xanh phục hồi Kết so sánh bảng cho thấy sử dụng số NDVI để phân tách trạng thái rừng đất chưa có rừng khu vực nghiên cứu thành phần riêng biệt Từ kết tính 46 TB -0,067 0,169 0,260 0,358 0,398 0,418 0,445 0,469 NDVI s 0,104 0,026 0,033 0,020 0,004 0,010 0,006 0,011 n 15 13 53 52 21 40 36 25 U 8,5 10,6 18,0 13,7 11,7 14,0 9,8 giá trị NDVI, đề tài xác định ngưỡng giá trị NDVI cho loại đất, loại rừng khố xác định trạng thái rừng giải đoán ảnh vệ tinh lưu vực (bảng 4) TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 3-2015 ... mềm ArcGIS 2.4.2 Phương pháp xây dựng đồ trạng lớp phủ thực vật lưu vực lựa chọn + Xác định ranh giới lưu vực từ điểm đầu lưu vực: tiến hành vẽ lưu vực tự động từ điểm đầu lưu vực xác định trước... khoanh vẽ ranh giới lưu vực lựa chọn - Có lưu vực lưu vực nhỏ nằm lưu vực lớn Cụ thể: lưu vực BK04 BK06 nằm lưu vực BK02, lưu vực BK05 nằm lưu vực BK04, hai lưu vực BK07 BK 08 nằm lưu vực BK06 Kết khoanh... 105 .87 9 22.201 5904 08 2455 480 BK04 105 .85 4 22.2 28 587 819 24 584 93 BK05 105 .83 7 22.244 586 002 2460176 BK06 105 .87 3 22.249 589 799 2460776 BK07 105 .87 8 22.276 590212 246 380 7 BK 08 105 .85 9 22. 281 588 328