Trn Quc Bình, Trng HKHTN, HQGHN
I HC QUC GIA HÀ NI
TRNG I HC KHOA HC T NHIÊN

_____


_____
Trn Quc Bình

Bài ging

GIS ng dng

(dùng cho hc viên cao hc)


phangdoc
SS 
vi  là góc nghiêng ca hình đa giác (so vi mt phng nm
ngang).  bt đu, hãy chng minh cho trng hp tam giác
đn gin nh hình v bên.
2. Bt Spatial Analyst, đt Cell Size bng 20m và Working Directory tùy ý.
3. To DEM có tên là DEM_DS t lp DiemDoCao bng phng pháp ni suy
Spline Tension, Cell size 20m (xem bài 13 trong cun bài ging ArcGIS v phng pháp
ni suy).
4. T lp DEM  bc trên, lp raster đ dc có tên là Doc_DS, đn v dc là đ
(xem bài 13 trong cun bài ging ArcGIS v tính đ dc).
5. Dùng công c Raster Calculator to raster vi các cell có giá tr là 1/cos() vi  là
đ dc. Chú ý v đn v đo ca góc. Xem v Raster Calculator trong phn tr giúp ca
ArcGIS.
6. i tên lp kt qu thành InvCosGama.
7. Tính tng


n
i
i
1
cos
1

theo tng đn v hành chính bng công c Zonal Statistics As
Table. t kt qu đu ra là bng KetQua.dbf. Hãy da vào tham s Cell Size tr li din
tích ca các đn v hành chính là bao nhiêu?
Các din tích trên s ch là gn đúng vì khi s dng công c Zonal Statistics có
nhng pixel nm c bên trong và bên ngoài ranh gii hành chính.  tính din tích chính


3
Bài 2. Spatial Join, Topology trong chun hóa d liu
phc v đánh giá bin đng s dng đt

D liu đu vào: 2 file bn đ hin trng s dng đt nm 2005 và 2007  th trn
Gôi, huyn V Bn, tnh Nam nh  đnh dng DGN ca Microstation.
Yêu cu: chuyn đi d liu sang mt c s d liu geodatabase ca ArcGIS. Sau
đó kim tra và chun hóa d liu sao cho các vùng hin trng:
Có màu sc khp vi nhãn mc đích s dng;
Không đè lên nhau và không có các khong trng.
Các vn đ cn gii quyt:
- Chuyn d liu đu vào t đnh dng *.dgn thành geodatabase;
- Lc b các đi tng không cn thit trên bn đ hin trng s dng đt;
- Khc phc các li chng đè, vùng h ca các vùng hin trng s dng đt;
- Gán và kim tra thông tin v hin trng s dng đt ca tng khoanh đt;
- Hin th ni dung hin trng s dng đt.
Kin thc nn: công c Quick Import và Feature Class to Feature Class, Topology
(xem trình chiu trên lp và đc phn tr giúp Topology  [described, creating,
validating, errors, common tasks]). Cách to mt style đ trình bày nhanh các lp d liu
bn đ (xem hng dn s b trên lp và đc phn tr giúp Styles  [described, creating,
match to symbol in a style]). Các bc thc hin:
1. M các bn đ trong Microstation ri xem các thông tin cn thit (vùng hin
trng, nhãn hin trng) nm  nhng lp nào. Ghi thông tin này ra giy.
2. To mi mt Geodatabase và Feature Dataset có h ta đ VN-2000 trong múi 48.
3. i vi bn đ ca tng nm (Goi2005.dgn và Goi2007.dgn), làm các bc sau:
- Bng công c Feature Class to Feature Class, hãy Import các đi tng Polygon

theo th t tng dn ri đi chiu Fill vi Text_ theo chun
tô màu cho trong bng di (nhng đi tng có Distance > 0 thì cha cn đi chiu).
Nhng đi tng nào có màu và nhãn không khp? i vi nhng đi tng này hãy sa
màu theo nhãn cho khp.
Các loi đt có trong bn đ:
Màu Loi đt Màu Loi đt Màu Loi đt
5 LUC 53 ANI 72 DYT
12 BHK 56 SKC 75 DGD
26 RPT 58 SKX 81 DCH
36 TSN 60 DGT 85 RAC
42 ODT 63 DTL 88 TIN
45 DTS 66 DNT 89 NTD
52 QPH 69 DVH 91 SON

9. Gán MDSD cho các tha đt nh: Select
tt c nhng đi tng có Distance > 0, sau đó s
dng Field Calculator tính cho trng Text_ bng
lnh sau (Trong ArcGIS 10 cn b dòng DIM đi,
nu không s b li):
Select case [Fill]
case 5
HT = "luc"
case 36
HT = "tsn"
case 42
HT = "odt"
case 56
HT = "skc"
case 60
HT = "dgt"

