Ứng dụng địa thống kê trong phân tích và mô hình hóa dữ liệu dị thường độ cao vùng tây nguyên việt nam

i Bộ GIáO DụC Và ĐàO TạO TRƯờNG ĐạI HọC Mỏ - ĐịA CHấT VY TH THU TấN TI ỨNG DỤNG ĐỊA THỐNG KÊ TRONG PHÂN TÍCH VÀ MƠ HÌNH HỐ DỮ LIỆU DỊ THƢỜNG ĐỘ CAO VÙNG TÂY NGUYÊN - VIỆT NAM Ngành: Kỹ thuật trắc địa - đồ Mã số: 60520503 LUẬN VĂN THẠC SĨ KỸ THUẬT NGƢỜI HƢỚNG DẪN KHOA HỌC PGS.TS ĐẶNG NAM CHINH HÀ NỘI - 2015 ii LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi Tất q trình nghiên cứu; số liệu tính tốn; báo cáo kết nghiên cứu đƣợc trình bày luận văn trung thực, xác chƣa đƣợc cơng bố cơng trình khác Tơi xin chịu trách nhiệm nội dung đƣợc trình bày luận văn Hà Nội, tháng 10 năm 2015 Tác giả luận văn Vy Thị Thùy iii MỤC LỤC LỜI CAM ĐOAN i MỤC LỤC iii DANH MỤC CÁC BẢNG vi DANH MỤC CÁC HÌNH vii MỞ ĐẦU .1 Chƣơng 1: DỮ LIỆU ĐỊA KHÔNG GIAN VÀ LÝ THUYẾT ĐỊA THỐNG KÊ 1.1 Dữ liệu địa không gian (Geospatial Data) 1.1.1 Khái quát liệu quan trắc phân bố không gian 1.1.2 Hệ quy chiếu trắc địa liệu địa không gian 1.1.3 Vấn đề mơ hình hóa liệu địa khơng gian .12 1.2 Khái lƣợc lý thuyết địa thống kê 13 1.2.1 Tính quy luật tính ngẫu nhiên liệu địa không gian 13 1.2.2 Trƣờng ngẫu nhiên không gian phân loại 14 1.2.3 Phân tích sơ liệu địa không gian .15 1.2.4 Lập đồ địa thống kê cho liệu địa không gian 19 1.3 Một số đặc trƣng địa thống kê liệu địa không gian 21 1.3.1 Phƣơng sai, hiệp phƣơng sai bán phƣơng sai .21 1.3.2 Hiệp phƣơng sai thực nghiệm hàm hiệp phƣơng sai lý thuyết 22 1.3.3 Bán phƣơng sai thực nghiệm hàm bán phƣơng sai lý thuyết .26 1.3.4 Mối quan hệ hiệp phƣơng sai bán phƣơng sai .31 1.4 Thuật toán Kriging 32 1.4.1 Hệ phƣơng trình Kriging phân loại Kriging .32 1.4.2 Quy trình tính tốn nội suy Kriging 36 1.5 Thuật toán Collocation .38 1.5.1 Công thức Collocation 38 1.5.2 Quy trình tính tốn nội suy Collocation 42 iv Chƣơng 2: CÁC HỆ THỐNG ĐỘ CAO VÀ DỊ THƢỜNG ĐỘ CAO 45 2.1 Các hệ thống độ cao 45 2.1.1 Hệ thống độ cao 45 2.1.2 Hệ thống độ cao chuẩn 46 2.1.3 Hệ độ cao trắc địa 47 2.2 Dị thƣờng độ cao, phân loại phƣơng pháp xác định .48 2.2.1 Dị thƣờng độ cao phân loại 48 2.2.2 Phƣơng pháp xác định dị thƣờng độ cao 49 2.2.3 Geoid Quasigeoid, mối liên hệ chúng 51 2.2.4 Mơ hình trọng trƣờng Trái đất EGM2008 56 2.3 Phƣơng pháp đo cao GPS 58 2.3.1 Nguyên lý đo cao GPS 58 2.3.2 Vai trị mơ hình dị thƣờng độ cao (mơ hình Geoid) đo cao GPS .63 2.3.3 Triển vọng ứng dụng phƣơng pháp đo cao GPS 63 Chƣơng 3: DỮ LIỆU DỊ THƢỜNG ĐỘ CAO GPS-THỦY CHUẨN VÙNG TÂY NGUYÊN 65 3.1 Dữ liệu dị thƣờng độ cao GPS-thủy chuẩn vùng Tây nguyên 65 3.1.1 Giới thiệu chung khu vực Tây Nguyên 65 3.1.2 Dữ liệu dị thƣờng độ cao GPS-thủy chuẩn khu vực Tây Nguyên 66 3.1.3 Dị thƣờng độ cao EGM2008 vùng Tây Nguyên .69 3.1.4 Dữ liệu dị thƣờng độ cao phần dƣ vùng Tây Nguyên .70 3.2 Các đặc trƣng địa thống kê liệu dị thƣờng độ cao phần dƣ 71 3.2.1 Phƣơng sai, bán phƣơng sai hiệp phƣơng sai thực nghiệm 71 3.2.2 Các tham số hàm bán phƣơng sai lý thuyết (hàm cầu, hàm mũ) .73 3.2.3 Các tham số hàm hiệp phƣơng sai lý thuyết (hàm Markov bậc 3) 74 Chƣơng 4: SỬ DỤNG PHƢƠNG PHÁP KRIGING ĐỂ MƠ HÌNH HĨA DỮ LIỆU DỊ THƢỜNG ĐỘ CAO VÙNG TÂY NGUYÊN 77 4.1 Mơ hình hóa liệu dị thƣờng độ cao nội suy Kriging 77 v 4.1.1 Mơ hình hóa liệu dạng lƣới chuẩn (Grid) 77 4.1.2 Quy trình tính tốn nội suy xây dựng mơ hình dị thƣờng độ cao dựa liệu dị thƣờng độ cao phần dƣ 77 4.2 Khảo sát nội suy Kriging (OK SK) .83 4.2.1 Phân tích sơ liệu số dƣ dị thƣờng độ cao 83 4.2.2 Nội suy Kriging (SK) sử dụng hàm bán phƣơng sai cầu theo R thay đổi 87 4.3 Nội suy Collocation so sánh với kết nội suy Kriging 91 4.3.1 Nội suy collocation sử dụng hàm hiệp phƣơng sai Markov bậc 91 4.3.2 So sánh kết nội suy Kriging sử dụng hàm bán phƣơng sai cầu theo R thay đổi nội suy collocation 94 KẾT LUẬN VÀ KIẾN NGHỊ 97 TÀI LIỆU THAM KHẢO 99 PHỤ LỤC 100 vi DANH MỤC CÁC BẢNG Bảng 2.1 Thống kê liệu g (5’ 5’) mGal 57 Bảng 4.1 Bảng số dƣ dị tƣờng độ cao theo thống kê sơ .86 Bảng 4.2 Bảng kết tính bán phƣơng sai thực nghiệm .89 Bảng 4.3 Bảng kết hàm xấp xỉ cầu (trọng số nhau) 90 Bảng 4.4 Bảng số liệu quan trắc (với sai số đơn vị 004) 90 Bảng 4.5 Bảng tham số hàm bán phƣơng sai cầu 90 Bảng 4.6 Bảng kết tính hiệp phƣơng sai thực nghiệm (với dung sai khoảng cách 2.0km) .92 Bảng 4.7 Bảng kết xấp xỉ hàm Markov bậc (trọng số 93 Bảng 4.8 Bảng số liệu quan trắc (sai số đơn vị: 39.612) 93 Bảng 4.9 Bảng tham số hàm Markov bậc .93 Bảng 4.10 Bảng so sánh kết tính tốn nội suytheo hai phƣơng pháp Kriging collocation 95 vii DANH MỤC CÁC HÌNH Hình 2.1 Các hệ thống độ cao 48 Hình 2.2 Mơ hình Geoid EGM2008 – Việt Nam .58 Hình 3.1 Khu vực Tây Nguyên duyên hải Nam Trung Bộ ảnh chụp từ vệ tinh 66 Hình 3.2: Phạm vi phân bố điểm song trùng 68 Hình 3.3 Mơ hình Geoid EGM2008 vùng Tây Nguyên 70 Hình 4.1 Ngun lý vịng trịn chuyển động bán kính thay đổi .80 Hình 4.2 Tạo mơ hình Geoid cục (grid) từ mơ hình tồn cầu .81 Hình 4.3 Quy trình tính tốn xác hóa dị thƣờng độ cao 82 Hình 4.4 Biểu đồ tần suất thực nghiệm 87 Hình 4.5 Các điểm đƣợc chọn để nội suy Kriging 88 Hình 4.6 Đồ thị hàm bán phƣơng sai cầu 91 Hình 4.7 Đồ thị hàm hiệp phƣơng sai Markov bậc .94 Hình 4.8 Mơ hình EGM08TN vùng Tây Ngun sau xử lý địa thống kê phần dƣ dị thƣờng độ cao 96 MỞ ĐẦU Tính cấp thiết đề tài Trong thời đại công nghệ GNSS (Global Navigation Satellite System), vị trí điểm bề mặt Trái Đất đƣợc giải cách hiệu quả, nhƣng độ cao điểm lại chƣa thể giải trọn vẹn nhƣ khơng có thơng tin vị trí Geoid hay Quasigeoid hệ quy chiếu Trái Đất Chính vấn đề nghiên cứu trọng trƣờng xây dựng mơ hình Geoid/Quasigeoid phạm vi tồn cầu nhƣ cục ln nhiệm vụ quan trọng ngành Trắc địa nói riêng ngành khoa học Trái Đất nói chung Đồng thời áp dụng phƣơng pháp địa thống kê để phân tích xử lý số loại số liệu trắc địa nhƣ số liệu trọng lực, số liệu đo cao vệ tinh, dị thƣờng độ cao xây dựng mơ hình Geoid, mơ hình số địa hình DTM vv…Thuật tốn Kriging – cung cấp giải pháp nội suy tối ƣu, với ƣớc lƣợng tuyến tính tốt vị trí, sử dụng tới mơ hình bán phƣơng sai xác định vị trí cho liệu không gian Hiện nay, nhiều nhà khoa học tự nhiên hƣớng nghiên cứu vào lĩnh vực thống kê địa thống kê để tìm lời giải tối ƣu cho nhiều toán phức tạp khoa học Trái Đất Cùng với xu lựa chọn đề tài luận văn là: “Ứng dụng địa thống kê phân tích mơ hình hóa liệu dị thường độ cao vùng Tây Nguyên – Việt Nam” Mục đích nghiên cứu luận văn: Nghiên cứu phát triển ứng dụng phƣơng pháp địa thống kê xử lý số liệu phục vụ xây dựng mơ hình Geoid Quasigeoid khu vực Đối tƣợng phạm vi nghiên cứu: - Địa thống kê ứng dụng chúng phân tích, mơ hình hóa liệu dị thƣờng độ cao - Phƣơng pháp Kriging số phƣơng pháp nội suy khác - Nghiên cứu liệu dị thƣờng độ cao vùng Tây Nguyên – Việt Nam Nội dung nghiên cứu: - Tìm hiểu lý thuyết địa thống kê; - Dữ liệu dị thƣờng độ cao ứng dụng đo cao GPS; - Các đặc trƣng địa thống kê liệu dị thƣờng độ cao vùng Tây Nguyên; - Lập chƣơng trình để tính tốn thực nghiệm nội suy Kriging so sánh với kết tính theo phƣơng pháp Collocation Phƣơng Pháp nghiên cứu: - Phƣơng pháp thu thập tài liệu, số liệu: Thu thập tài liệu có; cập nhập thơng tin mạng Internet; tìm kiếm số liệu tọa độ GPS, độ cao thủy chuẩn có đủ độ xác tin cậy phục vụ cho nghiên cứu - Phƣơng pháp phân tích: Tìm hiểu lý thuyết địa thống kê; hàm phƣơng sai, bán phƣơng sai thuật toán nội suy Kriging, Collocation - Phƣơng pháp so sánh: So sánh ƣu nhƣợc điểm thuật toán phƣơng pháp sử dụng nghiên cứu, để tìm phƣơng án tối ƣu - Phƣơng pháp tổng hợp: Tập hợp kết nghiên cứu, tìm đƣợc thuật tốn xác hóa dị thƣờng độ cao Phƣơng pháp ứng dụng cơng nghệ tin học: Lập chƣơng trình tính tốn cho máy tính thực quy trình mơ hình hóa liệu dị thƣờng độ cao Ý nghĩa khoa học thực tiễn đề tài: Phân tích đặc trƣng thống kê liệu số dƣ dị thƣờng độ cao tính từ mơ hình Geoid EGM2008 dị thƣờng độ cao GPS-thủy chuẩn vùng Tây Nguyên Sử dụng thuật tốn Kriging để mơ hình hóa Geoid đƣợc làm khớp liệu Chính xác hóa mơ hình Geoid phục vụ cho cơng tác Trắc địa – Bản đồ nói chung đo cao GPS nói riêng Cấu trúc luận văn - Lời mở đầu - Chƣơng 1: Dữ liệu địa không gian lý thuyết địa thống kê - Chƣơng 2: Các hệ thống độ cao dị thƣờng độ cao - Chƣơng 3: Dữ liệu dị thƣờng độ cao GPS – Thủy chuẩn vùng Tây Nguyên - Chƣơng 4: Sử dụng phƣơng pháp Kriging để mơ hình hóa liệu dị thƣờng độ cao vùng Tây Nguyên - Kết luận kiến nghị Tơi xin cảm ơn hƣớng dẫn nhiệt tình PGS.TS Đặng Nam Chinh thầy cô mơn trắc địa cao cấp giúp đỡ tơi hồn thành luận văn Tôi xin chân thành cảm ơn! 131 III(XH-SL)8 14 20 21.4621 108 28 40.2240 1588074.078 875220.622 -2.488 132 III(SL-NT)8 14 21 34.2667 108 41 55.7768 1590685.626 899055.180 -2.648 133 III(CH-IAR)8 14 22 6.3756 107 32 6.9839 1589976.201 773423.090 -6.788 134 II(TT-KT)29 14 22 23.2678 108 45.6461 1591234.505 833937.797 -4.437 135 I(VL-HT)103 14 23 8.5935 109 26.1592 1594177.801 934128.100 -3.016 136 I(DN-BMT)59 14 24 28.2744 107 58 16.9928 1594902.508 820435.102 -5.144 137 III(CH-PQ)7 14 25 26.5055 107 47 21.5246 1596447.353 800764.282 -5.744 138 III(CH-IAR)4 14 27 14.6278 107 36 12.3887 1599536.841 780671.986 -6.616 139 III(AL-DT)5 14 28 46.2736 108 55 13.4839 1604380.484 922766.418 -3.266 140 II(TT-KT)24 14 29 31.0694 108 13 30.1361 1604586.716 847682.394 -3.866 141 III(XH-SL)3 14 31 3.8818 108 26 50.0359 1607792.455 871618.123 -2.886 142 III(CH-PQ)3 14 31 6.8084 107 44 7.6664 1606842.827 794828.945 -6.072 143 I(VL-HT)98 14 33 55.4058 109 54.6835 1614139.561 936432.576 -3.637 144 III(PLK-CH)4 14 36 44.0121 107 40 58.4810 1617145.262 789039.045 -6.558 145 I(DN-BMT)52 14 37 13.0164 107 51 15.9235 1618262.906 807518.398 -5.897 146 II(TT-KT)18 14 38 47.9431 108 23 47.6300 1621989.594 865936.933 -3.570 147 III(TN-AL)3 14 39 38.1245 108 44 28.5263 1624120.562 903092.000 -3.932 148 III(MB-CN)7 14 41 55.8662 108 14 31.8646 1627527.035 849205.580 -4.078 149 I(DN-BMT)47A 14 42 18.4785 107 41 11.4295 1627435.428 789304.654 -6.831 150 I(VL-HT)95 14 42 40.8804 109 19.4035 1630299.918 935085.890 -4.217 151 III(MT-AL)2 14 43 16.4735 108 56 37.5104 1631214.167 924821.329 -4.226 152 II(TT-KT)11 14 45 42.8237 108 32 43.8115 1635003.642 881800.082 -4.481 153 II(TT-KT)4 14 48 37.6828 108 45 21.1565 1640754.691 904391.797 -4.785 154 III(TMR-DT)2 14 52 26.8789 107 57 50.0825 1646523.672 818952.636 -5.394 155 III(TN-LT)3 14 53 47.7955 108 31 50.0095 1649902.822 879953.250 -5.465 156 III(MR-MB)9 14 55 40.2254 108 39.0192 1652737.618 838283.330 -5.109 157 I(VL-HT)87 14 56 39.1386 108 53 28.4154 1655825.679 918725.211 -5.302 158 I(DN-BMT)40 14 58 32.6194 107 43 37.7029 1657446.207 793317.553 -7.093 159 III(HT-TN)2 15 50.6070 108 33 50.2592 1666665.360 883283.229 -6.220 160 III(HD-MR)11 15 14.9478 108 57.7714 1672297.934 840359.174 -6.409 161 I(VL-HT)83 15 46.9056 108 48 12.4520 1674370.320 908945.978 -6.077 162 I(DN-BMT)35 15 18.5936 107 44 55.5310 1673650.986 795442.675 -7.317 163 III(LV-ML)9 15 10 56.5214 108 39 41.6282 1681793.295 893542.536 -6.578 164 III(HD-BHB)8 15 15 33.4764 108 30 26.6040 1690044.173 876816.532 -6.806 165 III(LV-ML)3 15 16 13.2624 108 43 6.7876 1691646.261 899508.649 -6.742 166 III(KD-DC)8 15 17 8.6550 107 53 35.1722 1691999.696 810727.820 -7.618 167 III(HD-BHB)3-1 15 17 17.4152 108 19 3.3406 1692921.921 856356.676 -7.299 168 I(VL-HT)78 15 17 27.3060 108 45 42.3169 1694005.464 904116.450 -6.733 169 I(DN-BMT)28 15 18 28.6161 107 43 48.2678 1694232.000 793175.252 -8.215 170 III(BHA-HD)9 15 20 30.4450 108 10 36.4228 1698634.062 841130.832 -7.928 Diem tay nam: x = 1294150.946; y = 717836.918 Diem dong bac: x = 1698634.062; y = 977841.969 DelN(tb0) = 81514 m Phụ lục 3: TINH BAN PHUONG SAI THUC NGHIEM DIEM SONG TRUNG TAY NGUYEN So diem: 170 dN(tb)= 8151 m TT x(km) y(km) del(m) 1359.9010 782.2480 862 1294.1510 804.8520 618 1299.2220 790.4870 531 1301.9450 943.8850 864 1302.9090 866.5690 723 1302.9080 805.3120 644 1306.3240 910.2420 574 1306.7320 744.6790 606 1309.4970 889.2600 726 10 1314.4440 947.4550 925 11 1314.3030 818.6760 699 12 1319.3690 799.6460 676 13 1320.3710 866.2680 619 14 1321.6670 909.9440 463 15 1320.1570 764.5670 569 16 1321.4030 844.4480 709 17 1323.0070 880.2590 741 18 1323.3560 717.8370 938 19 1327.7800 934.6610 565 20 1329.1950 956.6410 1.003 21 1326.6870 732.5140 808 22 1327.3370 773.9360 589 23 1331.7510 726.0970 1.007 24 1333.9320 916.0310 450 25 1340.5220 866.0640 787 26 1340.0520 789.1280 668 27 1340.0060 769.9290 555 28 1342.3810 949.0180 876 29 1341.6700 805.7950 613 30 1348.4240 862.4710 801 31 1348.1390 815.3090 650 32 1347.6580 739.1240 1.040 33 1352.1880 840.0720 768 34 1356.6630 752.6050 897 35 1362.0010 956.0330 1.063 36 1361.0600 827.0240 782 37 1361.2920 841.4640 813 38 1364.0840 779.7500 900 39 1372.5780 929.6140 697 40 1371.6140 813.1510 1.136 41 1374.2720 953.8000 946 42 1382.1060 888.3260 872 detb(m) 047 -.197 -.284 049 -.092 -.171 -.241 -.209 -.089 110 -.116 -.139 -.196 -.352 -.246 -.106 -.074 123 -.250 188 -.007 -.226 192 -.365 -.028 -.147 -.260 061 -.202 -.014 -.165 225 -.047 082 248 -.033 -.002 085 -.118 321 131 057 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 1382.1510 1382.7380 1385.8130 1387.7790 1389.7220 1392.4800 1394.4160 1396.7530 1394.5600 1404.1340 1403.6050 1405.5250 1406.1660 1410.5610 1412.1420 1415.8460 1418.1090 1417.3700 1418.4970 1420.4290 1421.9580 1426.9160 1428.4450 1427.3930 1429.6800 1435.7510 1437.2020 1442.7860 1445.2150 1446.3640 1448.8690 1456.9770 1459.9210 1460.3200 1460.5270 1462.3350 1465.8740 1465.4390 1466.4940 1468.7860 1469.0870 1471.8760 1473.2580 1478.6430 1480.8320 1481.9170 1487.0900 785.5170 797.8440 898.1230 866.2360 844.5560 808.9670 919.3890 954.2790 781.1630 904.3450 824.6600 957.8480 867.6130 907.8030 781.4600 816.8460 884.4550 792.1760 873.2890 831.1500 809.9510 977.8420 934.0320 805.5810 854.6650 904.9570 781.1580 956.1520 908.1490 849.2000 931.2570 777.5380 805.5250 819.5960 779.0150 829.1830 962.3800 857.2820 839.1920 929.5470 892.8020 783.1380 841.6360 802.7990 879.5730 929.8090 956.4480 1.076 1.084 737 812 1.076 1.130 741 859 1.068 862 1.160 852 1.099 909 1.097 1.057 1.054 1.001 1.094 1.148 954 931 762 1.057 1.150 884 1.011 749 634 1.134 779 1.111 1.035 981 1.075 1.049 888 887 1.097 792 648 930 1.116 909 684 419 896 261 269 -.078 -.003 261 315 -.074 044 253 047 345 037 284 094 282 242 239 186 279 333 139 116 -.053 242 335 069 196 -.066 -.181 319 -.036 296 220 166 260 234 073 072 282 -.023 -.167 115 301 094 -.131 -.396 081 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 1491.5160 1496.0650 1494.7850 1495.6260 1499.7620 1501.1580 1502.9880 1503.8120 1505.6700 1508.0500 1507.2090 1515.6560 1516.6570 1517.6960 1523.3560 1523.4280 1525.2870 1527.0330 1527.2740 1529.2270 1538.3990 1537.9110 1540.2580 1545.1450 1544.4640 1545.8660 1546.6580 1547.0910 1547.9050 1549.2320 1549.2840 1553.3800 1558.3640 1565.3050 1566.7050 1569.7820 1570.5890 1576.4790 1577.8450 1577.1390 1576.5660 1588.0740 1590.6860 1589.9760 1591.2350 1594.1780 1594.9030 835.5830 928.2430 797.7950 863.8480 873.7680 806.8510 835.5490 849.1440 783.4990 933.0150 815.9810 902.4560 952.0510 880.7550 854.3880 827.3600 949.0500 883.0800 789.7940 796.3120 931.4580 886.7720 824.9440 940.9060 805.7540 895.2170 913.6450 876.3700 815.4000 835.9940 776.1780 865.5450 896.5480 818.6790 765.5290 887.8880 937.1210 879.2850 949.5090 889.0870 835.5910 875.2210 899.0550 773.4230 833.9380 934.1280 820.4350 1.015 630 881 774 618 877 974 895 837 635 844 685 1.008 869 792 909 951 783 668 703 757 866 951 1.050 682 855 702 888 867 982 608 886 897 811 172 892 906 903 1.186 873 685 799 780 169 768 928 743 200 -.185 066 -.041 -.197 062 159 080 022 -.180 029 -.130 193 054 -.023 094 136 -.032 -.147 -.112 -.058 051 136 235 -.133 040 -.113 073 052 167 -.207 071 082 -.004 -.643 077 091 088 371 058 -.130 -.016 -.035 -.646 -.047 113 -.072 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 1596.4470 1599.5370 1604.3800 1604.5870 1607.7920 1606.8430 1614.1400 1617.1450 1618.2630 1621.9900 1624.1210 1627.5270 1627.4350 1630.3000 1631.2140 1635.0040 1640.7550 1646.5240 1649.9030 1652.7380 1655.8260 1657.4460 1666.6650 1672.2980 1674.3700 1673.6510 1681.7930 1690.0440 1691.6460 1692.0000 1692.9220 1694.0050 1694.2320 1698.6340 800.7640 780.6720 922.7660 847.6820 871.6180 794.8290 936.4330 789.0390 807.5180 865.9370 903.0920 849.2060 789.3050 935.0860 924.8210 881.8000 904.3920 818.9530 879.9530 838.2830 918.7250 793.3180 883.2830 840.3590 908.9460 795.4430 893.5430 876.8170 899.5090 810.7280 856.3570 904.1160 793.1750 841.1310 482 -.333 243 -.572 726 -.089 676 -.139 611 -.204 343 -.472 952 137 390 -.425 555 -.260 680 -.135 665 -.150 720 -.095 450 -.365 1.049 234 940 125 650 -.165 663 -.152 712 -.103 732 -.083 900 085 970 155 568 -.247 834 019 908 093 911 096 762 -.053 870 055 686 -.129 948 133 593 -.222 663 -.152 974 159 851 036 692 -.123 dmax= 371;dmin= -.646;dl= 1.017;del= 102 Phụ lục 4: TINH HIEP PHUONG SAI THUC NGHIEM DIEM SONG TRUNG TAY NGUYEN So diem:170 dN(tb)= 000118 m TT x(km) y(km) del(m) 1359.901 782.248 4.686 1294.151 804.852 -19.714 1299.222 790.487 -28.414 1301.945 943.885 4.886 1302.909 866.569 -9.214 1302.908 805.312 -17.114 1306.324 910.242 -24.114 1306.732 744.679 -20.914 1309.497 889.260 -8.914 10 1314.444 947.455 10.986 11 1314.303 818.676 -11.614 12 1319.369 799.646 -13.914 13 1320.371 866.268 -19.614 14 1321.667 909.944 -35.214 15 1320.157 764.567 -24.614 16 1321.403 844.448 -10.614 17 1323.007 880.259 -7.414 18 1323.356 717.837 12.286 19 1327.780 934.661 -25.014 20 1329.195 956.641 18.786 21 1326.687 732.514 -.714 22 1327.337 773.936 -22.614 23 1331.751 726.097 19.186 24 1333.932 916.031 -36.514 25 1340.522 866.064 -2.814 26 1340.052 789.128 -14.714 27 1340.006 769.929 -26.014 28 1342.381 949.018 6.086 29 1341.670 805.795 -20.214 30 1348.424 862.471 -1.414 31 1348.139 815.309 -16.514 32 1347.658 739.124 22.486 33 1352.188 840.072 -4.714 34 1356.663 752.605 8.186 35 1362.001 956.033 24.786 36 1361.060 827.024 -3.314 37 1361.292 841.464 -.214 38 1364.084 779.750 8.486 39 1372.578 929.614 -11.814 40 1371.614 813.151 32.086 41 1374.272 953.800 13.086 42 1382.106 888.326 5.686 43 1382.151 785.517 26.086 44 1382.738 797.844 26.886 45 1385.813 898.123 -7.814 46 1387.779 866.236 -.314 47 1389.722 844.556 26.086 48 1392.480 808.967 31.486 49 1394.416 919.389 -7.414 50 1396.753 954.279 4.386 51 1394.560 781.163 25.286 52 1404.134 904.345 4.686 53 1403.605 824.660 34.486 54 1405.525 957.848 3.686 55 1406.166 867.613 28.386 56 1410.561 907.803 9.386 57 1412.142 781.460 28.186 58 1415.846 816.846 24.186 59 1418.109 884.455 23.886 60 1417.370 792.176 18.586 61 1418.497 873.289 27.886 62 1420.429 831.150 33.286 63 1421.958 809.951 13.886 64 1426.916 977.842 11.586 65 1428.445 934.032 -5.314 66 1427.393 805.581 24.186 67 1429.680 854.665 33.486 68 1435.751 904.957 6.886 69 1437.202 781.158 19.586 70 1442.786 956.152 -6.614 71 1445.215 908.149 -18.114 72 1446.364 849.200 31.886 73 1448.869 931.257 -3.614 74 1456.977 777.538 29.586 75 1459.921 805.525 21.986 76 1460.320 819.596 16.586 77 1460.527 779.015 25.986 78 1462.335 829.183 23.386 79 1465.874 962.380 7.286 80 1465.439 857.282 7.186 81 1466.494 839.192 28.186 82 1468.786 929.547 -2.314 83 1469.087 892.802 -16.714 84 1471.876 783.138 11.486 85 1473.258 841.636 30.086 86 1478.643 802.799 9.386 87 1480.832 879.573 -13.114 88 1481.917 929.809 -39.614 89 1487.090 956.448 8.086 90 1491.516 835.583 19.986 91 1496.065 928.243 -18.514 92 1494.785 797.795 6.586 93 1495.626 863.848 -4.114 94 1499.762 873.768 -19.714 95 1501.158 806.851 6.186 96 1502.988 835.549 15.886 97 1503.812 849.144 7.986 98 1505.670 783.499 2.186 99 1508.050 933.015 -18.014 100 1507.209 815.981 2.886 101 1515.656 902.456 -13.014 102 1516.657 952.051 19.286 103 1517.696 880.755 5.386 104 1523.356 854.388 -2.314 105 1523.428 827.360 9.386 106 1525.287 949.050 13.586 107 1527.033 883.080 -3.214 108 1527.274 789.794 -14.714 109 1529.227 796.312 -11.214 110 1538.399 931.458 -5.814 111 1537.911 886.772 5.086 112 1540.258 824.944 13.586 113 1545.145 940.906 23.486 114 1544.464 805.754 -13.314 115 1545.866 895.217 3.986 116 1546.658 913.645 -11.314 117 1547.091 876.370 7.286 118 1547.905 815.400 5.186 119 1549.232 835.994 16.686 120 1549.284 776.178 -20.714 121 1553.380 865.545 7.086 122 1558.364 896.548 8.186 123 1565.305 818.679 -.414 124 1566.705 765.529 -64.314 125 1569.782 887.888 7.686 126 1570.589 937.121 9.086 127 1576.479 879.285 8.786 128 1577.845 949.509 37.086 129 1577.139 889.087 5.786 130 1576.566 835.591 -13.014 131 1588.074 875.221 -1.614 132 1590.686 899.055 -3.514 133 1589.976 773.423 -64.614 134 1591.235 833.938 -4.714 135 1594.178 934.128 11.286 136 1594.903 820.435 -7.214 137 1596.447 800.764 -33.314 138 1599.537 780.672 -57.214 139 1604.380 922.766 -8.914 140 1604.587 847.682 -13.914 141 1607.792 871.618 -20.414 142 1606.843 794.829 -47.214 143 1614.140 936.433 13.686 144 1617.145 789.039 -42.514 145 1618.263 807.518 -26.014 146 1621.990 865.937 -13.514 147 1624.121 903.092 -15.014 148 1627.527 849.206 -9.514 149 1627.435 789.305 -36.514 150 1630.300 935.086 23.386 151 1631.214 924.821 12.486 152 1635.004 881.800 -16.514 153 1640.755 904.392 -15.214 154 1646.524 818.953 -10.314 155 1649.903 879.953 -8.314 156 1652.738 838.283 8.486 157 1655.826 918.725 15.486 158 1657.446 793.318 -24.714 159 1666.665 883.283 1.886 160 1672.298 840.359 9.286 161 1674.370 908.946 9.586 162 1673.651 795.443 -5.314 163 1681.793 893.543 5.486 164 1690.044 876.817 -12.914 165 1691.646 899.509 13.286 166 1692.000 810.728 -22.214 167 1692.922 856.357 -15.214 168 1694.005 904.116 15.886 169 1694.232 793.175 3.586 170 1698.634 841.131 -12.314 Phụ lục 5: THONG KE KET QUA NOI SUY KRIGING SU DUNG HAM BAN PHUONG SAI CAU Theo so lieu:170 DIEM SONG TRUNG TAY NGUYEN Ty so kr = 1.000 r0 = 147.07 km Deltb = 8151 m TT DIEM H N Delz DEL h(tinh) h(tc) delh nk I(BMT-APD)35 294.751 -2.777 -.142 673 296.855 296.987 -.132 36 III(LS-BN)21 19.831 4.371 009 824 14.636 14.769 -.133 50 III(DX-DL)3 516.170 1.910 -.042 773 513.487 513.423 064 86 I(BMT-APD)1-2 431.389 -.998 339 1.154 431.233 431.204 029 104 II(MT-TH)21 127.255 526 -.162 653 126.076 126.039 037 III(DM-CR)7 162.938 -1.638 -.147 668 163.908 163.925 -.017 114 III(DC-NB)4 304.938 -5.587 -.087 728 309.797 309.713 084 III(XH-SL)16 496.065 -2.386 082 897 497.554 497.527 027 100 III(AL-DT)1 33.897 -3.723 -.014 801 36.819 36.866 -.047 74 10 III(MR-HT)3 640.588 -6.129 030 845 645.872 645.737 135 62 109 84 n = 10; SSTP(h) = 084 m THONG KE KET QUA NOI SUY KRIGING SU DUNG HAM BAN PHUONG SAI CAU Theo so lieu:DIEM SONG TRUNG TAY NGUYEN Ty so kr = 1.250 r0 = 183.84 km Deltb = 8151 TT DIEM H N Delz DEL h(tinh) I(BMT-APD)35 294.751 -2.777 -.180 III(LS-BN)21 19.831 4.371 029 844 III(DX-DL)3 516.170 1.910 -.004 I(BMT-APD)1-2 431.389 -.998 II(MT-TH)21 127.255 635 296.893 h(tc) delh nk 296.987 -.094 53 14.616 14.769 -.153 77 811 513.449 513.423 026 105 317 1.132 431.255 431.204 051 127 526 -.212 603 26.126 126.039 087 134 III(DM-CR)7 162.938 -1.638 -.100 715 163.861 163.925 -.064 142 III(DC-NB)4 304.938 -5.587 -.051 764 309.761 309.713 048 127 III(XH-SL)16 496.065 -2.386 070 885 497.566 497.527 039 119 III(AL-DT)1 33.897 -3.723 021 836 36.784 36.866 -.082 96 10 III(MR-HT)3 640.588 -6.129 053 868 645.849 645.737 112 80 n = 10; SSTP(h) = 084 m THONG KE KET QUA NOI SUY COLLOCATION SU DUNG HAM HIEP PHUONG SAI MARKOV BAC Theo so lieu:DIEM SONG TRUNG TAY NGUYEN SO LIEU NGUYEN DANG-CHUAN Gia tri N(MIN): -64.61cm; Gia tri N(MAX): 37.09cm Khoang rong N(MAX)-N(MIN): 101.70cm TT Ten H(m) N(m) Restb Noisu h(tinh) h(tc) sai Dmin Dmax I(BMT-APD)35 294.751 -2.777 -.150 665 296.863 296.987 -.124 16.6 434.1 III(LS-BN)21 19.831 4.371 007 822 14.638 III(DX-DL)3 516.170 1.910 -.029 787 513.473 14.769 513.423 I(BMT-APD)1-2 431.389 -.998 360 1.175 431.212 431.204 II(MT-TH)21 526 127.255 -.228 587 126.142 -.131 11.3 398.8 050 10.9 351.0 008 4.9 302.2 126.039 103 8.6 267.9 III(DM-CR)7 162.938 -1.638 -.178 637 163.939 163.925 014 7.0 229.9 III(DC-NB)4 304.938 -5.587 -.097 718 309.807 309.713 094 9.2 269.3 III(XH-SL)16 496.065 -2.386 025 299.1 III(AL-DT)1 33.897 -3.723 -.042 773 10 III(MR-HT)3 640.588 -6.129 084 899 497.552 497.527 36.847 36.866 -.019 2.8 13.5 357.6 082 898 645.819 645.737 082 14.3 Gia tri noi suy MIN: -.2279; noi suy MAX: 3603 Sai lech trung phuong: 0790 m Sai(min) = -.1308 m; Sai(max) = 1028 m 377.1 ... ? ?Ứng dụng địa thống kê phân tích mơ hình hóa liệu dị thường độ cao vùng Tây Nguyên – Việt Nam? ?? Mục đích nghiên cứu luận văn: Nghiên cứu phát triển ứng dụng phƣơng pháp địa thống kê xử lý số liệu. .. cứu liệu dị thƣờng độ cao vùng Tây Nguyên – Việt Nam Nội dung nghiên cứu: - Tìm hiểu lý thuyết địa thống kê; - Dữ liệu dị thƣờng độ cao ứng dụng đo cao GPS; - Các đặc trƣng địa thống kê liệu dị. .. hệ thống độ cao dị thƣờng độ cao - Chƣơng 3: Dữ liệu dị thƣờng độ cao GPS – Thủy chuẩn vùng Tây Nguyên - Chƣơng 4: Sử dụng phƣơng pháp Kriging để mơ hình hóa liệu dị thƣờng độ cao vùng Tây Nguyên