Nghiên cứu các giải pháp nâng cao độ chính xác của bình đồ ảnh tỷ lệ lớn từ ảnh vệ tinh độ phân giải cao

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT ĐỖ THỊ HOÀI NGHIÊN CỨU CÁC GIẢI PHÁP KỸ THUẬT NÂNG CAO ĐỘ CHÍNH XÁC CỦA BÌNH ĐỒ ẢNH TỶ LỆ LỚN TỪ ẢNH VỆ TINH ĐỘ PHÂN GIẢI CAO LUẬN ÁN TIẾN SĨ KỸ THUẬT Hà Nội, Năm 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT ĐỖ THỊ HOÀI NGHIÊN CỨU CÁC GIẢI PHÁP KỸ THUẬT NÂNG CAO ĐỘ CHÍNH XÁC CỦA BÌNH ĐỒ ẢNH TỶ LỆ LỚN TỪ ẢNH VỆ TINH ĐỘ PHÂN GIẢI CAO Ngành : Kỹ thuật Trắc địa - Bản đồ Mã số : 62.52.05.03 LUẬN ÁN TIẾN SĨ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TS TRẦN ĐÌNH TRÍ TS NGUYỄN XUÂN LÂM Hà Nội, Năm 2016 i LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng thân Toàn trình nghiên cứu tiến hành cách khoa học, số liệu, kết trình bày luận án xác, trung thực chưa công bố công trình Tác giả Đỗ Thị Hoài ii MỤC LỤC LỜI CAM ĐOAN i MỤC LỤC ii DANH MỤC CÁC CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG BIỂU vii DANH MỤC CÁC HÌNH VẼ viii DANH MỤC PHỤ LỤC x MỞ ĐẦU 1 Tính cấp thiết đề tài Mục đích nghiên cứu Đối tượng nghiên cứu Nội dung nghiên cứu Phương pháp nghiên cứu Ý nghĩa khoa học ý nghĩa thực tiễn Các luận điểm bảo vệ Những điểm luận án Cơ sở tài liệu viết luận án 10 Khối lượng kết cấu luận án 11 Lời cảm ơn CHƢƠNG TỔNG QUAN VỀ TÌNH HÌNH NGHIÊN CỨU 1.1 Tình hình nghiên cứu sử dụng ảnh vệ tinh độ phân giải cao để thành lập bình đồ ảnh phục vụ cho công tác thành lập chỉnh đồ 1.1.1 Trên giới 1.1.2 Ở Việt Nam 12 1.2 Kết luận 19 CHƢƠNG ĐẶC TÍNH KỸ THUẬT CỦA ẢNH VỆ TINH ĐỘ PHÂN GIẢI CAO 21 iii 2.1 Các đặc tính kỹ thuật vệ tinh độ phân giải cao 21 2.1.1 Các hệ thống vệ tinh độ phân giải cao 21 2.1.2 Đặc điểm quỹ đạo vệ tinh 24 2.1.3 Đặc điểm cảm biến 26 2.2 Một số đặc tính hình học ảnh vệ tinh độ phân giải cao 30 2.2.1 Đặc tính hình học 30 2.2.2 Sai số méo hình ảnh vệ tinh 33 2.2.3 Độ phân giải khả chiết tách thông tin từ ảnh vệ tinh 35 2.3 Thành lập đồ tỷ lệ lớn từ ảnh vệ tinh độ phân giải cao 38 2.3.1 Khả chiết tách thông tin từ ảnh vệ tinh độ phân giải cao 38 2.3.2 Yêu cầu độ xác ảnh vệ tinh độ phân giải cao cho việc thành lập đồ 42 2.4 Kết luận chương 45 CHƢƠNG BÌNH ĐỒ ẢNH VỆ TINH VÀ QUY TRÌNH THÀNH LẬP 46 3.1 Khái niệm bình đồ ảnh vệ tinh phương pháp thành lập bình đồ ảnh 46 3.1.1 Khái niệm bình đồ ảnh vệ tinh 46 3.1.2 Cơ sở khoa học phương pháp nắn ảnh 48 3.2 Quy trình thành lập bình đồ ảnh vệ tinh 49 3.3 Các phương pháp nắn ảnh vệ tinh 55 3.3.1 Phương pháp nắn theo mô hình vật lý 56 3.3.2 Phương pháp nắn ảnh theo hàm đa thức 60 3.4 Các yếu tố ảnh hưởng đến độ xác nắn ảnh vệ tinh 71 3.4.1 Ảnh hưởng độ phân giải ảnh vệ tinh 72 3.4.2 Ảnh hưởng mô hình toán học sử dụng để nắn ảnh 73 3.4.3 Ảnh hưởng điểm khống chế ảnh 73 iv 3.4.4 Ảnh hưởng mô hình số độ cao (DEM) 76 3.5 Trộn ảnh 78 3.5.1 Phương pháp biến đổi hệ màu RGB-IHS-RGB 78 3.5.2 Phương pháp phân tích thành phần - PCA 80 3.5.3 Phương pháp tổ hợp kênh 81 3.5.4 Phương pháp nhân ảnh (Multiplicative method) 82 3.5.5 Phương pháp trộn ảnh lọc tần số cao (HPF- High Pass Filter 82 3.6 Tiêu chuẩn đánh giá độ xác bình đồ ảnh 83 3.7 Kết luận chương 85 CHƢƠNG CÁC GIẢI PHÁP NÂNG CAO ĐỘ CHÍNH XÁC CỦA BÌNH ĐỒ TRỰC ẢNH TỶ LỆ LỚN TỪ ẢNH VỆ TINH ĐỘ PHÂN GIẢI CAO 86 4.1 Phân tích liệu khu vực thực nghiệm 86 4.1.1 Dữ liệu thực nghiệm 86 4.1.2 Đặc điểm địa lý khu vực thực nghiệm 90 4.2 Giải pháp số lượng, đồ hình bố trí điểm khống chế ảnh 94 4.2.1 Sơ đồ phương án bố trí điểm khống chế ảnh 94 4.2.2 Kết đánh giá độ xác 97 4.3 Lựa chọn mô hình nắn với trường hợp khối ảnh 100 4.4 Độ xác bình đồ ảnh từ ảnh vệ tinh độ phân giải cao sử dụng mô hình DEM 102 4.5 Lựa chọn phương pháp trộn ảnh ảnh độ phân giải cao 108 4.5 Xây dựng phần mềm đánh giá chất lượng bình đồ ảnh 114 4.5.1 Sơ đồ khối chương trình 116 4.5.2 Công cụ chương trình 117 4.5.3 Tính 117 4.5.4 Giao diện chương trình 117 v 4.5.5 Các bước thực chương trình 118 KẾT LUẬN VÀ KIẾN NGHỊ 123 DANH MỤC CÁC CÔNG TRÌNH KHOA HỌC CỦA TÁC GIẢ ĐÃ CÔNG BỐ CÓ LIÊN QUAN TỚI LUẬN ÁN 126 TÀI LIỆU THAM KHẢO PHỤ LỤC vi DANH MỤC CÁC CHỮ VIẾT TẮT Chữ viết tắt CCD Tiếng Anh Charge-couple device Tiếng Việt Bộ cảm biến Cơ sở liệu CSDL DBM Digital Building Model Mô hình nhà DEM Digital Elevation Model Mô hình số độ cao DLT Direct Linear Transformation Mô hình chuyển đổi tuyến tính trực tiếp DSM Digital Surface Model Mô hình số bề mặt DTM Digital Terrain Model Mô hình số địa hình GIS Geographic Information SysHệ thống thông tin địa lý tem GPS Global Positioning System Hệ thống định vị toàn cầu GSD Ground Sampling Distance Độ phân giải mặt đất HPF High Pass Filter Lọc tần số cao LiDAR Light Detection And Ranging Công nghệ đo Laser OLI Operational Land Imager Bộ thu nhận ảnh mặt đất PT Projective Transform Mô hình phép chiếu RFCs Rational Function Coefficients Các hệ số hàm hữu tỷ RFM Rational Function Model Mô hình hàm hữu tỷ Sai số trung phương SSTP TDI Time Delay and Integration Bộ cảm biến tích hợp làm trễ thời gian TIN Triangulated Irregular Network Cấu trúc mạng tam giác không mô hình số độ cao TIRS Thermal Infrared Sensor Bộ cảm biến hồng ngoại nhiệt vii DANH MỤC CÁC BẢNG BIỂU Bảng 2.1 Bảng tổng hợp đặc tính kỹ thuật ảnh vệ tinh độ phân giải cao22 Bảng 2.2 Nội dung thông tin đoán đọc ảnh IKONOS đối chiếu với yêu cầu thành lập đồ 1:5 000 theo [27] .39 Bảng 2.3 Nội dung thông tin nhận biết từ ảnh QuickBird .41 Bảng 3.1: Số lượng điểm không chế tối thiểu cho mô hình hàm đa thức 61 Bảng 3.2 Giá trị xê dịch vị trí điểm ảnh nắn vệ tinh ảnh hưởng sai số mô hình số độ cao 77 Bảng 4.1 Thông số ảnh WorldView - dùng thử nghiệm 87 Bảng 4.2 Các hệ số mô hình RPC 89 Bảng 4.3 SSTP vị trí điểm kiểm tra mô hình nắn chỉnh hình học (đơn vị mét) 98 Bảng 4 Chỉ số đánh giá chất lượng phổ ảnh sau trộn 113 Bảng 4.5 Chỉ số tương quan (Corr) ảnh trước sau trộn 114 viii DANH MỤC CÁC HÌNH VẼ Hình 1.1 Ảnh GeoEye - với điểm khống chế điểm kiểm tra 10 Hình 1.2 Ảnh Geoeye với điểm khống chế dùng để nắn ảnh 11 Hình 1.3 Bình đồ ảnh Quickbird [19] 14 Hình 1.4 Mối quan hệ hình học ảnh bề mặt trái đất 16 Hình 2.1 Quỹ đạo đồng Mặt trời .25 Hình 2.2 Phạm vi trạm thu ảnh vệ tinh 26 Hình 2.3 Cách bố trí phần tử CCD cảm biến thu chụp 27 Hình 2.4 Cách bố trí hàng CCD so le với cảm biến toàn sắc 28 Hình 2.5 Bộ cảm biến TDI chế độ quét theo hai hướng Bắc sang Nam với dòng quét nét đứt ngược lại với hướng từ Nam sang Bắc 30 Hình 2.6 Nguyên tắc quét ảnh dạng chổi quét 31 Hình 2.7 Nguyên tắc quét ảnh dạng máy quét chổi đẩy 33 Hình 2.8 Méo hình nguồn sai số cảm biến 34 Hình 2.9 Méo hình nguồn sai số cảm biến 34 Hình 2.10 Ảnh toàn sắc ảnh pan-sharpened QuickBird 42 Hình 3.1 Xê dịch vị trí điểm ảnh chênh cao địa hình 47 Hình 3.2 Sơ đồ quy trình công nghệ thành lập bình đồ ảnh vệ tinh 50 độ phân giải cao .50 Hình 3.3: Sai số xác định điểm khống chế ảnh .74 Hình 3.4: Sai số đoán nhận điểm ảnh ảnh vệ tinh 75 Hình 3.5 Xê dịch vị trí điểm ảnh ảnh hưởng chênh cao điểm địa hình mặt phẳng nắn ảnh, sai số DEM 76 Hình 4.1 Sự khác biệt độ phân giải kênh toàn sắc kênh đa phổ 91 Hình 4.2 Vị trí có chênh cao địa hình ảnh toàn sắc - WorldView2 .92 Hình 4.3 Phân bố khu vực địa hình ảnh WorldView-2 – Bắc Giang 93 Hình 4.4 Nội dung sơ liệu .95 - Sai số vị trí điểm khống chế sau bình sai theo mô hình hàm đa thức hữu tỷ qua phương án sử dụng số lượng điểm khống chế khác điểm (±m) 1.17 1.01 0.77 0.56 0.55 0.45 0.38 0.18 13 điểm (±m) 1.17 0.96 0.57 0.56 0.55 0.53 0.51 0.58 0.25 0.22 0.21 0.19 20 điểm (±m) 1.21 1.15 0.93 0.87 0.72 0.65 0.63 0.61 0.47 0.39 0.34 0.32 0.28 0.25 0.23 0.21 0.19 0.11 0.06 21 điểm (±m) 1.17 1.05 0.87 0.84 0.64 0.53 0.51 0.41 0.4 0.39 0.38 0.37 0.36 0.35 0.32 0.32 0.29 0.28 0.19 0.09 25 điểm (±m) 0.85 0.85 0.84 0.74 0.72 0.67 0.67 0.66 0.66 0.61 0.54 0.53 0.52 0.51 0.45 0.45 0.43 0.4 0.39 0.31 0.31 0.25 0.11 0.02 30 điểm 40 điểm (±m) (±m) 0.98 0.72 0.97 0.71 0.91 0.55 0.84 0.52 0.81 0.45 0.79 0.45 0.68 0.44 0.65 0.44 0.64 0.37 0.61 0.35 0.57 0.34 0.55 0.33 0.53 0.32 0.51 0.31 0.45 0.29 0.42 0.28 0.41 0.27 0.38 0.26 0.35 0.26 0.33 0.20 0.31 0.20 0.29 0.18 0.24 0.16 0.22 0.15 0.21 0.14 0.17 0.14 0.13 0.13 0.09 0.12 0.05 0.10 0.10 0.08 0.06 0.05 0.04 0.04 0.01 0.01 - Kết bình sai khối ảnh Quickbird bậc Standard áp dụng mô hình hệ số hàm hữu tỷ bậc + Sai số vị trí điểm khống chế điểm kiểm tra sau bình sai Sai số vị trí điểm sau bình sai khối ảnh (±pixel) Cảnh ảnh 6100 Cảnh ảnh 7010 Cảnh ảnh 8170 Điểm Điểm Điểm Điểm Điểm Điểm k.chế k.tra k.chế k.tra k.chế k.tra 0.77 1.11 1.12 1.21 1.05 1.32 0.65 1.02 1.03 1.08 0.91 1.07 0.49 0.89 0.98 1.05 0.82 1.04 0.46 0.71 0.84 0.97 0.78 0.97 0.43 0.53 0.67 0.89 0.61 0.94 0.41 0.23 0.58 0.85 0.44 0.82 0.38 0.11 0.49 0.68 0.32 0.53 0.21 0.05 0.41 0.52 0.18 0.32 0.20 0.37 0.15 0.17 0.31 0.12 0.15 0.18 0.10 0.14 0.12 0.09 0.11 0.11 0.07 0.09 0.10 0.05 Cảnh ảnh 8180 Điểm Điểm k.chế k.tra 1.15 1.02 1.07 1.08 0.99 1.02 0.91 0.97 0.89 0.79 0.85 0.62 0.75 0.21 0.45 0.19 0.27 0.25 0.19 0.13 0.08 0.02 + Sai số trung phương vị trí điểm khống chế điểm kiểm tra sau bình sai Cảnh ảnh 6100 Điểm Điểm k.chế k.tra 0.41 0.79 Cảnh ảnh 7010 Cảnh ảnh 8170 Điểm Điểm Điểm Điểm k.chế k.tra k.chế k.tra 0.65 0.83 0.55 Cảnh ảnh 8180 Điểm Điểm k.chế k.tra 0.93 0.72 0.81 + Sai số vị trí điểm ảnh nắn Số hiệu điểm Tọa độ gốc Loại Tọa độ ảnh nắn X (m) Y (m) điểm X (m) Y (m) Cảnh ảnh 6100 DT01 582190.4 KT01 583452.1 DT02 584761.9 DT15 587351.7 1150531 1150518 1150625 1144401 KC KT KC KC 582190.82 583452.83 584762.26 587351.59 1150531.28 1150518.76 1150625.83 1144400.11 dx 0.41 0.75 0.32 -0.15 dx 0.51 0.31 0.67 -0.45 dxy 0.65 0.81 0.74 0.47 DT03 KT02 DT04 KT03 DT06 DT19 DT07 KT04 KT07 DT08 KT09 DT18 KT10 DT05 KT11 DT11 DT12 Cảnh ảnh 7010 DT18 DT19 DT15 KT10 KT11 DT25 KT17 DT28 DT29 DT35 DT36 DT40 DT39 KT14 KT15 KT20 KT22 587698 586227.4 583555.1 582525.6 582535.5 585899.5 584936.1 587520.8 587438.6 587489.4 582652.9 582757.7 583713 586414.4 586645.8 584014.1 586436.3 582757.7 585899.5 587351.7 583713 586645.8 583950.8 582331.5 585679.1 587203.6 583807 586951.1 585421.1 582423.8 582172.9 587522.6 586971 582231 1150637 1150384 1149023 1148891 1147622 1143933 1147355 1149041 1146299 1147477 1145327 1143948 1144172 1149201 1144146 1145807 1145634 1143948 1143933 1144401 1144172 1144146 1142171 1141763 1140697 1141267 1139036 1138765 1137411 1137269 1142490 1142626 1140117 1139141 KT KT KC KT KC KC KC KT KT KC KT KC KT KC KT KC KC KC KC KC KT KT KC KT KC KC KC KC KC KC KT KT KT KT 587697.85 586227.52 583555.36 582525.81 582535.85 585899.79 584936.85 587521.47 587438.87 587489.00 582653.27 582757.92 583713.18 586413.94 586646.40 584014.95 586436.33 582757.7 585899.8 587351.7 583713.3 586645.8 583951.3 582331.5 585679.4 587203.6 583807.6 586951.1 585421.9 582423.8 582173.6 587522.6 586971.9 582231 1150637.65 1150384.62 1149023.89 1148891.75 1147621.31 1143933.37 1147355.23 1149041.24 1146299.79 1147476.18 1145327.11 1143947.92 1144172.53 1149201.23 1144146.10 1145807.07 1145634.59 1143948.1 1143933.3 1144400.9 1144172.4 1144146.1 1142170.6 1141763.4 1140696.2 1141267.1 1139036.9 1138765.2 1137411 1137269.6 1142489.9 1142626.3 1140117.8 1139140.9 -0.17 0.15 0.23 0.21 0.31 0.32 0.77 0.68 0.28 -0.42 0.35 0.19 0.22 -0.42 0.58 0.87 0.04 0.21 0.35 0.52 0.33 0.41 0.47 -0.2 0.31 0.91 0.68 -0.6 0.87 0.95 0.68 0.35 0.86 0.36 0.64 0.36 0.51 0.34 -0.51 0.31 0.65 0.36 0.36 -0.34 0.31 0.12 0.72 0.21 0.41 0.32 0.37 0.31 0.24 0.31 0.63 0.41 -0.1 0.36 -0.3 0.58 0.39 0.63 0.38 0.87 0.04 0.62 0.67 0.31 0.66 0.39 0.56 0.40 0.60 0.45 1.01 0.77 0.46 0.54 0.47 0.22 0.75 0.47 0.71 0.93 0.37 0.37 0.42 0.61 0.71 0.58 0.49 0.42 0.45 1.08 0.78 0.84 0.95 1.29 0.68 0.71 1.09 0.48 KT23 DT26 DT27 583877.8 1137211 KT 586762.3 1141936 KC 582425 1141010 KC 583877.1 586762.3 582424.4 1137211.3 1141935.3 1141009.6 Cảnh ảnh8170 DT08 DT09 DT10 DT13 DT14 DT15 DT16 DT17 DT23 DT30 DT29 DT32 KT05 KT06 KT07 KT08 KT19 KT15 KT16 KT18 587489.4 589816.7 592229 588560.1 591146.9 587351.7 589428.9 592020 588618 589906 587203.6 592343.1 588695.8 591118.4 587438.6 592214 591223.1 587522.6 592319.5 588859.2 1147477 1147423 1147474 1145798 1145890 1144401 1144527 1144579 1142657 1141179 1141267 1141608 1147468 1147705 1146299 1145883 1141194 1142626 1142866 1141438 KC KC KC KC KC KC KC KC KC KC KC KC KT KT KT KT KT KT KT KT 587489.7 589817 592229.5 588560.4 591147.1 587352.1 589428.6 592020.2 588618.6 589907 587204.2 592342.4 588696.1 591119.1 587438.8 592213.3 591223.4 587522.9 592320.4 588857.5 1147476.81 1147423.83 1147474.74 1145797.76 1145890.31 1144400.42 1144527.11 1144579.23 1142657.92 1141179.84 1141266.65 1141608.33 1147468.08 1147704.42 1146299.11 1145883.43 1141194.07 1142625.55 1142865.37 1141438.14 Cảnh ảnh 8180 DT30 DT29 DT32 DT33 DT34 DT38 DT36 DT37 DT42 DT44 589906 587203.6 592343.1 588424.6 591449.9 592199.2 586951.1 589707.7 587152.7 592088.6 1141179 1141267 1141608 1139771 1139991 1138678 1138765 1138890 1136039 1135969 KC KC KC KC KC KC KC KC KC KC 589905.34 587204.01 592343.36 588425.06 591450.67 592199.03 586950.21 589707.45 587152.89 592088.60 1141179.31 1141266.34 1141608.48 1139771.25 1139991.49 1138677.78 1138764.44 1138890.07 1136039.50 1135969.61 -0.6 0.36 -0.6 0.28 -0.2 0.08 0.70 0.42 0.62 0.31 0.29 0.27 0.43 0.56 0.42 0.35 -0.15 0.21 0.20 0.31 -0.14 -0.3 0.16 0.23 0.25 0.56 0.45 0.98 0.61 0.58 0.13 -0.7 0.25 0.32 -0.3 0.67 -0.5 0.21 -0.3 -0.7 0.06 0.32 -0 0.23 -0.1 0.87 -0.3 -1.8 -0 0.42 0.51 0.70 0.38 0.29 0.34 0.38 0.34 0.72 1.15 0.59 0.72 0.43 0.85 0.38 0.68 0.32 0.25 0.91 1.78 -0.65 0.38 0.26 0.45 0.78 -0.16 -0.89 -0.21 0.16 0.01 0.08 -0.18 0.40 0.48 0.27 -0.04 -0.09 -0.25 0.45 0.78 0.65 0.42 0.48 0.66 0.83 0.16 0.89 0.33 0.48 0.78 DT45 DT41 DT44 KT19 KT20 KT18 KT21 KT26 KT27 KT24 KT25 588714.2 591022.2 592088.6 591223.1 586971 588859.2 593018.2 588289.6 591008.8 587161.5 592137 1137532 1137532 1135969 1141194 1140117 1141438 1140052 1136141 1136138 1137291 1137323 KC KC KC KT KT KT KT KT KT KT KT 588714.62 591021.59 592088.97 591223.36 586971.49 588860.01 593017.82 588289.34 591007.35 587161.36 592137.73 1137531.53 1137531.60 1135967.38 1141193.81 1140117.04 1141438.14 1140052.77 1136140.25 1136137.35 1137291.00 1137322.20 0.41 -0.65 0.38 0.26 0.45 0.78 -0.38 -0.31 -1.45 -0.09 0.73 -0.38 -0.31 -1.45 -0.30 -0.07 -0.04 0.90 -0.68 -0.17 0.20 -0.30 0.56 0.72 1.50 0.40 0.46 0.78 0.98 0.75 1.46 0.22 0.79 Phụ lục 5: Kết sai số mặt phẳng điểm khống chế sử dụng DEM khác - Với DEM10 Y TT Hb1 Hb2 Hb3 Hb4 Hb5 Hb6 Hb7 Hb8 Hb9 Hb10 Hb11 Hb12 Hb13 Hb14 Hb15 Hb16 Hb17 Hb18 Hb19 Hb20 Hb21 Hb22 Hb23 Hb24 Hb25 Hb26 Hb27 Xg(m) 564,694.100 563,821.200 563,919.500 564,075.057 562,828.423 558,590.451 555,724.824 552,208.204 558,668.087 561,039.986 556,047.637 560,142.066 559,424.861 562,202.546 561,926.252 560,365.057 559,939.904 560,823.561 564,911.706 560,055.026 552,192.290 555,949.342 557,600.798 556,059.324 555,167.100 553,586.054 554,277.440 Yg(m) 2,315,764.800 2,313,346.600 2,309,526.500 2,306,368.995 2,303,172.362 2,308,386.508 2,308,652.886 2,302,306.735 2,305,832.329 2,305,282.999 2,311,387.582 2,312,111.193 2,314,941.397 2,312,851.516 2,306,816.764 2,307,546.237 2,309,807.427 2,310,399.864 2,302,751.414 2,304,853.911 2,308,097.122 2,309,334.924 2,315,161.598 2,305,248.811 2,308,051.089 2,308,680.371 2,315,849.226 Xn(m) 564,695.024 563,821.446 563,919.515 564,075.967 562,828.532 558,591.466 555,724.571 552,208.882 558,667.677 561,039.554 556,049.291 560,141.890 559,425.258 562,203.003 561,926.837 560,365.281 559,939.759 560,824.017 564,911.727 560,054.868 552,192.441 555,950.200 557,601.326 556,059.448 555,167.123 553,586.057 554,277.954 Yn(m) X (m) (m) 2,315,765.724 0.924 0.322 2,313,346.846 0.246 0.003 2,309,526.515 0.015 0.432 2,306,369.905 0.910 0.484 2,303,172.471 0.109 0.106 2,308,387.523 1.015 0.803 2,308,652.633 -0.253 0.321 2,302,307.413 0.678 -0.512 2,305,831.919 -0.410 0.986 2,305,282.567 -0.432 0.436 2,311,389.236 1.654 0.723 2,312,111.017 -0.176 -0.542 2,314,941.794 0.397 -0.368 2,312,851.973 0.457 0.109 2,306,817.349 0.585 0.459 2,307,546.461 0.224 -0.019 2,309,807.282 -0.145 -0.764 2,310,400.320 0.456 0.012 2,302,751.435 0.021 0.431 2,304,853.753 -0.158 0.003 2,308,097.273 0.151 -0.020 2,309,335.782 0.858 0.210 2,315,162.126 0.528 1.411 2,305,248.935 0.124 0.002 2,308,051.112 0.023 0.432 2,308,680.374 0.003 0.328 2,315,849.740 0.514 0.342 S (m) 0.978 0.246 0.432 1.031 0.152 1.294 0.409 0.850 1.068 0.614 1.805 0.570 0.541 0.470 0.744 0.225 0.778 0.456 0.432 0.158 0.152 0.883 1.507 0.124 0.433 0.328 0.617 - Với DEM5 TT Hb1 Hb2 Hb3 Hb4 Hb5 Hb6 Hb7 Hb8 Hb9 Hb10 Hb11 Hb12 Hb13 Hb14 Hb15 Hb16 Hb17 Hb18 Hb19 Hb20 Hb21 Hb22 Hb23 Hb24 Hb25 Hb26 Hb27 X (m) Y (m) S (m) 564,694.100 2,315,764.800 564,694.379 2,315,765.162 0.279 0.362 0.457 563,821.200 2,313,346.600 563,821.573 2,313,345.789 0.373 -0.811 0.893 563,919.500 2,309,526.500 563,920.151 2,309,526.730 0.651 0.230 0.690 564,075.057 2,306,368.995 564,075.185 2,306,369.507 0.128 0.512 0.528 562,828.423 2,303,172.362 562,828.802 2,303,172.523 0.379 0.161 0.412 558,590.451 2,308,386.508 558,590.988 2,308,386.718 0.537 0.210 0.577 555,724.824 2,308,652.886 555,725.202 2,308,653.232 0.378 0.346 0.512 552,208.204 2,302,306.735 552,208.258 2,302,306.164 0.054 -0.571 0.574 558,668.087 2,305,832.329 558,668.523 2,305,832.536 0.436 0.207 0.483 561,039.986 2,305,282.999 561,040.043 2,305,283.435 0.057 0.436 0.440 556,047.637 2,311,387.582 556,047.893 2,311,387.705 0.256 0.123 0.284 560,142.066 2,312,111.193 560,142.189 2,312,111.463 0.123 0.270 0.297 559,424.861 2,314,941.397 559,424.651 2,314,941.629 -0.210 0.232 0.313 562,202.546 2,312,851.516 562,202.648 2,312,851.870 0.102 0.354 0.368 561,926.252 2,306,816.764 561,926.793 2,306,816.933 0.541 0.169 0.567 560,365.057 2,307,546.237 560,365.303 2,307,546.608 0.246 0.371 0.445 559,939.904 2,309,807.427 559,940.121 2,309,807.850 0.217 0.423 0.475 560,823.561 2,310,399.864 560,823.578 2,310,400.209 0.017 0.345 0.345 564,911.706 2,302,751.414 564,911.867 2,302,751.173 0.161 -0.241 0.290 560,055.026 2,304,853.911 560,055.292 2,304,854.424 0.266 0.513 0.578 552,192.290 2,308,097.122 552,192.590 2,308,096.801 0.300 -0.321 0.439 555,949.342 2,309,334.924 555,949.800 2,309,335.134 0.458 0.210 0.504 557,600.798 2,315,161.598 557,600.944 2,315,162.009 0.146 0.411 0.436 556,059.324 2,305,248.811 556,059.390 2,305,248.936 0.066 0.125 0.141 555,167.100 2,308,051.089 555,167.201 2,308,051.303 0.101 0.214 0.237 553,586.054 2,308,680.371 553,586.466 2,308,680.586 0.412 0.215 0.465 554,277.440 2,315,849.226 554,277.781 2,315,849.328 0.341 0.102 0.356 Xg(m) Yg(m) Xn(m) Yn(m) - Với DEM2 TT Hb1 Hb2 Hb3 Hb4 Hb5 Hb6 Hb7 Hb8 Hb9 Hb10 Hb11 Hb12 Hb13 Hb14 Hb15 Hb16 Hb17 Hb18 Hb19 Hb20 Hb21 Hb22 Hb23 Hb24 Hb25 Hb26 Hb27 X (m) Y (m) S (m) 564,694.100 2,315,764.800 564,693.955 2,315,764.489 -0.145 0.311 0.343 563,821.200 2,313,346.600 563,821.656 2,313,346.480 0.456 0.120 0.472 563,919.500 2,309,526.500 563,919.521 2,309,526.285 0.021 0.215 0.216 564,075.057 2,306,368.995 564,074.899 2,306,369.393 -0.158 -0.398 0.428 562,828.423 2,303,172.362 562,828.574 2,303,172.310 0.151 0.052 0.160 558,590.451 2,308,386.508 558,590.683 2,308,386.385 0.232 0.123 0.263 555,724.824 2,308,652.886 555,725.178 2,308,652.287 0.354 0.599 0.696 552,208.204 2,302,306.735 552,208.373 2,302,306.957 0.169 -0.222 0.279 558,668.087 2,305,832.329 558,668.458 2,305,832.008 0.371 0.321 0.491 561,039.986 2,305,282.999 561,040.409 2,305,282.856 0.423 0.143 0.447 556,047.637 2,311,387.582 556,047.982 2,311,387.903 0.345 -0.321 0.471 560,142.066 2,312,111.193 560,142.522 2,312,110.747 0.456 0.446 0.638 559,424.861 2,314,941.397 559,425.293 2,314,941.562 0.432 -0.165 0.462 562,202.546 2,312,851.516 562,202.704 2,312,851.619 0.158 -0.103 0.189 561,926.252 2,306,816.764 561,926.404 2,306,817.180 0.152 -0.416 0.443 560,365.057 2,307,546.237 560,365.179 2,307,546.090 0.122 0.147 0.191 559,939.904 2,309,807.427 559,940.009 2,309,807.017 0.105 0.410 0.423 560,823.561 2,310,399.864 560,823.684 2,310,399.975 0.123 -0.111 0.166 564,911.706 2,302,751.414 564,911.713 2,302,751.676 0.007 -0.262 0.262 560,055.026 2,304,853.911 560,055.382 2,304,853.824 0.356 0.087 0.366 552,192.290 2,308,097.122 552,192.721 2,308,097.594 0.431 -0.472 0.639 555,949.342 2,309,334.924 555,949.499 2,309,334.917 0.157 0.007 0.157 557,600.798 2,315,161.598 557,601.032 2,315,161.715 0.234 -0.117 0.262 556,059.324 2,305,248.811 556,059.447 2,305,248.810 0.123 0.001 0.123 555,167.100 2,308,051.089 555,167.176 2,308,050.898 0.076 0.191 0.206 553,586.054 2,308,680.371 553,586.221 2,308,680.159 0.167 0.212 0.270 554,277.440 2,315,849.226 554,277.529 2,315,849.638 0.089 -0.412 0.422 Xg(m) Yg(m) Xn(m) Yn(m) - DEM TT Hb1 Hb2 Hb3 Hb4 Hb5 Hb6 Hb7 Hb8 Hb9 Hb10 Hb11 Hb12 Hb13 Hb14 Hb15 Hb16 Hb17 Hb18 Hb19 Hb20 Hb21 Hb22 Hb23 Hb24 Hb25 Hb26 Hb27 X (m) Y (m) S (m) 564,694.100 2,315,764.800 564,694.002 2,315,765.034 0.098 0.234 0.254 563,821.200 2,313,346.600 563,821.410 2,313,346.745 -0.210 0.145 0.255 563,919.500 2,309,526.500 563,919.506 2,309,526.945 -0.006 0.445 0.445 564,075.057 2,306,368.995 564,074.823 2,306,369.109 0.234 0.114 0.260 562,828.423 2,303,172.362 562,828.465 2,303,172.575 -0.042 0.213 0.217 558,590.451 2,308,386.508 558,590.106 2,308,386.841 0.345 0.333 0.479 555,724.824 2,308,652.886 555,725.145 2,308,653.131 -0.321 0.245 0.404 552,208.204 2,302,306.735 552,208.080 2,302,306.522 0.124 -0.213 0.246 558,668.087 2,305,832.329 558,668.351 2,305,832.741 -0.264 0.412 0.489 561,039.986 2,305,282.999 561,039.854 2,305,283.578 0.132 0.579 0.594 556,047.637 2,311,387.582 556,047.322 2,311,387.384 0.315 -0.198 0.372 560,142.066 2,312,111.193 560,141.984 2,312,111.327 0.082 0.134 0.157 559,424.861 2,314,941.397 559,424.896 2,314,941.464 -0.035 0.067 0.076 562,202.546 2,312,851.516 562,202.247 2,312,851.767 0.299 0.251 0.390 561,926.252 2,306,816.764 561,925.819 2,306,816.517 0.433 -0.247 0.498 560,365.057 2,307,546.237 560,364.955 2,307,546.558 0.102 0.321 0.337 559,939.904 2,309,807.427 559,940.154 2,309,807.760 -0.250 0.333 0.416 560,823.561 2,310,399.864 560,823.228 2,310,400.098 0.333 0.234 0.407 564,911.706 2,302,751.414 564,911.692 2,302,751.835 0.014 0.421 0.421 560,055.026 2,304,853.911 560,054.614 2,304,853.919 0.412 0.008 0.412 552,192.290 2,308,097.122 552,192.570 2,308,097.356 -0.280 0.234 0.365 555,949.342 2,309,334.924 555,949.275 2,309,335.141 0.067 0.217 0.227 557,600.798 2,315,161.598 557,600.504 2,315,161.892 0.294 0.294 0.416 556,059.324 2,305,248.811 556,059.323 2,305,248.937 0.001 0.126 0.126 555,167.100 2,308,051.089 555,167.153 2,308,051.494 -0.053 0.405 0.408 553,586.054 2,308,680.371 553,586.218 2,308,680.798 -0.164 0.427 0.457 554,277.440 2,315,849.226 554,277.228 2,315,848.916 0.212 -0.310 0.376 Xg(m) Yg(m) Xn(m) Yn(m) Phụ lục 6: Nguồn chƣơng trình “ Đánh giá độ xác bình đồ ảnh” //Su dung cac thu vien cua MDL MicroStation #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include ******************************************************************* #ifdef unix char sepchar = '/'; #else char sepchar = 92; #endif // [1] Open file dialog Private int File_getFile(char *filename, char *defaultExtension, char *prompt) { int stringLength; char dir[MAXDIRLENGTH], localSpec[MAXFILELENGTH]; //Thay doi vung mdlFile_getcwd (dir, MAXDIRLENGTH); stringLength = strlen (dir); if (dir[stringLength-1] != sepchar) { dir[stringLength] = sepchar; dir[stringLength+1] = 0; } //Mo dialog if (mdlDialog_fileOpen (localSpec,NULL, 0, NULL, defaultExtension, dir, prompt)) return SUCCESS; else strcpy (filename, localSpec); return SUCCESS; } #include #include #define DIALOGID_Noi #define HOOKDIALOGID_DialogBoxNoi #define HOOKITEMID_ToggleButtonNoi #define HOOKITEMID_PushButtonNoi #define PUSHBUTTONID_RunNoi #define PUSHBUTTONID_QuitNoi #define PUSHBUTTONID_NOI #define TEXTID_NOI #define RTYPE_File 'Noidlg' #define RSCID_File typedef struct MAIN_Params { char asciiFile[MAXFILELENGTH]; char DS[200]; }MAIN_Params; typedef struct Toado { double x; double y; }Toado; #define PUSHBUTTONID_OPEN #define PUSHBUTTONID_THUCHIEN #define PUSHBUTTONID_EXIT #define PUSHBTNID_Button 200000001 #define PUSHBUTTONID_THUCHIEN #define PUSHBUTTONID_EXIT #define PUSHBUTTONID_PHUNDIEM #define TEXTID_Text2 200000002 #define TEXTID_Text 200000001 #define RADIOBUTTONID_radio1 200000001 #define RADIOBUTTONID_radio11 200000002 #define RADIOBUTTONLISTID_200000003 200000003 #define TXT_DialogNoi "Danh gia chat luong binh anh" #define TXT_File "File toa khong che ngoai nghiep" #define TXT_pBtn_Select " " #define TXT_pBtn_Noi "~Thuc hien" #define TXT_pBtn_Cancel "~Xuat bao cao" #define TXT_SelectFile "Chon file " #define TXT_FileNotOpened "Khong the mo file " #define TXT_MoFile "Ban phai nhap file truoc" #define TXT_Noi " TXT " #define TXT_PBTN_THOAT "~Xuat bao cao" #define TXT_PBTN_THUCHIEN "~Thuc hien" #define TXT_PBTN_BUTTON "~Hien diem" #define TXT_DTXT_FALSECOORDINATEXY "False (X, Y)" #define TXT_DTXT_TEXT "X:" #define TXT_DTXT_TEXT2 "Y:" #define TXT_GBOX_GROUPBOX "XY False (m)" #define TXT_RBTN_RADIO12 "Nan anh" #define TXT_RBTN_RADIO1 "Tiep bien" #define TXT_DTXT_CHATLUONGBINHDOANH_FINAL "chatluongbinhdoanh_final\" #define TXT_DTXT_0 "0" *************************************************************