Nghiên cứu chính xác hoá dị thường độ cao egm 2008 dựa trên số liệu gps thuỷ chuẩn trên phạm vi cục bộ ở việt nam

bộ giáo dục đào tạo trờng đại học mỏ - địa chất nguyễn Đô nghiên cứu xác hóa dị thờng độ cao EGM2008 dựa số liệu gps - thủy chuẩn phạm vi cục việt nam luận án tiến sĩ kĩ thuật Hà nội - 2012 giáo dục đào tạo trờng đại học mỏ - địa chất nguyễn Đô nghiên cứu xác hóa dị thờng độ cao EGM2008 dựa số liệu gps - thủy chuẩn phạm vi cục bé ë viÖt nam Chuyên ngành : Trắc địa cao cấp Mã số : 62.52.85.10 luËn ¸n tiÕn sÜ kÜ thuật ngời Hớng dẫn khoa học PGS TS Đặng Nam Chinh TS Lê Minh Tá Hà nội - 2012 i LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu riêng thân 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 trình bày luận án xác, trung thực chưa cơng bố cơng trình khác Tác giả Nguyễn Duy Đô ii MỤC LỤC Lời cam đoan ……………………………………………………………………………… Mục lục …………………………………………………………………………………… Danh mục hình ………………………………………………………………………… Danh mục bảng ……………………………………………………………………… MỞ ĐẦU ………………………………………………………………………………… Chương TỔNG QUAN ………………………………………………………………… 1.1 Tình hình nghiên cứu giới 1.2 Tình hình nghiên cứu Việt nam Chương MƠ HÌNH GEOID/QUASIGEOID VÀ VAI TRỊ CỦA NĨ TRONG CÔNG TÁC ĐO CAO BẰNG GPS ……………………………………………………………… 2.1 Các sở gốc độ cao loại độ cao …………………………………………… 2.1.1 Các sở gốc độ cao ………………………………………………………… 2.1.2 Các loại độ cao ……………………………………………………………… 2.2 Mơ hình Geoid/Quasigeoid phương pháp xây dựng …………………… 2.2.1 Mơ hình Geoid phân loại mơ hình Geoid ………………………………… 2.2.2 Một số phương pháp xây dựng mơ hình Geoid ……………………………… 2.2.3 Một số mơ hình Geoid tồn cầu cục …………………………………… 2.3 Đo cao GPS vai trò mơ hình Geoid đo cao GPS ……………………… 2.3.1 Đo cao GPS …………………………………………………………………… 2.3.2 Vai trị mơ hình Geoid/Quasigeoid đo cao GPS …………………… Chương ĐÁNH GIÁ ĐỘ CHÍNH XÁC MỘT SỐ MƠ HÌNH GEOID TỒN CẦU TRÊN LÃNH THỔ VIỆT NAM ………………………………………………………… 3.1 Phương pháp khảo sát đánh giá độ xác mơ hình Geoid tồn cầu dạng lưới… 3.1.1 Mơ hình Geoid tồn cầu lãnh thổ Việt Nam …………………………… 3.1.2 Phương pháp đánh giá độ xác mơ hình Geoid tồn cầu dạng lưới …… 3.2 Khảo sát đánh giá độ xác số mơ hình Geoid lãnh thổ Việt Nam …… 3.2.1 Giới thiệu số mơ hình Geoid tồn cầu để khảo sát ……………………… 3.2.2 Khảo sát đánh giá số mơ hình Geoid tồn cầu giới thiệu …………… Chương 4, NGHIÊN CỨU PHƯƠNG PHÁP CHÍNH XÁC HĨA DI THƯỜNG ĐỘ CAO EGM2008 KHU VỰC TÂY NGUYÊN VÀ DUYÊN HẢI NAM TRUNG BỘ ………… 4.1 Phần lý thuyết ……………………………………………………………………… 4.1.1 Quá trình ngẫu nhiên bề mặt trái đất …………………………………… 4.1.2 Sự thay đổi hiệu độ cao trắc địa tính chuyển tọa độ ……………………… 4.2 Phần thực nghiệm …………………………………………………………………… 4.2.1 Giới thiệu khu vực thực nghiệm ……………………………………………… 4.2.2 Quy trình tính tốn chương trình máy tính ………………………………… 4.2.3 Tính tốn thực nghiệm ………………………………………………………… 4.2.4 Tính tốn thử nghiệm để so sánh …………………………………………… KẾT LUẬN VÀ KIẾN NGHỊ …………………………………………………………… Các cơng trình nghiên cứu khoa học ……………………………………………………… Tài liệu tham khảo ……………………………………………………………………… Phụ lục …………………………………………………………………………………… i ii iii iv 6 12 16 16 16 20 26 26 27 30 37 37 41 42 42 42 42 46 46 53 63 63 63 82 86 86 90 97 112 123 125 126 131 iii DANH MỤC CÁC HÌNH Hình 1.1 Bản đồ Geoid với khoảng cao 10m (1964) ……………………… Hình 1.2 Bản đồ Geoid GRIM-3 ………………………………………………… Hình 1.3 Bản đồ Geoid EGM96 ………………………………………………… Hình 2.1 Quasigeoid Geoid …………………………………………………… 17 Hình 2.2 Hai sở gốc cao A B …………………………………………… 18 Hình 2.3 Độ cao …………………………………………………………… 22 Hình 2.4 Độ cao chuẩn …………………………………………………………… 23 Hình 2.5 Mối quan hệ độ cao độ cao chuẩn ……………………… 25 Hình 2.6 Mặt Geoid tồn cầu mực nước trung bình khu vực ………………… 26 Hình 2.7 Mơ hình Quasigeoid Cộng Hịa Liên Bang Đức GCG05 ……………… 32 Hình 2.8 Mơ hình Quasigeoid Ba Lan ………………………………… …… 33 Hình 2.9 Mơ hình WWGEOID04 34 Hình 2.10 Mơ hình EGMGEOID05 34 Hình 2.11 Bản đồ Geoid KOGD2003 Hàn Quốc ……………………….…… 35 Hình 2.12 Mơ hình Geoid AUSGEOID98 Australia …………………… … 36 Hình 2.13 Mối quan hệ độ cao ellipsoid độ cao chuẩn …………………… 37 Hình 2.14 Hiệu độ cao trắc địa hiệu độ cao thủy chuẩn ……………………… 38 Hình 3.1 Mơ hình Geoid DMA10 ………………………………………………… 48 Hình 3.2 Mơ hình Geoid USU91A ……………………………………………… 49 Hình 3.3 Mơ hình Geoid EGM96 ………………………………………………… 50 Hình 3.4 Mơ hình Geoid EGM2008 ……………………………………………… 52 Hình 4.1 Dạng hàm Hirvonen …………………………………………………… 67 Hình 4.2 Phân chia vùng để xác định đặc trưng hiệp phương sai dị thường trọng lực … 78 Hình 4.3 Sơ đồ khối chương trình tính tốn khảo sát …………………………… 84 Hình 4.4 Một phần vùng đất Tây Nguyên ……….………………………………… 87 Hình 4.5 Phạm vi phân bố điểm song trùng …………………………………… 89 Hình 4.6 Tạo mơ hình Geoid cục (grid) từ mơ hình tồn cầu ….……………… … 92 Hình 4.7 Tạo mơ hình EGM08 cục vùng Tây Ngun ………………………… 93 Hình 4.8 Quy trình tính tốn xác hóa dị thường độ cao …………………… 94 Hình 4.9 Đồ thị hàm hiệp phương sai Markov bậc ……………………………… 99 Hình 4.10 Đồ thị hàm hiệp phương sai Markov bậc …………………………… 104 Hình 4.11 Sơ đồ 2D số hiệu chỉnh δζ i …………………………………… 106 Hình 4.12 Mơ hình xác hóa EGM08C …………………………………… Hình 4.13 Hiệu việc xác hóa dị thường độ cao Geoid …………………… Hình 4.14 Đồ thị hàm hiệp phương sai thực nghiệm …………………………… 107 112 114 iv Hình 4.15 Sơ đồ 2D Quasigeoid GPS-TC xác lập từ 163 điểm song trùng … Hình 4.16 Phân bố 110 điểm GPS-TC ………………………………………… Hình 4.17 Hiệp phương sai dị thường độ cao GPS-TC xác định từ 110 điểm Hình 4.18 Sơ đồ 2D Quasigeoid GPS-TC xác lập từ 110 điểm song trùng … Hình 4.19 Hiệp phương sai số chênh dị thường độ cao ……………………… Hình 4.20 Sơ đồ 2D Quasigeoid xác hóa …………………… Hình 4.21 Phân bố 60 điểm GPS-TC ………………………………………… Hình 4.22 Hiệp phương sai số chênh dị thường độ cao GPS-TC ……………… Hình 4.23 Sơ đồ 2D Quasigeoid GPS-TC xác lập từ 60 điểm song trùng … Hình 4.24 Đồ thị hiệp phương sai số chênh dị thường độ cao ………………… Hình 4.25 Sơ đồ 2D Quasigeoid ………………………………………… 114 116 117 117 118 118 119 120 121 121 122 v DANH MỤC CÁC BẢNG Bảng 2.1.Tổng hợp số đặc điểm độ cao độ cao chuẩn …… 25 Bảng 2.2 Sự khác ζ N trường hợp khác biệt độ cao ……… 26 Bảng 2.3 Một số mơ hình geoid tồn cầu ……………………………………… 31 Bảng 3.1 Một số mơ hình Geoid tồn cầu có Việt Nam 46 Bảng 3.2 Số liệu thống kê giá trị trung bình dị thường trọng lực chân không theo ô 30’x30’ mơ hình EGM96 47 Bảng 3.3 Thống kê liệu ∆g (5’x5’) mGal 51 Bảng 4.1 Hiệp phương sai thực nghiệm tính theo 180 điểm 99 Bảng 4.2 Các tham số hàm hiệp phương sai Markov bậc 99 Bảng 4.3 Tọa độ loại độ cao điểm song trùng 100 Bảng 4.4 Hiệp phương sai thực nghiệm tính theo 163 điểm 103 Bảng 4.5 Các tham số hàm hiệp phương sai Markov bậc 104 Bảng 4.6 Giá trị dị thường độ cao nội suy 163 điểm song trùng 108 Bảng 4.7 Giá trị độ cao Geoid nội suy từ mơ hình EGM08C 111 Bảng 4.8 Tổng hợp kết tính nội suy theo phương pháp 112 Bảng 4.9 Tổng hợp kết tính tốn trường hợp 163 điểm 115 Bảng 4.10 Tổng hợp kết tính tốn trường hợp 110 điểm 119 Bảng 4.11 Tổng hợp kết tính tốn trường hợp 60 điểm 122 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 giải cách hiệu quả, độ cao điểm lại chưa thể giải trọn vẹn 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 cục 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 Trên quy mơ tồn cầu, nhà khoa học nghiên cứu xây dựng số mơ hình trọng trường, tiêu biểu số mơ hình EGM2008, có bậc hạng tới 2190 Tương ứng xác định Geoid theo mơ hình với độ phân giải 2,5’x2,5’ độ xác khoảng dm So với mơ hình trọng trường tồn cầu trước đó, độ xác mơ hình EGM 2008 tăng lên đáng kể Tuy vậy, mơ hình EGM2008 chưa đáp ứng yêu cầu nghiên cứu khoa học số ứng dụng thực tiễn phạm vi cục bộ, đặc biệt ứng dụng công nghệ đo cao GPS Trên quy mơ cục bộ, mơ hình Geoid quốc gia, châu lục liên tục hoàn thiện nâng cấp thỏa mãn nhu cầu ngày cao người sử dụng Độ xác Geoid cục số quốc gia đạt mức cm Tiêu biểu Canada (2.5cm), Thụy sỹ (2.5cm), Balan (2cm) vv… Ở Việt Nam, từ đầu năm 90 thập kỷ trước, GS.TSKH Phạm Hồng Lân xây dựng mơ hình grid trọng lực 5’x5’ sở sử dụng giá trị trọng lực chi tiết mặt đất cho vùng lãnh thổ sử dụng hệ số hàm điều hịa cầu mơ hình Rapp có bậc hạng (n=m=180) để cải cho vùng xa tính tốn độ cao geoid Độ xác mơ hình geoid đạt khoảng 2,0m ÷3,0m Tiếp đó, giai đoạn 19982000, Tổng Cục Địa giao cho GS.TSKH Đặng Hùng Võ làm trưởng ban số nhà khoa học trắc địa xây dựng mô hình Geoid VnGeo 96R VnGeo 96U với độ xác từ 0.5m ÷0,8m [28] Giai đoạn 2002 -2005, có nghiên cứu xây dựng mơ hình Geoid lãnh thổ với độ xác 0,22m Một vài nét sơ lược cho thấy, Việt nam, có hai lựa chọn để khai thác Geoid phục vụ cho việc giải nhiệm vụ khoa học thực tiễn Thứ nhất, sử dụng mơ hình EGM2008 Thứ hai sử dụng Geoid cục xây dựng nhà khoa học nước Cả hai phương án cho dị thường độ cao đạt cỡ 2dm Độ xác chưa thỏa mãn nhu cầu ứng dụng công nghệ đo cao GPS số ứng dụng khác Trong bối cảnh chúng tơi lựa chọn đề tài “Nghiên cứu xác hóa dị thường độ cao EGM2008 dựa số liệu GPS - Thủy chuẩn phạm vi cục Việt Nam” với mục tiêu xây dựng mơ hình Quasigeoid cục có độ xác cao mức 2dm, nhằm đáp ứng tốt yêu cầu nghiên cứu khoa học thực tiễn ngành Trắc địa nói riêng ngành khoa học trái đất nói chung nước ta Mục đích nghiên cứu luận án Nghiên cứu đánh giá độ xác mơ hình Geoid tồn cầu lãnh thổ Việt Nam xây dựng mơ hình Quasigeoid cục dạng lưới (đã xác hóa) cho vùng Tây Nguyên duyên hải Nam Trung Bộ, hoàn toàn đáp ứng yêu cầu đo cao GPS tương đương thủy chuẩn hạng IV vùng núi Đối tượng nghiên cứu - Mơ hình Geoid/Quasigeoid phương pháp xác định chúng hệ quy chiếu trái đất - Phương pháp Collocation bình phương nhỏ số phương pháp nội suy khác - Sử dụng mô hình Geoid/Quasigeoid cục phục vụ cho cơng tác đo cao GPS Phạm vi nghiên cứu Nghiên cứu phương pháp đánh giá độ xác mơ hình Geoid; lựa chọn mơ hình Geoid tiên nghiệm xử lý số liệu GPS-TC để xác hóa dị thường độ cao vùng Tây Nguyên Duyên hải Nam Trung Bộ Nội dung nghiên cứu - Tìm hiểu lý thuyết mơ hình Geoid, phương pháp nội suy dị thường độ cao độ cao Geoid - Nghiên cứu lý thuyết thực nghiệm số phương pháp đánh giá độ xác mơ hình Geoid tồn cầu - Nghiên cứu tính tốn thực nghiệm thay đổi hiệu độ cao trắc địa tính chuyển tọa độ sai số tọa độ khởi tính - Nghiên cứu phương pháp quy trình xác hóa dị thường độ cao cho khu vực nhỏ; lập chương trình cho máy tính, cho phép tự động xử lý xác hóa dị thường độ cao mơ hình Geoid tiên nghiệm dạng lưới dựa số liệu tọa độ, độ cao điểm song trùng 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ật 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 Geoid/Quasigeoid; phương pháp xây dựng đánh giá mơ hình Geoid/Quasigeoid - Phương pháp so sánh: So sánh ưu điểm 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 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: Viết chương trình tính tốn cho máy tính thực quy trình xác hóa dị thường độ cao Ý nghĩa khoa học thực tiễn luận án - Việc nghiên cứu phương pháp xác hóa dị thường độ cao mơ hình để đạt mục đích nghiên cứu luận án tiến sĩ kỹ thuật góp phần bổ sung vào lý thuyết xây dựng mơ hình Quasigeoid cục có độ xác cao - Xác định độ cao thủy chuẩn điểm lưới GPS đo cao GPS khu vực có mơ hình Quasigeoid cục xác hóa Các luận điểm bảo vệ luận điểm luận án a) Các luận điểm bảo vệ Đánh giá chất lượng mô hình Geoid tiên nghiệm theo phương pháp so sánh hiệu dị thường độ cao cặp điểm cho kết tin cậy Các điểm song trùng GPS-Thủy chuẩn (GPS-TC) phải xác định hệ WGS84 với độ xác cần thiết Độ xác phụ thuộc vào diện tích khu vực cần xác hóa dị thường độ cao Dị thường độ cao EGM2008 (mơ hình tiên nghiệm) xác hóa sở xác định số chênh dị thường độ cao chuẩn hóa xác định điểm GPS-TC để hiệu chỉnh vào điểm mắt lưới mơ hình tiên nghiệm Kết cho mơ hình Quasigeoid cục giữ phù hợp tổng thể với Geoid toàn cầu với mơ hình cục vùng tiếp biên khu vực lân cận sau xác hóa Dựa phương pháp xấp xỉ hàm thuật toán Collocation, xây dựng chương trình máy tính gồm modul có tên PRECOV.EXE, COVFUN.EXE COLLO2.EXE, cho phép tự động xử lý xác hóa dị thường độ cao Geoid/Quasigeoid tiên nghiệm dạng lưới dựa số liệu tọa độ, độ cao điểm song trùng; chương trình GEOINT.EXE cho người sử dụng muốn nội suy dị thường độ cao điểm khu vực 185 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 II(MT-TH)21 II(MT-TH)21 II(MT-TH)21 II(MT-TH)21 II(MT-TH)21 II(MT-TH)21 III(DM-CR)7 III(DM-CR)7 III(DM-CR)7 III(DM-CR)7 III(DM-CR)7 III(DM-CR)7 III(DM-CR)7 I(VL-HT)123 I(VL-HT)123 I(VL-HT)123 I(VL-HT)123 I(VL-HT)123 I(VL-HT)123 III(DC-NB)4 III(DC-NB)4 III(DC-NB)4 III(DC-NB)4 III(DC-NB)4 III(XH-SL)16 III(XH-SL)16 III(XH-SL)16 III(XH-SL)16 III(CH-PQ)3 III(CH-PQ)3 III(CH-PQ)3 III(AL-DT)1 III(AL-DT)1 III(MR-HT)3 Del(min)= III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(MR-HT)3 III(HD-BHB)3-1 III(HD-BHB)3-1 -.038; Del(max)= 141.8 120.0 192.1 165.6 221.6 245.5 82.3 92.2 78.9 138.4 134.2 175.9 200.9 167.0 86.0 181.3 107.3 180.1 200.6 112.6 88.0 167.1 163.4 187.6 97.3 58.9 101.8 125.5 125.1 84.8 105.8 81.7 97.3 25.4 297 347 322 345 347 73 286 325 20 353 355 270 309 299 342 326 331 62 53 24 22 291 35 342 346 84 44 35 306 319 42 48 57 10 14 47 42 28 40 56 50 15 58 30 50 22 29 54 20 44 24 41 30 33 42 37 35 54.92 59.93 48.78 54.96 38.33 41.55 27.02 3.68 43.03 40.98 16.46 9.72 57.54 25.41 31.67 22.28 33.17 9.56 8.80 47.03 44.12 38.10 18.01 18.65 47.39 12.72 55.46 53.12 41.76 50.95 23.32 47.84 39.33 22.48 177.683 368.810 487.078 -93.358 513.333 29.365 70.259 142.000 333.127 451.395 -129.041 477.650 -6.318 71.741 262.868 381.136 -199.300 407.391 -76.577 191.127 309.395 -271.041 335.650 -148.318 118.268 -462.168 144.523 -339.445 -580.436 26.255 -457.713 606.691 122.723 -483.968 177.716 368.830 487.120 -93.339 513.371 29.406 70.246 142.021 333.134 451.425 -129.034 477.676 -6.290 71.775 262.888 381.179 -199.281 407.430 -76.536 191.114 309.404 -271.055 335.655 -148.311 118.291 -462.169 144.541 -339.424 -580.459 26.251 -457.715 606.710 122.745 -483.965 033 020 042 019 038 041 -.013 021 007 030 007 026 028 034 020 043 019 039 041 -.013 009 -.014 005 007 023 -.001 018 021 -.023 -.004 -.002 019 022 003 066 Sai khác có giá trị nhỏ là: - 0.038m (ở khoảng cách 215,4km) Sai khác có giá trị lớn : + 0.066 m (ở khoảng cách 319,1km) Như cho thấy, sai lệch vị trí tọa độ tuyệt đối hệ WGS84 1m dẫn đến sai lệch hiệu độ cao điểm khu vực đến cỡ 6,6cm, giá trị không nhỏ so với yêu cầu độ xác đo cao GPS 2) Giả thiết tham số dịch gốc có giá trị 0,1m Trong trường hợp ta giả thiết dX=dY=dZ=0,1m, tính sai khác hiệu đọ cao cặp điểm sau: TT DAU I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 10 I(BMT-APD)12 11 I(BMT-APD)12 12 I(BMT-APD)12 CUOI I(BMT-APD)35 III(QK-LT)5 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 D(km)PHUONGVI(0 ’ ’’) 81.5 61.5 167.1 68.9 118.7 62.4 68.0 158.8 161.1 231.1 159.5 234.7 217 157 101 100 77 48 24 53 34 47 26 30 58 46 23 29 54 24 17 14.52 3.99 42.48 16.19 41.27 31.84 22.35 10.41 48.38 16.81 26.86 44.68 ∆H(m) ∆H’(m) -612.927 -7.749 -887.847 -391.508 -415.266 -476.289 -705.204 -780.423 -744.740 -674.481 -602.740 -411.613 -612.927 -7.750 -887.851 -391.510 -415.268 -476.289 -705.204 -780.424 -744.740 -674.482 -602.738 -411.612 δH(m) 000 -.001 -.004 -.002 -.002 000 000 -.001 000 -.001 002 001 186 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 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 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(KRKM-MDR)10 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(QK-LT)5 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 247.3 293.0 316.0 341.2 73.1 215.4 128.7 188.6 143.2 148.4 238.1 242.5 311.7 230.2 315.0 317.7 373.6 391.6 416.5 142.4 63.5 124.4 101.2 119.1 183.7 201.7 259.3 217.4 279.5 304.1 335.9 368.0 393.3 98.3 76.7 139.1 166.5 132.5 182.1 191.0 251.7 252.4 319.1 294.2 354.0 377.9 61.2 57.4 84.2 121.7 146.7 197.2 183.4 225.6 264.8 279.6 319.8 345.1 70.8 95.3 68.8 109.8 141.5 175.0 185.4 244.0 27 13 12 83 81 65 61 42 31 48 35 44 13 28 11 29 18 17 80 44 48 13 34 19 34 354 16 358 19 7 282 321 302 305 344 336 320 346 331 354 344 346 51 338 331 29 31 339 347 14 282 292 10 346 22 319 356 334 54 54 32 59 51 51 18 16 23 21 27 44 58 17 59 16 45 50 13 31 13 57 41 28 39 36 51 35 27 59 29 51 16 32 30 51 17 44 49 56 17 44 27 33 49 14 58 57 21 41 35 17 25 55 23 45 57 59.34 23.02 12.99 45.03 29.17 14.89 6.97 43.37 3.24 29.19 29.14 5.62 30.80 42.58 45.57 14.89 25.99 4.25 53.45 27.54 26.11 35.21 57.56 52.35 26.96 44.46 10.94 13.72 46.08 28.38 8.59 6.90 23.10 49.75 52.85 2.58 53.95 22.15 21.16 8.98 8.85 50.13 43.53 51.70 19.14 48.49 8.56 13.02 51.63 20.37 29.95 56.07 33.24 10.97 39.49 7.51 10.22 47.81 40.60 24.93 15.07 36.77 58.23 16.20 50.31 6.83 -293.345 -873.781 -267.090 -751.058 605.178 -274.920 221.419 197.661 136.638 -92.277 -167.496 -131.813 -61.554 10.187 201.314 319.582 -260.854 345.837 -138.131 -880.098 -383.759 -407.517 -468.540 -697.455 -772.674 -736.991 -666.732 -594.991 -403.864 -285.596 -866.032 -259.341 -743.309 496.339 472.581 411.558 182.643 107.424 143.107 213.366 285.107 476.234 594.502 14.066 620.757 136.789 -23.758 -84.781 -313.696 -388.915 -353.232 -282.973 -211.232 -20.105 98.163 -482.273 124.418 -359.550 -61.023 -289.938 -365.157 -329.474 -259.215 -187.474 3.653 121.921 -293.342 -873.780 -267.088 -751.055 605.177 -274.924 221.417 197.659 136.637 -92.277 -167.497 -131.813 -61.556 10.189 201.315 319.585 -260.854 345.839 -138.128 -880.101 -383.759 -407.518 -468.539 -697.454 -772.674 -736.990 -666.732 -594.988 -403.862 -285.592 -866.030 -259.337 -743.305 496.341 472.583 411.561 182.647 107.426 143.111 213.368 285.113 476.238 594.509 14.070 620.763 136.795 -23.758 -84.780 -313.694 -388.915 -353.231 -282.973 -211.229 -20.103 98.167 -482.271 124.422 -359.546 -61.021 -289.936 -365.156 -329.472 -259.214 -187.470 3.656 121.926 003 001 002 003 -.001 -.004 -.002 -.002 -.001 000 -.001 000 -.002 002 001 003 000 002 003 -.003 000 -.001 001 001 000 001 000 003 002 004 002 004 004 002 002 003 004 002 004 002 006 004 007 004 006 006 000 001 002 000 001 000 003 002 004 002 004 004 002 002 001 002 001 004 003 005 187 79 III(KRKM-MDR)10 III(AL-DT)1 80 III(KRKM-MDR)10 III(MR-HT)3 81 III(KRKM-MDR)10 III(HD-BHB)3-1 82 I(BMT-APD)1-2 III(BD-BMT)4 83 I(BMT-APD)1-2 II(MT-TH)21 84 I(BMT-APD)1-2 III(DM-CR)7 85 I(BMT-APD)1-2 I(VL-HT)123 86 I(BMT-APD)1-2 III(DC-NB)4 87 I(BMT-APD)1-2 III(XH-SL)16 88 I(BMT-APD)1-2 III(CH-PQ)3 89 I(BMT-APD)1-2 III(AL-DT)1 90 I(BMT-APD)1-2 III(MR-HT)3 91 I(BMT-APD)1-2 III(HD-BHB)3-1 92 III(BD-BMT)4 II(MT-TH)21 93 III(BD-BMT)4 III(DM-CR)7 94 III(BD-BMT)4 I(VL-HT)123 95 III(BD-BMT)4 III(DC-NB)4 96 III(BD-BMT)4 III(XH-SL)16 97 III(BD-BMT)4 III(CH-PQ)3 98 III(BD-BMT)4 III(AL-DT)1 99 III(BD-BMT)4 III(MR-HT)3 100 III(BD-BMT)4 III(HD-BHB)3-1 101 II(MT-TH)21 III(DM-CR)7 102 II(MT-TH)21 I(VL-HT)123 103 II(MT-TH)21 III(DC-NB)4 104 II(MT-TH)21 III(XH-SL)16 105 II(MT-TH)21 III(CH-PQ)3 106 II(MT-TH)21 III(AL-DT)1 107 II(MT-TH)21 III(MR-HT)3 108 II(MT-TH)21 III(HD-BHB)3-1 109 III(DM-CR)7 I(VL-HT)123 110 III(DM-CR)7 III(DC-NB)4 111 III(DM-CR)7 III(XH-SL)16 112 III(DM-CR)7 III(CH-PQ)3 113 III(DM-CR)7 III(AL-DT)1 114 III(DM-CR)7 III(MR-HT)3 115 III(DM-CR)7 III(HD-BHB)3-1 116 I(VL-HT)123 III(DC-NB)4 117 I(VL-HT)123 III(XH-SL)16 118 I(VL-HT)123 III(CH-PQ)3 119 I(VL-HT)123 III(AL-DT)1 120 I(VL-HT)123 III(MR-HT)3 121 I(VL-HT)123 III(HD-BHB)3-1 122 III(DC-NB)4 III(XH-SL)16 123 III(DC-NB)4 III(CH-PQ)3 124 III(DC-NB)4 III(AL-DT)1 125 III(DC-NB)4 III(MR-HT)3 126 III(DC-NB)4 III(HD-BHB)3-1 127 III(XH-SL)16 III(CH-PQ)3 128 III(XH-SL)16 III(AL-DT)1 129 III(XH-SL)16 III(MR-HT)3 130 III(XH-SL)16 III(HD-BHB)3-1 131 III(CH-PQ)3 III(AL-DT)1 132 III(CH-PQ)3 III(MR-HT)3 133 III(CH-PQ)3 III(HD-BHB)3-1 134 III(AL-DT)1 III(MR-HT)3 135 III(AL-DT)1 III(HD-BHB)3-1 136 III(MR-HT)3 III(HD-BHB)3-1 Del(min)= -.004; Del(max)= 234.0 285.3 309.9 28.1 96.9 101.7 168.7 126.0 178.6 208.4 235.6 267.5 292.9 105.4 95.0 170.8 100.6 166.8 185.5 225.2 249.6 274.9 54.2 75.6 141.8 120.0 192.1 165.6 221.6 245.5 82.3 92.2 78.9 138.4 134.2 175.9 200.9 167.0 86.0 181.3 107.3 180.1 200.6 112.6 88.0 167.1 163.4 187.6 97.3 58.9 101.8 125.5 125.1 84.8 105.8 81.7 97.3 25.4 007 351 352 317 57 25 46 339 18 350 22 5 72 41 56 345 26 355 29 10 316 33 297 347 322 345 347 73 286 325 20 353 355 270 309 299 342 326 331 62 53 24 22 291 35 342 346 84 44 35 306 319 12 15 23 23 49 50 36 21 34 33 20 21 36 48 19 39 49 19 43 31 17 42 48 57 10 14 47 42 28 40 56 50 15 58 30 50 22 29 54 20 44 24 41 30 33 42 37 35 37.73 44.28 19.55 52.12 11.06 37.28 33.28 6.00 24.38 1.80 42.11 23.81 41.72 51.26 6.13 10.76 59.99 58.22 27.22 2.60 17.32 6.12 42.27 3.71 54.92 59.93 48.78 54.96 38.33 41.55 27.02 3.68 43.03 40.98 16.46 9.72 57.54 25.41 31.67 22.28 33.17 9.56 8.80 47.03 44.12 38.10 18.01 18.65 47.39 12.72 55.46 53.12 41.76 50.95 23.32 47.84 39.33 22.48 -458.515 148.176 -335.792 -228.915 -304.134 -268.451 -198.192 -126.451 64.676 182.944 -397.492 209.199 -274.769 -75.219 -39.536 30.723 102.464 293.591 411.859 -168.577 438.114 -45.854 35.683 105.942 177.683 368.810 487.078 -93.358 513.333 29.365 70.259 142.000 333.127 451.395 -129.041 477.650 -6.318 71.741 262.868 381.136 -199.300 407.391 -76.577 191.127 309.395 -271.041 335.650 -148.318 118.268 -462.168 144.523 -339.445 -580.436 26.255 -457.713 606.691 122.723 -483.968 -458.512 148.180 -335.787 -228.914 -304.135 -268.451 -198.193 -126.449 64.677 182.947 -397.491 209.202 -274.766 -75.221 -39.536 30.721 102.466 293.591 411.862 -168.577 438.116 -45.852 35.684 105.942 177.686 368.812 487.082 -93.356 513.337 29.369 70.258 142.002 333.128 451.398 -129.040 477.653 -6.315 71.744 262.870 381.140 -199.298 407.395 -76.573 191.126 309.396 -271.042 335.650 -148.317 118.270 -462.168 144.525 -339.443 -580.438 26.255 -457.713 606.693 122.725 -483.968 Sai khác có giá trị nhỏ là: - 0.004m (ở khoảng cách 215,4 km) Sai khác có giá trị lớn : + 0.007 m (ở khoảng cách 319,1 km) 003 004 005 001 -.001 000 -.001 002 001 003 001 003 003 -.002 000 -.002 002 000 003 000 002 002 001 000 003 002 004 002 004 004 -.001 002 001 003 001 003 003 003 002 004 002 004 004 -.001 001 -.001 000 001 002 000 002 002 -.002 000 000 002 002 000 188 Như cho thấy, sai lệch vị trí tọa độ tuyệt đối hệ WGS84 0.1m dẫn đến sai lệch hiệu độ cao điểm khu vực đến cỡ 0,7cm, giá trị nhỏ, chấp nhận Phụ lục 14 c 412 111 NOI SUY COLLOCATION THEO HAM HIEP PHUONG SAI (MARKOV3) integer n,i,j,k,l,m1,m,ites,itest character*12 ten1,BANG character*30 tct,tct2,dong3,dong1 character*1 c1,c2 character*15 ten(180),ten2(190) double precision r1,pi,ag,bg,fg,x(180),y(180),g(180),b(180,180), * cs(180),du(180),f(180),s1,s2,d12,u1,u2,x0,y0,dx,dy,xp,yp,de,df, * sn,z1,z2,up,rmax,c0,b0,sr,s0,demax,ua,vk(190),xt,yt,sumx,sumy, * sume,e1,e2,e3,emin,emax,a0,a1,a2,a3,h1,h2,ge,delslop,um,ge1,sk, * vmax,vmin,m0,x1,y1,gb1,ul1,ul2,gama,tyle,he1,he2,he3,truc,de1, * slomin,slomax,des,restb,e11,dmin,dmax,dh0,hm common b WRITE(*,*)' ======================================' write(*,*)' HANOI UNIVERSITY OF MINING AND GEOLOGY' WRITE(*,*)' Copyright by Geoid-Group ' WRITE(*,*)' NOI SUY COLLOCATION-Markov3' WRITE(*,*)' =================2010=================' open(5,file='collo2.par',status='old',err=866) READ(5,412)ITEC FORMAT(I2) open(7,file='Kernel.par',status='old',err=867) READ(7,412)ITEC2 IF(ITEC.NE.ITEC2) GOTO 139 IF(ITEC.EQ.0) THEN WRITE(*,*)' NOI SUY THEO SO LIEU CHUAN HOA TRUNG BINH' ENDIF IF(ITEC.EQ.1) THEN WRITE(*,*)' NOI SUY THEO SO LIEU CHUAN HOA HOI QUY' ENDIF write(*,*)' Hay viet ten FILE so lieu (goc)(*.xyn)' read(*,111)ten1 format(a) open(1,file=ten1,status='old') write(*,*)' Hay viet ten FILE so lieu (nsuy)(*.col)' read(*,111)ten1 open(2,file=ten1,status='old') write(*,*)' Hay viet ten FILE thong ke ket qua (.kq)' read(*,111)ten1 open(3,file=ten1,status='new') open(4,file='collo2.dat',status='new') IF(ITEC.EQ.1) THEN open(6,file='slop.dat',status='new') ENDIF open(8,file='zzz.yyy',status='new') open(9,file='collo2.vao',status='new') open(10,file='collo2.pos',status='new') open(11,file='colcheck.txt',status='new') write(*,*)' Co su dung pham vi gioi han S0 ? (0=khong, 1: co)' read(*,*)ites write(*,*)' Co TEN diem khong ? (0=khong(GRID), 1: co(DIEM))' 189 69 699 c c c 70 55 698 697 c c read(*,*)itest READ(1,111)TCT read(1,111)dong1 FORMAT(A30) r1=0.206264806247096D+6 pi=0.314159265358979D+1 ag=0.6378137D+7 fg=1./298.257223563 bg=ag-(ag*fg) m0=0.9996D+0 write(3,*)' THONG KE KET QUA NOI SUY COLLOCATION' write(3,*)' SU DUNG HAM HIEP PHUONG SAI MARKOV BAC 3' write(3,699)tct format(' Theo so lieu:',a30) write(3,111)dong1 read(1,*) m m: diem read(5,*)u1,u2,sr,s0 u1,u2 la cac tham so ham hiep phuong sai sr: khoang cach lien he, s0: khoang cach ket thuc sumx=0 sumy=0 emin=0 emax=0 i=1,m read(1,70)k,c1,ten(i),c2,x(i),y(i),g(i) format(i3,a1,a15,a1,2(f12.3),f8.1) write(*,55)i format(' Doc diem goc:',i3) if(g(i).lt.emin) emin=g(i) if(g(i).gt.emax) emax=g(i) sumx=sumx+x(i) sumy=sumy+y(i) continue de1=emax-emin write(3,698)emin,emax format(' Gia tri N(MIN):',f9.2,' cm; Gia tri N(MAX):',f9.2,' *cm') write(3,697)de1 format(' Khoang rong N(MAX)-N(MIN):',f9.2,' cm') if(itec.eq.1) then read(7,*)ido,iphu,giay truc=(ido*3600+iphu*60+giay)/r1 read(7,*)a0,a1,a2,a3 read(7,*)restb endif if(itec.eq.0) then read(7,*)restb endif xt=sumx/m yt=sumy/m doc xong toa cac diem goc g(i) k1=m*(m-1)/2 k=0 i=1,m b(i,i)=1 continue bat dau thiet lap ma tran hiep phuong sai 10 i=1,(m-1) 190 50 11 10 14 13 c c 144 c 740 623 624 11 j=(i+1),m k=k+1 dx=(x(j)-x(i))/1000 dy=(y(j)-y(i))/1000 sn=DSQRT(dx*dx+dy*dy) z1=-sn/u2 z2=1.-z1-0.5*z1*z1 b(i,j)=z2*EXP(z1) format(2(i3),f13.5) continue continue 13 i=1,m 14 j=1,m b(j,i)=b(i,j) continue continue ==== thiet lap xong ma tran hiep phuong read(2,111)tct2 CALL DMT2(m,b) write(*,*)' Dao xong MT C' format(/) ======= DOC TEP 2=========== read(2,*)n demax=0 sume=0 if(itest.eq.1) then if(itec.eq.1) then write(10,*)' KET QUA NOI SUY DE write(11,*)' KET QUA NOI SUY DE write(10,623) write(11,624) write(3,*)' TT Ten H(m) *) sai Delsl Dmin Dmax' endif if(itec.eq.0) then write(9,*)' KET QUA NOI SUY' write(10,*)' KET QUA NOI SUY DE write(11,*)' KET QUA NOI SUY DE write(9,740)restb format(f7.3) write(9,623) format('TT TEN X(m) * delN N(hc) h(tc)') write(10,623) write(11,624) format(' TT TEN DIEM H(m) * h(tinh) h(m) del') write(3,*)' TT Ten H(m) * sai Dmin Dmax' endif endif if(itest.eq.0) then if(itec.eq.0) then write(3,*)' TT X(m) Y(m) *c) Dmin Dmax' endif if(itec.eq.1) then write(3,*)' TT X(m) Y(m) *c) Dmin Dmax Delop' sai======= KIEM TRA' KIEM TRA' N(m) Delslop Noisu h(tinh KIEM TRA' KIEM TRA' Y(m) N(m) N(m) H(m) dN(m) N(m) N(hc) Restb Noisu h(tinh) N(08) HC(m) N(h/ N(08) HC(m) N(h/ 191 833 56 endif endif emin=0 emax=0 vmax=0 vmin=0 slomin=9999 slomax=0 i=1,n if(itest.eq.0) then read(2,*)k,xp,yp,ge endif if(itest.eq.1) then read(2,833)k,c1,ten2(i),xp,yp,h1,h2,ge format(i3,a1,a15,2(f12.3),2(f9.3),f8.3) endif write(*,56)i format(' Noi suy diem:',i5) x1=xp y1=yp-500000 CALL TINHBL(ag,bg,x1,y1,m0,gb1,ul1,gama,tyle) ul2=ul1+truc he1=DCOS(gb1)*DCOS(ul2) he2=DCOS(gb1)*DSIN(ul2) he3=DSIN(gb1) delslop=a0+a1*he1+a2*he2+a3*he3 IF(ITEC.EQ.1) THEN write(6,71)yp,xp,delslop ENDIF if(delslop.lt.slomin) slomin=delslop if(delslop.gt.slomax) slomax=delslop des=slomax-slomin dmin=999 dmax=0 j=1,m dx=(xp-x(j))/1000 dy=(yp-y(j))/1000 sp=DSQRT(dx*dx+dy*dy) if(sp.lt.dmin) dmin=sp if(sp.gt.dmax) dmax=sp if(ites.eq.1) then if(sp.gt.s0) then cs(j)=0 endif if(sp.le.s0) then z1=-sp/u2 z2=1.-z1-0.5*z1*z1 cs(j)=z2*EXP(z1) endif endif if(ites.eq.0) then z1=-sp/u2 z2=1.-z1-0.5*z1*z1 cs(j)=z2*EXP(z1) endif continue 30 l=1,m z1=0 31 k=1,m 192 31 30 32 c c 71 79 77 857 z1=z1+b(l,k)*g(k) continue f(l)=z1 continue z1=0 32 k=1,m z1=z1+cs(k)*f(k) continue up=z1 um=up/100 if(um.lt.vmin) vmin=um if(um.gt.vmax) vmax=um if(itec.eq.1) then e11=delslop-restb endif if(itec.eq.0) then e11=0 endif hieu chinh e11 la HChinh nghieng ====== ge1=ge+um+e11 ge1=ge+um write(4,71)yp,xp,ge1 format(2(f12.3),f10.4) dx=(xp-xt)/1000 dy=(yp-yt)/1000 sk=DSQRT(dx*dx+dy*dy) if(itest.eq.1) then if(itec.eq.0) then z2=h1-ge-restb-um endif if(itec.eq.1) then z2=h1-(ge+um)-delslop z3=ge+um endif e1=z2-h2 if(e1.lt.emin) emin=e1 if(e1.gt.emax) emax=e1 if(i.gt.m) then sume=sume+e1*e1 endif IF(N.EQ.M) SUME=0 if(itec.eq.1) then write(3,79)i,ten2(i),h1,ge,delslop,um,z2,e1,e11,dmin,dmax format(i3,1x,a15,f9.3,f7.3,2(f6.3),f9.3,2(f6.3),f6.1,f6.1) if(i.gt.m) then k=i-m e2=ge+um write(10,857)k,ten2(i),xp,yp,h1,ge,um,e2,h2 write(11,858)k,ten2(i),h1,ge,um,e2,z2,h2,e1 endif endif if(itec.eq.0) then z1=h1-ge-restb-um e1=z2-h2 write(3,77)i,ten2(i),h1,ge,restb,um,z2,e1,dmin,dmax format(i3,1x,a15,f9.3,f7.3,2(f6.3),f9.3,f6.3,2(f6.1)) e2=ge+um write(9,857)i,ten2(i),xp,yp,h1,ge,um,e2,h2 format(i3,1x,a15,2(f12.3),f9.3,3(f7.3),f9.3) 193 858 80 840 933 934 935 896 767 866 867 139 138 55 101 45 if(i.gt.m) then k=i-m write(10,857)k,ten2(i),xp,yp,h1,ge,um,e2,h2 write(11,858)k,ten2(i),h1,ge,um,e2,z1,h2,e1 format(i3,1x,a15,f10.3,3(f8.3),2(f9.3),f7.3) endif endif endif if(itest.eq.0) then if(itec.eq.0) then write(3,80)i,xp,yp,ge,um,ge1,dmin,dmax format(i5,1x,2(f12.3),3(f9.4),2(f6.1)) endif if(itec.eq.1) then write(3,840)i,xp,yp,ge,um,ge1,dmin,dmax,e11 format(i5,1x,2(f12.3),3(f9.4),2(f6.1),f7.4) endif endif continue write(3,933)vmin,vmax format(' Gia tri noi suy MIN:',f9.4,'; noi suy MAX:',f9.4) if(itec.eq.1) then write(3,934)slomin,slomax format(' SLOP MIN:',f9.4,'; SLOP MAX:',f9.4) write(3,935)des format(' Hieu so SLOMAX-SLOMIN:',f10.4,' m') endif if(itest.eq.1) then if(n.gt.m) then e3=DSQRT(sume/(n-m)) endif if(n.eq.m) e3=0 write(3,896)e3 format(' Sai lech trung phuong:',f8.4,' m') write(3,767)emin,emax format(' Sai(min)=',f8.4,' m; Sai(max)=',f8.4,' m') endif goto 138 write(*,*)' Thieu tep COLLO2.PAR' goto 138 write(*,*)' Thieu tep Kernel.PAR' GOTO 138 WRITE(*,*)' KHONG PHU HOP GIUA COLLO2.PAR VA KERNEL.PAR' stop end SUBROUTINE DMT2(N,b) double precision b(180,180),E1,E2 INTEGER I,J,L,N 43 l=1,n e1=0 if(l.eq.1) goto 100 45 i=1,l-1 e2=0 if(i.eq.1) goto 101 55 j=1,i-1 e2=e2+b(j,l)*b(j,i) b(i,l)=(b(i,l)-e2)/b(i,i) e1=e1+b(i,l)*b(i,l) 194 100 43 c 58 57 56 c 102 111 110 95 95 97 96 b(l,l)=DSQRT(b(l,l)-e1) continue write(*,*)' TINH XONG R=SQRT(N)' 56 l=1,n b(l,l)=1./b(l,l) if(l.eq.n) goto 102 57 j=l+1,n e1=0 58 i=l,j-1 e1=e1+b(l,i)*b(i,j) b(l,j)=-e1/b(j,j) continue write(*,*)' TINH XONG q=R-1' 95 i=1,n 110 j=i,n e1=0 111 l=j,n e1=e1+b(i,l)*b(j,l) b(i,j)=e1 b(J,I)=b(I,J) continue RETURN END SUBROUTINE AZYMUT(a1,a2,a13) double precision a1,a2,a13,pi pi=.314159265358979D+1 if(a1.ne.0.) goto 95 a13=pi/2 if(a2.lt.0.) a13=a13+pi goto 96 a13=DATAN(a2/a1) if(a1.lt.0.) goto 97 if(a2.lt.0.) a13=a13+2.*pi goto 96 a13=a13+pi return end SUBROUTINE DOIGOC(iu1,iu2,hp3,b33) double precision r1,b33,hp3 integer iu1,iu2 r1=0.206264806247096D+6 b33=(iu1*3600.+iu2*60.+hp3)/r1 return end SUBROUTINE HOANGOC(u1,is1,is2,s3) double precision w1,w2,u1,r1,s3 integer is1,is2 r1=0.206264806247096D+6 w1=u1*r1/3600 is1=int(w1) w2=(w1-is1)*60 is2=int(w2) s3=(w2-is2)*60 return end 195 SUBROUTINE CUNGKT(ra,rb,b1,ss) double precision ra,rb,b1,z1,z2,z3,z4,ea,eb,h1,h2,h3,h4,h5,ss double precision rc,cb1,cb2 z1=ra*ra z2=rb*rb ea=(z1-z2)/z1 eb=(z1-z2)/z2 z1=(3.*eb)/4 z2=(45.*eb*eb)/64 z3=(175.*eb*eb*eb)/256 z4=(11025.*eb*eb*eb*eb)/16384 h1=1.-z1+z2-z3+z4 h2=h1-1 h3=15.*eb*eb/32.-175.*eb*eb*eb/384.+3675.*eb*eb*eb*eb/8192 h4=-35.*eb*eb*eb/96.+735.*eb*eb*eb*eb/2048 h5=315.*eb*eb*eb*eb/1024 z1=1.-ea rc=ra/DSQRT(z1) cb1=DCOS(b1) cb2=cb1*cb1 z2=h1*b1+(h2+h3*cb2+h4*cb2*cb2+h5*cb2*cb2*cb2)*cb1*DSIN(b1) ss=rc*z2 return end 333 444 609 553 SUBROUTINE TINHBL(ra,rb,x2,y2,m0,b2,zl2,gama1,tyle1) double precision ra,rb,x2,y2,b2,zl2,db,b0,z1,z2,z3,ea,eb,h1,h5, * gama1,tyle1,rm2 double precision sb1,sb2,cb1,cb2,et2,t1,t2,rn,rm,s3,dk,m0,s1,s2, * s4,t4,t6,m2,m4,m6,u1,u2,u3,u4,e2,e4,e6,e8,rn2,rn4,rn6,z4,z5,eps eps=0.1D-6 z1=ra*ra z2=rb*rb ea=(z1-z2)/z1 eb=(z1-z2)/z2 z3=eb*eb h1=1.-3.*eb/4.+45.*z3/64.-175.*z3*eb/256.+11025.*z3*z3/16384 db=0 z3=ra*h1 b0=x2*SQRT(1.-ea)/(z3*m0) b0=b0+db CALL CUNGKT(ra,rb,b0,s3) dk=x2-s3*m0 if(dabs(dk).lt.eps) goto 444 db=dk*SQRT(1.-ea)/(ra*h1) goto 333 write(8,609)b0 format(' bo=',f12.9) sb1=DSIN(b0) sb2=sb1*sb1 cb1=DCOS(b0) cb2=cb1*cb1 t1=DTAN(b0) write(8,553)t1 format(' t0=',f13.7) t2=t1*t1 t4=t2*t2 t6=t4*t2 m2=m0*m0 196 551 552 c c m4=m2*m2 m6=m4*m2 e2=y2*y2 e4=e2*e2 e6=e4*e2 e8=e4*e4 s1=(1.-ea*sb2)/(1.-ea) s2=s1*s1 s4=s2*s2 rn=ra/DSQRT(1.-ea*sb2) write(8,551)rn format(' N=',f13.3) rn2=rn*rn rn4=rn2*rn2 rn6=rn4*rn2 rm=rn*(1.-ea)/(1.-ea*sb2) write(8,552)rm format(' M=',f13.3) rm2=rm*rm z1=t1/(m0*rm) z2=e2/(2.*m0*rn) u1=z1*z2 z2=e4/(24.*m0*m2*rn2*rn) z3=-4.*s2+9.*s1*(1.-t2)+12.*t2 u2=z1*z2*z3 z2=e6/(720.*m4*m0*rn4*rn) z3=8.*s4*(11.-24.*t2)-12.*s2*s1*(21.-71.*t2) z4=15.*s2*(15.-98.*t2+15.*t4) z5=180.*s1*(5.*t2-3.*t4)+360.*t4 u3=z1*z2*(z3+z4+z5) z2=e8/(40320.*m6*m0*rn6*rn) z3=1385.+3633.*t2+4095.*t4+1575*t6 u4=z1*z2*z3 b2=b0-u1+u2-u3+u4 Tinh xong vi trac dia B============= z1=1./cb1 z2=y2/(m0*rn) u1=z1*z2 z2=e2*y2/(6.*m2*m0*rn2*rn) z3=s1+2.*t2 u2=z1*z2*z3 z2=e4*y2/(120.*m4*m0*rn4*rn) z3=-4.*s1*s2*(1.-6.*t2)+s2*(9.-68.*t2)+72.*s1*t2+24.*t4 u3=z1*z2*z3 z2=e6*y2/(5040.*m6*m0*rn6*rn) z3=61.+662*t2+1320*t4+720*t6 u4=z1*z2*z3 zl2=u1-u2+u3-u4 Tinh xong kinh trac dia l'=L-Lo=========== z1=y2/(m0*rn) u1=t1*z1 z1=e2*y2/(3.*m0*m2*rn*rn2) z2=-2.*s2+3.*s1+t2 u2=t1*z1*z2 z1=e4*y2/(15.*m4*m0*rn4*rn) z2=s4*(11.-24.*t2)-3.*s1*s2*(8.-23.*t2) z3=5.*s2*(3.-14.*t2)+30.*s1*t2+3.*t4 u3=t1*z1*(z2+z3) z1=e6*y2/(315.*m6*m0*rn6*rn) 197 c  z2=17.+77.*t2+105.*t4+45.*t6 u4=t1*z1*z2 gama1=u1-u2+u3-u4 tinh xong goc hoi tu kinh tuyen =========== u1=e2/(2.*m2*rm*rn) z1=e4/(24.*m4*rm2*rn2) z2=4.*s1*(1.-6.*t2)-3.*(1.-16.*t2)-24.*t2/s1 u2=z1*z2 u3=e6/(720.*m6*rm*rm2*rn*rn2) tyle1=m0*(1.+u1+u2+u3) return end 198 Phụ lục 15: KẾT QUẢ TÍNH TỐN CHO 136 TUYẾN GIỮA 17 ĐIỂM KIỂM TRA TT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 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 DAU I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)12 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 I(BMT-APD)35 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(QK-LT)5 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(LS-BN)21 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 III(DX-DL)3 CUOI I(BMT-APD)35 III(QK-LT)5 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(QK-LT)5 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(LS-BN)21 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(DX-DL)3 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 III(KRKM-MDR)10 I(BMT-APD)1-2 III(BD-BMT)4 II(MT-TH)21 III(DM-CR)7 I(VL-HT)123 III(DC-NB)4 III(XH-SL)16 III(CH-PQ)3 III(AL-DT)1 III(MR-HT)3 III(HD-BHB)3-1 D(km) ∆hNS 81.5 -610.956 61.5 -9.685 167.1 -893.405 68.9 -394.645 118.7 -418.993 62.4 -476.784 68.0 -704.172 158.8 -781.882 161.1 -744.091 231.1 -675.155 159.5 -598.239 234.7 -410.463 247.3 -287.940 293.0 -871.137 316.0 -262.204 341.2 -744.707 73.1 601.271 215.4 -282.449 128.7 216.311 188.6 191.963 143.2 134.172 148.4 -93.216 238.1 -170.926 242.5 -133.135 311.7 -64.199 230.2 12.717 315.0 200.493 317.7 323.016 373.6 -260.181 391.6 348.752 416.5 -133.751 142.4 -883.720 63.5 -384.960 124.4 -409.308 101.2 -467.099 119.1 -694.487 183.7 -772.197 201.7 -734.406 259.3 -665.470 217.4 -588.554 279.5 -400.778 304.1 -278.255 335.9 -861.452 368.0 -252.519 393.3 -735.022 98.3 498.760 76.7 474.412 139.1 416.621 166.5 189.233 132.5 111.523 182.1 149.314 191.0 218.250 251.7 295.166 252.4 482.942 319.1 605.465 294.2 22.268 354.0 631.201 377.9 148.698 61.2 -24.348 57.4 -82.139 84.2 -309.527 121.7 -387.237 146.7 -349.446 197.2 -280.510 183.4 -203.594 225.6 -15.818 264.8 106.705 279.6 -476.492 319.8 132.441 345.1 -350.062 ∆htc -610.989 -9.734 -893.207 -394.553 -418.977 -476.772 -704.045 -781.937 -744.051 -675.286 -598.263 -410.449 -287.914 -871.110 -262.239 -744.720 601.255 -282.218 216.436 192.012 134.217 -93.056 -170.948 -133.062 -64.297 12.726 200.540 323.075 -260.121 348.750 -133.731 -883.473 -384.819 -409.243 -467.038 -694.311 -772.203 -734.317 -665.552 -588.529 -400.715 -278.180 -861.376 -252.505 -734.986 498.654 474.230 416.435 189.162 111.270 149.156 217.921 294.944 482.758 605.293 22.097 630.968 148.487 -24.424 -82.219 -309.492 -387.384 -349.498 -280.733 -203.710 -15.896 106.639 -476.557 132.314 -350.167 δh -.033 -.049 198 092 016 012 127 -.055 040 -.131 -.024 014 026 027 -.035 -.013 -.016 231 125 049 045 160 -.022 073 -.098 009 047 059 060 -.002 020 247 141 065 061 176 -.006 089 -.082 025 063 075 076 014 036 -.106 -.182 -.186 -.071 -.253 -.158 -.329 -.222 -.184 -.172 -.171 -.233 -.211 -.076 -.080 035 -.147 -.052 -.223 -.116 -.078 -.066 -.065 -.127 -.105 fgh(3) fgh(4)fgh(KT) 108 226 451 dat 094 196 392 dat 155 323 646 100 207 415 dat 131 272 545 dat 095 197 395 dat 099 206 412 151 315 630 dat 152 317 635 dat 182 380 760 dat 152 316 632 dat 184 383 766 dat 189 393 786 dat 205 428 856 dat 213 444 889 dat 222 462 924 dat 103 214 427 dat 176 367 734 136 284 567 dat 165 343 687 dat 144 299 598 dat 146 305 609 185 386 772 dat 187 389 779 dat 212 441 883 dat 182 379 759 dat 213 444 887 dat 214 446 891 dat 232 483 966 dat 237 495 989 dat 245 510 1.020 dat 143 298 597 096 199 398 134 279 558 dat 121 251 503 dat 131 273 546 163 339 678 dat 170 355 710 dat 193 403 805 dat 177 369 737 dat 201 418 836 dat 209 436 872 dat 220 458 916 dat 230 480 959 dat 238 496 992 dat 119 248 496 dat 105 219 438 142 295 590 155 323 645 dat 138 288 575 162 337 675 dat 166 345 691 190 397 793 191 397 794 dat 214 447 893 dat 206 429 858 dat 226 470 941 233 486 972 dat 094 196 391 dat 091 189 379 dat 110 229 459 dat 132 276 552 145 303 606 dat 169 351 702 163 339 677 dat 180 376 751 dat 195 407 814 dat 201 418 836 dat 215 447 894 dat 223 464 929 dat KL III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III III 199 71 III(KRKM-MDR)10I(BMT-APD)1-2 70.8 72 III(KRKM-MDR)10III(BD-BMT)4 95.3 73 III(KRKM-MDR)10II(MT-TH)21 68.8 74 III(KRKM-MDR)10III(DM-CR)7 109.8 75 III(KRKM-MDR)10I(VL-HT)123 141.5 76 III(KRKM-MDR)10III(DC-NB)4 175.0 77 III(KRKM-MDR)10III(XH-SL)16 185.4 78 III(KRKM-MDR)10III(CH-PQ)3 244.0 79 III(KRKM-MDR)10III(AL-DT)1 234.0 80 III(KRKM-MDR)10III(MR-HT)3 285.3 81 III(KRKM-MDR)10III(HD-BHB)3-1 309.9 82 I(BMT-APD)1-2 III(BD-BMT)4 28.1 83 I(BMT-APD)1-2 II(MT-TH)21 96.9 84 I(BMT-APD)1-2 III(DM-CR)7 101.7 85 I(BMT-APD)1-2 I(VL-HT)123 168.7 86 I(BMT-APD)1-2 III(DC-NB)4 126.0 87 I(BMT-APD)1-2 III(XH-SL)16 178.6 88 I(BMT-APD)1-2 III(CH-PQ)3 208.4 89 I(BMT-APD)1-2 III(AL-DT)1 235.6 90 I(BMT-APD)1-2 III(MR-HT)3 267.5 91 I(BMT-APD)1-2 III(HD-BHB)3-1 292.9 92 III(BD-BMT)4 II(MT-TH)21 105.4 93 III(BD-BMT)4 III(DM-CR)7 95.0 94 III(BD-BMT)4 I(VL-HT)123 170.8 95 III(BD-BMT)4 III(DC-NB)4 100.6 96 III(BD-BMT)4 III(XH-SL)16 166.8 97 III(BD-BMT)4 III(CH-PQ)3 185.5 98 III(BD-BMT)4 III(AL-DT)1 225.2 99 III(BD-BMT)4 III(MR-HT)3 249.6 100 III(BD-BMT)4 III(HD-BHB)3-1 274.9 101 II(MT-TH)21 III(DM-CR)7 54.2 102 II(MT-TH)21 I(VL-HT)123 75.6 103 II(MT-TH)21 III(DC-NB)4 141.8 104 II(MT-TH)21 III(XH-SL)16 120.0 105 II(MT-TH)21 III(CH-PQ)3 192.1 106 II(MT-TH)21 III(AL-DT)1 165.6 107 II(MT-TH)21 III(MR-HT)3 221.6 108 II(MT-TH)21 III(HD-BHB)3-1 245.5 109 III(DM-CR)7 I(VL-HT)123 82.3 110 III(DM-CR)7 III(DC-NB)4 92.2 111 III(DM-CR)7 III(XH-SL)16 78.9 112 III(DM-CR)7 III(CH-PQ)3 138.4 113 III(DM-CR)7 III(AL-DT)1 134.2 114 III(DM-CR)7 III(MR-HT)3 175.9 115 III(DM-CR)7 III(HD-BHB)3-1 200.9 116 I(VL-HT)123 III(DC-NB)4 167.0 117 I(VL-HT)123 III(XH-SL)16 86.0 118 I(VL-HT)123 III(CH-PQ)3 181.3 119 I(VL-HT)123 III(AL-DT)1 107.3 120 I(VL-HT)123 III(MR-HT)3 180.1 121 I(VL-HT)123 III(HD-BHB)3-1 200.6 122 III(DC-NB)4 III(XH-SL)16 112.6 123 III(DC-NB)4 III(CH-PQ)3 88.0 124 III(DC-NB)4 III(AL-DT)1 167.1 125 III(DC-NB)4 III(MR-HT)3 163.4 126 III(DC-NB)4 III(HD-BHB)3-1 187.6 127 III(XH-SL)16 III(CH-PQ)3 97.3 128 III(XH-SL)16 III(AL-DT)1 58.9 129 III(XH-SL)16 III(MR-HT)3 101.8 130 III(XH-SL)16 III(HD-BHB)3-1 125.5 131 III(CH-PQ)3 III(AL-DT)1 125.1 132 III(CH-PQ)3 III(MR-HT)3 84.8 133 III(CH-PQ)3 III(HD-BHB)3-1 105.8 134 III(AL-DT)1 III(MR-HT)3 81.7 135 III(AL-DT)1 III(HD-BHB)3-1 97.3 136 III(MR-HT)3 III(HD-BHB)3-1 25.4 Dat:112 Khong dat: 24 Dmin= 25.4 km ; Dmax= 416.5 km m1= 0000m/km; m2= 0088 m/km -57.791 -285.179 -362.889 -325.098 -256.162 -179.246 8.530 131.053 -452.144 156.789 -325.714 -227.388 -305.098 -267.307 -198.371 -121.455 66.321 188.844 -394.353 214.580 -267.923 -77.710 -39.919 29.017 105.933 293.709 416.232 -166.965 441.968 -40.535 37.791 106.727 183.643 371.419 493.942 -89.255 519.678 37.175 68.936 145.852 333.628 456.151 -127.046 481.887 -.616 76.916 264.692 387.215 -195.982 412.951 -69.552 187.776 310.299 -272.898 336.035 -146.468 122.523 -460.674 148.259 -334.244 -583.197 25.736 -456.767 608.933 126.430 -482.503 -57.795 -285.068 -362.960 -325.074 -256.309 -179.286 8.528 131.063 -452.133 156.738 -325.743 -227.273 -305.165 -267.279 -198.514 -121.491 66.323 188.858 -394.338 214.533 -267.948 -77.892 -40.006 28.759 105.782 293.596 416.131 -167.065 441.806 -40.675 37.886 106.651 183.674 371.488 494.023 -89.173 519.698 37.217 68.765 145.788 333.602 456.137 -127.059 481.812 -.669 77.023 264.837 387.372 -195.824 413.047 -69.434 187.814 310.349 -272.847 336.024 -146.457 122.535 -460.661 148.210 -334.271 -583.196 25.675 -456.806 608.871 126.390 -482.481 -.004 111 -.071 024 -.147 -.040 -.002 010 011 -.051 -.029 115 -.067 028 -.143 -.036 002 014 015 -.047 -.025 -.182 -.087 -.258 -.151 -.113 -.101 -.100 -.162 -.140 095 -.076 031 069 081 082 020 042 -.171 -.064 -.026 -.014 -.013 -.075 -.053 107 145 157 158 096 118 038 050 051 -.011 011 012 013 -.049 -.027 001 -.061 -.039 -.062 -.040 022 101 117 100 126 143 159 163 187 184 203 211 064 118 121 156 135 160 173 184 196 205 123 117 157 120 155 163 180 190 199 088 104 143 131 166 154 179 188 109 115 107 141 139 159 170 155 111 162 124 161 170 127 113 155 153 164 118 092 121 134 134 110 123 108 118 060 210 244 207 262 297 331 340 390 382 422 440 133 246 252 325 281 334 361 384 409 428 257 244 327 251 323 340 375 395 415 184 217 298 274 347 322 372 392 227 240 222 294 290 332 354 323 232 337 259 335 354 265 235 323 320 342 247 192 252 280 280 230 257 226 247 126 421 488 415 524 595 661 681 781 765 845 880 265 492 504 649 561 668 722 767 818 856 513 487 653 501 646 681 750 790 829 368 435 595 548 693 643 744 783 454 480 444 588 579 663 709 646 464 673 518 671 708 531 469 646 639 685 493 384 505 560 559 460 514 452 493 252 dat dat dat dat III III III III dat dat dat dat dat dat III III III III III III dat dat dat dat dat dat dat dat dat III III III III III III III III III dat III dat dat dat dat dat III III III III III dat dat dat dat dat dat dat III III III III III III III dat dat dat dat dat dat dat III III III III III III III dat III dat dat dat dat dat dat dat dat dat dat dat dat dat dat dat dat dat III III III III III III III III III III III III III III III III III ... dụng công nghệ đo cao GPS số ứng dụng khác Trong bối cảnh chúng tơi lựa chọn đề tài ? ?Nghiên cứu xác hóa dị thường độ cao EGM2 008 dựa số liệu GPS - Thủy chuẩn phạm vi cục Vi? ??t Nam? ?? với mục tiêu... Geoid (độ cao chính) so với mặt Quasigeoid (độ cao chuẩn) Trên giới nói chung Vi? ??t Nam nói riêng, sử dụng độ cao chuẩn Chỉ có số nước sử dụng độ cao Vi? ??c tính độ cao Geoid N hay dị thường độ cao. .. OSU91A, EGM9 6, EGM2 008 vv… để xử lý lưới GPS, xác định độ cao thủy chuẩn cho điểm lưới Tuy nhiên, độ xác dị thường độ cao (hoặc độ cao Geoid ) hiệu dị thường độ cao (hoặc hiệu độ cao Geoid) xác định