A new set distance and its application to shape registration

We propose a new distance measure, called Complement weighted sum of minimal distances, between finite sets in ${\mathbb Z }^n$ and evaluate its usefulness for shape registration and matching. In this set distance the contribution of each point of each set is weighted according to its distance to the complement of the set. In this way, outliers and noise contribute less to the new similarity measure. We evaluate the performance of the new set distance for registration of shapes in binary images and compare it to a number of often used set distances found in the literature. The most extensive evaluation uses a set of synthetic 2D images. We also show three examples of real problems: registering a set of 2D images extracted from synchrotron radiation micro-computed tomography (SR $\upmu $ CT) volumes depicting bone implants; the difficult multi-modal registration task of finding the exact location of a 2D slice of a bone implant, as imaged by a light microscope, within a 3D SR $\upmu $ CT volume of the same implant; and finally recognition of handwritten characters. The evaluation shows that our new set distance performs well for all tasks and outperforms the other observed distance measures in most cases. It is therefore useful in many image registration and shape comparison tasks.     

Extracting 3D information on bone remodeling in the proximity of titanium implants in SRμCT image volumes

Bone-implant integration is measured in several ways. Traditionally and routinely, 2D histological sections of samples, containing bone and the biomaterial, are stained and analyzed using a light microscope. Such histological section provides detailed cellular information about the bone regeneration in the proximity of the implant. However, this information reflects the integration in only a very small fraction, a 10 μm thick slice, of the sample. In this study, we show that feature values quantified on 2D sections are highly dependent on the orientation and the placement of the section, suggesting that a 3D analysis of the whole sample is of importance for a more complete judgment of the bone structure in the proximity of the implant. We propose features describing the 3D data by extending the features traditionally used for 2D-analysis. We present a method for extracting these features from 3D image data and we measure them on five 3D SRμCT image volumes.We also simulate cuts through the image volume positioned at all possible section positions. These simulations show that the measurement variations due to the orientation of the section around the center line of the implant are about 30%.     

Segmentation of virus particle candidates in transmission electron microscopy images

In this paper, we present an automatic segmentation method that detects virus particles of various shapes in transmission electron microscopy images. The method is based on a statistical analysis of local neighbourhoods of all the pixels in the image followed by an object width discrimination and finally, for elongated objects, a border refinement step. It requires only one input parameter, the approximate width of the virus particles searched for. The proposed method is evaluated on a large number of viruses. It successfully segments viruses regardless of shape, from polyhedral to highly pleomorphic.     

A refined circular template matching method for classification of human cytomegalovirus capsids in TEM images

An automatic image analysis method for describing, segmenting, and classifying human cytomegalovirus capsids in transmission electron micrograph (TEM) images of host cell nuclei has been developed. Three stages of the capsid assembly process in the host cell nucleus have been investigated. Each class is described by a radial density profile, which is the average grey-level at each radial distance from the center. A template, constructed from the profile, is used to find possible capsid locations by correlation based matching. The matching results are further refined by size and distortion analysis of each possible capsid, resulting in a final segmentation and classification.     

On Reversible Skeletonization Using Anchor-Points from Distance Transforms

In many applications thinning of objects is of great interest. We here present a skeletonization algorithm that is based on the idea of iteratively thinning the distance transform of an object layer by layer until either an anchor-point is reached or the connectivity breaks. Our algorithm is general in the sense that any metric and any connectivity can be used. Also, it is based on ideas that are not specific for 2D. The properties of the resulting skeletons are evaluated according to the “Lee–Lam–Suen properties.”     

Distance transforms for three-dimensional grids with non-cubic voxels

Distance transforms on the face-centered cubic (fcc) grid and the body-centered cubic (bcc) grid are examined. Since the voxels on the fcc and bcc grids are better approximations of a Euclidean ball than the cube, the distance transforms (DTs) on these grids can be less rotation dependent than those in , which is a desirable feature. Optimal (according to the error function) weights are calculated and integer approximations of these weights are found. Also, the two-dimensional city block distance is generalized to the fcc and bcc grids by considering a unit distance between gridpoints whose corresponding voxels share a face. A method to compute the DTs is presented. The results are evaluated both theoretically and by actually computing some DTs.     

