Struct JgdTransformer 

pub struct JgdTransformer { /* private fields */ }

JGD2011座標変換器

Tokyo Datum/JGD2000からJGD2011への変換を行う。

• Step 1: Tokyo Datum → JGD2000 (TKY2JGD)
• Step 2: JGD2000 → JGD2011 (PatchJGD統合)
• Step 3: JGD2011 → WGS84 (proj)

Implementations

impl JgdTransformer

pub fn with_survey_date<P: AsRef<Path>>( tky_path: P, patch_dir: P, survey_date: SurveyStartDate, ) -> Result<Self, TransformError>

パラメータファイルを読み込んで変換器を初期化（調査日指定あり）

Arguments

• tky_path - TKY2JGD.par のパス
• patch_dir - PatchJGD系ファイルが入ったディレクトリのパス
• survey_date - 調査日（地震補正の適用判定に使用）

地震補正の適用ルール

• 調査日が地震発生日より前の場合のみ、その地震の補正を適用
• 調査日が地震発生日より後の場合、補正は適用しない（既に変動後の位置）
• ジオイド補正（hyokorev2014_geoid2011_h.par）は常に適用

Example

let survey_date = SurveyStartDate::new(2010, 5, 15);
let transformer = JgdTransformer::with_survey_date(
    "params/TKY2JGD.par",
    "params/patch_files",
    survey_date
)?;

Examples found in repository

examples/bench.rs (lines 18-22)

fn main() {
    // ファイル読み込み
    let xml_bytes = std::fs::read("107282 (1).xml").unwrap();
    let (utf8, _, _) = encoding_rs::SHIFT_JIS.decode(&xml_bytes);
    let xml_str = utf8.to_string();
    
    // パース計測
    let start = Instant::now();
    let boring = Boring300::parse_from_str(&xml_str).unwrap();
    println!("Parse: {:?}", start.elapsed());
    
    // 変換器初期化計測
    let start = Instant::now();
    let transformer = JgdTransformer::with_survey_date(
        "params/TKY2JGD.par",
        "params/patch_files",
        SurveyStartDate::new(2010, 5, 15),
    ).unwrap();
    println!("Transformer init: {:?}", start.elapsed());
    
    // 座標変換計測
    let start = Instant::now();
    let loc = boring.geo_location();
    let result = loc.to_wgs84(&transformer);
    println!("Coordinate transform: {:?}", start.elapsed());
    println!("Result: {:?}", result);
    
    // 100回ループ計測
    let start = Instant::now();
    for _ in 0..100 {
        let loc = boring.geo_location();
        let _ = loc.to_wgs84(&transformer);
    }
    println!("100x transform: {:?}", start.elapsed());
}

pub fn new<P: AsRef<Path>>( tky_path: P, patch_dir: P, ) -> Result<Self, TransformError>

パラメータファイルを読み込んで変換器を初期化（調査日なし）

注意: 調査日が不明な場合、地震補正は適用されません。 ジオイド補正のみが適用されます。

Arguments

• tky_path - TKY2JGD.par のパス
• patch_dir - PatchJGD系ファイルが入ったディレクトリのパス

pub fn new_tky_only<P: AsRef<Path>>(tky_path: P) -> Result<Self, TransformError>

TKY2JGDのみで変換器を初期化 (PatchJGDファイルなし)

JGD2000→JGD2011変換が不要な場合に使用

pub fn transform_horizontal( &self, lat: f64, lon: f64, system_code: &str, ) -> Option<(f64, f64)>

座標をJGD2011に変換 (水平座標のみ)

Arguments

• lat - 緯度 (度, 10進数)
• lon - 経度 (度, 10進数)
• system_code - 測地系コード (“00”=Tokyo, “01”=JGD2000, “02”=JGD2011)

Returns

変換後の (緯度, 経度) または変換失敗時にNone

pub fn transform_full( &self, lat: f64, lon: f64, alt: f64, system_code: &str, ) -> Result<(f64, f64, f64), TransformError>

座標をJGD2011に変換 (水平+垂直)

Arguments

• lat - 緯度 (度, 10進数)
• lon - 経度 (度, 10進数)
• alt - 標高 (m)
• system_code - 測地系コード (“00”=Tokyo, “01”=JGD2000, “02”=JGD2011)

Returns

変換後の (緯度, 経度, 標高)

pub fn to_wgs84( &self, lat: f64, lon: f64, system_code: &str, ) -> Option<(f64, f64)>

座標をWGS84に変換 (水平座標のみ)

Tokyo Datum/JGD2000 → JGD2011 → WGS84 の変換を行う。

Arguments

• lat - 緯度 (度, 10進数)
• lon - 経度 (度, 10進数)
• system_code - 測地系コード (“00”=Tokyo, “01”=JGD2000, “02”=JGD2011)

Returns

変換後の (経度, 緯度) - WGS84座標