### Description¶

Clip geometry using a 2D texture map and two implicit functions. The technique is described in : Geometric Clipping Using Boolean Textures.

Other languages

See (Python)

Question

### Code¶

#include <vtkActor.h>
#include <vtkBooleanTexture.h>
#include <vtkCamera.h>
#include <vtkDataSetMapper.h>
#include <vtkImplicitTextureCoords.h>
#include <vtkNamedColors.h>
#include <vtkNew.h>
#include <vtkProperty.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkSmartPointer.h>
#include <vtkSphereSource.h>
#include <vtkTexture.h>

namespace {

double positions[][3]{{-4, 4, 0},  {-2, 4, 0},  {0, 4, 0},  {2, 4, 0},
{-4, 2, 0},  {-2, 2, 0},  {0, 2, 0},  {2, 2, 0},
{-4, 0, 0},  {-2, 0, 0},  {0, 0, 0},  {2, 0, 0},
{-4, -2, 0}, {-2, -2, 0}, {0, -2, 0}, {2, -2, 0}};

unsigned char solid[]{255, 255};
unsigned char clear[]{255, 0};
unsigned char edge[]{0, 255};

vtkSmartPointer<vtkBooleanTexture> MakeBooleanTexture(int, int, int);
} // namespace

int main(int /* argc */, char* /* argv */[])
{
vtkNew<vtkNamedColors> colors;

vtkNew<vtkRenderWindow> renWin;

vtkNew<vtkRenderWindowInteractor> iren;
iren->SetRenderWindow(renWin);

vtkNew<vtkRenderer> aren;

// define two elliptical cylinders
quadric1->SetCoefficients(1, 2, 0, 0, 0, 0, 0, 0, 0, -.07);

quadric2->SetCoefficients(2, 1, 0, 0, 0, 0, 0, 0, 0, -.07);

// create a sphere for all to use
vtkNew<vtkSphereSource> aSphere;
aSphere->SetPhiResolution(21);
aSphere->SetThetaResolution(21);

// create texture coordinates for all
vtkNew<vtkImplicitTextureCoords> tcoords;
tcoords->SetInputConnection(aSphere->GetOutputPort());

vtkNew<vtkDataSetMapper> aMapper;
aMapper->SetInputConnection(tcoords->GetOutputPort());

// create a mapper, sphere and texture map for each case
for (int i = 0; i < 16; i++)
{
auto aBoolean = MakeBooleanTexture(i, 64, 0);

vtkNew<vtkTexture> aTexture2;
aTexture2->SetInputConnection(aBoolean->GetOutputPort());
aTexture2->InterpolateOff();
aTexture2->RepeatOff();

vtkNew<vtkActor> anActor2;

anActor2->SetMapper(aMapper);
anActor2->SetTexture(aTexture2);
anActor2->SetPosition(&positions[i][0]);
anActor2->SetScale(2.0, 2.0, 2.0);
anActor2->GetProperty()->SetColor(
colors->GetColor3d("MistyRose").GetData());

}

aren->SetBackground(colors->GetColor3d("SlateGray").GetData());
renWin->SetSize(500, 500);

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

iren->Start();

return EXIT_SUCCESS;
}

