Hideaki Haneishi

Independent-component analysis of skin color image

The spatial distributions of melanin and hemoglobin in human skin are separated by independent-component analysis of a skin color image. The analysis is based on the skin color model with three assumptions: (1) Spatial variation of color in the skin is caused by two pigments, melanin and hemoglobin; (2) the quantities of the two pigments are mutually independent spatially; and (3) linearity holds among the quantities and the observed color signals in the optical density domain. The results of the separation agree well with physiological knowledge. The separated components are synthesized to simulate the various facial color images by changing the quantities of the two separated pigments.

System design for accurately estimating the spectral reflectance of art paintings

Accurately estimating the spectral reflectance of art paintings from low-dimensional multichannel images requires that both image-acquisition hardware with appropriate spectral characteristics and appropriate estimation software be applied to the captured multichannel image. In this study, a system that incorporates both factors is designed and developed on the basis of the minimum-mean-square error criterion. The accuracy of spectral estimation by use of this system is evaluated, and the systems high performance is demonstrated.

Development of a small prototype for a proof-of-concept of OpenPET imaging

The OpenPET geometry is our new idea to visualize a physically opened space between two detector rings. In this paper, we developed the first small prototype to show a proof-of-concept of OpenPET imaging. Two detector rings of 110 mm diameter and 42 mm axial length were placed with a gap of 42 mm. The basic imaging performance was confirmed through phantom studies; the open imaging was realized at the cost of slight loss of axial resolution and 24% loss of sensitivity. For a proof-of-concept of PET image-guided radiation therapy, we carried out the in-beam tests with (11)C radioactive beam irradiation in the heavy ion medical accelerator in Chiba to visualize in situ distribution of primary particles stopped in a phantom. We showed that PET images corresponding to dose distribution were obtained. For an initial proof-of-concept of real-time multimodal imaging, we measured a tumor-inoculated mouse with (18)F-FDG, and an optical image of the mouse body surface was taken during the PET measurement by inserting a digital camera in the ring gap. We confirmed that the tumor in the gap was clearly visualized. The result also showed the extension effect of an axial field-of-view (FOV); a large axial FOV of 126 mm was obtained with the detectors that originally covered only an 84 mm axial FOV. In conclusion, our initial imaging studies showed promising performance of the OpenPET.

Transaxial system models for jPET-D4 image reconstruction

A high-performance brain PET scanner, jPET-D4, which provides four-layer depth-of-interaction (DOI) information, is being developed to achieve not only high spatial resolution, but also high scanner sensitivity. One technical issue to be dealt with is the data dimensions which increase in proportion to the square of the number of DOI layers. It is, therefore, difficult to apply algebraic or statistical image reconstruction methods directly to DOI-PET, though they improve image quality through accurate system modelling. The process that requires the most computational time and storage space is the calculation of the huge number of system matrix elements. The DOI compression (DOIC) method, which we have previously proposed, reduces data dimensions by a factor of 1/5. In this paper, we propose a transaxial imaging system model optimized for jPET-D4 with the DOIC method. The proposed model assumes that detector response functions (DRFs) are uniform along line-of-responses (LORs). Then each element of the system matrix is calculated as the summed intersection lengths between a pixel and sub-LORs weighted by a value from the DRF look-up-table. 2D numerical simulation results showed that the proposed model cut the calculation time by a factor of several hundred while keeping image quality, compared with the accurate system model. A 3D image reconstruction with the on-the-fly calculation of the system matrix is within the practical limitations by incorporating the proposed model and the DOIC method with one-pass accelerated iterative methods.

Color digital halftoning taking colorimetric color reproduction into account

Taking colorimetric color reproduction into account, the conventional error diffusion method is modified for color digital half-toning. Assuming that the input to a bilevel color printer is given in CIE-XYZ tristimulus values or CIE-LAB values instead of the more conventional RGB or YMC values, two modified versions based on vector operation in (1) the XYZ color space and (2) the LAB color space were tested. Experimental results show that the modified methods, especially the method using the LAB color space, resulted in better color reproduction performance than the conventional methods. Spatial artifacts that appear in the modified methods are presented and analyzed. It is also shown that the modified method (2) with a thresholding technique achieves a good spatial image quality.

Nonlinear motion correction of respiratory-gated lung SPECT images

