r3.cross.rast.html 2.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475
  1. <h2>DESCRIPTION</h2>
  2. This module creates a cross section 2D map from one 3D raster volume
  3. map based on a 2D elevation map. It checks if the value of the
  4. elevation map is located in the z-coordinate space of the 3D map. If
  5. so, the 3D voxel value for this position is transferred to the related
  6. cross section output map cell, otherwise the NULL value is set.
  7. <center>
  8. <img src="r3.cross.rast.png" border="0"><br>
  9. <table border="0" width="700">
  10. <tr><td><center>
  11. <i>How r3.cross.rast works</i>
  12. </center></td></tr>
  13. </table>
  14. </center>
  15. <p>If the 2D and 3D region settings are different,
  16. the 2D resolution will be adjust to the 3D resolution.
  17. <h2>NOTES</h2>
  18. To create a cut plane elevation map use <em>r.mapcalc</em>. Some examples:
  19. <ul>
  20. <li>To create a cut plane elevation map in x direction type
  21. <br><i>r.mapcalc "cutplane = col()*x"</i>,<br> x be the value for
  22. the elevation. If the range of col() is 1 ... 10, the elevation map
  23. has the range 1 ... 10 if x == 1 and if x == 10 the range 10
  24. ... 100</li>
  25. <li>To create a cut plane elevation map in y direction type
  26. <br><i>r.mapcalc "cutplane = row()*x"</i>,<br> x be the value for
  27. the elevation. If the range of col() is 1 ... 10, the elevation map
  28. has the range 1 ... 10 if x == 1 and if x == 10 the range 10
  29. ... 100</li>
  30. <li>The user can also make a cut in y and x direction with <em>r.mapcalc</em> by
  31. using <br><i>r.mapcalc "cutplane = (row()+col())*x"</i></li>
  32. </ul>
  33. <h2>EXAMPLE</h2>
  34. <h3>Simple Spearfish example</h3>
  35. <div class="code"><pre>
  36. g.region -d
  37. g.region res=150 res3=150 t=1000 b=0 tbres=100
  38. # synthetic data, could be geological structures:
  39. r3.mapcalc "map3d = sin(row())+sin(col())+sin(depth()*depth())"
  40. #create a cutplane map
  41. r.mapcalc "cutplane = col()*10"
  42. #create the cross section map
  43. r3.cross.rast input=map3d elevation=cutplane output=crosssection
  44. </pre></div>
  45. <h2>SEE ALSO</h2>
  46. <em>
  47. <a href="g.region.html">g.region</a>,
  48. <a href="r.mapcalc.html">r.mapcalc</a>,
  49. <a href="r3.mapcalc.html">r3.mapcalc</a>,
  50. <a href="r3.to.rast.html">r3.to.rast</a>
  51. </em>
  52. <h2>AUTHOR</h2>
  53. S&ouml;ren Gebbert
  54. <!--
  55. <p>
  56. <i>Last changed: $Date$</i>
  57. -->