Kết hợp giữa thế mạnh giải đoán, xử lý ảnh viễn thám của phần mềm Envi và thế mạnh xử lý số liệu, biên tập bản đồ của Arcgis đã cho hiệu quả cao và khách quan để giải quyết bài toán biến động sử dụng đất.
Trang 1PHAN 1 TAI ANH VE TINH
Dữ liệu ảnh vệ tỉnh cĩ thể được lấy từ nhiều nguồn khác nhau tuy nhiên hiện nay cách được sử dụng phổ biến nhất là tải từ các trang web tải ảnh vệ tỉnh thơng qua mạng Internet Tùy thuộc vào loại ảnh vệ tỉnh, chất lượng và khả năng ứng dụng mà người dùng cĩ thể tải miễn phí hay trả phí cho từng loại dữ liệu Hiện nay cĩ một số trang web thơng dụng như:
http://earthexplore.usgs.gov http://glovis.usgs gov
Để tải ảnh vệ tỉnh chúng ta truy cập một trong số các trang web trên
Đây là giao điện của trang chủ của http://earthexplore.usgs.gov
Muốn download, người dùng cần cĩ một tài khoản dé đăng nhập vào hệ thống Tài khoản này cĩ thê đăng ký miễn phí bằng cách chọn Regited và nhập các thơng tin cân thiết theo yêu câu, sau khi đăng ký xong cân vào Email đăng ký để kích hoạt tài khoản Khi đã cĩ một tài khoản, click vào login sau đĩ điền các thơng tin đăng nhập đề đăng nhập và trang web
Trang 2
vẻ tên, tính chất bức ảnh, tìm bằng cách khoanh vùng hay chọn điểm cụ thể trên map ngồi ra trang web này cịn hỗ trợ lọc về khoảng thời gian, ngày-đêm, lượng mây che phủ, senso Sc ee ype in an address or place ` "Tim bằng cách khoarh vững cụ thể G ì } Cae cr 0 28H an EEE fp tae w 2527 K ln 1 aE ae 10 96 OF W Lan 5 av ar Sau khi lựa chọn xong | ®JECAMSites Lựa chọn dữ liệu vệ tịnh
các yêu cầu về ảnh, chọn Data : #3 Land Cover UPAR
Sets để tùy chọn loại vệ tinh| Và [Z1anasatrcnwe uo kiểu dữ liệu ảnh vệ tỉnh tương #'@ L8 oLrls
ứng '® | L8 OLUTIRS Pre-WRS-2 Aiwa hai te, ` L7 ETM+ SLC-off(2003-present) Thun ung #/@IEl L7 ETM+ SLC-on (1999-2003)
chọn cách khoanh vùng TX ©1517 ETM+ int Ground stations (Search Only)
Hồng Lĩnh sau đĩ lựa chọn 2 #'@ L4-5TM
OB u1-5mss
moc thdi gian là 2003 và 2013
Ứng với 2 mốc thời gian đĩ và
Landsat COR Gi
ý Landsat Legacy
yêu câu cia bai lam téilya chon | Landsat MRLC
anh Landsat 5 ứng với năm (NASA LPDAAC Collections
2003 va Landsat 8 ứng với năm | ` 9'2% Radar Tuy chọn lọc nâng cao
2013 Trong mục Landsat predic Q.20) & Hiển thị kết quả
Archiver của Data Sets tích vào Nnibisaibiiouel các loại dữ liệu ảnh ve tinh EROS
Trang 3tương ứng Sau đĩ chọn Additional Criteria để chọn các tiêu chuẩn lọc nâng cao hoặc chọn Results để hiễn thị các kết quả tìm kiếm
Kết quả sẽ tơng hợp cho chúng ta tất cả các bức ảnh thỏa mãn yêu cau đặt
ra và sắp xếp theo thời gian gan nhất
‘Show Result Controls x Data Set enero expat youresuts » mm 7 ny 10: LT512604720032088C00 AcqusifønDat: 11OCT03 PaN: Tào Row 47 )sdd/v® nut 10: LT512604720032226,¢00 ‘Acquistion Date: 1-AUG03 Jid7LaWSe nut 0: LT512604720022478,C00 ‘Acquistion Date: 1-06C-02 ath 126 Row: 47 Danh sách ảnh L5 năm 2003 'ENHW ULCST281452013297LNUU Coordinates: 7.3443 108 05788 ‘Acquistion Date:24-OCT-13 Path: 126 Row 48 Sudq7lao nity ID: LC81280472013297L6N00 Coordhates: 1878739 108 38843 ‘Acquistion Date: 24-OCT-13 Path: 126 Row 47 )sdzv8a© Entity 10: Lc8 126048201328 1LGNO0 Coordinates: 1734436, 105 04825 ‘Acquistion Date: 08-OCT-13 ath: 128 Row 48 sd⁄v8a© Entity 10: LOB 26047207328 TUGNOO Coordinates: 1878732 105 37688 ‘Acquistion Date: 08-OCT-13 Path: 126 Row 47 Jud ao rst «Provious (5 ¥] Next> La Danh sách ảnh L8 năm 2013
Cĩ thẻ căn cứ vào các cơng cụ như ảnh xem trước, Show footprint, Show
Browse Overlay, Show Metadata and Browse để biết các thơng tin về bức ảnh
đĩ như vị trí, thời gian chụp, chất lượng ảnh đẻ đảm bảo tải đúng dữ liệu mình cẩn Sau khi lựa chọn được ảnh, click vào|biểu tượng tải về để tải ảnh Sẽ cĩ nhiều lựa chọn để tải ảnh như ảnh màu tự nhiên, ảnh nhiệt Lựa chọn Download Lever
Trang 4[CƠ rcsat.ook stra Color image (22 mB)
LandsatLook Thermal Image (2.6 MB)
LandsatLook Images with Geographic Reference (4.9 MB)
Ea Ea) Level 1 Product (106.8 MB)
Tải ảnh Landsat 5 Hà Tĩnh năm 2003
‹aicseLook “Natural Color” image (2.2 MB)
LandsatLook "Thermar” Image (961.4 KB)
LandsatLook "Quality" Image (190.6 KB) (7 Click dé tai anh dinh dang TIFF
LandstLook images win Geographic Reterg
EES Level 1 GeoTlFF Data Product (765.4 MB)
Tải ảnh Landsat 8 Hà Tĩnh năm 2013
Sau khi quá trình download hồn tất, chúng ta được 2 file ảnh ở 2 thời kỳ năm
2003 và 2013 của khu vực Hà Tĩnh ở đẹ:g nén tar
BOs WOK Dt Baw DỰ Ban Ree
omy Pate oy am Tac cnre arial PORE vustory GB wer selection
© © + th [SE HGwHTO say, vi GÌ | Semen gusty
tame Duemodied Type sice
“ 10/31/2014 9:55PM File folder
Trang 5Giải nén và đổi tên thư mụt ta được nhiễu band ảnh tương ứng như sau:
GL] Steere eee oe con Tos [Ợt Deeeesee =)
Cactep rong Anh
= > compar» HOCTAO) › goay ' NMESĐEDHDXUC HIẾU GHI5 1 downloade ^ 1 Recent places Lives ước (TSI2E0.72HB)19 (T5LXQG200J LTSIGỜSGĐĐI9 LTEBEOOAODS (TỚIMSOZAOĐĐĐ TEIS6OGEDMS B Videos 2% Homesrour
Computer (TSI2EƠ/720019 IEIMĐQƠAOMS LIAIĨƠGMM0Đl2 LTSI2GOKTIONIG ƯIIMQGAOOMA - HEADMEGIE ao) GĐCPS.87.If /08006CPue O6(C0QMTLb€ G8CSOVEpg - G8UCĐĐVEReL ca SETUP 0) ca HOCTAP () aout fa Removable Dik
TS on mT] [OF Deemer sos
Ca : “Các tệp tong ảnh Landsat 8 tài về
Trang 6|PHẢN 2 GOI ANH
- Bước 1 Mở các kênh ảnh cần ghép
Khởi động phần mêm Envi 4.5 Xuất hiện thanh Menu chính của Envi
p Basic Tools Clasification Transform Fier Spectral Map Vector Topographic Radar Window Help Trên thanh menu của phần mềm Envi chon: File\Open External File\
Landsat\ Geo TIFF dé mé file theo dinh dạng TIF của ảnh Landsat
[File | Basic Tools Classification Trans Open Image File
Open Vector File Mo anh theo dinh dang Geo TIFF
Open Remote File
Open External File >) Landsat
Open Previous File >| SPOT Geol
Launch ENVI Zoom QuickBird ron NI › NLAPS
Edit ENVI Header Worldview › wae
Generate Test Data OrbView-3 II
Data Viewer IRS › ESA CEOS
Save File As › AVHRR k
R SeaWIFS ›
Import from IDL Variable fos ,
Export to IDL Variable = ,
Compile IDL Module ewvisat ,
IDL CPU Parameters ALoS ,
Tape Utilities >) CARTOSAT-1 ›
‘Scan Directory List ADM0
Change Output Directory ATSR
DMSP (NOAA) Save Session to Script FoRMOSAT.2
Execute Startup Script
eee Thermal ›
Radar >
ENVI Queue Manager = ,
ENVI Log Manager EXMexin 0
Close All Files uses »
Preferences LASLIDAR »
Trang 7+ File name: chỉ ra đường dẫn đến file ảnh cân gĩi, đuơi *TIE
+ Chọn các kênh ảnh cần gĩi đã giải nén/ Open € ©` + 4|[Ư « HGCTAP& › guay » 2003 w/e | | Search 2003 2 Organize > Newfoldef BF Videos ˆ XÃ Homegroup ihe Local Disk (C9 LT5126047200319| '0BIC00 B2,TIE |LT5I26047200319| | LT5126047200319| | LT5126047200319 08JC00,B3TIF 08JC00 B4.TIE ea SETUP (0) ca HOCTAP (E) ca GIẢI TRỊ (E) = Removable Disk LT5126047200319 | | LT5126047200319 | | LT5126047200319 oBICO0 BS TIF | | 08UC0.B6TF | | oBic Chọn các kênh ảnh LT5 cần gĩi File name: |°LT512604720031908)C00_61.TIF* "LTS126087 ¥ GH Network MB NHAT-BLACK pC.PC - Bước 2 Ghép các kênh ảnh
Trên thanh Menu của phần mềm Envi lựa chọn:
Basic Tools\Layer Stacking, xuat hién hop thoai Layer Stacking Parmeters
+ Chon Import File/ Chọn các Band cân gĩi sau đĩ chon OK
Resize Data (Spatial/Spectral)
Trang 8Giao dign chinh cia Layer Stacking Parameter nhw sau: Độ phân giải Output Resuitto © Fle C Enter Output Filename Choose, I 0K | cance |
+ Lựa chọn Import File để thêm các kênh ảnh Trong hộp thoại Layer
Stacking Input File lựa chọn các kênh cần gĩi sau đĩ chọn OK Thơng thường
Trang 9+ Chọn Reođer Files ấn chuột trái để sắp xếp các band theo thứ tự tăng dan sau đĩ bấm OK + 9 Reorder Files ‘Selected Files for Layer Stackang: ([751260472003 190B/C00_B7 TIF [Band 7] ip xếp các kênh ảnh theo thứ tự nhất định -2: LT512604720031908JC00_B2 TIF [Band 1] 3: LT512604720031908JC00_B3.TỊF [Band 1] 4: LT6126047200319084C00, B4 TIF [Band 1 | ‘5: LT51260472003190BJC00_B5.TIF [Band 1} 6: LT5126047200319081C00_86 TIF [Band 1] | 7: LT512604720031908JC00_B7.TIF [Band 1} = ¬ OưputResultto @ Fle CC Memoy
Enter Output Filename Choose |
XX Pxel Sge|3U.U0000000 Meters Meters 1: 7512604720031908iC00_B1.TIF [Bend 1) | Resamping [Nearest Neighbor _~| _øK |_cane |
Lựa chọn các thơng tin giống với hình đưới đây sau đĩ chọn Choose\chỉ ra
đường dẫn lưu kết quả, Ok để thực hiện gĩi ảnh
Trang 10
Tùy thuộc vào số lượng kênh ảnh cũng như
dung lượng ảnh mà quá trình gĩi ảnh cĩ thể mất | "= 9ts= ` — ape and 5160472003 một vài phút tye and 1LT51280072003 Layer (Band 1751260472003, Cape Sond 1475126047003 tape va 1115126017003 Ly Bạnd 1LT5I280072003 a Input Fie : E:\giualy\2003\LT512604720031908/C00_85.TIF Fle: E`\guaky\anhgo2003 Layer Gand (LT51260472003 Ta) Ls126047200 1 aCOO_E7 TF JOutput Input Fie 5 of 7
Các kênh ảnh được gối trong "ảnh gĩi 2003" |] LT612604720021908ico0_ © Band 1 s5 TIE v
| © Gray Scale © RGB Cols
Sau khi gĩi xong, xuất hiện hộp thoại ‘Selected Band
Available Bands List anh g6i sé hién thị ngay |=" venice
trên hộp thoai Available Bands list voi day du các band đã gĩi, ảnh sau khi gĩi sẽ được đưa về cùng độ phân giải [p=Pmmmmama | @
File Basic Tools Classification Transform Filter
Trang 11PHAN 3 CAT ANH
Để cắt anh ta cdn cé file vecter địa gi
hành chính của khu vực đĩ, dữ
này cĩ thể xin tại các cơ quan hoặc download trên các trang web, sau đây là quy trình tách file địa giới hành chính của khu vực TX Hồng Lĩnh từ file vecter hành
chính cấp huyện của Việt Nam
Trên thanh Menu cia phan mém Envi chon File\Open vector File dé mé file
vecter hanh chinh Viét Nam, hép thoai Select Vecter Filenames xuat hién, lya chon file vecter céin mở rồi chọn Open
Hộp thoại Import Vecter Files
Parameters xuất hiện, lựa chọn file vecter, đường dẫn lưu các thay đổi của file, thay đổi lưới chiếu, phép chiếu và đơn vị hiển
thị Trong trường hợp này ta thường giữ
nguyên các giá trị, chọn ØK Hộp thoại Available vectors List xuất hiện ta chon
file vecter can hién thj sau đĩ chon Load
Selected H6p thoai Load vecter xuat
hiện, lựa chọn cửa số hiền thị và chọn OK để hiển thị file
& _ import Vector Files Parameters MEH
Selected input Fies:
Trang 12File vecter hành chính| cấp huyện, thị xã của Việt Nam Edit Options Help
được hiển thị như sau:
Trên thanh Menu của cửa số Vector Window #1: Cursor Query chọn lệnh: Edit/View/ Query Attributes
Vector Window #1 : Cursor Query
Mode [Edit] Options Help Edit Layer Properties
View/Edit Query Atibutes Query Attributes
Hộp thoại Leyer Attributes xuat hién, chứa đây đủ thơng tin các lớp cĩ trong
file vecter đĩ, tìm đến và chọn lớp Hồng Lĩnh, trên thanh Menu của hộp thoại Leyer
Attributes chon File/Save Selected Records to New Layer Hơp thoại Save
Selected Records Layer xuất hiện, nhập tên và thư mục lưu, chọn OK
Fle Options
Save Selected Record to New Laer Tae NET
Save i Record to ASC
Trang 13Sau khi cĩ file địa giới hành chính của Hồng Lĩnh, tiến hành cắt ảnh theo các bước như sau:
- Bước 1 Mở ảnh cần cắt (ảnh Hà Tĩnh đã gĩi)
Trên thanh Menu của phân mềm Envi chọn: File\Open Image File dé mé file
ảnh Hà Tĩnh đã gĩi
- Bước 2 Mở tiếp file vector ranh giới hành chính Hồng Lĩnh (đuơi evf)
File\Open vector File
- Bước 3 Chọn file Veeter\ Load Selected\Display #1\OK trên hộp thoại Vecter' Parameters chọn Apply để hiển thị chồng file vecter lên file ảnh Hà Tĩnh đã #1 Vector Parameters:C — File Mode Edit Options Help Window: Image © Scol © Zoom © Of
‘Avalable Vector Layers
File vecter duge hién thi chéng
lên file ảnh Hà Tĩnh như hình bên
Trén Menu của hộp thoại Available
Vecter List chọn FileExport Active Layer to ROIs ,
@ Available Vectors List - © MEM
[File Options ‘Open Vector File
Export Layers to ROI Export Layers to Shapefile
Export Layers to ArcGiS Geodatabase ‘Open Layers in ENVI Zoom
Trang 14Xuất hiện hộp thoại Expolrt EVE
Leyer to ROI chọn Convert all records of an EVF layer to one ROI \OK Chuyên tồn bộ thơng tin của u ROI - Bước 4 Trên cửa số ảnh Envi: chọn lệnh Overlay\Region of Interest
Hộp thoại ROI Tool xuất hiện, trên
thanh cơng cụ của hộp thoại ROI Tool Vao File\ Subset Data via ROIs, xuất hiện hộp thoại Select Input Files to
Subset via ROI
B===ls=i=is==is=] IEEIIEEEenni|
"Fie Information:
= [EE-EsEoanhoo2p03 [Dene 7821x6951 27 (8S) [r1s1 260472003 1906sC00-86 TIF L151 2604720031 908NC00-85 TEL a CON
Trang 15Chon File vecter lam đữi liệu
chồng xép 6 muc Select Input ROIs, Chon Yes 6 Mask pixels
ouside of ROI?
Chon Choose\chi ra duong dan
lưu kết quả, nhập tên mới của ảnh cắt
'Open\OK để thực hiện cắt ảnh
Sau khi cắt ảnh xong, ảnh sẽ
hiển thị lên trên cùng của hộp thoại
Available Band list File Options Efssrzr
ROI Mask (Layer (Band 1:LT51) RO! Mask (Layer (Band 1:LT51 ROI Mask (Layer (Band 1:LT51 RO! Mask (Layer (Band 1:LT51 RO! Mask (Layer (Band 1:LT51 ROI Mask (Layer (Band 1:LT51
Trang 17[PHAN 4 GIAI DOAN ANH TREN ENVI
- Bước 1 Mở ảnh cần phân loại File\Open Image File
- Bước 2 Trên thanh Menu của ảnh hiển thị chọn:
n Overlay\Region of Interest
[Fie [Overy Enhance Tools Window TT uớon Xuất hiện hộp thoại RƠI Tòt Hoặc Basic Tools\Region of
Classification Interest \ROI Tool
ad 'Xuất hiện hộp thoại ROI Tool
Chọn New Region đề thêm các
Region of Interest q lớp chọn mẫu cĩ trong ảnh
ROI Type: chọn Polygon, Multi Part: Off
Fill: chon None
Lựa chọn cửa số Zoom để chọn mẫu Sau đĩ tiến hành chọn mẫu theo các lớp đã đặt tên Densiy Sice Lop dang
Sree lua cle, Chọn cửa số chọn mẫu
Trang 18- Bước 3 Chọn vùng mẫu |
hấp chuột vào lớp cẩn chọn mẫu Sử dụng cửa số Zoom để chọn, ta phĩng to cửa số này lên để thấy rõ các pixel, nhấp chuột để bắt điểm, hoặc nhấn giữ chuột trái và khoanh vùng để chọn vùng mẫu, nháy chuột phải lần 1 để khép vịng và lần 2 để chọn mẫu, muốn xĩa vùng mẫu đã chọn ta di chuột vào ơ vùng mẫu và nháy chuột giữa Trong cột Color nên chọn các màu tương phản với chính lớp
đĩ để đễ quan sát và thuận tiện khi chọn mẫu Ta cĩ thể thay đổi các kiểu tổ họp màu hiễn thị khác nhau để đễ dàng giải đốn
Trang 19Một số mẫu lựa chọn đối tượng lrên ảnh Landsat 8 té hop màu theo kênh 5-4-3:
Đất trống
Đất đơ thị Đất rừng
Các vùng được chọn sẽ hiện thị trên ảnh, cứ làm tương tự như vậy ta sẽ lần
lượt chọn hết được tất cả vùng mẫu cho các lớp
'Kết quả của việc giải đốn phụ thuộc vào kết quả của việc chọn mẫu vì vậy việc chọn mẫu phải hết sức cân thận và tỉ mi, cần am hiều về các khĩa giải đốn và các đặc trưng phản xạ của đối tượng được chọn Ngồi ra cân lưu ý về thời gian chụp ảnh để biết mức phản xạ phổ của một số đối tượng cĩ tính thay đơi theo mùa, ví dụ như các đối tượng về nơng nghiệp, hoa màu cố tính chất ngắn ngày Các mẫu cần được lựa chọn ở
nhiều vùng khác nhau và rải đều trên cả ảnh, số lượng điểm được chọn càng nhiêu, càng
chỉ tiết thì xác suất giải đốn chính xác càng cao Đồi với yêu câu lớn sau khi chọn mẫu
Trang 20Kiểm tra độ phân tách phổ giữa|các đối tượng: ~ _ Trên hộp thoại ROI Tool, Option\Compute Roi bu Separabili parabilty —_— KỈ:
- Xuất hiện hộp thoại Repo ea RO > mn a ea
Select Input File For Roi pị = —=e-ta 6 em 6 0U Separability: chon anh\OK seonctenon
Reconcile ROS va Map v — Xuất hiện hộp thoại Roi Separability Calarlation\ Select All Items\OK fund Teo 0L | ete rage tom Create Biter Zone rm mm TS Đan ee, ane" (lettries-Watusite, Tranetoraed Divergence) 00) CẮT BÙI [White] 214 poiate: (1 99997201 2 0000000) PAT BONG (Groen) 711 point MAT WGC’ [Hied} 615 points: (1, 99999925 2 00000000)
DAT MAO (Blue) és) points (1 93367184 2 00000000) DAT WONG NOELEF (Yel1ov)"709 pointe (1 33333388 3 00000000) DAT TRONG (Cyan} 73 points: (199999998 200000000)
GAY BUT (vhite] 216 points: (1189303921 199999999) PAT Mia 0, (Blue 617 points ind) 618 points (1 29503002 2.00000000)
{Green} 711 points’ (1.99957184, 2 00000000) HEP, felon] 709 points (1 9034541 2 20000000) DAT Teowio (Cyan) 73 poiate' ( 33211132 1 33339333)
CẮT BÙI [White] 214 polate: (1 88572418 1 33629491) XANG NGHẾP (Jelle) 70 po
IAT woe fed} E13 pounes’” (196677572 2.00000000) DAT RONG Green} 711 roi DAT Miko (Blue) 617 points (1 33896541 2 00000000) DAT THONG (Cyan), 73 points’ (1 $9998988 200000000) CẮT BỞI [White] 214 poiate: (1 99940412 1 99998002) DAT TONG {Cyan} 73 points: MAT HvGc’ (Red) 13 points: (1, 99966850 2, 00000000)
DAT HONG (Green) 731 points (1 33333398 2.00000000) DAT Miko (Blue) 617 pointe, (1 33711197 1 33983388) DAT WONG NGHỆP (Yellow) 709 points (1 99388988 2 00000000) GY BUT (White) 216 pointe: (189869003 200000000)
BUT (White) 214 points WUGC filed} 613 points: (1, 99997201 2 00000000) Dit mmc freen) 721 poinee’ (1 353033) 1 3333330) DAT MHA © (Blue) 617 pointe (1 68972018 1 34623481) DAT WONG NGHIEP (Yel1ov) 709 poiate (1.39840412 1 33338002) DAT TRENG [C/an) 73 polate: (1 93509004 2.00000000)
[Pair Separation (least to seo)
PAT Mia 0 {Blve] 617 points sod CAY BUT (White) 214 poines — 1 0057264 T MOQC [Red] 612 points and DAT HONG NGHIEP (Yellov] 709 pointe 1
Trang 21
- Bước 4 Phân loại vùng mẫu|
Để tiền hành phân loại với các mẫu đã chọn ta làm như sau:
Trên thanh menu chính của ENVI chon Classification\Supervised va chon phương pháp phân loại phù hợp Ở đây ta chọn Phương pháp Maximum
Likelihood
‘Neural Net
‘Support Vector Machine
Xuất hiện hộp thoai Classification Input File
© ROI Mask (Layer Band 1:1C8126087
8 Meide
Chon ảnh cân phân loại, chọn 1 band bất kỳ ở Select Mask Band\OK Xuất hiện hộp thoại Maximum Likelihood Parameters trong phan Select
Classes from Regions ta chỉ chọn các lớp cĩ trong chọn mẫu phân loại sau đĩ
chọn Choose, chỉ ra đường dẫn lưu kết quả\OK đề thực hiện phân loại ảnh
Trang 22
Kết quả ảnh phân loại của 2 thời kỳ như sau:|
#1 Max Lke (honglinh2003)anhpha - © MER B x
File Overy Enhance Tools Window Fe Overly Enhance Tools Window
Ảnh phân loại Hồng Lĩnh năm 2003 Ảnh phân loại Hồng Lĩnh năm 2013 Bước 5 Lập ma trận sai sơ phân loại
Trên thanh menu chính của Envi, Classification\Post Classification\Confusion Matrix\Using ground Truth Image, chon anh vita phan loai\OK
Trang 23
c Classification Input File =
“Select Input Fie: 283
Xuất hiện hộp thoại € Confusion Matrix Parameters
Classification Input File\OK | Output Confusion Matrix in F7 Paxels I% Percent
Report Accuracy Assessment © Yes No Xuất hiện hộp thoại Confusion [s2 2= # Matrix Parameters\Chosse, chỉ Output Result to @ Fie CO Memory ra đường dẫn lưu kết quả\OK | cme ce JONKIEMTRA\matran2003—C—C—C——C—C~*S Bảng ma trận ea os Class Confusion Matrix - ra loverati Accuracy = (1261707126170) 100 0000% lappa Somtticiene = “178000 "` 3 3 1 i 3 xeosi ị : i Ệ 3 ị Ệ m8 Tey SERED ố rota ĐẠI BA ST 3 i 3 wat N8 on 3 3 Hee bài xi 3 bộ egies cold
Trang 24- Bước 6 Lọc nhiễu, làm mịn ảnh
Trên thanh menu chính của Envi chon Classification\ Post Classification\ ification Input File
Majority/Minority/Parameters, Xuất hiện hộp thoại Cla chọn ảnh đã phân loại\OK
Xuất hiện hộp thoại Majority/Minority/Parameters, chỉ chọn các lớp cĩ trong
chon mau sau dé chon Choose để đặt tên cho ảnh lọc nhiễu và lựa chọn vị trí lưu
ảnh, Click OK để tiến hành lọc nhiễu
Kết quả ảnh phân loại sau lọc nhiễu hiễn thị như sau:
.€ #1 Majoriy Analysis (anhphanioaizo - © ME File Overly Enhance Tools Window
DẠT RUNG [Green] 711 ponls T NHAO [Bue]
Select Alltems | Clear Al tems
h phân loại 2003 đã lọc nhiễu
Analysis Method (* Majorty Minority #2 Majority Analysis (anhphanioaizo - > EM
Trang 25- Bước 7 Chuyển đữ liệu änh sang định dạng vecter (.evf) Trên thanh menu chính của Envi chon:
Classification\Post Classification\ Classification to Vector ‘Assign Clas Colors Rule Cacifier Clas Statistics Change Detection Statistics Confusion Matic › ROC Curves › Generate Random Sample > Majority/Minoiy Analysic Chump Classes Sieve Clases Combine Classes Overay Classes Xuất hiện hộp thoại Raster to Vecter Input Band chon file anh cần chuyển sau đĩ bam OK,
Xuất hiện hộp thoại Raster To vecter Parameters, chọn các lớp của ảnh cĩ trong
chọn mẫu sau đĩ chọn Choose chỉ ra đường dẫn đề lưu kết quả\OK Q Raster to Vector Input Band [Descpton Maơt Anas FesitSatNov Ot 001223 2014)
Trang 26- Bước 8 Đỗi từ đuơi*.evf sang *.shp Øf@ƒƒƒW؃ƒ{Œ “2
Trên cửa số ảnh vector Available œ Optio
Vecter List chon File\ Export Active Open Vector File
Layer to Shapefile |
Xuất hiện hộp thoai Output EVF
Layer to Shapefile chon Choose dé dit
Export Layers to RO tên cho file Shapefile và chỉ ra đường dẫn để lưu kết quảOK © Output EVE Layer to Shapefile HE Chuyên định dạng từ EVF sang SHP [EVF Layer: RTV Gocnhiew2003) ey |Se: 358,308 bytes |Current Output Directory: C:\Program Files (x86)\ITT\IDL7O\sb\hook a ¥ |Records: 722 [19.936 nodes] IPtoj: UTM Zone 48 North Meters} |Datum: WGS-84 ltibute: Yes Enter Output Filename [ shp]| Choose ï Select Al Layers | Deselect AIL ==} a a Load Selected | Remove Selected | OK | Cancel |
Kết quả chuyển đổi của 2 ảnh vecter như sau:
@ VedorWindow#1:Cursor Query - 8 Vector Window #2: Cursor Query - Ð MEE
File Mode Edt Options Help
ve “ATV foci Jct.“ Voce2037
Trang 27PHẢN 5 XỬ LY SO LIEU VA BIEN TAP BAN DO TREN PHAN MEM
ARCGIS
- Bước 1 Khởi động phần mềm ArcGIS, ‘le [tion fotmans je ấ©esen Teer Medex Hĩp
Owes vex oe, 1] e@H»w
- Bước 2 Mở file ảnh bản đồ cần cần biên tập (đuơi *.shp)
Trên thanh cơng cụ chọn biểu tượng Add đata + hoặc thao tac File\Add data
theo hinh sau:
| |File Edit View Bookmarks Insert Selectic
Add Data From Resource Center
[J Page and Print Setup
Trang 28Hộp thoại Add đata xuất hiện chọn đường dẫn vào thư mục chứa ảnh sau đĩ
chon 2 file vecter ảnh Hồng Lĩnh năm 2003 và 2013 sau đĩ chọn để
thêm đữ liệu vào Arcgis
Sau khi add dữ liệu sẽ hiển thị như sau:
Trang 29Bước 3: Gộp tên và Gộp lớp trong bảng thuộc tính
Nhap chuột phải vào từng file vecter 2003 và 2013 rồi chon Open
Attribute_Table
Label Features
Convert Labels to Annotation
@ Convert Features to Graphics
Convert Symbology to Representation
Data »
Bảng thuộc tính xuất hiện, để them trường đữ liệu Gộp lớp và Gộp tên ta chọn:
Option\Add field
m Attributes of 2003 5) switch Selection mm
‘Shape ‘Class_Name_ Class_1d | Part () Select Alt ˆ
Poygon | MATNUOC [Red] 374 pons Poygon | MAT NUOC [Red] 374 pons Potygon Poygon Polygon
WAT NUOC [Red] 374 ponts
WAT NUOC [Red] 374 ponts
WAT NUOC [Red] 374 ponts
Polygon | WAT NUOC [Red] 374 ponts Polygon | WAT NUOC [Red] 374 ponts Polygon | WAT NUOC [Red] 374 points Polygon | WAT NUOC [Red] 274 points Polygon | WAT NUOC [Red] 374 pons Polygon | WAT NUOC [Red] 374 ponts Polygon | WAT NUOC [Red] 374 ponts 72 | Polygon | WAT UOC [Red] 376 pomts Polygon Polygon Polygon Polygon Polygon Polygon Polygon
WAT NUOC [Red] 374 ponts WAT NUOG [Red] 374 pons WAT NUOG [Red] 374 ponts WAT NUOC [Red] 374 ponts WAT NUOC [Red] 374 ponts WAT NUOC [Red] 374 ponts
WAT NUOC [Red] 374 ponts rae
Trang 30Đặt tên và khai báo kiểu dữ liệu mới cho trường dữ liệu được thêm sau 46 chon OK
Trường đữ liệu mới đã được thêm, để gộp lớp và gộp tên vào trường đữ liệu mới
ta chon Option\Select by Attribute ‘cai Ỷ tšZãx:: xšz3 Ề 833533 šššSs š33813333ã533s Pompe WAT MVOC fe so Poy | WAT MIDE fe pn it eee i
Trong hộp thoại Select by Attributes chọn trường “Class_Name”= “DAT
NHA O ” chon Apply Dé chon tat cả các đối tượng mà trong trường
%Class Name” cĩ giá trị là “DAT NHA O ”
Trang 31
Select by Attributes 7 za
'DAT NHÀ O [Bue] 661 pons
[DAT NONG NGHIEP [Yelow] 1170 points’ "DAT RUNG [Green] 757 points’
[DAT TRONG [Cyen 250 pons Ce) & |
lề gộp lớp (gan lớp mới) cho các đối tượng được chọn ta thực hiện:
Trang 32Trong bảng chọn giá trị ta gõ số 1 và chon OK Field Calculator Fels: 3 Tre : Lee name © mmber [ia 6
S2, a © sero com EO lee somo orrex eo) BỘ v
Guy anh lớp mới sis) E[EIE
oe lớp mới [am [*|[-||=|
tực
được chon mm) BE
m
Để gộp tên (gán tên mới) cho các đối tượng đã được chọn ta thao tác tương tự:
Trên cột gopten của bảng Attributeb click chuột phải sau đĩ chọn lệnh Field Caculator trong bang dién gid tri ta nhap gia tri “DAT DO THI” Field Calculator ede: Type [x (Ciaes_name œ Ea C sưng Tản © Date "DAT OO THT
Như vậy, tất cả các đối tượng “DAT NHA O ” đã được chọn sẽ được gán
Trang 33Dé gộp các lớp khác ta làm tương tự, trong trường hợp gộp 2 lớp trở lên
thành 1 lớp mới ta thực hiện như sau:
Trong bang Attribute ta chon Option\Select by Attribute Trong hộp
thoai Select by Attributes nhap gia tri “Class Name”= “DAT RUNG” OR
“Class_Name”= “DAT NONG NGHIEP” chon Apply Dé chon tat ca các đối
tượng mà trong trường “Class_Name” cĩ giá tri la “DAT RUNG” hodc “DAT
NONG NGHIEP”
Để gộp lớp (gán lớp mới) cho 2 lớp này ta cũng thực hiện :
Trên cột goplop của bang Attributes click chuột phải sau đĩ chọn lệnh Field
Caculator để biên tập lệnh Trong bảng chọn giá trị ta gõ số 2 và chọn OK
Va để gộp tên (gán tên mới) cho các đối tượng đã được chọn của 2 lớp này ta thao tác tương tự: Trên cột gopten của bảng Attributes click chuột phải sau đĩ
chọn lệnh Field Caculator trong bảng điên giá trị ta nhập giá trị “DAT SX” Tương tự ta gộp 3 lớp cịn lại theo lệnh:
“Class_Name”= “DAT TRONG” OR “Class_Name”= “ MAT NUOC” OR
“Class_Name”=“CAY BUI” Dé chọn ra tất cả các đối tượng cĩ giá trị là “DAT
'TRONG"” hoặc “ MAT NUOC” hoặc “CAY BUI” trong truéng Class_Name Trong cột goplop nhập giá trị = 3 và cột gopten nhập giá tri = “DAT KHAC”
Như vậy sau khi gộp lớp xong sẽ cĩ 3 lớp nh sau: lp 1 âđDAT DO THI,
lớp 2 ® DAT SX va lop 3 © DAT KHAC m Attributes of 2003,
FiO | Shape Cua hame [cioss-ta | Parts | tenot | Area|Gontop|S0PTM]] ^- 208 | Poygen_[ DAT RUNG [Green] 757 ponte 180 1800) DATSX
209 | Poygen [OAT RUNG [Green] 757 pons 880| 18906) DaTSx
210 | Polygon | OAT RUNG [Green] 757 points 211 | Polygon | DAT RUNG [Green] 757 points 212 | Polygon | DAT RUNG [Green] 757 points 213 | Polygon | DAT RUNG [Green] 757 points 214 | Polygon | DAT RUNG [Green] 757 points 215 | Polygon | DAT RUNG [Green] 757 points [ 216 | Polygon | DAT RUNG [Green] 757 points 217 | Polygon | DAT RUNG [Green] 757 points 218 | Polygon | DAT RUNG [Green] 757 points nfl || ||] fl | || oo 5 2
219 | Polygon | DAT RUNG [Green] 757 points DATSx 220 | Polygon | DAT NHA 0 [Blue] 651 pons 1 DATOOT 221 | Polygon | DAT NHA 0 [Bue] 681 pons + DATOOT 222 | Polygon | DAT NHA 0 [Blue] 681 points 1Ị— taƑ %@| ĐATDOT 223 | Polygon | DAT NHA 0 [Blue] 651 points 1| — 180, 1809] DATOOT 224 | Polygon | DAT NHA O [Blue] 651 points 1Ị — 129, %6 DATOOT 225 | Polygon | OAT NHA O [Blue] 651 points 1Ị — %0| 2706 ĐATDOT 22 roypn | OAT A OB pon 1Ị — 120, 900 ĐATDOT 227 | Polygon | DAT NHA 0 [Blue] 651 points 11120] 800} paroorl| v
recor: «| <ff c9 |M[ show: [a selected | Record (out of 38 SORE
Các giá trị của 2 trường mới trong bảng thuộc tính
Trang 34Trong hộp thoại Ar Toolbox ta chọn lệnh Analysis Tools\Overlay\intersect B® Analysis Tools 4% Conversion “Toole
(@- Data Interoperability Tools
Data Management Tools
Geocoding Tools
ti ĐƯ Geostatistical Analyst Tools
[6-€@ Linear Referencing Tools Mobile Tools: Muttidimension Tools Network Analyst Tools ‘Schematics Tools Server Tools,
pf Gratial Onalvet Tonle ‹ >
Trong hộp thoại Intersect ở muc Input Features chon 2 file cin chéng xếp và chon Add > intersect mm Thư mực và tên dữ liệu mới được lưu ‘Tro oto) = | conf reroment >
Trong muc Output Feature Class ta chon muc luu va dat tén cho file chong xếp
chon save sau d6 chon OK
Trang 36
Ảnh chồng xếp sẽ chứa đữ
liệu của 2 bức ảnh Hơng Lĩnh năm
2003 và năm 2013, để thành lập
được bản đơ biến động và đưa ra diện tích biến động ta cần tạo thêm
1 trường dữ liệu mã chuyên đơi thể
hiện sự chuyển đơi diện tích giữa các loại hình sử dụng đất và điện tích chuyển đổi giữa các loại hình
sử dụng đất đĩ Trên cửa sơ Table
of contens nhấp chuột phải vào leyer ANHCHONGXEP chọn Open Attribute Table để mở bảng thuộc tính ZỞ ConvertEestures to Giaphics Convert Symbology to Representation Data › Save As Layer File (SF Propertis
Trong bảng thuộc tính ta chon
Options\Add Field dé them trường mới ste ty states Sa r8] = See reese sonen eetion y= Thơm mới † tường dữ hộu, —= BE mm [=4 ted te —= =
Trong hộp thoại Add Field chon ; ma tin như sau:
Trang 37Để tạo mã chuyên đồi ta thao tác như sau:
Trên cột MA_CD click chuột phải chọn Field Caeulator trong bảng Field Caculator nhập [goplop]&[GOPLOP'_1] sau đĩ chọn OK HS: s«án —— Ê SoftD&scending — Â| AdhancelSoing Summarize Z Statistics
Bẻ thống kê diện tích chuyên đổi, trên cột DT_CD trong bang
Attributes nhap chuét phai chon Calculate Geometry hộp thoại Calculate
Geometry xuất hiện, lựa chọn thành phần muốn thống kê, ở đây là diện tích nên
chọn là Area, sau đĩ chọn đơn vị và click OK ˆ Calculate Geometry 7a Sort Ascending
F Sort Descending Property: [EE ~]
4} Advanced Sorting Coordinate System
Summarize Use coordinate system ofthe datasource:
E sti [PCS GS 1984UTM Zone aN
Trang 38'Như vậy, đã tạo xong mã chuyển đổi và thơng kê được diện tích chuyển
lơi kết quả được hiền thị trong bảng thuộc tính như sau: Cues it | or [ tenons Are corcon | corres FRO] OBERT], = + mr =m 1 = T806 : at #46 3 [17100 : i = = es 00 at 00 5 = 1888 : = Tang S[ THỊ — 3êBh| 6744600 1|BATSANXUAT|— 31 — T88 %[ —TTỊ — SsêBh| 6744600 1[BATSANXUAT|-— sĩ 300 %[ —TTỊ sea | 6744600 1[BATSANXUA|— sr] 22500 %[ —TTỊ seas | 6744600 1[BATSANXUAT|-— sĩ 2700 SÍ TTIT[ — 348Bh| 8raseo6 1[BATSANXUAI|— š1— 18888 a[—a7 30280 | 700 1[BATSaNXUAI|— ãĩ S46 a] rene | ores 1[BATSaNXUAI|-— ãĩ 2h : TỊ — 6 2700 z|Barwua6 |—Z 2700 [ree | Brea T[oarsan aT] —or[ 20000 : TL MB| 388 >[Bareue = 600 |i i 2| — T889} 7966 3eae | 6744600 HE 3 [oar sanxva 31 = 3700 5400 i TỊ — 268|— S4 3|BaTieae = TI
~ Để thuận tiện cho việc tính tốn, ta Export dữ liệu sang bảng Exel, các bước thực hiện như sau:
- Trén hop céng cu Are Toolbox chon lénh Data Interoperoperability Tools\ Quick Export
~ Trong mục Input Layer của hộp thoại Quick
Export click chuột vào biểu tượng _s&i đĩ chọn
đường dẫn đến file ANH CHONG XEP.shp để chon file cần Export le lo ele le a! Record: 14] «ff T9 p[M| - show[AH Selected | Records (0 out of 1817 Selected) œ AreToolbox tạ ĐƯ 3D Analst Tools BY Analysis Tools 8 Cartography Tools f8- Conversion Tools -&% Data Interoperability Tools > AP Quick Import 4 Data Management Tools 2 Geocoding Tools
5 Geostatistical Analyst Tools A Linear Referencing Tools 2% Mobile Tools 8 Multidimension Tools Network Analyst Tools Ч Samples fa Schematics Tools f5- Server Tools
Trang 39Trong muc Output Dataset click vao biéu trong dé chon dinh dang,
Trang 40Bước §: Thành lập bản đồ
Sau khi xử lý xong số liệu và các bước chỉnh sửa cân thiết chúng ta tiến hành thành lập bản đồ trên phân mềm Arcgis
Để thành lập bản đơ, trước hết cần mở các file vecter dữ liệu sau đĩ thực hiện như sau:
Click chuột phải vao Layer trén cita sé Table Of Content, chon Properties
a, Bar
5H 2003 X Remove
m Open Attribute Table
o Joins and Relates › @ Zoom To Layer
@ Zoom To Make Visit
Visible Scale Range ›
Use Symbol Levels Selection › Label Features & Convert Features to Graphics Chỉnh sử thơng tin dữ liệu
- Hộp thoại Layer Layer Properties ox |
Properties xuat hién, chuyén #==ls=lseelss»[ESE]selxe=aelwelsetsee|reses| sang tab Symbology iw gir i an oe a wo
+ Tại cửa số Show chọn
Categories — Unique values
+ Cửa số Value Field chọn
trường Class name
+ Color Ramp cé thé chon tùy ý
- Sau khi chọn xong thì chọn vào Add All Values
Sau 46 chon Ok
Sinh vién: Lé Duy Nhat 40 MSSV: 1153076018