MIA - a free and open source software for gray scale medical image analysis

Background Gray scale images make the bulk of data in bio-medical image analysis, and hence, the main focus of many image processing tasks lies in the processing of these monochrome images. With ever improving acquisition devices, spatial and temporal image resolution increases, and data sets become very large. Various image processing frameworks exists that make the development of new algorithms easy by using high level programming languages or visual programming. These frameworks are also accessable to researchers that have no background or little in software development because they take care of otherwise complex tasks. Specifically, the management of working memory is taken care of automatically, usually at the price of requiring more it. As a result, processing large data sets with these tools becomes increasingly difficult on work station class computers. One alternative to using these high level processing tools is the development of new algorithms in a languages like C++, that gives the developer full control over how memory is handled, but the resulting workflow for the prototyping of new algorithms is rather time intensive, and also not appropriate for a researcher with little or no knowledge in software development. Another alternative is in using command line tools that run image processing tasks, use the hard disk to store intermediate results, and provide automation by using shell scripts. Although not as convenient as, e.g. visual programming, this approach is still accessable to researchers without a background in computer science. However, only few tools exist that provide this kind of processing interface, they are usually quite task specific, and don’t provide an clear approach when one wants to shape a new command line tool from a prototype shell script. Results The proposed framework, MIA, provides a combination of command line tools, plug-ins, and libraries that make it possible to run image processing tasks interactively in a command shell and to prototype by using the according shell scripting language. Since the hard disk becomes the temporal storage memory management is usually a non-issue in the prototyping phase. By using string-based descriptions for filters, optimizers, and the likes, the transition from shell scripts to full fledged programs implemented in C++ is also made easy. In addition, its design based on atomic plug-ins and single tasks command line tools makes it easy to extend MIA, usually without the requirement to touch or recompile existing code. Conclusion In this article, we describe the general design of MIA, a general purpouse framework for gray scale image processing. We demonstrated the applicability of the software with example applications from three different research scenarios, namely motion compensation in myocardial perfusion imaging, the processing of high resolution image data that arises in virtual anthropology, and retrospective analysis of treatment outcome in orthognathic surgery. With MIA prototyping algorithms by using shell scripts that combine small, single-task command line tools is a viable alternative to the use of high level languages, an approach that is especially useful when large data sets need to be processed.

[the Modified Us7 Score In The Assessment Of Synovitis In Early Rheumatoid Arthritis].

To evaluate the modified US7 score (MUS7 score SYN) in the assessment of patients with early rheumatoid arthritis (ERA). In addition, dorsal and palmar recesses of the wrists as well as of small joints of the hands and feet were examined for the presence of synovitis by means of a global assessment of joints. The study sample comprised 32 patients treated for arthritis, with an average disease duration of 13 months. An ultrasound machine with high frequency transducer was used. Hands were also X-rayed and analysed by Larsen score. Out of the 832 examined joints, synovitis was detected in 173 (20,79%), tenosynovitis in 22 (4,91%), and erosions in 3 (1,56%). Synovitis was predominantly detected in the dorsal recess (73,38%) of MCP and PIP joints, when compared with palmar recess (26%). The presence of synovitis in the joints evaluated correlated with clinical (HAQ-DI, DAS28), laboratory (ACPA, RF, CRP), and ultrasound results (r = 0,37 to r = 0,42; p = 0,04 to p = 0,003). We found correlation of the MUS7 score SYN of the gray scale US or of the power Doppler US with DAS28 (PCR) values (r = 0,38; p = 0,0332), and with CRP results (r = 0,39; p = 0,0280), respectively. The dorsal recess, the wrist, and small joints can be considered as important sites to detect synovitis by the MUS7 score SYN in patients with ERA.

WOOD DENSITY VARIATION AND TREE RING DEMARCATION IN Gmelina arborea TREES USING X-RAY DENSITOMETRY

