Kích “Show List” đ xem danh sách nh ng đi m GCP đã ch n... Sauđó nh p Load Selected... Ti p theo trên menu ch n Basic Tools\ Masking\Apply Mask... Ph ng pháp này áp d ng phù h p cho nh c
Trang 1X 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
20
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
Hà N i, 2008
6/17, Ngõ 139 Nguy n Ng c V - C u Gi y
Hà n i – VI T NAM
Tel./Fax: +84.4.5564 710 Email: geoviet@gmail.com Web: http://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
21
M C L C
II N N CH NH VÀ X LÝ NH V I ENVI 22
2.1 T ng quan v module này 22
2.2 Các files s d ng 22
2.3 Các h to đ thông d ng và nh đ c n n ch nh trong ENVI 22
2.4 N n ch nh hình h c nh v i nh 23
Bài t p 2 N n ch nh nh SPOT v i nh ASTER (n n) vùng Nghĩa Đô HN 23
1 M và hi n th nh SPOT & ASTER 23
2 Kh i đ ng n n nh và đ c các đi m kh ng ch (GCP) 23
3 Làm vi c v i các GCP 25
4 Th c hi n n n nh 25
5 So sánh k t qu n n 27
6 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 3 N n ch nh hình h c nh SPOT v i b n đ n n vùng Nghĩa Đô 27
1 M và hi n th nh 27
2 M các file d li u Vector GIS trong shape file format 27
3 M t s b c trong 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 chính 34
2.7 Ghép nh (mosaicking) 37
2.8 K thu t tăng đ nét và 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 tđi u ch nh và bi n đ i nh 39
2.9.1 T l c a các kênh nh (band ratios) 39
2.9.2 T o kênh ch s th c v t (NDVI) 40
2.9.3 Các phép tính trên kênh (Band Math) 42
Trang 2II N N CH NH VÀ X LÝ NH V I ENVI
2.1 T ng quan v module này
Module th c hành này g m 2 ph n: n n ch nh và 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 và 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 là b n đã có 1 s hi u bi t nh t đ nh v n n ch nh
qua bài gi ng lý thuy t Ph n sau cung c p m t s ch c năng x lý / bi n đ i nh
c s đ giúp cho vi c hi n th và gi i đoán nh
2.2 Các files s d ng
nh dùngđ th c hi n các bài th c hành trong module này g m nh ng nh đ c
cung c p trong th m c \RS_Training_Dataset\DATA nh đã gi i thi u
module 1.Đ 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 và nh đ c n n ch nh trong ENVI
ENVI h tr vi cđ c và x lý các nh đã đ c n n ch nh v nhi u h chi u to đ
thông d ng nh đã gi i thi u trên Ngoài ra ENVI còn có th h tr nh ng h
t ađ do ng i dùng t l p b ng cách đ a vào thay th nh ng thông s đ c tr
trong 6 nhóm h quy chi u, 35 d ng ellipsoids khác nhau và h n 100 datumsđ
đáp ng h u h t nh ng yêu c u v đo v b n đ
“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 1 ph n m m so n
th ođ n gi n) Nh ng thông tin trong file này s đ c s d ng trong file header
g n v i m i d li u nh và 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 năng c a ENVI có th dùng thông
tin nàyđ làm vi c v i d li u nh trong 1 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ý và s a đ i đ thích h p v i hình h c nh c s Nh ng đi m kh ng
ch (GCP)đ c ch n trên c a s v i đ đ phân gi i không gian (c a s Chính)
và c a s Zoomđ i v i c ph ng th c n n nh v i nh và n n nh v i b n đ
n n To đ đ c hi n th cho các đi m GCP trên c nh g c và 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 1 đ n n)
ho c n n tam giác Delaunay Vi c resampling (tính toán và ghi d li u t ng pixel
lên khung đã đ c n n) h tr các ph ng pháp nearest neighbor (pixel g n
nh t), n i suy tuy n tính b c 2 và n i suy b c 3 (cubic convolution) So sánh gi a
nh g c (n n) và nh n n s d ng kh năng Dynamic Overlay c a ENVI cho phép đánh giá đ chính xác c a phép n n m t cách nhanh chóng
2.4 N n ch nh hình h c nh v i nh
Ph n này c a module s d n d t các 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 cùng 1 khung hình h c
Bài t p 2 N n ch nh nh SPOT v i nh ASTER (n n) vùng Ngh a ô HN
- Th c hi n các b cđã đ c gi i thi u module “D li u vi n thám và ph n
m m ENVI”đ m và hi n th nh nghiado_spot95 trong hi n th (Display) #1
& nghiado_aster123geo t i hi n th #2;
- Hi n th v trí / giá tr con tr đ quan sát to đ c a các pixel trong nh SPOT
& nh ASTER;
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) khi ch n Base Image ( nh chu n) và Display #1 ( nh SPOT) khi ch n Warp Image ( nh c n n n);
Hình 2.1: C a s l a ch n nh n n và nh g c
Kích vào “OK”đ b t đ n quá trình n n nh
Trang 3c) H p tho i Ground Control Points Selection s xu t hi n và t ng đi m kh ng
ch (GCP) s đ c nh p vào b ng cách đ t v trí con tr trên c 2 nh vào đúng
cùng 1 v trí trên m tđ t (giao đi m c a đ ng, đi m đ c bi t ) (Hình 2.2)
d) Di chuy n con tr trong 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 con tr trong 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í con tr trong c a s Zoom và đi 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 năng 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 đi m GCP đó vào danh sách Kích “Show List”
đ xem danh sách nh ng đi m GCP đã ch n Th làm v i vài đi m GCP đ có
đ c k năng ch n GCP
g) Chú ý vào danh sách c a nh ng đi m th c và đi m d tính (predicted), khi s
đi m đã ch n l n h n 4, sai s RMS s đ c tính và hi n th
h) Ch n Options\ Clear All Points trong h p tho i Registration GCP List đ xoá đi
nh ngđi m b n đã ch n;
i) Ch n File\ Restore GCPs from ASCII trong h p tho i Ground Control Point
Selection và kích vào tên file nghiado_spot95_toAster123geo.pts
j) Kích vào “OK” đ đ c danh sách nh ng đi m GCP đã l u tr c;
k) Kích vào t ng đi m GCP trong h p tho i Image to Image GCP List và ki m tra
v trí c a cácđi m đó trên 2 nh, nh ng t a đ th c t và t a đ d tính, và sai
s RMS Quan sát sai s RMS t ng trên h p tho i Ground Control Point
Selection
25
Hình 2.2: H p tho i Ground Control Point Selection và h p tho i Image to Image GCP List
khi n n ch nh hình h c nh v i nh
Nh ng miêu t d iđây có th đ c dùng nh thông tin tham kh o:
- V trí c a m i đi m GCP có th đ c s a b ng cách ch n chúng trên danh sách h p tho i GCP List và s a chúng trong h p tho i GCP Selection Ho c
là b n có th nh p vào m t v trí m i c a pixel ho c là 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 trong h p tho i;
- Kích vào nút “On/Off” trên h p tho i GCP List đ lo i b đi m GCP đ c
ch n ra kh i ph ng trình tính toán h s n n và công th c tính sai s Nh ng
đi m GCP đó ch b b qua khi tính toán ch ch a b xoá kh i danh sách và
b n có th s d ng l i chúng khi kích vào nút “On/Off” 1 l n n a
- Kích vào nút “Delete” đ xoá đi m GCP ra kh i danh sách
- Đ t v trí con tr trong 2 c a s zoom và kích vào nút “Update” đ thay đ i
đi m GCP đ c ch n v v trí m i c a con tr
- Nút “Predict” cho phép tính toán nh ng đi m GCP m i d a trên mô hình n n đang dùng
a) Th đ t con tr vào m t v trí m i trên nh ASTER và kích vào nút “Predict”, khiđó v trí con tr trên 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í chính xác sau đó có th đ c tinh ch nh b ng cách di chuy n m t chút v trí pixel trên nh SPOT
c) Kích vào nút “Add Point” đ nh p thêm đi m GCP cho danh sách
nh có th đ c n n ch cho nh ng kênh đang hi n th ho c cho t t c các kênh
c a nh cùng 1 lúc, đây là ví d c a vi c ch n n nh ng kênh đang hi n th :
a Ch n Options\ Warp File… trong h p tho i GCP Selection H p tho i Input Wrap Image hi n ra ch n nh c n n n (nghiado_Spot95) r i nh p OK
b Khi h p tho i Registration Parameters xu t hi n, ch n “RST” cho Warp Method và ph ng th c resampling là “Nearest Neighbor”
c Ch n tên file nh n n là nghiado_SPOTtoASTER123geo b ng cách tích vào File
và nh p vào chooseđ l a ch n th m c ch a file d li u và kích “OK”
nhđ c n n s đ c hi n th trong Available Bands List sau khi quá trình
n n ch nh k t thúc
d Th l p l i v i ph ng pháp resampling Bilinear Interpolation và Cubic Convolution
e Đ t tên các file k t qu là nghiado_SPOTtoASTER123geo_b &
nghiado_SPOTtoASTER123geo_c
Trang 4f Th c hi n vi c n n đa th c b c 1 và n n Delaunay Triangulation v i Cubic
Convolution
g Đ t tên các file k t qu là nghiado_SPOTtoASTER123geo_pc &
nghiado_SPOTtoASTER123geo_dtc
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 và t nh ti n nh
Œ Ph ng pháp Polynomial – Hàmđa th c: ph ng pháp này cho k t qu t t h n
ph ng pháp RST, v i yêu c u s v s đi 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 các đi m kh ng ch làm các đ nh c a các
tam giác khôngđ u và ti n hành n i suy
Đ ti n hành n n nh ta cũng ph i l a ch n m t trong ba ph ng pháp tái chia
m u – Resampling sao 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 tuy n: ti n hành n i suy tuy n tính s d ng giá tr c a b n
pixel
pixelđ ti n hành n i suy
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 và ch n File\ Close Selected File trên c a s Available Bands List
b) Đ c các file nh đã n n trên và l n l t hi n th lên các Display m i
c) Đ t con tr vào c a s Chính và kích vào phím ph i chu t và sau đó ch n File\ Link\ Link Displays Kích vào “OKđ liên k t nh ASTER và 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 các ph ng pháp n n và ph ng th c resampling
Làm hi n lên h p tho i Cursor Location/ Value và di chuy n con chu t trong 1 c a
s và quan sát hi u ng c a các ph ng phán n n khác nhau 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 đ trong h t a đ UTM Zone 48, datum WGS 84 d a trên d li u GIS trong 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 trong th m c làm vi c) s đ c s d ng làm Base (d li u chu n) và to đ pixel trên nh nghiado_Spot95 s đ c n n ch nh v t a đ b n đ To đ c a t ng
đi m GCP trên b n đ s đ c đ c sau khi hi n th nh ng file đó trên ph n m m ArcView
Bài t p 3 N n ch nh hình h c nh SPOT v i b n đ n n vùng Ngh a ô
M và hi n th nh nghiado_spot95 trong hi n th (Display) #1
2 M các file d li u Vector GIS trong 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 3 file d li u trên và nh p vào Open H p tho i Avaiable Vectors List hi n ra:
Trang 5Hì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.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 trên 1 c a s khác c a s
display # 1 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
29
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 là 15 m trong h p tho i Image to Map Registration
Hình 2.7: H p tho i Image to Map Registration
c) Đ hi n th cho rõ nét các l p d li u vector trên c a s Vector Window #1: Cursor Query ch n Edit\ Edit Layer Properties H p tho i Edit Vector Layers
hi n ra cho phép thayđ i cách hi n th các l p d li u vector
Hình 2.8: H p tho i Edit Vector Layers
Trang 6d) Đ l a ch n l p d li u đ c dùng làm d li u chu n thì trê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 1 đi m GCP đ c nh p b ng cáchđ a con tr đ n v trí thích h p trên nh
và nh p to đ b n đ qua bàn phím vào ô “E” & “N”, sau đó kích vào nút
“Add” Đ ý trên h p tho i Avaiable Vectors List hi n th thêm 1 l p d li u
<Registration #1 GCPs>: L p d li u ch a cácđi m kh ng ch
Ngoài ra còn có th nh p cácđi m kh ng ch b ng cách t đ ng: Kích trái chu t
lên trên nh sauđó kích trái chu t lên v trí t ng trên d li u vector Kích chu t
ph i ch n Export Map Location
f) Sau khi đã ch n đ và hài lòng v i s GCPs & sai s RMS (thông th ng RMS <
1.0 pixel) H p tho i trong Image to map registration có d ng nh hình
Hình 2.9: H p tho i Ground Control Points Selection và Image to Map GCP List trong
n n nh v i b n đ n n.
g) B t đ u n n nh b ng cách trên 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 ra Nh p các thông
s nh hình d i:
Hình 2.10: H p tho i n n nh Registration Parameters
Đ t tên cho file n n là nghiado_spotToVector Nh p OK
i) Hi n th nhđ c n n lên m t Display #2 m i đ quan sát và đánh giá ch t
l ng n n.Đ c shapefile (vector) ch a đ ng và sông vào Display #2 m i này
b ng cách trên h p tho i Avaiable Vectors List ch n 3 file vectorđ ng và sông Sauđó nh p Load Selected C a s Load Vect… hi n ra 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
Nh p OK
Khiđó các l p d li u Vector s đ c hi n th trên c a s Display #2 và đ 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 này thì các l p Vector s t t kh i c a s Display #2
Trang 7Hình 2.12: H p tho i Vector Parameters
j) Đóng các file sau khi k t thúc
2.6 C t nh
N u trên m t t m nh vi n thám mà b n ch nghiên c u m t khu v c nào c a t m
nh thì b n s d ng ch c năng Resize Data (Spatial\Spectral) c a ENVI C t nh có
th d a trên 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 các thông s nh hình d i:
Hình 2.13: H p tho i Import Vector Files
33
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) khi 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.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 các h p tho i và ch nđ ng d n l u nh v i tên AnhCauGiay
Trang 8Hình 2.16: 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 chính
M file nh AnhCauGiay
M File vector CauGiay.evf và 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.17: L a ch n c a s hi n th Mask.
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đ ng d n và l u k t qu v i tên AnhCauGiay_hanhchinh
Hình 2.19: H p tho i Mask Definition
Nh p Apply khiđó trên c a s Available Bands List xu t hi n l p Mask Band trong file AnhCauGiay_hanhchinh
Ti p theo trên menu ch n Basic Tools\ Masking\Apply Mask H p tho i Apply Mask Input File hi n ra Nh p ch n AnhCauGiay sauđó nh p ch n Select Mask Band C a s Select Mask Input Band hi n ra nh p ch n Mask Band
Trang 9Hình 2.20: H p tho i Apply Mask
Input File hi n ra Hình 2.21: H p tho i l a ch n Mask
Nh p OKđ đóng hai h p tho i Select Mask Input Band và Apply Mask Input File
H p tho i Apply Mask Parameters hi n ra 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.22: nh khu v c C u Gi y
37
2.7 Ghép nh (mosaicking)
Ghép nh là ngh thu tđ k t h p nhi u nh thành m t nh ghép Ghép có th d a trên to đ pixel nh ho c d a trên to đ b n đ đ t o thành 1 nh ph vùng đ a
lý r ng h n
Các b c ghép nhđã n n (georeferenced mosaicking) nh sau:
a) Đ c và 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 chính 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 trên h p tho i: đ u tiên đ c vào file
nh AnhCauGiay , sauđó file nh Khuvucnoithanh
Hình 2.23: H p tho i Map Based Mosaic
d) Kích vào “Apply” trong menu File trên thanh menu c a h p tho i đ ti n thành ghép nh.Đ t tên nh ghép là Hanoi_mosaic và hi n th trên m t Display m i đ
ki m tra k t qu ghép
Trang 10Hình 2.24: nh đ c ghép t hai nh
2.8 K thu t t ng đ nét và 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
s Image, Zoom hay Scroll theo các ph ng pháp
Œ Linear Tuy n tính: s d ng giá tr nh nh t và l n nh t c a nh đ th c hi n
phép giãn tuy n tính Ph ng pháp này áp d ng phù h p cho nh có ít giá tr
Œ Linear 0 255 Tuy n tính 0 255: ph ng pháp này s hi n th các giá tr th c
pixel c a nh theo giá tr hi n th c a màn hình t 0đ n 255
Œ Linear 2% Tuy n tính 2%: ph ng pháp tăng c ng tuy n tính s c t b t 2%
c a 2đ u d li u đ tăng kh năng hi n th nh
Œ Gaussian: ph ng pháp này tăng c ng nh s d ng giá tr đ xám trung bình
là 127 vàđ l ch chu n c a d li u là 3 đ tăng c ng
Œ Equalization – Cân b ng: ph ng pháp này 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 này s tính căn b c hai c a đ th đ u
vào sauđó m i th c hi n giãn tuy n tính
ENVI còn 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 năng Histogram Matching hay cho ng i dùng t tăng c ng d a trên
đ th và theo các hàm toán h c đ nh s n thông qua ch c năng Interactive
Stretching.
K thu t tăng đ nét này ch thay đ i cách hi n th trên màn hình (tăng gi m đ
t ng ph n) mà không làm thay đ i các giá tr th t c a pixel Các k thu t stretching g m các ph ng pháp tuy n tính và không tuy n tính
Th v i các nh nghiado_ASTER & nghiado_SPOT95 và quan sát các hi u ng khác nhau c a các k thu tđó
2.8.2 L c nh (Filter)
L c nh cũng 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 các giá tr pixel và t o ra 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đang hi n th
b ng cách ch n enhance\ filter\ t thanh menu c a c a s hi n th chính Th
v i các b l c làm tr n, làm rõ nét (sharpening) và l c medianđ làm s c nét ho c làm m n nh
T menu chính c a ENVI ch n Filter, b n có th th c hi nđ c phép l c nh tr c
ti p trên file nh v i 1 s k thu t l c thông d ng trong Filter\ Convolutions (v i các thông s là kích th c c a s l c và d ng l c khác nhau) và các k thu t l c nâng cao nh texture, Fourier cho nh ng nh cóđ nhi u cao Th s d ng và thay
đ i các thông s c a l c convolutions và quan sát các hi u ng trên nh
2.9 K thu t đi u ch nh và bi n đ i nh
M t s ch c năng x lý / bi n đ i nh nâng cao s đ c gi i thi u thêm trên l p
n u th i gian cho phép vàđ 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 chính (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 các kênh nh (band ratios)
M File nh nghiado_spot95toAster123geo
