# RGrid

VTKExamples/Cxx/RectilinearGrid/RGrid

### Description

Creating a rectilinear grid dataset. The coordinates along each axis are defined using an instance of vtkDataArray.

Info

See Figure 5-20 in Chapter 5 the VTK Textbook.

Other Languages

See (Python)

### Code

RGrid.cxx

    // This example shows how to create a rectilinear grid.
    //
    #include <vtkActor.h>
    #include <vtkCamera.h>
    #include <vtkDoubleArray.h>
    #include <vtkNamedColors.h>
    #include <vtkPolyDataMapper.h>
    #include <vtkProperty.h>
    #include <vtkRectilinearGrid.h>
    #include <vtkRectilinearGridGeometryFilter.h>
    #include <vtkRenderWindow.h>
    #include <vtkRenderWindowInteractor.h>
    #include <vtkRenderer.h>
    #include <vtkSmartPointer.h>

    int main(int, char *[])
    {
      int i;
      static double x[47]={
        -1.22396, -1.17188, -1.11979, -1.06771, -1.01562, -0.963542,
        -0.911458, -0.859375, -0.807292, -0.755208, -0.703125, -0.651042,
        -0.598958, -0.546875, -0.494792, -0.442708, -0.390625, -0.338542,
        -0.286458, -0.234375, -0.182292, -0.130209, -0.078125, -0.026042,
        0.0260415, 0.078125, 0.130208, 0.182291, 0.234375, 0.286458,
        0.338542, 0.390625, 0.442708, 0.494792, 0.546875, 0.598958,
        0.651042, 0.703125, 0.755208, 0.807292, 0.859375, 0.911458,
        0.963542, 1.01562, 1.06771, 1.11979, 1.17188};
      static double y[33]={-1.25, -1.17188, -1.09375, -1.01562, -0.9375, -0.859375,
                           -0.78125, -0.703125, -0.625, -0.546875, -0.46875, -0.390625,
                           -0.3125, -0.234375, -0.15625, -0.078125, 0, 0.078125,
                           0.15625, 0.234375, 0.3125, 0.390625, 0.46875, 0.546875,
                           0.625, 0.703125, 0.78125, 0.859375, 0.9375, 1.01562,
                           1.09375, 1.17188, 1.25};
      static double z[44]={0, 0.1, 0.2, 0.3, 0.4, 0.5,
                           0.6, 0.7, 0.75, 0.8, 0.9, 1,
                           1.1, 1.2, 1.3, 1.4, 1.5, 1.6,
                           1.7, 1.75, 1.8, 1.9, 2, 2.1,
                           2.2, 2.3, 2.4, 2.5, 2.6, 2.7,
                           2.75, 2.8, 2.9, 3, 3.1, 3.2,
                           3.3, 3.4, 3.5, 3.6, 3.7, 3.75,
                           3.8, 3.9};

      // Create a rectilinear grid by defining three arrays specifying the
      // coordinates in the x-y-z directions.
      vtkSmartPointer<vtkDoubleArray> xCoords =
        vtkSmartPointer<vtkDoubleArray>::New();
      for (i=0; i<47; i++)
      {
        xCoords->InsertNextValue(x[i]);
      }
      vtkSmartPointer<vtkDoubleArray> yCoords =
        vtkSmartPointer<vtkDoubleArray>::New();
      for (i=0; i<33; i++)
      {
        yCoords->InsertNextValue(y[i]);
      }
      vtkSmartPointer<vtkDoubleArray> zCoords =
        vtkSmartPointer<vtkDoubleArray>::New();
      for (i=0; i<44; i++)
      {
        zCoords->InsertNextValue(z[i]);
      }
      // The coordinates are assigned to the rectilinear grid. Make sure that
      // the number of values in each of the XCoordinates, YCoordinates,
      // and ZCoordinates is equal to what is defined in SetDimensions().
      //
      vtkSmartPointer<vtkRectilinearGrid> rgrid =
        vtkSmartPointer<vtkRectilinearGrid>::New();
      rgrid->SetDimensions(47,33,44);
      rgrid->SetXCoordinates(xCoords);
      rgrid->SetYCoordinates(yCoords);
      rgrid->SetZCoordinates(zCoords);

      // Extract a plane from the grid to see what we've got.
      vtkSmartPointer<vtkRectilinearGridGeometryFilter> plane =
        vtkSmartPointer<vtkRectilinearGridGeometryFilter>::New();
      plane->SetInputData(rgrid);
      plane->SetExtent(0,46, 16,16, 0,43);

      vtkSmartPointer<vtkPolyDataMapper> rgridMapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
      rgridMapper->SetInputConnection(plane->GetOutputPort());

      vtkSmartPointer<vtkNamedColors> colors =
        vtkSmartPointer<vtkNamedColors>::New();

      vtkSmartPointer<vtkActor> wireActor =
        vtkSmartPointer<vtkActor>::New();
      wireActor->SetMapper(rgridMapper);
      wireActor->GetProperty()->SetColor(colors->GetColor3d("Banana").GetData());
      wireActor->GetProperty()->EdgeVisibilityOn();

      // Create the usual rendering stuff.
      vtkSmartPointer<vtkRenderer> renderer =
        vtkSmartPointer<vtkRenderer>::New();
      vtkSmartPointer<vtkRenderWindow> renWin =
        vtkSmartPointer<vtkRenderWindow>::New();
      renWin->AddRenderer(renderer);
      vtkSmartPointer<vtkRenderWindowInteractor> iren =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
      iren->SetRenderWindow(renWin);

      renderer->AddActor(wireActor);
      renderer->SetBackground(colors->GetColor3d("Beige").GetData());
      renderer->ResetCamera();
      renderer->GetActiveCamera()->Elevation(60.0);
      renderer->GetActiveCamera()->Azimuth(30.0);
      renderer->GetActiveCamera()->Zoom(1.0);

      renWin->SetSize(640, 480);

      // interact with data
      renWin->Render();
      iren->Start();

      return EXIT_SUCCESS;
    }

### CMakeLists.txt

    cmake_minimum_required(VERSION 2.8)

    PROJECT(RGrid)

    find_package(VTK REQUIRED)
    include(${VTK_USE_FILE})

    add_executable(RGrid MACOSX_BUNDLE RGrid.cxx )

    target_link_libraries(RGrid ${VTK_LIBRARIES})

### Download and Build RGrid

Click here to download RGrid and its CMakeLists.txt file. Once the tarball RGrid.tar has been downloaded and extracted,

    cd RGrid/build

If VTK is installed:

    cmake ..

If VTK is not installed but compiled on your system, you will need to specify the path to your VTK build:

    cmake -DVTK_DIR:PATH=/home/me/vtk_build ..

Build the project:

    make

and run it:

    ./RGrid

WINDOWS USERS

Be sure to add the VTK bin directory to your path. This will resolve the VTK dll's at run time.