Skip to Content

Nous adressons toutes nos pensées à la famille de notre ami Jérôme !

http://www.forumsig.org/showthread.php/43488-Disparition-de-Phoenix

calcul de l'indicateur d'écart RMSE

  1. import numpy as np
  2. # première série de points
  3. x1 = [194184.81,194286.47,194289.02,195292.21,196502.49]
  4. x2 = [194186.12,194283.64,194288.34,195294.12,196501.84]
  5. # deuxième série de points
  6. y1 = [469054.43,469550.98,469945.65,470539.12,470941.38]
  7. y2 = [469053.24,469549.89,469943.46,470542.21,470942.73]
  8. # calcul de des erreurs types (RMSE) sur x et y
  9. RMSEx = np.sqrt(np.sum([(x1[i] - x2[i])**2 for i in range(len(x1))])/float(len(x1)))
  10. RMSEy = np.sqrt(np.sum([(y1[i] - y2[i])**2 for i in range(len(y1))])/float(len(y1)))
  11. #calcul du RMSE planaire
  12. RMSEplanaire = np.sqrt((a1*a1)+(b1*b1))
  13.  
  14. # résultat avec 3 chiffres significatifs
  15. print("%.2f" % RMSEplanaire)
  16. 2.57