MỞ ĐẦU 1. Tính cấp thiết của đề tài Đất nước ta đang trong quá trình phát triển mạnh mẽ, sự ra đời của các khu đô thị, các công trình dân dụng và công nghiệp… là những minh chứng cho điều đó. Để xây dựng lên các công trình thành phố, thủy lợi – thủy điện, đường hầm, công trình dân dụng và công nghiệp ấy, đặc biệt là những công trình có quy mô lớn, thì vai trò của công tác trắc địa là không thể thiếu. Một trong những nhiệm vụ của ông tác trắc địa khi xây dựng các công trình này đó là việc xây dựng hệ tọa độ cho công trình, rồi vấn đề đo nối tọa độ của công trình với tọa độ Nhà nước. Khi thực hiện những nhiệm vụ trên thì công tác tính chuyển tọa độ là hết sức quan trọng. Để hiểu rõ thêm vai trò của công tác tính chuyển tọa độ trong xây dựng các công trình và để làm luận văn tốt nghiệp của mình, em đã lựa chọn đề tài: “Nghiên cứu phương pháp chuyển đổi tọa độ ứng dụng trong trắc địa công trình”.
-1LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi, hướng dẫn khoa học TS Phạm Quốc Khánh Các số liệu, kết Luận văn trung thực chưa cơng bố cơng trình khác Hà Nội, ngày tháng 04 năm 2017 Tác giả Luận văn Đặng Đình Thịnh -2- MỤC LỤC MỤC LỤC DANH MỤC CÁC BẢNG DANH MỤC HÌNH VẼ MỞ ĐẦU Chương - CÁC HỆ TỌA ĐỘ DÙNG TRONG TRẮC ĐỊA 1.1 GIỚI THIỆU VỀ PHÉP CHIẾU 1.1.1 Phép chiếu Gauss-Kriugơ 1.1.2 Phép chiếu UTM ( Universal transverse Mercator) 12 1.2 CÁC HỆ TỌA ĐỘ NÓI CHUNG 13 1.2.1 Khái niệm hệ tọa độ 14 1.2.2 Hệ tọa độ trắc địa (B,L,H) 15 1.2.3 Hệ tọa độ vng góc khơng gian địa tâm 16 1.2.4 Hệ tọa độ vng góc khơng gian địa diện 17 1.2.5 Hệ tọa độ vng góc phẳng 19 1.3 CÁC HỆ THỐNG TỌA ĐỘ DÙNG TRÊN KHU VỰC XÂY DỰNG CƠNG TRÌNH 19 1.3.1 Hệ tọa độ nhà nước 19 1.3.2 Hệ tọa độ cơng trình 21 1.3.3 Hệ tọa độ quy ước 24 Chương - MỘT SỐ BÀI TOÁN CHUYỂN ĐỔI TỌA ĐỘ 26 -32.1 BÀI TỐN CHUYỂN ĐỔI TỪ HỆ TỌA ĐỘ VNG GĨC KHƠNG GIAN SANG HỆ TỌA ĐỘ VNG GĨC PHẲNG 26 2.1.1 Nguyên tắc chung 26 2.1.2 Phương pháp tham số 29 2.1.3 Phương pháp chuyển qua hệ tọa địa diện chân trời 33 2.2 BÀI TOÁN CHUYỂN ĐỔI TỌA ĐỘ GIỮA HAI MÚI CHIẾU 36 2.3 BÀI TOÁN CHUYỂN ĐỔI TỌA ĐỘ VNG GĨC PHẲNG 39 2.3.1 Tính chuyển công thức Helmert 40 2.3.2 Tính chuyển theo cơng thức Aphin bậc 43 2.3.3 Tính chuyển theo công thức Aphin bậc mở rộng 44 2.3.4 Tính chuyển thao cơng thức Aphin bậc hai 46 Chương – PHẦN MỀM XỬ LÝ ỨNG DỤNG CHUYỂN ĐỔI TỌA ĐỘ TRONG TRẮC ĐỊA CƠNG TRÌNH 48 3.1 SỰ CẦN THIẾT VÀ YÊU CẦU BÀI TOÁN CHUYỂN ĐỔI TỌA ĐỘ 48 3.2 NỘI DUNG CƠNG TÁC TÍNH CHUYỂN TỌA ĐỘ 49 3.2.1 Nếu tọa độ ban đầu sử dụng tọa độ trắc địa B,L 49 3.2.2 Nếu tọa độ ban đầu sử dụng tọa độ vng góc phẳng x,y 50 3.2.3 Nếu tọa độ ban đầu điểm đo GPS 50 3.2.2 Tính chuyển tọa độ độ cao khu vực 56 3.3 THÀNH LẬP PHẦN MỀM TÍNH CHUYỂN TỌA ĐỘ 61 3.3.1 Chương trình nguồn 61 3.3.2 Sử dụng phần mềm “ Convert coordinates” 62 KẾT LUẬN 69 -4- TÀI LIỆU THAM KHẢO 71 PHỤ LỤC 73 PHỤ LỤC 77 -5- DANH MỤC CÁC BẢNG Bảng 1: Kích thước Ellipxoid WGS-84 20 Bảng 2: Các tham số tính chuyển hai hệ tọa độ WGS 84 VN 2000 32 Bảng 3: Tọa độ vng góc khơng gian địa tâm khu vực Yên Phong 54 Bảng 4: Tọa độ GPS sau tính chuyển khu vực Yên Phong 55 Bảng 5: Bảng tọa độ phẳng đô cao lưới thủy điện Quảng Trị 59 Bảng 6: Bảng thành tính chuyển tọa độ ( mặt phẳng Ho =300 m) 60 Bảng 7: Bảng so sánh chiều dài cạnh sau tính chuyển 61 Bảng 8: Tọa độ điểm GPS lưới khống chế thi công cầu Bãi Cháy 66 Bảng 9: Tọa độ điểm GPS múi 103(sau tính chuyển ) 67 -6- DANH MỤC HÌNH VẼ Hình 1.1: Phép chiếu Gauss - Kriugơ 10 Hình 1.2: Phép chiếu UTM 12 Hình 1.3: Hệ tọa độ Trắc địa 16 Hình 1.4: Hệ tọa độ vng góc khơng gian địa tâm 16 Hình 1.5: Hệ tọa độ địa diện xích đạo 18 Hình 1.6: Hệ tọa độ địa diện chân trời 18 Hình 1.7: Chọn mặt chiếu trắc địa cơng trình 22 Hình 2.1: Quan hệ hai hệ tọa độ phẳng 41 Hình 3.1: Sơ đồ lưới GPS khu công nghiệp Yên Phong – Bắc Ninh 53 Hình 3.2: Lưới khống chế thủy cơng thủy điện Quảng trị 58 Hình 3.3: Giao diện chương trình 62 Hình 3.4:Giao diện chương trình tính chuyển tọa độ 63 Hình 3.5:Giao diện chương trình tính chuyển BLH sang XYZ ngược lại 64 Hình 3.6: Giao diện chương trình XYZ (WGS-84) sang xy (VN-2000) 65 Hình 3.7:Giao diện chương trình tính chuyển tọa độ GPS 66 Hình 3.8: Giao diện chương trình tính chuyển tọa độ độ cao mặt chiếu 67 -7- MỞ ĐẦU Tính cấp thiết đề tài Đất nước ta trình phát triển mạnh mẽ, đời khu đô thị, cơng trình dân dụng cơng nghiệp… minh chứng cho điều Để xây dựng lên cơng trình thành phố, thủy lợi – thủy điện, đường hầm, cơng trình dân dụng cơng nghiệp ấy, đặc biệt cơng trình có quy mơ lớn, vai trị cơng tác trắc địa khơng thể thiếu Một nhiệm vụ ông tác trắc địa xây dựng cơng trình việc xây dựng hệ tọa độ cho cơng trình, vấn đề đo nối tọa độ cơng trình với tọa độ Nhà nước Khi thực nhiệm vụ cơng tác tính chuyển tọa độ quan trọng Để hiểu rõ thêm vai trị cơng tác tính chuyển tọa độ xây dựng cơng trình để làm luận văn tốt nghiệp mình, em lựa chọn đề tài: “Nghiên cứu phương pháp chuyển đổi tọa độ ứng dụng trắc địa công trình” Mục đích đề tài - Nghiên cứu tốn tính chuyển tọa độ ứng dụng trắc địa cơng trình - Thành lập phần mềm chuyển đổi tọa độ ứng dụng trắc địa cơng trình viết ngơn ngữ lập trình Visual Basic 6.0 Nội dung đề tài - Nghiên cứu lý thuyết tính chuyển tọa độ dùng trắc địa, đặc biệt ứng dụng trắc địa cơng trình - Xử lý tốn tính chuyển đổi tọa độ phần mềm viết ngơn ngữ lập trình - Tính tốn thực nghiệm Phương pháp nghiên cứu Phương pháp nghiên cứu từ lý thuyết đến thực tiễn, kèm theo kết -8nghiên cứu tính tốn thực nghiệm cụ thể: - Nghiên cứu lý thuyết hệ tọa độ Trắc địa, toán chuyển đổi tọa độ độ cao lĩnh vực cơng trình - Tiến hành thực nghiệm cụ thể để chứng minh lý thuyết, khẳng định tính đắn, khả thi phương án đặt đến kết luận Ý nghĩa khoa học - Ý nghĩa khoa học: Góp phần làm hồn thiện quy trình đo đạc xử lý số liệu GPS lĩnh vực trắc địa cơng trình - Ý nghĩa thực tiễn: Các kết nghiên cứu đề tài mở rộng thêm hướng nghiên cứu công tác chuyển đổi tọa độ trắc địa cơng trình tài liệu tham khảo cho tổ chức cá nhân quan tâm đến đề tài Cấu trúc luận văn Cấu trúc luận văn bao gồm: 1- Phần mở đầu 2- Phần nội dung gồm chương: Chương 1: Các hệ tọa độ dùng Trắc địa Chương 2: Một số tốn tính chuyển tọa độ Chương 3: Phần mềm xử lý ứng dụng chuyển đổi tọa độ trắc địa cơng trình 3- Phần kết luận Lời cảm ơn Em xin bày tỏ lòng biết ơn đến thầy giáo hướng dẫn khoa học TS.Phạm Quốc Khánh tận tình giúp đỡ trình thực luận văn Xin cảm ơn thầy, cô khoa Trắc địa, đặc biệt thầy cô Bộ mơn Trắc địa cơng trình bạn bè đồng nghiệp tận tình bảo, giúp đỡ để em hoàn thành luận văn -9- Chương - CÁC HỆ TỌA ĐỘ DÙNG TRONG TRẮC ĐỊA 1.1 GIỚI THIỆU VỀ PHÉP CHIẾU Như biết, mặt Ellipxoid trái đất mặt quy chiếu phù hợp với hình dạng đất Sau mạng lưới trắc địa sở chiếu lên mặt Ellipxoid việc xác định tọa độ điểm mạng lưới hồn tồn thực cách chặt chẽ xác Mặt Ellipxoid hệ tọa độ gắn với thuận lợi cho cơng tác nghiên cứu khoa học trắc địa cao cấp, cho việc giải tốn phạm vi lớn Tuy nhiên mặt cong Không thuận tiện cho công tác trắc địa cơng trình số cơng tác trắc địa khác Trước hết tọa độ trắc địa điểm mặt Ellipxoid xác định đơn vị góc, chiều dài cung mặt Ellipxoid tương ứng với đơn vị lại thay đổi theo vĩ độ Mặt khác kinh tuyến dùng để xác định góc phương vị trắc địa lại khơng song song với nhau, phải xét thêm tính khơng song song sử dụng tài liệu trắc địa Ngồi cơng thức để giải tốn trắc địa, dù khoảng cách ngắn phức tạp Như vậy, để phục vụ rộng rãi cho cơng tác trắc địa trắc địa cơng trình, người ta dùng hệ tọa độ vng góc phẳng Do chiếu mặt Ellipxoid lên mặt phẳng cần thiết Có nhiều phép chiếu mặt Ellipxoid lên mặt phẳng, nội dung đề tài em xin giới thiệu phép chiếu sử dụng phổ biến Việt Nam nhiều nước giới, phép chiếu Gauss – Kruigơ phép chiếu UTM 1.1.1 Phép chiếu Gauss-Kriugơ Đây phép chiếu đồng góc đối xứng , để hạn chế biến dạng chiều dài người ta chia mặt Ellipxoid trái đất thành 60 múi chiếu dọc theo - 10 kinh tuyến Kinh tuyến chia múi chiếu thành hai phần đối xứng với gọi kinh tuyến trục, hai kinh tuyến hai biên gọi kinh tuyến biên Hiệu kinh độ hai kinh tuyến biên múi chiếu 60 Đối với trắc địa công trình ta cịn dùng loại múi chiếu 30 1030’ Số thứ tự n múi chiếu 60 kinh tuyến trục liên hệ với theo biểu thức: Lo = 60.n- 30 (1.1) đó: L – kinh độ kinh tuyến trục múi chiếu 60 Công thức (1.1) ứng với phần mặt Ellipxoid đông bán cầu kinh độ có dấu dương, cịn tây bán cầu kinh độ âm phải trừ 3600 Từ công thức (1.1) ta thấy biết Lo xác định n ngược lại Kinh tuyến gốc Greenwich biên giới hai múi chiếu 60 Nước Việt Nam (lãnh thổ, lãnh hải hải đảo) nằm múi chiếu có số thứ tự 18,19 20 ứng với kinh độ kinh tuyến trục 1050, 1110, 1170 Kinh tuyến trục múi chiếu chuyển lên mặt phẳng thành trục x, xích đạo chuyển lên mặt phẳng thành trục y Hình chiếu kinh tuyến vĩ tuyến mặt phẳng cắt thành góc vng tính chất chiếu đồng góc (Hình 1.1) Hình 1.1: Phép chiếu Gauss - Kriugơ - 70 Cần tiếp tục nghiên cứu, thử nghiệm phươn pháp tính chuyển tọa độ độ cao cho nhiều dạng cơng trình vùng địa hình khác ( đặc biệt vùng núi) để tiếp tục hoàn thiện phương pháp tính chuyển tọa độ , độ cao cơng tắc trắc địa cơng trình Đây vấn đề tiếp tục phát triển nghiên cứu thực nghiệm hoàn chỉnh đem lại hiệu rõ rệt ứng dụng công nghệ lĩnh vực trắc địa cơng trình Cuối tơi xin chân thành cảm ơn thầy giáo Phạm Quốc Khánh toàn thể thầy mơn trắc địa cơng trình giúp đỡ tơi hồn thiện Luận văn - 71 - TÀI LIỆU THAM KHẢO Cục Đo đạc Bản đồ - Bộ Tài nguyên Môi trường (2007), Hướng dẫn sử dụng tham số tính chuyển từ Hệ tọa độ quốc tế WGS-84 sang Hệ tọa độ quốc gia VN-2000 ngược lại Đặng Nam Chinh, Đỗ Ngọc Đường (2004), Định vị vệ tinh, Nxb Khoa học kỹ thuật, Hà Nội 3.Đỗ Ngọc Đường, Đặng Nam Chinh (2003), Bài giảng công nghệ GPS, Trường đại học Mỏ - Địa Chất Phạm Hoàng Lân, Đặng Nam Chinh (1998), Giáo trình trắc địa cao cấp (học phần : Trắc địa mặt cầu), Trường đại học Mỏ - Địa Chất Phan Văn Hiến nnk (1999), Trắc địa cơng trình, Nxb Giao thơng vận tải, Hà Nội Đinh Cơng Hịa(2010), Lập trình tốn Trắc địa sở,Nxb Giao thông vận tải, Hà Nội Nguyễn Quang Phúc, Hoàng Thị Minh Hương (2009), Nghiên cứu phương pháp tính chuyển tọa độ lưới GPS hệ tọa độ thi cơng cơng trình Tuyển tập Báo cáo Hội nghị khoa học “ Đo đạc đồ Việt Nam nghiệp xây dựng bảo vệ Tổ quốc” , tr 164 -171 Nguyễn Quang Phúc (2008), “ Những lưu ý sử dụng tọa độ Nhà nước Trắc địa cơng trình”, Tập chí KHCN Xây dựng Trần Viết Tuấn (2007), “ Nghiên cứu ứng dụng cơng nghệ GPS trắc địa cơng trình Việt Nam”, Luận án tiến sĩ kĩ thuật, Thư viện trường Đại học Mỏ Địa Chất 10 Trần Bạch Giang (2004), Giới thiệu hệ quy chiếu hệ tọa độ quốc gia Việt Nam - 72 11 Công ty tư vấn xây dựng điện 1( 2003), Báo cáo xây dựng lưới tam giác thủy công lưới thủy chuẩn thủy cơng cơng trình thủy lợi – thủy điện Quảng Trị 12 TCVN 9401:2012, Công tác trắc địa xây dựng cơng trình – u cầu chung, Bộ Xây dựng, Hà Nội -2012 13 TCVN 9401:2012, Kỹ thuật đo xử lý số liệu GPS Trắc địa công trình, Bộ Xây dựng, Hà Nội -2012 14 Trung tâm Thơng tin – Lưu trữ tư liệu Địa (2001), Tài liệu hướng dẫn sử dụng chương trình Geotools 1.2 - 73 - PHỤ LỤC KẾT QUẢ TÍNH CHUYỂN TỌA ĐỘ - 74 A TINH CHUYEN TOA DO TU BL SANG XY & NGUOC LAI KET QUA TINH CHUYEN TOA DO -Tinh Chuyen BL-xy =============***============ Luoi chieu UTM_Elipxoid WGS-84 quoc te He thong goc: Mui chieu 6^, Hs bien dang chieu dai Ko=0.9996 Tinh chuyen toa Trac dia B,L toa X,Y +===========================================+ | S | TEN | TOADO | TOA DO | | T | | TRAC DIA | VUONG GOC PHANG | | T | DIEM | B-L | XY | | | | o , ,, | 105 0.0 | | -| | | -| | 1| GPS1| 21 55.8948| 2318654.895 | | | |105 16.3015| 505666.465 | | -| | | -| | 2| GPS2| 21 27.8032| 2317788.659 | | | |105 43.9556| 504733.010 | | -| | | -| | 3| GPS5| 21 37.7917| 2318099.115 | | | |104 54 1.2833| 489644.887 | | -| | | -| | 4| GPS6| 21 40.0000| 2325563.747 | | | |105 0.0000| 506924.993 | KET QUA TINH CHUYEN TOA DO -Tinh Chuyen Toa Do xy sang BL =============***============ Luoi chieu UTM_Elipxoid WGS-84 quoc te He thong goc: Mui chieu 6^, Hs bien dang chieu dai Ko=0.9996 Tinh chuyen toa X,Y toa Trac dia B,L ============================================= | S | TEN | TOADO | TOA DO | | T | | TRAC DIA | VUONG GOC PHANG | | T | DIEM | B-L | XY | | | | o , ,, | 105 0.0 | | -| | | -| | 1| GPS1| 21 27.4921| 2324838.329 | | | |105 16.3949| 505668.829 | | -| | | -| | 2| GPS2| 21 59.3888| 2323974.070 | | | |105 44.0336| 504734.984 | | -| | | -| | 3| GPS3| 21 18.5777| 2324720.768 | | | |104 18 15.4449| 427704.900 | | -| | | -| | 4| GPS4| 21 50.4776| 2323860.947 | | | |104 17 43.0880| 426767.058 | | -| | | -| | 5| GPS4| 21 0.4673| 2324241.514 | | | |104 0.2369| 411674.811 | | -| | | -| - 75 B TINH CHUYEN TOA DO QUA CAC MUI CHIEU KET QUA TINH CHUYEN TOA DO -Tinh Chuyen Toa Do Mui Chieu =============***============ Luoi chieu UTM_Elipxoid WGS-84 quoc te He thong goc: Mui chieu 6^, Hs bien dang chieu dai Ko=0.9996 He thong moi: Mui chieu 3^, Hs bien dang chieu dai Ko=0.9999 ================================================================= | S | TEN | TOADO | TOA DO VUONG GOC PHANG (XY) | | T | | TRAC DIA | -| | T | DIEM | B-L | HE THONG GOC | HE THONG MOI | | | | o , ,,| 105 0.0 | 106 0.0 | | -| | | | | | 1| GPS1 | 21 27.4921 | 2324838.329 | 2319918.771 | | | |105 16.3949 | 505668.829 | 401753.514 | | -| | | | | | 2| GPS2 | 21 59.3888 | 2323974.070 | 2319057.929 | | | |105 44.0336 | 504734.984 | 400814.164 | | -| | | | | | 3| GPS3 | 21 18.5777 | 2324720.768 | 2320288.810 | | | |104 18 15.4449 | 427704.900 | 323773.295 | | -| | | | | | 4| GPS4 | 21 50.4776 | 2323860.947 | 2319432.264 | | | |104 17 43.0880 | 426767.058 | 322829.785 | | -| | | | | | 5| GPS4 | 21 0.4673 | 2324241.514 | 2319908.324 | | | |104 0.2369 | 411674.811 | 307735.069 | | -| | | | | C.TINH CHUYEN XYZ (WGS-84) SANG XY (VN-2000) KET QUA TINH CHUYEN TOA DO GPS ============**================ TEN CONG TRINH: Cau Bai Chay 1/So luong diem can tinh: + -+ + + | TOA DO XYZ TRONG HE WGS-84 + XY TRONG HE VN-2000 + + -| -+ + | No |Ten diem | X(m) | Y(m) | Z(m) + x(m) | y(m) | + -+ -+ -+ + + + -+ | 1| GPS19 | -1748205.206| 5696493.076| 2267454.689| 2361936.277| 713906.459| | 2| BM29 | -1748255.585| 5696442.799| 2267527.835| 2362018.542| 713968.146| | 3| HL03 | -1748969.774| 5696355.963| 2267205.691| 2361676.525| 714679.195| | 4| GPS12 | -1748891.209| 5696404.603| 2267102.621| 2361569.260| 714591.435| + -+ -+ -+ + + + -+ - 76 D TINH CHUYEN TOA DO GPS VE TOA DO CONG TRINH KET QUA TINH CHUYEN TOA DO GPS ============**================ TEN CONG TRINH: Cau Bai Chay 1/So luong diem can tinh: 2/TOA DO VUONG GOC KHONG GIAN DIA TAM + -+ | No |Ten diem | X(m) | Y(m) | Z(m) | + -+ -+ -+ + + | 1| GPS19 | -1748205.206| 5696493.076| 2267454.689| | 2| BM29 | -1748255.585| 5696442.799| 2267527.835| | 3| HL03 | -1748969.774| 5696355.963| 2267205.691| | 4| GPS12 | -1748891.209| 5696404.603| 2267102.621| + -+ 3/TOA DO VUONG GOC KHONG GIAN DIA DIEN\n') + -+ | No |Ten diem | x(m) | y(m) | z(m) | + -+ -+ + -+ + | 1| GPS19 | 132.945| -385.572| 5.978| | 2| BM29 | 219.257| -316.045| 1.064| | 3| HL03 | -132.945| 385.572| 4.019| | 4| GPS12 | -231.481| 302.805| -10.962| + -+ 4/TOA DO THI CONG + + | No |Ten diem | x(m) | y(m) | + -+ -+ -+ + | 1| GPS19 | 2319133.959| 402321.422| | 2| BM29 | 2319219.873| 402391.460| | 3| HL03 | 2318863.538| 403091.087| | 4| GPS12 | 2318765.476| 403007.737| + + - 77 - PHỤ LỤC CHƯƠNG TRÌNH TÍNH ỨNG DỤNG - 78 PROGRAME TINH CHUYEN TOA DO GPS Option Explicit Dim FileSL$, FileKQ$, Ngoc%, Nst%, E2#, AA#, bb# Dim B0 As Double, L0 As Double, H0 As Double, x1() As Double, y1() As Double, Z1() As Double Dim x2() As Double, y2() As Double, Z2() As Double, X() As Double, Y() As Double, Z() As Double Dim Tenluoi$ Dim Nghiem() As Double, TS1() As Double, TS2() As Double Dim tendiem() As String, ttst() As Integer, Hs() As Double, TS() As Double, twords() As String Const Pi As Double = 3.14159265358979 Const pi2 As Double = Pi * 2# Const Ro As Double = (180# * 3600#) / Pi Const AW As Double = 6378137# Const EW As Double = 6.6943800042084E-03 Const EPW As Double = 6.73949675658691E-03 Const AK As Double = 7378245# Const EPK As Double = 0.00673852541468 Const EGK As Double = 0.00669342162296 Private Sub CmdConvert_Click() Dim txtSl$, k% Dim Xtb#, Ytb#, Ztb#, xtb1#, ytb1#, xtb2#, ytb2#, xtb3#, ytb3# Dim i%, j% Dim B1() As Double, L1() As Double, H1() As Double FileKQ = DoiTenFile(FileSL, ".KQ") Open FileSL For Input As #1 Open FileKQ For Output As #4 Input #1, Tenluoi Line Input #1, txtSl Ngoc = Val(txtSl) ReDim x1(1 To Ngoc) ReDim y1(1 To Ngoc) ReDim Z1(1 To Ngoc) ReDim x2(1 To Ngoc) ReDim y2(1 To Ngoc) ReDim Z2(1 To Ngoc) ReDim tendiem(1 To Ngoc) Xtb = 0# Ytb = 0# Ztb = 0# For i = To Ngoc Line Input #1, txtSl Words txtSl, " ", twords() k = Val(twords(1)) tendiem(i) = twords(2) x1(i) = twords(3) Xtb = Xtb + x1(i) y1(i) = twords(4) Ytb = Ytb + y1(i) Z1(i) = twords(5) Ztb = Ztb + Z1(i) Next Xtb = Xtb / Ngoc Ytb = Ytb / Ngoc Ztb = Ztb / Ngoc Line Input #1, txtSl Words txtSl, " ", twords() Nst = Val(twords(1)) - 79 ReDim ttst(1 To Nst) For i = To Nst ttst(i) = Val(twords(i + 1)) Next For i = To Nst Line Input #1, txtSl Words txtSl, " ", twords() x2(i) = Val(twords(2)) xtb1 = xtb1 + x2(i) y2(i) = Val(twords(3)) ytb1 = ytb1 + y2(i) Next xtb1 = xtb1 / Nst ytb1 = ytb1 / Nst If OptWGS84 = True Then AA = 6378137 bb = AA * (1 - / 298.257223563) E2 = (AA * AA - bb * bb) / (AA * AA) End If If OptKrasovski = True Then AA = AK bb = AA * (1 - / 298.257223563) E2 = (AA * AA - bb * bb) / (AA * AA) End If ReDim B1(1 To Ngoc) ReDim L1(1 To Ngoc) ReDim H1(1 To Ngoc) For i = To Ngoc Call ChuyenXYZBLH(x1(i), y1(i), Z1(i), B1(i), L1(i), H1(i)) Next Call ChuyenXYZBLH(Xtb, Ytb, Ztb, B0, L0, H0) xtb2 = 0# ytb2 = 0# ReDim X(1 To Ngoc) ReDim Y(1 To Ngoc) ReDim Z(1 To Ngoc) For i = To Ngoc Call ChuyenBLHXYZ(B1(i), L1(i), B0, L0, H0, x1(i), y1(i), Z1(i), X(i), Y(i), Z(i)) Next For i = To Nst xtb2 = xtb2 + X(ttst(i)) ytb2 = ytb2 + Y(ttst(i)) Next xtb2 = xtb2 / Nst ytb2 = ytb2 / Nst For i = To Nst X(ttst(i)) = (X(ttst(i)) - xtb2) Y(ttst(i)) = (Y(ttst(i)) - ytb2) Next For i = To Nst x2(i) = (x2(i) - xtb1) y2(i) = (y2(i) - ytb1) Next ' ' ''TinhtsHelmert; If OptHelMert.Value = True Then If Nst < Then - 80 MsgBox " So diem song trung phai lon hon " GoTo L: End If ReDim Hs(1 To * Ngoc, To 4) ReDim TS(1 To 4) ReDim Nghiem(1 To * Ngoc) Call Helmert(X(), Y(), x2(), y2(), ttst(), TS(), Nst) For i = To Ngoc Hs(i, 1) = Hs(i, 2) = Hs(i, 3) = X(i) Hs(i, 4) = -Y(i) 'Vx(i, j) = [1,0,x1(i),-y1(i)] Hs(i + Ngoc, 1) = Hs(i + Ngoc, 2) = Hs(i + Ngoc, 3) = Y(i) Hs(i + Ngoc, 4) = X(i) 'Vy(i+Nst,:)=[0,1,y1(tt2(i)),x1(tt2(i))]; Next Nhan_Matran1 Hs(), TS(), Nghiem(), * Ngoc, xtb3 = (x2(1) + xtb1) - Nghiem(1) ytb3 = (y2(1) + ytb1) - Nghiem(1 + Ngoc) For i = To Ngoc Nghiem(i) = Nghiem(i) + xtb1 Nghiem(i + Ngoc) = Nghiem(i + Ngoc) + ytb1 Next End If ' ' ''Tinhts APhine; If OptAphine.Value = True Then ReDim Hs(1 To Ngoc, To 4) ReDim TS1(1 To 4) ReDim TS1(1 To 4) ReDim Nghiem(1 To * Ngoc) For i = To Ngoc Hs(i, 1) = Hs(i, 2) = Hs(i, 3) = X(i) Hs(i, 4) = -Y(i) 'Vx(i, j) = [1,0,x1(i),-y1(i)] Next Call APhine(x1(), y1(), x2(), y2(), ttst(), TS1(), TS2(), Nst) End If Call Infile Open FileKQ For Input As #4 TxtDisplay.Text = "" Do While Not EOF(4) Line Input #4, txtSl TxtDisplay.Text = TxtDisplay.Text & txtSl & vbCrLf Loop Close #4 MsgBox ("Done") Close #1 Close #4 L: 'msgbox('Done'); - 81 End Sub Private Sub CmdExit_Click() FrmGPS.Hide End Sub Private Sub CmdOpenFile_Click() Dim txtSl As String FileSL = "" CdlFile.Filter = "Data(*.txt)|*.dat|So lieu(*.sl)|*.sl|All File(*.*)|*.*" CdlFile.ShowOpen FileSL = CdlFile.FileName FileKQ = DoiTenFile(FileSL, ".Kq") 'On Error GoTo L: Open FileSL For Input As #1 Do While Not EOF(1) Line Input #1, txtSl TxtDisplay.Text = TxtDisplay.Text & txtSl & vbCrLf Loop Close #1 'End Sub 'L: 'MsgBox " Ban lien he voi anh Thinh Trac dia B-k54 " End Sub Private Sub Helmert(x1() As Double, y1() As Double, x2() As Double, y2() As Double, tt2() As Integer, Thamso() As Double, Nst%) Dim V() As Double, L() As Double, LT() As Double Dim VT() As Double, PTC() As Double, QQ() As Double Dim i%, j% ReDim V(1 To * Nst, To 4) ReDim L(1 To * Nst) ReDim LT(1 To 4) ReDim VT(1 To 4, To * Nst) ReDim QQ(1 To 4, To 4) For i = To Nst V(i, 1) = V(i, 2) = V(i, 3) = x1(tt2(i)) V(i, 4) = -y1(tt2(i)) 'Vx(i, j) = [1,0,x1(tt2(i)),-y1(tt2(i))] V(i + Nst, 1) = V(i + Nst, 2) = V(i + Nst, 3) = y1(tt2(i)) V(i + Nst, 4) = x1(tt2(i)) 'Vy(i+Nst,:)=[0,1,y1(tt2(i)),x1(tt2(i))]; L(i) = -x2(i) 'L(i,:)=-x2(i); L(i + Nst) = -y2(i) 'L(i+Nst,:)=-y2(i); Next Matran_chuyen_vi V(), VT(), * Nst, Nhan_Matran VT(), V(), PTC(), 4, * Nst, Nghichdao PTC(), QQ(), Nhan_Matran1 VT(), L(), LT(), 4, * Nst For i = To LT(i) = -LT(i) Next Nhan_Matran1 QQ(), LT(), Thamso(), 4, 'Thamso=(-1)*inv((V')*V)*((V')*L); - 82 End Sub Private Sub APhine(x1() As Double, y1() As Double, x2() As Double, y2() As Double, tt2() As Integer, Thamso1() As Double, Thamso2() As Double, Nst%) Dim V() As Double, La() As Double, LT1() As Double, Lb() As Double, LT2() As Double Dim VT() As Double, PTC() As Double, QQ() As Double Dim i%, j% ReDim V(1 To Nst, To 4) ReDim La(1 To Nst) ReDim Lb(1 To Nst) ReDim LT1(1 To 4) ReDim LT2(1 To 4) ReDim VT(1 To 4, To * Nst) ReDim QQ(1 To 4, To 4) For i = To Nst V(i, 1) = V(i, 2) = x1(tt2(i)) V(i, 3) = y1(tt2(i)) V(i, 4) = x1(tt2(i)) * y1(tt2(i)) La(i) = -x2(i) 'L(i,:)=-x2(i); Lb(i) = -y2(i) 'L(i+Nst,:)=-y2(i); Next Matran_chuyen_vi V(), VT(), Nst, Nhan_Matran VT(), V(), PTC(), 4, Nst, Nghichdao PTC(), QQ(), Nhan_Matran1 VT(), La(), LT1(), 4, Nst Nhan_Matran1 VT(), Lb(), LT2(), 4, Nst For i = To LT1(i) = -LT1(i) LT2(i) = -LT2(i) Next Nhan_Matran1 QQ(), LT1(), Thamso1(), 4, Nhan_Matran1 QQ(), LT2(), Thamso2(), 4, 'Thamso=(-1)*inv((V')*V)*((V')*L); End Sub Private Sub ChuyenXYZBLH(X#, Y#, Z#, B#, L#, H#) Dim p#, Bgd#, No#, Ho#, EP#, Bs# p = Sqr(X * X + Y * Y) Bgd = (Z / p) * (1 - E2) ^ -1 No = AA * AA / (Sqr(AA * AA * Cos(Bgd) * Cos(Bgd) + bb * bb * Sin(Bgd) * Sin(Bgd))) Ho = (p / (Cos(Bgd))) - No B = Atn((Z / p) / (1 - (E2 * No / (No + Ho)))) EP = Abs(B - Bgd) Do While EP > (0.00001 * Pi / (3600# * 180#)) No = AA * AA / (Sqr(AA * AA * Cos(B) * Cos(B) + bb * bb * Sin(B) * Sin(B))) Ho = (p / (Cos(B))) - No Bs = Atn((Z / p) / (1 - (E2 * No / (No + Ho)))) EP = Abs(B - Bs) B = Bs Loop H = Ho L = Atn(Y / X) If L < Then L = L + Pi End If - 83 End Sub Private Sub ChuyenBLHXYZ(B#, L#, Bo#, Lo#, Ho#, Xx#, Yy#, ZZ#, XXX#, YYY#, ZZZ#) Dim No# Dim mtx(1 To 3) As Double, mty(1 To 3) As Double, mtz(1 To 3) As Double Dim mtb1#, mtb2#, mtb3# No = AA * AA / (Sqr(AA * AA * Cos(Bo) * Cos(Bo) + bb * bb * Sin(Bo) * Sin(Bo))) mtx(1) = -Sin(B) * Cos(L) mtx(2) = -Sin(B) * Sin(L) mtx(3) = Cos(B) mty(1) = -Sin(L) mty(2) = Cos(L) mty(3) = 0# mtz(1) = Cos(B) * Cos(L) mtz(2) = Cos(B) * Sin(L) mtz(3) = Sin(B) mtb1 = (Xx - (No + Ho) * Cos(Bo) * Cos(Lo)) mtb2 = Yy - (No + Ho) * Cos(Bo) * Sin(Lo) mtb3 = ZZ - (No * (1 - E2) + Ho) * Sin(Bo) XXX = mtx(1) * mtb1 + mtx(2) * mtb2 + mtx(3) * mtb3 YYY = mty(1) * mtb1 + mty(2) * mtb2 + mty(3) * mtb3 ZZZ = mtz(1) * mtb1 + mtz(2) * mtb2 + mtz(3) * mtb3 End Sub Private Sub Infile() Dim i As Integer Print #4, "KET QUA TINH CHUYEN TOA DO GPS" Print #4, "============**================" Print #4, "TEN CONG TRINH: " & Tenluoi Print #4, "1/So luong diem can tinh: " & Ngoc & vbCrLf & vbCrLf Print #4, "2/TOA DO VUONG GOC KHONG GIAN DIA TAM" Print #4, "+ -+" Print #4, "| No |Ten diem | X(m) | Y(m) | Z(m) |" Print #4, "+ -+ -+ -+ + +" For i = To Ngoc Print #4, "|" & Dinhdang1(i, 5) & "|" & Dinhdang1(tendiem(i), 8) & _ " |" & Fm(x1(i), 13, 3) & "|" & Fm(y1(i), 12, 3) & "|" & Fm(Z1(i), 12, 3) & "|" Next Print #4, "+ -+" & vbCrLf & vbCrLf Print #4, "3/TOA DO VUONG GOC KHONG GIAN DIA DIEN\n')" Print #4, "+ -+" Print #4, "| No |Ten diem | x(m) | y(m) | z(m) |" Print #4, "+ -+ -+ + -+ +" For i = To Ngoc Print #4, "|" & Dinhdang1(i, 5) & "|" & Dinhdang1(tendiem(i), 8) & _ " |" & Fm(X(i), 12, 3) & "|" & Fm(Y(i), 13, 3) & "|" & Fm(Z(i), 12, 3) & "|" Next Print #4, "+ -+" & vbCrLf & vbCrLf Print #4, "4/TOA DO THI CONG" Print #4, "+ +" Print #4, "| No |Ten diem | x(m) | y(m) |" Print #4, "+ -+ -+ -+ +" For i = To Ngoc Print #4, "|" & Dinhdang1(i, 5) & "|" & Dinhdang1(tendiem(i), 8) & _ - 84 " |" & Fm(Nghiem(i), 13, 3) & "|" & Fm(Nghiem(i + Ngoc), 12, 3) & "|" Next Print #4, "+ +" & vbCrLf & vbCrLf Print #4, " Ngay thuc hien:" & Date Close (4) End Sub ... dụng trắc địa cơng trình? ?? Mục đích đề tài - Nghiên cứu tốn tính chuyển tọa độ ứng dụng trắc địa cơng trình - Thành lập phần mềm chuyển đổi tọa độ ứng dụng trắc địa cơng trình viết ngơn ngữ lập trình. .. cơng trình, b, Tính chuyển tọa độ độ cao khu vực cơng trình 3.2.1 Nếu tọa độ ban đầu sử dụng tọa độ trắc địa B,L Trường hợp 1: Tọa độ trắc địa hệ quy chiếu với hệ tọa độ cơng trình Khi đó, tính chuyển. .. niệm hệ tọa độ nói 1.2.2 Hệ tọa độ trắc địa (B,L,H) Đây hệ tọa độ trắc địa Vị trí điểm T nằm trái đất xác định hệ tọa độ trắc địa theo ba tọa độ sau ( B,L,H) [2] -Vĩ độ trắc địa B: góc kẹp pháp