Bài giảng GIS ứng dụng (Dành cho học viên cao học)

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)

Hà Ni - 2010
Trn Quc Bình, Trng HKHTN, HQGHN

1
Bài 1. Tính din tích mt dc

D liu đu vào: feature class DiemDoCao và DonViHanhChinh  geodatabase
C:\ArcGIS_course\CaoHoc\Matdoc\DoSon.mdb
Yêu cu: Hãy tính din tích t nhiên ca 3 xã ng Vn, Hoàng Xá và  Xá.
Kin thc nn: các công c Raster Calculator, Zonal Statistics, Attribute Join (Join
 Add Join tool)

Các bc thc hin:
1. Hãy chng minh mi quan h gia din tích ca mt
hình đa giác phng bt k vi din tích ca hình chiu ca nó
trên mt phng nm ngang:

cos
1

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

xác hn, hãy làm các bc tip theo đây.
8. Join lp DonViHanhChinh vi bng KetQua.dbf va to ra đ gán kt qu Zonal
Statistics vi các đn v hành chính. S dng trng Ten_DVHC đ join.
9. M bng thuc tính ca DonViHanhChinh, to mt trng mi là DT_doc (đnh
dng Double). Tính toán cho trng DT_doc mi to theo công thc nh hình v di:

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

2
pixel
phang
thuc
phang
pixel
doc
thuc
doc
DT
DT
DTDT 
Công thc tng ng trong Field Calculator nh sau:
[DonViHanhChinh.Shape_Area] * [Ketqua.SUM] / [Ketqua.COUNT]
Hoc:
[DonViHanhChinh.Shape_Area] * [Ketqua.DT_doc] / [Ketqua.AREA]
Kt qu s đc nh sau (kt qu ca bn có th hi khác đôi chút do các tham s
ni suy có th khác)

Hãy so sánh din tích thu đc  bc 9 so vi kt qu  bc 7.

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

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

thành lp Poly_xx (xx là nm hin trng) ch gi li 2 trng Level_ và Fill_
- Import các đi tng Annotation thành lp Anno_xx, ch gi li trng Level_ và
Text_
4. M 4 lp trong ArcMap, xóa tt c các đi tng không phi là vùng hin trng 
lp Poly_xx, và các đi tng không phi là nhãn hin trng trong lp Anno_xx (chú ý s
dng công c Select By Attributes và kích hot phiên biên tp Editing).
Chú ý: bc 5 và 6 di đây thc hin ln lt cho tng feature class Poly_xx
5. To mt Topology có 2 quy tc Must not overlap và Must not have gaps. t
Tolerance bng ~0.5m. S dng phn tr giúp ca ArcGIS, hãy gii thích 2 quy tc
topology nói trên là gì, ý ngha ca tolerance khi thit lp topology. Mt feature class có
th tham gia vào 2 topology đc không?
Trn Quc Bình, Trng HKHTN, HQGHN

4
6. M topology va to trong mt bn đ mi 
ArcMap. Kim tra và sa li topology cho toàn b lp
Poly_xx
1
.
7. (Spatial) Join các lp Poly_xx vi các lp Anno_xx
ca nm tng ng. 2 lp đu ra đt tên là HT_xx, các tham
s khác đt nh trong hp thoi  hình v bên. Hãy tr li:
- Spatial Join là gì ?
- Gii thích ý ngha ca các tham s (thông tin) trong
hp thoi  trên.
- Mc tiêu ca bc 7 này là gì ?
Xem thêm tr giúp ca ArcGIS v Spatial Join  các
phn Join  [Described, Spatial Join].
8. i vi 2 lp HT_xx va to, m bng thuc tính
ca chúng, kéo trng Text_ vào gn trng Fill, sp xp Fill

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"

case 63

