libisg

库读取/写入 ISG 2.0 格式。

use std::fs;

use libisg;
use libisg::{Data, DataBounds, ISG};

let s = fs::read_to_string("Example 1.isg").unwrap();

let isg = libisg::from_str(&s).unwrap();


let (a_max, b_max, delta_a, delta_b) = match isg.header.data_bounds {
    DataBounds::GridGeodetic { lat_max, lon_max, delta_lat, delta_lon, .. } => {
        (lat_max, lon_max, delta_lat, delta_lon)
    },
    _ => unimplemented!("`Example 1.isg` is grid geodetic"),
};

match &isg.data {
    Data::Grid(data) => {
        for (nrow, row) in data.iter().enumerate() {
            for (ncol, value) in row.iter().enumerate() {
                let a = a_max - delta_a * nrow;
                let b = b_max - delta_b * ncol;
                // do something
            }
        }
    }
    Data::Sparse(data) => {
        for row in data {
            let (a, b, value) = row;
            // do something
        }
    }
}

特性

• 支持 ISG 格式的序列化/反序列化
• 支持 serde（需要 serde 功能）

许可证

MIT 或 Apache-2.0

参考

规范：https://www.isgeoid.polimi.it/Geoid/format_specs.html