Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình MỤC LỤC SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình LỜI MỞ ĐẦU - - - Trắc địa môn khoa học phát triển từ lâu, trình phát triển người không ngừng ứng dụng thành tựu khoa học kỹ thuật để hoàn thiện phương pháp thiết bị đo nhằm giảm bớt khó khăn công tác đo đạc đồng thời nâng cao độ xác thành trắc địa Trong khoảng 20 năm trở lại đây, thiết bị đo đạc đươc đổi chất Đó giai đoạn chuyển dần từ thiết bị quang sang thiết bị đo đạc điện tử, thời gian người ta ứng dụng khoa học vũ trụ vào lĩnh vực trắc địa Mỹ xây dựng hệ thống định vị toàn cầu (GPS), Nga xây dựng hệ thống vệ tinh dẫn đường toàn cầu riêng cho hệ thống GLONASS (hệ thống gồm 24 vệ tinh, bay mặt phẳng quĩ đạo độ cao từ 18800 đến 19900km, với chu kỳ 676 phút), nước châu Âu xây dựng hệ thống định vị toàn cầu GALILEO Cũng xin lưu ý, thời gian gần đây, cộng đồng quốc tế sử dụng Global Navigation Satellite System (GNSS-Hệ thống vệ tinh dẫn đường toàn cầu) thuật ngữ chung đề cập tới hệ thống vệ tinh hoạt động theo nguyên lý định vị Ngày nay, hệ thống GNSS ứng dụng rộng rãi trắc địa, nhiều lĩnh vực khác : Ứng dụng dẫn đường: với sai số vị trí điểm cỡ m (hoặc chục m), nhằm dẫn đường cho tàu thuyền biển, phương tiện không (máy bay, tên lửa), lĩnh vực quân sự, giao thông vận tải, du lịch … Ứng dụng trắc địa: với sai số xác định vị trí điểm từ mm đến dm hoàn toàn đáp ứng yêu cầu kỹ thuật toàn nhiệm vụ trắc địa như: xây dựng lưới khống chế, đo vẽ đồ, mặt cắt, định vị tàu thuyền đo đạc biển, ứng dụng trắc địa công trình (thành lập mạng lưới sở trắc địa công trình lưới thi công công trình v.v ) Ứng dụng nghiên cứu biến đổi môi trường toàn cầu: Với ưu vượt trội (độ xác cao, hoạt động liên tục quy mô toàn cầu) GNSS áp dụng xác định chuyển dịch đại vỏ Trái đất, nghiên cứu tầng khí (điện ly, khí tượng) tượng nước biển dâng Ở nước ta, công nghệ GPS đưa vào từ năm 1990, sau thời gian ngắn, sử dụng GPS để hoàn chỉnh mạng lưới thiên văn- trắc địa quốc gia nhiều ứng dụng khác có ích cho quốc gia.Ví dụ từ năm 1991 đến 1993 : Xây dựng lưới GPS hệ ngắn cho khu vực Minh Hải-Sông Bé- Tây Nguyên (gồm 117mốc).Năm 1992: Xây dựng lưới GPS hệ dài biển đất liền, lưới gồm 36 mốc Năm 1995: Xây dựng lưới cấp 0, gồm 71 mốc phủ toàn lãnh thổ Việt Nam làm sở để xây dựng hệ toạ độ VN-2000.vv… Từ đến công SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình nghệ GPS ứng dụng rộng rãi để xây dựng nhiều lưới, việc ứng dụng GPS để thành lập lưới khống chế phần quan trọng Xử lý số liệu khâu then chốt ứng dụng GNSS Với khối lượng số liệu khổng lồ (dung lượng ca đo 24 với tần xuất 15 giây điểm in ~ 1000 tờ A4), mô hình toán học vật lý phức tạp, số liệu GNSS xử lý phần mềm chuyên nghiệp Trên giới, công ty, viện, trường đại học tung nhiều phần mềm khác nhau, song chia làm hai loại: 1.phần mềm mang tính thương mại (TGO, GPS survey TTC, TBC…); 2.phần mềm chuyên dụng mang tính khoa học (Bernese, Gamit, Gipsy-Oasis…) Phần mềm mang tính thương mại thường cung cấp theo máy thu, việc sử dụng không phức tạp, đáp ứng nhiệm vụ trắc địa thường gặp phục vụ cho đông đảo khách hàng Ngược lại, phần mềm chuyên dụng mang tính khoa học thiết lập nhằm xử lý số liệu cho nhiệm vụ đặc biệt xác định chuyển dịch vỏ Trái đất Để đáp ứng yêu cầu độ xác cao quy mô không gian lớn (sai số cỡ mm không gian hàng ngàn km) phần mềm yêu cầu hàng loạt số liệu hỗ trợ trước xử lý Với khuôn khổ đồ án tốt nghiệp, khó xử lý số liệu phần mềm Bernese, đề tài với tiêu đề: “Chuẩn bị số liệu hỗ trợ cho công tác xử lý số liệu GPS độ xác cao phần mềm Bernese” lựa chọn Đây hội thuận lợi, không làm quen với công tác nghiên cứu mà tiếp cận thực tế, góp phần giải nhiệm vụ phức tạp Khái niệm số liệu hiểu số liệu đo điểm IGS, số liệu chung liên quan tới hệ thống GNSS, phần mềm (thuộc thư mục GEN) số liệu hỗ trợ xử lý số liệu đo (lịch vệ tinh, mô hình điện ly, tải trọng đại dương …) Với tinh thần trên, phần mở đầu kết luận, nôi dung đồ án gồm: Chương 1: Xây dựng lưới khống chế công nghệ GPS Chương 2: Quy trình xử lý số liệu GPS độ xác cao Chương 3: Chuẩn bị số liệu bổ trợ cho công tác xử lý số liệu GPS độ xác cao phần mềm Bernese Em xin chân thành cảm ơn TS Phạm Thị Hoa TS Vy Quốc Hải thầy cô môn Trắc địa cao cấp tận tình giúp đỡ để em hoành thành đồ án Do trình độ, kinh nghiệm có hạn nên tránh khỏi sai sót, em mong bảo vào đóng góp ý kiến thầy, cô bạn sinh viên hoàn thiện đồ án Em xin chân thành cảm ơn! Chương 1: XÂY DỰNG LƯỚI KHỐNG CHẾ BẰNG CÔNG NGHỆ GPS SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình 1.1.KHÁI QUÁT VỀ HỆ THỐNG GPS 1.1.1 Vài nét lịch sử phát triển Tháng 10 năm 1957, Liên Xô (cũ) phóng thành công công nghệ vệ tinh nhân tạo (VTNT) Trái Đất Sputnhic – I lên quỹ đạo, mở đầu cho kỷ nguyên người chinh phục không gian vũ trụ Cũng từ đó, người sử dụng VTNT vào giải nhiệm vụ nhiều lĩnh vực khác nhau, có Trắc địa cao cấp Ngay từ năm 1958, số liệu quan sát từ vệ tinh, người ta xác định độ dẹt cực Trái Đất f = 1/298,3 gần với kết từ đo đạc mặt đất Cũng thời gian này, quan sát VTNT người ta xác định số tham số vật lý Trái Đất Vào năm 1960, VTNT đưa lên quỹ đạo đóng vai trò mục tiêu cao, dùng thiết bị quang học để quan sát vệ tinh từ mặt đất phục vụ xây dựng lưới tam giác vệ tinh (còn gọi mạng lưới tam giác vũ trụ) cho phép chuyền tọa độ điểm cách xa bề mặt mặt đất Phương pháp tam giác vệ tinh chịu ảnh hưởng đáng kể điều kiện thời tiết, máy móc, thiết bị quan sát nặng nề, không thuận tiện cho công tác đo đạc thực địa Vì vậy, sau hoàn thành nhiệm vụ lịch sử nói trên, phương pháp không áp dụng thực tế Để khắc phục nhược điểm phương pháp tam giác vệ tinh, người ta thiết kế hệ thống đạo hàng vệ tinh làm việc điều kiện thời tiết liên tục 24 ngày Năm 1962, Mỹ thiết kế xây dựng hệ thống đạo hàng hải quân NNSS, gọi Transit Cũng thời gian này, Liên Xô (cũ) xây dựng hệ thống TSICADA có tính tương tự hệ thống Transit Mỹ Cả hai hệ thống hoạt động theo nguyên lý hiệu ứng Doppler, dựa theo tín hiệu từ vệ tinh phát xuống mặt đất Trong trường hợp vệ tinh đóng vai trò điểm gốc (có tọa độ), phương tiện truyền thông tin quỹ đạo vệ tinh, tạo trị đo Doppler để cung cấp cho máy thu thực toán định vị biển mặt đất Từ năm 1967, phương pháp Doppler vệ tinh không đột phá cho nhiệm vụ định vị biển mà mở khả xây dựng lưới khống chế tọa độ cho số quốc gia trước thập niên 80 kỷ trước Tuy vậy, hệ thống Transit có nhược điểm không đáp ứng yêu cầu định vị tức thời cần độ xác cao SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Năm 1973, hệ thống GPS thiết kế Từ năm 1978 đến 1985, người ta đưa lên quỹ đạo 11 vệ tinh khối I (block I) mang tính chất thực nghiệm Từ năm 1989 đến 1990, người ta đưa lên quỹ đạo vệ tinh thuộc khối II (block II) Các vệ tinh khối II khác vệ tinh khối I chỗ phát tín hiệu có nhiễu cố ý SA có kỹ thuật bảo mật AS Từ năm 1990 đến năm 1994, người ta đưa lên quỹ đạo 15 vệ tinh hệ II – A có khả liên hệ vệ tinh Từ sau năm 1995, hệ thống GPS tiếp tục trì bảo dưỡng thay vệ tinh già tuổi Năm 2000, số vệ tinh chòm GPS tăng lên 28 vệ tinh Những vệ tinh hệ GPS-IIR phóng lên để thay vệ tinh già tuổi Vệ tinh phóng lên ngày 16/9/2005 mang tên GPS-IIR-M1, vệ tinh thuộc hệ vệ tinh đại GPS-IIR-M Theo kế hoạch, vệ tinh phóng lên không gian vào tháng giêng năm (2006) Ngoài hệ thống GPS Mỹ, năm 1980, Liên Xô (cũ) triển khai xây dựng hệ thống định vị toàn cầu quân có tên gọi GLONASS Nguyên lý hoạt động hệ thống tương tự hệ thống GPS Để tăng cường độ xác định vị GPS GLONASS, từ cuối năm 2002, dịch vụ dẫn đường sử dụng vệ tinh phủ trùm Châu Âu EGNOS cung cấp khả định vị xác toàn Châu Âu vùng lân cận Để tăng cường độ xác cho hệ thống GPS, Mỹ xây dựng hệ thống định vị tăng cường diện rộng WAAS, Nhật Bản xây dựng hệ thống MSAS có tính tương tự WAAS EGNOS Các hệ thống cung cấp khả định vị tức thời toàn vùng phủ sóng với sai số không lớn 3m Từ tháng năm 2002, Liên minh Châu Âu bắt đầu đưa lên quỹ đạo vệ tinh hệ thống định vị toàn cầu GALILEO Hệ thống GALILEO đưa vào hoạt động thử nghiệm từ năm 2008 dự kiến hoàn thành vào năm 2013 2015 Năm 2007, Trung Quốc phát triển hệ thống định vị khu vực Bắc Đẩu – thành hệ thống định vị toàn cầu với tên gọi COMPASS hay Bắc Đẩu – Các hệ thống định vị vệ tinh GPS, GALILEO, GLONASS, COMPASS gọi chung hệ thống định vị dẫn đường toàn cầu GPS Để phục vụ cho công tác định vị điều kiện hạn chế thu liên tục tín hiệu từ vệ tinh, người ta SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình chế tạo thiết bị định vị tích hợp công nghệ định vị với hệ thống định vị quán tính INS Từ năm 2000 trở lại đây, người ta đưa vệ tinh CHAMP, GRACE, GOCE lên quỹ đạo phục vụ quan trắc gradien trọng lực trái đất với độ phân giải cao 1.1.2 Hệ thống GPS a Cấu trúc hệ thống GPS Toàn phần cứng hệ thống GPS có tên đầy đủ NAVSTAR GPS SYSTEM NAVSTAR viết tắt chữ NAVIGATION SYSTEM WITH TIME AND RANGING Phần cứng gồm phần: phần điều khiển (Control Segment), phần không gian (Space Segment) phần sử dụng (User Segment) Hình 1.1.1- Sơ đồ hoạt động hệ thống định vị toàn cầu GPS SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình b Phần điều khiển (Control Segment): Hình 1.1.2- Trạm điều khiển hệ thống GPS Phần điều khiển gồm trạm mặt đất có trạm theo dõi (Monitor Station): Diego Garcia, Ascension, Kwajalein Hawai; trạm điều khiển trung tâm (Master Control Station) trạm hiệu chỉnh số liệu (Upload Station) Lưới trắc địa đặt trạm xác định phương pháp giao thoa đường đáy dài (VLBI) Trạm trung tâm làm nhiệm vụ tính toán lại tọa độ vệ tinh theo số liệu trạm theo dõi thu từ vệ tinh Sau tính toán, số liệu gửi từ trạm trung tâm tới trạm hiệu chỉnh số liệu từ gửi tiếp tới vệ tinh Như vậy, vòng vệ tinh có số liệu hiệu chỉnh để phát cho máy thu c Phần không gian (Space Segment): SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Hình 1.1.3- Phân bố vệ tinh quỹ đạo Hình 1.1.4- Vệ tinh GPS * Chòm vệ tinh GPS: Bao gồm 24 vệ tinh bay quỹ đạo có độ cao đồng 20 200 km, chu kỳ 12 giờ, phân phối mặt phẳng quỹ đạo nghiêng với xích đạo góc 55 Việc bố trí nhằm mục đích để thời điểm vị trí trái đất quan sát vệ tinh Mỗi vệ tinh phát tần số sóng mang với tần số cao L1=1575.42 MHz L2=1227.60 MHz Loại sóng phát sở dãy số tựa ngẫu nhiên bao gồm số Mã gọi tên mã P (Precise) Bên cạnh mã P, sóng mang mã C/A (Clear/Acquisition) sóng L1 Mã C/A phát với tần số 10.23 MHz 1.023 MHz Ngoài mã vệ tinh phát mã phụ có tần số 50 Hz chứa thông tin lịch vệ tinh Các vệ tinh trang bị đồng hồ nguyên tử với độ xác cao Các vệ tinh NAVSTAR có trạng thái: "hoạt động khỏe" ( Healthy) "hoạt động không khoẻ ( Unhealthy) Hai trạng thái vệ tinh định trạm điều khiển mặt đất Chúng ta sử dụng tín hiệu vệ tinh hai trạng thái "hoạt động khỏe" "hoạt động không khỏe" d Phần sử dụng (User Segment): Phần sử dụng bao gồm máy thu tín hiệu từ vệ tinh đất liền, máy bay tàu thủy Các máy thu phân làm loại: máy thu tần số máy thu tần số Máy SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình thu tần số nhận mã phát với sóng mang L1 Các máy thu tần số nhận sóng mang L1 L2 Các máy thu tần số phát huy tác dụng đo tọa độ tuyệt độ xác 10 m tọa độ tương độ xác từ đến cm khoảng cách nhỏ 50 km Với khoảng cách lớn 50 km độ xác giảm đáng kể (độ xác cỡ dm) Để đo khoảng cách dài đến vài nghìn km phải sử dụng máy tần số để khử ảnh hưởng tầng ion khí trái đất Toàn phần cứng GPS hoạt động hệ thống tọa độ WGS-84 với kích thước elipsoid a=6378137.0 m α=1:298,5722 e Nguyên lý hoạt động hệ thống GPS Như biết nguyên lý hoạt động hệ thống DOPPLER, nguyên lý thay đổi tần số tín hiệu nơi phát tín hiệu chuyển động Hệ thống GPS hoạt động nguyên lý hoàn toàn khác Để xác định tọa độ tuyệt đối điểm mặt đất sử dụng kỹ thuật "tựa khoảng cách" Kỹ thuật mô tả công thức: C.t + C.∆ t = ( x s − x p ) + ( y s − y p ) + ( z s − z p ) (1) đây: s =[xs ys zs] - Tọa độ vệ tinh; p=[xp yp zp] - Tọa độ điểm mặt đất; c - Tọa độ sóng; t - Thời gian sóng từ vệ tinh tới máy thu ∆t - Số hiệu chỉnh thời gian Tập hợp phương trình đo dạng (1) ta có hệ thống phương trình sai số có ẩn số t, xp yp zp xs ys zs biết từ mã lịch vệ tinh (tần số 50Hz), t xác định theo đồng hồ vệ tinh máy thu theo mã C/A, c số tốc độ truyền sóng điện từ Theo kỹ thuật xác định tọa độ với độ xác 10 m Nếu kết gửi tới trạm điều khiển trung tâm, có tọa độ tuyệt đối mặt đất với độ xác m Sở dĩ độ xác tăng lên đáng kể máy thu thu lịch vệ tinh dự báo, trạm điều khiển trung tâm có lịch vệ tinh xác Qua thấy tọa độ tuyệt đối điểm mặt đất xác định có độ xác phương pháp DOPPLER Sở dĩ vệ tinh hệ thống GPS có SV: Trần Đại An Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình độ cao gấp đôi hệ thống DOPPLER Tọa độ tuyệt độ xác 10 m hệ thống GPS dùng để đáp ứng mục đích: - Đạo hàng ( định vị cho đối tượng chuyển động tàu biển, máy bay ) - Cung cấp tọa độ gần cho phương pháp đo tọa độ tương đối GPS Ngược lại với độ xác tọa độ tuyệt đối, công nghệ GPS đạt thành tựu đáng kể việc xác định tọa độ tương đối Nguyên lý đo tọa độ tương đối xác định pha sóng mang L1 (với máy thu tần số) hay L1 L2 (với máy thu tần số) Chúng ta có công thức: S = Nλ + ϕλ (2) Trong đó: λ - Bước sóng (λ = c/f) f: Tần số sóng; N: Số nguyên lần bước sóng; ϕ: Pha sóng; S: Khoảng cách vệ tinh - máy thu Từ công thức (2) có: ϕ = (f/c).S - N (3) Xét công thức (3) từ phía khác viết: ϕ(t) = φs(ts ) - φp(t) + Nsp (4) φs(ts ) - Pha sóng thời điểm ts vệ tinh bắt đầu phát tín hiệu; φp(t) - Pha sóng thời điểm t máy thu nhận tín hiệu; Nsp - Số nguyên lần bước sóng Từ công thức ta suy ra: ϕ(t) = φs(t) - (f/c).Ssp - φp(t) + Nsp (5) Kết hợp thành phần vế phải công thức (5) biểu diễn dạng: ϕ(t) = - (f/c).Ssp - αp(t) + βs(t) + γsp (6) Trong đó: αp(t) - Thành phần ảnh hưởng hệ thống pha (t) máy thu gây (chủ yếu số hiệu chỉnh đồng hồ máy thu) SV: Trần Đại An 10 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Observation Data RINEX VERSION / TYPE CS Pre 4.3 (5) Leica Inc Sun 02 Jun 2013 PGM / RUN BY / DATE TCMS MARKER NAME 23604S002 MARKER NUMBER CMS ITRI OBSERVER / AGENCY 80057 LEICA RS500 5.10 REC # / TYPE / VERS 720 LEIAT504 NONE ANT # / TYPE -2982782.2543 4966650.1111 2658815.4481 APPROX POSITION XYZ 0.0000 0.0000 0.0000 ANTENNA: DELTA H/E/N 1 WAVELENGTH FACT L1/2 L1 C1 L2 P2 # / TYPES OF OBSERV 30 INTERVAL 16 LEAP SECONDS 2013 0 0.000000 TIME OF FIRST OBS 2013 23 59 30.000000 TIME OF LAST OBS 30 # OF SATELLITES 728 728 725 725 PRN / # OF OBS 975 975 964 964 PRN / # OF OBS 837 837 794 794 PRN / # OF OBS 977 977 935 935 PRN / # OF OBS 987 987 960 960 PRN / # OF OBS 716 716 700 700 PRN / # OF OBS 853 853 828 828 PRN / # OF OBS 959 959 894 894 PRN / # OF OBS 1037 1037 961 961 PRN / # OF OBS 10 838 838 789 789 PRN / # OF OBS 11 780 780 729 729 PRN / # OF OBS 12 1085 1085 1041 1041 PRN / # OF OBS 13 1040 1040 935 935 PRN / # OF OBS 14 925 925 851 851 PRN / # OF OBS 15 917 917 915 915 PRN / # OF OBS 16 956 956 905 905 PRN / # OF OBS 17 895 895 885 885 PRN / # OF OBS 18 839 839 747 747 PRN / # OF OBS 19 848 848 834 834 PRN / # OF OBS 20 951 951 795 795 PRN / # OF OBS 21 847 847 673 673 PRN / # OF OBS 22 758 758 740 740 PRN / # OF OBS 23 1011 1011 959 959 PRN / # OF OBS 24 708 708 704 704 PRN / # OF OBS 25 871 871 754 754 PRN / # OF OBS 26 710 710 686 686 PRN / # OF OBS 28 761 761 715 715 PRN / # OF OBS 29 1012 1012 998 998 PRN / # OF OBS 31 1015 1015 994 994 PRN / # OF OBS 32 1032 1032 962 962 PRN / # OF OBS BIT OF LLI FLAGS DATA COLLECTED UNDER A/S CONDITION COMMENT SNR is mapped to RINEX snr flag value [2-9] COMMENT L1&L2: < 25dBHz -> 1; 25-27dBHz -> 2; 28-31dBHz -> COMMENT 32-35dBHz -> 4; 36-38dBHz -> 5; 39-41dBHz -> COMMENT 42-44dBHz -> 7; 45-48dBHz -> 8; 49-55dBHz -> COMMENT END OF HEADER 13 00 00 0.0000000 12 14 15 18 21 22 24 25 31 114324001.78908 21755134.207 89083638.47609 21755134.792 122536221.70347 23317862.214 95482785.29347 23317863.346 122778631.22307 23363997.105 95671688.28908 23363999.119 107159035.85248 20391681.618 83500529.66849 20391680.927 114998547.39247 21883500.766 89609239.92848 21883500.903 SV: Trần Đại An 49 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình 2.11 OBSERVATION DATA M (MIXED) RINEX VERSION / TYPE teqc 2013May1 NOAA/NOS/NGS/CORS 20130612 23:50:38UTCPGM / RUN BY / DATE Solaris x86 5.10|AMD64|cc SC5.8 -xarch=amd64|=+|=+ COMMENT teqc 2013May1 NOAA/NOS/NGS/CORS 20130612 23:50:36UTCCOMMENT 0000.000 (antenna height) COMMENT BIT OF LLI FLAGS DATA COLLECTED UNDER A/S CONDITION COMMENT WUHN MARKER NAME 21602M001 MARKER NUMBER Prof Guigen Nie WHU OBSERVER / AGENCY 4918K35018 TRIMBLE NETR8 NP 4.16 / SP 3.59 REC # / TYPE / VERS 4852A62817 TRM59800.00 ENCL ANT # / TYPE -2267750.2660 5009156.1420 3221291.8980 APPROX POSITION XYZ 2.3612 -0.0022 -0.0094 ANTENNA: DELTA H/E/N 1 WAVELENGTH FACT L1/2 11 L1 L2 L5 C1 P1 C2 P2 C5 S1# / TYPES OF OBSERV S2 S5 # / TYPES OF OBSERV 30.0000 INTERVAL SNR is mapped to RINEX snr flag value [0-9] COMMENT L1 & L2: min(max(int(snr_dBHz/6), 0), 9) COMMENT teqc windowed: start @ 2013 Jun 00:00:00.000 COMMENT teqc windowed: delta = 86400.000 sec COMMENT teqc edited: all SBAS satellites excluded COMMENT teqc edited: all Galileo satellites excluded COMMENT teqc edited: all Compass satellites excluded COMMENT teqc edited: all QZSS satellites excluded COMMENT 16 LEAP SECONDS Forced Modulo Decimation to 30 seconds COMMENT 2013 0 0.0000000 GPS TIME OF FIRST OBS END OF HEADER 13 0 0.0000000 14G24R10G25G18R11G22R08G21R06G12G15R09 G14R07 113142468.552 88163077.210 21530313.883 21530324.328 21530324.379 40.300 48.600 102570724.334 79777268.160 19242010.473 19242003.188 53.500 19242002.484 44.100 121788466.17715 23175596.332 94900194.05116 23175596.777 23175583.531 31.800 41.400 106473205.057 82966321.22144 20261225.113 20261218.313 43.100 28.300 117865603.042 91673259.121 22056944.227 22056937.008 45.100 22056936.316 36.400 112592863.727 87734792.76544 21425728.234 21425722.273 40.700 24.700 110794036.660 86173148.766 20690004.480 20689999.484 49.900 20689998.652 45.300 115826524.504 90254476.74342 22041071.359 22041064.195 37.300 16.900 124504881.220 96837148.243 23332166.262 23332160.008 37.300 23332158.758 32.500 117868156.412 22429566.195 91845334.927 22429566.063 22429557.906 35.400 41.300 124588414.482 23708417.242 97081987.541 23708417.805 23708407.328 35.200 41.100 111098380.662 86410091.088 20805230.047 20805224.313 50.400 20805223.570 45.600 3.3.2 Lịch vệ tinh xác Lịch vệ tinh xác gồm tệp IGS17426.SP3, IGS17430.SP3, IGS17431.SP3, IGS17432.SP3 phần số liệu tệp trình bày bảng 3.4 SV: Trần Đại An 50 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Bảng 3.4 Một phần số liệu lịch tinh xác #cP2013 0 0.00000000 96 ORBIT IGb08 HLM IGS ## 1742 518400.00000000 900.00000000 56444 0.0000000000000 + 32 G01G02G03G04G05G06G07G08G09G10G11G12G13G14G15G16G17 + G18G19G20G21G22G23G24G25G26G27G28G29G30G31G32 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 ++ 2 2 2 2 2 2 2 2 ++ 2 2 2 2 2 2 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 %c G cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %f 1.2500000 1.025000000 0.00000000000 0.000000000000000 %f 0.0000000 0.000000000 0.00000000000 0.000000000000000 %i 0 0 0 0 %i 0 0 0 0 /* RAPID ORBIT COMBINATION FROM WEIGHTED AVERAGE OF: /* cod emr esa gfz jpl ngs sio usn whu /* REFERENCED TO IGS TIME (IGST) AND TO WEIGHTED MEAN POLE: /* PCV:IGS08_1740 OL/AL:FES2004 NONE Y ORB:CMB CLK:CMB * 2013 0 0.00000000 PG01 13512.466402 -12220.845980 19279.970808 32.308429 7 115 PG02 -13529.778371 -14013.560316 -18060.113029 434.861765 110 PG03 25108.907229 8764.536122 1594.215224 214.316866 127 PG04 -8436.586785 -22468.562303 -11354.726472 233.348709 98 PG05 -20342.363349 -2295.696102 -17022.732229 -396.004817 130 PG06 22958.912550 13487.166525 -2613.173916 207.731072 10 127 PG07 6273.812878 -24171.053768 -8545.782808 220.061714 9 97 PG08 1376.666647 -26096.325695 2373.243135 6.485850 10 119 PG09 -7626.593013 -22427.261286 11354.841038 244.723538 10 133 PG10 -5743.405798 -14726.057211 -21612.033966 -81.866456 108 PG11 15851.970057 -4685.557838 20350.563067 -387.283203 115 PG12 -24071.361244 10443.179863 3252.210215 136.854077 9 95 PG13 8583.607621 -12083.936613 -22180.457305 113.877300 11 132 PG14 10707.341363 19440.529945 14902.047726 220.656683 106 SV: Trần Đại An 51 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình #cP2013 0 0.00000000 96 ORBIT IGb08 HLM IGS ## 1743 0.00000000 900.00000000 56445 0.0000000000000 + 32 G01G02G03G04G05G06G07G08G09G10G11G12G13G14G15G16G17 + G18G19G20G21G22G23G24G25G26G27G28G29G30G31G32 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 ++ 2 2 2 2 2 2 2 ++ 2 2 2 2 2 2 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 %c G cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %f 1.2500000 1.025000000 0.00000000000 0.000000000000000 %f 0.0000000 0.000000000 0.00000000000 0.000000000000000 %i 0 0 0 0 %i 0 0 0 0 /* RAPID ORBIT COMBINATION FROM WEIGHTED AVERAGE OF: /* cod emr esa gfz jpl ngs usn whu /* REFERENCED TO IGS TIME (IGST) AND TO WEIGHTED MEAN POLE: /* PCV:IGS08_1740 OL/AL:FES2004 NONE Y ORB:CMB CLK:CMB * 2013 0 0.00000000 PG01 13548.719155 -11617.470856 19626.058441 32.681916 8 106 PG02 -13504.382657 -14580.256507 -17642.894318 435.020216 113 PG03 25127.721219 8860.976364 835.248321 214.753350 139 PG04 -8374.991090 -22794.726816 -10709.954363 234.240326 11 129 PG05 -19899.059470 -2614.351601 -17491.441018 -396.079423 128 PG06 22840.642475 13513.951718 -3369.466505 206.280890 107 PG07 6327.750771 -23891.322308 -9265.648814 220.425320 7 11 113 PG08 1428.560873 -26149.348880 1578.071995 6.506430 6 88 PG09 -7499.972290 -22879.519972 10544.687587 244.956907 98 PG10 -5157.939906 -15042.849489 -21530.162109 -81.985548 123 PG11 16052.900516 -4049.517326 20336.350450 -387.593420 114 PG12 -23998.667584 10328.199915 4053.493403 137.033904 11 11 96 PG13 9129.362910 -11700.412142 -22171.141619 113.562713 94 PG14 10312.321950 19218.296514 15462.558463 220.635036 7 90 PG15 -22872.813743 -583.240860 13673.991635 -121.561664 92 #cP2013 0 0.00000000 96 ORBIT IGb08 HLM IGS ## 1743 86400.00000000 900.00000000 56446 0.0000000000000 + 32 G01G02G03G04G05G06G07G08G09G10G11G12G13G14G15G16G17 + G18G19G20G21G22G23G24G25G26G27G28G29G30G31G32 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 ++ 2 2 2 2 2 2 2 2 ++ 2 2 2 2 2 2 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 %c G cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %f 1.2500000 1.025000000 0.00000000000 0.000000000000000 %f 0.0000000 0.000000000 0.00000000000 0.000000000000000 %i 0 0 0 0 %i 0 0 0 0 /* RAPID ORBIT COMBINATION FROM WEIGHTED AVERAGE OF: /* cod emr esa gfz jpl ngs sio usn whu /* REFERENCED TO IGS TIME (IGST) AND TO WEIGHTED MEAN POLE: /* PCV:IGS08_1740 OL/AL:FES2004 NONE Y ORB:CMB CLK:CMB * 2013 0 0.00000000 PG01 13594.216095 -11004.050777 19947.188864 33.056546 7 94 PG02 -13486.124444 -15133.972689 -17202.598344 435.178457 81 PG03 25126.143601 8948.998504 74.999599 215.172295 10 10 146 PG04 -8317.124261 -23102.681245 -10051.614102 235.136419 11 108 PG05 -19448.518625 -2946.011060 -17936.974922 -396.153400 107 PG06 22701.665946 13532.898795 -4121.909586 204.836591 124 PG07 6386.100145 -23590.241006 -9973.778438 220.789144 11 112 PG08 1477.688789 -26178.684110 780.575489 6.521192 7 122 PG09 -7377.508913 -23301.617987 9714.710845 245.196479 90 PG10 -4579.836756 -15365.902274 -21420.985701 -82.102857 90 PG11 16260.751966 -3416.102499 20294.892960 -387.901695 7 91 PG12 -23905.555411 10199.595949 4848.945271 137.213781 11 10 123 PG13 9680.417629 -11325.603384 -22134.788353 113.250392 8 102 PG14 9900.212409 18993.108189 16004.408704 220.611965 115 SV: Trần Đại An 52 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình #cP2013 0 0.00000000 96 ORBIT IGb08 HLM IGS ## 1743 172800.00000000 900.00000000 56447 0.0000000000000 + 32 G01G02G03G04G05G06G07G08G09G10G11G12G13G14G15G16G17 + G18G19G20G21G22G23G24G25G26G27G28G29G30G31G32 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 ++ 2 2 2 2 2 2 2 2 ++ 2 2 2 2 2 2 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 ++ 0 0 0 0 0 0 0 0 %c G cc GPS ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %c cc cc ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc %f 1.2500000 1.025000000 0.00000000000 0.000000000000000 %f 0.0000000 0.000000000 0.00000000000 0.000000000000000 %i 0 0 0 0 %i 0 0 0 0 /* RAPID ORBIT COMBINATION FROM WEIGHTED AVERAGE OF: /* cod emr esa gfz jpl ngs sio usn whu /* REFERENCED TO IGS TIME (IGST) AND TO WEIGHTED MEAN POLE: /* PCV:IGS08_1740 OL/AL:FES2004 NONE Y ORB:CMB CLK:CMB * 2013 0 0.00000000 PG01 13649.205239 -10381.266458 20243.079678 33.431964 PG02 -13474.392444 -15674.085087 -16739.773577 435.336425 PG03 25103.933349 9029.156059 -685.332778 215.594198 PG04 -8262.176433 -23392.051997 -9380.357096 236.027035 10 PG05 -18991.598623 -3290.841466 -18358.844853 -396.226743 PG06 22541.868013 13544.631561 -4869.326691 203.384844 PG07 6449.542225 -23268.196856 -10669.353017 221.155543 PG08 1524.885181 -26184.155624 -18.355450 6.530190 10 PG09 -7258.288615 -23692.813907 8866.437427 245.412921 PG10 -4009.687604 -15694.838324 -21284.569761 -82.219990 PG11 16475.078673 -2786.082826 20226.302889 -388.210214 PG12 -23792.632047 10056.796818 5637.520318 137.393829 PG13 10235.747725 -10960.165940 -22071.446949 112.937582 PG14 9471.417729 18765.819864 16526.800861 220.588955 PG15 -23537.264698 -1032.334178 12482.190534 -121.835291 10 8 9 9 7 7 8 6 106 131 132 110 121 101 121 110 112 122 112 142 89 104 66 3.3.3 Các số liệu thư mục GEN Phần lớn có phần mềm, bổ sung theo thời điểm đo: sat_yyyy.crx; C04_yyyy.erp Một phần số liệu trình bày bảng 3.5 SV: Trần Đại An 53 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Bảng 3.5 Một phần số liệu thư mục GEN SATELLITE PROBLEMS: MANEUVERS OR BAD OBSERVATION INTERVALS -SATELLITE *** PROBLEM * 18 27 27 30 30 30 30 30 30 30 37 102 108 114 120 125 126 126 126 322 -27 -108 26 -27 26 -27 -27 -27 -27 -27 -27 -27 -27 -27 -27 -108 -27 -108 108 108 -108 -27 -27 -27 -27 -5 -5 -27 -27 -27 -27 -27 -27 120 -20 20 20 -20 -27 -27 -27 123 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 ACTION * 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 FROM YYYY MM DD HH MM SS TO YYYY MM DD HH MM SS 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2011 2013 2013 2013 2099 2013 2012 2013 2113 2113 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 11 03 05 03 05 08 09 12 12 01 27 30 31 12 22 18 04 18 23 23 23 23 23 23 23 23 23 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 2099 2013 2013 2013 2099 2099 2013 2013 2199 2199 12 07 01 03 08 12 01 02 12 12 31 01 29 13 02 31 05 22 31 31 23 09 23 23 23 23 23 23 23 23 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 59 11 03 05 03 05 08 09 10 12 12 07 07 01 03 08 07 03 01 02 02 12 12 12 12 12 12 12 12 12 12 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 01 14 09 29 07 08 17 08 18 18 13 01 26 08 02 30 15 06 26 26 25 25 22 26 23 27 28 29 30 31 01 01 02 03 04 05 05 06 06 06 06 07 07 08 09 10 10 10 11 11 12 13 14 15 16 16 16 16 16 17 17 14 18 15 19 19 16 00 00 23 00 00 00 00 00 00 12 00 00 00 00 06 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 15 00 00 00 00 00 00 00 00 00 22 22 00 00 00 00 00 00 00 00 00 00 59 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 03 00 00 00 00 00 00 00 00 00 48 48 00 00 00 00 00 00 00 00 IERS C04 POLE SV: Trần Đại An 00 00 59 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 01 00 00 00 00 00 00 53 00 00 00 00 00 00 00 00 00 08 08 00 00 00 00 00 00 00 00 ?? ?? ?? ?? 2013 01 06 00 00 01 2013 01 16 23 59 59 11-OCT-13 00:55 54 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình NUTATION MODEL : NO SUBDAILY POLE MODEL: RAY DATE TIME X-POLE Y-POLE UT1-UTC GPS-UTC RMS XP RMS YP RMS DT DE-CPO DP-CPO RMS EP RMS PS YYYY MM DD HH MM (") (") (S) (S) REM (") (") (S) (") (") (") (") 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 2013 0.00000 0 0.07531 0.00004 0.00006 0 0.07330 0.00004 0.00007 0 0.07143 0.00005 0.00008 0 0.07010 0.00003 0.00004 0 0.06904 0.00003 0.00004 0 0.06820 0.00004 0.00005 0 0.06720 0.00005 0.00007 0 0.06629 0.00006 0.00008 0 0.06535 0.00007 0.00011 10 0 0.06441 0.00009 0.00014 11 0 0.06323 0.00011 0.00017 0.28986 0.277091 16 C04 0.00002 0.00003 0.000005 0.00000 0.28990 0.276137 16 C04 0.00002 0.00003 0.000003 0.00000 0.28970 0.275037 16 C04 0.00002 0.00003 0.000001 0.00000 0.28958 0.273858 16 C04 0.00002 0.00003 0.000000 0.00000 0.28957 0.272657 16 C04 0.00002 0.00003 0.000012 0.00000 0.28967 0.271455 16 C04 0.00002 0.00003 0.000010 0.00000 0.29007 0.270379 16 C04 0.00002 0.00003 0.000006 0.00000 0.29053 0.269461 16 C04 0.00002 0.00003 0.000001 0.00000 0.29131 0.268601 16 C04 0.00002 0.00003 0.000009 0.00000 0.29218 0.267717 16 C04 0.00002 0.00003 0.000005 0.00000 0.29300 0.266780 16 C04 0.00002 0.00003 0.000002 0.00000 3.3.4 Mô hình điện ly Mô hình điện ly gồm file COD17426.ION, COD17430.ION, COD17431.ION, COD 17432.ION Một phần số liệu trình bày bảng 3.6 Bảng 3.6 Mô hình điện ly CODE'S GLOBAL IONOSPHERE MAPS FOR DAY 152, 2013 05-JUN-13 20:48 -MODEL NUMBER / STATION NAME : 1520-01 MODEL TYPE (1=LOCAL,2=GLOBAL,3=STATION) : MAXIMUM DEGREE OF SPHERICAL HARMONICS : 15 MAXIMUM ORDER : 15 DEVELOPMENT WITH RESPECT TO GEOGRAPHICAL (=1) OR GEOMAGNETIC (=2) FRAME : MEAN (=1) OR TRUE (=2) POSITION OF THE SUN : MAPPING FUNCTION (0=NONE,1=COSZ,2=MSLM,3=ESM) : HEIGHT OF SINGLE LAYER AND ITS RMS ERROR (KM) : 450.00 0.00 COORDINATES OF EARTH-CENTERED DIPOLE AXIS LATITUDE OF NORTH GEOMAGNETIC POLE (DEGREES) : 80.22 EAST LONGITUDE (DEGREES) : -72.50 PERIOD OF VALIDITY FROM EPOCH / REFERENCE EPOCH (Y,M,D,H,M,S) : 2013 06 01 00 00 00 TO EPOCH : LATITUDE BAND COVERED MINIMUM LATITUDE (DEGREES) : -89.73 MAXIMUM LATITUDE (DEGREES) : 89.69 ADDITIONAL INFORMATION NUMBER OF CONTRIBUTING STATIONS : 289 NUMBER OF CONTRIBUTING SATELLITES : 56 ELEVATION CUT-OFF ANGLE (DEGREES) : 10 MAXIMUM TEC AND ITS RMS ERROR (TECU) : 66.69 0.42 COMMENT / WARNING : COEFFICIENTS DEGREE ORDER VALUE (TECU) RMS (TECU) 0 22.66247690 0.0312 7.66702120 0.0302 1 8.67648540 0.0325 SV: Trần Đại An 55 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp 2 2 3 3 3 4 4 4 4 5 5 5 5 5 6 6 6 6 6 6 7 7 -1 -1 -2 -1 -2 -3 -1 -2 -3 -4 -1 -2 -3 -4 -5 -1 -2 -3 -4 -5 -6 -1 SV: Trần Đại An Bộ môn: trắc địa Cao cấp – Công trình 6.50632401 -4.96019554 -0.56057276 3.24986854 0.80392265 1.37821099 -3.66069475 -2.32887629 -0.83315728 -0.76009924 0.14166464 -0.26953405 -0.69593028 0.33542654 -0.56144167 -1.10261111 0.21496871 -0.94275589 0.22293862 0.54773993 0.16197710 -0.08743816 0.57311987 1.21228502 -0.38227440 0.21851714 -1.16533448 0.29949337 0.54973834 0.28684273 0.40495678 0.56423709 0.08995775 -0.25269054 0.20860471 0.12001072 0.47895950 0.36991313 0.43839780 0.24411103 -0.27659339 0.46511894 0.03678153 -0.16619137 0.17986592 -0.07572805 1.03583893 -0.07925163 -0.06334417 -0.75179400 0.0289 0.0291 0.0313 0.0290 0.0293 0.0301 0.0263 0.0319 0.0290 0.0288 0.0294 0.0295 0.0286 0.0285 0.0280 0.0261 0.0311 0.0306 0.0275 0.0283 0.0283 0.0293 0.0265 0.0279 0.0268 0.0281 0.0275 0.0302 0.0300 0.0268 0.0270 0.0281 0.0290 0.0269 0.0268 0.0254 0.0277 0.0280 0.0281 0.0282 0.0297 0.0296 0.0259 0.0265 0.0285 0.0285 0.0260 0.0271 0.0248 0.0260 56 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình CODE'S GLOBAL IONOSPHERE MAPS FOR DAY 153, 2013 06-JUN-13 21:18 -MODEL NUMBER / STATION NAME : 1530-01 MODEL TYPE (1=LOCAL,2=GLOBAL,3=STATION) : MAXIMUM DEGREE OF SPHERICAL HARMONICS : 15 MAXIMUM ORDER : 15 DEVELOPMENT WITH RESPECT TO GEOGRAPHICAL (=1) OR GEOMAGNETIC (=2) FRAME : MEAN (=1) OR TRUE (=2) POSITION OF THE SUN : MAPPING FUNCTION (0=NONE,1=COSZ,2=MSLM,3=ESM) : HEIGHT OF SINGLE LAYER AND ITS RMS ERROR (KM) : 450.00 0.00 COORDINATES OF EARTH-CENTERED DIPOLE AXIS LATITUDE OF NORTH GEOMAGNETIC POLE (DEGREES) : 80.22 EAST LONGITUDE (DEGREES) : -72.50 PERIOD OF VALIDITY FROM EPOCH / REFERENCE EPOCH (Y,M,D,H,M,S) : 2013 06 02 00 00 00 TO EPOCH : LATITUDE BAND COVERED MINIMUM LATITUDE (DEGREES) : -89.63 MAXIMUM LATITUDE (DEGREES) : 89.69 ADDITIONAL INFORMATION NUMBER OF CONTRIBUTING STATIONS : 286 NUMBER OF CONTRIBUTING SATELLITES : 56 ELEVATION CUT-OFF ANGLE (DEGREES) : 10 MAXIMUM TEC AND ITS RMS ERROR (TECU) : 56.95 0.49 COMMENT / WARNING : COEFFICIENTS DEGREE ORDER VALUE (TECU) RMS (TECU) 0 18.79324656 0.0313 3.89743585 0.0303 1 7.24918970 0.0326 -1 3.66127553 0.0291 -6.43779497 0.0292 -1.51402139 0.0315 -1 0.66827104 0.0291 2 0.44442208 0.0294 -2 -0.38938783 0.0302 -1.38008635 0.0264 -2.08496249 0.0321 -1 -1.31895456 0.0290 -1.48924176 0.0289 -2 -0.37998128 0.0296 3 -0.36883627 0.0297 -3 -0.42776731 0.0287 1.86658805 0.0286 0.28190454 0.0282 -1 -0.58818214 0.0262 -0.22210715 0.0311 -2 0.03736080 0.0308 0.13428035 0.0277 -3 0.54823946 0.0285 CODE'S GLOBAL IONOSPHERE MAPS FOR DAY 154, 2013 07-JUN-13 20:49 -MODEL NUMBER / STATION NAME : 1540-01 MODEL TYPE (1=LOCAL,2=GLOBAL,3=STATION) : MAXIMUM DEGREE OF SPHERICAL HARMONICS : 15 MAXIMUM ORDER : 15 DEVELOPMENT WITH RESPECT TO GEOGRAPHICAL (=1) OR GEOMAGNETIC (=2) FRAME : MEAN (=1) OR TRUE (=2) POSITION OF THE SUN : MAPPING FUNCTION (0=NONE,1=COSZ,2=MSLM,3=ESM) : HEIGHT OF SINGLE LAYER AND ITS RMS ERROR (KM) : 450.00 0.00 COORDINATES OF EARTH-CENTERED DIPOLE AXIS LATITUDE OF NORTH GEOMAGNETIC POLE (DEGREES) : 80.22 EAST LONGITUDE (DEGREES) : -72.50 PERIOD OF VALIDITY FROM EPOCH / REFERENCE EPOCH (Y,M,D,H,M,S) : 2013 06 03 00 00 00 SV: Trần Đại An 57 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình TO EPOCH LATITUDE BAND COVERED MINIMUM LATITUDE (DEGREES) MAXIMUM LATITUDE (DEGREES) ADDITIONAL INFORMATION NUMBER OF CONTRIBUTING STATIONS NUMBER OF CONTRIBUTING SATELLITES ELEVATION CUT-OFF ANGLE (DEGREES) MAXIMUM TEC AND ITS RMS ERROR (TECU) COMMENT / WARNING COEFFICIENTS DEGREE ORDER VALUE (TECU) RMS (TECU) 0 19.46207812 0.0306 5.24349622 0.0295 1 7.95420551 0.0320 -1 3.43997351 0.0283 -6.28767162 0.0285 -1.18139055 0.0307 -1 0.96225069 0.0283 2 1.94199895 0.0288 -2 1.36529104 0.0295 -1.59407238 0.0258 -2.55398783 0.0313 -1 -0.83584247 0.0281 3 3 4 4 4 -2 -3 -1 -2 -3 SV: Trần Đại An 0.16695399 1.07126481 -0.56907327 0.50683527 1.93883129 0.32013896 -0.18118316 -0.18695722 0.07558568 -0.29710879 1.10725840 : : : -89.53 89.68 : : : : : 286 56 10 64.14 0.44 0.0283 0.0288 0.0291 0.0282 0.0279 0.0275 0.0255 0.0303 0.0299 0.0270 0.0279 58 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình CODE'S GLOBAL IONOSPHERE MAPS FOR DAY 155, 2013 08-JUN-13 20:54 -MODEL NUMBER / STATION NAME : 1550-01 MODEL TYPE (1=LOCAL,2=GLOBAL,3=STATION) : MAXIMUM DEGREE OF SPHERICAL HARMONICS : 15 MAXIMUM ORDER : 15 DEVELOPMENT WITH RESPECT TO GEOGRAPHICAL (=1) OR GEOMAGNETIC (=2) FRAME : MEAN (=1) OR TRUE (=2) POSITION OF THE SUN : MAPPING FUNCTION (0=NONE,1=COSZ,2=MSLM,3=ESM) : HEIGHT OF SINGLE LAYER AND ITS RMS ERROR (KM) : 450.00 0.00 COORDINATES OF EARTH-CENTERED DIPOLE AXIS LATITUDE OF NORTH GEOMAGNETIC POLE (DEGREES) : 80.22 EAST LONGITUDE (DEGREES) : -72.50 PERIOD OF VALIDITY FROM EPOCH / REFERENCE EPOCH (Y,M,D,H,M,S) : 2013 06 04 00 00 00 TO EPOCH : LATITUDE BAND COVERED MINIMUM LATITUDE (DEGREES) : -89.41 MAXIMUM LATITUDE (DEGREES) : 89.68 ADDITIONAL INFORMATION NUMBER OF CONTRIBUTING STATIONS : 287 NUMBER OF CONTRIBUTING SATELLITES : 56 ELEVATION CUT-OFF ANGLE (DEGREES) : 10 MAXIMUM TEC AND ITS RMS ERROR (TECU) : 70.27 0.45 COMMENT / WARNING : COEFFICIENTS DEGREE ORDER VALUE (TECU) RMS (TECU) 0 20.54146150 0.0303 6.09832781 0.0293 1 8.56705469 0.0315 -1 5.09185885 0.0280 -6.75350057 0.0283 -0.54468781 0.0305 -1 1.71246408 0.0282 2 1.52745580 0.0285 -2 1.67290024 0.0291 -2.59239443 0.0256 -2.60305933 0.0310 -1 -1.44470189 0.0280 -0.38610538 0.0281 -2 0.40369763 0.0286 3 -0.75139660 0.0287 -3 0.24071574 0.0277 1.38980847 0.0277 -0.37943693 0.0272 3.3.5.Tải trọng đại dương Tải trọng đại dương khai thác file FES2004.BLQ file gồm nhiều số liệu cho điểm IGS Từ số liệu điểm cần xử lý dự án (BAKO, WUHN,…) lựa chọn trình bày bảng 3.7 SV: Trần Đại An 59 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Bảng 3.7 Một phần số liệu tải trọng đại dương $$ Ocean loading displacement $$ $$ Calculated on holt using olfg/olmpp of H.-G Scherneck $$ $$ COLUMN ORDER: M2 S2 N2 K2 K1 O1 P1 Q1 MF MM SSA $$ $$ ROW ORDER: $$ AMPLITUDES (m) $$ RADIAL $$ TANGENTL EW $$ TANGENTL NS $$ PHASES (degrees) $$ RADIAL $$ TANGENTL EW $$ TANGENTL NS $$ $$ Displacement is defined positive in upwards, South and West direction $$ The phase lag is relative to Greenwich and lags positive The $$ Gutenberg-Bullen Greens function is used In the ocean tide model the $$ deficit of tidal water mass has been corrected by subtracting a uniform $$ layer of water with a certain phase lag globally $$ $$ Complete : No interpolation of ocean model was necessary $$ _PP : Ocean model has been interpolated near the station $$ (PP = Post-Processing) $$ $$ Ocean tide model: FES2004 $$ $$ Loading type: displacement $$ $$ CMC: NO (corr.tide centre of mass) $$ $$ CMC file : http://holt.oso.chalmers.se/loading/CMC/FES2004.cmc $$ CMC start : center of mass coefficient file content $$ CMC format: (a,1p,t42,3(2x,2e12.4)) $$ CMC frequ : M2 NCDF_FES2004 -1.2661E-03 -1.4298E-03 -1.3724E-03 8.2077E-04 1.1479E-03 2.3005E-04 $$ CMC frequ : S2 NCDF_FES2004 -1.7763E-04 -5.7273E-04 -5.3350E-04 -3.1591E-04 -5.1370E-05 2.8184E-04 $$ CMC frequ : N2 NCDF_FES2004 -3.2372E-04 -2.8986E-04 -2.7121E-04 1.9849E-04 2.6018E-04 -1.4302E-04 $$ CMC frequ : K2 NCDF_FES2004 -1.1814E-04 -1.5250E-04 -1.1223E-04 -1.0889E-05 -1.5751E-05 1.2367E-04 $$ CMC frequ : K1 NCDF_FES2004 -1.1370E-03 4.4839E-03 -1.8539E-03 -8.6426E-04 -9.1022E-04 -1.7823E-03 $$ CMC frequ : O1 NCDF_FES2004 -1.6802E-04 2.9702E-03 -1.3985E-03 -2.2975E-04 -8.8858E-04 -6.4989E-04 $$ CMC frequ : P1 NCDF_FES2004 -3.6495E-04 1.4941E-03 -6.1436E-04 -2.9129E-04 -2.9261E-04 -5.7461E-04 $$ CMC frequ : Q1 NCDF_FES2004 3.0709E-05 4.5472E-04 -2.7831E-04 -2.9313E-05 -2.1734E-04 -4.1637E-05 $$ CMC frequ : Mf NCDF_FES2004 -5.0643E-04 -7.3040E-05 -2.2065E-04 4.1472E-04 -1.0212E-04 8.2276E-05 $$ CMC frequ : Mm NCDF_FES2004 -2.7885E-04 2.0596E-05 4.6882E-05 1.8399E-04 -7.4897E-06 1.3209E-05 $$ CMC frequ : Ssa NCDF_FES2004 -1.4899E-04 2.6146E-06 1.3687E-04 3.5475E-05 -2.4093E-05 3.1666E-07 $$ CMC end : $$ BAKO $$ FES2004_PP ID: 2013-08-02 03:13:54 $$ Computed by OLMPP by H G Scherneck, Onsala Space Observatory, 2013 $$ BAKO, RADI TANG lon/lat: 106.8489 -6.4911 158.141 01023 00288 00246 00083 00741 00558 00254 00122 00112 00060 00050 SV: Trần Đại An 60 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình 00325 00134 00061 00035 00055 00033 00019 00007 00009 00003 00202 00068 00055 00014 00214 00115 00067 00026 00004 00003 160.4 -147.7 137.2 -150.8 -12.3 -25.5 -11.4 -42.9 -169.8 -173.5 -102.8 -64.5 -116.9 -66.6 38.1 -24.2 35.0 -15.6 -31.1 -52.7 -1.6 77.0 -37.1 75.1 158.8 151.7 155.9 142.9 -155.1 -175.0 00000 00002 -178.6 -120.1 -178.2 $$ PIMO $$ FES2004_PP ID: 2013-08-02 08:16:01 $$ Computed by OLMPP by H G Scherneck, Onsala Space Observatory, 2013 $$ PIMO, RADI TANG lon/lat: 121.0777 14.6357 00977 00350 00211 00078 00746 00834 00256 00171 00074 00047 00466 00185 00091 00048 00355 00238 00116 00042 00006 00002 00144 00060 00024 00017 00174 00146 00061 00031 00003 00001 133.6 145.4 114.7 141.7 -30.6 -54.6 -29.1 -66.2 -164.9 -169.3 102.1 119.1 92.7 111.7 -136.5 -168.8 -136.7 -176.2 -56.4 -101.3 92.1 128.7 90.4 129.2 -164.4 168.5 -166.3 155.9 1.8 17.4 $ TCMS $$ FES2004_PP ID: 2013-08-02 09:57:23 $$ Computed by OLMPP by H G Scherneck, Onsala Space Observatory, 2013 $$ TCMS, RADI TANG lon/lat: 120.9874 24.7980 01217 00288 00279 00066 00863 00801 00286 00168 00034 00025 00768 00262 00139 00064 00275 00194 00090 00035 00007 00003 00390 00113 00065 00019 00133 00124 00044 00023 00003 00002 -150.6 -153.3 174.2 179.8 -74.8 -90.6 -73.0 -101.0 -158.4 -164.3 104.0 129.2 91.0 126.8 -131.0 -156.2 -131.3 -165.9 -55.1 -103.7 -107.1 -80.7 -126.9 -85.3 177.0 144.1 176.2 131.6 23.2 26.8 $$ WUHN $$ Complete FES2004 $$ Computed by OLFG, H.-G Scherneck, Onsala Space $$ WUHN, RADI TANG lon/lat: 00329 00113 00084 00025 00334 00320 00109 00155 00069 00038 00017 00131 00094 00043 00093 00035 00019 00009 00109 00099 00036 145.6 141.3 134.8 135.5 -86.7 -98.5 -85.5 140.2 137.1 133.5 130.4 -104.6 -122.5 -104.2 -4.2 20.7 -23.7 20.7 153.0 124.7 154.4 $$ SV: Trần Đại An 61 95.523 00042 00001 00002 -177.7 -164.2 2.7 77.244 00024 00002 00003 -176.8 -162.8 3.9 Observatory 2013-Aug- 114.3573 30.5317 25.859 00067 00011 00009 00009 00018 00007 00002 00001 00019 00001 00001 00001 -107.7 -157.7 -163.0 -176.8 -133.9 -44.1 -86.2 -156.8 111.4 132.5 89.7 7.3 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình KẾT LUẬN Xử lý số liệu GNSS phần Bernese đáp ứng độ xác cao không gian lớn, số liệu đo yêu cầu số lượng lớn tệp số liệu hỗ trợ Ngoài kiến thức chung công nghê GNSS, đồ án sâu tìm hiểu trình bày cách ngắn gọn có hệ thống loại tệp hỗ trợ Số liệu hỗ trợ sản phẩm quốc tế, đa dạng, phong phú nên sở tài liệu phần mềm Bernese 5.0, đồ án đã: -Thống kê loại tệp số liệu hỗ trợ, khả giải thích vai trò xử lý -Đề cập tới số khái niệm liên quan tới công nghệ GNSS xử lý số liệu GNSS, làm sở nhận biết nhằm khai thác tệp số liệu -Tập hợp địa Website, nơi khai thác tệp số liệu -Áp dụng hiểu biết, đồ án tập hợp đầy đủ số liệu hỗ trọ cần thiết, hoạt động thiết thực đảm bảo cho việc xử lý số liệu đợt đo thành công Với điều trình bày trên, cho đồ án hoàn thành đầy đủ mục tiêu đặt ra, tiền đề thuận lợi quan trọng cho việc tìm hiểu, xử lý số liệu phần mềm Bernese tương lai TÀI LIỆU THAM KHẢO [1] Rolf Dach, Urs Hugentober, Peter Walser (2008): Bernese GPS Software Version 5.0 Astronomical Institute, University of Bern, 2008 [2] Werner Gurtner, Lou Estey (2009): RINEX: The Receiver Independent Exchange Format Version 3.01 www.ngs.noaa.gov/CORS/ SV: Trần Đại An 62 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình [3] Vy Quốc Hải (2004): So sánh kết xử lý số liệu GPS lưới địa động lực phần mềm GPSurvey 2.35 Bernese 4.2 Tạp trí Các khoa học Trái đất – PS 4(T26) tr.426-431 [4] Bài Giang Công Nghệ GPS: PGS.TS Đỗ Ngọc Đường PGS.TS Đặng Nam Chinh [5] Xây Dựng Lưới Trắc Địa: PGS.TS Đỗ Ngọc Đường [6] Giáo Trình Trắc Địa Cao Cấp: GS.TS Phạm Hoang Lân (chu biên) PGS.TS Đặng Nam Chinh (biên soạn) SV: Trần Đại An 63 Lớp:LĐH2TĐ2 [...]... các hệ số phụ thuộc vào kích thước (a, a) của ellipxoid SV: Trần Đại An 32 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình Chương 2: QUY TRÌNH XỬ LÝ SỐ LIỆU GPS ĐỘ CHÍNH XÁC CAO 2.1 GIỚI THIỆU CHUNG VỀ PHẦN MỀM XỬ LÝ SỐ LIỆU ĐỘ CHÍNH XÁC CAO BERNESE 2.1.1 Khái quát về sự phát triển của phần mềm Bernese Phần mềm BERNESE là một trong bốn phần mềm có độ ổn định và khả năng xử lý các... đến độ chính xác của lưới, nhưng không nên thiết kế góc kẹp quá nhỏ 1.3 CƠ SỞ LÝ THUYẾT XỬ LÝ SỐ LIỆU ĐO LƯỚI GPS 1.3.1 Nguyên lý cơ bản xử lý số liệu Nguyên tắc xử lý toán học của các trị đo GPS là áp dụng phương pháp “bình phương nhỏ nhất” và phép toán phân tích xác suất thống kê để tính toán số liệu và đánh giá chính xác kết quả đo Tuỳ theo thiết kế máy và kỹ thuật đo, phần mềm FIROMWARE tự động... Chỉ số DOP chia ra làm các loại: - PDOP chỉ số phân tán độ chính xác về vị trí (Positional DOP), - TDOP chỉ số phân tán độ chính xác về thời gian (Time DOP), - HDOP chỉ số phân tán độ chính xác về mặt phẳng (Horizontal DOP), - VDOP chỉ số phân tán độ chính xác về độ cao (Verical DOP), - GDOP chỉ số phân tán độ chính xác về hình học (Geometric DOP), Đồ hình phân bố vệ tinh được thiết kế sao cho chỉ số. .. khả năng xử lý tự động và tăng các lựa chọn để đưa ra những kết quả hợp lý nhất a Những trường hợp ứng dụng Bernese 5.0: - Xử lý nhanh một mạng lưới nhỏ đo máy GPS 1 hoặc 2 tần số Báo cáo nghiên cứu sử dụng phần mềm Bernese 5.0 - Tự động xử lý một mạng lưới lớn đo liên tục - Xử lý dữ liệu trong một mạng lưới có nhiều máy thu - Xử lý kết hợp nhiều loại máy thu, nhiều loại antenna khác nhau - Xử lý kết... năng xử lý tự động và tăng các lựa chọn để đưa ra kết quả hợp lý nhất Phần mềm Bernese đã được công nhận là phần mềm đa năng và đạt độ chính xác cao, đã có trên 400 các trường đại học, viện nghiên cứu sử dụng với nhiều mục đích khác nhau như: + Xử lí nhanh lưới qui mô nhỏ, sử dụng các máy thu 1 tần số và 2 tần số + Tự động xử lí các mạng lưới cố định SV: Trần Đại An 33 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ. .. môn: trắc địa Cao cấp – Công trình thuật khác nhau Song tất cả các phần mềm đều có 3 chức năng xử lý chủ yếu đó là: nhập số liệu từ máy thu sang máy tính, xử lý cặp điểm để tính ra hiệu tọa độ giữa chúng, bình sai mạng lưới gồm nhiều cặp điểm đo Việc xử lý có thể được tiến hành theo chương trình tự động gồm nhiều công đoạn nối tiếp nhau cho ra kết quả cuối cùng nhưng cũng có cả chế độ xử lý bằng tay trong... Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ môn: trắc địa Cao cấp – Công trình và xác định sự trượt chu kỳ trong tín hiệu GPS Một số máy thu có thể nhận biết được sự trượt chu kỳ và thêm vào số hiệu chỉnh tương ứng khi xử lý số liệu Mặt khác khi tính toán bình sai xử lý số liệu GPS có thể dùng sai phân bậc ba để nhận biết và xử lý sự trượt chu kỳ f Sự suy giảm độ chính xác (DOPs) do đồ hình các vệ tinh Hình 1.1.12- Khoảng... tức thời của vệ tinh với độ chính xác cỡ 20 ÷ 100 m Ngoài lịch vệ tinh quảng bá còn có lịch vệ tinh chính xác (Precise Ephemerit) Lịch vệ tinh này được thành lập từ kết quả hậu xử lý số liệu quan trắc ở các thời điểm trong khoảng thời gian quan trắc, có độ chính xác toạ độ vệ tinh cỡ 10 ÷ 50 m Sai số vị trí điểm của vệ tinh chịu ảnh hưởng gần như trọn vẹn đến độ chính xác toạ độ điểm định vị tuyệt đối... Sai số trung phương trọng số đơn vị: δo = V T WV n−t (1.4.11) Trong đó: n: Là trị quan sát t: Là các tham số (ẩn số) 1.3.2 Quy trình xử lý số liệu Ta chỉ xét các bước xử lý chủ yếu với trường hợp đo tĩnh Đễ xử lý dữ kiện đo GPS có các chương trình phần mềm lập sẵn bán kèm với máy thu Mỗi hãng sản xuất có phần mềm riêng với nhiều đặc điểm, tính năng, kỹ SV: Trần Đại An 27 Lớp:LĐH2TĐ2 Đồ án tốt nghiệp Bộ. .. chức năng về phần xử lý quỹ đạo chính xác của vệ tinh, mô hình tầng điện ly, Modul ADDNEQ Phiên bản mới nhất hiện nay (Bernese 5.0) tiếp tục kế thừa và phát huy các tính năng ưu việt của các phiên bản trước đó là: tính hiệu quả cao, độ chính xác cao và xử lý linh hoạt dữ liệu từ nhiều trạm quan trắc GPS/ GLONASS Ngoài ra, phiên bản 5.0 còn chú trọng vào các phần như: sử dụng các nguồn dữ liệu đã được ... nghiệp Bộ môn: trắc địa Cao cấp – Công trình Chương 3: CHUẨN BỊ BỘ SỐ LIỆU HỖ TRỢ CHO CÔNG TÁC XỬ LÝ SỐ LIỆU ĐO GPS ĐỘ CHÍNH XÁC CAO BẰNG PHẦN MỀM BERNESE Việc chuẩn bị số liệu hỗ trợ đủ số lượng,... THIỆU CHUNG VỀ PHẦN MỀM XỬ LÝ SỐ LIỆU ĐỘ CHÍNH XÁC CAO BERNESE 2.1.1 Khái quát phát triển phần mềm Bernese Phần mềm BERNESE bốn phần mềm có độ ổn định khả xử lý trị đo GPS với độ xác cao (tương đương... PHƯƠNG PHÁP CHUẨN BỊ CÁC LOẠI SỐ LIỆU HỖ TRỢ CHO CÔNG TÁC XỬ LÝ SỐ LIỆU ĐO GPS ĐỘ CHÍNH XÁC CAO BẲNG PHẦN MỀM BERNESE Với đắc tính toàn cầu, nói, tất số liệu hỗ trợ tổ chức quốc tế tập hợp, công bố;