Due to its relationship with other properties, wood density is the main wood quality parameter. Modern, accurate methods - such as X-ray densitometry - are applied to determine the spatial distribution of density in wood sections and to evaluate wood quality. The objectives of this study were to determinate the influence of growing conditions on wood density variation and tree ring demarcation of gmelina trees from fast growing plantations in Costa Rica. The wood density was determined by X-ray densitometry method. Wood samples were cut from gmelina trees and were exposed to low X-rays. The radiographic films were developed and scanned using a 256 gray scale with 1000 dpi resolution and the wood density was determined by CRAD and CERD software. The results showed tree-ring boundaries were distinctly delimited in trees growing in site with rainfall lower than 25 10 mm/year. It was demonstrated that tree age, climatic conditions and management of plantation affects wood density and its variability. The specific effect of variables on wood density was quantified by for multiple regression method. It was determined that tree year explained 25.8% of the total variation of density and 19.9% were caused by climatic condition where the tree growing. Wood density was less affected by the intensity of forest management with 5.9% of total variation.

A critical analysis of Doppler velocimetry in the differential diagnosis of malignant and benign ovarian masses

Objectives: To evaluate the intratumoral reliability of color Doppler parameters and the contribution of Doppler sonography to the gray-scale differential diagnosis of ovarian masses. Methods: An observational study was performed including 67 patients, 15 (22.4%) with malignant ovarian neoplasm and 52 (77.6%) with benign ovarian diseases. We performed the Doppler evaluation in two distinct vessels selected after decreasing the Doppler gain to sample only vessels with higher velocity flow. Doppler measurements were obtained from each identified vessel, and resistive index (RI), pulsatility index (PI), peak systolic velocity (PSV), and end-diastolic velocity (EDV) were measured. Intraclass coefficient of correlation (ICC), sensitivity, specificity, and potential improvement in gray-scale ultrasound performance were calculated. Results: The general ICC were 0.60 (95% CI 0.42- 0.73) for RI, 0.65 (95% CI 0.49- 0.77) for PI, 0.07 (95% CI- 0.17-0.30) for PSV, and 0.19 (95% CI -0.05-0.41) for EDV. The sensitivity and specificity were respectively 84.6% and 86.7% for RI, 69.2% and 93.3% for PI, 80.0% and 65.4% for gray-scale sonography, and 93.3% and 65.4% for gray-scale plus RI (p = 0.013). Conclusions: Gynecologists must be careful in interpreting results from Doppler evaluation of ovarian masses because PSV and EDV present poor intratumoral reliability. The lower RI value, evaluated in at least two distinct sites of the tumor, was able to improve the performance of gray-scale ultrasound in differential diagnosis of ovarian masses.

Low-Level Laser Therapy in Chronic Autoimmune Thyroiditis: A Pilot Study

