Characterization and classification of tumor lesions using computerized fractal-based texture analysis and support vector machines in digital mammograms

Objective: This paper presents a detailed study of fractal-based methods for texture characterization of mammographic mass lesions and architectural distortion. The purpose of this study is to explore the use of fractal and lacunarity analysis for the characterization and classification of both tumor lesions and normal breast parenchyma in mammography. Materials and methods: We conducted comparative evaluations of five popular fractal dimension estimation methods for the characterization of the texture of mass lesions and architectural distortion. We applied the concept of lacunarity to the description of the spatial distribution of the pixel intensities in mammographic images. These methods were tested with a set of 57 breast masses and 60 normal breast parenchyma (dataset1), and with another set of 19 architectural distortions and 41 normal breast parenchyma (dataset2). Support vector machines (SVM) were used as a pattern classification method for tumor classification. Results: Experimental results showed that the fractal dimension of region of interest (ROIs) depicting mass lesions and architectural distortion was statistically significantly lower than that of normal breast parenchyma for all five methods. Receiver operating characteristic (ROC) analysis showed that fractional Brownian motion (FBM) method generated the highest area under ROC curve (A z = 0.839 for dataset1, 0.828 for dataset2, respectively) among five methods for both datasets. Lacunarity analysis showed that the ROIs depicting mass lesions and architectural distortion had higher lacunarities than those of ROIs depicting normal breast parenchyma. The combination of FBM fractal dimension and lacunarity yielded the highest A z value (0.903 and 0.875, respectively) than those based on single feature alone for both given datasets. The application of the SVM improved the performance of the fractal-based features in differentiating tumor lesions from normal breast parenchyma by generating higher A z value. Conclusion: FBM texture model is the most appropriate model for characterizing mammographic images due to self-affinity assumption of the method being a better approximation. Lacunarity is an effective counterpart measure of the fractal dimension in texture feature extraction in mammographic images. The classification results obtained in this work suggest that the SVM is an effective method with great potential for classification in mammographic image analysis.

Investigation of support vector machine for the detection of architectural distortion in mammographic images

This paper investigates detection of architectural distortion in mammographic images using support vector machine. Hausdorff dimension is used to characterise the texture feature of mammographic images. Support vector machine, a learning machine based on statistical learning theory, is trained through supervised learning to detect architectural distortion. Compared to the Radial Basis Function neural networks, SVM produced more accurate classification results in distinguishing architectural distortion abnormality from normal breast parenchyma.

Sensitivity of feature-based analysis methods of storm tracks to the form of background field removal

For the tracking of extrema associated with weather systems to be applied to a broad range of fields it is necessary to remove a background field that represents the slowly varying, large spatial scales. The sensitivity of the tracking analysis to the form of background field removed is explored for the Northern Hemisphere winter storm tracks for three contrasting fields from an integration of the U. K. Met Office's (UKMO) Hadley Centre Climate Model (HadAM3). Several methods are explored for the removal of a background field from the simple subtraction of the climatology, to the more sophisticated removal of the planetary scales. Two temporal filters are also considered in the form of a 2-6-day Lanczos filter and a 20-day high-pass Fourier filter. The analysis indicates that the simple subtraction of the climatology tends to change the nature of the systems to the extent that there is a redistribution of the systems relative to the climatological background resulting in very similar statistical distributions for both positive and negative anomalies. The optimal planetary wave filter removes total wavenumbers less than or equal to a number in the range 5-7, resulting in distributions more easily related to particular types of weather system. For the temporal filters the 2-6-day bandpass filter is found to have a detrimental impact on the individual weather systems, resulting in the storm tracks having a weak waveguide type of behavior. The 20-day high-pass temporal filter is less aggressive than the 2-6-day filter and produces results falling between those of the climatological and 2-6-day filters.

Adaptive constraints for feature tracking

In this paper extensions to an existing tracking algorithm are described. These extensions implement adaptive tracking constraints in the form of regional upper-bound displacements and an adaptive track smoothness constraint. Together, these constraints make the tracking algorithm more flexible than the original algorithm (which used fixed tracking parameters) and provide greater confidence in the tracking results. The result of applying the new algorithm to high-resolution ECMWF reanalysis data is shown as an example of its effectiveness.

Feature-point detection using distance transforms: Application to tracking tropical convective complexes

