CHƯƠNG 1: TỔNG QUAN VỀ CÔNG NGHỆ GPS 1.1. LỊCH SỬ PHÁT TRIỂN CỦA HỆ THỐNG GPS Từ những năm 1960 cùng với sự tiến bộ của kỹ thuật điện tử, chế tạo tên lửa và lý thuyết định vị vệ tinh, người ta đã xây dựng được các hệ thống định vị vệ tinh đầu tiên. Trước khi có hệ thống định vị toàn cầu Mỹ đã dùng hệ thống định vị vệ tinh khu vực (thuộc lãnh thổ Mỹ) như hệ thống STARFIX OMINITRACS, Ở châu Âu có hệ thống định vị vệ tinh Euteltracs gồm các vệ tinh địa tĩnh. Trước năm 1960, Bộ quốc phòng Mỹ và cơ quan hàng không NASA đã triển khai hệ thống đạo hàng mang tên TRANSIT.Hệ thống này đã sớm đạt được các ưu điểm của một hệ thống đạo hàng hoàn chỉnh và chính thức trở thành dịch vụ dẫn đường từ năm 1976, Hệ thống TRANSIT hoạt động trên nguyên lý Dopler. Các vệ tinh của hệ thống này phân bố trên 6 quỹ đạo có mặt phẳng quỹ đạo vuông góc với mặt phẳng xích đạo có độ cao 1100km. Ở độ cao này sự chuyển động của các vệ tinh khá nhanh và chịu ảnh hưởng khá lớn của lực hút trái đất. Các vệ tinh TRANSIT phát hai tín hiệu ở hai tần số là150MHz và 400MHz. Ở những tần số này các tín hiệu truyền từ vệ tinh dễ bị tầng điện ly làm chậm và dễ bị nhiễu. Để phục vụ cho nhiều công tác khác nhau cần phải xây dựng một hệ thống định vị toàn cầu với các yêu cầu mới như : Có thể định vị 2424 giờ trong một ngày ở mọi nơi trên Trái Đất và trên không gian, thời gian định vị nhanh và có độ chính xác cao. Đứng trước yêu cầu cấp thiết đó, Bộ quốc phòng Mỹ đã giao cho không quân kết hợp cùng hải quân nghiên cứu thiết lập một hệ thống định vị hiện đại, có khả năng khắc phục những nhược điểm về thời gian và độ chính xác của hệ thống TRANSIT. Đến năm 1973 sáu vệ tinh NAVASATS đầu tiên đã được sản xuất phóng lên quỹ đạo từ bang Kasas để phục vụ cho chương trình thử nghiệm. Sau kết quả thử nghiệm tốt đẹp, hệ thống này tiếp tục được xây dựng và hoàn thành vào năm 1986. Trước năm 1980, hệ thống GPS chỉ được dùng cho mục đích quân sự do Bộ Quốc phòng Mỹ quản lý.Từ những năm đầu thập kỷ 80, quân đội Mỹ cho phép sử dụng hệ thống GPS cho dân sự. Các vệ tinh đầu tiên của hệ thống được phóng lên quỹ đạo vào tháng 2 năm 1978.Toàn bộ hệ thống được đưa vào hoạt động hoàn chỉnh vào tháng 5 năm 1994.Chi phí cho việc thiết lập hệ thống này cỡ 12 tỷ USD.Ngày nay, GPS được ứng dụng vào nhiều lĩnh vực khác nhau trên hầu hết các quốc gia trên thế giới. 1.2. CẤU TRÚC HỆ THỐNG GPS Hệ thống định vị tàn cầu GPS bao gồm 3 bộ phận cấu thành: Đoạn không gian (Space Segment) Đoạn điều khiển (Control Segment) Đoạn sử dụng (User Segment) 1.2.1. Đoạn không gian Đoạn này gồm 24 vệ tinh, trong đó có 3 vệ tinh dự trữ, chuyển động trên 6 mặt phẳng quỹ đạo cách đều nhau và có góc nghiêng 55o so với mặt phẳng xích đạo của trái đất. Quỹ đạo của vệ tinh gần như tròn, vệ tinh bay ở độ cao xấp xỉ 20.200 km so với bề mặt Trái đất chuyển động trên quỹ đạo gần tròn với chu kỳ 718 phút. Do sự phân bố vệ tinh như vậy mà bất kỳ thời điểm nào, ở bất cứ vị trí nào trên Trái đất cũng có thể quan trắc được ít nhất 4 vệ tinh.
BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯƠNG HÀ NỘI KHOA: TRẮC ĐỊA – BẢN ĐỒ BÁO CÁO TỔNG KẾT ĐỀ TÀI NGHIÊN CỨU KHOA HỌC SINH VIÊN NĂM HỌC 2014 - 2015 TÊN ĐỀ TÀI KHẢO SÁT ĐỘ CHÍNH XÁC XÁC ĐỊNH TỌA ĐỘ TUYỆT ĐỐI THEO ĐỘ DÀI CA ĐO Thuộc nhóm ngành khoa học: Tài nguyên Môi trường HÀ NỘI, THÁNG – NĂM 2015 BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG TRƯỜNG ĐẠI HỌC TÀI NGUYÊN VÀ MÔI TRƯƠNG HÀ NỘI KHOA: TRẮC ĐỊA – BẢN ĐỒ BÁO CÁO TỔNG KẾT ĐỀ TÀI NGHIÊN CỨU KHOA HỌC SINH VIÊN NĂM HỌC 2014 – 2015 TÊN ĐỀ TÀI KHẢO SÁT ĐỘ CHÍNH XÁC XÁC ĐỊNH TỌA ĐỘ TUYỆT ĐỐI THEO ĐỘ DÀI CA ĐO Thuộc nhóm ngành khoa học: Tài Nguyên Môi Trường Sinh viên thực hiện: Dân tộc: Lớp: Trần Thị Mai Anh Nam, Nữ: Nữ Kinh ĐH2TĐ2 Năm thứ: Khoa: Trắc địa-Bản đồ Số năm đào tạo:4 Ngành học: Trắc địa – Bản đồ Người hướng dẫn:TS Vy Quốc Hải TS Phạm Thị Hoa HÀ NỘI, THÁNG – NĂM 2015 THÔNG TIN VỀ SINH VIÊN CHỊU TRÁCH NHIỆM CHÍNH THỰC HIỆN ĐỀ TÀI I SƠ LƯỢC VỀ SINH VIÊN: Ảnh 4x6 Họ tên: Lê Chí Hiếu Sinh ngày: 04 tháng 03 năm 1995 Nơi sinh: Hà Nội Lớp: ĐH5TĐ Khóa: 2015-2019 Khoa: Trắc địa-Bản đồ Địa liên hệ: Phú Châu – Ba Vì – Hà Nội Điện thoại: 0969026845 Email: lechihieu39@gmail.com II QUÁ TRÌNH HỌC TẬP (kê khai thành tích sinh viên từ năm thứ đến năm học): * Năm thứ 1: Ngành học: Kỹ thuật Trắc địa-Bản đồ Khoa: TĐ-BĐ Kết xếp loại học tập: TB Sơ lược thành tích: * Năm thứ 2: Ngành học: Kỹ thuật Trắc địa-Bản đồ Khoa: TĐ-BĐ Kết xếp loại học tập: K Sơ lược thành tích: Ngày tháng năm 2018 Xác nhận trường đại học Sinh viên chịu trách nhiệm (ký tên đóng dấu) thực đề tài THÔNG TIN KẾT QUẢ NGHIÊN CỨU CỦA ĐỀ TÀI Thông tin chung: - Tên đề tài: “TỰ ĐỘNG HĨA BÌNH SAI LƯỚI TỰ DO TRONG TRẮC ĐỊA CƠNG TRÌNH” - Nhóm sinh viên thực hiện: - Lớp: ĐH5TĐ - Năm thứ: Lê Chí Hiếu Mai Văn Hùng Phùng Ngọc Hưng Khoa: Trắc địa – Bản đồ Số năm đào tạo: - Người hướng dẫn: TS Đinh Xuân Vinh Mục tiêu đề tài: Tự động hóa bình sai lưới tự trắc địa cơng trình Tính sáng tạo: Là đề tài nghiên cứu ứng dụng tự động hóa bình sai lưới tự trắc địa cơng trình Kết nghiên cứu: Tự động hóa tìm nghiệm bình phương nhỏ phương trình sai số phương trình điều kiện phụ Đóng góp mặt kinh tế - xã hội,giáo dục đào tạo, an ninh, quốc phòng khả áp dụng đề tài: Cung cấp tài liệu chương trình giúp tự động hóa bình sai lưới tự trắc địa cơng trình 6.Cơng bố khoa học sinh viên từ kết nghiên cứu đề tài (ghi rõ tên tạp chí có) nhận xét, đánh giá sở áp dụng kết nghiên cứu (nếu có): Ngày tháng năm 2018 Sinh viên chịu trách nhiệm thực đề tài Nhận xét người hướng dẫn đóng góp khoa học sinh viên thực đề tài Qua đề tài nhóm sinh viên tiến bước từ sâu tìm hiểu lý thuyết đến ứng dụng vào thực tế, làm chủ công nghệ Đây tiền đề cho sinh viên nghiên cứu nhà trường (đề tài nghiên cứu khoa học, khóa luận tốt nghiệp,…) quan công tác tương lai Ngày tháng năm 2015 Xác nhận trường đại học Người hướng dẫn (ký tên đóng dấu) (ký, họ tên) LỜI CẢM ƠN LỜI NÓI ĐẦU CHƯƠNG 1: TỔNG QUAN VỀ CÔNG NGHỆ GPS .4 1.1 LỊCH SỬ PHÁT TRIỂN CỦA HỆ THỐNG GPS 1.2 CẤU TRÚC HỆ THỐNG GPS 1.2.1 Đoạn không gian 1.2.2 Đoạn điều khiển 1.2.3 Đoạn sử dụng 1.3 CÁC PHƯƠNG PHÁP ĐỊNH VỊ GPS 1.3.1 Định vị tuyệt đối 1.3.2 Định vị tương đối 11 1.4 CÁC PHƯƠNG PHÁP ĐO GPS .12 1.4.1 Phương pháp đo tĩnh (Static) 12 1.4.2 Phương pháp đo động (Kinematic) 12 1.4.3 Phương pháp đo giả động 13 1.5 CÁC NGUỒN SAI SỐ TRONG ĐO GPS .14 1.5.1 Sai số đồng hồ 14 1.5.2 Sai số quỹ đạo vệ tinh 14 1.5.3 Sai số tầng điện ly tầng đối lưu .14 1.5.4 Sai số nhiễu xạ tín hiệu vệ tinh .16 1.5.5 Sai số người đo 16 CHƯƠNG 2: GIỚI THIỆU PHẦN MỀM BERNESE 17 2.1.KHÁI QUÁT VỀ SỰ PHÁT TRIỂN CỦA PHẦN MỀM BERNESE 17 2.2 TỔNG QUAN VỀ PHẦN MỀM BERNESE 5.0 .17 2.3 KHẢ NĂNG ỨNG DỤNG PHẦN MỀM BERNESE 5.0 .18 CHƯƠNG 3: KẾT QUẢ ỨNG DỤNG 20 3.1 CÔNG TÁC THU THẬP SỐ LIỆU 20 3.2 QUY TRÌNH XỬ LÝ SỐ LIỆU 22 3.3.TOẠ ĐỘ THEO CÁC CA ĐO 31 3.4.TỔNG HỢP ĐỘ LỆCH TỌA ĐỘ .35 KẾT LUẬN VÀ KIẾN NGHỊ 43 TÀI LIỆU THAM KHẢO .44 PHỤ LỤC 45 Danh mục chữ viết tắt: IGS-International GPS Service for Geodynamics Tổ chức dịch vụ GPS quốc tế ITRF- International Terrestrial Reference Frame Khung quy chiếu mặt đất quốc tế GNSS-Global Navigation Satellite Systems Hệ thống vệ tinh dẫn đường toàn cầu GPS - Global Positioning System Hệ thống định vị toàn cầu RINEX - Receiver Independent Exchange Format Định dạng số liệu độc lập với máy thu Danh mục hình Hình 1.2 Chuyển động vệ tinh quỹ đạo Hình 1.2 Vị trí đặt trạm điều khiển GPS mặt đất .7 Hình 1.2 Phần thiết bị sử dụng dẫn đường GPS .8 YHình 1.3 Định vị tuyệt đối đo khoảng cách giả 10 Hình 1.3 Địnhvịtươngđối .11 YHình 1.5 Sai số tầng điện ly .15 Hình 1.5 Sai số tầng khí sai số đa đường dẫn 15 YHình 3.2 Quy trình xử lý số liệu phần mềm Bernese 5.0 .22 Hình 3.2 Thiết lập Project 23 Hình 3.2 Chuyển số liệu sang dạng Bernese 24 Hình 3.2 Hình ảnnh file RXOBV3.L00 26 Hình 3.2 Chuyển lịch vệ tinh 27 Hình 3.2 Phần đầu file ORB12364.OUT 28 Hình 3.2 Xác định toạ độ phần mềm Bernese 5.0 29 Hình 3.2 Hình ảnh phần file CODSPP.L00 31 Danh mục bảng: Bảng 3.1 Các tệp số liệu chung .2 Bảng 3.3 Tọa độ theo ca đo ngày 364 31 Bảng 3.3 Tọa độ theo ca đo ngày 365 32 Bảng 3.3 Tọa độ theo ca đo ngày 366 32 Bảng 3.3 Tọa độ theo ca đo ngày 001 33 Bảng 3.3 Tọa độ theo ca đo ngày 002 33 Bảng 3.3 Tọa độ theo ca đo ngày 003 34 Bảng 3.3 Tọa độ theo ca đo ngày 004 34 YBảng 3.4 Kết tính tốn ngày 364 35 Bảng 3.4 Kết tính toán ngày 365 36 Bảng 3.4 Kết tính tốn ngày 366 37 Bảng 3.4 Kết tính tốn ngày 001 38 Bảng 3.4 Kết tính tốn ngày 002 39 Bảng 3.4 Kết tính tốn ngày 003 40 Bảng 3.4 Kết tính tốn ngày 004 41 Bảng 3.4 Toạ độ gốc điểm 42 Bảng 3.4 SAI SỐ THỰC CỦA CÁC CA ĐO 24 GIỜ 42 Hình 3.2.8Hình ảnh phần file CODSPP.L00 Kết trình tọa độ điểm sai số liên quan xác định Như viết, ca đo, tọa độ xác định theo độ dài ca đo 2,4,6 22, 24 Kết xử lý tập hợp theo ca đo bảng 3.3.TOẠ ĐỘ THEO CÁC CA ĐO Bảng 3.3.1 Tọa độ theo ca đo ngày 364 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.713 8.709 8.705 8.705 8.702 8.703 8.702 8.704 8.704 8.705 8.704 8.703 109 109 109 109 109 109 109 109 109 109 109 109 32 0 0 0 0 0 0 Độ cao (m) 36.089 36.090 36.088 36.089 36.091 36.090 36.088 36.090 36.090 36.090 36.088 36.087 12.56 12.23 12.27 12.28 12.16 12.32 12.43 12.48 12.38 12.28 12.24 12.23 Bảng 3.3.2 Tọa độ theo ca đo ngày 365 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.708 8.702 8.701 8.703 8.701 8.699 8.699 8.701 8.701 8.702 8.701 8.701 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 Độ cao (m) 36.088 36.086 36.088 36.089 36.091 36.090 36.088 36.089 36.089 36.089 36.087 36.086 13.04 12.71 12.49 12.35 12.22 12.23 12.25 12.25 12.10 12.00 11.98 11.98 Bảng 3.3.3 Tọa độ theo ca đo ngày 366 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.717 8.699 8.697 8.699 8.698 8.699 8.700 8.702 8.702 8.702 8.702 8.702 109 109 109 109 109 109 109 109 109 109 109 109 33 0 0 0 0 0 0 Độ cao (m) 36.085 36.077 36.079 36.081 36.082 36.080 36.080 36.082 36.083 36.083 36.082 36.081 12.80 12.21 11.72 11.67 11.76 11.81 11.71 11.71 11.65 11.61 11.66 11.64 Bảng 3.3.4 Tọa độ theo ca đo ngày 001 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.717 8.707 8.704 8.707 8.703 8.702 8.702 8.704 8.705 8.705 8.704 8.703 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 Độ cao (m) 36.083 36.088 36.086 36.086 36.088 36.087 36.087 36.089 36.090 36.089 36.088 36.087 12.55 11.96 11.56 11.40 11.32 11.44 11.54 11.70 11.60 11.60 11.59 11.61 Bảng 3.3.5 Tọa độ theo ca đo ngày 002 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.708 8.705 8.703 8.704 8.702 8.700 8.701 8.703 8.703 8.703 8.702 8.701 109 109 109 109 109 109 109 109 109 109 109 109 34 0 0 0 0 0 0 Độ cao (m) 36.089 36.080 36.081 36.082 36.086 36.084 36.084 36.086 36.088 36.087 36.086 36.085 12.79 12.95 12.34 12.09 11.91 11.98 11.97 12.04 11.93 11.90 11.88 11.89 Bảng 3.3.6 Tọa độ theo ca đo ngày 003 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.718 8.711 8.707 8.708 8.704 8.703 8.702 8.703 8.703 8.703 8.702 8.702 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 Độ cao (m) 36.092 36.084 36.084 36.084 36.087 36.086 36.087 36.089 36.089 36.088 36.086 36.085 13.32 13.31 12.28 11.93 11.63 11.72 11.75 11.81 11.71 11.71 11.68 11.64 Bảng 3.3.7 Tọa độ theo ca đo ngày 004 Thời gian(giờ ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.712 8.706 8.703 8.707 8.704 8.703 8.703 8.705 8.705 8.704 8.703 8.703 109 109 109 109 109 109 109 109 109 109 109 109 35 0 0 0 0 0 0 Độ cao (m) 36.083 36.081 36.084 36.084 36.086 36.085 36.087 36.089 36.090 36088 36.087 36.086 12.91 12.38 11.66 11.24 11.18 11.32 11.37 11.44 11.38 11.41 11.42 11.46 3.4.TỔNG HỢP ĐỘ LỆCH TỌA ĐỘ Từ tọa độ thống kê bảng tính độ lệch tọa độ xác định từ ca đo thành phần với ca đo 24 theo thành phần tọa độ Cũng xin lưu ý, theo kết phần mềm xử lý, độ vĩ, độ kinh có đơn vị độ - phút – giây, với cơng thức thích hợp đề tài chuyển độ lệch sang đơn vị chiều dài (vì đơn vị có tính trực quan hơn) Các đồ thị thiết lập sở bảng tương ứng Bảng 3.4.1 Kết tính tốn ngày 364 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.713 8.709 8.705 8.705 8.702 8.703 8.702 8.704 8.704 8.705 8.704 8.703 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 36.089 36.090 36.088 36.089 36.091 36.090 36.088 36.090 36.090 36.090 36.088 36.087 Độ cao (m) 12.560 12.230 12.270 12.280 12.160 12.320 12.430 12.480 12.380 12.280 12.240 12.230 dLa (m) 0.303 0.182 0.061 0.061 -0.030 0.000 -0.030 0.030 0.030 0.061 0.030 0.000 dLo dh (m) (m) 0.062 0.33 0.031 0.00 0.093 0.04 0.062 0.05 0.124 -0.07 0.093 0.09 0.031 0.20 0.093 0.25 0.093 0.15 0.093 0.05 0.031 0.01 0.000 0.00 0.400 0.350 0.300 0.250 0.200 dLa (m) 0.150 dLo (m) 0.100 dh (m) 0.050 0.000 -0.050 -0.100 Hình 3.4.1 Đồ thị độ lệch toạ độ ngày 364 36 Bảng 3.4.2Kết tính tốn ngày 365 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 8.708 8.702 8.701 8.703 8.701 8.699 8.699 8.701 8.701 8.702 8.701 8.701 Kinh độ 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 36.088 36.086 36.088 36.089 36.091 36.090 36.088 36.089 36.089 36.089 36.087 36.086 Độ cao (m) 13.040 12.710 12.490 12.350 12.220 12.230 12.250 12.250 12.100 12.000 11.980 11.980 dLa (m) 0.212 0.030 0.000 0.061 0.000 -0.061 -0.061 0.000 0.000 0.030 0.000 0.000 dLo dh (m) (m) 0.062 1.06 0.062 0.73 0.000 0.51 0.093 0.37 0.155 0.24 0.124 0.25 0.062 0.27 0.093 0.27 0.093 0.12 0.093 0.02 0.031 0.00 0.000 0.00 1.200 1.000 0.800 0.600 dLa (m) dLo (m) 0.400 dh (m) 0.200 0.000 -0.200 Hình 3.4.2 Đồ thị độ lệch toạ độ ngày 365 37 Bảng 3.4.3Kết tính tốn ngày 366 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.717 8.699 8.697 8.699 8.698 8.699 8.700 8.702 8.702 8.702 8.702 8.702 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 36.085 36.077 36.079 36.081 36.082 36.080 36.080 36.082 36.083 36.083 36.082 36.081 Độ cao (m) 12.800 12.210 11.720 11.670 11.760 11.810 11.710 11.710 11.650 11.610 11.660 11.640 dLa (m) 0.454 -0.091 -0.151 -0.091 -0.121 -0.091 -0.061 0.000 0.000 0.000 0.000 0.000 dLo dh (m) (m) 0.124 1.16 -0.062 0.57 -0.124 0.08 0.000 0.03 0.031 0.12 -0.031 0.17 -0.031 0.07 0.031 0.07 0.062 0.01 0.062 -0.03 0.031 0.02 0.000 0.00 1.4 1.2 0.8 #REF! dLa (m) dLo (m) dh (m) 0.6 0.4 0.2 -0.2 -0.4 Hình 3.4.3 Đồ thị độ lệch toạ độ ngày 366 38 Bảng 3.4.4 Kết tính tốn ngày 001 Thời gian(giờ) Vĩ độ 11 36 10 12 14 16 18 20 22 24 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 8.717 8.707 8.704 8.707 8.703 8.702 8.702 8.704 8.705 8.705 8.704 8.703 Kinh độ 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 36.083 36.088 36.086 36.086 36.088 36.087 36.087 36.089 36.090 36.089 36.088 36.087 Độ cao (m) 12.55 11.960 11.560 11.400 11.320 11.440 11.540 11.700 11.600 11.600 11.590 11.610 dLa (m) dLo (m) 0.424 -0.124 0.121 -0.031 0.030 0.031 0.121 -0.031 0.000 0.031 -0.030 0.000 -0.030 0.000 0.030 0.062 0.061 0.093 0.061 0.062 0.030 0.031 0.000 0.000 dh (m) 0.94 0.35 -0.05 -0.21 -0.29 -0.17 -0.07 0.09 -0.01 -0.01 -0.02 0.00 1.200 1.000 0.800 0.600 dLa (m) 0.400 dLo (m) dh (m) 0.200 0.000 -0.200 -0.400 Hình 3.4.4 Đồ thị độ lệch toạ độ ngày 001 39 Bảng 3.4.5Kết tính tốn ngày 002 Thời gian(giờ) Vĩ độ 11 36 11 36 11 36 11 36 10 12 14 11 11 11 36 36 36 16 11 36 18 20 22 24 11 11 11 11 36 36 36 36 Kinh độ 8.708 109 36.089 8.705 109 36.080 8.703 109 36.081 8.704 8.702 8.700 8.701 109 109 109 109 0 0 36.082 36.086 36.084 36.084 8.703 8.703 8.703 8.702 8.701 109 109 109 109 109 0 0 36.086 36.088 36.087 36.086 36.085 Độ dLa cao (m) (m) 12.79 0.212 12.95 0.121 12.34 0.061 12.09 0.091 11.910 0.030 11.980 -0.030 11.970 0.000 12.04 0.061 11.930 0.061 11.900 0.061 11.880 0.030 11.890 0.000 dLo (m) dh (m) 0.124 0.90 -0.124 1.06 -0.155 0.45 -0.093 0.031 -0.031 -0.031 0.20 0.02 0.09 0.08 0.031 0.093 0.062 0.031 0.000 0.15 0.04 0.01 -0.01 0.00 1.200 1.000 0.800 0.600 dLa (m) 0.400 dLo (m) dh (m) 0.200 0.000 -0.200 -0.400 Hình 3.4.5 Đồ thị độ lệch toạ độ ngày 002 40 Bảng 3.4.6Kết tính tốn ngày 003 Thời gian(giờ) 10 12 14 16 18 20 22 24 Vĩ độ 11 11 11 11 11 11 11 11 11 11 11 11 36 36 36 36 36 36 36 36 36 36 36 36 Kinh độ 8.718 8.711 8.707 8.708 8.704 8.703 8.702 8.703 8.703 8.703 8.702 8.702 109 109 109 109 109 109 109 109 109 109 109 109 0 0 0 0 0 0 36.092 36.084 36.084 36.084 36.087 36.086 36.087 36.089 36.089 36.088 36.086 36.085 Độ cao (m) 13.320 13.310 12.280 11.930 11.630 11.720 11.750 11.810 11.710 11.710 11.680 11.640 dLa (m) dLo (m) 0.485 0.273 0.151 0.182 0.061 0.030 0.000 0.030 0.030 0.030 0.000 0.000 0.216 -0.031 -0.031 -0.031 0.062 0.031 0.062 0.124 0.124 0.093 0.031 0.000 dh (m) 1.68 1.67 0.64 0.29 -0.01 0.08 0.11 0.17 0.07 0.07 0.04 0.00 1.800 1.600 1.400 1.200 1.000 dLa (m) 0.800 dLo (m) 0.600 dh (m) 0.400 0.200 0.000 -0.200 Hình 3.4.6 Đồ thị độ lệch toạ độ ngày 003 41 Bảng 3.4.7Kết tính tốn ngày 004 Thời gian(giờ) Vĩ độ 11 36 11 36 11 36 11 36 10 11 36 12 11 36 14 11 36 16 11 36 18 11 36 20 11 36 22 11 36 24 11 36 Kinh độ 8.71 8.70 8.70 8.70 8.70 8.70 8.70 8.70 8.70 8.70 8.70 8.70 109 109 109 109 109 109 109 109 109 109 109 109 36.08 36.08 36.08 36.08 36.08 36.08 36.08 36.08 36.09 36.08 36.08 36.08 42 Độ cao (m) 12.91 12.38 dLa (m) dLo (m) 0.273 -0.093 1.45 0.091 -0.062 0.92 11.660 0.000 -0.155 0.20 11.240 0.121 -0.062 -0.22 11.180 0.030 0.000 -0.28 11.320 0.000 -0.031 -0.14 11.370 0.000 0.031 -0.09 11.440 0.061 0.093 -0.02 11.380 0.061 0.124 -0.08 11.410 0.030 0.062 -0.05 11.420 0.000 0.031 -0.04 11.460 0.000 0.000 0.00 dh (m) 1.60 1.40 1.20 1.0 0 80 dLa (m) dLo (m) dh (m) 60 0 40 0 20 0 0 - 20 - 40 Hình 3.4.7 Đồ thị độ lệch toạ độ ngày 004 Bên cạnh khảo sát mối liên hệ độ xác xác định tọa độ với độ dài ca đo, đề tài tính sai số thực toạ độ ca đo 24 với toạ độ gốc Giá trị đại lượng xác định theo phương pháp tương đối đo nối với điểm IGS, sai số vị trí điểm cỡ vài mm Toạ độ gốc điểm: Bảng 3.4.8 Toạ độ gốc điểm Vĩ độ 11 36 Kinh độ 8.683 109 36.08 Độ cao (m) 11.1 Bảng 3.4.9SAI SỐ THỰC CỦA CÁC CA ĐO 24 GIỜ Thời gian(giờ) Vĩ độ 364 11 36 365 11 36 366 11 36 Kinh độ 8.703 109 8.701 8.702 109 109 0 43 36.08 36.08 36.08 Độ cao (m) dLa (m) 12.23 0.606 -0.031 1.06 11.98 11.64 0.545 -0.216 0.576 -0.062 0.81 0.47 dLo (m) dh (m) 001 11 36 002 11 36 003 11 36 004 11 36 8.703 109 8.701 109 8.702 109 8.703 109 36.08 36.08 36.08 36.08 11.61 0.606 -0.031 0.44 11.89 0.545 -0.093 0.72 11.64 0.576 -0.093 0.47 11.46 RMS 0.606 -0.062 0.508 -0.084 0.29 0.60 3.5 NHẬN XÉT, ĐÁNH GIÁ Với số liệu bảng đồ thị trên, thấy rằng: -Tuy bảng đồ thị có đơi chút khác biệt song khẳng định xu thế: độ xác toạ độ điểm tăng dần theo thời gian, hay nói cách khác thời gian ca đo dài toạ độ xác định phương pháp tuyệt đối xác -Độ lệch tọa độ có dao động, theo thứ tự từ lớn đến nhỏ là: độ cao, độ kinh, độ vĩ, điều phản ảnh đặc điểm xác định vị trí điểm cơng nghệ GPS -Tọa độ điểm giai đoạn đầu có dao động cỡ m, sau (6-8 giờ) ổn định dần từ khoảng 12 giở trở gần thay đổi (độ lệch nhỏ 0,2 m) Đặc điểm liên quan tới chu kỳ vệ tinh ~ 12 (vệ tinh chuyển động vòng/ ngày) -So với toạ độ gốc (bảng 3.4.8), sai số thực toạ độ theo vĩ độ ~ 0.5m, kinh độ ~ 0.08m, độ cao ~ 0.6 m, sai số vị trí điểm xấp xỉ 0,8 m -So với kết xử lý GPSurvey 2.35 (Vy Quốc Hải nnk, 2004), đặc điểm biến thiên sai số thực có khác đơi chút song thống xu thế, bậc độ lớn 44 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận: Đề tài thực đầy đủ nội dung hoàn thành mục tiêu mà đề cương xác định, thể qua: -Thu thập số lượng số liệu đủ lớn (7 ca đo 24/24 giờ) Số liệu quan trắc theo phương pháp đo tĩnh, máy thu Trimble 4000ssi, anten có vành chống phản xạ, đảm bảo chất lượng -Thu thập đầy đủ lịch vệ tinh xác, số liệu hỗ trợ cho ngày đo, cho Project -Bước đầu tiếp cận, khai thác phần mềm Bernese 5.0 Cho dù việc xử lý Bernese 5.0 không đơn giản, song đề tài lựa chọn sử dụng modul thích hợp xác định tọa độ điểm -Số liệu xử lý theo độ dài ca đo: 2,4,6 22,24 Bên cạnh đó, sử dụng tọa độ xác (xác định theo phương pháp tương đối từ điểm IGS) đánh giá sai số thực từ ca đo 24 -Các kết xử lý tổng hợp, thống kê xu thế: ca đo dài, tọa độ điểm xác Sai số thực xấp xỉ 0,8 m Với sở phương pháp tường minh, số liệu có dung lượng đủ lớn đảm bảo chất lượng, xử lý phần mềm Bernese 5.0, hồn tồn tin kết xử lý có độ xác độ tin cậy cao Kiến nghị: Với nội dung phần 3.5.và phần kết luận cho rằng: -Nếu xử lý số liệu có độ xác cao phần mềm Bernese 5.0, tọa độ khởi đầu nên xác định phương pháp tương đối, đo nối với điểm IGS, theo [5] tọa độ khởi đầu so với giá trị thực phải nhỏ 0,3 m -Với nhiệm vụ trắc địa-bản đồ (không u cầu độ xác cao nhất) hồn tồn sử dụng phương pháp để xác định tọa độ điểm khởi đầu Nên ưu tiên lưới (khống chế, tăng dày) cóđiểm quan trắc 4-6 giờ, nhằm xác định tọa độ khởi đầu xác mà đáp ứng yếu tố kinh tế Trong trường hợp này, hồn tồn xác định tọa độ phần mềm thương mại (GPSurvey tương đương) 45 TÀI LIỆU THAM KHẢO [1] Đặng Nam Chinh, Đỗ Ngọc Đường (2007), Bài giảng công nghệ GPS [2] Đặng Nam Chinh, Đỗ Ngọc Đường (2012), Định vị vệ tinh [3] Vy Quốc Hải, Trần Đình Tơ, Dương Chí Cơng (2004): Xác định toạ độ tuyệt đối GPS WGS-84 Tạp chí khoa học Trái đất, 1/2004, tr.76-81, Hà Nội [4] Báo cáo xây dựng hệ quy chiếu hệ tọa độ Quốc gia.Tổng cục Địa chính, Hà Nội, 1998 [5] Rolf Dach, Urs Hugentober, Peter Walser (2008): Bernese GPS Software Version 5.0.Astronomical Institute, University of Bern, 2008 46 ... trình tọa độ điểm sai số liên quan xác định Như viết, ca đo, tọa độ xác định theo độ dài ca đo 2,4,6 22, 24 Kết xử lý tập hợp theo ca đo bảng 3.3.TOẠ ĐỘ THEO CÁC CA ĐO Bảng 3.3.1 Tọa độ theo ca đo. .. ca đo 24 phải trích xuất theo ca đo thành phần có độ dài 2, 4, 22, 24 Sau đó, xác định tọa độ điểm đo theo ca đo thành phần, từ tổng hợp so sánh với kết ca đo 24 để mối liên quan độ xác độ dài. .. 3.3 Tọa độ theo ca đo ngày 365 32 Bảng 3.3 Tọa độ theo ca đo ngày 366 32 Bảng 3.3 Tọa độ theo ca đo ngày 001 33 Bảng 3.3 Tọa độ theo ca đo ngày 002 33 Bảng 3.3 Tọa