Background and Objectives: Chronic autoimmune thyroiditis (CAT) remains the most common cause of acquired hypothyroidism There is currently no therapy that is capable of regenerating CAT-damaged thyroid tissue The objective of this study was to gauge the value of applying low-level laser therapy (LLLT) in CAT patients based on both ultrasound studies (USs) and evaluations of thyroid function and thyroid autoantibodies. Study Design/Materials and Methods: Fifteen patients who had hypothyroidism caused by CAT and were undergoing levothyroxine (LT4) treatment were selected to participate in the study Patients received 10 applications of LLLT (830 nm, output power 50 mW) in continuous mode, twice a week, using either the punctual technique (8 patients) or the sweep technique (7 patients), with fluence in the range of 38-108 J/cm(2) USs were performed prior to and 30 days after LLLT USs included a quantitative analysis of echogenicity through a gray-scale computerized histogram index (El). Following the second ultrasound (30 days after LLLT), LT4 was discontinued in all patients and, if required, reintroduced Truodothyronine, thyroxine (T4), free T4, thyrotropin, thyroid peroxidase (TPOAb) and thyroglobulin (TgAb) antibodies levels were assessed before LLLT and then 1, 2, 3, 6, and 9 months after LT4 withdrawal. Results: We noted all patients` reduced LT4 dosage needs, including 7 (47%) who did not require any LT4 through the 9-month follow-up The LT4 dosage used pre-LLLT (96 +/- 22 mu g/day) decreased in the 9th month of follow-up (38 23 mu g/day; P<0.0001) TPOAb levels also decreased (pre-LLLT = 982 +/- 530 U/ml, post-LLLT = 579 454 U/ml, P = 0 016) TgAb levels were not reduced, though we did observe a post-LLLT increase in the EI (pre-LLLT = 0 99 +/- 0.09, post-LLLT= 1.21 +/- 0.19, P=0.001) Conclusion: The preliminary results indicate that LLLT promotes the improvement of thyroid function, as patients experienced a decreased need for LT4, a reduction in TPOAb levels, and an increase in parenchymal echogenicity Lasers Surg. Med. 42:589-596, 2010. (C) 2010 Wiley-Liss, Inc

Abnormal left ventricular filling with increasing age reflects abnormal myocardial characteristics independent of ischemia or hypertrophy

Abnormal left ventricular (IV) filling may occur with increasing age despite apparently normal IV size and function, and is usually attributed to IV hypertrophy and coronary artery disease. The purpose of this study was to determine whether myocardial abnormalities could be identified in 67 such patients (36 men, mean age 57 +/- 9 years) whose IV hypertrophy and coronary artery disease were excluded by dobutamine echocardiography. All patients underwent gray scale and color tissue Doppler imaging from 3 apical views, which were stored and analyzed off line. Disturbances in structure and function were assessed by averaging the cyclic variation of integrated backscatter, strain rate, and peak systolic strain from each myocardial segment. Calibrated integrated backscatter (corrected for pericardial backscatter intensity) was measured in the septum and posterior wall from the parasternal long-axis view. Abnormal IV filling was present in 36 subjects (54%). Subjects with and without abnormal IV filling had similar IV mass, but differed in age (p <0.01), cyclic variation (p = 0.001), strain rate (p <0.01), and peak systolic strain (p <0.001). Multivariate logistic regression analysis demonstrated that age (p = 0.016) and cyclic variation (p = 0.042) were the most important determinants of abnormal IV filling in these apparently normal subjects. (C) 2003 by Excerpta Medica, Inc.

Comparison between MDCT and Grayscale IVUS in a Quantitative Analysis of Coronary Lumen in Segments with or without Atherosclerotic Plaques

Background: The diagnostic accuracy of 64-slice MDCT in comparison with IVUS has been poorly described and is mainly restricted to reports analyzing segments with documented atherosclerotic plaques. Objectives: We compared 64-slice multidetector computed tomography (MDCT) with gray scale intravascular ultrasound (IVUS) for the evaluation of coronary lumen dimensions in the context of a comprehensive analysis, including segments with absent or mild disease. Methods: The 64-slice MDCT was performed within 72 h before the IVUS imaging, which was obtained for at least one coronary, regardless of the presence of luminal stenosis at angiography. A total of 21 patients were included, with 70 imaged vessels (total length 114.6 ± 38.3 mm per patient). A coronary plaque was diagnosed in segments with plaque burden > 40%. Results: At patient, vessel, and segment levels, average lumen area, minimal lumen area, and minimal lumen diameter were highly correlated between IVUS and 64-slice MDCT (p < 0.01). However, 64-slice MDCT tended to underestimate the lumen size with a relatively wide dispersion of the differences. The comparison between 64-slice MDCT and IVUS lumen measurements was not substantially affected by the presence or absence of an underlying plaque. In addition, 64-slice MDCT showed good global accuracy for the detection of IVUS parameters associated with flow-limiting lesions. Conclusions: In a comprehensive, multi-territory, and whole-artery analysis, the assessment of coronary lumen by 64-slice MDCT compared with coronary IVUS showed a good overall diagnostic ability, regardless of the presence or absence of underlying atherosclerotic plaques.

Persistence of ultrasound synovitis in patients with rheumatoid arthritis fulfilling the DAS28 and/or the new ACR/EULAR RA remission definitions: results of an observational cohort study.

OBJECTIVE: The primary aim of the study was to evaluate whether rheumatoid arthritis (RA) patients considered to be in remission according to clinical criteria sets still had persisting ultrasound (US) synovitis. We further intended to evaluate the capacity of our US score to discriminate between the patients with a clinically active disease versus those in remission. METHODS: This is an observational study nested within the Swiss Clinical Quality Management in Rheumatic Diseases (SCQM) rheumatoid arthritis cohort. A validated US score (SONAR score) based on a semi-quantitative B-mode and Doppler (PwD) score as part of the regular clinical workup by rheumatologists in different clinical settings was used. To define clinically relevant synovitis, the same score was applied to 38 healthy controls and the 90st percentile was used as cut-off for 'relevant' synovitis. RESULTS: Three hundred and seven patients had at least one US examination and concomitant clinical information on disease activity. More than a third of patients in both DAS28 and ACR/EULAR remission showed significant gray scale synovitis (P=0.01 and 0.0002, respectively) and PwD activity (P=0.005 and 0.0005, respectively) when compared to controls. The capacity of US to discriminate between the two clinical remission groups and patients with active disease was only moderate. CONCLUSION: This observational study confirms that many patients considered to be in clinical remission according the DAS and the ACR/EULAR definitions still have residual synovitis on US. The prognostic significance of US synovitis and the exact place of US in patients reaching clinical remission need to be further evaluated.

Plaque characteristics and arterial remodeling in coronary and peripheral arterial systems.

BACKGROUND: Few studies have examined plaque characteristics among multiple arterial beds in vivo. The purpose of this study was to compare the plaque morphology and arterial remodeling between coronary and peripheral arteries using gray-scale and radiofrequency intravascular ultrasound (IVUS) at clinical presentation. METHODS AND RESULTS: IVUS imaging was performed in 68 patients with coronary and 93 with peripheral artery lesions (29 carotid, 50 renal, and 14 iliac arteries). Plaques were classified as fibroatheroma (VH-FA) (further subclassified as thin-capped [VH-TCFA] and thick-capped [VH-ThCFA]), fibrocalcific plaque (VH-FC) and pathological intimal thickening (VH-PIT). Plaque rupture (13% of coronary, 7% of carotid, 6% of renal, and 7% of iliac arteries; P = NS) and VH-TCFA (37% of coronary, 24% of carotid, 16% of renal, and 7% of iliac arteries; P = 0.02) were observed in all arteries. Compared with coronary arteries, VH-FA was less frequently observed in renal (P < 0.001) and iliac arteries (P < 0.006). Lesions with positive remodeling demonstrated more characteristics of VH-FA in coronary (84% vs. 25%, P < 0.001), carotid (72% vs. 20%, P = 0.001), and renal arteries (42% vs. 4%, P = 0.001) compared with those with intermediate/negative remodeling. There was positive relationship between remodeling index and percent necrotic area in all four arteries. CONCLUSIONS: Atherosclerotic plaque phenotypes were heterogeneous among four different arteries; renal and iliac arteries had more stable phenotypes compared with coronary artery. In contrast, the associations of remodeling pattern with plaque phenotype and composition were similar among the various arterial beds.

A post-processing pipeline to reconstruct the developing fetal brain using low-resolution MRI

Introduction. Development of the fetal brain surfacewith concomitant gyrification is one of the majormaturational processes of the human brain. Firstdelineated by postmortem studies or by ultrasound, MRIhas recently become a powerful tool for studying in vivothe structural correlates of brain maturation. However,the quantitative measurement of fetal brain developmentis a major challenge because of the movement of the fetusinside the amniotic cavity, the poor spatial resolution,the partial volume effect and the changing appearance ofthe developing brain. Today extensive efforts are made todeal with the âeurooepost-acquisitionâeuro reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution (Rousseau, F., 2006;Jiang, S., 2007). We here propose a framework devoted tothe segmentation of the basal ganglia, the gray-whitetissue segmentation, and in turn the 3D corticalreconstruction of the fetal brain. Method. Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences in fetuses aged from 29 to 32gestational weeks (slice thickness 5.4mm, in planespatial resolution 1.09mm). For each fetus, 6 axialvolumes shifted by 1 mm were acquired (about 1 min pervolume). First, each volume is manually segmented toextract fetal brain from surrounding fetal and maternaltissues. Inhomogeneity intensity correction and linearintensity normalization are then performed. A highspatial resolution image of isotropic voxel size of 1.09mm is created for each fetus as previously published byothers (Rousseau, F., 2006). B-splines are used for thescattered data interpolation (Lee, 1997). Then, basalganglia segmentation is performed on this superreconstructed volume using active contour framework witha Level Set implementation (Bach Cuadra, M., 2010). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed (Bach Cuadra, M., 2009). Theresulting white matter image is then binarized andfurther given as an input in the Freesurfer software(http://surfer.nmr.mgh.harvard.edu/) to provide accuratethree-dimensional reconstructions of the fetal brain.Results. High-resolution images of the cerebral fetalbrain, as obtained from the low-resolution acquired MRI,are presented for 4 subjects of age ranging from 29 to 32GA. An example is depicted in Figure 1. Accuracy in theautomated basal ganglia segmentation is compared withmanual segmentation using measurement of Dice similarity(DSI), with values above 0.7 considering to be a verygood agreement. In our sample we observed DSI valuesbetween 0.785 and 0.856. We further show the results ofgray-white matter segmentation overlaid on thehigh-resolution gray-scale images. The results arevisually checked for accuracy using the same principlesas commonly accepted in adult neuroimaging. Preliminary3D cortical reconstructions of the fetal brain are shownin Figure 2. Conclusion. We hereby present a completepipeline for the automated extraction of accuratethree-dimensional cortical surface of the fetal brain.These results are preliminary but promising, with theultimate goal to provide âeurooemovieâeuro of the normal gyraldevelopment. In turn, a precise knowledge of the normalfetal brain development will allow the quantification ofsubtle and early but clinically relevant deviations.Moreover, a precise understanding of the gyraldevelopment process may help to build hypotheses tounderstand the pathogenesis of several neurodevelopmentalconditions in which gyrification have been shown to bealtered (e.g. schizophrenia, autismâeuro¦). References.Rousseau, F. (2006), 'Registration-Based Approach forReconstruction of High-Resolution In Utero Fetal MR Brainimages', IEEE Transactions on Medical Imaging, vol. 13,no. 9, pp. 1072-1081. Jiang, S. (2007), 'MRI of MovingSubjects Using Multislice Snapshot Images With VolumeReconstruction (SVR): Application to Fetal, Neonatal, andAdult Brain Studies', IEEE Transactions on MedicalImaging, vol. 26, no. 7, pp. 967-980. Lee, S. (1997),'Scattered data interpolation with multilevel B-splines',IEEE Transactions on Visualization and Computer Graphics,vol. 3, no. 3, pp. 228-244. Bach Cuadra, M. (2010),'Central and Cortical Gray Mater Segmentation of MagneticResonance Images of the Fetal Brain', ISMRM Conference.Bach Cuadra, M. (2009), 'Brain tissue segmentation offetal MR images', MICCAI.

Hardcopy quality parameters to ensure structures detection at digital mammography

Objective To develop procedures to ensure consistency of printing quality of digital images, by means of hardcopy quantitative analysis based on a standard image. Materials and Methods Characteristics of mammography DI-ML and general purpose DI-HL films were studied through the QC-Test utilizing different processing techniques in a FujiFilm®-DryPix4000 printer. A software was developed for sensitometric evaluation, generating a digital image including a gray scale and a bar pattern to evaluate contrast and spatial resolution. Results Mammography films showed maximum optical density of 4.11 and general purpose films, 3.22. The digital image was developed with a 33-step wedge scale and a high-contrast bar pattern (1 to 30 lp/cm) for spatial resolution evaluation. Conclusion Mammographic films presented higher values for maximum optical density and contrast resolution as compared with general purpose films. The utilized digital processing technique could only change the image pixels matrix values and did not affect the printing standard. The proposed digital image standard allows greater control of the relationship between pixels values and optical density obtained in the analysis of films quality and printing systems.

Reconstruction and analysis of surface variation using photometric stereo

In many industrial applications, accurate and fast surface reconstruction is essential for quality control. Variation in surface finishing parameters, such as surface roughness, can reflect defects in a manufacturing process, non-optimal product operational efficiency, and reduced life expectancy of the product. This thesis considers reconstruction and analysis of high-frequency variation, that is roughness, on planar surfaces. Standard roughness measures in industry are calculated from surface topography. A fast and non-contact method to obtain surface topography is to apply photometric stereo in the estimation of surface gradients and to reconstruct the surface by integrating the gradient fields. Alternatively, visual methods, such as statistical measures, fractal dimension and distance transforms, can be used to characterize surface roughness directly from gray-scale images. In this thesis, the accuracy of distance transforms, statistical measures, and fractal dimension are evaluated in the estimation of surface roughness from gray-scale images and topographies. The results are contrasted to standard industry roughness measures. In distance transforms, the key idea is that distance values calculated along a highly varying surface are greater than distances calculated along a smoother surface. Statistical measures and fractal dimension are common surface roughness measures. In the experiments, skewness and variance of brightness distribution, fractal dimension, and distance transforms exhibited strong linear correlations to standard industry roughness measures. One of the key strengths of photometric stereo method is the acquisition of higher frequency variation of surfaces. In this thesis, the reconstruction of planar high-frequency varying surfaces is studied in the presence of imaging noise and blur. Two Wiener filterbased methods are proposed of which one is optimal in the sense of surface power spectral density given the spectral properties of the imaging noise and blur. Experiments show that the proposed methods preserve the inherent high-frequency variation in the reconstructed surfaces, whereas traditional reconstruction methods typically handle incorrect measurements by smoothing, which dampens the high-frequency variation.

Color ad spectral image assessment using novel quality and fidelity techniques

The ongoing development of the digital media has brought a new set of challenges with it. As images containing more than three wavelength bands, often called spectral images, are becoming a more integral part of everyday life, problems in the quality of the RGB reproduction from the spectral images have turned into an important area of research. The notion of image quality is often thought to comprise two distinctive areas – image quality itself and image fidelity, both dealing with similar questions, image quality being the degree of excellence of the image, and image fidelity the measure of the match of the image under study to the original. In this thesis, both image fidelity and image quality are considered, with an emphasis on the influence of color and spectral image features on both. There are very few works dedicated to the quality and fidelity of spectral images. Several novel image fidelity measures were developed in this study, which include kernel similarity measures and 3D-SSIM (structural similarity index). The kernel measures incorporate the polynomial, Gaussian radial basis function (RBF) and sigmoid kernels. The 3D-SSIM is an extension of a traditional gray-scale SSIM measure developed to incorporate spectral data. The novel image quality model presented in this study is based on the assumption that the statistical parameters of the spectra of an image influence the overall appearance. The spectral image quality model comprises three parameters of quality: colorfulness, vividness and naturalness. The quality prediction is done by modeling the preference function expressed in JNDs (just noticeable difference). Both image fidelity measures and the image quality model have proven to be effective in the respective experiments.

Wrist ultrasound analysis of patients with early rheumatoid arthritis

In the present study, we evaluated 42 wrists using the semi-quantitative scales power Doppler ultrasound (PDUS) and gray scale ultrasound (GSUS) with scores ranging from 0 to 3 and correlated the results with clinical, laboratory and radiographic data. Twenty-one patients (17 women and 4 men) with rheumatoid arthritis according to criteria of the American College of Rheumatology were enrolled in the study from September 2008 to July 2009 at Universidade Estadual de Campinas (UNICAMP). The average disease duration was 14 months. The patients were 66.6% Caucasians and 33.3% non-Caucasians, with a mean age of 42 and 41 years, respectively. A dorsal longitudinal scan was performed by ultrasound on the radiocarpal and midcarpal joints using GE LOGIQ XP-linear ultrasound and a high frequency (8-10 MHz) transducer. All patients were X-rayed, and the Larsen score was determined for the joints, with grades ranging from 0 to V. This study showed significant correlations between clinical, sonographic and laboratory data: GSUS and swollen right wrist (r = 0.546), GSUS of right wrist and swelling of left wrist (r = 0.511), PDUS of right wrist and pain in left wrist (r = 0.436), PDUS of right wrist and C-reactive protein (r = 0.466). Ultrasound can be considered a useful tool in the diagnosis of synovitis in early rheumatoid arthritis mainly when the anti-cyclic citrullinated peptide and rheumatoid factor are negative, and can lead to an early change in the therapeutic decision.