Xử lý và phân tích ảnh viễn thám với ENVI tài liệu hướng dẫn thực hành phần II

C¤NG TY TNHH T¦ VÊN GEOVIÖT X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h 6/17, Ngõ 139 Nguy n Ng c V - C u Gi y Hà n i – VI T NAM ng d n th c hành M CL C Tel./Fax: +84.4.5564 710 Email: geoviet@gmail.com Web: http://www.geoviet.vn II N N CH NH VÀ X LÝ NH V I ENVI 22 2.1 T ng quan v module 22 2.2 Các files s d ng 22 2.3 Các h to đ thông d ng nh đ c n n ch nh ENVI 22 TÀI LI U H NG D N TH C HÀNH X LÝ VÀ PHÂN TÍCH D LI U VI N THÁM V I PH N M M ENVI Biên t p: TS Tr n Hùng KS Ph m Quang L i 2.4 N n ch nh hình h c nh v i nh 23 Bài t p N n ch nh nh SPOT v i nh ASTER (n n) vùng Nghĩa Đô HN 23 M hi n th nh SPOT & ASTER 23 Kh i đ ng n n nh đ c m kh ng ch (GCP) 23 Làm vi c v i GCP 25 Th c hi n n n nh 25 So sánh k t qu n n 27 Ki m tra t a đ b n đ 27 2.5 N n ch nh hình h c nh v i b n đ n n 27 Bài t p N n ch nh hình h c nh SPOT v i b n đ n n vùng Nghĩa Đô 27 M hi n th nh 27 M file d li u Vector GIS shape file format 27 M t s b c quy trình n n nh v i b n đ n n 29 2.6 C t nh 32 2.6.1 C t nh theo m t khu v c nghiên c u 32 2.6.2 Hi n th khu v c Qu n C u Gi y theo đ n v hành 34 2.7 Ghép nh (mosaicking) 37 2.8 K thu t tăng đ nét l c nh (convolutions) 38 2.8.1 Tăng đ nét 38 2.8.2 L c nh (Filter) 39 2.9 K thu 2.9.1 2.9.2 2.9.3 t u ch nh bi n đ i nh 39 T l c a kênh nh (band ratios) 39 T o kênh ch s th c v t (NDVI) 40 Các phép tính kênh (Band Math) 42 Hà N i, 2008 20 @2008 T v n GeoVi t – www.geoviet.vn 21 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h II N N CH NH VÀ X ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h lên khung đ c n n) h tr ph ng pháp nearest neighbor (pixel g n nh t), n i suy n tính b c n i suy b c (cubic convolution) So sánh gi a nh g c (n n) nh n n s d ng kh Dynamic Overlay c a ENVI cho phép đánh giá đ xác c a phép n n m t cách nhanh chóng LÝ NH V I ENVI 2.1 T ng quan v module Module th c hành g m ph n: n n ch nh x lý / bi n đ i nh Ph n đ u cung c p nh ng thông tin c s v n n ch nh hình h c nh v i ENVI (n n ch nh nh v i nh n n ch nh nh v i b n đ n n) Quy trình n n ch nh t ng b c s đ c trình bày v i gi thuy t b n có s hi u bi t nh t đ nh v n n ch nh qua gi ng lý thuy t Ph n sau cung c p m t s ch c x lý / bi n đ i nh c s đ giúp cho vi c hi n th gi i đoán nh 2.2 Các files s 2.4 N n ch nh hình h c nh v i nh Ph n c a module s d n d t b n qua t ng b c c a quy trình n n ch nh nh v i nh nh nghiado_aster123geo đ c dùng làm nh chu n đ n n ch nh nh nghiado_spot95 cho v khung hình h c Bài t p N n ch nh nh SPOT v i nh ASTER (n n) vùng Ngh a ô HN d ng nh dùng đ th c hi n th c hành module g m nh ng nh đ c cung c p th m c \RS_Training_Dataset\DATA nh gi i thi u module Đ làm vi c theo nhóm mang tính ch t g n v i nh ng ng d ng th c ti n, m t s nh v tinh ph m t ph n Hà n i, Kiên giang & Cà mau đ c ch n dùng 2.3 Các h to đ thông d ng nh đ ng d n th c hành - c n n ch nh ENVI ENVI h tr vi c đ c x lý nh đ c n n ch nh v nhi u h chi u to đ thông d ng nh gi i thi u Ngoài ENVI có th h tr nh ng h t a đ ng i dùng t l p b ng cách đ a vào thay th nh ng thông s đ c tr nhóm h quy chi u, 35 d ng ellipsoids khác h n 100 datums đ đáp ng h u h t nh ng yêu c u v đo v b n đ M hi n th - nh SPOT & ASTER Th c hi n b c đ c gi i thi u module “D li u vi n thám ph n m m ENVI” đ m hi n th nh nghiado_spot95 hi n th (Display) #1 & nghiado_aster123geo t i hi n th #2; Hi n th v trí / giá tr tr đ quan sát to đ c a pixel nh SPOT & nh ASTER; Kh i đ ng n n nh đ c m kh ng ch (GCP) a) Ch n Map\ Registration\ Select GCPs: Image to Image b) Khi h p tho i Image to Image Registration xu t hi n, kích vào Display #2 ( nh ASTER) ch n Base Image ( nh chu n) Display #1 ( nh SPOT) ch n Warp Image ( nh c n n n); Nh ng thông s h to đ ENVI đ c l u gi file ASCII “map_proj.txt”, có th d dàng thay đ i / b xung b ng nh ng công c ENVI Map Projection (h quy chi u) ho c có th đ c s a tr c ti p (b ng ph n m m so n th o đ n gi n) Nh ng thông tin file s đ c s d ng file header g n v i m i d li u nh cho phép s liên k t v trí “Magic Pixel” v i nh ng to đ c a h quy chi u xác đ nh M t s ch c c a ENVI có th dùng thông tin đ làm vi c v i d li u nh không gian c a d li u có t a đ Nh ng công c n n ch nh hình h c c a ENVI cho phép b n n n nh ng nh v h t a đ đ a lý s a đ i đ thích h p v i hình h c nh c s Nh ng m kh ng ch (GCP) đ c ch n c a s v i đ đ phân gi i không gian (c a s Chính) c a s Zoom đ i v i c ph ng th c n n nh v i nh n n nh v i b n đ n n To đ đ c hi n th cho m GCP c nh g c nh đ c n n ch nh v i sai s cho t ng thu t toán n n ch nh Vi c n n đ c th c hi n b ng ph ng th c RST (resampling, scaling & translation), n n đa th c (b c đ n n) ho c n n tam giác Delaunay Vi c resampling (tính toán ghi d li u t ng pixel Hình 2.1: C a s l a ch n nh n n nh g c Kích vào “OK” đ b t đ n trình n n nh 22 @2008 T v n GeoVi t – www.geoviet.vn 23 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành c) H p tho i Ground Control Points Selection s xu t hi n t ng m kh ng ch (GCP) s đ c nh p vào b ng cách đ t v trí tr c nh vào v trí m t đ t (giao m c a đ ng, m đ c bi t ) (Hình 2.2) d) Di chuy n tr nh ASTER đ n 201, 200 b ng cách nh p nh ng giá tr vào h p tho i Ground Control Point Selection ô Base X & Y; Di chuy n tr nh SPOT đ n 1261.84, 1834.58 b ng cách nh p nh ng giá tr t ng t vào h p tho i ô Warp X & Y; e) Ki m tra v trí tr c a s Zoom u ch nh v trí n u c n b ng cách kích phím trái chu t m i c a Zoom đ n v trí yêu c u; Chú ý: c a s Zoom có th h tr cho vi c tinh ch nh v trí đ n giá tr hàng ph n trăm c a pixel (sub pixel) Kh tinh ch nh tùy thu c vào nhân t khuy ch đ i c a c a s Zoom f) Kích vào “Add” đ có th nh p m GCP vào danh sách Kích “Show List” đ xem danh sách nh ng m GCP ch n Th làm v i vài m GCP đ có đ c k ch n GCP g) Chú ý vào danh sách c a nh ng m th c m d tính (predicted), s m ch n l n h n 4, sai s RMS s đ c tính hi n th h) Ch n Options\ Clear All Points h p tho i Registration GCP List đ xoá nh ng m b n ch n; i) Ch n File\ Restore GCPs from ASCII h p tho i Ground Control Point Selection kích vào tên file nghiado_spot95_toAster123geo.pts j) Kích vào “OK” đ đ c danh sách nh ng m GCP l u tr c; k) Kích vào t ng m GCP h p tho i Image to Image GCP List ki m tra v trí c a m nh, nh ng t a đ th c t t a đ d tính, sai s RMS Quan sát sai s RMS t ng h p tho i Ground Control Point Selection Hình 2.2: H p tho i Ground Control Point Selection h p tho i Image to Image GCP List n n ch nh hình h c nh v i nh Làm vi c v i GCP Nh ng miêu t d i có th đ c dùng nh thông tin tham kh o: - V trí c a m i m GCP có th đ c s a b ng cách ch n chúng danh sách h p tho i GCP List s a chúng h p tho i GCP Selection Ho c b n có th nh p vào m t v trí m i c a pixel ho c b n có th thay đ i v trí theo t ng pixel m t b ng cách nh n vào nút có hình mũi tên c nh giá tr X, Y h p tho i; - Kích vào nút “On/Off” h p tho i GCP List đ lo i b m GCP đ c ch n kh i ph ng trình tính toán h s n n công th c tính sai s Nh ng m GCP ch b b qua tính toán ch ch a b xoá kh i danh sách b n có th s d ng l i chúng kích vào nút “On/Off” l n n a - Kích vào nút “Delete” đ xoá m GCP kh i danh sách - Đ t v trí tr c a s zoom kích vào nút “Update” đ thay đ i m GCP đ c ch n v v trí m i c a tr - Nút “Predict” cho phép tính toán nh ng m GCP m i d a mô hình n n dùng a) Th đ t tr vào m t v trí m i nh ASTER kích vào nút “Predict”, v trí tr nh SPOT s t đ ng chuy n đ n v trí t ng ng v i v trí d tính theo mô hình n n b) V trí xác sau có th đ c tinh ch nh b ng cách di chuy n m t chút v trí pixel nh SPOT c) Kích vào nút “Add Point” đ nh p thêm m GCP cho danh sách Th c hi n n n nh a b c d e nh có th đ c n n ch cho nh ng kênh hi n th ho c cho t t c kênh c a nh lúc, ví d c a vi c ch n n nh ng kênh hi n th : Ch n Options\ Warp File… h p tho i GCP Selection H p tho i Input Wrap Image hi n ch n nh c n n n (nghiado_Spot95) r i nh p OK Khi h p tho i Registration Parameters xu t hi n, ch n “RST” cho Warp Method ph ng th c resampling “Nearest Neighbor” Ch n tên file nh n n nghiado_SPOTtoASTER123geo b ng cách tích vào File nh p vào choose đ l a ch n th m c ch a file d li u kích “OK” nh đ c n n s đ c hi n th Available Bands List sau trình n n ch nh k t thúc Th l p l i v i ph ng pháp resampling Bilinear Interpolation Cubic Convolution Đ t tên file k t qu nghiado_SPOTtoASTER123geo_b & nghiado_SPOTtoASTER123geo_c 24 @2008 T v n GeoVi t – www.geoviet.vn ng d n th c hành 25 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành f Th c hi n vi c n n đa th c b c n n Delaunay Triangulation v i Cubic Convolution g Đ t tên file k t qu nghiado_SPOTtoASTER123geo_pc & nghiado_SPOTtoASTER123geo_dtc ng d n th c hành So sánh k t qu n n S d ng Dynamic Overlay đ so sánh k t qu : a) Đóng file nghiado_SPOT b ng cách ch n kích vào tên file ch n File\ Close Selected File c a s Available Bands List b) Đ c file nh n n l n l t hi n th lên Display m i c) Đ t tr vào c a s Chính kích vào phím ph i chu t sau ch n File\ Link\ Link Displays Kích vào “OK đ liên k t nh ASTER nh SPOT n n d) So sánh nh SPOT & ASTER v i Dynamic Overlay b ng cách nháy nút trái chu t e) Quan sát hi u qu c a ph ng pháp n n ph ng th c resampling Ki m tra t a đ b n đ Làm hi n lên h p tho i Cursor Location/ Value di chuy n chu t c a s quan sát hi u ng c a ph ng phán n n khác lên giá tr c a pixel 2.5 N n ch nh hình h c nh v i b n đ n n Quy trình n n ch nh hình h c nh v i b n đ n n (chu n) t ng t nh quy trình n n nh v i nh T a đ b n đ h t a đ UTM Zone 48, datum WGS 84 d a d li u GIS shape file format v i đ ng giao thông (roadpoly.shp & roadlines.shp) & sông su i (hydro.shp) c a vùng nghiên c u (đã đ c cung c p th m c làm vi c) s đ c s d ng làm Base (d li u chu n) to đ pixel nh nghiado_Spot95 s đ c n n ch nh v t a đ b n đ To đ c a t ng m GCP b n đ s đ c đ c sau hi n th nh ng file ph n m m ArcView Hình 2.3: H p tho i Registration Parameters Các ph ng pháp n n nh: Œ Ph ng pháp RST – Rotating, Scaling, Translation: ch th c hi n nh ng chuy n d ch đ n gi n: xoay, xác đ nh t l t nh ti n nh Œ Ph ng pháp Polynomial – Hàm đa th c: ph ng pháp cho k t qu t t h n ph ng pháp RST, v i yêu c u s v s m kh ng ch N t ng ng v i b c c a hàm n nh sau: N > (n+1)2 Œ Ph ng pháp Triangulation – l i tam giác: ENVI s d ng nguyên lý tam giác Delaunay đ n n nh b ng cách ch n m kh ng ch làm đ nh c a tam giác không đ u ti n hành n i suy Đ ti n hành n n nh ta ph i l a ch n m t ba ph ng pháp tái chia m u – Resampling cho đ t đ c k t qu mong mu n (Hình 2.8) Œ Nearest Neighbor – ng i láng gi ng g n nh t s d ng giá tr c a pixel g n nh t mà không c n ti n hành n i suy Œ Bilinear – hàm song n: ti n hành n i suy n tính s d ng giá tr c a b n pixel Cubic Convolution – xo n l p ph ng: s d ng hàm l p ph ng v i giá tr c a 16 pixel đ ti n hành n i suy Bài t p N n ch nh hình h c nh SPOT v i b n đ n n vùng Ngh a M hi n th M hi n th nh nh nghiado_spot95 hi n th (Display) #1 M file d li u Vector GIS shape file format Shape file format v i đ ng giao thông (roadpoly.shp & roadlines.shp) & sông su i (hydro.shp) s đ c s d ng làm Base (d li u chu n): Trên Menu nh p ch n Vector\ Open Vector File tìm đ n th m c …\RS_Training_Dataset\ DATA nh p ch n đ ng th i file d li u nh p vào Open H p tho i Avaiable Vectors List hi n ra: 26 @2008 T v n GeoVi t – www.geoviet.vn ô 27 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành M ts b ng d n th c hành c quy trình n n nh v i b n đ n n a) Ch n Map\ Registration\ Select GCPs: Image to Map b) Ch n h quy chi u UTM, Zone 48, datum WGS 84 & pixel size 15 m h p tho i Image to Map Registration Hình 2.4: H p tho i l a ch n File d li u Vector Nh p vào Select All Layers sau nh p Load Selected H p tho i Load Vector… hi n ra: Hình 2.7: Hình 2.5: Ch n c a s hi n th File d li u Vector Nh p ch n New Vector Window (hi n th d li u vector c a s khác c a s display # c a s nh) Nh p OK Khi nh n đ c h p tho i Vector Window #1: Cursor Query : Hình 2.6: C a s hi n th File Vector H p tho i Image to Map Registration c) Đ hi n th cho rõ nét l p d li u vector c a s Vector Window #1: Cursor Query ch n Edit\ Edit Layer Properties H p tho i Edit Vector Layers hi n cho phép thay đ i cách hi n th l p d li u vector Hình 2.8: 28 @2008 T v n GeoVi t – www.geoviet.vn H p tho i Edit Vector Layers 29 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành ng d n th c hành d) Đ l a ch n l p d li u đ c dùng làm d li u chu n c a s Vector Window #1: Cursor Query ch n Option\ Select Active Layer Ch n l p d li u c n làm vi c e) M i m GCP đ c nh p b ng cách đ a tr đ n v trí thích h p nh nh p to đ b n đ qua bàn phím vào ô “E” & “N”, sau kích vào nút “Add” Đ ý h p tho i Avaiable Vectors List hi n th thêm l p d li u : L p d li u ch a m kh ng ch Ngoài có th nh p m kh ng ch b ng cách t đ ng: Kích trái chu t lên nh sau kích trái chu t lên v trí t ng d li u vector Kích chu t ph i ch n Export Map Location f) Sau ch n đ hài lòng v i s GCPs & sai s RMS (thông th ng RMS < 1.0 pixel) H p tho i Image to map registration có d ng nh hình Hình 2.10: H p tho i n n nh Registration Parameters Đ t tên cho file n n nghiado_spotToVector Nh p OK i) Hi n th nh đ c n n lên m t Display #2 m i đ quan sát đánh giá ch t l ng n n Đ c shapefile (vector) ch a đ ng sông vào Display #2 m i b ng cách h p tho i Avaiable Vectors List ch n file vector đ ng sông Sau nh p Load Selected C a s Load Vect… hi n ch n c a s hi n th Display #2 Hình 2.11: L a ch n c a s hi n th File Vector Hình 2.9: H p tho i Ground Control Points Selection Image to Map GCP List n n nh v i b n đ n n g) B t đ u n n nh b ng cách h p tho i Ground Control Points Selection ch n Option\ Wap File… H p tho i Input Warp Image hi n ra: Nh p ch n nh c n n n (nghiado_spot95) Nh p OK h) H p tho i n n nh Registration Parameters theo vector hi n Nh p thông s nh hình d i: Nh p OK Khi l p d li u Vector s đ c hi n th c a s Display #2 đ ng th i xu t hi n thêm h p tho i #2 Vector Parameters: Cursor…N u t t h p tho i l p Vector s t t kh i c a s Display #2 30 @2008 T v n GeoVi t – www.geoviet.vn 31 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành Nh p OK H p tho i Available Vectors List nh p ch n Layer: CauGiay.shp ti p theo nh p Load Selected Sau ch n m hi n th file Vector lên c a s nh Trên menu ch n Basic Tools\ Resize Data (Spatial\Spectral) xu t hi n h p tho i Resize Data Input File Nh p ch n File nh sau nh p vào Spatial Subset H p tho i Select Spatial Subset hi n ra: Hình 2.12: H p tho i Vector Parameters j) Đóng file sau k t thúc 2.6 C t nh N u m t t m nh vi n thám mà b n ch nghiên c u m t khu v c c a t m nh b n s d ng ch c Resize Data (Spatial\Spectral) c a ENVI C t nh có th d a m t khu v c, m t ranh gi i … 2.6.1 C t nh theo m t khu v c nghiên c u M File nh nghiado_spot95toAster123geo M file d li u vector Caugiay.shp Trên menu ch n Vector\ Open vector file nh p ch n file d li u Caugiay.shp Khi xu t hi n h p tho i Import Vector Files Parameters Ch n thông s nh hình d i: Hình 2.14: H p tho i l a ch n không gian nh Nh p ch n ROI\EVF H p tho i Subset Image by ROI\EVF Extent xu t hi n: Hình 2.15: H p tho i l a ch n theo ROI\EVF Nh p ch n l p d li u EVF:Layer: CauGiay.shp Nh p OK đ đóng h p tho i ch n đ ng d n l u nh v i tên AnhCauGiay Hình 2.13: H p tho i Import Vector Files 32 @2008 T v n GeoVi t – www.geoviet.vn 33 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành Xu t hi n h p tho i #2 Mask Definition: Nh p ch n Option\ Import EVFs… Nh p ch n select l p d li u Layer: CauGiay.shp Hình 2.18: L a ch n File EVF Quay tr l i h p tho i #2 Mask Definition Ch n đ AnhCauGiay_hanhchinh Hình 2.16: ng d n l u k t qu v i tên nh c t khu v c C u Gi y 2.6.2 Hi n th khu v c Qu n C u Gi y theo đ n v hành M file nh AnhCauGiay M File vector CauGiay.evf cho hi n th lên c a s Display #2 (C a s khu v c v a c t) Trên menu ch n Basic Tools\ Masking\ Build Mask Nh p ch n c a s hi n th Display #2 Hình 2.19: H p tho i Mask Definition Nh p Apply c a s Available Bands List xu t hi n l p Mask Band file AnhCauGiay_hanhchinh Ti p theo menu ch n Basic Tools\ Masking\Apply Mask H p tho i Apply Mask Input File hi n Nh p ch n AnhCauGiay sau nh p ch n Select Mask Band C a s Select Mask Input Band hi n nh p ch n Mask Band Hình 2.17: L a ch n c a s hi n th Mask 34 @2008 T v n GeoVi t – www.geoviet.vn 35 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành ng d n th c hành 2.7 Ghép nh (mosaicking) Ghép nh ngh thu t đ k t h p nhi u nh thành m t nh ghép Ghép có th d a to đ pixel nh ho c d a to đ b n đ đ t o thành nh ph vùng đ a lý r ng h n Hình 2.20: H p tho i Apply Mask Input File hi n Hình 2.21: H p tho i l a ch n Mask Các b c ghép nh n n (georeferenced mosaicking) nh sau: a) Đ c hi n th nh AnhCauGiay; Khuvucnoithanh b) Kh i đ ng vi c ghép nh b ng cách ch n Map\ Mosaicking\ Georerefenced (ghép nh theo t a đ ) t menu c a ENVI H p tho i Map Based Mosaic xu t hi n c) Đ c vào file nh c n ghép v i menu Import h p tho i: đ u tiên đ c vào file nh AnhCauGiay , sau file nh Khuvucnoithanh Nh p OK đ đóng hai h p tho i Select Mask Input Band Apply Mask Input File H p tho i Apply Mask Parameters hi n cho phép đ t tên d li u nh( đ t v i tên AnhCauGiay_spot_hanhchinh) Nh p OK Nh n đ c k t qu : Hình 2.23: H p tho i Map Based Mosaic d) Kích vào “Apply” menu File menu c a h p tho i đ ti n thành ghép nh Đ t tên nh ghép Hanoi_mosaic hi n th m t Display m i đ ki m tra k t qu ghép Hình 2.22: nh khu v c C u Gi y 36 @2008 T v n GeoVi t – www.geoviet.vn 37 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành ng d n th c hành ENVI cho phép ta tăng c ng nh d a theo m t nh đ c tăng c ng s d ng ch c Histogram Matching hay cho ng i dùng t tăng c ng d a đ th theo hàm toán h c đ nh s n thông qua ch c Interactive Stretching K thu t tăng đ nét ch thay đ i cách hi n th hình (tăng gi m đ t ng ph n) mà không làm thay đ i giá tr th t c a pixel Các k thu t stretching g m ph ng pháp n tính không n tính Th v i nh nghiado_ASTER & nghiado_SPOT95 quan sát hi u khác c a k thu t ng 2.8.2 L c nh (Filter) L c nh nh m m c đích tăng đ t ng ph n c a nh đ th hi n rõ h n m t s đ i t ng c n quan tâm, nh ng làm thay đ i giá tr pixel t o file nh l c m i L c nh thông d ng có th th c hi n m t cách t ng tác cho nh hi n th b ng cách ch n enhance\ filter\ t menu c a c a s hi n th Th Hình 2.24: nh đ v i b l c làm tr n, làm rõ nét (sharpening) l c median đ làm s c nét ho c c ghép t hai nh làm m n nh 2.8 K thu t t ng đ nét l c nh (convolutions) 2.8.1 Tăng đ nét Đ tăng đ nét c a nh m t cách t ng tác t c a s nh đ c m , ch n Enhance, m t danh sách s s cho ta ch n di n tích đ c tăng c ng c a s Image, Zoom hay Scroll theo ph ng pháp Œ Linear Tuy n tính: s d ng giá tr nh nh t l n nh t c a nh đ th c hi n phép giãn n tính Ph ng pháp áp d ng phù h p cho nh có giá tr Œ Linear 255 Tuy n tính 255: ph ng pháp s hi n th giá tr th c pixel c a nh theo giá tr hi n th c a hình t đ n 255 Œ Linear 2% Tuy n tính 2%: ph ng pháp tăng c ng n tính s c t b t 2% c a đ u d li u đ tăng kh hi n th nh Œ Gaussian: ph ng pháp tăng c ng nh s d ng giá tr đ xám trung bình 127 đ l ch chu n c a d li u đ tăng c ng Œ Equalization – Cân b ng: ph ng pháp s kéo giãn cân b ng đ th c a d li u đ c hi n th Œ Square Root – Căn b c hai: ph ng pháp s tính b c hai c a đ th đ u vào sau m i th c hi n giãn n tính T menu c a ENVI ch n Filter, b n có th th c hi n đ c phép l c nh tr c ti p file nh v i s k thu t l c thông d ng Filter\ Convolutions (v i thông s kích th c c a s l c d ng l c khác nhau) k thu t l c nâng cao nh texture, Fourier cho nh ng nh có đ nhi u cao Th s d ng thay đ i thông s c a l c convolutions quan sát hi u ng nh 2.9 K thu t u ch nh bi n đ i nh M t s ch c x lý / bi n đ i nh nâng cao s đ c gi i thi u thêm l p n u th i gian cho phép đ c đ a vào ph n d án m u (Ph n V d i đây) - Chu n hoá c p đ xám giá tr pixels (calibration) ph i màu - Bi n đ i nh theo ch đ màu HSV - K thu t thành ph n (PCA) Các k thu t bi n đ i t o kênh m i thông d ng bao g m: 2.9.1 T l c a kênh nh (band ratios) M File nh nghiado_spot95toAster123geo 38 @2008 T v n GeoVi t – www.geoviet.vn 39 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành Ch n Transform\ Band ratios ch n kênh cho t s m u s đ nh p vào h p tho i Hình 2.26: H p tho i tính ch s NDVI Hi n th kênh m i tính đ c quan sát giá tr pixel v i vùng đ i t ng bi t Sau hi n th c a s nh p ch n Tool\ Color Mapping\ ENVI color Tables Hình 2.25: H p tho i Band Ratio Input Bands Hình 2.27: H p tho i ENVI color Tables Ch n d i màu phù h p Nh n đ Nh p OK ch n l u File d li u c k t qu : 2.9.2 T o kênh ch s th c v t (NDVI) Ch n Transform\ NDVI (Vegetation Index) ch n file nh nghiado_SPOT95toAster123geo đ tính kênh ch s th c v t H p tho i NDVI Calculation Parameters hi n ch n th m c tên file d li u Hình 2.28: 40 @2008 T v n GeoVi t – www.geoviet.vn nh NDVI hi n ch nh màu 41 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h 2.9.3 ng d n th c hành Các phép tính kênh (Band Math) Các kênh t l kênh NDVI có th đ c tính t Band Math Ch n Basic Tools\ Band Math nh p vào công th c tính đ th c hi n vi c tính toán t o kênh m i Ví d , kênh NDVI cho nh SPOT ho c ASTER s tính đ c theo công th c NDVI = (Kênh – Kênh 1) / Kênh + Kênh 1) Hình 2.29: H p tho i tính kênh nh Hình 2.30: H p tho i Variables to Bands Pairings 42 @2008 T v n GeoVi t – www.geoviet.vn [...]...X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành Ch n Transform\ Band ratios và ch n kênh cho t s và m u s đ nh p vào h p tho i Hình 2.26: H p tho i tính ch s NDVI Hi n th kênh m i tính đ c và quan sát giá tr pixel v i các vùng đ i t ng đã bi... GeoVi t – www.geoviet.vn nh NDVI khi hi n ch nh màu 41 @2008 T v n GeoVi t – www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h 2.9.3 ng d n th c hành Các phép tính trên kênh (Band Math) Các kênh t l và kênh NDVI trên có th đ c tính t Band Math Ch n Basic Tools\ Band Math và nh p vào công th c tính đ th c hi n vi c tính toán t o ra các kênh m i Ví d , kênh NDVI cho nh SPOT... Color Mapping\ ENVI color Tables Hình 2.25: H p tho i Band Ratio Input Bands Hình 2.27: H p tho i ENVI color Tables Ch n d i màu phù h p Nh n đ Nh p OK và ch n l u File d li u c k t qu : 2.9.2 T o kênh ch s th c v t (NDVI) Ch n Transform\ NDVI (Vegetation Index) và ch n file nh nghiado_SPOT95toAster123geo đ tính kênh ch s th c v t H p tho i NDVI Calculation Parameters hi n ra ch n th m c và tên file ...X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h II N N CH NH VÀ X ng d n th c hành X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h lên khung đ... www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành ng d n th c hành ENVI cho phép ta... www.geoviet.vn X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h X lí & phân tích d li u vi n thám v i ph n m m ENVI – Tài li u h ng d n th c hành ng d n th c hành d) Đ l a ch n l