We propose a method for correcting the motion of the lungs between different phase images obtained by respiratory-gated single photon emission computed tomography (SPECT). This method is applied to SPECT images that show a preserved activity distribution in the lungs such as 99m-Tc macro aggregated albumin (99m-Tc-MAA) perfusion images and 99m-Tc-Technegas ventilation images. In the proposed method, an objective function, which consists of both the degree of similarity between a reference image and a deformed image, and the smoothness of deformation is defined and optimized using a simulated annealing algorithm. For the degree of similarity term in the objective function, an expansion ratio, defined as the ratio of change in local volume due to deformation, is introduced to preserve the total activity during the motion correction process. This method was applied to data simulated from computer phantoms, data acquired from a physical phantom, and 17 sets of clinical data. In all cases, the motion correction between inspiration and expiration phase images was successfully achieved.

DOI-PET image reconstruction with accurate system modeling that reduces redundancy of the imaging system

A high-performance positron emission tomography (PET) scanner, which measures depth-of-interaction (DOI) information, is under development at the National Institute of Radiological Sciences in Japan. Image reconstruction methods with accurate modeling of the system response functions have been successfully used to improve PET image quality. It is, however, difficult to apply these methods to the DOI-PET scanner because the dimension of DOI-PET data increases in proportion to the square of the number of DOI layers. In this paper, we propose a compressed imaging system model for DOI-PET image reconstruction, in order to reduce computational cost while keeping image quality. The basic idea of the proposed method is that the DOI-PET imaging system is highly redundant. First, DOI-PET data is transformed into compact data so that data bins with highly correlating sensitivity functions are combined. Then image reconstruction methods based on accurate system modeling, such as the maximum likelihood expectation maximization (ML-EM), are applied. The proposed method was applied to simulated data for the DOI-PET scanner operated in 2-D mode. Then the tradeoff between the background noise and the spatial resolution was investigated. Numerical simulation results show that the proposed method followed by ML-EM reduces computational cost effectively while keeping the advantages of the accurate system modeling and DOI information.

Preliminary resolution performance of the prototype system for a 4-Layer DOI-PET scanner: jPET-D4

We are developing a high-performance brain PET scanner, jPET-D4, which provides 4-layer depth-of-interaction (DOI) information. The scanner is designed to achieve not only high spatial resolution but also high scanner sensitivity with the DOI information obtained from multi-layered thin crystals. The scanner has 5 rings of 24 detector blocks each, and each block consists of 1024 GSO crystals of 2.9 mm/spl times/2.9 mm/spl times/7.5 mm, which are arranged in 4 layers of 16/spl times/16 arrays. At this stage, a pair of detector blocks and a coincidence circuit have been assembled into an experimental prototype gantry. In this paper, as a preliminary experiment, we investigated the performance of the jPET-D4s spatial resolution using the prototype system. First, spatial resolution was measured from a filtered backprojection reconstructed image. To avoid systematic error and reduce computational cost in image reconstruction, we applied the DOI compression (DOIC) method followed by maximum likelihood expectation maximization that we had previously proposed. Trade-off characteristics between background noise and resolution were investigated because improved spatial resolution is possible only when enhanced noise is avoided. Experimental results showed that the jPET-D4 achieves better than 3 mm spatial resolution over the field-of-view.

Multiple current dipole estimation using simulated annealing

A method for estimating electrical current distribution in the human brain using a multiple current dipole model is presented. A cost function for estimating multiple dipoles is proposed and a simulated annealing algorithm is used to obtain an acceptable solution. Computer simulation is used to evaluate the effectiveness of this method.<<ETX>>

Details of simulated annealing algorithm to estimate parameters of multiple current dipoles using biomagnetic data

The details of the simulated annealing algorithm proposed to estimate the parameters of multiple current dipoles using biomagnetic data are described. The effects of the choices of such numerical conditions as the amount of estimate transitions, the equilibrium criterion, and the temperature decrement on the algorithms performance are discussed. Incorrect results from the computer simulation obtained with inappropriate choices of such conditions are presented. The verification of a near-optimum convergence by reheating and reannealing is demonstrated. A modification of this algorithm is proposed for implementation by parallel computer, and the validity of this parallel algorithm is demonstrated by computer simulation.