The identification, tracking, and statistical analysis of tropical convective complexes using satellite imagery is explored in the context of identifying feature points suitable for tracking. The feature points are determined based on the shape of complexes using the distance transform technique. This approach has been applied to the determination feature points for tropical convective complexes identified in a time series of global cloud imagery. The feature points are used to track the complexes, and from the tracks statistical diagnostic fields are computed. This approach allows the nature and distribution of organized deep convection in the Tropics to be explored.

Feature tracking on the unit-sphere

Techniques used in a previous study of the objective identification and tracking of meteorological features in model data are extended to the unit sphere. An alternative feature detection scheme is described based on cubic interpolation for the sphere and local maximization. The extension of the tracking technique, used in the previous study, to the unit sphere is described. An example of the application of these techniques to a global relative vorticity field from a model integration are presented and discussed.

A comparison of recent reanalysis datasets using objective feature tracking: Storm tracks and tropical easterly waves

Data from four recent reanalysis projects [ECMWF, NCEP-NCAR, NCEP - Department of Energy ( DOE), NASA] have been diagnosed at the scale of synoptic weather systems using an objective feature tracking method. The tracking statistics indicate that, overall, the reanalyses correspond very well in the Northern Hemisphere (NH) lower troposphere, although differences for the spatial distribution of mean intensities show that the ECMWF reanalysis is systematically stronger in the main storm track regions but weaker around major orographic features. A direct comparison of the track ensembles indicates a number of systems with a broad range of intensities that compare well among the reanalyses. In addition, a number of small-scale weak systems are found that have no correspondence among the reanalyses or that only correspond upon relaxing the matching criteria, indicating possible differences in location and/or temporal coherence. These are distributed throughout the storm tracks, particularly in the regions known for small-scale activity, such as secondary development regions and the Mediterranean. For the Southern Hemisphere (SH), agreement is found to be generally less consistent in the lower troposphere with significant differences in both track density and mean intensity. The systems that correspond between the various reanalyses are considerably reduced and those that do not match span a broad range of storm intensities. Relaxing the matching criteria indicates that there is a larger degree of uncertainty in both the location of systems and their intensities compared with the NH. At upper-tropospheric levels, significant differences in the level of activity occur between the ECMWF reanalysis and the other reanalyses in both the NH and SH winters. This occurs due to a lack of coherence in the apparent propagation of the systems in ERA15 and appears most acute above 500 hPa. This is probably due to the use of optimal interpolation data assimilation in ERA15. Also shown are results based on using the same techniques to diagnose the tropical easterly wave activity. Results indicate that the wave activity is sensitive not only to the resolution and assimilation methods used but also to the model formulation.

The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

Effect of combined heat and high-pressure treatments on the texture of chicken breast muscle (Pectoralis fundus)

Commercially supplied chicken breast muscle was subjected to simultaneous heat and pressure treatments. Treatment conditions ranged from ambient temperature to 70 °C and from 0.1 to 800 MPa, respectively, in various combinations. Texture profile analysis (TPA) of the treated samples was performed to determine changes in muscle hardness. At treatment temperatures up to and including 50 °C, heat and pressure acted synergistically to increase muscle hardness. However, at 60 and 70 °C, hardness decreased following treatments in excess of 200 MPa. TPA was performed on extracted myofibrillar protein gels that after treatment under similar conditions revealed similar effects of heat and pressure. Differential scanning calorimetry analysis of whole muscle samples revealed that at ambient pressure the unfolding of myosin was completed at 60 °C, unlike actin, which completely denatured only above 70 °C. With simultaneous pressure treatment at >200 MPa, myosin and actin unfolded at 20 °C. Unfolding of myosin and actin could be induced in extracted myofibrillar protein with simultaneous treatment at 200 MPa and 40 °C. Electrophoretic analysis indicated high pressure/temperature regimens induced disulfide bonding between myosin chains.

Chromatographic alignment of LC-MS and LC-MS/MS datasets by genetic algorithm feature extraction

