Ứng dụng gis trong quản lý tài nguyên và môi trường tổng hợp bài tập thực hành

BÀI TẬP THUC HANH 2: PHAN TÍCH MẠNG VÀ 3D Bài tập 1: Tao network dataset Tai Catalog, click chuot phai vao Transportation -> New -> Network dataset Requires the Network © Feature Cli

TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA THANH PHO HO CHI MINH

KHOA MOI TRUONG VA TAI NGUYEN

s

eo

UNG DUNG GIS TRONG QUAN LY TAI NGUYEN VA MOI TRUONG

TONG HOP BAI TAP THUC HANH

Giảng viên hướng dẫn: Ths Lưu Đình Hiệp

Họ và tên sinh viên: Huỳnh Tân Dũng MSSYV: 2110955

Thành phố Hồ Chí Minh — 2024

MỤC LỤC

BÀI TẬP THỰC HÀNH 1: CÁCH SỬ DỤNG PHẪN MÈM ENVI 2222 2222225 22c 3

BAI TAP THUC HANH 2: PHAN TICH MẠNG VẢ 3D -Q Q TH Hee 13

Bài tập I; TẠO HGMUOTĂE (ÏÍ(HSCÍ G001 TH HH TH ng re 13

Bài tập 2: Tìm đường ngắn nhất thông qua Network dafqset .-‹ l8

Bài tập 3: Tìm trạm cứu hỏa, bệnh viện 20

Bài tập 4: Vùng dich vu 24

Bài tập 5: Nội suy các đim quan trắc 37

Bài tập 6: Nội suy độ cao đỊA ÏHHÍU SG ch KH cá c0 4 0 37

BAI TAP THUC HANH 3: THUC HANH PHAN MEM ENVI VA VIEN THAM 46

Bài 1: Thực ảnh gộp ảnh và cắt ảnh Landsat 08 46

Bai 2: Phan lớp không giảm sắt (tự động phân logi) 50

801010580) 1108.G)/1.-) 1.4 50