1
Khi sa li topology trong ArcGIS 10.0, s báo li "Can not create feature". ây là bug trong bn 10.0 và
cách khc phc là cài đt SP1. Link v vn đ.  khc phc tm thi hãy dùng các công c v / biên tp đi
tng trong Editor, ví d nh Edit Vertices đ sa li Must not have gaps.
Trn Quc Bình, Trng HKHTN, HQGHN

5
HT = "dtl"
case 88
HT = "tin"
Case else
HT = "";
End Select

Text_=HT
Các tham s khác đt nh hình bên, chú ý đánh du ô Calculate Selected Records only.
Hãy gii thích cách thc gán mc đích s dng đt trong lnh tính toán trên.
10. To Style đ mô t Symbology cho các loi hình s dng đt. Khi to Style hãy
s dng quy đnh v màu sc RGB đi vi tng loi đt có trong cun "QUY NH K
THUT XÂY DNG BN  HIN TRNG S DNG T VÀ BN  QUY HOCH S DNG
T DNG S TRÊN PHN MM MICROSTATION"
(có th tìm cun này trong Google ri ti
v t Internet). Áp dng Style mi to cho c 2 lp HT_2005 và HT_2007 đ th hin ni
dung hin trng s dng đt lên màn hình.
11. t Alias là HTxx (xx là nm hin trng) cho trng Text_ ca c 2 lp HT_2005
và HT_2007 (đ có th phân bit chúng sau khi Union) bng cách bm chut phi vào tên
trng trong hp thoi Attribute Table ri chn Properties.

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

* Chú ý: th t và tên ca các trng Text_ và Text1 có th ngc li nu trong Union thay đi th
t các lp đu vào. Quan trng là trong trng LoaiBD phi th hin đúng th t ca hin trng các nm.
Hãy gii thích lnh tính toán trên dùng đ làm gì và làm nh th nào ?

3. Lc nhng bin đng gi (quá nh): đi vi nhng đi tng ca lp Biendong có
din tích nh hn 50 ri gán cho trng BD là "" (không bin đng). Hãy t suy ngh cách
làm.

4. Khái quát hóa và trình bày: chp các vùng bin đng cùng loi vi nhau, đt tên
lp kt qu là Biendong05_07. To Symbology và nhãn cho lp Biendong05_07. Hãy t
suy ngh cách thc hin bc này.

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

7
Bài 4. T đng hóa quy trình
phân tích không gian bng Model Builder


D liu đu vào: các d liu v hin trng s dng đt đã s dng  bài 2 (các file
Microstation hin trng s dng đt th trn Gôi các nm 2005 và 2007).
Yêu cu: T đng hóa mt s bc chuyn đi đnh dng d liu, chun hóa bn đ
hin trng s dng đt và thành lp bn đ bin đng s dng đt.
Kin thc nn: xem phn tr giúp ca ArcGIS v Model Builder  [Overview of,
Introduction to] và phn trình chiu minh ha trên lp.


Các bc thc hin: ni dung bài này gm 3 phn. Phn A, B hãy 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 (liên quan đn bài

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

dng công c Select By Attribute đc).
Hãy tìm công c Make Feature Layer bng ca s Search ri kéo nó vào ca s
Model Builder. Nháy đúp vào nó ri đt Input Feature là Poly_05, Output Layer là
Poly_05_Layer.
Tìm công c Select Layer by Attributes bng ca s Search ri kéo nó vào ca s
Model Builder, đt đu vào là Poly_05_Layer (biu tng màu xanh), bm vào nút SQL
đ đt lnh truy vn tìm các đi tng có [Level_] <> 30 (chú ý tên trng cho vào du
ngoc vuông).

1.5. Xóa các đi tng ngoài level 30 đã chn  bc trc: s dng công c Delete
Features. t Input Feature là lp kt qu ca bc trc.
1.6. To topology và xác đnh các feature class tham gia vào topology: thêm công
c Create Topology vào model, đt Input Feature Dataset là VN2000 đã to  bc chun
b, Output topology là HT_05_topo và Cluster Tolerance bng 0.5.
Thêm công c Add Feature Class to Topology vào model, đt Input Topology là
HT_05_topo  bc trên, Input Feature Class là lp đu ra ca công c Delete Features.
Trn Quc Bình, Trng HKHTN, HQGHN

