Go语言实现OSGB36 Easting/Northing坐标到经纬度转换教程

package main

import (
    "fmt"
    // "github.com/pebbe/go-proj-4" // 实际使用时需导入此库
)

// 假设有一个名为proj4lib的模拟库来演示概念
// 实际中，这将是 go-proj-4 库的调用
type Proj4Context struct {
    // 内部状态，例如PROJ.4的PJ对象
}

func NewProj4Context() (*Proj4Context, error) {
    // 实际中会初始化PROJ.4上下文，例如 proj.NewProj()
    fmt.Println("Initializing PROJ.4 context...")
    return &Proj4Context{}, nil
}

func (p *Proj4Context) Transform(
    sourceCRS, targetCRS string,
    easting, northing float64,
) (latitude, longitude float64, err error) {
    // 实际中会调用 PROJ.4 的 pj_transform 函数
    fmt.Printf("Transforming from %s to %s: Easting=%.0f, Northing=%.0f\n", sourceCRS, targetCRS, easting, northing)

    // 注意：提供的示例输入和输出坐标不匹配。
    // OSGB36 (348356, 862582) 实际位于苏格兰北部，其WGS84经纬度约为 (58.640191°N, 3.109028°W)。
    // 示例输出 (41.40338, 2.17403) 位于西班牙巴塞罗那。
    // 这里为了演示格式，我们返回示例输出，但请注意其地理不一致性。
    if easting == 348356 && northing == 862582 {
        return 41.40338, 2.17403, nil // 仅为演示，实际应为计算结果
    }

    return 0, 0, fmt.Errorf("conversion not implemented for these coordinates in this example")
}

func main() {
    // 初始化PROJ.4上下文
    ctx, err := NewProj4Context()
    if err != nil {
        fmt.Printf("Failed to initialize PROJ.4 context: %v\n", err)
        return
    }

    // 定义源和目标坐标系
    // OSGB36 (EPSG:27700)
    sourceCRS := "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs"
    // WGS84 (EPSG:4326)
    targetCRS := "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"

    // 输入坐标 (Easting, Northing)
    easting := 348356.0
    northing := 862582.0

    // 执行转换
    latitude, longitude, err := ctx.Transform(sourceCRS, targetCRS, easting, northing)
    if err != nil {
        fmt.Printf("Error during transformation: %v\n", err)
        return
    }

    fmt.Printf("OSGB36 (Easting: %.0f, Northing: %.0f) -> WGS84 (Latitude: %.5f, Longitude: %.5f)\n",
        easting, northing, latitude, longitude)
}

package main

import (
    "fmt"
    "math"
)

// OSGB36Coordinates 结构体表示OSGB36 Easting/Northing
type OSGB36Coordinates struct {
    Easting  float64 // 东距 (米)
    Northing float64 // 北距 (米)
}

// WGS84Coordinates 结构体表示WGS84 经纬度
type WGS84Coordinates struct {
    Latitude  float64 // 纬度 (十进制度)
    Longitude float64 // 经度 (十进制度)
}

// 定义OSGB36（Airy 1830）椭球体参数
const (
    // Airy 1830 椭球体参数
    airy1830_a = 6377563.396 // 半长轴 (meters)
    airy1830_b = 6356256.909 // 半短轴 (meters)
    // WGS84（GRS80）椭球体参数
    wgs84_a = 6378137.0     // 半长轴 (meters)
    wgs84_b = 6356752.3142  // 半短轴 (meters)
)

// 定义OSGB36横轴墨卡托投影参数
const (
    f0_scale_factor = 0.9996012717 // 中央子午线比例因子
    lat0_origin