Ì3j0111):13ð1)ì-318L`1):L:1+- RE aaaiaiadđdiiiỞỎỐỒ 51

Bài 3: Phân lớp có giữm sát 32

BAI THUC HANH 04: THUC HIEN PHẦN TÍCH KHÔNG GIAN S -cc- c2 59

BÀI TẬP THỰC HÀNH 1: CÁCH SỬ DỤNG PHAN MEM ENVI

Fle [dt Dingley Plocemarts Views Help

ôn 9 ©iÑMeeg2 “2 r6 LẦU 8 “ á Ko + Eheramere [Ge To

5ì @—+ b2) Ø fTO| @ Asus v @ ø @-) ole) eo} objec & = Hels

Boureeheeocs 22

eỞ)Ở

ane el

Ộ # + v'Kmageas [GTe 5] @ỞỞỞỞ đệ}] ụ fỉo| lạ Ux vụ ụ' 0-ỞỞỞ Gf9)ụ ẹ + Of Jo oe

Layer Manager -

4 & ỘBayer stacking input Fite XI ì =_

Select out Fle Fle rhomancr EỞỞ| /Raeer Meểsgenert/1asye Sang

x&

8

Bowe eReegs 2ouencus tr -leuege - = -

ệỞ+Ởcs) ụ @Ở+ỞỞỞ Ư3] ụ ft ẹ| @ Gating) ẹ & 0 -+Ở_ 90] 2 ẹ ec ce mialk

Bước 4: Vào Data Manaper đề tao band mau cho anh anhdem

Layer Manager M Teobox

lÑ L6 L22 129062 24

@ hse Layer Bard 1.008 L257 _1

Fle Edt Dinglay Placemarks Views Help

[Mê e6 £ £ c#r36 tít sue gee + vỗ Arddess [Ge To

Bước 6: Cắt ảnh khu vực nghiên cứu

Mé file shapefile BinhDuong.shp

Fle E42 Dagley Placemark Views Help

Boureeheeos 2eo ema ele Ameideres [Ge To

lẾ Lc09 L250 125052 20221213

LC08.LZ5P 120522027

{1055583 Le 107173 LUTM Zone 48 Noệt WGS-34

Bước 7: Mở Region of Interest (ROI) Tool -> Open -> Import Vector

i i

A reconts to 8 ergle FO!

7) Each mond to 2 separate RO!

© Record subset to 9 angle PO!

Select Record Subset Number of Selected Records [1/1 Number of Qượư: ROk Osos To @ New FO!

Fi Oaptay remit

Lat: 111766 N Lon: 1072655 91°E JTM Zone 48 Moth, WGS 84

Sow

Fle Edt Dinpley Placements Views Help

Boureheecs 2 210 148) ele SK D> Vetom= ol teretstore +

(®) select input File to Subset via ROI %

Select Input File: File Information

File: D:\ArcGIS\anhdem Dims: 7561 x 7721 x 7 [BSQ]

Size: [Unsigned Int] 817,298,734 bytes

File Type : ENVI Standard Sensor Type: Unknown Byte Order : Host (Inte!) Projection : UTM, Zone 48 North Pixel -: 30 Meters

Nemberof teen snlected 1 Select Aitoms Gear Al heme Mask pois outut of FOI 7| Yes Ms: Bay Vake 0

[Lat 1057150N.Ln 9722/80 UTM, Zone 48 Noth WGS-S

Buéc 9: Thé hién anh bang cach diéu chinh RGB Band: Mé ENVI Classic 5.3 >

File — Open — Chọn “catanhdem” —› Điều chỉnh các band RGB lần lượt như sau:

œQ@ Available Bands List =

File Options

ROI! Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B2.TIF):anhdem) ROI Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B3.TIF):anhdem) ROI! Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B4.TIF):anhdem) RO! Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B5.TIF):anhdem) RO! Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B6.TIF):anhdem) ROI Mask (Layer (Band 1:LC08_L2SP_125052_20221213_20221228_02_T1_SR_B7.TIF):anhdem) 4-8 Map Info

Dims | 2318 x 2343 (Unsigned Int) [BSQ]

Load Band _No Display v

1 Kết quả thực hiện

- Kết quả sau khi cắt ảnh:

ck & flee

Lư 1141 87M Lớn W649$ 9E UTM Zene 48 Noệt: WGS3⁄

- Kết quả sau khi chỉnh màu ảnh:

Tedbor fiegons of reevest Subset Data tom PO,

»

OGay Seae WAGE Cor

@p_ [ROI Mask (Layer and 11.008 12SP_128082 20201213 22122802

Aq [ROI Mack (ayer Gard 1LCO8_L25P 24021228 _02

UTM, Zone 48 Nosh WGS-84

Weoas xố Amadews [Ge To

Moamieny

BÀI TẬP THUC HANH 2: PHAN TÍCH MẠNG VÀ 3D

Bài tập 1: Tao network dataset

Tai Catalog, click chuot phai vao Transportation -> New -> Network dataset

Requires the Network

©) Feature Cli Analyst extension Disabled LÔ ơi

KR Relationsh if that extension is not >

Terrain checked in the Customize > ,

|) New Network Dataset x

Hinh 2 Chon Streets

é ĐAU 2AINGIVSE * = HWA:

3.8) New Network Dataset x

dt ‘Do you ware io mode! rns in ts network? On =- Over

È New Network Dataset x

‘The default connectivity settings for nework datasets establish connectivity only at

Coincident endpoints of line features curing the bud process Hyouwantto use diferent connectuty settings, click the Connectwity button below

‘You can change the connectivity settings now, of you can change them star the

network dataset has been created

Connectity

Hinh 4 Kiểm tra kết nỗi cua cde con dong tai Connectivity -> Next

New Network Dataset x

‘How would you like to model fhe elevation of your network features?

‘Using Z Coordinate Values from Geometry

OUsing Elavation Ficids

streets From Eno F_FLEV

Click in fhe Field column to set elevaton fields

<Back Cancel

Hinh 5 Chon Using Elevation Fields

ONG

Oves Traffic Tables:

Type

Time Slice Field Time Slice Start Time Time Slice Finish Time

Do you want to use Watfic data with this network dataset?

Speed

‘SpeedFactor_0000

12 AM

12AM Streets_DailyProfiles

Hinh 6 Chon Yes -> Next

"` SOHNE tert ous Lun seers pus

