Tutorial on ParaView
Purpose
This is a tutorial that will demonstrate how to use ParaView to examine AIM files, and how to use vtkski and python scripting to construct finite element models with complicated coundary conditions.
This tutorial was done on OS-X, but should be essentially similar for Windows or Linux, with some minor differences.
Prerequisites
You'll need ParaView built with the bonelab extensions. You should be able to find in on the server mike . On OS-X, click on this link: afp://mike.kin.ucalgary.ca/Data0/Software/ParaView_with_extensions/ . On Windows, enter \\mike\Data0\Software\ParaView_with_extensions in an Explorer address bar. In either case, you'll need to be on the university network.
You'll also need python. Windows users can get it from [[1]] ; everyone else has it already.
You will need a compiled version of the vtkski library. To compile it yourself, refer to Bone lab libraries#Build_Instructions .
You can find the data file used in this tutorial, SAWBONE51_BC_CSTR.AIM, as well as the python script that we will construct, on the server mike at /Data0/Tutorials/ParaView_and_vtkski .
Background Material
- If you're not familiar with Python, you might want to go through the tutorial [[2]] .
- There are a couple of books on VTK floating about the lab, "The Visualization Toolkit" and "VTK User's Guide".
- The API documentation for VTK is at [[3]] . This is where you go to find out exactly what classes are available in VTK and exactly how to use them.
- The API documentation for vtkski is at [[4]] .
Model Overview
We want to construct a finite element model of an AIM file that contains a scan of a bone with a screw in it. Physical tests were done in which the bone was constrained with a ring and the screw was pulled out. We want to model this with a finite element model.
Examining the AIM file with ParaView
Open ParaView. Click File->Open, and browse to SAWBONE51_BC_CSTR.AIM . Click OK. Our custom version of ParaView with bonelab extensions is able to open AIMs, but it may inquire about the file type, as shown in the following dialog. Just select "AIM File" and proceed.
Click the Apply button in the Object Inspector. You should see the following, which is so far only the outline of the bounds of the data.
In order to convert this data, which is 3D image data to something that can be rendered, a number of different filters can be used. We are going to use the Threshold filter. You can find this at the menu item Filters->Common->Threshold . There is also a button for it on the toolbar. Click it. Now in the Object Inspector, set the Lower Threshold to 1. (If you keep the default value of 0, you just end up rendering a box with the size of the bounds of the data.) Now hit Apply.
File:ParaView Threshold AIM.png
We still don't see much! Why not? Well, this is such a large data set, that ParaView has decided to render just the outline for now, to save time. To see the object, go to the Display tab of the Object Inspector, scroll down to the Style box, and change the Representation to Surface. Voila!
File:ParaView Change Representation.png
You can click and drag on the image to rotate it. You can also right-click and hold down shift and control, with different effects. We are now going to show only the screw, which can be identified as material 127. Make sure the Threshold filter is selected, then select the Properties tab in the Object Inspector, and set both Lower and Upper Threshold to 127, as shown. Also in the Display tab, I have scrolled to the Color box, and clicked on "Set Solid Color" to make the screw blue.
Now we are going to show the bone as well. Make sure that SAWBONES51_BC_CSTR.AIM is selected in the Pipeline browser, then click on the Threshold button to create a second Threshold filter. Change its Lower and Upper Thresholds to 126, then click Apply.
