倾斜摄影osgb格式文件，进行坐标转换

倾斜摄影OSGB格式的文件，很棘手，今天需要把osgb放到UE中渲染。碰到的问题如下：

1、osgb文件导进去后，Z轴不想上，是歪的，小人放进去后，就斜站在马路上。

2、根本原因是坐标系，UE的插件cesium for UE支持WGS84坐标（wkid:4326）。

怎么解决问题呢？

1、当然是问GPT呀，搞搞问答，看它出什么答案，结果，很拉跨。

2、直接转换：通过osgbLab软件进行osgb坐标系变换，用了，能成功，但是，是假性的成功。用其他软件打开看，会发现转换的根本不对，跟极坐标似的。

3、GDAL，一个不错的库，然而，也不行。

问题来了：怎么才行？

当然是间接投影呀！

1、转换osgb为obj格式。使用osgbLab工具，把倾斜摄影的成功转成一个obj，大几百MB。

2、obj进行格式转换，见代码。注意修改input.obj与output.obj。核心功能：找到obj代码里的坐标部分，进行坐标系变换。

# -*- coding: utf-8 -*-
import pyproj
X0, Y0 = 1111111,111111111 #此处需要通过metadata.xml获取，我发现obj里的xy坐标，与metadata.xml结合才能使用
# 初始化坐标系转换器
transformer = pyproj.Transformer.from_crs("EPSG:4528", "EPSG:4326", always_xy=True)

with open('input.obj', 'r') as infile, open('output.obj', 'w') as outfile:
    for line in infile:
        if line.startswith('v '):  # 处理顶点信息
            parts = line.split()
            x, y, z = map(float, parts[1:4])
            # 执行坐标系转换

            x, y = transformer.transform(x+X0, y+Y0)

#            x, y = transformer.transform(x, y)
            outfile.write(f"v {x} {y} {z}\n")
        else:
            outfile.write(line)

3、使用osgbLab进行格式转换，obj转osgb。至此，曲线救国完成。

相关说明：

obj格式简介 https://www.cnblogs.com/dzmsjs/p/18389082