Planetiler中GeoPackageReader处理EPSG:3067坐标系的坐标交换问题分析

问题背景

在使用Planetiler处理来自芬兰国家土地调查局(NLS)的GeoPackage数据时，发现了一个坐标系统处理的问题。该数据使用的是ETRS-TM35FIN投影系统(对应EPSG:3067)，其坐标顺序(北向、东向)与常见的EPSG:4326(WGS84)坐标系(纬度、经度)顺序相反。

问题现象

当使用Planetiler的GeoPackageReader读取这种特殊坐标系的数据时，会导致生成的矢量切片中出现经纬度坐标交换的情况。具体表现为：

1. 地图上的经纬度坐标被交换
2. 地图整体呈现90度逆时针旋转
3. 随后又进行了垂直轴的镜像翻转
4. 最终地图在垂直方向上被拉伸

这种坐标转换错误导致最终生成的地图显示异常，无法正常使用。

技术分析

问题的根源在于GeoPackageReader在处理不同坐标系转换时的逻辑。在原始代码中，当从EPSG:3067转换到EPSG:4326时，没有正确处理坐标轴的顺序差异。

ETRS-TM35FIN(EPSG:3067)是芬兰使用的横轴墨卡托投影，其坐标顺序为：

• 第一个值：北向坐标(Northing)
• 第二个值：东向坐标(Easting)

而WGS84(EPSG:4326)的坐标顺序为：

• 第一个值：纬度(Latitude)
• 第二个值：经度(Longitude)

解决方案

通过修改GeoPackageReader.java中的坐标转换逻辑可以解决此问题。关键修改包括：

1. 明确指定目标坐标系为EPSG:4326
2. 使用CRS.findMathTransform方法进行坐标转换时，设置lenient参数为true以允许宽松的转换

修改后的核心代码如下：

var latLonCRS = CRS.decode("EPSG:4326", true);
coordinateTransform = CRS.findMathTransform(sourceCRS, latLonCRS, true);

后续发展

根据项目维护者的反馈，此问题已在Planetiler 0.9版本中得到修复。修复可能涉及对坐标转换逻辑的全面改进，而不仅仅是针对EPSG:3067这一特定情况。

经验总结

处理地理空间数据时，坐标系的差异是需要特别注意的关键点。开发者和用户应当：

1. 了解数据源使用的具体坐标系及其特性
2. 注意不同坐标系之间坐标轴顺序的差异
3. 在坐标转换时明确指定源坐标系和目标坐标系
4. 测试生成的矢量切片是否正确保持了原始数据的空间关系

对于使用Planetiler处理非标准坐标系数据的用户，建议：

1. 升级到最新版本(0.9或更高)
2. 如果仍遇到类似问题，可以检查坐标转换逻辑
3. 必要时可以手动指定坐标转换参数

通过正确处理坐标系转换，可以确保生成的矢量切片保持准确的地理空间关系，为地图应用提供可靠的基础数据。