Khi bạn add các dữ liệu lên khung dữ liệu thì một hộp thoại xuất hiện cho bạn biết một hay vài lớp dữ liệu thiếu thông tin tham chiếu không gian cho nên dữ liệu từ lớp này không thể tham chiếu được
Trang 1ĐẠI HỌC KINH TẾ QUỐC DÂN
KHOA KINH TẾ - QUẢN LÝ TÀI NGUYÊN, MÔI TRƯỜNG VÀ ĐÔ THỊ
************
BÀI TẬP LỚN MÔN HỌC GIS
Số:2
Đầu bài:THÔNG BÁO GIẢI PHÓNG MẶT BẰNG ĐỂ MỞ RỘNG ĐƯỜNG
Họ và tên sinh viên: Lương Thị Đức
Lớp: Kinh tế và Quản lý Môi trường 47
Tôi xin cam đoan tự mình làm các nội dung trong Bài tập lớn này,
không sao chép, gian lận dưới bất kỳ hình thức nào
Chữ ký của sinh viên
Phần chấm điểm
Câu Điểm Ghi chú
1
2
3
4
5
6
7
8
Tổng số
Trang 2BÀI LÀM
Khởi động ArcMap, mở A new emty map, Add các dữ liệu lên khung dữ liệu
Câu hỏi 1:
Khi bạn add các dữ liệu lên khung dữ liệu thì một hộp thoại xuất hiện cho bạn biết một hay vài lớp dữ liệu thiếu thông tin tham chiếu không gian cho nên dữ liệu từ lớp này không thể tham chiếu được.Chúng ta cần phải tham chiếu cho lớp dữ liệu này
Đóng ArcMap lại không Save những thay đổi Bây giờ sử dụng ArcToolbox để xác định hệ toạ độ cho
cả 2 lớp dữ liệu Làm việc này bằng cách sử dụng “ Define Projection Winzard( Shapefile, Geodatabase)”
Click đúp vào Define Projection Wizard ( Shapefiles, geodatabase) một hộp thoại xuất hiện, ta cần tìm đường dẫn đến bài tập lớn 2 này ( bằng cách click vào ) sau đó click vào Add, ta có hộp thoại sau:
1
Trang 3Click vào Next, tiếp theo một hộp thoại khác xuất hiện chọn Select Coordinate System, bấm vào Select chọn Projected Coordinate System, National Grids, Autralia, chọn AGD1984 AMG Zone 56.Click vào Apply, OK Một hộp thoại khác xuất hiện click vào Next, Finish
Trang 4Khi hệ toạ độ đã được xác định thì có các file có đuôi (*prj) được thêm vào cơ sở dữ liệu Dùng chương trình Notepad để mở file dcdb.prj ( bằng cách bấm chuột phải lên file cần mở chọn Open With….Notepad Nội dung của nó là:
PROJCS["AGD_1984_AMG_Zone_56",GEOGCS["GCS_Australian_1984",DATUM["D_Austral ian_1984",SPHEROID["Australian",6378160.0,298.25]],PRIMEM["Greenwich",0.0],UNIT["Deg ree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easti ng",500000.0],PARAMETER["False_Northing",10000000.0],PARAMETER["Central_Meridian" ,153.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT[" Meter",1.0]]
Câu hỏi 2:
Mở ArcMap, Add các dữ liệu lên khung dữ liệu
Vì đường Sir Fred Schonell Drive gồm nhiều đoạn phố ( street) cùng tên Vì vậy, chúng ta phải tiến hành lọc các đường này ra Bằng cách: Mở bảng thuộc tính của lớp sl_Street_amg, bấm vào Option chọn Select By Attributes và thực hiện lệnh “STREET_NAME”=’S IR FRED SCHONELL’ Click vào Apply,
Sau khi đã lọc xong các đường Sir Fred Schonell ta tiến hành thêm trường mới vào bảng dữ liệu ( Bấm vào Option, chọn Add Field…, đặt tên cho trường là Lenghth, chọn Type là Double.)Click OK
3
Trang 5Trước hết ta phải tìm đoạn mã Visual Basic tính độ dài.Vào Help => ArcGis Destop help, gõ Lenghth, vào ô tìm kiếm trong thẻ Index, chọn Caculating for lines, mở rộng How to make fiele caculations chọn Updating lenghth for a shapefile, copy đoạn mã trong đó ra:
Dim dblLength as double
Dim pCurve as ICurve
Set pCurve = [shape]
dblLength = pCurve.Length
Trong bảng thuộc tính của sl_street_mag, bấm chuột phải vào trường Lenghth chọn Caculates Values…, trong hộp thoại Field Caculator Click vào Advanced, trong hộp Pre-logic VBA script Code nhập đoạn
mã trên vào,
Trong hộp lenghth nhập dblLength Click OK.
Trang 6Sau đó Click chuột phải vào trường Lenghth chọn Statictis…xuất hiện một hộp thoại qua đó ta biết được
độ dài của con đường Sir fred schonell là: 1835,702498m.
Câu hỏi 3:
Xác định các mảnh đất nằm hoàn toàn hay một phần trong phạm vi 200m từ tâm đường Sir fred schonell drive Có bao nhiêu mảnh đất thuộc quyền sở hữu tư nhân có Tenure là FH được chọn? Tổng diện tích của những mảnh đất tư nhân được chọn này là bao nhiêu?
Đây là câu hỏi yêu cầu phải có truy vấn không gian, sử dụng công cụ Select by location để thực hiện truy vấn Vào Selection chọn Select by location:
Trong ô thứ nhất chọn “Select features from”
Trong ô thứ hai chọn lớp bạn muốn truy vấn, chọn “dcdb”
Trong ô thứ ba, nhập tiêu chí không gian mà truy vấn không gian phải tuân theo Trong trường hợp này, chọn “are within a distance of” (nằm hoàn toàn bên trong hoặc một phần)
5
Trang 7 Trong ô thứ 4, chọn lớp’sl_street_amg’ Đánh dấu vào ô Use Selected Features để đảm bảo rằng chỉ có các mảnh đất hoàn toàn hoặc một phần nằm trong phạm vi 200m tính từ tâm đường Sir fred schonell drive
Đánh dấu vào vùng đệm (Buffer) Apply a buffer to the features in sl_street_amg, đánh vào ô
Of là 200meters
Click Apply
Vậy tổng số các mảnh đất nằm hoàn toàn hoặc một phần trong phạm vi 200m từ tâm đường Sir fred schonell drive là 785 mảnh
Tiến hành xuất bảng thuộc tính của 785 mảng ra (bằng cách vào Option chọn Export….và add lên khung dữ liệu
Trang 8Từ bảng thuộc tính Export_Output_3 vào Select By Attributes, một hộp thoại xuất hiện thực hiện lệnh “TENURE”= ‘FH’.click Apply
Như vậy có 649 mảnh thuộc quyền sở hữu tư nhân
Sau khi lọc xong bảng dữ liệu, tiến hành xuất bảng thuộc tính này ra ( Vào Option chọn Export…) tạo thành bảng thuộc tính Export_Output_7
Trong bảng thuộc tính Export_Output_7 click chuột phải vào trường Area chọn Statictis
7
Trang 9Vậy tổng diện tích của những mảnh đất thuộc quyền sở hữu tư nhân là: 1080992,593180m2.
Câu hỏi 4: Tạo một bản đồ chỉ bao gồm những mảnh đất tư nhân được chọn trong câu 3 Không
chỉ có phần cắt ( Clip) từ vùng đệm.
Theo kết quả như câu 3, chúng ta chỉ quan tâm đến lớp Export_Output_7, vì vậy tắt tất cả các lớp khác
đi, chúng ta đã tạo ra được 1 bản đồ chỉ gồm những mảnh đất thuộc quyền sở hữu tư nhân