Kỹ thuật làm trơn đường lô sau giải đoán ảnh vệ tinh

Công nghệ viễn thám được chứng minh có nhiều ưu điểm và triển vọng trong điều tra rừng, nhất là trong việc xây dựng bản đồ hiện trạng rừng và theo dõi diễn biến diện tích các trạng thái và chất lượng rừng ở nước ta. Một trong những bước công việc mất nhiều thời gian sau giải đoán là làm trơn đường lô để biên tập bản đồ hiện trạng rừng. Bài báo sử dụng một số phần mềm chuyên dụng, lớp bản đồ sau giải đoán từ ảnh vệ tinh SPOT6 tại xã Quảng Thành, thị xã Gia Nghĩa, tỉnh Đắk Nông, các bản đồ và tài liệu phù trợ... để làm trơn đường bao lô bằng 3 phương pháp: 1) Bằng công cụ có sẵn trên phần mềm eCognition; 2) Bằng công cụ có sẵn trên phần mềm ArcGIS; 3) Bằng tổ hợp các công cụ trên phần mềm ArcGIS. Kết quả cho thấy sai số về diện tích xác định theo phương pháp bình phương nhỏ nhất của lớp bản đồ làm trơn theo phương pháp 3 so với lớp bản đồ sau giải đoán ở cấp lô rừng là 9,4 và ở cấp trạng thái rừng là 0,03. Các giá trị này nhỏ hơn so với sai số khi làm trơn bằng phương pháp 1 (Sai số cấp lô: 665,8 ; cấp trạng thái: 0,11) hoặc phương pháp 2 (Sai số cấp lô: 12,9; cấp trạng thái: 0,06). Từ kết quả nghiên cứu đã xác định các bước kỹ thuật làm trơn đường lô: 1) Chuyển lớp bản đồ sau giải đoán từ dạng vùng thành dạng đường; 2) Làm trơn lớp bản đồ dạng đường; 3) Chuyển lớp đường sau làm trơn thành lớp vùng và cập nhật dữ liệu để sử dụng.

