2. 海洋数据中的坐标系基础

在海洋数据处理中,坐标是连接“数值”与“真实海域位置”的桥梁。同一组温度、盐度或水深数据,如果坐标参照不同,落到地图上可能相差几十米甚至数百米。理解坐标系,核心是分清两类表达:地理坐标系用经纬度描述球面上的点,投影坐标系用平面直角坐标描述展开后的地图上的点。本节围绕这两类坐标系的概念、常用海洋坐标系、坐标表示形式及其在海洋数据中的意义展开。

一、地理坐标系:以椭球和角度定位

地理坐标系(Geographic Coordinate System,GCS)由“基准面(椭球体+定位定向)+ 角度单位”构成。它用经度 λ、纬度 φ 表示点位,单位是度(或度分秒)。经度以本初子午线为 0°,向东为正、向西为负,范围 −180°~180°;纬度以赤道为 0°,向北为正、向南为负,范围 −90°~90°。这一正负号约定在后续所有图示与计算中保持一致。

需要特别强调:经纬度必须依附于某个基准面才有意义。同一个物理点,在 WGS84 椭球下和在 CGCS2000 椭球下算出的经纬度会有细微差异。海洋数据中常见的基准面包括:

  • WGS84:GPS 的默认基准,也是全球海洋观测、AIS、浮标、卫星高度计数据最通用的坐标参照,长半轴约 6378137 m,扁率 1/298.257223563。
  • CGCS2000:我国 2000 国家大地坐标系,与 WGS84 在厘米级上接近,是我国海洋调查、海域管理和基础测绘的法定基准。
  • ITRF 框架下的历元坐标:高精度 GNSS 与部分科研数据使用,需注明参考历元(如 ITRF2014, 2010.0)。

实际处理中,若数据源标注为“WGS84 经纬度”,而业务底图使用 CGCS2000,二者在多数海洋应用中可以近似等同,但在高精度岸线、岛礁测绘中应做基准转换并记录转换参数,不能默认“经纬度一样就是同一个点”。

二、投影坐标系:把球面展开成平面

投影坐标系(Projected Coordinate System,PCS)= 地理坐标系 + 投影方法 + 线性单位。它把椭球面上的点按某种数学规则映射到平面上,用 (x, y) 或 (E, N) 表示,单位通常为米。海洋数据常用的投影有:

  • 墨卡托投影(Mercator):等角,航向为直线,是海图与航海导航的经典选择,但高纬度面积变形大。
  • 横轴墨卡托 / UTM:分带投影,每带 6°,适合区域尺度海洋工程与调查数据,变形小。
  • 高斯-克吕格投影:我国大比例尺地形图与海域管理常用,按 3° 或 6° 分带,需注明带号与中央经线。
  • 兰勃特等角圆锥投影:适合中纬度东西向延伸的海区。

投影坐标的关键参数是中央经线、原点纬度和假东/假北。例如高斯-克吕格 3° 带第 40 带,中央经线为 120°E,为避免负值,东坐标通常加 500000 m 的假东;我国还常在 y 坐标前加带号。读取一个平面坐标时,必须同时知道“哪个投影、哪个带、什么单位”,否则数值无法还原为地理位置。

地理坐标系 → 投影 → 平面坐标:海洋点位表达的转换关系 P(φ, λ) 赤道 0° 本初子午线 0° 地理坐标系(GCS) 基准面:WGS84 / CGCS2000 单位:度(°),经度东正西负,纬度北正南负 投影变换 墨卡托 / UTM / 高斯-克吕格 P(x, y) 投影坐标系(PCS) 投影方法 + 中央经线 + 假东/假北 单位:米(m),如 x=500000+ΔE 同一物理点:经纬度 (φ, λ) 依附基准面;平面坐标 (x, y) 依附投影与带号。 海洋数据入库时必须同时记录:坐标系名称、基准面、投影参数、单位、参考历元(如适用)。

三、坐标表示形式及其在海洋数据中的意义

海洋数据中常见的坐标表示有三类:经纬度(度或度分秒)、平面直角坐标(米)、以及栅格行列号。经纬度便于全球统一表达和跨数据集拼接,是 NetCDF、卫星产品、AIS 报文的主流形式;平面坐标便于距离、面积计算和工程放样,常用于海图、海域使用、管线与平台定位;栅格行列号则是遥感与模型网格的内部索引,必须通过仿射变换参数才能还原为地理坐标。

以一份示例海温格点数据为例(以下为教学假设参数):数据为 WGS84 经纬度,经度范围 118°E~124°E,纬度 22°N~28°N,分辨率 0.05°,共 121×121 个格点。若要把某观测站 (121.50°E, 25.30°N) 对应到格点,可按 col = (121.50−118)/0.05 = 70,row = (28−25.30)/0.05 = 54 计算(行号自北向南递增),得到行列号 (54, 70)。这里的关键是明确格点代表的是像元中心还是角点,以及行序方向,否则会整体偏移半个格点。

若业务底图采用高斯-克吕格 3° 带第 40 带(中央经线 120°E),则需把上述经纬度投影为平面坐标,再与底图叠加。此时应检查:带号是否为 40、假东是否已加 500000 m、单位是否为米。若把未加带号的 y 值与加了带号的底图混用,会出现数百公里的错位,这是海洋 GIS 数据叠加中最常见的错误之一。

判断坐标是否“对得上”的实用方法:取一个已知地标(如某灯塔或浮标),分别用两套坐标计算其位置,若平面距离偏差在米级且方向一致,说明基准与投影匹配;若偏差呈系统性平移,多半是基准面或假东设置不一致。

四、实际处理中的关键约定

在海洋数据流程中,建议遵循以下约定:第一,元数据先行,每个数据集必须写明坐标系名称(如 WGS84 / CGCS2000 / CGCS2000_3_Degree_GK_Zone_40)、单位与参考历元;第二,转换留痕,记录从源坐标到目标坐标所用的转换参数与工具版本;第三,方向一致,经度东正西负、纬度北正南负,平面坐标明确 x 向东、y 向北;第四,精度匹配,全球尺度可用 WGS84 经纬度,区域工程用投影坐标,高精度岸线需做基准转换。

坐标系的正确使用,直接决定海洋数据能否在空间上“对齐”。把地理坐标系理解为球面上的角度定位,把投影坐标系理解为平面上的米制定位,并在每个环节保留完整的坐标元数据,才能让温度、盐度、水深、航迹等数据在同一张海图上可靠地叠加与分析。