3. 海洋坐标系与投影选择
海洋数据有一个绕不开的前提:地球是曲面,而屏幕和纸张是平面。把海表温度、船舶轨迹、海底地形或风暴潮淹没范围放到QGIS里时,必须先回答两个问题——用哪个坐标参照系(CRS)描述位置,以及用哪种投影把曲面“摊平”。这两个选择直接决定了图上量算的距离、面积和方向是否可信。本节围绕地理坐标系与投影坐标系的区别、海洋常用投影的适用场景、QGIS中的定义与转换操作,以及投影变形对海洋数据的实际影响展开。
地理坐标系与投影坐标系:描述位置与描述平面是两件事
地理坐标系(GCS)用经纬度描述点位,单位是度,基准是椭球体。常见的WGS 84(EPSG:4326)就是地理坐标系,GPS定位、AIS报文、多数海洋再分析产品的原始坐标都基于它。地理坐标系的问题是:经纬度不是平面长度单位,一度经度在赤道约111公里,在纬度60°处只剩约55.5公里。因此直接在EPSG:4326下用QGIS测量距离或面积,结果会随纬度严重失真,不能作为工程量算依据。
投影坐标系(PCS)则通过数学变换把椭球面映射到平面,坐标单位是米,可以直接量算。它由“地理坐标系+投影方法+参数”共同定义。理解这一点很关键:同一个投影方法配不同中央经线,就是不同的CRS,EPSG代码也不同。海洋工作中常见的组合包括基于WGS 84的UTM分区、Web Mercator(EPSG:3857)以及各类Lambert投影。
海洋常用投影及其适用边界
Mercator(墨卡托)是等角投影,航向线(等角航线)在图上为直线,这正是航海图长期采用它的原因。但它的面积变形随纬度急剧放大,纬度80°附近的面积放大倍数可达数十倍,因此不适合做极区海冰范围统计或大范围面积对比。Web Mercator(EPSG:3857)是网络地图常用的球面近似版本,方便与在线底图叠加,但精度低于椭球版墨卡托,不宜用于高精度量算。
UTM(通用横轴墨卡托)按6°经差分带,每带一条中央经线,中央经线附近变形极小,是区域尺度海洋工程的首选。我国近海大致跨UTM 48N—52N带,例如渤海、黄海部分区域常用EPSG:32650(UTM 50N),南海部分区域用EPSG:32649(UTM 49N)。UTM的局限是跨带拼接时会出现明显错位,大范围成图需谨慎。
Lambert等角圆锥投影用两条标准纬线控制变形,适合中纬度东西延伸的大区域,如我国近海大范围专题图、海洋功能区划图。它保持角度不变形,面积变形在标准纬线之间较小,但离开标准纬线越远变形越大。
QGIS中的坐标系定义与转换
在QGIS中,每个图层都有自己的CRS,项目(Project)也有一个项目CRS。图层CRS描述数据本身“是什么坐标”,项目CRS决定“画在屏幕上用什么坐标”。二者不一致时,QGIS会实时进行动态投影(on-the-fly reprojection),让不同CRS的图层能叠在一起显示。但要注意:动态投影只改变显示,不改变数据本身。若要在某个CRS下做精确量算或导出成果,应使用“另存为”并指定目标CRS,把数据真正重投影。
实际操作路径:右键图层 →“导出”→“要素另存为”,在“CRS”下拉框选择目标坐标系(如EPSG:32650),勾选“添加保存的文件到地图”,即可得到重投影后的新图层。也可用“处理工具箱”中的“重投影图层”(Reproject layer)批量处理。判断是否成功,可查看图层属性中的CRS字段,并用测量工具在同一段已知距离上对比重投影前后的量算值。
示例:将南海某海域AIS轨迹从WGS 84重投影到UTM 49N
源CRS:EPSG:4326(经纬度)
目标CRS:EPSG:32649(UTM zone 49N,单位:米)
操作:图层右键 → 导出 → 要素另存为 → CRS选EPSG:32649
验证:对同一段航迹测量长度,重投影后单位应为米,且数值与球面距离接近
投影变形对海洋数据的实际影响
投影变形不是抽象概念,它会直接改变分析结论。以一次示例性对比说明:在EPSG:4326下直接计算某纬度40°附近一块1°×1°海域的面积,QGIS会按“度”的平面单位处理,得到的结果没有物理意义;若先重投影到UTM 50N再计算,面积单位变为平方米,数值才可用于统计。再如,用Web Mercator统计高纬度海冰覆盖面积,会因面积放大而系统性高估,纬度越高偏差越大。
另一个常见问题是跨带数据拼接。若把分别位于UTM 49N和50N的两块海底地形数据直接叠加,接缝处会出现明显错位或拉伸,因为两带的中央经线不同。正确做法是先统一到同一CRS(如Lambert或某一UTM带),或在大范围分析时改用适合的圆锥投影。此外,方向也有影响:在墨卡托图上量取的方位角是等角航线方位,与大地线方位不同,长距离导航规划时需注意区分。
把上述关系落到操作上,可以形成一条稳定的工作习惯:拿到海洋数据先确认其原始CRS(多数为EPSG:4326),明确本次任务是展示、量算还是成图;展示可保留地理坐标并叠加在线底图,量算则重投影到对应UTM带或Lambert;导出成果前检查图层CRS与项目CRS是否一致,避免“看起来对、算起来错”。投影选择没有唯一正确答案,只有与任务尺度、位置和精度要求相匹配的答案。