Qu n lý tài nguyên r ng & Môi tr TH ng NGHI M PH NG PHÁP LÀM TR N SAU GI I ỐN T NH V TINH NG LƠ Ph m V n Du n, V Th Thìn, 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 đ c ch ng minh có nhi u u m tri n v ng u tra r ng, nh t vi c xây d ng b n đ hi n tr ng r ng theo dõi di n bi n di n tích tr ng thái ch t l ng r ng n c ta M t nh ng b c công vi c m t nhi u th i gian sau gi i đoán làm tr n đ ng lô đ biên t p b n đ hi n tr ng r ng Tác gi s d ng m t s ph n m m chuyên d ng, l p b n đ sau gi i đoán t nh v tinh SPOT6 t i xã Qu ng Thành, th xã Gia Ngh a, t nh k Nông, b n đ tài li u phù tr đ làm tr n đ ng bao lô b ng ph ng pháp: 1) B ng cơng c có s n ph n m m eCognition; 2) B ng cơng c có s n ph n m m ArcGIS; 3) B ng t h p công c ph n m m ArcGIS K t qu cho th y sai s v di n tích xác đ nh theo ph ng pháp bình ph ng nh nh t c a l p b n đ làm tr n theo ph ng pháp so v i l p b n đ sau gi i đoán c p lô r ng 9,4 c p tr ng thái r ng 0,03 Các giá tr nh h n so v i sai s làm tr n b ng ph ng pháp (Sai s c p lô: 665,8 ; c p tr ng thái: 0,11) ho c ph ng pháp (Sai s c p lô: 12,9; c p tr ng thái: 0,06) T k t qu nghiên c u xác đ nh b c k thu t làm tr n đ ng lô: 1) Chuy n l p b n đ sau gi i đoán t d ng vùng thành d ng đ ng; 2) Làm tr n l p b n đ d ng đ ng; 3) Chuy n l p đ ng sau làm tr n thành l p vùng c p nh t d li u đ s d ng T khoá: Hi n tr ng r ng, ph n m m ArcGIS, sau gi i đoán, SPOT6, tr n đ I TV N B n đ hi n tr ng r ng m t nh ng công c quan tr ng cho qu n lý tài nguyên r ng M t nh ng ph ng pháp tri n v ng nh t hi n đ xây d ng b n đ hi n tr ng r ng đ m b o khoa h c, chi phí h p lý v i đ xác c n thi t đ c c p nh t th ng xuyên ng d ng công ngh vi n thám h th ng thơng tin đ a lý đ gi i đốn tr ng thái, tr l ng c a khu r ng Công tác xây d ng b n đ hi n tr ng r ng t nh v tinh s li u th c đ a g m b c chính: B c 1: gi i đốn nh; B c 2: biên t p b n đ hi n tr ng Tr c đây, vi c gi i đoán nh ch y u khoanh v tr c ti p n n nh b i k thu t viên đốn đ c thơng qua ph n m m chuyên d ng Cách gi i đoán ph thu c vào kinh nghi m c a ng i gi i đốn, th ng có s khơng đ ng nh t v k t qu gi a nh ng cán b gi i đoán khác nhau, m t nhi u th i gian, nh ng đ ng lô l p b n đ k t qu gi i đốn tr n có th s d ng 38 ng lô đ biên t p b n đ hi n tr ng Hi n nay, vi c gi i đoán ch y u đ c th c hi n t đ ng b ng ph n m m chuyên d ng V i ph ng pháp v a ti t ki m đ c th i gian l i b nh h ng c a y u t kinh nghi m c a ng i gi i đốn nh ng đ ng lơ b n đ k t qu gi i đoán khơng tr n mà dích d c theo đ phân gi i c a nh v tinh, nên sau gi i đốn ph i m t nhi u cơng đ làm tr n đ ng bao lô M t khác, ph n m m gi i đoán nh chuyên d ng nh : Erdas, Ecognition, ArcGIS đ u có s n cơng c đ làm tr n đ ng lơ sau gi i đốn Tuy nhiên, theo đánh giá cơng c cho hình nh đ ng lô làm tr n th ng không phù h p v i n n nh v tinh theo mong mu n Vì v y,“Th nghi m ph ng pháp làm tr n đ ng lô sau gi i đoán t nh v tinh” đ c th c hi n nh m đ xu t b c k thu t làm tr n đ ng bao lô r ng đ t khơng có r ng sau gi i đoán, nâng cao hi u su t c a công tác xây d ng b n đ hi n tr ng t k t qu gi i đoán nh v tinh T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 Qu n lý tài nguyên r ng & Môi tr II PH NG PHÁP NGHIÊN C U 2.1 M c tiêu nghiên c u Xác đ nh b c k thu t làm tr n đ ng bao lô r ng đ t khơng có r ng nh m nâng cao hi u qu công tác xây d ng b n đ hi n tr ng r ng t k t qu gi i đoán nh v tinh 2.2 N i dung nghiên c u - Nghiên c u hi u qu làm tr n đ ng lơ sau gi i đốn m t s ph n m m thông d ng - Nghiên c u k thu t làm tr n đ ng lô b ng t h p công c ph n m m ArcGIS 2.3 V t li u nghiên c u th c hi n n i dung nghiên c u trên, v t li u nghiên c u nh sau: - L p b n đ sau gi i đoán t nh SPOT6 t i xã Qu ng Thành, th xã Gia Ngh a, t nh k Nơng - B n đ rà sốt quy ho ch ba lo i r ng xã Qu ng Thành, Th xã Gia Ngh a, t nh k Nông n m 2013 - Các ph n m m: eCognition, ArcGIS, Mapinfo, Excel 2.4 Ph ng pháp nghiên c u 2.4.1 Ph ng pháp nghiên c u hi u qu làm tr n đ ng lô sau gi i đốn m t s ph n m m thơng d ng Trên ph n m m gi i đoán nh c ng nh ph n m m thu c h th ng thông tin đ a lý đ u có ch c n ng làm tr n đ ng bao lô Trong n i dung này, báo nghiên c u hi u qu làm tr n đ ng bao lơ sau gi i đốn t i xã Qu ng Thành, th xã Gia Ngh a, t nh k Nông ph n m m hay đ c s d ng gi i đoán biên t p b n đ thành qu sau gi i đoán: 1) Ph n m m Ecognition; 2) Ph n m m ArcGIS C th : - S d ng l p b n đ sau gi i đoán ch a làm tr n đ ng bao lô c a xã Qu ng Thành, Th xã Gia Ngh a, t nh k Nông (l p b n đ 1) - S d ng ch c n ng Export results/Polygon smoothed ph n m m Ecognition đ làm ng tr n đ ng bao lô c a l p b n đ đ c l p b n đ làm tr n b ng cơng c có s n ph n m m eCognition (l p b n đ 2) - S d ng công c ArcToolbox/Cartography Tools/Generalization/Smooth Polygon ph n m m ArcGIS đ làm tr n đ ng bao lô c a l p b n đ đ c l p b n đ làm tr n b ng cơng c có s n ph n m m ArcGIS (l p b n đ 3) Hai l p b n đ (l p b n đ l p b n đ 3) t o sau làm tr n đ ng bao lô đ c s d ng đ so sánh v i l p b n đ ch a làm tr n ban đ u (l p b n đ 1) đ đánh giá hi u qu thông qua lo i sai s : 1) Sai s di n tích theo lơ; 2) Sai s di n tích theo tr ng thái - Sai s di n tích theo lơ Tính di n tích cho t ng lô l p b n đ (1,2,3) b ng hàm CartesianArea(obj, "sq m") ph n m m Mapinfo Ch n ng u nhiên 10 lô r ng l p b n đ ghi thành l p riêng theo b c: 1) ánh s th t lô c a l p b n đ t đ n h t t ng s 11.503 lô b ng hàm Rowid ph n m m Mapinfo; 2) Ch n ng u nhiên lô 11.503 lô l p b n đ b ng hàm ng u nhiên Random ph n m m Excel; 3) Ch n lơ cịn l i theo ph ng pháp h th ng, c 900 m l y m, tr ng h p ch n đ n cu i danh sách v n ch a đ s m ti n hành ch n ng c l i t m có s th t 1; 4) Các lơ r ng l a ch n đ c ghi thành l p riêng Ch n 10 lô r ng t ng ng l p b n đ 10 lô r ng t ng ng l p b n đ Xác đ nh sai s theo lô c a l p b n đ v i l p b n đ 1, c a l p b n đ v i l p b n đ b ng ph ng pháp bình ph ng nh nh t (2.1) Trong đó: X2 sai s ; n s c p lô so sánh (10 c p); yi di n tích lơ i l p b n đ 1; f(xi) di n tích lơ t ng ng l p b n đ ho c T P CHÍ KHOA H C VÀ CÔNG NGH LÂM NGHI P S 2-2015 39 Qu n lý tài nguyên r ng & Môi tr ng - Sai s di n tích theo tr ng thái r ng đ t khơng có r ng theo 10 lô r ng ch n l p b n đ n i dung nghiên c u Tính di n tích cho t ng lô l p b n đ (1,2,3) b ng hàm CartesianArea(obj, "hectare") ph n m m Mapinfo Xác đ nh sai s theo lô c a l p b n đ v i l p b n đ b ng ph ng pháp bình ph ng nh nh t theo công th c (2.1) Xác đ nh t ng di n tích c a t ng tr ng thái r ng đ t khơng có r ng l p b n đ Trong đó: X2 sai s ; n s c p lô so sánh (10 c p); yi di n tích lơ l p b n đ 1; f(xi) di n tích lơ t ng ng l p b n đ Xác đ nh t ng di n tích c a t ng tr ng thái r ng đ t khơng có r ng t ng ng l p b n đ Xác đ nh sai s di n tích theo tr ng thái c a l p b n đ v i l p b n đ 1, c a l p b n đ v i l p b n đ b ng ph ng pháp bình ph ng nh nh t thơng qua cơng th c (2.1) Trong đó: X sai s ; n s c p tr ng thái so sánh; yi di n tích tr ng thái r ng ho c đ t khơng có r ng l p b n đ 1; f(xi) di n tích tr ng thái t ng ng l p b n đ ho c 2.4.2 Ph ng pháp nghiên c u k thu t làm tr n đ ng lô b ng t h p công c ph n m m ArcGIS T k t qu nghiên c u n i dung th nh t k t h p v i ph ng pháp chun gia tìm hi u cơng c ph n m m ArcGIS đ th c hi n làm tr n đ ng bao lô ánh giá hi u qu làm tr n đ theo k thu t m i t o ng bao lô L p b n đ t o sau làm tr n đ ng bao lô theo k thu t m i (L p b n đ 4) đ c s d ng đ so sánh v i l p b n đ ch a làm tr n ban đ u (L p b n đ 1) đ đánh giá hi u qu thông qua lo i sai s : 1) Sai s di n tích theo lơ; 2) Sai s di n tích theo tr ng thái - Sai s di n tích theo lơ Tính di n tích cho t ng lô l p b n đ b ng hàm CartesianArea(obj, "sq m") ph n m m Mapinfo Ch n 10 lô r ng t 40 ng ng l p b n đ - Sai s theo tr ng thái r ng đ t khơng có r ng Tính di n tích cho t ng lô l p b n đ b ng hàm CartesianArea(obj, "hectare") ph n m m Mapinfo Xác đ nh t ng di n tích c a t ng tr ng thái r ng đ t khơng có r ng l p b n đ Xác đ nh sai s theo tr ng thái c a l p b n đ v i l p b n đ b ng ph ng pháp bình ph ng nh nh t theo cơng th c (2.1) Trong đó: X2 sai s ; n s c p tr ng thái so sánh; yi di n tích tr ng thái r ng ho c đ t khơng có r ng l p b n đ 1; f(xi) di n tích tr ng thái t ng ng l p b n đ III K T QU NGHIÊN C U VÀ TH O LU N 3.1 Hi u qu làm tr n đ ng lơ sau gi i đốn m t s ph n m m thông d ng Trên ph n m m gi i đốn nh thơng d ng đ u có ch c n ng làm m m đ ng lơ sau gi i đốn nh m m c đích h tr cơng tác xây d ng b n đ hi n tr ng sau gi i đoán Nh ng t ng ph n m m gi i đoán nh khác nhau, ch c n ng làm tr n có s khác nh t đ nh Trong n i dung này, báo trình bày hi u qu làm tr n đ ng lô b ng cơng c có s n ph n m m thông d ng là: eCognition ArcGIS v i l p b n đ sau gi i đoán t ph n m m eCognition ch a làm tr n t i xã Qu ng Thành, Th xã Gia Ngh a, t nh k Nông Hi n tr ng r ng đ t ch a có r ng đ a bàn xã đ c t p h p b ng 3.1 T P CHÍ KHOA H C VÀ CÔNG NGH LÂM NGHI P S 2-2015 Qu n lý tài nguyên r ng & Môi tr ng B ng 3.1 K t qu gi i đoán hi n tr ng r ng đ t ch a có r ng xã Qu ng Thành TT Maldlr Ldlr I 10 II Tên ldlr Di n tích (ha) Di n tích r ng đ t ch a có r ng quy ho ch ba lo i r ng 14 txg R ng g t nhiên núi đ t r ng th ng xanh gi u 15 txb R ng g t nhiên núi đ t r ng th ng xanh trung bình 16 txn R ng g t nhiên núi đ t r ng th ng xanh nghèo 17 txk R ng g t nhiên núi đ t r ng th ng xanh nghèo ki t 52 tnk R ng tre n a khác núi đ t 60 rtg R ng tr ng g núi đ t 72 dtr t tr ng núi đ t 78 dt2 t tr ng có g tái sinh 82 dt1 t tr ng khơng có g tái sinh 88 nn t nông nghi p núi đ t Di n tích đ t ngồi quy ho ch ba lo i r ng T ng 3.565,3 3,9 1.449,9 64,2 10,6 4,2 55,6 0,3 62,8 61,4 1.852,4 4.205,3 7.770,6 Hình nh lơ r ng đ t ch a có r ng sau gi i đốn n n nh v tinh SPOT6 đ c minh h a t i hình 3.1 Hình 3.1 Hình nh lơ r ng đ t ch a có r ng sau gi i đoán K t qu s d ng công c Polygon smoothed đ làm tr n đ ng bao lơ r ng đ t khơng có r ng ph n m m eCognition ch ng x p nh v tinh SPOT6 đ c minh h a t i hình 3.2 T P CHÍ KHOA H C VÀ CÔNG NGH LÂM NGHI P S 2-2015 41 Qu n lý tài nguyên r ng & Mơi tr ng Hình 3.2 Hình nh lơ r ng đ t ch a có r ng làm tr n b ng ph n m m eCognition T hình 3.2 th y: V m t hình nh đ ng lô sau làm tr n h u nh lo i b đ c hi n t ng dích d c so v i đ ng lô sau gi i đốn Tuy nhiên, đ ng lơ t o v n nhi u đo n g p khúc đ t ng t d n đ n s không phù h p v i n n nh bên d i K t qu s d ng công c Smooth Polygon ph n m m ArcGIS đ làm tr n đ ng bao lô r ng đ t r ng ch ng x p nh v tinh SPOT6 đ c minh h a t i hình 3.3 Hình 3.3 Hình nh lơ r ng đ t ch a có r ng làm tr n b ng công c Smooth Polygon ph n m m ArcGIS T hình 3.3 th y: v m t hình nh đ ng lơ sau làm tr n lo i b đ c hi n t ng dích d c so v i đ ng lô sau gi i đốn, đ ng lơ bám sát v i ranh gi i hi n tr ng nh v tinh Tuy nhiên, khu v c ti p giáp gi a lô r ng l i b h ho c ch ng đè nh t khu v c ti p giáp gi a lô tr lên, đ kh c ph c hi n 42 t ng đòi h i ph i m t nhi u công s c trình biên t p b n đ sau b s h đánh giá hi u qu làm tr n đ ng ph ng pháp nêu trên, tác gi s lo i sai s : 1) Sai s di n tích theo lơ; di n tích theo tr ng thái K t qu đ p b ng 3.2 b ng 3.3 T P CHÍ KHOA H C VÀ CÔNG NGH LÂM NGHI P S 2-2015 ng lô d ng 2) Sai ct p Qu n lý tài nguyên r ng & Môi tr ng B ng 3.2 Sai s di n tích theo lơ làm tr n b ng ph n m m thông d ng TT 10 T ng Sai s Di n tích lơ gi i đốn (Yi ) (m2) 4.418,3 6.112,0 13.233,4 15.985,5 19.679,0 7.303,1 10.842,3 7.759,6 7.628,5 8.141,7 101.103,4 Di n tích lơ làm tr n ph n m m ecognition (f(xi)) - (m2) 4.409,9 6.138,5 13.206,9 16.003,5 19.662,2 7.281,6 1.0811 7.781,2 7.610,5 8.116,5 101.021,8 ((Yi -f(xi))2 Di n tích lơ làm tr n b ng cơng c Smooth Polygon ph n m m ArcGIS (f(xj)) - (m2) 70,56 702,25 702,25 324,00 282,24 462,25 979,69 466,56 324,00 635,04 6.658,56 665,8 ((Yi -f(xj))2 4.417,7 6.111,7 13.232,5 15.984,1 19.678,5 7.299,1 10.841,3 7.759,0 7.627,5 8.140,6 101.092,0 0,36 0,09 0,81 1,96 0,25 16,00 1,00 0,36 1,00 1,21 129,96 12,9 B ng 3.3 Sai s di n tích theo tr ng thái làm tr n b ng ph n m m thông d ng TT mald lr 10 14 15 16 17 52 60 72 78 82 88 11 T ng ldlr txg txb txn txk tnk rtg dtr dt2 dt1 nn Ngoài LN 3,90 1.449,90 64,20 10,60 4,20 55,60 0,30 62,80 61,40 1.852,40 Di n tích làm tr n ph n m m ecognition (f(xi)) - (ha) 4,20 1.449,00 64,20 10,50 4,20 55,90 0,30 62,60 61,70 1.852,70 4.205,30 7.770,60 Di n tích tr ng thái gi i đốn (Yi) (m2) 0,09 0,81 0,00 0,01 0,00 0,09 0,00 0,04 0,09 0,09 Di n tích làm tr n b ng cơng c Smooth Polygon ph n m m ArcGIS (f(xj)) - (ha) 3,90 1.449,90 64,20 10,70 4,20 55,60 0,30 62,80 61,90 1.851,80 4.205,30 0,00 4.205,30 0,00 7.770,60 1,22 7.770,60 0,62 Sai s ((Yi-f(xi))2 0,11 T k t qu xác đ nh sai s di n tích theo lơ theo tr ng thái làm tr n đ ng lô b ng m t s ph n m m thông d ng th y: Sai s di n tích theo lơ sai s di n tích theo tr ng thái c a ph ng pháp làm tr n đ ng ((Yif(xj))2 0,00 0,00 0,00 0,01 0,00 0,00 0,00 0,00 0,25 0,36 0,06 bao lô thông qua s d ng công c Smooth Polygon ph n m m ArcGIS đ u nh h n so v i ph ng pháp làm tr n đ ng bao lô s d ng ph n m m eCognition T P CHÍ KHOA H C VÀ CÔNG NGH LÂM NGHI P S 2-2015 43 Qu n lý tài nguyên r ng & Môi tr ng 3.2 K thu t làm tr n đ ng lô b ng t h p công c ph n m m ArcGIS T k t qu nghiên c u t i n i dung 3.1 th y: làm tr n đ ng lô b ng công c Smooth Polygon ph n m m ArcGIS cho k t qu đ ng lô sau làm tr n bám sát hi n tr ng nh sai s nh h n so v i th c hi n ph n m m eCognition Tuy nhiên, làm tr n b ng cơng c có nh c m lơ c nh khơng trùng khít v đ ng biên chung Nguyên nhân c a hi n t ng cách qu n lý d li u đ h a ph n m m ArcGIS t o Cách qu n lý d li u đ h a d ng vùng d ng đ ng ph n m m ArcGIS nói riêng ph n m m thu c h th ng thơng tin đ a lý nói chung đ c minh h a t i hình 3.4 Trên hình 3.4 th hi n lơ c nh là: Lô (ACBEA) lô (ADBFA), cách qu n lý d li u đ h a ph n m m ArcGIS nh sau: - D ng vùng: M i lô m t đ i t ng riêng bi t Nh v y làm tr n, khu v c ti p giáp gi a lô (AEB c a lô AFB c a lô 2) đ c nh n d ng thu c đ i t ng nên s làm tr n khác nhau, d n đ n lô sau làm tr n khơng có s phù h p (trùng nhau) khu v c ti p biên - D ng đ ng: M i Polyline m t đ i t ng đ c chia tách khu v c giao nhau, lô hình 3.4 n u d ng đ ng s g m đ i t ng: Polyline (ACB), Polyline (ADB) Polyline (AEB ho c AFB đ c qu n lý đ i t ng) Nh v y làm tr n, khu v c ti p giáp gi a lô (AEB c a lô AFB c a lô 2) đ c nh n d ng đ i t ng nên đ ng bao lô sau làm tr n trùng khu v c ti p biên Do đó, tác gi nh n th y c n ph i chuy n l p lô d ng vùng thành l p lơ d ng đ ng sau làm tr n l p lô d ng đ ng t o chuy n l i thành d ng vùng đ s d ng Hình 3.4 Cách qu n lý d li u d ng vùng d ng đ 44 ng ph n m m ArcGIS T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 Qu n lý tài nguyên r ng & Môi tr chuy n l p ranh gi i lơ sau gi i đốn t d ng vùng thành d ng đ ng ph n m m ArcGIS có th s d ng m t s cơng c sau: QT1) S d ng công c Feature To Line (ArcToobox\Data Management Tools\Features\Feature To Line); QT2) S d ng công c Polygon To Line (ArcToobox\Data Management Tools\Features\Polygon To Line); QT3) S Coverage d ng công c Feature Class To (ArcToobox\Conversion Tools\To Coverage\Feature Class To Coverage) ng L p lô d ng đ ng t o s đ c làm tr n b ng: QT4) S d ng công c Smooth Line (ArcToobox\Cartography Tools\Generalization\Smooth Line) Th c hi n h t QT4, đ ng bao lô đ c làm tr n, nhiên l p b n đ v n d ng đ ng, c n ph i chuy n thành d ng vùng đ s d ng chuy n l p b n đ d ng đ ng thành d ng vùng ArcGIS s d ng công c : Feature To Polygon (QT5) (ArcToobox\Data Management Tools\Features\Feature To Polygon) K t qu th c hi n đ c minh h a t i hình 3.5 Hình 3.5 Hình nh lơ r ng đ t ch a có r ng sau làm tr n b ng t h p công c ph n m m ArcGIS T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 45 Qu n lý tài ngun r ng & Mơi tr ng T hình 3.5 th y: V m t hình nh đ ng lô sau làm tr n lo i b đ c hi n t ng dích d c so v i đ ng lơ sau gi i đốn, đ ng lô bám sát v i ranh gi i hi n tr ng nh v tinh, khơng cịn hi n t ng tr ng đè ho c h t i khu v c ti p giáp gi a lô r ng đánh giá hi u qu làm tr n đ ng lô b ng t h p công c ph n m m ArcGIS, tác gi s d ng lo i sai s : 1) Sai s di n tích theo lơ; 2) Sai s di n tích theo tr ng thái K t qu đ c t p h p b ng 3.4 b ng 3.5 B ng 3.4 Sai s di n tích theo lô làm tr n b ng t h p công c ph n m m ArcGIS TT Di n tích lơ gi i đốn (Yi) (m2) 10 T ng Sai s Di n tích lơ làm tr n b ng t h p công c ph n m m ArcGIS (f(xj)) - (m2) 4.418,3 6.112,0 13.233,4 15.985,5 19.679,0 7.303,1 10.842,3 7.759,6 7.628,5 8.141,7 101.103,4 ((Yi-f(xj))2 4.418,3 6.114,2 13.233,9 15.987,1 19.678,5 7.303,7 10.844,7 7.760,7 7.630,5 8.141,5 101.113,1 4,84 0,25 2,56 0,25 0,36 5,76 1,21 0,04 94,09 9,4 B ng 3.5 Sai s di n tích theo tr ng thái làm tr n b ng t h p công c ArcGIS TT maldlr ldlr 10 11 T ng Sai s 14 15 16 17 52 60 72 78 82 88 txg txb txn txk tnk rtg dtr dt2 dt1 nn Ngồi LN Di n tích tr ng thái gi i đoán (Yi) (m2) 3,90 1.449,90 64,20 10,60 4,20 55,60 0,30 62,80 61,40 1.852,40 4.205,30 7.770,60 Nh v y, k t h p công c ph n m m ArcGIS đ làm tr n đ ng lơ sau gi i đốn cho k t qu sai s di n tích theo lơ sai s di n tích theo tr ng thái nh h n so v i 46 Di n tích làm tr n b ng t h p công c ArcGIS (f(xi)) - (ha) 3,9 1.450,3 64,2 10,6 4,2 55,6 0,3 62,8 61,4 1852 4.205,3 7.770,6 ((Yi-f(xi))2 0,16 0 0 0 0,16 0,32 0,03 ph ng pháp làm tr n đ ng lô b ng công c Smooth Polygon ph n m m ArcGIS ho c ph ng pháp làm tr n đ ng bao lô s d ng ph n m m eCognition T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 Qu n lý tài nguyên r ng & Môi tr IV K T LU N - S d ng ph n m m eCognition đ làm tr n đ ng lơ sau gi i đốn cho k t qu không t t v m t đ h a đ ng lơ sau làm tr n khơng bám sát v i hi n tr ng c a nh - S d ng công c Smooth Polygon ph n m m ArcGIS đ làm tr n đ ng lơ sau gi i đốn s xu t hi n hi n t ng ch ng đè h gi a lơ g n có chung đ ng biên - K t h p công c ph n m m ArcGIS đ làm tr n đ ng lơ sau gi i đốn cho k t qu sai s di n tích theo lơ (9,4) sai s di n tích theo tr ng thái (0,03) nh h n so v i vi c s d ng cơng c có s n eCognition ho c ArcGIS đ ng lô sau làm tr n bám sát v i hi n tr ng c a nh - Quy trình làm tr n đ ng bao lơ b ng t h p công c ph n m m ArcGIS + Chuy n l p b n đ sau gi i đoán t d ng vùng thành d ng đ ng b ng m t công c : Feature To Line (QT1); Polygon To Line (QT2) ho c; Feature Class To Coverage (QT3) + Làm m m l p b n đ d ng đ công c Smooth Line (QT4) ng b ng ng + Chuy n l p đ ng thành l p vùng b ng công c Feature To Polygon (QT5) TÀI LI U THAM KH O Tr n Quang B o, Nguy n V n Th , Ph m V n Du n (2014), ng d ng GIS qu n lý tài nguyên thiên nhiên, Nhà xu t b n Nông nghi p, Hà N i Ph m V n Du n, Phùng V n Khoa (2013), Th nghi m ph ng pháp xây d ng b n đ ki m kê r ng l u v c t nh v tinh SPOT5, T p chí Khoa h c Lâm nghi p s 1, n m 2013, trang 2619 – 2630 Luciana F Alves, Simone A Vieira, Marcos A Scaranello, Plinio B Camargo, Flavio A.M Santos, Carlos A Joly, Luiz A Martinelli (2010) Forest structure and live aboveground biomass variation along an elevational gradient of tropical Atlantic moist forest (Brazil) Forest Ecology and Management, pp 679-691 R.J Hall, R.S Skakun, E.J Arsenault, B.S Case (2006), Modeling forest stand structure attributes using Landsat ETM+ data: Application to mapping of aboveground biomass and stand volume Forest Ecology and Management, No 225, pp 378-390 Peter T Wolter, Philip A Townsend, Brian R Sturtevant (2009), Estimation of forest structural parameters using and 10 meter SPOT5 satellite data Remote Sensing of Environment, No 113, pp 2019-2036 Foody, G.M., Boyd, D.S and Cutler, M.E.J (2003), Predictive relation of tropical forest biomass from Landsat TM data and their transferability between regions Remote Sensing of Environment, No 85, pp 463-474 T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 47 Qu n lý tài nguyên r ng & Môi tr ng TEST METHOD SMOOTHING ROAD OF PLOTS AFTER INTERPRETATION FROM IMAGE SATELLITE Pham Van Duan, Vu Thi Thin, Nguyen Quang Giap SUMMARY Remote sensing technology has proven advantages and prospects for forest inventory, especially mapping and monitoring of the quality forest in our country One of the jobs lost time after interpretation is smooth plots to editor mapping the status of the forest The article uses some specialized software, map layer after interpreted from satellite imagery SPOT6 in Quang Thanh commune , Gia Nghia town, Dak Nong province, maps and supporting documents the plots smoothed by three methods: 1) By tool available on eCognition software; 2) By tool available on ArcGIS software; 3) By combining tools in ArcGIS software The results show the wrong number of identified by the method of least squares of mapping method with the after map interpreted at 9.4 for forest plots and 0.03 at the status forest These values are smaller than the wrong number when using method (wrong number plot: 665.8; wrong number status of forest: 0.11) or method (wrong number plot: 12.9; wrong number status of forest: 0,06) From the research results have identified the technical steps smoothing plots: 1) Moving map layer after interpretation from regions to lines ; 2) To smoothing lines; 3) Move the line layer after smoothed into region and update the data to use Keywords: After interpretation, ArcGIS software, Smooth plot line, SPOT6, State Forest Ng i ph n bi n Ngày nh n Ngày ph n bi n Ngày quy t đ nh đ ng 48 : PGS.TS Tr n Quang B o : 18/3/2015 : 17/4/2015 : 09/6/2015 T P CHÍ KHOA H C VÀ CƠNG NGH LÂM NGHI P S 2-2015 ... hi n lô c nh là: Lô (ACBEA) lô (ADBFA), cách qu n lý d li u đ h a ph n m m ArcGIS nh sau: - D ng vùng: M i lô m t đ i t ng riêng bi t Nh v y làm tr n, khu v c ti p giáp gi a lô (AEB c a lô AFB... v y làm tr n, khu v c ti p giáp gi a lô (AEB c a lô AFB c a lô 2) đ c nh n d ng đ i t ng nên đ ng bao lô sau làm tr n trùng khu v c ti p biên Do đó, tác gi nh n th y c n ph i chuy n l p lô d... K thu t làm tr n đ ng lô b ng t h p công c ph n m m ArcGIS T k t qu nghiên c u t i n i dung 3.1 th y: làm tr n đ ng lô b ng công c Smooth Polygon ph n m m ArcGIS cho k t qu đ ng lô sau làm tr