Trn Quc Bình, Trng HKHTN, HQGHN

6
Bài 3. Thành lp bn đ bin đng s dng đt D liu đu vào: 2 lp d liu hin trng s dng đt th trn Gôi (HT_2005 và
HT_2007) đã đc chun hóa trong bài 2.
Yêu cu: Thành lp bn đ bin đng s dng đt giai đon 2005-2007 ca th trn
Gôi.
Kin thc nn: các công c Geoprocessing đã hc  bài 12 cun ArcGIS. Các bc thc hin:

1. Chng xp 2 lp HT_2005 và HT_2007 bng công c Union trong ArcToolBox,
đt tên lp đu ra là BienDong di dng Feature Class trong Geodatabase. Trong hp thoi
Union nh b không đánh du ô Gaps Allowed (nu không thì ArcToolBox có th s báo li).

2. Biên tp bin đng bng cách m bng thuc tính ca lp Biendong, to mt
trng mi có tên là LoaiBD dng Text(10) ri tính cho trng đó bng lnh nh sau
(Trong ArcGIS 10 cn b dòng DIM đi, nu không s b li):

if [Text_] = [Text1] then
BD=""
Else
BD =[Text_] & "->" & [Text1]
End if

LoaiBD = BD

2), phn C hãy t suy ngh to model đ t đng hóa quy trình thành lp bn đ bin đng
s dng đt (liên quan đn bài 3)

Phn A: thc hin theo hng dn chi tit  di đ t đng hóa quy trình chun hóa
d liu

1. T đng hóa các bc: chuyn đi d liu, xóa các đi tng không cn thit, to
quy tc topology cho các đi tng dng vùng
1.1. Chun b: trc khi làm bc này, trong th
mc C:\ArcGIS_Course\Caohoc\BiendongSDD hãy to
sn mt Geodatabase có tên là CH_model vi mt
Feature DataSet là VN2000, h ta đ VN2000 vi múi
chiu 48N.
1.2. To model: ArcCatalog, tìm nhánh
Toolbox.tbx  danh sách ri bm chut phi vào đó, chn
New  Model. Trên màn hình s xut hin ca s Model
Builder.
1.3. Thc thi chc nng chuyn đi d liu t
Microstation sang ArcGIS. Trong ca s Search, tìm
công c Feature Class to Feature Class ri kéo nó vào ca
s Model Builder. Nháy đúp vào công c này ri đt lp đu vào là Goi2005Polygon,
Output location là feature dataset VN2000 trong CSDL đã to, lp đu ra đt là Poly_05,
xóa ht các trng chuyn đi, ch gi li trng Level_ và Fill. Kt qu s đc nh hình
di bên phi.
Trn Quc Bình, Trng HKHTN, HQGHN

8

