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

lecture et traitement du worldfile

  1. from __future__ import with_statement #nécessaire pour les versions < 2.6
  2. from PIL import Image
  3. tfw = "montif.tfw"
  4. raster = Image.open("montif.tif")
  5.  
  6. # dimension du raster
  7. rasterx = raster.size[0]
  8. rastery = raster.size[1]
  9.  
  10. #paramètres du fichier tfw
  11. #lecture du fichier
  12. with open(tfw) as f:
  13. data = [line.strip('\n').strip() for line in f]
  14. #paramètres
  15. taillepixelx = float(data[0])
  16. taillepixely = float(data[3])
  17. topcentregauchex = float(data[4])
  18. topcentregauchey = float(data[5])
  19.  
  20. #calculs
  21. xsupgauche = topcentregauchex - (taillepixelx * .5)
  22. ysupgauche = topcentregauchey - (taillepixely * .5)
  23. print xsupgauche,ysupgauche
  24. 162012.677881 108172.869385
  25. #taille du raster en m
  26. lx = rasterx * taillepixelx
  27. print "dimension en x: ", lx
  28. 8002.3097846
  29. ly = rastery * taillepixely
  30. print "dimension en y: ", ly
  31. -10002.8872308
  32. #calcul des coordonnées du point inférieur droit
  33. xinfdroit = xsupgauche + lx
  34. yinfdroit = ysupgauche + ly
  35. print xinfdroit,yinfdroit
  36. 170014.987666 98169.9821547