محول WGS84 → BJ2000

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

// Step 1: Convert from WGS84 to geocentric Cartesian
// Using source ellipsoid: a = 6378137.0, 1/f = 298.257223563
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 (WGS84 → BJ2000)
// Parameters: dx, dy, dz (m), rx, ry, rz (rad), s (ppm)
// 
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 CGCS2000/GRS80 ellipsoid
// a = 6378137.0, 1/f = 298.257222101
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')

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

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