WGS-84与SK-42（EPSG:4284）坐标转换结果偏差问题求助

我目前拥有一套经测试可正常运行的Python代码，能实现WGS-84与SK-42（米单位）坐标参考系统（CRS）间的双向转换。

为验证转换准确性，我使用公开参考资源进行对比，以下EPSG代码对应的WGS-84 ↔ SK-42（米单位）转换结果符合预期：

• WGS-84：EPSG: 4326
• SK-42（米单位）：EPSG: 28405、28406、28407（适用于乌克兰区域）

但当使用GDAL或PYPROJ库实现WGS-84与SK-42（度单位）的双向转换时，对应以下EPSG代码的转换结果与参考源存在显著差异：

• WGS-84：EPSG:4326
• SK-42（度单位）：EPSG:4284

我怀疑问题出在EPSG:4284的使用上，但缺乏相关专业知识，特此附上测试代码及对比结果，寻求解决思路。

测试代码 #

import pandas as pd

dataset_file = './dataset.xlsx'
dataset_sheet = 'WGS84'
dataset = pd.read_excel(dataset_file, sheet_name=dataset_sheet)

sk42_code = 4284 # Pulkovo 42
wgs84_code = 4326 # WGS-84

#### CONVERSION WITH HELP OF PYPROJ

from pyproj import Transformer

transformer = Transformer.from_crs(wgs84_code, sk42_code)
    
res_pyproj = [transformer.transform(row.lon_WGS84, row.lat_WGS84) \
              for _, row in dataset.iterrows()]
res_pyproj = pd.DataFrame(res_pyproj, columns=['lon_SK42', 'lat_SK42'])

转换结果对比 #

各列含义：

• lon_WGS84、lat_WGS84：原始WGS-84数据
• lon_SK42、lat_SK42：上述代码的转换结果
• lon_SK42_ref、lat_SK42_ref：参考服务的转换结果

内容的提问来源于stack exchange，提问作者Olegh Bondarenko