Minimizing scanning electron microscope artefacts by filter design

A new type of non‐linear filter for digital images has been developed. By using distance transforms we estimate the average point spread function for a set of fibre cross‐sectional images. Then a fast filter technique, based on lookup tables for distance layers, attenuates the uneven background response from the scanning electron microscope. Compared to the convolution‐based techniques that we tried, this approach caused less blurring effects on our fibre images and also made the background pixels more homogeneous. The only assumption we make is that we can roughly segment the background pixels by using a pixel‐wise classifier. Although the assumption that the uneven background response can be described by a circular point spread function is only approximately true in the case discussed here, this method greatly attenuates the effect and provides a fast and general filtering method that can also be of use for other applications.     

Comparison of histomorphometrical data obtained with two different image analysis methods

A common way to determine tissue acceptance of biomaterials is to perform histomorphometrical analysis on histologically stained sections from retrieved samples with surrounding tissue, using various methods. The “time and money consuming” methods and techniques used are often “in house standards”. We address light microscopic investigations of bone tissue reactions on un-decalcified cut and ground sections of threaded implants. In order to screen sections and generate results faster, the aim of this pilot project was to compare results generated with the in-house standard visual image analysis tool (i.e., quantifications and judgements done by the naked eye) with a custom made automatic image analysis program. The histomorphometrical bone area measurements revealed no significant differences between the methods but the results of the bony contacts varied significantly. The raw results were in relative agreement, i.e., the values from the two methods were proportional to each other: low bony contact values in the visual method corresponded to low values with the automatic method. With similar resolution images and further improvements of the automatic method this difference should become insignificant. A great advantage using the new automatic image analysis method is that it is time saving—analysis time can be significantly reduced.     

Hierarchical chamfer matching: a parametric edge matching algorithm

The algorithm matches edges by minimizing a generalized distance between them. The matching is performed in a series of images depicting the same scene with different resolutions, i.e. in a resolution pyramid. Using this hierarchical structure reduces the computational load significantly. The algorithm is reasonably simple to implement and is insensitive to noise and other disturbances. The algorithm has been tested in several applications. Two of them are briefly presented. In the first application the outlines of common tools are matched to gray-level images of the same tools, with overlapping. In the second application lake edges from aerial photographs are matched to lake edges from a map, with translation, rotation, scale, and perspective changes. The hierarchical chamfer matching algorithm gives correct results using a reasonable amount of computational resources in all tested applications     

Automatic measurement of compression wood cell attributes in fluorescence microscopy images

This paper presents a new automated method for analyzing compression wood fibers in fluorescence microscopy. Abnormal wood known as compression wood is present in almost every softwood tree harvested. Compression wood fibers show a different cell wall morphology and chemistry compared to normal wood fibers, and their mechanical and physical characteristics are considered detrimental for both construction wood and pulp and paper purposes. Currently there is the need for improved methodologies for characterization of lignin distribution in wood cell walls, such as from compression wood fibers, that will allow for a better understanding of fiber mechanical properties. Traditionally, analysis of fluorescence microscopy images of fiber cross-sections has been done manually, which is time consuming and subjective. Here, we present an automatic method, using digital image analysis, that detects and delineates softwood fibers in fluorescence microscopy images, dividing them into cell lumen, normal and highly lignified areas. It also quantifies the different areas, as well as measures cell wall thickness. The method is evaluated by comparing the automatic with a manual delineation. While the boundaries between the various fiber wall regions are detected using the automatic method with precision similar to inter and intra expert variability, the position of the boundary between lumen and the cell wall has a systematic shift that can be corrected. Our method allows for transverse structural characterization of compression wood fibers, which may allow for improved understanding of the micro-mechanical modeling of wood and pulp fibers.