Skip to Content

traitement

  1. """coordonnées du premier sommet"""
  2. debligne =grass.read_command("v.to.db", flags="p", map="testgrass", type="line", option="start", units="meters" , quiet=True)
  3. debx=float(debligne.split("|")[1])
  4. deby=float(debligne.split("|")[2])
  5. debz=float(debligne.split("|")[3])
  6.  
  7. """coordonnées du deuxième sommet"""
  8. finligne=grass.read_command("v.to.db", flags="p", map="testgrass",type="line", option="end", units="meters", quiet=True)
  9. finx=float(finligne.split("|")[1])
  10. finy=float(finligne.split("|")[2])
  11. finz=float(finligne.split("|")[3])
  12.  
  13. """création d'une couche vide testpt"""
  14. grass.read_command("v.edit",tool="create",map="testpt")
  15.  
  16. """création d'un fichier temporaire au format standard:
  17. P  1 1
  18. 186139.12370417 53082.65419398
  19. 1 1        
  20. P  1 1
  21. 188199.1227985 53467.75855873
  22. 1 2"""
  23. fichtemp = grass.tempfile()
  24. pfich = open(fichtemp, 'w')
  25. print>> pfich, "P", "1", "1"
  26. print>> pfich, debx, deby, debz
  27. print>> pfich, "1", 1
  28. print>> pfich, "P", "1", "1"
  29. print>> pfich, finx, finy, finz
  30. print>> pfich, "1", 2
  31. pfich.close()
  32.  
  33. """importation des points dans la couche testpt"""
  34. grass.run_command("v.edit", input=fichtemp, flags="n", tool="add",map="testpt", snap="node", thresh="1")
  • 3394 lectures


Theme by Dr. Radut.