Usage Units Data Type Add

Remove Minutes Cost Double Remove All

RoacClass Descriptor Integer Rename

© TravelTime Cost Double -

'WeekendFallbackTravelTime Cost Double Tom

Thursday ProfilelD Field: PROFILE_5

Friday ProfilelD Field: PROFILE_6

Saturday ProfilelD Field: PROFILE_7

Hinh 9 Chon Finish

SWS toi setup > Witershed Delineator = HRU Analysis Vikite pur Taties - Edit SWAT Input ~ SWAT Simulation + ¿

Bài tập 2: Tìm đường ngăn nhat thong qua Network dataset

Tại thanh công cụ chọn Customize -> Extensions -> Network Analyst

File Edit View Bookmarks Insert Seiection Geoprocessing Customize Windows Help

Hinh I Bat Network Analyst

Chon New route

Hinh 2 New route

Chon Create Network Location Tool -> Chon cac tram dtrng -> Solve

Network Anayst» EER HD sreets ND ~ ae

Seated Cot (© Polygon Barriers (Restriction

2D Shoreline

3 Cioyent facliqy(TGE+73⁄2:- sities

Click vao Restriction 6 Network Analyst table -> Create Network Location

Tool -> Chọn các điêm bị câm không được đi qua -> Solve

Ngwok Asahzte ƑE)PBV, | sweets NO HN Steps (4)

O ceaptic Pk 23 ra TOs 8 tridenh

Point Barriers (1) TẾ]tvcone Lecenes

Deine Barriers

‘Added Cost (0) Được

Une (0

Restriction (0) addes Com

‘Scaled Cost (0) Rot

Polygon Barriers (0) te

40 Bi Une Baniers

‘Scales Cost (0) restriction

Scaled C Polygon Banners restriction scaled C

3 Seváce Area(S68AE7B.16 Facies

e

@ Locores

®Quwecnes

1 Point Bares trợ Restrict Added Cos

Bài tập 3: Từm trạm cứu hỏa, bệnh viện

File Et View Rockmart: ice Selection GesgGoeting Cuiomae Windows Hep

O#8ô+98xia= ®-(0854 -Ig)=zs8п nays Wie pt neu

+R ON OE RALS OR, [thd < 2|BAIØ B

Hình 1 Các điểm trạm cứu hóa

Ở Network Analyst -> New Closet Facility

= Gi Route:(08932E9A-AAED-4DE The closest facility solver

EI Stops finds the cost of traveling

Binet between incidents and

Hinh 2 New Closest Facility

Chuột phai vao Facilities -> Load Location -> Load From: Fire Stations

File Edit View Bookmarks Inset Selection Geoprocessing Customize Windows

O#8@+9Bxia+e-[ãe —— >8 @8n1+g lglc= wre put out i

AB oOunes #- ROS B As OR,

Eater-| > 1- #JEših < 2|8ã|B

Nghượrt Anant u |_ Giuse ox

Hinh 3 Tai cac trạm cứu hỏa

Chon vi tri xay ra tai nan -> Solve

20

‘Sealed Cost (0) ?fBunocated

— Pelygen Barriers (0) BB Point Barriers

Hình 5 Kết qua sau khi Solve

Chon biéu tugng Properties bén canh “Closest Facility 2”

Routes (1) 3 D\GIS TH2\DATA\SanFrancisco_v

& Point Barriers (0) & (3 Route:{08932E9A-AAED-4D8A-84

Restriction (0) & (3 Closest Facility:(70E17348-E5D3-4 Added Cost (0) 3 @ ees

© Line Barriers (0) @ Error

Restrictinn (0) @ | arated

Hình 6 Chọn nhiều trạm cứu hỏa

Doi Facility to find theo nhu cầu, ở day em chon 5

21

Gi Closest Facility 2(8043C9AE-7273-4-

GB DAGSS TH2\XULY SanFrancisoo\ Sant}

+3 %rwcc Aren:(56BAE781-7B£7-$AB4-

Gi Vehicie Routing Probiem‡542AB1BF

Hinh 8 Két qua sau khi chính lại số trạm cứu hỏa gần nhất

Chon Create Network Location Tool -> Chon 1 điểm Incidents -> Solve

Hình 9 Kết qua thu duoc sau khi Solve

22

Chon Select By Attributes Tai Layer chon Streets, tai bang Select nhập NAME =

“Van Ness Ave’ OR NAME = ‘Turk St’

1 (_] Only show selectable layers in this list Method: Create a new selection v

| OBJECTID

lịp

FEATTYP

FT F_JNCTID

| = <> Like "Van Buren St’

l "Van Dyke Ave'

> >= And "Van Fleet Ave'

"Van Keuran Ave' I

< <= Or "Van Ness Ave'

ob

Hinh 12 Két quả sau khi chỉnh lại s

dich vu tip 4: Vang Bai

->S

Chon 1 diém Facilities

24

sweese —_~ (a) [9 OB

© Feces (1) _— BO Sweets

s ao

‘Lines (0) BG Service Ares

© Line Barriers (0) ?QUniocated

Restriction Point Barriers

‘Scaled Cost (0) Error

© Added Cost

=Lines Line Barriers

8 OD Cost Matrix

1B & Vehicle Routing Problem

© Locstion-Allocation

2 Analysis inputs (0 Restricted Tums

Hinh 1 Két qua sau khi Solve

Click chuét phai o Polygons, chon Properties va sua Default Breaks thanh

Graphic Pick 210-1 @ D Closest Facility

Lines (0) s ice Ares

= Point Barriers (0) 8 Facilities

Restriction (0) @ Error

© Line Barriers (0) ?QUnocated

49) 2 Point Barriers

Cost (0) Q Error

© Polygon Barriers (0)

© Bee, Scaled ni B Unes

Network Analyst ax Contents ax =

Service Area xifm) BỊ

© Facilities (1) BO Streets

al BO Closest Facility

Graphic Pick2:5-10 || 5H Service Aves

© Graphic Pick 2:3-5 2 Facilities

Graphic Pick 2:1 -3 Error

Graphic Pick 2:0-1

Lines (0) 72?Quniocated

© Point Barriers (0) 4 63 Bont Bacriers

40) Q Error Added Cost (0) Restriction

© Line Barriers (0) Added Cost

EY Scaled Cost

8 OD Cost Matrix

6 Vehicle Routing Problem

© B Location-Allocation 1D Analysis inputs

BO Resticted Tums

BM Streets ND

1 &2 Basemap Layers

Hinh 3 Két qua thu duoc

Chuột phai vao Facilities -> Load Location -> Load From: CentralDeposts

M6é Layer Properties, d6i Default Breaks thanh 7 -> Apply

Lines (0) BD Cosest Facility

© Point Barriers (0) Sm Service Area

Restriction (0) Gri Faciities

© Polygon Barriers (0) Error

@ 0 OD Cost Matrix

@ 0 Venice Routing Problem

@ Dl Location-Allocation

© Ol Anaiysis inputs jeamalsted 3 -Facieee | GO Restricted Tums

26

Network Analyst + Table Of Contents ox

Service Area x ae81

‘Downtown Depot Ea

East Bay Depot

@Nortn Bay Depot 0 Streets

SENorth Bay Depot:0- 7) BI Service Area

- Depot: 0 LÁC

= lúne Barriers (0) Enor

5 Két qua thu duoc

Mé Layer Properties, déi Default Breaks thanh 2 4 6 -> Apply

@ East Bay Depot °

'North Bay Depot =O Streets

Polygons (9) BO Route

mEast Bay Depot: 4-6 | Closest Facility

Nowth Bay Depot: 4-€ 1 Facilities

North Bay Depot 2-4 ® ae

© Downtown Depot: 2 ®Qunecated

sm East Bay Depot: 0-2 © Point Barriers

™@ Downtown Depot: 0 aaa

Hinh 6 Két qua thu duoc

Ta xuất vùng vừa thực hiện thành một layer khác và đặt tên thành “Vungdichvu”

27

Hinh 6 Két qua thu duoc sau khi a@ export data

Click chuét phai tai bang chon Add Field -> Nhap MoTa, type chon Text -> OK

B-| Sat © oh x

Add Field Name:

Hình 7 Thêm trường dữ liệu

Chọn Select By Attributes trong Table -> Nhập như hình bên dưới -> Apply

28

._ Entera WHERE clause to selectrecords ¡n the table window

Method: Create a new selection ~ jer [ER] at, HE streets ND v RWS

ats *FaciliylD" ° ox

3 “Name”

= *FromBreak"

ID | *MoTa" Name | FromBreak | ToBreak |

—= 1 East Bay Depot : 4 - 6 4 6

East Bay Depot - 0-2 0 2

ay i=) LO Not Downtown Depot : 0-2 0 2

iit Is In Null GetUnique Values Go To:

MC SELECT * FROM Export_Outout WHERE:

9¢ *FromBreak" = 0 AND *"ToBreak" = 2| a

Hinh 8 Select by Attributes

Click chuột phải tại cột Mo Ta vừa tạo chon Field Calculator -> Nhập “Vùng dịch

Bs Type: Functions:

Shape ObjectID Ostnng FacilityID :

Dat Name ©pate FromBreak

0 » »! [B]*I @outof9 Selected)

Tình 10 Kết quả thu được

Chon Select By Location -> Chon Vungdichvu ở Source layer -> Apply

Select By Location x

relation to the features in the source layer

(Juse selected features (3 features selected) Spatial selection method for target layer feature(s):

(Apply 2 search distance

Hinh 12 Field Calculator

Bài tập 5: Nội suy các điêm quan trắc

Tao XY data cho cac diém quan trac như sau:

Choose a table from the map or browse for another table:

‘Specify the fields for the X, Y and Z coordinates:

Projected Coordinate System:

Name: TPHCM_VN2000 Geographic Name: GCS_VN_2000

v

(CO) Show Details Edit

@ War me if the resutting layer will have restricted

2 Ta được các điểm quan trắc

Trên thanh công cu chon Customize -> Extensions -> 3D Analyst

4 - SWAT Project Set

Extensions Select the ArcGIS Desktop

extensions you want to use

Extensions provide extended capabilities and

usually require that you have a license to use them

The dialog lists the extensions that are currently installed on your system and which work with the application you are currently using

Hinh 3 Bat 3D Analyst

Trén ArcToolBox -> Spatial Analyst Tools -> Interpolation -> IDW

32

Hình

tư © Serves Too's |

mm ^ © Soatial Aratys: Tous

s — * -Ÿ Spine win Bames,

‘Seach Radvs Sefteas A Topo te Rater

eee 2 J Á Nhohobeeetyfie

Hình 4 Nội suy IDIV

= Table Of Contents ax

Tại layer MauQuanTrac, click chuột phai chon Open Attribute Table -> Chon 1

diém bén ngoai nhw hinh:

33

H

ĐT TP CER Tp

og Faye Si

a bi Praga Co

ym Tm 8t

ra nag lùn

Ma Tang: gC Hư

p His, x8 TSn Th yi CO On Tp HOM

a eg hạ GD cw Hồ Gs Tp HO ONG Lon? 18 Tây An

He Mũ Lên 2 Xã Tân nộ GỌI Tạ H2 yi COCA Te HCM

Hình 7 Các điểm được chọn

Trén ArcToolBox -> Spatial Analyst Tools -> Interpolation -> IDW

34