9


1.7. Thêm 2 quy tc topology và kim tra topology:
Thêm công c Add Rule to Topology vào model, đt Input Topology là topology
đu ra ca bc trc, Rule Type là "Must not have gaps", Input Feature Class là Poly_05.
Tng t nh vy hãy thêm quy tc "Must not overlap".
Thêm công c Validate Topology vào model, đt Input Topology là lp đu ra ca
quy tc “Must not overlap” va làm  trên.

1.9. Thit lp tham s đu vào:
Bm chut phi vào lp Polygon (lp đu vào đu tiên) ri chn Model Parameter.

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

11
- Vào ArcCatalog ri nháy đúp chut vào model (Nhap du lieu), hãy đt các tham s
cn thit ri bm nút Run.
Hãy chy th vi d liu ca nm 2005, sau đó chy ln na vi d liu ca nm 2007.

Phn B: 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 cho các đi tng dng đim
Cách làm tng t nh trc nhng thay vì các lp Polygon ta cn các lp
Annotation, tên đu ra đt là Anno_05, các trng thuc tính cn chuyn đi là Level_ và
Text_. Bn cn xóa các đi tng không phi  Level 33.
i vi các đi tng dng đim, không cn thit lp Topology.

Phn C: To model đ t đng hóa quy trình thành lp bn đ bin đng s dng đt
Hãy to model đ t đng hóa các bc thành lp bn đ bin đng s dng đt 
bài 3.


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

12
Bài 5. Thành lp bn đ vùng giá tr đt đai

D liu đu vào: các d liu đu vào ca phng Nguyn Du, qun Hai Bà Trng,
TP. Hà Ni, nm  feature dataset VunggiatriDD trong geodatabase Bando_NguyenDu 
th mc VungGiaTri, trong đó có 3 feature class:
- GiaDat: lp đi tng dng đim cha thông tin thu thp đc v giá tr ca (mt
s) tha đt có giao dch trên th trng.
- Thua: lp đi tng dng vùng mô t các tha đt. Lp Thua đã có 2 trng thuc
tính mô t các h s v trí (HS_VT) và h s hình th (HS_Hinhthe).
- TenDT : tên đa danh, ch dùng đ trình bày bn đ.

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:

GiaTT = Giá quy chun × h s v trí × h s hình th
Trn Quc Bình, Trng HKHTN, HQGHN

13
Khi chy lnh tính toán có th xut hin báo li Null field do mt s tha không có
h s v trí (các tha không nm trong din có th giao dch đc, ví d đt c quan). Bm
OK đ chp nhn

7. Trình bày bn đ: Remove Join cho lp Thua. i vi nhng tha có
GiaTT=Null, hãy thay giá tr này (null) bng 0. To Symbology cho lp Thua theo cách
thc thích hp nht.
Trn Quc Bình, Trng HKHTN, HQGHN

14
NG DNG GIS TRONG THÀNH LP
BN  VÙNG GIÁ TR T AI
Trn Quc Bình, Lê Phng Thúy,  Th Minh Tâm
Khoa a lý, Trng i hc Khoa hc T nhiên, HQG Hà Ni

1. M đu
Trong bng xp hng tính minh bch ca th trng bt đng sn do Jones Lang
LaSalle - tp đoàn kinh doanh đa c và bt đng sn ln nht th gii - công b nm
2006, Vit Nam đc đng cui bng xp hng gm 56 nc có tham gia đánh giá [2].
n nm 2008, Vit Nam đã thoát khi v trí cui bng song mi ch đc nâng lên vài
bc. ánh giá ca Jones Lang LaSalle v th trng bt đng sn Vit Nam không phi là
duy nht, đa s các ý kin chuyên gia trong và ngoài nc đu khng đnh s kém minh
bch là mt trong nhng yu t gây cn tr ln nht ti s phát trin ca th trng bt
đng sn Vit Nam.  gii quyt vn đ này, mt trong nhng hng nghiên cu nhiu
trin vng là áp dng h thông tin đa lý (GIS) nhm đánh giá và công khai hóa giá tr bt
đng sn di dng các bn đ vùng giá tr đt đai. ( đã lc bt )

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