1.4. Chuyn feature class Poly_05 (đu ra ca bc trc) thành mt feature layer
ri chn các vùng  các level khác 30 (ta phi chuyn thành feature layer thì mi có th s

S xut hin ch P bên cnh đó. Bm chut phi vào lp Polygon mt ln na ri chn
Trn Quc Bình, Trng HKHTN, HQGHN

10
Rename và đi tên nó thành Input Polygon (tên này s hin ra trong hp thoi tham s khi
chy model)
Bm chut phi vào công c Feature Class to Feature Class ri chn Make Variable
 From Parameter  Output Feature Class, trên ca s Model Builder s xut hin thêm
1 ellips có tên là Output Feature. Bm chut phi vào nó ri chn Model Parameter. Bng
cách này chúng ta s có th đt tên cho lp d liu đu ra ca công c Feature Class to
Feature Class nh mt tham s.
Làm tng t nh vy vi công c Create Topology đ đt tên ca HT_05_topo nh
mt tham s ca model.
1.10. Thêm các lp đu ra vào ArcMap
Bm chut phi vào lp cui cùng ca model (có tên là HT_05_topo (5) trong hình
di, tên có th khác trong mô hình ca bn) ri chn Model Parameter và Add to Display
(bn phi đt đu ra là Parameter thì chc nng Add to Display mi có tác dng).
Bm chut phi vào lp đu ra ca công c Delete Features (Poly_05_layer(3)) ri
cng chn Model Parameter và Add to Display. Nu cn thit hãy đi tên nó thành mt
tên khác d nhn bit hn.

1.11. t tên, ghi và chy th model
Trong thc đn ca ca s Model Builder, vào Model  Model Properties ri đt
Name là Nhapdulieu (vit lin) và Label là Nhap du lieu.
Bm vào nút Save đ ghi model.
 chy model có 2 cách:
- Vào thc đn Model  Run Entire model đ chy model (trong trng hp này
bn không đt tham s đc).
Trn Quc Bình, Trng HKHTN, HQGHN


Các lp d liu trên đã đc đa sn vào bn đ có tên là VungGiaTri.mxd.
Yêu cu: Thành lp bn đ vùng giá tr đt đai theo giá th trng ca phng
Nguyn Du.
Kin thc nn: đc bài báo (đã lc bt) v thành lp bn đ vùng giá tr đt đai
theo giá th trng ca phng Láng Thng. Các bc thc hin:

1. M bn đ VungGiaTri.mxd

2. Tính giá quy chun cho các đim thu thp giá đt: to trng mi GiaQuyChuan
dng Float (hay Double) trong bng thuc tính ca lp GiaDat. Dùng Field Calculator tính
giá tr cho trng này nh sau:
Giá quy chun = giá thu thp đc / (h s v trí × h s hình th x 10
6
)

3. Ni suy giá quy chun: Bt và đt thông s cho Spatial Analyst: working
directory, Extent là Same as GiaDat, Cell size=1m.
Ni suy bng phng pháp IDW theo trng GiaQuyChuan. Các thông s ni suy t
la chn cho thích hp nht.

5. Tính giá quy chun cho tng tha đt: S dng công c Zonal Statistics đ tính
giá quy chun trung bình ca các tha đt theo giá tr ca các pixel t lp raster là kt qu
ni suy  bc 3. Nh đánh du ô Join output table trong hp thoi Zonal Statistics.
Giá quy chun trung bình ca các tha đt nm  trng nào ca bng kt qu ?

6. Tính giá d báo cho tng tha đt: to mt trng mi trong bng thuc tính ca
lp Thua vi tên là GiaTT, dng Double. Tính trng GiaTT theo công thc sau:

2. Thành lp bn đ vùng giá tr đt đai theo giá th trng
Khi thành lp bn đ vùng giá tr đt đai theo giá th trng, ngi ta phi tin hành
thu thp d liu v giá tr trng ca các tha đt nm trong khu vc nghiên cu. Rõ ràng
là tp hp d liu này không th ph kín toàn b khu vc cn thành lp bn đ bi ch có
nhng tha đt đc giao dch trên th trng mi có d liu v giá. Vn đ là trên c s
tp hp hu hn các d liu thu thp đc cn tính toán, d báo giá đt cho nhng tha
đt cha có d liu. Gi s giá đt là mt hàm s
G
ph thuc vào rt nhiu yu t
mn
vvvuuu , ,,,, ,
2121
, trong đó u là các yu t xác đnh (ví d nh nh hng ca v trí,
hình th, din tích, ca các tha đt), còn v là các yu t cha xác đnh (ví d nh nh
hng ca môi trng, dân trí, an ninh, ). Ta có th vit G mt cách gn đúng di dng
tích ca các hàm s nh sau:
), ,,(
), ,,()( )(
), ,,,, ,(
211
2111
2121
mn
mnn
mn
vvvGkk
vvvGufuf
vvvuuuGG



i
k (xem công thc 2). T giá quy chun tính
đc, tin hành ni suy thành 1 lp raster trong đó giá tr ca mi ô (cell) là
G tính cho
cell đó. Giá tr
G cho tha đt cn đnh giá, ví d nh tha s 19, đc coi là giá tr trung
bình ca các ô nm bên trong tha đt đó. Công vic này có th đc thc hin d dàng
bng các chc nng ca GIS, ví d nh công c Zonal Statistics ca phn mm ArcGIS. k
(19)
Líp néi suy
:k
(17)

:k
(21)
Zonal
Statistics
Líp gi¸ thùc
Líp gi¸ quy chuÈn
Trung
b×nh
17
21
19

đnh giá đnh k hàng nm da trên cn c là v trí và mc đích s dng. Vy có cách nào
đ cp nht d liu cho hàng nghìn tha đt mt cách nhanh chóng nht? Có th nhn thy
rng trong bng khung giá nhà nc thng ch có giá đt là thay đi, còn h thng phân
loi v trí ca các tha đt thng rt ít thay đi. Mc đích s dng ca các tha đt cng
không bin đng thng xuyên (hàng nm). Bi vy, gii pháp đa ra là không gán trc
tip khung giá nhà nc cho tng tha đt mà thit lp mt c s d liu khung giá riêng
ri liên kt nó vi c s d liu tha đt thông qua mã giá tr đc hình thành t mã
đng, mã mc đích s dng và mã v trí ca tha đt (hình 2), trong đó:
- Mã đng là s ca tuyn đng trong bng giá.
- Mã mc đích s dng đt (đt , đt sn xut kinh doanh phi nông nghip, ).
- Mã v trí ca tha đt gm có các v trí 1, 2, 3, 4 (v trí đc tính theo các tuyn
đng giao thông, đng ph, ngõ, ngách, ).  xác đnh nhanh mã v trí ca các tha
đt trên bn đ có th s dng các công c tìm kim không gian ca GIS, ví d nh Select
by Location ca phn mm ArcGIS.
Nh vy mi tha đt đu có 1 mã giá tr. Mi giá tr s có 1 loi giá tng ng
trong bng khung giá nhà nc hin hành. Bt c s thay đi nào trong khung giá thì giá
tr ca các tha đt trong c s d liu cng đc cp nht t đng. Bi vì mi quan h
gia c s d liu tha đt và c s d liu khung giá là Nhiu - Mt nên vic cp nht c
s d liu khung giá s nhanh hn nhiu so vi cp nht giá đt cho tng tha trong c s
d liu tha đt.

311 20 Chùa Láng VT1
312 12 Chùa Láng VT2
313 10.1 Chùa Láng VT3
314 9.1 Chùa Láng VT4
CSDL
thöa ®Êt
M· gi¸ trÞ
M· gi¸ trÞ = M· ®−êng & M· môc ®Ých sö dông ®Êt & M· vÞ trÝ thöa ®Êt
* 1

5. ê La Thành 511 512 513 514 521 522 523 524
6. Pháo đài Láng 611 612 613 614 621 622 623 624

4.2. Thành lp bn đ vùng giá tr đt đai theo giá th trng nm 2007
 có đc thông tin v giá th trng, các tác gi đã điu tra thu thp s liu t
nhiu ngun khác nhau nh: ch s dng đt, báo chí, công ty, vn phòng dch v nhà đt,
các trang web v bt đng sn. Các thông tin đc chia ra thành tng mc v nhà, đt, các
điu kin c s h tng, đ d x lý s liu. Sau khi lc nhng thông tin không tin cy
nh không phù hp vi giá trung bình trong khu vc, không có giao dch c th, thì còn
li 47 mu giá. Mt s mu giá này đc trình bày trong bng 2.
Bng 2. Mt s d liu thu thp đc v giá đt theo th trng  phng Láng Thng
Khu
vc
V trí
Din
tích (m
2
)
Hình th
Giá thu thp
(triu đng)
Ghi chú
Nguyn
Chí
Thanh
3 65 Vuông 30 Khu tp th
3 44 Vuông 40 Tha kim tra (TKT)
3 85 Vuông 35
Chùa
Láng

không ging nhau cho tng khu vc (ngay c trong phm vi mt phng). Trên đa bàn
phng Láng Thng, theo điu tra thì có th chia s b thành 5 vùng cho các h s nh
hng nh trong bng 3. Giá tr ca các h s đc xác đnh bng phng pháp so sánh
hay tham kho ý kin chuyên gia.
Bng 3. Phân vùng cho các h s nh hng
Vùng
Xung quanh
đng ph
H s v trí
VT
k H s hình th
HT
k
H s din tích
DT
k
(theo m
2
)
VT1 VT2 VT3 VT4
Vuông
vn
Hình
thang
a
giác
<40

40-
60

phng pháp ni suy, đng thi đ kim tra kt qu tính toán, 5 trong s 47 d liu thu
thp đc v giá đt th trng đã đc tách riêng (không tham gia ni suy) đ làm d liu
kim tra.
Các kt qu thu đc cho thy c 3 phng pháp ni suy cho kt qu gn tng
đng nhau: chênh lch gia giá tính đc và giá thu thp nm trong khong t 3 đn
22%. So vi 2 phng pháp khác, phng pháp ni suy tr trung bình trng s cho kt qu
tt hn mt chút. Tuy nhiên, cha th khng đnh đc phng pháp này v bn cht là tt
hn các phng pháp còn li hay ch đn gin là nó thích hp nht đi vi b d liu ca
phng Láng Thng. Các nghiên cu  nhng khu vc khác cn phi tip tc đc tin
hành nhm làm rõ vn đ này.

Kt qu thành lp bn đ vùng giá tr đt đai theo giá th trng nm 2007 t l
1:2000 cho phng Láng Thng đc th hin trên hình 3. ( đã lc bt )