Liquid chromatography-mass spectrometry (LC-MS) datasets can be compared or combined following chromatographic alignment. Here we describe a simple solution to the specific problem of aligning one LC-MS dataset and one LC-MS/MS dataset, acquired on separate instruments from an enzymatic digest of a protein mixture, using feature extraction and a genetic algorithm. First, the LC-MS dataset is searched within a few ppm of the calculated theoretical masses of peptides confidently identified by LC-MS/MS. A piecewise linear function is then fitted to these matched peptides using a genetic algorithm with a fitness function that is insensitive to incorrect matches but sufficiently flexible to adapt to the discrete shifts common when comparing LC datasets. We demonstrate the utility of this method by aligning ion trap LC-MS/MS data with accurate LC-MS data from an FTICR mass spectrometer and show how hybrid datasets can improve peptide and protein identification by combining the speed of the ion trap with the mass accuracy of the FTICR, similar to using a hybrid ion trap-FTICR instrument. We also show that the high resolving power of FTICR can improve precision and linear dynamic range in quantitative proteomics. The alignment software, msalign, is freely available as open source.

Texture development in ice cream - effects of adding cocoa solids

Dynamic rheological techniques can aid the understanding of the factors contributing to ice cream structure, though the data obtained differs from that deduced from destructive techniques. Studies have shown that ice cream systems are both strain- and frequency-dependent. Chocolate ice cream is normally more viscous than the equivalent vanilla ice cream during mix preparation and has more body on freezing. Ice creams were prepared with and without cocoa solids and frequency sweeps were made from 0.1 to 100 Hz at 0.1% strain. With rapidly frozen ice creams, both G' and G" increased in the presence of cocoa solids. Comparison of mixes made with and without low-fat cocoa powder or non-gelatinizing starch demonstrated a similar relationship, with higher apparent viscosities in those mixes containing either cocoa powder or the starch. The results were consistent with the cocoa particles adding to the effect of the fat globules in increasing viscosity.

Composition, color, and texture of Thai indigenous and broiler chicken muscles

Chemical compositions and physical properties of mixed-sex Thai indigenous (Gallus domesticus) and broiler (commercial breed, CP707) chicken biceps femoris and pectoralis muscles were determined. Indigenous chicken muscles contained higher protein contents but lower fat and ash contents compared to broiler muscles (P < 0.001). The amino acid profile of the indigenous chicken muscles was similar to that of the broiler muscles except they were slightly richer in glutamic acid (P < 0.05). The indigenous chicken muscles contained more saturated and less polyunsaturated fatty acids than the broiler muscles. There were no differences in the monounsaturated fatty acid contents between the breeds. The total collagen contents of indigenous pectoralis and biceps femoris muscles were 5.09 and 12.85 mg/g, respectively, which were higher than those found in broiler pectoralis (3.86 mg/g) and biceps femoris muscles (8.70 mg/g) (P < 0.001). Soluble collagen contents were lower for indigenous pectoralis and biceps femoris muscles, 22.16 vs. 31.38% and 26.06 vs. 33.87%, respectively. The CIE system values of lightness (L*), redness (a*), and yellowness (b*) of indigenous chicken muscles were higher than those of broiler muscles. The shear values of indigenous chicken muscles either raw or cooked were higher than those of broiler muscles (P < 0.05). After cooking, the shear values decreased for broiler biceps femoris and pectoralis muscles (P < 0.05), whereas no change was observed for indigenous chicken biceps femoris muscle (P > 0.05). Shear values increased for indigenous chicken pectoralis muscle (P < 0.05).

Effect of heat treatment on changes in texture, structure and properties of Thai indigenous chicken muscle

Changes in texture, microstructure, colour and protein solubility of Thai indigenous and broiler chicken Pectoralis muscle stripes cooked at different temperatures were evaluated. The change in shear value of both chicken muscles was a significant increase from 50 to 80 degrees C but no change from 80 to 100 degrees C. A significant decrease in fibre diameter was obtained in samples heated to an internal temperature of 60 degrees C and the greatest shrinkage of sarcomeres was observed with internal temperatures of 70-100 and 80-100 C for broiler and indigenous chicken muscles, respectively (P < 0.05). Cooking losses of indigenous chicken muscles increased markedly in the temperature range 80-100 C and were significantly higher than those of the broiler (P < 0.001). With increasing temperature, from 50 to 70 degrees C, cooked chicken muscle became lighter and yellower. Relationships between changes in sarcomere length, fibre diameter, shear value, cooking loss and solubility of muscle proteins were evaluated. It was found that the solubility of muscle protein was very highly correlated with the texture of cooked broiler muscle while sarcomere length changes and collagen solubility were important factors influencing the cooking loss and texture of cooked indigenous chicken muscle. (c) 2004 Elsevier Ltd. All rights reserved.