محول ED50 → WGS84

// ED50 → WGS84 (7-parameter Helmert transformation)
// Between WGS84 ellipsoid and Hayford 1909 ellipsoid

// Step 1: Convert from Hayford 1909 to geocentric Cartesian
// Using source ellipsoid: a = 6378388.0, 1/f = 297.0
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 (ED50 → 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')

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

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