Giải pháp xử lý dữ liệu tích hợp phục vụ xây dựng mô hình mặt biển

Bài viết Giải pháp xử lý dữ liệu tích hợp phục vụ xây dựng mô hình mặt biển phân tích hiện trạng cơ sở toán học của các nguồn dữ liệu khác nhau (dữ liệu toàn cầu và dữ liệu quốc gia) được thu thập phục vụ xây dựng một số mô hình mặt biển.

Nghiên cứu GIẢI PHÁP XỬ LÝ DỮ LIỆU TÍCH HỢP PHỤC VỤ XÂY DỰNG MƠ HÌNH MẶT BIỂN LƯƠNG THANH THẠCH Trường Đại học Tài nguyên Môi trường Hà Nội Tóm tắt: Bài báo phân tích trạng sở toán học nguồn liệu khác (dữ liệu toàn cầu liệu quốc gia) thu thập phục vụ xây dựng số mơ hình mặt biển Trên sở phân tích ưu - nhược điểm số giải pháp, tác giả đề xuất lựa chọn hệ tọa độ quốc tế WGS84 hệ độ cao quốc gia Hòn Dấu (HP72) làm sở toán học để xây dựng mơ hình mặt biển Kết thực nghiệm cho thấy phương án đề xuất phù hợp với thực tiễn mang tính khả thi cao Đặt vấn đề Để xây dựng mơ hình mặt biển phục vụ điều tra bản, khai thác sử dụng quản lý tổng hợp tài nguyên môi trường biển cần sử dụng nhiều nguồn số liệu khác Số liệu toàn cầu: Các mơ hình mặt biển trung bình động lực MDT, mơ hình mặt biển thấp tồn cầu, số liệu đo cao vệ tinh,… Số liệu nước quan thuộc chuyên ngành khác nhau, bao gồm: Mơ hình số độ cao DEM, đồ địa hình tỷ lệ 1:25.000 (của Cục Bản đồ, Bộ Tổng tham mưu); hải đồ, độ cao mặt biển thấp trạm quan trắc mực nước (của Đoàn Đo đạc Biên vẽ hải đồ Nghiên cứu biển); độ cao mặt biển (trung bình, trung bình cao nhiều năm, trung bình thấp nhiều năm, cao thấp nhiều năm) trạm quan trắc mực nước (của Trung tâm Hải văn); đồ địa hình đáy biển (của Trung tâm Trắc địa Bản đồ biển),… Mỗi loại số liệu, tuỳ theo nhu cầu sử dụng nguồn gốc xuất xứ dựa sở toán học khác Như vậy, cần phải thống sở toán học tất số liệu hệ quy chiếu hệ tọa độ, đảm bảo phù hợp với thực tế khai thác sử dụng Bài báo khoa học trình bày sở tốn học loại liệu, đề xuất lựa chọn hệ quy chiếu hệ tọa độ phù hợp với phạm vi lãnh thổ Việt Nam thực chuyển đổi dạng liệu nhằm thống sở toán học nguồn số liệu, phục vụ nội dung nghiên cứu để xây dựng mơ hình mặt biển phục vụ điều tra bản, khai thác sử dụng quản lý tổng hợp tài nguyên môi trường biển Giải vấn đề 2.1 Cơ sở toán học loại liệu 2.1.1 Mơ hình trọng trường Trái Đất EGM2008 Mơ hình trọng trường Trái Đất EGM2008 Cơ quan Tri thức - Địa không gian quốc gia (NGA) Mỹ xây dựng Các hệ số điều hòa phép khai triển địa theo hàm cầu Ngày nhận bài: 11/8/2022, ngày chuyển phản biện: 15/8/2022, ngày chấp nhận phản biện: 19/8/2022, ngày chấp nhận đăng: 28/8/2022 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu mơ hình có bậc mức đến 2159, số hệ số bổ sung khai triển đến mức 2190 Dữ liệu sử dụng để xây dựng mơ hình EGM2008 bao gồm: liệu trọng trường dự án vệ tinh CHAMP, GRACE với mơ hình trọng trường tồn cầu GGM2S; mơ hình địa hình đại dương động lực (DOT); mơ hình trọng trường kết hợp EGM96, GGM02C, EIGEN-05C; mơ hình sử dụng vệ tinh EIGEN05S, GGM03S, ITG-GRACE03S; mơ hình mặt đất dư (RTM) mơ hình trọng trường ban đầu PGM2007B Nhóm làm việc quốc tế bảo trợ Hội Trắc địa quốc tế (IAG) Trên đại dương sử dụng giá trị dị thường trọng lực xác định kết đo cao từ vệ tinh altimetry tổ chức DNSC SIO/NOAA với mô hình DOT tương ứng [6] Các hệ số điều hịa hệ khơng phụ thuộc triều mơ hình EGM2008 xác định tương ứng với ellipsoid trung bình TFS2008 dựa giá trị dị thường trọng lực phạm vi toàn cầu Tổ chức NGA chuyển đổi giá trị dị thường độ cao toàn cầu xác định theo hệ số điều hịa mơ hình EGM2008 từ ellipsoid trung bình TFS2008 ellipsoid quy chiếu chung WGS84 phục vụ ứng dụng trắc địa tồn cầu Mơ hình độ cao quasigeoid tồn cầu EGM2008-WGS84 cung cấp hai dạng: tiện ích “1x1-Minute Geoid Undulation Grid in WGS84 - BIG ENDIAN” theo địa chỉ: http:/earthinfo.nga.mil/GandG/wgs84/gravitymod/egm 2008/und_min 1x1_egm2008_isw=82_WGS84_TideFree.gz  Mạng lưới (grid) ô chuẩn hình vng độ phân giải 2.5’ x 2.5’ bao phủ toàn cầu Các liệu đỉnh ô chuẩn bao gồm tọa độ trắc địa B, L (tương ứng với ellipsoid quy chiếu chung WGS84) dị thường độ cao tồn cầu  n hệ khơng phụ thuộc triều Dữ liệu nêu truy nhập từ Website “EGM2008-WGS84 Version” với tiện ích “2.5x2.5-Minute Geoid Undulation Grid in WGS84 - BIG ENDIAN” theo địa chỉ: http:/earthinfo.nga.mil/GandG/wgs84/gravitymod/egm 2008/und_min 2.5x2.5_egm2008_isw=82_WGS84_Tid eFree.gz Đối với ứng dụng hải dương học, tổ chức NGA cung cấp CSDL mơ hình dị thường độ cao toàn cầu EGM2008 dạng mạng lưới chuẩn hình vng độ phân giải 1’ x 1’ bao phủ toàn cầu Các liệu đỉnh ô chuẩn bao gồm tọa độ trắc địa B, L (tương ứng với ellipsoid trái đất trung  Mạng lưới (grid) chuẩn hình vng độ phân giải 1’ x 1’ bao phủ toàn cầu Các liệu đỉnh ô chuẩn bao gồm tọa độ trắc địa B, L (tương ứng với bình TOPEX/Poseidon) dị thường độ cao ellipsoid quy chiếu chung WGS84) dị Mơ hình địa hình động lực trung bình MDT thực chất mơ hình địa hình mặt biển trung bình biển đại dương giới với giá trị MDT điểm mặt thường độ cao toàn cầu  n hệ khơng phụ thuộc triều tồn cầu  m hệ triều trung bình 2.1.2 Mơ hình địa hình mặt biển trung bình tồn cầu Dữ liệu nêu truy nhập từ Website “EGM2008-WGS84 Version” với TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu biển trung bình độ cao điểm so với mặt geoid toàn cầu Trung tâm Vũ trụ quốc gia Đan Mạch (DNSC) đổi tên thành Viện Vũ trụ Quốc gia (NSI) trực thuộc Trường Đại học Tổng hợp Kỹ thuật Đan Mạch (DTU) thành lập mơ hình địa hình động lực trung bình DNSC08MDT (Mean Dynamic Topography) dựa mơ hình mặt biển trung bình DNSC08MSS (Mean Sea Surface) mơ hình trọng trường EGM2008 [1], [2] Mơ hình DNSC08MSS xây dựng hệ triều trung bình từ dự án vệ tinh TOPEX/Poseidon, JASON-1, GEOSAT GM, GFO ERM, ERS-1 GM, ERS-2 ERM, ENVISAT ERM, ICESAT giai đoạn từ 1993 - 2004 Các giá trị SSH mơ hình MSS thực chất độ cao trắc địa H điểm mặt biển trung bình tồn cầu (mặt địa hình biển) so với ellipsoid trung bình TOPEX/Poseidon (hoặc GRIM), thêm vào độ cao trắc địa H xác định dựa liệu altimetry thu nhận từ dự án vệ tinh altimetry Mơ hình DTU10MSS xây dựng tương ứng với ellipsoid trung bình TOPEX/Poseidon hệ triều trung bình Mơ hình địa hình động lực trung bình DTU10MDT với độ phân giải 1’ x 1’ bao phủ toàn cầu xây dựng dựa mơ hình mặt biển trung bình DTU10MSS mơ hình EGM2008 [1] Mơ hình phát triển từ mơ hình DNSC08MSS theo cơng thức MDTn = H M −  M , H M độ cao trắc địa điểm M mặt biển trung bình tồn cầu so với mặt ellipsoid trung bình TOPEX/Poseidon lấy từ mơ hình DNSC08MSS DTU10MSS, cịn dị thường độ cao tồn cầu  M hệ không phụ thuộc triều điểm M tính từ mơ hình geoid tồn cầu tương ứng với mơ hình EGM2008 (xem hình 1) Hình 1: Độ cao mặt địa hình biển MDTn Các giá trị MDTn điểm M mặt địa hình biển từ mơ hình MDT khơng phụ thuộc vào ellipsoid, mà phụ thuộc vào mặt geoid tồn cầu Bởi mơ hình DNSC08MDT, DTU10MDT xác định tương ứng với mơ hình EGM2008, nên chúng xây dựng tương ứng với mặt geoid tồn cầu trọng trường thực W0 = 62636856 ,0 m s −2 Về phần mình, giá trị MDTn độ cao điểm M mặt biển trung bình tồn cầu so với mặt geoid tồn cầu MDTn hệ khơng phụ thuộc triều độ Trong năm 2013, 2015 2019, Tổ chức DTU tiếp tục phát triển mô hình địa hình động lực trung bình tồn cầu DTU13MDT, DTU15MDT DTU19MDT dựa geoid toàn cầu EIGEN6C1 EIGEN6C4 cao điểm mặt địa hình biển so với mặt geoid tồn cầu mơ hình DNSC08MDT DTU10MDT tính Mơ hình DTU13MDT xây dựng năm 2013 sở sử dụng mô hình DTU13MSS geoid EIGEN6C1 Mơ hình Theo tài liệu [1], [2], giá trị TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu DTU13MSS xây dựng sở sử dụng số liệu đo cao vệ tinh thời gian 20 năm Độ phân giải mơ hình DTU13MDT 1’ x 1’, độ bao phủ từ 900 Nam đến 900 Bắc [3] Mơ hình DTU15MDT xây dựng từ mơ hình DTU15MSS mơ hình geoid tồn cầu EIGEN6C4 So với mơ hình DTU13MSS, mơ hình DTU15MSS sử dụng thêm số liệu vệ tinh Cryosat-2 để làm tăng độ phân giải Ngồi ra, mơ hình geoid tồn cầu EIGEN6C4 sử dụng thêm số liệu gradient trọng lực vệ tinh GOCE nên độ xác tốt [4] Năm 2021, phiên mơ hình MDT DTU19MDT mơ hình DTU21MSS, DTU21LAT, Bathymetry (Gebco 14) giáo sư Ole Baltazar Andersen thuộc Viện Vũ trụ Quốc gia, Trường Đại học Tổng hợp kỹ thuật Đan Mạch (DTU) cơng bố Các mơ hình xây dựng dựa ellipsoid quy chiếu Topex/Poseidon mơ hình geoid EGM2008 2.1.3 Mơ hình số độ cao DEM Cục Bản đồ, Bộ Tổng tham mưu Mơ hình số độ cao (DEM) tỷ lệ 1:25.000 phủ trùm phần đất liền lãnh thổ Việt Nam Cục Bản đồ, Bộ Tổng tham mưu xây dựng từ năm 2014-2018 phương pháp đo vẽ ảnh lập thể trạm ảnh số Mô hình xậy dựng Hệ tọa độ VN-2000, Hệ độ cao quốc gia HP72 Kích thước mắt lưới 0,3” x 0,3”, số mảnh đặc biệt có kích thước mắt lưới 0,15” x 0,15” Mức độ chi tiết thành lập DEM phụ thuộc vào độ dốc tính phức tạp địa hình Độ xác DEM đáp ứng yêu cầu thành thành lập đồ địa hình tỷ lệ 1:25.000 tồn lãnh thổ DEM quản lý theo mảnh đồ tỷ lệ 1:25.000 có tên phiên hiệu tương ứng, mảnh có định dạng Geotiff BIL 2.1.4 Bản đồ địa hình đáy biển Trung tâm Trắc địa Bản đồ biển Bản đồ địa hình đáy biển Trung tâm Trắc địa Bản đồ biển thành lập dựa hệ tọa độ VN2000, phép chiếu UTM, hệ độ cao quốc gia HP72 dãy tỷ lệ 1:5.000, 1:10.000 cho đảo khu vực kinh tế trọng điểm gần bờ, 1:50.000 cho vùng biển ven bờ (phạm vi phủ đến hết vùng tiếp giáp lãnh hải) Bản đồ địa hình đáy biển tỷ lệ 1:100.000 biên tập từ đồ tỷ lệ 1:50.000 phương pháp tổng qt hóa Bản đồ địa hình đáy biển tỷ lệ thành lập cho khu vực gần bờ (trên phạm vi khoảng 100 km tính từ bờ biển) 2.1.5 Hải đồ Đoàn Đo đạc Biên vẽ hải đồ Nghiên cứu biển Hải đồ vùng biển Việt Nam Đoàn Đo đạc Biên vẽ hải đồ Nghiên cứu biển (Đoàn 6) thành lập dựa hệ tọa độ quốc gia VN2000 (2002 - 2009) quốc tế WGS84 (2010 đến nay), phép chiếu Mercator, hệ độ cao tính theo mực nước triều thiên văn thấp cho dãy tỷ lệ: 1:1.000, 1:2.000, 1:5.000, 1:10.000 1:25.000 khu vực cảng biển, đảo/đá; 1:100.000, 1:1:200.000, 1:300.000, 1:400.000, 1:500.000, 1:1.000.000, 1:1.500.000 1:2.500.000 cho toàn vùng biển Việt Nam 2.1.6 Độ cao mặt biển trạm quan trắc mực nước Độ cao mặt biển trạm quan trắc mực nước Trung tâm Hải văn, Tổng cục Biển hải đảo Việt Nam Đoàn Đo đạc Biên vẽ hải đồ Nghiên cứu biển đo đạc tính tốn phương pháp điều hịa thủy triều, phương pháp Vladimirsky phương pháp Peresipkin TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu 2.2 Đề xuất phương án sử dụng thống hệ tọa độ, độ cao để xử lý liệu Từ kết phân tích nguồn liệu mục 2.1, thấy rằng: - Dữ liệu có nguồn gốc toàn cầu: Được qui chiếu dựa hệ tọa độ, độ cao toàn cầu mà đặc trưng hệ tọa độ tồn cầu WGS84 mơ hình EGM2008 - Dữ liệu có nguồn gốc nước: Chủ yếu qui chiếu dựa Hệ tọa độ VN2000 hệ độ cao quốc gia HP72 Chỉ có liệu hải đồ quy chiếu dựa hệ tọa độ quốc gia VN2000 (2002 - 2009) quốc tế WGS84 (2010 đến nay) độ cao mực nước triều thiên văn thấp Như vậy, để thống sở toán học liệu, cần phải lựa chọn hệ quy chiếu tọa độ phù hợp Các phương án lựa chọn hệ quy chiếu tọa độ để xử lý liệu sau:  Phương án Lựa chọn Hệ quy chiếu tọa độ quốc gia VN2000 hệ độ cao quốc gia HP72 Hệ quy chiếu tọa độ quốc gia VN2000 hệ độ cao HP72 quen thuộc với chuyên ngành Việt Nam thuận tiện cho người sử dụng Đồng thời có tham số chuyển đổi tọa độ hệ tọa độ quốc gia VN2000 hệ tọa độ quốc tế WGS84 Tuy nhiên, tham số chuyển đổi tọa độ xây dựng dựa điểm tọa độ Nhà nước phân bố lãnh thổ đất liền Việt Nam Tuy nhiên, chuyển đổi tọa độ hai hệ thống nêu cho liệu thuộc phạm vi vùng biển Việt Nam kết khơng cịn nữa, đặc biệt cho vùng biển xa bờ Độ cao điểm lãnh thổ Việt Nam xác định chủ yếu phương pháp thủy chuẩn hình học hay GNSS dựa mơ hình VNGeo Cục Đo đạc Bản đồ Thông tin địa lý xây dựng Mơ hình VNGeo xây dựng từ mơ hình tiên nghiệm EGM2008, làm khớp với điểm GNSS-Thủy chuẩn phân bố lãnh thổ Việt Nam phương pháp sóng phương pháp phần dư Vì vậy, sử dụng mơ hình VNGeo để chuyển đổi độ cao tồn cầu dựa mơ hình EGM2008 độ cao chuẩn Việt Nam cho liệu thuộc vùng biển Việt Nam khơng đảm bảo độ tin cậy Vì vậy, sử dụng Hệ tọa độ quốc gia VN2000 hệ độ cao HP72 để xử lý liệu phục vụ xây dựng mơ hình mặt biển phạm vi vùng biển Việt Nam không phù hợp  Phương án Lựa chọn Hệ quy chiếu tọa độ toàn cầu WGS84 hệ độ cao toàn cầu Hệ tọa độ tồn cầu WGS84 đảm bảo độ xác đồng phạm vi lãnh thổ Việt Nam (cả đất liền biển) Tuy nhiên, kết nối với độ cao chuẩn phạm vi lãnh thổ đất liền Việt Nam lại không phù hợp Mặt khác, hệ độ cao toàn cầu chưa sử dụng để sản xuất tư liệu địa hình Việt Nam Vì thế, lựa chọn phương án có nhiều hạn chế áp dụng kết nghiên cứu vào thực tế  Phương án Lựa chọn Hệ quy chiếu tọa độ toàn cầu WGS84 hệ độ cao quốc gia HP72 Khi sử dụng Hệ tọa độ toàn cầu WGS84 gắn với hệ độ cao quốc gia HP72 khu vực đất liền lãnh thổ Việt Nam có phương án chuyển đổi tọa độ từ Hệ tọa độ toàn cầu WGS84 hệ tọa độ quốc gia VN2000 (theo tham số Bộ Tà ngun Mơi trường cơng bố năm 2007), cịn độ cao chuyển đổi Như vậy, sử dụng Hệ tọa độ toàn cầu WGS84 gắn với hệ độ cao quốc gia HP72 để xử lý liệu đảm bảo độ xác đồng tồn vùng biển Việt Nam, đồng thời thuận lợi việc áp dụng kết nghiên cứu đề tài vào thực tế Đây TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu phương pháp phù hợp để xử lý liệu để xây dựng mơ hình mặt biển đó, phương pháp chuyển đổi độ sâu Hải đồ (dựa mực nước triều thiên văn thấp nhất) hệ độ cao quốc gia HP72 công bố công trình [5], [7] Tính tốn thực nghiệm Số lượng nguồn liệu phục vụ xây dựng số mơ hình mặt biển lớn, chuyển đổi hệ tọa độ WGS84, hệ độ cao HP72 để xây dựng số mơ hình mặt biển tác giả cơng bố số cơng trình [5], [7], nên lựa chọn hai loại liệu điển hình Bản đồ địa hình đáy biển Hải đồ để thực chuyển đổi toạ độ, độ sâu theo Phương án “Lựa chọn Hệ quy chiếu tọa độ toàn cầu WGS84 hệ độ cao quốc gia HP72” Các thông tin thông số đầu vào đầu liệu thực nghiệm thể hình 3.1 Khái quát số liệu thực nghiệm * Toạ độ điểm đo sâu đồ địa hình đáy biển (phiên hiệu E-49-85-A) tỷ lệ 1:50.000 khu vực Vụng Chân Mây phạm vi từ 16015’N÷16030’N 0 108 00’E÷108 15’E Bản đồ thành lập hệ toạ độ VN2000, phép chiếu UTM, kinh tuyến trục 1110, múi chiếu 60, Hệ độ cao HP72 Tổng số điểm toạ độ chuyển đổi Hệ quy chiếu toạ độ (WGS84, phép chiếu UTM) 779 điểm * Toạ độ, độ sâu điểm đo sâu Hải đồ tỷ lệ 1:25.000 khu vực Vụng Chân Mây phạm vi từ 16016’18”N÷16025’36”N 107050’48”E÷108003’54”E Hải đồ thành lập hệ toạ độ VN2000, phép chiếu Mercator, vĩ tuyến chuẩn 160, độ sâu theo mực nước triều thiên văn thấp Tổng số điểm toạ độ chuyển đổi Hệ quy chiếu toạ độ WGS84, phép chiếu UTM, hệ độ cao quốc gia HP72 1668 điểm Hình 2: Sơ đồ chuyển đổi toạ độ theo phương án lựa chọn Theo phân tích trình bày mục 2, giá trị độ sâu đồ địa hình đáy biển tính dựa vào hệ độ cao HP72, giá trị độ sâu đồ địa hình đáy biển không cần 3.2 Kết thực nghiệm thực chuyển đổi báo Trong Bảng 1: Kết chuyển đổi toạ độ điểm độ sâu đồ địa hình đáy biển STT … 100 101 102 103 VN2000, phép chiếu UTM, KTT 111 độ, MC độ X (m) Y (m) h_HP72 (m) 1799591.700 1800558.000 1801963.300 1800894.800 1802544.400 … 188480.854 191741.400 199727.800 203446.100 191764.500 … 14 19 28 33 20 1821089.400 1821080.900 1821053.500 1821008.800 191440.100 193392.300 195044.000 196995.900 42.6 45.8 48.4 51.3 WGS84, phép chiếu UTM B (0) L (0) 16.25606046 16.26520198 16.27889454 16.26970522 16.28313950 … 108.0875054 108.1178529 108.1923244 108.2272204 108.1178068 … 16.4505329 16.45070716 16.45067089 16.45051533 108.1123101 108.1305758 108.1460329 108.1643012 h_HP72 (m) 14 19 28 33 20 … 42.6 45.8 48.4 51.3 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Nghiên cứu STT 104 … 201 202 203 204 205 … 775 776 777 778 779 VN2000, phép chiếu UTM, KTT 111 độ, MC độ X (m) Y (m) h_HP72 (m) WGS84, phép chiếu UTM B (0) L (0) h_HP72 (m) 1820994.300 198647.600 … … 53.9 16.45059304 … 108.1797569 … 53.9 … 1825201.100 1825203.800 1825209.500 1823656.200 1823658.000 … 186064.600 185163.400 184262.300 189020.600 188119.600 … 40.1 38.9 37.7 42.8 41.7 16.48695424 16.48685991 16.48679234 16.47339344 16.47329223 … 108.0614591 108.0530261 108.0445937 108.0893278 108.0808972 … 40.1 38.9 37.7 42.8 41.7 … 1808224.200 1805672.900 1805738.900 1804664.900 1803699.600 180070.500 184711.100 183895.900 184448.500 184397.400 12.3 12.7 12.5 8.4 8.1 16.33289666 16.3104747 16.31096395 16.30134007 16.29261857 108.0077187 108.0514516 108.0438216 108.049133 108.0487859 12.3 12.7 12.5 8.4 8.1 Bảng 2: Kết chuyển đổi toạ độ điểm độ sâu Hải đồ STT VN2000, Phép chiếu Mercator, Vĩ tuyến chuẩn 16 độ X (m) Y (m) 1771060.491 1769412.028 1768877.537 1770891.922 1770732.784 … 1767032.391 1768064.148 1766811.743 1766581.675 1766286.197 … 826160.953 827598.983 827723.148 826360.575 826611.037 … 825532.461 824792.309 825854.634 826172.527 826597.461 … … 61 62 63 64 65 … 441 1766683.027 814120.59 442 1770468.82 810140.539 443 1769219.73 811846.28 444 1769206.01 810578.49 445 1769976.98 809721.14 … … … 1664 1758129.324 805928.915 Độ sâu theo MBTN h (m) WGS84, phép chiếu UTM B (0) L (0) 16.42438885 16.40952383 16.40470384 16.42286883 16.42143383 … 16.38806376 16.39736879 16.38607374 16.38399875 16.38133379 … 108.0490738 108.0625089 108.0636689 108.0509388 108.0532788 … 108.0432016 108.0362867 108.0462116 108.0491816 108.0531516 … 16.38491381 16.41905492 16.40779099 16.40766737 16.41461982 … 0.1 16.30775373 107.9365833 107.8993989 107.9153351 107.9034905 107.8954805 … 107.8600497 29.5 29.5 29.5 29.5 29.5 … 25.5 25.5 25.5 25.5 25.5 … 18.2 18.8 17.8 17.8 17.8 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 Hệ độ cao HP-72 (m) Độ sâu Độ sâu theo theo MBTN MBTB -0.474 0.033 -0.481 0.032 -0.482 0.032 -0.474 0.033 -0.481 0.032 … … -0.479 0.032 -0.475 0.032 -0.480 0.032 -0.482 0.032 -0.484 0.032 … … -0.459 -0.449 -0.452 -0.451 -0.450 … -0.475 0.034 0.037 0.036 0.036 0.037 … 0.035 Nghiên cứu STT VN2000, Phép chiếu Mercator, Vĩ tuyến chuẩn 16 độ X (m) Y (m) 1665 1757387.86 806948.955 1666 1760763.998 806541.663 1667 1762332.619 814921.399 1668 1756790.667 825880.198 Độ sâu theo MBTN h (m) WGS84, phép chiếu UTM B (0) L (0) 0.1 16.30106364 107.8695796 0.1 16.33152368 107.8657747 0.5 16.3456737 107.9440647 0.7 16.29567367 108.0464496 * Sau chuẩn hoá chuyển toạ độ điểm đo sâu địa hình đáy biển hải đồ hệ quy chiếu toạ độ WGS84, phép chiếu UTM hệ độ cao quốc gia HP72, tiến hành biên tập lớp độ sâu hai loại đồ với phạm vi ranh giới mảnh hải đồ tỷ lệ 1:25.000 Sản phẩm đồ lớp độ sâu sau biên tập trình bày hình Hệ độ cao HP-72 (m) Độ sâu Độ sâu theo theo MBTN MBTB -0.477 0.034 -0.468 0.035 -0.469 0.034 -0.508 0.032 Trên sở phân tích ưu - nhược điểm lựa chọn sở toán học để xử lý liệu, báo đề xuất lựa chọn Hệ tọa độ toàn cầu WGS84 làm sở để tính tốn toạ độ mặt Hệ độ cao quốc gia HP72 để làm sở tính tốn độ cao. Tài liệu tham khảo [1] Andersen, O.B., Knudsen, P., 2008 The DNSC08MDT Mean Dynamic Topography EGU 2008 Meeting, Vienna, Austria, April 14-18 [2] Andersen, O.B., Knudsen, P., 2009 DNSC08 mean sea surface and mean dynamic topography models Journal of geophysical Research, Vol 114, C11001, 12 pp., 2009, doi:10.1029/2008JC005179 Hình 3: Tọa độ, độ sâu đồ địa hình đáy biển tọa độ, độ sâu hải đồ hệ thống toạ độ WGS84, phép chiếu UTM, hệ độ cao HP72 Kết luận Bài báo nghiên cứu sở toán học liệu, mối quan hệ liệu toàn cầu liệu quốc gia, đồng thời đưa số phương án lựa chọn sở toán học phù hợp để xử lý liệu phục vụ cơng tác nghiên cứu xây dựng mơ hình mặt biển [3] Andersen O.B., Knudsen P., Stenseng L., 2013 The DTU13 Global mean sea surface from 20 years satellite altimetry DTU Space Technical University of Denmark [4] Nguyễn Văn Sáng, 2016 Nghiên cứu phương pháp xác định độ cao địa hình mặt biển số liệu đo cao vệ tinh Biển Đông, Đề tài khoa học công nghệ cấp Bộ, mã số: B2014-02-18 (Xem tiếp trang 17) TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 53-9/2022 ... nguồn liệu phục vụ xây dựng số mơ hình mặt biển lớn, chuyển đổi hệ tọa độ WGS84, hệ độ cao HP72 để xây dựng số mô hình mặt biển tác giả cơng bố số cơng trình [5], [7], nên lựa chọn hai loại liệu. .. vậy, sử dụng Hệ tọa độ quốc gia VN2000 hệ độ cao HP72 để xử lý liệu phục vụ xây dựng mô hình mặt biển phạm vi vùng biển Việt Nam không phù hợp  Phương án Lựa chọn Hệ quy chiếu tọa độ toàn cầu WGS84... liệu, mối quan hệ liệu toàn cầu liệu quốc gia, đồng thời đưa số phương án lựa chọn sở toán học phù hợp để xử lý liệu phục vụ công tác nghiên cứu xây dựng mơ hình mặt biển [3] Andersen O.B., Knudsen