محول WGS84 → ED50

// WGS84 → ED50 (7-parameter Helmert transformation)
// Between WGS84 ellipsoid and Hayford 1909 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 → ED50)
// 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 Hayford 1909 ellipsoid
// a = 6378388.0, 1/f = 297.0
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 إلى ED50 يتطلب تحويل هلمرت كاملًا ذي 7 معاملات لأن النظامين يستخدمان مجسمات إهليلجية مختلفة وأصول datum مختلفة. تعتمد دقة التحويل على جودة قيم 7 معاملات الإقليمية المستخدمة.

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