namespace {

vtkSmartPointer<vtkBooleanTexture>
MakeBooleanTexture(int caseNumber, int resolution, int thickness)
{
vtkNew<vtkBooleanTexture> booleanTexture;

booleanTexture->SetXSize(resolution);
booleanTexture->SetYSize(resolution);
booleanTexture->SetThickness(thickness);

switch (caseNumber)
{
case 0:
booleanTexture->SetInIn(solid);
booleanTexture->SetOutIn(solid);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(solid);
booleanTexture->SetOnIn(solid);
booleanTexture->SetOnOut(solid);
booleanTexture->SetInOn(solid);
booleanTexture->SetOutOn(solid);
break;
case 1: /* cut inside 1 */
booleanTexture->SetInIn(clear);
booleanTexture->SetOutIn(solid);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(solid);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(solid);
break;
case 2: /* cut outside 1 */
booleanTexture->SetInIn(solid);
booleanTexture->SetOutIn(clear);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(solid);
booleanTexture->SetInOn(solid);
booleanTexture->SetOutOn(edge);
break;
case 3: /* cut all 1 */
booleanTexture->SetInIn(clear);
booleanTexture->SetOutIn(clear);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(clear);
booleanTexture->SetOnOut(solid);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(edge);
break;
case 4:
booleanTexture->SetInIn(solid);
booleanTexture->SetOutIn(solid);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(solid);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(solid);
break;
case 5:
booleanTexture->SetInIn(clear);
booleanTexture->SetOutIn(solid);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(clear);
booleanTexture->SetOutOn(solid);
break;
case 6:
booleanTexture->SetInIn(solid);
booleanTexture->SetOutIn(clear);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(edge);
break;
case 7:
booleanTexture->SetInIn(clear);
booleanTexture->SetOutIn(clear);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutOut(solid);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(clear);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(clear);
booleanTexture->SetOutOn(edge);
break;
case 8:
booleanTexture->SetInIn(solid);
booleanTexture->SetOutIn(solid);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(solid);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(solid);
booleanTexture->SetOutOn(edge);
break;
case 9:
booleanTexture->SetInIn(clear);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutIn(solid);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(edge);
break;
case 10:
booleanTexture->SetInIn(solid);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutIn(clear);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(solid);
booleanTexture->SetOutOn(clear);
break;
case 11:
booleanTexture->SetInIn(clear);
booleanTexture->SetInOut(solid);
booleanTexture->SetOutIn(clear);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(clear);
booleanTexture->SetOnOut(edge);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(clear);
break;
case 12:
booleanTexture->SetInIn(solid);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutIn(solid);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(solid);
booleanTexture->SetOnOut(clear);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(edge);
break;
case 13:
booleanTexture->SetInIn(clear);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutIn(solid);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(clear);
booleanTexture->SetInOn(clear);
booleanTexture->SetOutOn(edge);
break;
case 14:
booleanTexture->SetInIn(solid);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutIn(clear);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(edge);
booleanTexture->SetOnIn(edge);
booleanTexture->SetOnOut(clear);
booleanTexture->SetInOn(edge);
booleanTexture->SetOutOn(clear);
break;
case 15:
booleanTexture->SetInIn(clear);
booleanTexture->SetInOut(clear);
booleanTexture->SetOutIn(clear);
booleanTexture->SetOutOut(clear);
booleanTexture->SetOnOn(clear);
booleanTexture->SetOnIn(clear);
booleanTexture->SetOnOut(clear);
booleanTexture->SetInOn(clear);
booleanTexture->SetOutOn(clear);
break;
}

return booleanTexture;
}
} // namespace


### CMakeLists.txt¶

cmake_minimum_required(VERSION 3.12 FATAL_ERROR)

find_package(VTK COMPONENTS
CommonColor
CommonCore
CommonDataModel
FiltersTexture
ImagingHybrid
InteractionStyle
RenderingContextOpenGL2
RenderingCore
RenderingFreeType
RenderingGL2PSOpenGL2
RenderingOpenGL2
)

if (NOT VTK_FOUND)
message(FATAL_ERROR "TextureCutQuadric: Unable to find the VTK build folder.")
endif()

# Prevent a "command line is too long" failure in Windows.
set(CMAKE_NINJA_FORCE_RESPONSE_FILE "ON" CACHE BOOL "Force Ninja to use response files.")
target_link_libraries(TextureCutQuadric PRIVATE ${VTK_LIBRARIES} ) # vtk_module_autoinit is needed vtk_module_autoinit( TARGETS TextureCutQuadric MODULES${VTK_LIBRARIES}
)


cd TextureCutQuadric/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:

./TextureCutQuadric


WINDOWS USERS

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