tµi liÖu híng dÉn
Th¸ng 12 n¨m 2005
Môc lôc Th¸ng 12 n¨m 2005..................................................................................1 Giíi thiÖu chung......................................................................................3 PhÇn 1: Thèng kª vµ ph©n tÝch sù biÕn ®éng t×nh tr¹ng sö dông ®Êt ®ai 4
0.1 Gi¶i ®o¸n ¶nh sè vÖ tinh ®a phæ b»ng phÇn mÒm ERDAS IMAGINE 8.4 trong x©y dùng b¶n ®å hiÖn tr¹ng rõng víi hai thêi kú:...................................... .......................................5
0.1.1 Bíc 1: NhËp ¶nh tõ CD vµo m¸y tÝnh............5 0.1.2 Bíc 2: Tæng hîp mµu gi¶..............................7 0.1.3 Bíc 3: N¾n chØnh to¹ ®é.............................8 0.1.4 Bíc 4: T¨ng cêng chÊt lîng ¶nh.....................9 0.1.5 Bíc 5: C¾t ¶nh..........................................11 0.1.6 Bíc 6: LÊy mÉu trªn ¶nh ®Ó ph©n lo¹i........12 0.1.7 Bíc 7: Ph©n lo¹i cã gi¸m ®Þnh...................15 0.1.8 Bíc 8: Läc ¶nh sau ph©n lo¹i......................16 0.1.9 Bíc 9: KiÓm tra, chØnh söa sau ph©n lo¹i...16 0.1.10 Bíc 10: Gép nhãm mÉu- nhãm ®èi tîng ... .18 0.1.11 ¶nh sau ph©n lo¹i kiÓm tra, gép nhãm.....18 0.1.12 ChuyÓn ®æi format ®Ó chång xÕp, xö lý b¶n ®å .............................................................19 0.1.13 ChuyÓn raster file.grid sang vecter file.shp ........................................................................20 0.1.14 ChuyÓn file.shp sang file.tab ®Ó xö lý b»ng Map/Infor..........................................................21 0.2 Chång xÕp b¶n ®å hiÖn tr¹ng cña hai thêi kú b»ng phÇn mÒm ARC VIEW:................................................ ...............21 1 PhÇn II: Ph©n tÝch, thèng kª c¸c lo¹i ®Êt ®ai, tµi nguyªn rõng theo ®ai cao, ®é dèc .....................................................................................30
1.1 X©y dùng m« h×nh sè ho¸ ®Þa h×nh DTM:.................30 1
1.2 Chång xÕp b¶n ®å ®ai cao vµ b¶n ®å hiÖn tr¹ng:........43 1.3 X©y dùng b¶n ®å ®é dèc:.................................... .......45 1.4 Chång xÕp b¶n ®å ®é dèc vµ b¶n ®å hiÖn tr¹ng:........48 2 PhÇn III: øng dông c«ng nghÖ viÔn th¸m vµ GIS vµo viÖc x¸c ®Þnh khu vùc u tiªn trång rõng..............................................................................48
2.1 X¸c ®Þnh c¸c chØ tiªu tham gia vµo viÖc x¸c ®Þnh vïng u tiªn cho trång rõng:..................................................... .....48 2.2 X¸c ®Þnh träng sè cho c¸c nh©n tè tham gia:...............50 2.3 Cho ®iÓm cho tõng b¶n ®å ........................................53 2.4 Chång xÕp c¸c líp b¶n ®å:........................ ...................57 2.5 Ph©n vïng u tiªn trång rõng:.......................................57
2
Giíi thiÖu chung
Kü thuËt ViÔn th¸m øng dông trong c«ng t¸c ®iÒu tra rõng hiÖn nay ®ang ®îc ¸p dông réng r·i trªn ph¹m vi c¶ níc bëi nh÷ng u thÕ râ rÖt cña nã. Trong ®ã, ®Æc biÖt lµ gi¶m ®¸ng kÓ chi phÝ tµi chÝnh vµ c«ng lao ®éng trong c«ng t¸c ®iÒu tra, x©y dùng b¶n ®å hiÖn tr¹ng rõng còng nh theo dâi, ®¸nh gi¸ diÔn biÕn tµi nguyªn rõng. Bªn c¹nh ®ã, kh¶ n¨ng cËp nhËt th«ng tin tõ t liÖu ¶nh ViÔn th¸m lµ rÊt nhanh chãng. Ngoµi ra, c¸c lo¹i ¶nh ViÔn th¸m còng kh«ng ngõng ®îc n©ng cao vÒ chÊt lîng, ®é ph©n gi¶i kh«ng gian vµ gi¸ c¶ còng gi¶m ®¸ng kÓ, phï hîp víi kh¶ n¨ng tµi chÝnh cña nhiÌu ®Þa ph¬ng. Song song víi kh¶ n¨ng øng dông cña ViÔn th¸m, hiÖn nay c¸c s¶n phÈm c«ng nghÖ th«ng tin chuyªn nghiÖp nh HÖ thèng th«ng tin ®Þa lý (GIS) ph¸t triÓn rÊt m¹nh mÏ, gãp phÇn vµo viÖc gi¶i quyÕt c¸c vÊn ®Ò vÒ nhËp, lu tr÷, ph©n tÝch, sö lý d÷ liÖu theo c¸c môc ®Ých, nhu cÇu kh¸c nhau. Tuy nhiªn, bªn c¹nh sù ph¸t triÓn m¹nh cña c«ng nghÖ ViÔn th¸m vµ GIS, trong thùc tÕ vÉn cã nh÷ng kho¶nh c¸ch rÊt lín gi÷a c¸c nhµ chuyªn m«n vµ viÖc øng dông c«ng nghÖ trong s¶n suÊt, hiÖn nay nhiÒu nhµ chuyªn m«n cha cã nhiÒu c¬ së sö dông hîp lý vµ ph¸t huy thÕ m¹nh cña c«ng nghÖ míi nµy. Do ®ã, viÖc n©ng cao n¨ng lùc kü thuËt viÔn th¸m, GIS vµ øng dông chóng trong viÖc ®iÒu tra vµ quy ho¹ch rõng lµ môc tiªu rÊt quan träng cña kho¸ huÊn luyÖn nµy §Ó ®¹t ®îc môc tiªu cña khãa häc, chóng t«i sÏ giíi thiÖu c¸c néi dung chÝnh sau: 1. Theo dâi ®¸nh gi¸ sù biÕn ®éng t×nh tr¹ng sö dông ®Êt vµ tµi nguyªn rõng b»ng ¶nh sè vÖ tinh ®a phæ vµ ®a thêi gian hoÆc b»ng b¶n ®å hiÖn trang sö dông ®Êt qua nhiÒu thêi kú. 2. Ph©n tÝch vµ thèng kª c¸c lo¹i ®Êt ®ai, tµi nguyªn rõng theo c¸c cÊp ®é cao vµ ®é dèc lµm c¬ së cho c«ng t¸c quy ho¹ch sö dông ®Êt ®ai. 3. øng dông c«ng nghÖ viÔn th¸m vµ GIS vµo viÖc ph©n vïng u tiªn trång rõng
3
PhÇn 1: Thèng kª vµ ph©n tÝch sù biÕn ®éng t×nh tr¹ng sö dông ®Êt ®ai
Nh chóng ta ®· biÕt, muèn theo dâi ®¸nh gi¸ sù biÕn ®éng t×nh tr¹ng sö dông ®Êt ®ai vµ tµi nguyªn rõng cña mét khu vùc trong giai ®o¹n thêi gian nµo ta cÇn ph¶i tiÕn hµnh x©y dùng chuçi c¸c b¶n ®å hiÖn tr¹ng sö dông ®Êt ®ai cña khu vùc theo c¸c thêi kú trong giai ®o¹n ®ã b»ng ph¬ng ph¸p gi¶i ®o¸n ¶nh sè vÖ tinh ®a phæ, ®a thêi gian. C¸c häc viªn ®· ®îc häc vÒ quy tr×nh x©y dùng b¶n ®å hiÖn tr¹ng sö dông ®Êt ®ai b»ng ph¬ng ph¸p gi¶i ®o¸n ¶nh sè vÖ tinh ®a phæ, nªn cã thÓ tãm t¾t c¸c bíc trong ph¬ng ph¸p theo dâi ®¸nh gi¸ sù biÕn ®éng t×nh tr¹ng sö dông ®Êt ®ai vµ tµi nguyªn rõng b»ng ¶nh vÖ tinh ®a phæ vµ ®a thêi gian nh sau :
4
0.1 Gi¶i ®o¸n ¶nh sè vÖ tinh ®a phæ b»ng phÇn mÒm ERDAS IMAGINE 8.4 trong x©y dùng b¶n ®å hiÖn tr¹ng rõng víi hai thêi kú: Trong kho¸ huÊn luyÖn nµy c¸c häc viªn sÏ ®îc híng dÉn sö dông phÇn mÒm ERDAS IMAGINE 8.4 gi¶i ®o¸n ¶nh sè vÖ tinh Landsat ®Ó thµnh lËp b¶n ®å hiÖn tr¹ng rõng vµ sö dông ®Êt. Khu vùc nghiªn cøu ®îc chän lµ huyÖn Th¹ch An. 0.1.1 Bíc 1: NhËp ¶nh tõ CD vµo m¸y tÝnh ¶nh gèc Landsat ETM ®îc lu tr÷ trªn ®Üa CD díi format *.HDF nhng cha thÓ xö lý ngay trªn m¸y tÝnh,v× vËy tríc hÕt ta ph¶i nhËp c¸c kªnh ¶nh K1, K2, K3 vµ K4 tõ ®Üa CD vµo trong m¸y tÝnh, tøc lµ chuyÓn tõ Format *. HDF sang *.IMG ¶nh gèc cã 8 kªnh nhng ®Ó sö dông cho b¶n ®å rõng chñ yÕu dïng K1, k2, k3 vµ k4 nªn chØ nhËp c¸c kªnh nµy. Thùc hiÖn tr×nh tù nh sau: Cho ®Üa vµo æ CD Rom -> më æ E -> file MTL, file REF ®Ó xem sè hµng (sample) vµ sè cét (colums) cña ¶nh, ghi sè hµng cét nµy ®Ó khai b¸o khi nhËp (Import).
5
Khëi ®éng phÇn mÒm ERDAS Trªn Menu chÝnh chon: Import Type: Generic binary Media: File Ghi input file name (®a tªn file ®Çu vµo cña kªnh ¶nh cÇn nh¹p) Ghi out put file name (®Æt tªn file ®Çu ra cña kªnh ¶nh t¬ng øng)
VÝ dô: Kªnh cÇn nhËp lµ 127-46-K3, khai th«ng sè hµng cét -> unsigned 8bit ->Type: Billary ->OK Lµm lÇn lît cho tÊt c¶ c¸c kªnh ¶nh cÇn nhËp b»ng c¸ch t¬ng tù. Trong trêng hîp cña chóng ta lµ nhËp 4 kªnh ¶nh K1, k2, k3 Vµ k4.
6
0.1.2 Bíc 2: Tæng hîp mµu gi¶ C¸c kªnh ¶nh ®¬n lÎ k1, k2, k3, k4 ®Òu lµ nh÷ng ¶nh toµn s¾c nªn cha thÓ gi¶i ®o¸n ¶nh ngay trªn c¸c kªnh ®¬n lÎ ®ã ®îc. §Ó m¾t ta nhËn biÕt ®èi tîng râ vµ chÝnh x¸c h¬n, ta cÇn tiÕn hµnh tæ hîp mÇu b»ng c¸ch kÕt hîp 3 kªnh cïng mét lóc. C¸c bíc thùc hiÖn trªn m¸y vi tÝnh nh sau: Trªn Menu chÝnh chän: Interpreter Utilities Layer Stack Input File: Ên con trá vµo biÓu tîng më file ®Ó lÇn lît ®a tªn file cña c¸c kªnh ¶nh ®¬n vµo. Sau mçi lÇn ®a 1 kªnh vµo, Ên Nót “Add” Output File: ®Æt tªn file tæng hîp mµu. OK
7
0.1.3 Bíc 3: N¾n chØnh to¹ ®é VÒ mÆt vÞ trÝ, ¶nh tæng hîp mµu míi cã c¸c gi¸ trÞ hµng cét, cha cã to¹ ®é vµ bÞ lÖch so víi b¶n ®å ®Þa h×nh, v× vËy cÇn g¾n to¹ ®é cho ¶nh vµ n¾n chØnh ®Ó cã thÓ sö dông. Cã nhiÒu ph¬ng ph¸p n¾n chØnh h×nh häc. C¸ch thø nhÊt lµ n¾n ¶nh dùa vµo file b¶n ®å (®Þa h×nh) ®· ®îc sè ho¸. C¸ch lµm nh sau: Trong Menu chÝnh chän Datapresentation Image Geometric corection From image File View Open Raster Tªn file ¶nh cÇn n¾n Ok Polinominal OK Projection Add change projection UTM WGS 84 North Zone 48(49) OK Map unit: m Set Projection from GCP Tool Vector layer Tªn File s«ng suèi cã to¹ ®é chuÈn (format *.Shp) Out put file name: ®Æt tªn file ®· n¾n (VÝ dô: 127-45nan.IMG) Nearest neighbour BÊm con trá vµo biÓu tîng h×nh b¸nh chng trªn thanh c«ng cô, 8
BÊm con trá vµo h×nh trßn g¹ch chÐo trªn thanh c«ng cô ®Ó lÊy ®iÓm khèng chÕ. BÊm vµo ®iÓm khèng chÕ trªn ¶nh sau ®ã bÊm con trá vµo ®iÓm t¬ng øng trªn file s«ng suèi. Sau khi cã ba ®iÓm khèng chÕ ®Çu tiªn, sö dông biÓu tîng h×nh trßn g¹ch chÐo trªn b¶ng thèng kª ®Ó lÊy ®iÓm khèng chÕ tiÕp theo. Lóc nµy bÊm vµo 1 ®iÓm trªn ¶nh th× ®iÓm ®ã sÏ xuÊt hiÖn t¬ng øng trªn b¶n ®å s«ng suèi. Dïng con trá ®Ó di chuyÓn ®iÓm khèng chÕ vÒ ®óng vÞ trÝ t¬ng øng. Chó ý: - Nªn phãng to ¶nh tæng hîp mµu vµ file s«ng suèi ®Ó lÊy ®iÓm khèng chÕ chÝnh x¸c. - Nªn ®æi mµu nÒn b¶n ®å s«ng suèi sang mµu tèi ®Ó thÊy râ c¸c ®iÓm khèng chÕ (mµu tr¾ng). - LÊy c¸c ®iÓm khèng chÕ ph©n bè ®Òu trªn tê ¶nh, sè lîng cµng nhiÒu cµng tèt.
0.1.4 Bíc 4: T¨ng cêng chÊt lîng ¶nh 0.1.4.1 T¨ng cêng ®é t¬ng ph¶n: T¨ng cêng chÊt lîng ¶nh lµ sù gi·n ®é t¬ng ph¶n, ®é s¸ng tèi. ¶nh sau khi tæng hîp mµu cã ®é t¬ng ph¶n cßn thÊp, do vËy cÇn d·n ®é t¬ng ph¶n ®Ó dÔ ph©n biÖt c¸c lo¹i ®èi tîng kh¸c nhau. C¸ch lµm nh sau: Trªn mµn h×nh Menu chÝnh chän: View File Open Raster Tªn File ¶nh OK Raster Fit image to Window Raster Contract 9
General contract Brightness/ Contract Apply Save change? Yes 0.1.4.2T¨ng cêng ®é ph©n gi¶i kh«ng gian: §èi víi ¶nh Landsat7-ETM+, riªng kªnh toµn s¾c Panchromatic (K8), ®é ph©n gi¶i kh«ng gian lµ 15 m, trong khi ®ã c¸c kªnh ®a phæ kh¸c l¹i cã ®é ph©n gi¶i kh«ng gian lµ 30 m. ChÝnh v× vËy, cÇn ph¶i t¨ng cêng ®é ph©n gi¶i cho ¶nh tæng hîp mµu cã ®é ph©n gi¶i 30 m lªn 15 m. Intepreter Spatial Enhancement Resolution Merge High Resolution Image File: Chän kªnh toµn s¾c Panchromatic (K8) cã ®é ph©n gi¶i kh«ng gian lµ 15 m. Multispectral Input File: Chän ¶nh tæng hîp mµu Output File: §Æt tªn file ¶nh ®îc t¨ng cêng lªn 15 m. Method: Chän Principle component Resampling Technique: chän Nearest Neighbor Output: Chän Unsigned 8 bit OK
10
0.1.5 Bíc 5: C¾t ¶nh Mét c¶nh ¶nh Landsat thêng cã kÝch thíc lµ 170x185 km vµ dung lîng ¶nh vµo kho¶ng 500-600 MB v× vËy g©y khã kh¨n trong viÖc xö lý vµ m¸y tÝnh hay bÞ trôc trÆc. Trong trêng hîp nÕu ta kh«ng cÇn dïng c¶ c¶nh ¶nh v× khu vùc ta quan t©m nhá h¬n mét c¶nh ¶nh th× ta cã thÓ chØ c¾t lÊy khu vùc ta cÇn. ViÖc c¾t ¶nh rÊt cã lîi v× tèc ®é xö lý vµ gi¶i ®o¸n ¶nh nhanh h¬n vµ m¸y tÝnh cïng Ýt bÞ trôc trÆc. C¸ch c¾t ¶nh nh sau: Tõ Menu chÝnh chän: View File Open Raster Tªn file ¶nh gèc OK Fit image to Window
11
Utility Inquire box (Type: File) Fit to AOI Apply Interpreter Utility Sub set In put (File ¶nh gèc, ®¸nh dÊu vßng trßn File) Out put (tªn file ¶nh c¾t ra) From Inquire box – OK
0.1.6 Bíc 6: LÊy mÉu trªn ¶nh ®Ó ph©n lo¹i C«ng viÖc lÊy mÉu trªn ¶nh phôc vô cho viÖc ph©n lo¹i cã gi¸m ®Þnh. §©y lµ c«ng ®o¹n rÊt quan träng ¶nh hëng ®Õn ®é chÝnh x¸c cña viÖc gi¶i ®o¸n nªn cÇn thùc hiÖn hÕt søc cÈn thËn. §Ó lµm ®îc ®iÒu nµy, gi¶i ®o¸n viªn kh«ng chØ ®¬n thuÇn biÕt sö dông phÇn mÒm vi tÝnh, cã kho¸ mÉu ¶nh qua kh¶o s¸t thùc tÕ, mµ cßn cÇn cã kiÕn thøc l©m nghiÖp. §Æc biÖt lµ hiÓu biÕt vÒ hÖ thèng ph©n lo¹i ®Êt ®ai, tµi 12
nguyªn rõng cña ViÖt Nam vµ cã kiÕn thøc ®Þa ph¬ng. C¸ch lÊy mÉu ¶nh nh sau: Më ¶nh cÇn ph©n lo¹i b»ng c¸ch: Trªn Menu chÝnh chän: View File Open Chän tªn file ¶nh cÇn ph©n lo¹i C¸ch lÊy mÉu ¶nh: Classifer Signature editor AOI (area of interest- vïng quan t©m) Tool Dïng biÓu tîng h×nh ®a gi¸c dÝch d¾c ®Ó chän mÉu BÊm con trá vµo biÓu tîng + ( add) ®Ó nhËp vµo b¶ng mÉu Class name: ®Æt tªn lo¹i ®Êt, lo¹i rõng. KÝch chuét ph¶i ®Ó chän mµu cho tõng mÉu, nhãm mÉu (vÝ dô nhãm ®Êt n«ng nghiÖp chän mµu vµng, nhãm ®Êt cã rõng giµu chän mµu xanh ®Ëm...)
13
Chó ý : Nhãm ®Êt cã rõng trung b×nh, dµy: Ký hiÖu mÉu: RG1, RG2,....Mµu xanh l¸ c©y ®Ëm Rõng tha (phôc håi, non): RN1, RN2...., Mµu xanh l¸ c©y nh¹t §Êt trèng cã c©y (Ic, Ib) : Ký hiÖu T1, T2....Mµu n©u nh¹t §Êt trèng kh«ng c©y (Ia) ®Êt NN, D©n c, §Êt kh¸c: Ký hiÖu N1, N2...Mµu vµng chanh Bãng nói ®en, M©y, bãng m©y: Ký hiÖu K1, K2...Mµu tr¾ng S«ng hå: Ký hiÖu S1, S2...Mµu xanh níc biÓn ®Ëm §©y lµ nhãm hiÖn tr¹ng c¬ së, sau khi ph©n lo¹i xong, tiÕn hµnh läc, kiÓm tra tõng líp mÉu ®Ó chØnh söa vµ gép nhãm sau. Chän xong nhí ghi file AOI vµ file *.sig vµo th môc ®ang xö lý ¶nh ®Ó sö dông l¹i khi cÇn. §Æt tªn file theo tªn tê ¶nh. VÝ dô: 127-45.aoi; 12745.sig LÊy tèi thiÓu mçi lo¹i ®èi tîng rõng trªn b¶n ®å kiÓm kª 10 mÉu trë lªn. Sè mÉu lÊy tû lÖ víi diÖn tÝch cña ®èi tîng trªn ¶nh. Kh«ng lÊy mÉu rõng trång 14
MÉu kh«ng nªn qu¸ lín, cÇn lÊy ®ång nhÊt trong 1 mÉu. 0.1.7 Bíc 7: Ph©n lo¹i cã gi¸m ®Þnh Classifer Supervised classification Input file (®a tªn ¶nh cÇn ph©n lo¹i) Out put file (®Æt tªn ¶nh sÏ ph©n lo¹i) Signature file (®a tªn file mÉu ®· ®Æt ë phÇn 6. VÝ dô 12447.sig) Attribute option: Mean Parametric rule: Maximum likelihood OK
Muèn xem ¶nh ®· ph©n lo¹i lµm nh sau: View file Open (raster) 15
Tªn file ®· ph©n lo¹i OK BÊm chuét ph¶i Fit image to window
0.1.8 Bíc 8: Läc ¶nh sau ph©n lo¹i Khi ¶nh ph©n lo¹i song thêng cã d¹ng “muèi tiªu”, tøc lµ gåm c¸c pixel nhá cña c¸c lo¹i ®Êt ®ai, c¸c lo¹i rõng xen kÏ nhau. Sù xen kÏ nµy cã thÓ lµ ®óng trong thùc tÕ nhng l¹i rÊt khã kiÓm tra vµ sö dông v× kh«ng thÓ x¸c ®Þnh ®îc ranh giíi gi÷a c¸c l« nhá ®ã v× thÕ cÇn lo¹i bít chóng. C¸ch thùc hiÖn nh sau: Trªn Menu chÝnh chän: Interpreter GIS analysis Neighbourhood Input file (®a tªn ¶nh sè cÇn läc) Out put file (®Æt tªn file sau läc) §¸nh dÊu chÊm vµo chç File; Out put: unsigned 8 bit Funtion: Majority Size: 3 x 3 -> OK
0.1.9 Bíc 9: KiÓm tra, chØnh söa sau ph©n lo¹i 0.1.9.1 KiÓm tra tæng thÓ §Ó kiÓm tra tæng thÓ ¶nh ®· ph©n lo¹i, tiÕn hµnh më ®ång thêi ¶nh ®· ph©n lo¹i vµ läc víi ¶nh cha ph©n lo¹i chång lªn nhau ®Ó kiÓm tra. C¸ch lµm nh sau: Trªn Menu chÝnh chon: View Open
16
File ¶nh cha ph©n lo¹i Ok Chuét ph¶i Fit image to window Open File ¶nh ®· ph©n lo¹i vµ ®· läc Display option: bá dÊu ë clear display OK Utility Swipe Dïng con trá di chuyÓn thanh c«ng cô Swipe ®Ó so s¸nh ¶nh sè vµ ¶nh gèc trªn cïng vÞ trÝ xem ®· ®¹t yªu cÇu cha, cßn sai ë nh÷ng ®èi tîng nµo cÇn bæ sung mÉu ®Ó ph©n lo¹i l¹i. 0.1.9.2 KiÓm tra tõng líp mÉu Thùc chÊt phÇn nµy lµ kiÓm tra tõng tr¹ng th¸i rõng hoÆc tõng lo¹i ®Êt ®ai ®îc gi¶i ®o¸n cã ®óng so víi ¶nh gèc hay kh«ng. Theo chñ quan cña ngêi gi¶i ®o¸n, nÕu thÊy ¶nh gi¶i ®o¸n cha tin cËy ®îc th× cÇn tiÕn hµnh lÊy mÉu bæ xung hoÆc lÊy mÉu l¹i cho tr¹ng th¸i cô thÓ ®ã. C¸ch lµm nh sau: Trªn Menu chän: View Raster Attribute M¸y sÏ hiÖn b¶ng mµu c¸c líp ®· lÊy mÉu vµ ph©n lo¹i. Tuy cã nhiÒu líp (mçi mÉu t¬ng ®¬ng 1 líp) nhng cã mµu theo nhãm nªn ®Ó kiÓm tra tõng líp tiÕn hµnh kÝch chuét vµo líp cÇn kiÓm tra theo vÞ trÝ rows (ngoµi cïng cña b¶ng), sau ®ã nh¸y chuét vµo mµu cò, ®æi thµnh mµu ®á. Trªn ¶nh mµu ®á cña líp nµy sÏ xuÊt hiÖn. KiÓm tra xem viÖc ph©n lo¹i líp nµy ®· hîp lý cha (theo c¸ch ë môc a). NÕu ph©n lo¹i cña líp mÉu cÇn kiÓm tra ®¹t yªu cÇu th× chuyÓn l¹i mµu cò ®Ó kiÓm tra mÉu kh¸c. NÕu mÉu cha ®¹t yªu cÇu, xÐt xem nã thuéc vÒ líp nµo kh¸c th× ®æi tªn vµ mµu vÒ líp ®ã. TËp trung kiÓm tra nh÷ng ®èi tîng n»m ë ranh giíi cã rõng, kh«ng rõng vµ c¸c chç cã bãng nói v× ë ®ã hay bÞ nhÇm lÉn, ...VÝ dô c¸c mÉu Ic, rõng phôc håi, rõng nghÌo... 17
0.1.9.3LÊy mÉu bæ sung, ph©n lo¹i l¹i C¸h lµm nh sau: Trªn Menu chän View Open file ¶nh gèc Open file *.SIG Open file *.AOI Chän: Append OK TiÕn hµnh lÊy mÉu bæ sung ë nh÷ng líp cÇn thiÕt, nh÷ng líp cßn sai sè lín. Sau khi lÊy mÉu bæ sung, tiÕn hµnh ph©n lo¹i nh môc bíc 7 vµ kiÓm tra l¹i nh bíc 8. 0.1.10Bíc 10: Gép nhãm mÉu- nhãm ®èi tîng ¶nh sau ph©n lo¹i vµ kiÓm tra tuy chØ cã mét sè mµu ®Æc trng cho tõng nhãm mÉu nhng thùc chÊt trong tõng nhãm l¹i cã rÊt nhiÒu m· hiÖu kh¸c nhau (mçi mÉu mét m· hiÖu), do vËy cÇn gép chóng l¹i ®Ó gi¶m sè polygon khi chång xÕp, tÝnh to¸n. M· hiÖu míi cã thÓ theo m· quy ®Þnh cña c¸c lo¹i ®Êt ®ai do ViÖn §iÒu tra Quy Ho¹ch Rõng quy ®Þnh hoÆc do ngêi gi¶i ®o¸n cô thÓ quy ®Þnh. C¸ch lµm nh sau: Interpreter GIS analysis Record: ®a tªn file cÇn record Out put: ®Æt tªn file sau khi record OK 0.1.11¶nh sau ph©n lo¹i kiÓm tra, gép nhãm Sau khi ®· läc ¶nh råi nhng mét sè l« vÉn cßn nhá, cßn khã sö dông trong thùc tÕ. §Ó phï hîp víi nh÷ng quy ®Þnh trong viÖc x©y dùng b¶n ®å ta cÇn gép nh÷ng l« nhá vµo víi nh÷ng l« lín h¬n ë bªn c¹nh. ViÖn §iÒu tra Quy ho¹ch Rõng quy ®Þnh l« nhá nhÊt trªn b¶n ®å cã kÝch thíc lµ 2x2mm, t¬ng ®¬ng víi 4 ha trªn b¶n ®å tû lÖ 1/100.000. Trªn Menu chän: Eliminated 18
Input file: tªn file ®· Clump Out put File: ®Æt tªn file sÏ lîc bít Minimum: 45 pixel OK. Chó ý: sau khi Recod-Clump- Eliminate cÇn ®æi l¹i mµu c¸c ®èi tîng nh cò. C¸ch lµm nh sau: Trªn cöa sæ View chän: View Raster Atribute
KÝch chuét vµo ®Çu mçi dßng ®Ó ®æi mµu cho nhãm. CÇn ghi l¹i m· cña tõng nhãm ë phÇn Recod ®Ó ®æi mµu cho ®óng. 0.1.12 ChuyÓn ®æi format ®Ó chång xÕp, xö lý b¶n ®å ChuyÓn ¶nh sè tõ format *.IMG sang *.Grid Menu chÝnh chän: EXPORT Type: Grid In put file: ®a tªn file ®· eliminate Out put name (®Æt tªn file, th môc ) OK Export option: File OK (®îc file d¹ng *.grid)
19
0.1.13 ChuyÓn raster file.grid sang vecter file.shp Khëi ®éng Acr/view File Extension Spatial analysic OK New view Add theme Tªn file th môc chøa file *.GRID: Data source type: grid data source T×m file *.grid cÇn chuyÓn Theme Convert to shape file -> (chän tªn, th môc ®Æt file kÕt qu¶ d¹ng *.Shp) OK Add shape file as theme to view Yes 20
0.1.14ChuyÓn file.shp sang file.tab ®Ó xö lý b»ng Map/Infor Khëi ®éng Map/Infor Tool Universal translator From format: *shp Data source: ®a tªn file cÇn chuyÓn To format: *tab Destination: ®a ®êng dÉn vµo th môc cÇn thiÕt 0.2 Chång xÕp b¶n ®å hiÖn tr¹ng cña hai thêi kú b»ng phÇn mÒm ARC VIEW: §Ó cã ®îc mét b¶n ®å diÔn biÕn hiÖn tr¹ng tµi nguyªn rõng vµ sö dông ®Êt ®ai gi÷a hai thêi kú, ta sÏ ph¶i x©y dùng b¶n ®å hiÖn tr¹ng riªng cho tõng thêi kú vµ sau ®ã chång xÕp hai b¶n ®å víi nhau. Sau khi cã ®îc b¶n ®å hiÖn tr¹ng cña hai thêi kú b»ng ph©n lo¹i ¶nh sè, ta chuyÓn sang file d¹ng *.shp Ta tiÕn hµnh chång xÕp hai b¶n ®å b»ng phÇn mÒm ARC VIEW nh sau:
Khëi ®éng phÇn mÒm ARC VIEW
Untitled
Views
Add theme ( HoÆc kÝch vµo dÊu + ) ®Ó më hai file b¶n ®å hiÖn tr¹ng d¹ng *.shp
21
Sau khi më hai file b¶n ®å ta cã:
22
23
24
X¸c ®Þnh ®¬n vÞ lµm viÖc cho b¶n ®å:
View
Properties
Map units: meters
Distance units: meters
TiÕp theo vµo:
File
Extensions
Geoprocessing
View
GeoProcessing Wizard
Intersect two theme
25
Next
Chän c¸c file input vµ output
Finish
Sau khi ch¹y xong ta ®îc file b¶n ®å chång xÕp cña hai b¶n ®å víi hai thêi ®iÓm kh¸c nhau.
26
TiÕp theo ta tÝnh l¹i diÖn tÝch cho b¶n ®å biÕn ®éng hai thêi kú :
Më file biÕn ®éng
Theme
Table
Table - Start Editing
Edit - Add Field §a vµo tªn trêng (DT_98_02), kiÓu trêng, ®é réng …
Field - Calculate {Shape}.ReturnArea/10000 ha)
( TÝnh cho ®¬n vÞ
27
TiÕp theo ta x¸c ®Þnh trêng tr¹ng th¸i míi cho b¶n ®å biÕn ®éng (thùc chÊt lµ x¸c ®Þnh xem tr¹ng th¸i l« biÕn ®æi tõ thêi kú I sang thêi kú II)
Edit - Add Field
Table - Query (Chän lÇn lît c¸c m· tr¹ng th¸i cña thêi kú I biÕn ®æi thµnh c¸c m· tr¹ng th¸i cña thêi kú II)
New Set
Calculate (§a m· tr¹ng th¸i míi t¬ng øng vµo)
28
KÕt thóc tÝnh to¸n b»ng:
Table - Stop Editing
Save edits - Yes
29
Ta ®îc b¶n ®å biÕn ®éng hiÖn tr¹ng sö dông ®Êt vµ tµi nguyªn rõng nh sau: Sau ®ã chuyÓn sang biªn tËp b¶n ®å biÕn ®éng b»ng phÇn mÒm
Map/Info vµ tÝnh to¸n b¶ng biÓu thèng kª b»ng Excel
1
PhÇn II: Ph©n tÝch, thèng kª c¸c lo¹i ®Êt ®ai, tµi nguyªn rõng theo ®ai cao, ®é dèc
1.1 X©y dùng m« h×nh sè ho¸ ®Þa h×nh DTM: §Ó x©y dùng c¸c b¶n ®å ®ai cao, ®é dèc tríc hÕt ta ph¶i tiÕn hµnh x©y dùng m« h×nh sè ho¸ ®Þa h×nh DTM ( Digital Terrain Model) Do ARC VIEW lµm viÖc trªn trêng ID, nªn ta ph¶i t¹o l¹i gi¸ trÞ cña trêng nµy b»ng chÝnh c¸c gi¸ trÞ ®é cao
Më ARC View
Më file ®Þa h×nh d¹ng *.shp
View
Properties 30
Map Unit: m
Distance Unit: m
T¹o trêng id víi c¸c gi¸ trÞ ®é cao:
Table
Start Editing
Add field
Name field: id
Calculate (trêng id = trêng value)
Stop editing
Save
Yes
1.1.1.1TiÕn hµnh néi suy DTM A, X©y dùng DTM b»ng thuËt to¸n TIN - Khëi ®éng extention 3D analysis:
File 31
-
-
Extentions
Chän 3D analysis
OK
Më File ®Þa h×nh:
Nh¸y chuét vµo biÎu tîng “+”
T×m ®êng dÉn ®Õn file ®Þa h×nh
TiÕn hµnh néi suy DTM
Surface Create TIN from Features
Height source: chän “id”
OK
X©y dùng b¶n ®å ®ai cao:
KÝch ho¹t file TIN
Them
Convert to Grid
§Æt tªn file
OK
X©y dùng b¶n ®å ®é dèc:
KÝch ho¹t file TIN
Surface
Derive Slope
Them
Convert to Grid 32
§Æt tªn file
OK
b, X©y dùng DTM b»ng thuËt to¸n Spline T¹o file sè ho¸ ®Þa h×nh b»ng c¸c ®iÓm cã gi¸ trÞ ®é cao:
File
Extention
3D
Scrip
Load text file Pl2intpt1.ave (C:\esri\av_gis30\arcview\samples\scrips\pl2int~1.ave)
Compile
Run
Enter point seperation distance in map unit: 20 ( kho¶ng c¸ch nµy tuú thuéc vµo sè lîng ®êng ®ång møc dµy hay tha)
Enter input field name for polyline elevation or id …: id 33
Enter output field name for point indentification: value
§Æt tªn file ra
Sau khi ch¹y xong ta ®îc mét file sè ho¸ ®Þa h×nh gåm toµn c¸c ®iÓm cã gi¸ trÞ ®é cao vµ tiÕn hµnh néi suy trªn file nµy.
§a file ranh giíi cña khu vùc vµo ®Ó giíi h¹n
KÝch ho¹t file cÇn néi suy
Surface
Interpolate Grid …
Output grid extent: Tªn file ranh giíi
Output grid cell size: 10
OK
Method : Spline
Z value field: id
Weight: 50
No of point: 50
Type: Tension
OK
34
Sau khi ch¹y néi suy xong, ph¶i ghi vµo ®Üa file DTM:
Theme
Save Data Set…
C¾t file DTM theo file ranh giíi:
Më file polygon ranh giíi d¹ng *.shp
KÝch ho¹t file ranh giíi
Theme
Start Editing Dïng mòi tªn chän ranh giíi
Edit
Coppy feature 35
Theme
Stop Editing
No Save
Edit
Past KÝch ho¹t file grid DTM
Transformation
Clip Grid
Inside
Theme
Save Data set
Sau khi kÕt thóc phÇn néi suy ta sÏ ®îc file DTM vµ cã thÓ kiÓm tra ®îc ®Þa h×nh, ®Þa thÕ cña khu vùc b»ng ph¬ng ph¸p Bãng nói (Hillshade) hoÆc h×nh häc kh«ng gian 3 chiÒu (3D).
36
Nh×n theo Hillshade:
37
38
39
Nh×n theo 3D:
40
X©y dùng b¶n ®å ®ai cao:
Tõ m« h×nh sè ho¸ ®Þa h×nh DTM ta tiÕn hµnh x©y dùng b¶n ®å ®ai cao nh sau:
Më file DTM
KÝch ho¹t file DTM
§æi gi¸ trÞ ®é cao tõ sè thËp ph©n sang sè nguyªn
Analysis
Map Calculator (File DTM).int
Evaluate
Theme
Save Data set…
Dùa theo c¸c tiªu chuÈn ®ai cao ®Ó x©y dùng b¶n ®å ®ai cao:
KÝch ho¹t file DTM 41
Analysis
Reclassify
Classification field: value
Classify Type: Equal Interval Number of classes: 5 Round value at: d Thay ®æi c¸c gi¸ trÞ ®ai cao vµo Old values
Theme
Save Data set…
Ta ®îc b¶n ®å ®ai cao nh sau:
42
1.2 Chång xÕp b¶n ®å ®ai cao vµ b¶n ®å hiÖn tr¹ng: Sau khi ®· x©y dùng ®îc b¶n ®å ®ai cao vµ b¶n ®å hiÖn tr¹ng ta sÏ thùc hiÖn chång xÕp hai b¶n ®å víi nhau theo c¸c bíc nh ë phÇn chång xÕp b¶n ®å biÕn ®éng gi÷a hai thêi kú Më hai file b¶n ®å hiÖn tr¹ng vµ b¶n ®å ®ai cao:
X¸c ®Þnh ®¬n vÞ lµm viÖc cho b¶n ®å:
View
Properties
Map units: meters
Distance units: meters
TiÕp theo vµo:
File
Extensions
Geoprocessing
View 43
GeoProcessing Wizard
Intersect two theme
Next
Chän c¸c file input vµ output
Finish
Sau khi ch¹y xong ta ®îc file b¶n ®å chång xÕp hiÖn tr¹ng - ®ai cao: 44
TiÕp theo ta tÝnh l¹i diÖn tÝch cho b¶n ®å hiÖn tr¹ng - ®ai cao:
Më file b¶n ®å
Theme
Table
Table - Start Editing
Edit - Add Field §a vµo tªn trêng (DT_dc_r), kiÓu trêng, ®é réng …
Field - Calculate {Shape}.ReturnArea/10000
(TÝnh cho ®¬n vÞ ha)
KÕt qu¶ ta còng ®îc mét b¶n ®å hiÖn tr¹ng ph©n bè theo ®ai cao vµ b¶ng thèng kª diÖn tÝch c¸c lo¹i ®Êt ®ai, lo¹i rõng theo ®ai cao 1.3 X©y dùng b¶n ®å ®é dèc: Ta x©y dùng b¶n ®å ®é dèc tõ file m« h×nh sè ho¸ ®Þa h×nh DTM:
45
Më file DTM
KÝch ho¹t file DTM
Surface
Derive Slope
Theme
Save Data set… (File ®é dèc d¹ng sè thËp ph©n)
Sau ®ã ®æi gi¸ trÞ ®é dèc tõ sè thËp ph©n vÒ sè nguyªn:
Më file ®é dèc
KÝch ho¹t file ®é dèc
Analysis
Map Calculator (File ®é dèc).int
Evaluate
Theme
Save Data set… 46
Dùa theo c¸c tiªu chuÈn ®é dèc ®Ó ph©n lo¹i l¹i b¶n ®å c¸c cÊp dé dèc:
KÝch ho¹t file ®é dèc
Analysis
Reclassify
Classification field: value
Classify Type : Equal Interval Number of classes: 4 Round value at: d Thay ®æi c¸c gi¸ trÞ ®é dèc vµo Old values
Theme
Save Data set…
47
Ta ®îc b¶n ®å cÊp ®é dèc:
1.4 Chång xÕp b¶n ®å ®é dèc vµ b¶n ®å hiÖn tr¹ng: Sau khi ®· x©y dùng ®îc b¶n ®å cÊp ®é dèc vµ b¶n ®å hiÖn tr¹ng ta sÏ thùc hiÖn chång xÕp hai b¶n ®å víi nhau qua c¸c bíc nh ®· lµm víi b¶n ®å hiÖn tr¹ng - ®ai cao KÕt qu¶ ta còng ®îc mét b¶n ®å hiÖn tr¹ng ph©n bè theo cÊp ®é dèc vµ b¶ng thèng kª diÖn tÝch t¬ng øng. 2
2.1
PhÇn III: øng dông c«ng nghÖ viÔn th¸m vµ GIS vµo viÖc x¸c ®Þnh khu vùc u tiªn trång rõng
X¸c ®Þnh c¸c chØ tiªu tham gia vµo viÖc x¸c ®Þnh vïng u tiªn cho trång rõng:
ViÖc x¸c ®Þnh c¸c chØ tiªu tham gia vµo viÖc x¸c ®Þnh vïng u tiªn cho trång rõng mang ý nghÜa rÊt quan träng. §©y lµ bíc ®ßi hái tÝnh chuyªn gia rÊt cao, ngoµi ra cßn phô thuéc vµo ®Æc tÝnh cña tõng khu 48
vùc. Díi sù hç trî cña phÇn mÒm Geo-Choice, viÖc lùa chän c¸c chØ tiªu nµy sÏ mang tÝnh kh¸ch quan trªn c¬ së kÕt hîp nhiÒu ý kiÕn cña c¸c chuyªn gia. KÕt qu¶ nµy sÏ h¹n chÕ yÕu ®iÓm cña c¸c ph¬ng ph¸p lùa chän tríc ®ay lµ chØ dùa vµo ý kiÕn chñ quan cña 1 c¸ nh©n ®Ó ®a ra kÕt qu¶ lùa chän.
C¸c bíc tiÕn hµnh: - S¬ bé x¸c ®Þnh c¸c nh©n tè tham gia - §a c¸c nh©n tè ®ã vµo Table trong ARC/VIEW
-
Më file Geo_choice
Më Table: nh¸y chuét vµo biÓu tîng h×nh d¹ng nh 1 b¶ng
Table
Start Editting
Edit
Add Field
§¸nh tªn nh©n tè tham gia
OK
TiÕp tôc cho hÕt c¸c nh©n tè.
File
Export
Chän Delimited Text
OK
§Æt tªn ra
OK
Chän c¸c nh©n tè tham gia:
Khëi ®éng Choice explorer 1.2
File
Import Delimited Text File
Chän file ®· export tríc ®ã
OK
Chän nh÷ng nh©n tè theo ý kiªn riªng råi Ên “
“
49
OK
OK
Vote
Criteria Selection
OK
§Æt tªn file
OK
- Khëi ®éng Choice Perspectivies 1.2
File
Open
Chän tÊt c¶ c¸c file
OK
View
Multiple files
X¸c ®Þnh c¸c nh©n tè ®îc c¸c c¸ nh©n chän nhiÒu.
2.2 X¸c ®Þnh träng sè cho c¸c nh©n tè tham gia: Do møc ®é ¶nh hëng cña c¸c nh©n tè tham gia vµo viÔc x¸c ®Þnh khu vùc u tiªn trång rõng sÏ kh¸c nhau. V× vËy, cÇn thiÕt ph¶i x¸c ®Þnh träng sè (møc ®é quan träng) cho tõng nh©n tè. Trong kho¸ häc nµy, träng sè cho tõng nh©n tè sÏ ®îc x¸c ®Þnh b»ng ph¬ng ph¸p “So s¸nh tõng cÆp” – Pairwise Comparision díi sù hç trî cña phÇn mÒm Choice explorer.
Khëi ®éng phÇn mÒm Choice explorer 50
File
Import Delimited Text File
Chän file chØ cã c¸c nh©n tè ®îc lùa chän ®· export tríc ®ã
OK
Chän tÊt c¶ nh÷ng nh©n tè råi Ên “
X¸c ®Þnh kh¶ n¨ng ¶nh hëng cña c¸c nh©n tè
OK
Tool
Asign Criterion Weight
Pairwise Comparision
“
51
TiÕn hµnh so s¸nh tõng cÆp.
Apply
File
Save: ®Æt tªn
Vote
Criterion Weight
§Æt tªn
- Tæng hîp c¸c kÕt qu¶: Khëi ®éng Choice Perspectivies 1.2
File
Open
Chän tÊt c¶ c¸c file
OK
View
52
Multiple files
X¸c ®Þnh träng sè b×nh qu©n cho tõng nh©n tè.
2.3 Cho ®iÓm cho tõng b¶n ®å ViÖc cho ®iÓm cho tõng b¶n ®å chuyªn ®Ò sÏ phô thuéc vµo sù ¶nh hëng tÝch cùc hay tiªu cùc cña c¸c nh©n tè ®ã ®Õn viÖc x¸c ®Þnh møc ®é u tiªn cho trång rõng §èi víi c¸c nh©n tè cã ¶nh hëng tiªu cùc nh §ai cao: Cµng lªn cao (gi¸ trÞ ®é cao cµng lín) th× kh¶ n¨ng trång rõng cµng gi¶m, viÖc tÝnh to¸n cho ®iÓm cho nh©n tè nµy sÏ ®îc tÝnh nh sau: Gi¸ tri míi =1 – [Vt – Vmin)/(Vmax- Vmin)] §èi víi c¸c nh©n tè cã ¶nh hëng tÝch cùc nh ®é dµy tÇng ®Êt: tÇng ®Êt cµng dÇy th× kh¶ n¨ng trång rõng cµng tèt, viÖc tÝnh to¸n cho ®iÓm cho nh©n tè nµy sÏ ®îc tÝnh nh sau: Gi¸ tri míi = Vt – Vmin)/(Vmax- Vmin) Trong ®ã:
Vt: gi¸ trÞ thùc Vmin: gi¸ trÞ nhá nhÊt Vmax: gi¸ trÞ lín nhÊt
2.3.1.1Cho ®iÓm b¶n ®å ®é cao:
Më b¶n ®å ®é cao
Analysis
Map Calculator
§¸nh c«ng thøc: ( ([Docao] - 180)/(1120-180))
Evaluate
Më map 1
Analysis
Map Calculator
§¸nh c«ng thøc: ( [1] - [Map Calculation 1])
Evaluate
Them
Convert to Grid
§Æt tªn file
OK 53
2.3.1.2Cho ®iÓm b¶n ®å ®é dèc
Më b¶n ®å ®é dèc
Analysis
Map Calculator
§¸nh c«ng thøc: ( [Dodoc] / 89.403)
Evaluate
Më map 1
Analysis
Map Calculator
§¸nh c«ng thøc: ( [1] - [Map Calculation 1])
Evaluate
Them
Convert to Grid
§Æt tªn file
OK
2.3.1.3T¹o b¶n ®å kho¶ng c¸ch ®Õn ®êng:
Më File ®êng
Analysis
Find Distance
Output cell Size: 15
Output Grid Extend: Same as Road_polyline
OK
Analysis
Map Calculator
§¸nh c«ng thøc: ( [Distance to Road_polyline.shp] / 11535.966)
Evaluate
Mowr File 1
Analysis
Map Calculator
§¸nh c«ng thøc: ( [1] - [Map Calculation 1]) 54
Evaluate
Them
Convert to Grid
§Æt tªn file
- T¹o b¶n ®å kho¶ng c¸ch ®Õn khu d©n c vµ s«ng suèi lµm t¬ng tù nh trªn. 2.3.1.4T¹o b¶n ®å §Êt: ViÖc x¸c ®Þnh møc ®é u tiªn cho trång rõng ®îc chia thµnh 3 møc: cÊp cao, cÊp trung b×nh, cÊp thÊp. Nh©n tè §Êt dîc ph©n cÊp thµnh 5 cÊp: + CÊp 1: RÊt tèt, hÇu nh kh«ng cã h¹n chÕ trong sö dông (T) + CÊp 2: Tèt, cã Ýt h¹n chÕ trong sö dông + CÊp 3: Trung b×nh + CÊp 4: XÊu, nhiÒu h¹n chÕ trong sö dông + CÊp 5: RÊt xÊu, trªn nói ®¸ v«i §iÓm sè cho c¸c chØ tiªu t¬ng øng víi c¸c cÊp 1,2,3,4,5 lµ ®iÓm 5,4,3,2,1 CÊp ®Êt
M« t¶
Ký hiÖu
§iÓm sè
1
RÊt tèt, hÇu nh kh«ng cã h¹n chÕ sö T dông
5
2
Tèt, cã Ýt h¹n chÕ trong sö dông
Fk
4
3
Trung b×nh
Fa, Fc
3
4
XÊu, nhiÒu h¹n chÕ trong sö dông
Fs
2
5
RÊt xÊu, trªn nói ®¸ v«i
K
1
Më file GRID "lo¹i ®Êt"
Më biÓu
Thªm trêng "Diem"
Cho gi¸ trÞ c¸c lo¹i ®Êt theo b¶ng trªn 55
Thªm trêng "Dat_val"
Nh¸y vµo biÓu tîng c¸i m¸y tÝnh (Calculate)
§¸nh c«ng thøc: [Diem] / 5
OK
2.3.1.5 T¹o b¶n ®å §é dÇy tÇng ®Êt: Ta sÏ cho ®iÓm trªn c¸c chØ tiªu cña c¸c tiªu chÝ theo b¶ng sau:
CÊp ®é dÇy
M« t¶
§iÓm sè
1
TÇng ®Êt dÇy trªn 100 cm
3
2
TÇng ®Êt trung b×nh tõ 50 - 100 cm
2
3
TÇng ®Êt máng díi 50 cm
1
Më File GRID : tangdat
Më biÓu
Thªm trêng "Diem"
Cho ®iÓm nh b¶ng trªn
Thªm trêng "DD_val"
Nh¸y biÓu tîng c¸i m¸y tÝnh (Calculate)
§¸nh c«ng thøc: [Diem] / 3
OK
Save
56
2.4 Chång xÕp c¸c líp b¶n ®å: Sau khi ta ®· cã ®Çy ®ñ 6 lo¹i b¶n ®å cña 6 tiªu chÝ cÇn thiÕt ®Ó x¸c ®Þnh vïng u tiªn trång rõng, ta sÏ tiÕn hµnh céng ®iÓm nh sau:
Më lÇn lît 6 file b¶n ®å
Analysis
Map calculater (Ta lµm phÐp céng cña 6 trêng ®iÓm sè cña 6 file b¶n ®å):
Evaluate
Theme
Save Data set... (ta ®îc file b¶n ®å céng ®iÓm cña 6 file d÷ liÖu – chÝnh lµ file b¶n ®å tæng hîp)
2.5 Ph©n vïng u tiªn trång rõng: Khu vùc trång rõng ®îc ph©n thµnh 3 møc ®é u tiªn nh sau: cao, trung b×nh, thÊp ViÖc x¸c ®Þnh ®iÓm cña 3 vïng u tiªn ®îc x¸c ®Þnh nh sau: Vtb = (Vmax - Vmin)/3 - Møc u tiªn thÊp cã gi¸ trÞ tõ
Vmin ®Õn Vmin + Vtb
- Møc u tiªn trung b×nh cã gi¸ trÞ tõ 2Vtb
Vmin + Vtb ®Õn Vmin +
- Møc u tiªn cao cã gi¸ trÞ tõ
Vmin + 2Vtb ®Õn Vmax
Trong ®ã: Vtb:
Gi¸ trÞ trung b×nh cña 1 cÊp
Vmin: Gi¸ trÞ nhá nhÊt trong b¶n ®å tæng hîp Vmax: Gi¸ trÞ lín nhÊt trong b¶n ®å tæng hîp Tõ c¬ së ph©n c¸c møc ®é u tiªn nh vËy ta sÏ tiÕn hµnh x©y dùng b¶n ®å ph©n vïng u tiªn trång rõng nh sau:
Më file tæng hîp
KÝch ho¹t file tæng hîp
Analysis
Reclassify (Ph©n l¹i thµnh 3 møc u tiªn) 57
Classification field: value
Classify Type: Equal Interval Number of classes: 3 Round value at: d Thay ®æi c¸c gi¸ trÞ cÊp tiÒm n¨ng ®Êt vµo Old values
Theme
Save Data set
Ta ®îc b¶n ®å ph©n vïng u tiªn trång rõng theo 3 cÊp
58
59