1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
//! GPS→ENU, Sim3 정렬, 공유 3D 점 정렬, 점군 후처리.
//!
//! - [`geodesy`]: WGS84 LLA ↔ ECEF ↔ ENU.
//! - [`umeyama`](mod@umeyama): Umeyama Sim3, LO-RANSAC 견고 추정, 반복 재가중(견고 Umeyama).
//! - [`model_aligner`]: 등록 영상 투영 중심 ↔ GPS ENU 정렬(강건 Sim3 추정).
//! - [`shared`]: 두 재구성 간 공유 3D 점 대응과 그 정렬.
//! - [`kdtree`], [`cloud`]: 점군 변환·근접 마스킹·추출·스냅샷 합성.
//! - [`reanchor`]: 새 정밀 모델 기준으로 이전 구역 좌표계를 다시 잇는 연쇄 Sim3.
//!
//! # 사용 예
//!
//! ```
//! use cumulus3d_align::{compose_snapshot, umeyama, EnuFrame, SnapshotOptions};
//! use cumulus3d_core::io::PointCloud;
//! use cumulus3d_core::{Quat, Sim3, Vec3};
//!
//! // GPS(WGS84) ↔ ENU: 원점 기준 동쪽·북쪽·위 미터 좌표.
//! let enu = EnuFrame::new(37.5, 127.0, 10.0);
//! let p = enu.lla_to_enu(37.5001, 127.0001, 10.0);
//! let (lat, lon, _alt) = enu.enu_to_lla(&p);
//! assert!((lat - 37.5001).abs() < 1e-9 && (lon - 127.0001).abs() < 1e-9);
//!
//! // 대응점 집합 사이의 Sim3(src → dst) 추정.
//! let truth = Sim3::new(2.0, Quat::from_axis_angle(&Vec3::z(), 0.3), Vec3::new(1.0, 2.0, 3.0));
//! let src: Vec<Vec3> = (0..10).map(|i| Vec3::new(i as f64, (i * i) as f64 * 0.1, (i as f64).sin())).collect();
//! let dst: Vec<Vec3> = src.iter().map(|x| truth.transform_point(x)).collect();
//! let est = umeyama(&src, &dst, true).unwrap();
//! assert!((est.scale - 2.0).abs() < 1e-9);
//!
//! // 스냅샷: 정밀 점과 1.5 m 안에 있는 초벌 점을 지우고 합친 뒤 1/stride 추출.
//! let fine = PointCloud { positions: vec![[0.0, 0.0, 0.0]], ..Default::default() };
//! let coarse = PointCloud { positions: vec![[0.5, 0.0, 0.0], [10.0, 0.0, 0.0]], ..Default::default() };
//! let snap = compose_snapshot(&[&fine], &[&coarse], &SnapshotOptions { mask_radius: 1.5, stride: 1 });
//! assert_eq!(snap.len(), 2);
//! ```
//!
//! ```no_run
//! use cumulus3d_align::{align_to_gps_file, ModelAlignerOptions};
//! let mut rec = cumulus3d_core::interop::read_model("model/0").unwrap();
//! let a = align_to_gps_file(&mut rec, "gps_ref.txt", &ModelAlignerOptions::default()).unwrap();
//! println!("인라이어 {}/{}, 중앙 오차 {:.2} m", a.num_inliers, a.common.len(), a.median_error);
//! ```
pub use ;
pub use ;
pub use ;
pub use KdTree;
pub use ;
pub use AnchorChain;
pub use ;
pub use ;