(2)
vi )(
iii
ufk  là các h s nh hng ca các yu t xác đnh. Các h s này s đc xác
đnh theo đánh giá chuyên gia hay bng cách so sánh nhng tha đt có các yu t
n
uuu , ,
21
gn ging nhau, ch khác nhau  yu t
i
u . Vn đ còn li cn gii quyt đ tính
giá G là xác đnh giá tr ca hàm
), ,,(
21 m
vvvGG  ca các yu t cha xác đnh. Trong
các phn tip theo,
G s đc gi là giá quy chun (đã đc tách khi các yu t nh
hng đã xác đnh).
Giá tr ca
G cho tng tha đt cn đnh giá có th đc xác đnh bng phng
pháp ni suy t nhng tha đt lân cn đã bit giá th trng. Cách thc ni suy đc mô
t trên hình 1.
Trong s đ trên hình 1, t giá tr th trng thu thp đc ca các tha đt s 17,
Trn Quc Bình, Trng HKHTN, HQGHN

15
21 (biu din trong lp giá tr thc) ta tính giá quy chun
G cho các tha đt này bng
cách ly giá thu thp đc chia cho các h s

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

Gi¸ trÞ tõn
g
cell

Hình 1. Ni suy giá đt.
T giá quy chun đc ni suy
G ca tha đt cn đnh giá, có th c tính đc
giá th trng ca nó bng cách nhân ngc vi các h s
i
k đã xác đnh cho tha đt này.
Các kt qu tính toán đc trình bày trên nn bn đ đa chính đ to ra bn đ vùng giá
tr đt đai theo giá th trng.
3. Thành lp bn đ vùng giá tr đt đai theo khung giá nhà nc
Nu nh giá đt theo th trng có tính khách quan cao và hay đc tham kho
trong các giao dch v đt đai thì khung giá nhà nc có vai trò rt quan trng trong qun
lý nhà nc v đt đai bi đây là c s đ giao đt, cho thuê đt, thu hi đt, tính thu s
dng đt, Theo Lut đt đai 2003, hàng nm U ban nhân dân các tnh, thành ph trc
thuc Trung ng công b khung giá đt cho phm vi đn v hành chính ca mình. Khung
giá nhà nc đc cho di dng bng nên khó tra cu và tính trc quan kém. Vic
chuyn hóa các thông tin này di dng bn đ s làm cho vic tra cu thông tin đc
thun tin hn và ngi s dng s có mt cái nhìn trc quan hn.
Xét v khía cnh k thut, vic chuyn đi thông tin t bng khung giá sang bn đ
là mt công vic không khó bng cách gán giá cho tng tha đt da theo bng khung giá
hin hành. Tuy nhiên,  đây chúng ta vn có th s dng GIS đ làm đn gin hóa công
vic này, và quan trng hn là t đng cp nht d liu khi mt bng khung giá mi đc
ban hành.
Trn Quc Bình, Trng HKHTN, HQGHN

16
Vn đ c bn đc đt ra là: đi vi khung giá nhà nc, các tha đt đc xác

đ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

Bng khung giá nhà nc

Hình 2. Liên kt CSDL tha đt vi bng khung giá nhà nc.
4. Th nghim và tho lun
( đã lc bt )
4.1. Thành lp bn đ vùng giá tr đt đai theo khung giá nhà nc nm 2007
C s đ thành lp bn đ vùng giá tr đt đai theo khung giá nhà nc là Quyt
đnh s 242/2006/Q-UBND ngày 29/12/2006 ca UBND Thành ph Hà Ni v vic ban
hành giá các loi đt trên đa bàn thành ph Hà Ni nm 2007 [1]. Theo quy tc đã trình
bày trong mc 3, mã giá tr ca các tha đt đc thit lp nh trong bng 1.
 gán mã giá tr cho các tha đt bng phng pháp bán t đng, công c tìm
kim không gian Select By Location ca phn mm ArcGIS đã đc s dng. Chng hn
nh tt c các tha đt có mã 211 (v trí 1 đng Nguyn Chí Thanh) đc xác đnh bng
lnh tìm kim nhng tha đt có chung cnh (share a line segment with) vi tuyn đng
Trn Quc Bình, Trng HKHTN, HQGHN

17
này. Nu mt tha đt thuc v 2 v trí khác nhau (ví d nh va mt ngõ, va mt ph)
thì nó s đc gán mã giá tr cao nht (mt ph). Bi vy, vic tìm kim các tha đt có
cùng mã giá tr đc thc hin t v trí thp nht (VT4) đn v trí cao nht (VT1). Kt qu
thành lp bn đ vùng giá tr đt đai theo khung giá nhà nc nm 2007 t l 1:2000 cho
phng Láng Thng đc th hin trên hình 3.
Bng 1. Mã giá tr ca các tha đt
ng ph
t  t SXKD phi nông nghip
VT 1 VT 2 VT 3 VT 4 VT 1 VT 2 VT 3 VT 4
1. Cu Giy 111 112 113 114 121 122 123 124
2. Nguyn Chí Thanh 211 212 213 214 221 222 223 224
3. Chùa Láng 311 312 313 314 321 322 323 324
4. Láng 411 412 413 414 421 422 423 424

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

1 46 Ch nht 100
2 50 Ch nht 80 Gn h Láng thng
2 50 Ch nht 76 Gn h (TKT)
ng
Láng
1 42 Vuông 80
1 64 Ch nht 76
2 55 Vuông 41
4 58 Ch nht 24
Pháo
ài
Láng
1 60 Ch nht 42
4 42 Ch nht 18 Tha kim tra
4 43 Thang 23
4 80 Vuông 15

Có 3 h s nh hng đn giá đt đc xác đnh trong quá trình th nghim là: h
s v trí
VT
k , h s hình th
HT
k
và h s din tích
DT
k
. Thang giá tr cho các h s này
Trn Quc Bình, Trng HKHTN, HQGHN

18

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

60-
80
>80
1 Nguyn Chí Thanh 1 0.5 0.34 0.26 1 0.98 - 0.88 1 0.9 0.88
2 ng Chùa Láng 1 0.43 0.34 0.2 1 - 0.88 0.99 1 0.95 0.85
3 ng Láng 1 0.51 0.39 0.3 1 0.96 0.92 0.93 1 0.95 0.74
4 Pháo ài Láng 1 - 0.71 0.57 1 0.96 - - 1 0.83 0.83
5 H Láng Thng 1 0.8 0.58 - 1 - 0.88 0.99 1 0.95 0.85

i vi các yu t nh hng ti giá đt cha xác đnh, cn phi ni suy cho các
tha đt mun đnh giá. Theo ý tng đã trình bày trong mc 2,  đây cn ni suy không
phi là giá đt thc (theo th trng)
G
mà là giá quy chun
G bng cách loi b các yu t
đã xác đnh. Trong th nghim ti phng Láng Thng,
G đc tính da theo công thc
(2) nh sau:
DTHTVT
kkk
G
G

 (3)
Sau khi ni suy, giá đt d báo ca các tha đt đc tính ngc li t công thc
(3):
DTHTVT
kkkGG  (4)
 ni suy, trong th nghim này chúng tôi đã s dng 3 phng pháp ni suy là:
tr trung bình trng s, ni suy Spline và ni suy Kriging. Nhm so sánh đ tin cy ca các

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 )