#include <itkParallelSparseFieldLevelSetImageFilter.h>
Inheritance diagram for itk::ParallelSparseFieldLevelSetImageFilter< TInputImage, TOutputImage >:
3. Starting with the first layers adjacent to the active set and moving outwards, reconstruct the distance transform by setting values in the layers according to their neighbors. At the very outer layers, add or remove indicies which have come into or moved out of the sparse field.
Definition at line 247 of file itkParallelSparseFieldLevelSetImageFilter.h.
|
Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 256 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The value type of the time step. This is distinct from PixelType because PixelType may often be a vector value, while the TimeStep is a scalar value. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 260 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 271 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Information derived from the image types. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 269 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
A type for a list of LayerPointerTypes Definition at line 291 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Memory pre-allocator used to manage layer nodes in a multithreaded environment. Definition at line 302 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Node type used in parallel sparse field layer lists. Definition at line 284 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 288 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
A list type used in the algorithm. Definition at line 287 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 304 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Some convenient typedefs. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 270 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The pixel type of the output image will be used in computations. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 275 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 255 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Standard class typedefs Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 253 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The type of the image used to index status information. Necessary for the internals of the algorithm. Definition at line 298 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Type used for storing status information Definition at line 294 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 254 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 277 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Typedefs from the superclass Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 259 of file itkParallelSparseFieldLevelSetImageFilter.h. Referenced by itk::ParallelSparseFieldLevelSetImageFilter< TInputImage, TOutputImage >::ApplyUpdate(). |
|
The data type used in numerical computations. Derived from the output image type. Definition at line 281 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
|
|
Definition at line 330 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Reserves memory in the update buffer Implements itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 414 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
This method is not implemented or necessary for this solver Implements itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 526 of file itkParallelSparseFieldLevelSetImageFilter.h. References itk::ParallelSparseFieldLevelSetImageFilter< TInputImage, TOutputImage >::TimeStepType. |
|
This method is not implemented or necessary for this solver Implements itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >. Definition at line 533 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Check if the load is fairly balanced among the threads. This is performed by just one thread while all other threads wait. This need NOT be performed every iteration because the level-set surface moves slowly and it is correct to believe that during an iteration the movement is small enough that the small gain obtained by load balancing (if any) does not warrant the overhead for calling this method. How often this is done is controlled by a parameter LOAD_BALANCE_ITERATION_FREQUENCY which is defined in the IterateThreaderCallback() function. A parameter that defines a degree of unbalancedness of the load among threads is MAX_PIXEL_DIFFERENCE_PERCENT which is defined in CheckLoadBalance(). |
|
Delete all nodes in a thread's own lists which its used to transfer nodes to neighboring threads. |
|
Delete all nodes in the List |
|
This performs the initial load distribution among the threads. Every thread gets a slab of the data to work on. The slabs created along a specific dimension. Load balancing is performed along the greatest numbered dimension (i.e. the 3rd dimension in the 3D case and the 2nd dimension in the 2D case). During the initializing of the sparse field layer an histogram is computed that stores the number of nodes in the active set for each index along the chosen dimension. This histogram is used to divide the work "equally" among threads so that each thread approximately get the same number of nodes to process. |
|
Constructs the active layer and initialize the first layers inside and outside of the active layer. The active layer defines the position of the zero level set by its values, which are constrained within a range around zero. |
|
Initializes a layer of the sparse field using a previously initialized layer. Builds the list of nodes in m_Layer[to] using m_Layer[from]. Marks values in the m_StatusImage. |
|
Copies the input to the output image. Processing occurs on the output image, so the data type of the output image determines the precision of the calculations (i.e. double or float). This method overrides the parent class method to do some additional processing. Implements itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >.
|
|
Make a copy of the nodes given to one thread by its neighbors to process and insert them into the thread's own list. |
|
Make a copy of the nodes in the FromList and insert them into the ToList. |
|
Delete the data and synchronization primitives used by the threads during iterations. |
|
Reimplement the GenerateData() function from FiniteDifferenceImageFilter for more effective multithreading Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >.
|
|
Definition at line 318 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Run-time type information (and related methods). Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >.
|
|
Set/Get the value of the isosurface to use in the input image. |
|
Set/Get the number of layers to use in the sparse field. Argument is the number of layers on ONE side of the active layer, so the total layers in the sparse field is 2 * NumberOfLayers + 1 |
|
Find the thread to which a pixel belongs |
|
Obtain a thread's region split as per the load balancing is done. |
|
Split the volume uniformly along the chosen dimension for post processing the output. |
|
Constructs the sparse field layers and initializes their values. Also creates data structures that are NOT local to a thread. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >.
|
|
Initializes the values of the active layer set. |
|
This method pre-processes pixels inside and outside the sparse field layers. The default is to set them to positive and negative values, respectively. This is not necessary as part of the calculations, but produces a more intuitive output for the user. |
|
This method calculates the change and does the update, i.e. one iteration of this iterative solver. A barrier class is used to synchronize execution and keep the CalculateChange and ApplyUpdate sections from executing simultaneously. |
|
This method calculates the change and does the update, i.e. one iteration of this iterative solver. A barrier class is used to synchronize execution and keep the CalculateChange and ApplyUpdate sections from executing simultaneously. |
|
|
|
Method for creation through the object factory. Reimplemented from itk::Object.
|
|
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes. Reimplemented from itk::FiniteDifferenceImageFilter< TInputImage, TOutputImage >.
|
|
|
|
Adjusts the values associated with all the index layers of the sparse field by propagating out one layer at a time from the active set. This method also takes care of deleting nodes from the layers which have been marked in the status image as having been moved to other layers. |
|
Adjusts the values in a single layer "to" using values in a neighboring layer "from". The list of indicies in "to" are traversed and assigned new values appropriately. Any indicies in "to" without neighbors in "from" are moved into the "promote" layer (or deleted if "promote" is greater than the number of layers). "InOrOut" == 1 indicates this propagation is inwards (more negative). "InOrOut" == 0 indicates this propagation is outwards (more positive). |
|
Set/Get the value of the isosurface to use in the input image. |
|
Set/Get the number of layers to use in the sparse field. Argument is the number of layers on ONE side of the active layer, so the total layers in the sparse field is 2 * NumberOfLayers + 1 |
|
Thread synchronization methods. |
|
Thread synchronization methods. |
|
Each thread allocates and initializes the data it will use by itself. This maintains the memory locality of data w.r.t. the thread in a shared memory environment. |
|
Does the actual work of updating the output from the UpdateContainer over an output region supplied by the multithreading mechanism. |
|
This method does the actual work of calculating change over a region supplied by the multithreading mechanism. |
|
This method allows a subclass to override the way in which updates to output values are applied during each iteration. The default simply follows the standard finite difference scheme of scaling the change by the timestep and adding to the value of the previous iteration. Definition at line 509 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Each thread allocates and initializes the data it will use by itself. This maintains the memory locality of data w.r.t. the thread in a shared memory environment. |
|
If child classes need an entry point to the start of every iteration step they can override this method. This method is defined but empty in this class. |
|
Redistribute an load among the threads to obtain a more balanced load distribution. This is performed in parallel by all the threads. |
|
Assign background pixels INSIDE the sparse field layers to a new level set with value less than the innermost layer. Assign background pixels OUTSIDE the sparse field layers to a new level set with value greater than the outermost layer. |
|
Performs two tasks. The situation here is that ThreadedProcessStatusList has been called just once after the active layer values have been updated and the UpLists and DownLists formed. Some nodes are now moving into the active layer. The two tasks performed are: 1. modify the status-image like it is performed by the ThreadedProcessStatusList. 2. Update the values in the output-image for those nodes that are moving IN the active layer. |
|
Push each index in the input list into its appropriate status layer (ChangeToStatus) and ... ... update the status image value at that index |
|
|
|
Push each index in the input list into its appropriate status layer (ChangeToStatus) and update the status image value at that index. Also examine the neighbors of the index, (with status SearchForStatus) to determine which need to go onto the output list. |
|
|
|
1. Updates the values (in the output-image) of the nodes in the active layer that are moving OUT of the active layer. These values are used in the ThreadedProcessFirstLayerStatusLists() method to assign values for new nodes that are moving IN the active layer. 2. This function also constructs the up/down lists for nodes that are moving out of the active layer. |
|
Thread synchronization methods. |
|
Thread synchronization methods. |
|
A global barrier used for synchronization between all threads. Definition at line 672 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The boundaries defining thread regions Definition at line 659 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
A boolean variable stating if the boundaries had been changed during CheckLoadBalance() Definition at line 656 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The constant gradient to maintain between isosurfaces in the spare-field of the level-set image. This value defaults to 1.0 Definition at line 339 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
An array storing the individual (local) data structures for each thread. Definition at line 719 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Histogram of number of pixels in each Z plane for the entire 3D volume Definition at line 662 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
This flag tells the solver whether or not to interpolate for the actual surface location when calculating change at each active layer node. By default this is turned on. Subclasses which do not sample propagation (speed), advection, or curvature terms should turn this flag off. Definition at line 729 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The value in the input which represents the isosurface of interest. Definition at line 396 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Storage for layer node objects. Definition at line 393 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
An array which contains all of the layers needed in the sparse field. Layers are organized as follows: m_Layer[0] = active layer, m_Layer[i:odd] = inside layer (i+1)/2, m_Layer[i:even] = outside layer i/2. This is used only during the initialization of the level set. Definition at line 377 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The mapping from a z-value to the thread in whose region the z-value lies Definition at line 665 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Connectivity information for examining neighbor pixels. Definition at line 335 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The number of layers to use in the sparse field. Sparse field will consist of m_NumberOfLayers layers on both sides of a single active layer. This active layer is the interface of interest, i.e. the zero level set. Definition at line 382 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The number of threads to use. Definition at line 646 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 386 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Definition at line 390 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
This image is a copy of the input with m_IsoSurfaceValue subtracted from each pixel. This way we only need to consider the zero level set in our calculations. Makes the implementation easier and more efficient. This is used only during the initialization of the level set. Definition at line 371 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The dimension along which to distribute the load. Definition at line 649 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Special status value which indicates a pending change to a more negative sparse field. Definition at line 353 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Special status value which indicates a pending change to a more positive sparse field. Definition at line 349 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Special status value which indicates a pixel is on the boundary of the image Definition at line 365 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Special status value which indicates pending change to another sparse field layer. Definition at line 361 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
An image of status values used internally by the algorithm. Definition at line 385 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Images used temporarily during the initialization of the thread data structures. Definition at line 389 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Special status value used as a default for indicies which have no meaningful status. Definition at line 357 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Used to check if there are too few pixels remaining. If yes, then we can stop iterating. Definition at line 723 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Multiplicative identity of the ValueType. Definition at line 342 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Additive identity of the ValueType. Definition at line 345 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
Cumulative frequency of number of pixels in each Z plane for the entire 3D volume Definition at line 669 of file itkParallelSparseFieldLevelSetImageFilter.h. |
|
The length of the dimension along which to distribute the load. Definition at line 652 of file itkParallelSparseFieldLevelSetImageFilter.h. |