Nested or intersecting surfaces are proven techniques for visualizing shape differences between static 3D objects (Weigle and Taylor II, IEEE Visualization, Proceedings, pp. 503-510, 2005). In this... Show moreNested or intersecting surfaces are proven techniques for visualizing shape differences between static 3D objects (Weigle and Taylor II, IEEE Visualization, Proceedings, pp. 503-510, 2005). In this paper we present an image-based formulation for these techniques that extends their use to dynamic scenarios, in which surfaces can be manipulated or even deformed interactively. The formulation is based on our new layered rendering pipeline, a generic image-based approach for rendering nested surfaces based on depth peeling and deferred shading. We use layered rendering to enhance the intersecting surfaces visualization. In addition to enabling interactive performance, our enhancements address several limitations of the original technique. Contours remove ambiguity regarding the shape of intersections. Local distances between the surfaces can be visualized at any point using either depth fogging or distance fields: Depth fogging is used as a cue for the distance between two surfaces in the viewing direction, whereas closest-point distance measures are visualized interactively by evaluating one surface's distance field on the other surface. Furthermore, we use these measures to define a three-way surface segmentation, which visualizes regions of growth, shrinkage, and no change of a test surface compared with a reference surface. Finally, we demonstrate an application of our technique in the visualization of statistical shape models. We evaluate our technique based on feedback provided by medical image analysis researchers, who are experts in working with such models. Show less
Kok, P.; Baiker, M.; Hendriks, E.A.; Post, F.H.; Dijkstra, J.; Lowik, C.W.G.M.; ... ; Botha, C.P. 2010
The analysis of multi-timepoint whole-body small animal CT data is greatly complicated by the varying posture of the subject at different timepoints. Due to these variations, correctly relating and... Show moreThe analysis of multi-timepoint whole-body small animal CT data is greatly complicated by the varying posture of the subject at different timepoints. Due to these variations, correctly relating and comparing corresponding regions of interest is challenging. In addition, occlusion may prevent effective visualization of these regions of interest. To address these problems, we have developed a method that fully automatically maps the data to a standardized layout of sub-volumes, based on an articulated atlas registration. We have dubbed this process articulated planar reformation, or APR. A sub-volume can be interactively selected for closer inspection and can be compared with the corresponding sub-volume at the other timepoints, employing a number of different comparative visualization approaches. We provide an additional tool that highlights possibly interesting areas based on the change of bone density between timepoints. Furthermore we allow visualization of the local registration error, to give an indication of the accuracy of the registration. We have evaluated our approach on a case that exhibits cancer-induced bone resorption. Show less
Krekel, P.R.; Valstar, E.R.; Groot, J. de; Post, F.H.; Nelissen, R.G.H.H.; Botha, C.P. 2010
Kinematics is the analysis of motions without regarding forces or inertial effects, with the purpose of understanding joint behaviour. Kinematic data of linked joints, for example the upper... Show moreKinematics is the analysis of motions without regarding forces or inertial effects, with the purpose of understanding joint behaviour. Kinematic data of linked joints, for example the upper extremity, i.e. the shoulder and arm joints, contains many related degrees of freedom that complicate numerical analysis. Visualisation techniques enhance the analysis process, thus improving the effectiveness of kinematic experiments. This paper describes a new visualisation system specifically designed for the analysis of multi-joint kinematic data of the upper extremity. The challenge inherent in the data is that the upper extremity is comprised of five cooperating joints with a total of fifteen degrees of freedom. The range of motion may be affected by subtle deficiencies of individual joints that are difficult to pinpoint. To highlight these subtleties our approach combines interactive filtering and multiple visualisation techniques. Our system is further differentiated by the fact that it integrates simultaneous acquisition and visual analysis of biokinematic data. Also, to facilitate complex queries, we have designed a visual query interface with visualisation and interaction elements that are based on the domain-specific anatomical representation of the data. The combination of these techniques form an effective approach specifically tailored for the investigation and comparison of large collections of kinematic data. This claim is supported by an evaluation experiment where the technique was used to inspect the kinematics of the left and right arm of a patient with a healed proximal humerus fracture, i.e. a healed shoulder fracture. Show less
Purpose Segmentation of rheumatoid joints from CT images is a complicated task. The pathological state of the joint results in a non-uniform density of the bone tissue, with holes and... Show morePurpose Segmentation of rheumatoid joints from CT images is a complicated task. The pathological state of the joint results in a non-uniform density of the bone tissue, with holes and irregularities complicating the segmentation process. For the specific case of the shoulder joint, existing segmentation techniques often fail and lead to poor results. This paper describes a novel method for the segmentation of these joints. Methods Given a rough surface model of the shoulder, a loop that encircles the joint is extracted by calculating the minimum curvature of the surface model. The intersection points of this loop with the separate CT-slices are connected by means of a path search algorithm. Inaccurate sections are corrected by iteratively applying a Hough transform to the segmentation result. Results As a qualitative measure we calculated the Dice coefficient and Hausdorff distances of the automatic segmentations and expert manual segmentations of CT-scans of ten severely deteriorated shoulder joints. For the humerus and scapula the median Dice coefficient was 98.9% with an interquartile range (IQR) of 95.8-99.4 and 98.5% (IQR 98.3-99.2%), respectively. The median Hausdorff distances were 3.06 mm (IQR 2.30-4.14) and 3.92 mm (IQR 1.96 5.92 mm), respectively. Conclusion The routine satisfies the criterion of our particular application to accurately segment the shoulder joint in under 2 min. We conclude that combining surface curvature, limited user interaction and iterative refinement via a Hough transform forms a satisfactory approach for the segmentation of severely damaged arthritic shoulder joints. Show less