محول BJ2000 → WGS84

// BJ2000 → WGS84 (7-parameter Helmert transformation)
// Between WGS84 ellipsoid and CGCS2000/GRS80 ellipsoid

// Step 1: Convert from CGCS2000/GRS80 to geocentric Cartesian
// Using source ellipsoid: a = 6378137.0, 1/f = 298.257222101
N = a_source / √(1 - e²_source × sin²(φ))
X = (N + h) × cos(φ) × cos(λ)
Y = (N + h) × cos(φ) × sin(λ)
Z = (N × (1 - e²_source) + h) × sin(φ)

// Step 2: Helmert 7-parameter transform (BJ2000 → WGS84)
// Parameters: dx, dy, dz (m), rx, ry, rz (rad), s (ppm)
// Reverse sign for inverse transformation
X' = ΔX + (1 + s)(X + Rz·Y - Ry·Z)
Y' = ΔY + (1 + s)(-Rz·X + Y + Rx·Z)
Z' = ΔZ + (1 + s)(Ry·X - Rx·Y + Z)

// Step 3: Convert back to geodetic using WGS84 ellipsoid
// a = 6378137.0, 1/f = 298.257223563
p = √(X'² + Y'²)
θ = atan2(Z' × a_target, p × (1 - f_target))
φ' = atan2(Z' + e²_target × (1 - f_target) × a_target × sin³(θ),
         p - e²_target × a_target × cos³(θ))
λ' = atan2(Y', X')

BJ2000 إلى WGS84 يتطلب تحويل هلمرت كاملًا ذي 7 معاملات لأن النظامين يستخدمان مجسمات إهليلجية مختلفة وأصول datum مختلفة. تعتمد دقة التحويل على جودة قيم 7 معاملات الإقليمية المستخدمة.

1. أدخل إحداثيات BJ2000 في حقل الإدخال (خط العرض، خط الطول، زوج واحد لكل سطر)
2. أدخل قيم هلمرت المناسبة ذات 7 معاملات (dx، dy، dz، rx، ry، rz، s) لمنطقتك
3. انقر فوق زر التحويل لتحويل الإحداثيات من BJ2000 إلى WGS84
4. راجع إحداثيات WGS84 المحولة في حقل الإخراج
5. انسخ النتائج أو احفظها كملف XLSX للاستخدام لاحقًا