Provided by: grass-doc_6.4.3-3_all bug

NAME

       r.out.vtk  - Converts raster maps into the VTK-Ascii format

KEYWORDS

       raster, export

SYNOPSIS

       r.out.vtk
       r.out.vtk help
       r.out.vtk    [-pstvoc]    input=name[,name,...]     [output=string]     [elevation=string]
       [null=float]               [elevation2d=float]               [rgbmaps=string[,string,...]]
       [vectormaps=string[,string,...]]      [elevscale=float]      [dp=integer]      [--verbose]
       [--quiet]

   Flags:
       -p
           Create VTK point data instead of VTK cell data (if no elevation map is given)

       -s
           Use structured grid for elevation (not recommended)

       -t
           Use polydata-trianglestrips for elevation grid creation

       -v
           Use polydata-vertices for elevation grid creation (to use with vtkDelauny2D)

       -o
           Scale factor effects the origin (if no elevation map is given)

       -c
           Correct the coordinates to fit the VTK-OpenGL precision

       --verbose
           Verbose module output

       --quiet
           Quiet module output

   Parameters:
       input=name[,name,...]
           Name of input raster map(s)

       output=string
           Name for VTK-ASCII output file

       elevation=string
           Elevation raster map

       null=float
           Value to represent no data cell
           Default: -10.0

       elevation2d=float
           Elevation (if no elevation map is specified)
           Default: 0.0

       rgbmaps=string[,string,...]
           Three (r,g,b) raster maps to create rgb values [redmap,greenmap,bluemap]

       vectormaps=string[,string,...]
           Three (x,y,z) raster maps to create vector values [xmap,ymap,zmap]

       elevscale=float
           Scale factor for elevation
           Default: 1.0

       dp=integer
           Number of significant digits (floating point only)
           Options: 0-20
           Default: 12

DESCRIPTION

       Outputs raster maps in VTK-ASCII format. Map's are  valid  raster  map's  in  the  current
       mapset.   output  is  the  name  of an VTK-ASCII file which will be written in the current
       working directory.  If output is not  specified  then  stdout  is  used.   The  module  is
       sensitive to region settings (set with g.region).

       Elevation,  scaling,  point/celldata, vector and RGB Data are supported.  If the map is in
       LL projection, the elevation values will automatically scaled to degrees. It  is  supposed
       that  the  elevation  values  are  provided  in  meters.  If the elevation values are in a
       different unit than meters, use the scale parameter to convert the units.

       If no elevation map is given, the user can set the hight of the map by one  value.   Point
       or  cell  data  are  available.  Also  scaling is supported for this elevation value.  The
       elevation value must be provided in meters.

       The RGB input requires three raster maps: red, green, blue  - in  this  order.   The  maps
       must have values between 0 and 255, otherwise you will get lots of warnings and the values
       are set to 0.  More than one RGB dataset (3 maps) is not supported.

       The vector input requires three raster maps: x, y, z -- defining the vector coordinates  -
       in this order.  More than one vector dataset (3 maps) is not supported.

NOTES

       This filter generates:

                     structured points with celldata or pointdata if no elevationfile is given

                     structured  grid  (not  recommendet)  with  pointdata if an elevationfile is
                     given

                     polydataset with pointdata if an elevationfile is given (default)
       and puts this in a simple VTK-ASCII file. Nor XML or binary output are  supported.  It  is
       possible  to  choose  more  then  one raster map to be written to the VTK-ASCII file. Each
       cell-/pointdata is named like the raster map it represents.  You can visualize  this  file
       with  the  VTK  Toolkit, Paraview and MayaVi which are based on VTK.  If you have a raster
       map with partly no data, use the threshold filter in paraview to visualize the valid data.
       Just  filter  all data which is greater/lesser than the chosen null value in the VTK-ASCII
       file.
       If elevation map is chosen, a polygonal grid is created  with  quads,  but  the  user  can
       choose also triangle strips or vertices.  These dataformats a documented at VTK Toolkit.

       If  the  "-c"  flag  is  used  and  the data should be visualised together with other data
       exported via *.out.vtk modules, be sure the "-c" flag was also set in these modules.   But
       this  will  only  work  with  data  from  the  SAME  location (The reference point for the
       coordinates transformation is based on the center point of the default region).

   Difference between point- and celldata
       r.out.vtk can export raster cells with different representations.

                     pointdata -- the cells/values are represented by the  center  of  the  cell.
                     Instead  of cells, points are created. Each point can hold different values,
                     but the user can only visualize one value at a time.  These  points  can  be
                     connected in different ways.

                     celldata  --  is  only provided if no elevation map is given.  The cells are
                     created with the same hight and width  as  in  GRASS.  Each  cell  can  hold
                     different values, but the user can only visualize one value at a time.

EXAMPLE

   Simple Spearfish example

       # set region
       g.region n=4926970 s=4914857 w=591583 e=607793 res=50 -p
       # export the data
       r.out.vtk input=elevation.10m,slope,aspect elevation=elevation.10m output=/tmp/out.vtk
       # visualize in Paraview or other VTK viewer:
       paraview --data=/tmp/out.vtk

   Spearfish example with RGB data

       #set the region
       g.region n=4926990 s=4914840 w=591570 e=607800 res=30 -p
       # using r.in.wms to create RGB data to get a satellite coverage
       r.in.wms layers=global_mosaic mapserver=http://wms.jpl.nasa.gov/wms.cgi \
                output=wms_global_mosaic
       # export the data to VTK
       r.out.vtk rgbmaps=wms_global_mosaic.red,wms_global_mosaic.green,wms_global_mosaic.blue \
                 elevation=elevation.10m output=/tmp/out.vtk
       # visualize in Paraview or other VTK viewer:
       paraview --data=/tmp/out.vtk
        Paraview RGB visualization notes
       To achieve proper RGB overlay:

                      In Paraview, click "Apply"

                      Select the "Display" tab and choose "Color by" to switch from input scalars
                     to rgb scalars

                      Disable the "Map Scalars" check button in the display tab to avoid the  use
                     of a lookup table

SEE ALSO

        r3.out.vtk, r.out.ascii, g.region
       GRASS and Paraview Wiki page

AUTHOR

       Soeren Gebbert

       Last changed: $Date: 2013-02-15 14:04:18 -0800 (Fri, 15 Feb 2013) $

       Full index

       © 2003-2013 GRASS Development Team