Đang tải... (xem toàn văn)
Trang 1 Bài 4: PHÂN TÍCH KHƠNG GIAN ĐỊA LÝ Advanced GIS operations: Xác định khu vực phù hợp để phát triển kinh tế dựa trên 1 số ràng buộc nhất định Link bộ dữ liệu: https://drive.google
IE402 – Hệ thống thông tin địa lý chiều Bài 4: PHÂN TÍCH KHƠNG GIAN ĐỊA LÝ Advanced GIS operations: Xác định khu vực phù hợp để phát triển kinh tế dựa số ràng buộc định Link liệu: https://drive.google.com/file/d/1QlS33NY_OVmVZEkxp4kqFeeejyampB7i/view?usp=sh aring Giới thiệu Lab: Phân tích lớp phủ có trọng số đa tiêu chí trình phân bổ khu vực dựa nhiều thuộc tính khác Nó thực tốt không gian raster cách sử dụng phương pháp dựa lưới (grid) Phân tích lớp phủ vector thực cách sử dụng cơng cụ xử lý địa lý, bị giới hạn câu trả lời nhị phân (phù hợp/không phù hợp) Phân tích lớp phủ raster cung cấp thứ hạng phù hợp cho phép kết hợp số lượng lớp đầu vào với trọng số khác Bài tập trình bày hướng dẫn quy trình làm việc để thực phân tích mức độ phù hợp địa điểm không gian raster Mục tiêu: Xác định khu vực phát triển kinh tế dựa yếu tố gần đường xá, tránh xa nguồn nước, khu vực bảo tồn, Bài 1: Xử lý liệu vùng liệu bị dư thừa gây mâu thuẫn (water polygon, water polyline) *Hướng dẫn bước thực hiện: Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Duyệt đến Assam.gpkg Bắt đầu cách mở Trình duyệt QGIS Nhấp vào nút "Mở" chọn tệp Assam.gpkg tải xuống Thêm lớp liệu Sau tệp mở, mở rộng để xem danh sách lớp liệu Kéo lớp vào khung đồ Bây giờ, bạn thấy năm lớp liệu bảng điều khiển lớp: boundary: Vùng biên giới Assam roads: Đường giao thông protected_regions: Các khu vực bảo vệ water_polygons: Các vùng nước dạng đa giác water_polylines: Các vùng nước dạng tuyến tính Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Bước phân tích lớp phủ chuyển đổi lớp liệu thành raster Điều giúp dễ dàng so sánh lớp liệu tính tốn giá trị kết hợp Một điều quan trọng cần cân nhắc tất raster phải có phạm vi Điều đảm bảo so sánh khu vực tương đương Trong ví dụ này, sử dụng lớp liệu "boundary" làm phạm vi cho tất raster Điều có nghĩa tất raster có kích thước độ phân giải Chúng ta sử dụng Rasterize (vector to raster) algorithm để thực thao tác Chọn lớp "roads" làm lớp Đầu vào để raster hóa đường giao thơng Nhập làm giá trị cố định A để ghi để pixel có đường giao thông gán giá trị Chọn đơn vị Geoferenced units làm đơn vị kích thước raster đầu để raster đầu có hệ tọa độ với lớp đầu vào Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Đặt độ phân giải raster đầu 15 mét để pixel đại diện cho khu vực rộng 15 mét Chọn ranh giới cho Phạm vi sử dụng để raster đầu có phạm vi với lớp đầu vào (boundary for Use extent) Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Cuộn xuống nhấp vào nút mũi tên phần Gán giá trị nút cụ thể cho dải đầu Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thơng tin địa lý chiều Sau raster hóa lớp "roads", chuyển đổi lớp vectơ lại thành raster cách sử dụng chức xử lý hàng loạt Nhấp chuột phải vào thuật toán "Rasterize (vector to raster)" chọn "Execute as Batch Process." Chọn tất lớp vectơ bạn muốn chuyển đổi thành raster Đặt tùy chọn raster hóa mong muốn Nhấp vào "OK" QGIS chuyển đổi tất lớp vectơ thành raster lúc Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Duyệt đến thư mục máy tính bạn đặt tên lớp raster_ Giao diện xử lý hàng loạt tự động điền tên với tên lớp điền vào tất hàng Đảm bảo hộp Tải lớp hoàn thành (load layers on completion (được chọn nhấp vào OK) Sau q trình xử lý kết thúc, bạn có lớp raster tải bảng điều khiển Layer Bạn nhận thấy có lớp liên quan đến nước - hai đại diện cho nước Chúng ta hợp chúng để có lớp đại diện cho vùng nước khu vực Sử dụng công cụ Raster calculator để xử lý vấn đề Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Tiến hành nhập biểu thức: "raster_water_polygons@1" + "raster_water_polylines@1" Nhằm hợp layer nước polygons polylines Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Raster hợp thu có pixel có giá trị cho tất khu vực có nước Nhưng bạn nhận thấy có số vùng có đa giác nước đường đa tuyến nước Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Những khu vực có pixel có giá trị - điều khơng xác Chúng ta sửa biểu thức đơn giản "raster_water_merged@1" > Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Kết quả: Bài 2: Proximity rasters - Euclidean distances Bây có lớp đại diện cho pixel đường nước, tạo trình quét vùng lân cận Chúng gọi khoảng cách Euclide - pixel raster đầu biểu thị khoảng cách đến pixel gần raster đầu vào Sau đó, raster kết sử dụng để xác định khu vực phù hợp nằm khoảng cách định từ đầu vào Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Sau trình xử lý kết thúc, lớp roads_proximitysẽ thêm vào bảng điều khiển Lớp Để hình dung rõ hơn, thay đổi kiểu dáng mặc định Nhấp vào nút Mở bảng điều khiển Kiểu lớp bảng điều khiển Lớp Thay đổi giá trị Max=5000 thành Color gradient Bài 3: Repeat the Proximity (Raster Distance) algorithm for the raster_water layer with same parameters and name the output water_proximity.tif Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Bài 4: Phân cấp điểm cho giá trị lân cận từ đường tới khu vực water (poly, line) Sau trình xử lý kết thúc, bạn áp dụng kiểu dáng tương tự trước để hình dung kết tốt Nếu bạn nhấp vào xung quanh raster kết quả, bạn thấy chuỗi giá trị liên tục từ đến 5000 Để sử dụng raster phân tích lớp phủ, trước tiên phải phân loại lại để tạo giá trị riêng biệt Gợi ý: 0-1000m → 100 1000-5000m → 50 >5000m → 10 Nhập biểu thức sau áp dụng tiêu chí vào đầu vào Nhấp vào nút … bên cạnh (các) lớp Tham chiếu chọn lớp ``roads_proximity` Đặt tên cho đầu roads_reclass.tif Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều nhấp vào Run Các Point có ý nghĩa tới khoảng cách cho điều kiện ban đầu tốn khu vực phát triển kinh tế 100*("roads_proximity@1"1000)*("roads_proximity@1"5000) Bài 5: Lặp lại trình phân loại lại cho water_proximity layer Ở sơ đồ ngược lại, pixel xa mặt nước có điểm cao 0-1000m → 10 1000 -5000m → 50 >5000m →> 100 Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web IE402 – Hệ thống thông tin địa lý chiều Bài 6: Thực phân tích lớp phủ cuối Hãy nhớ lại tiêu chí xác định phù hợp sau - gần đường, cách xa nguồn nước không nằm khu vực bảo vệ Gợi ý: ("roads_reclass@1" + "water_reclass@1")*("raster_protected_regions@1" )*"raster_boundary@1" Khoa Khoa học Kỹ thuật Thông tin Bộ môn Thiết bị di động Công nghệ Web !=