Abstract
Amplitudemodulation frequencymodulation (AMFM) decompositions represent images using spatiallyvarying sinusoidal waves and their spatiallyvarying amplitudes. AMFM decompositions use different scales and bandpass filters to extract the wide range of instantaneous frequencies and instantaneous amplitude components that may be present in an image. In the past few years, as the understanding of its theory advanced, AMFM decompositions have been applied in a series of medical imaging problems ranging from ultrasound to retinal image analysis, yielding excellent results. This article summarizes the theory of AMFM decompositions and related medical imaging applications.
Keywords:
multidimensional AMFM methods; digital image processing; medical imaging1. Introduction
In the field of computer aided detection and diagnostics (CAD), recent advances in image processing techniques have brought a wide array of applications into the field. Many existing CAD methods rely on fixed basis functions based on wavelet decompositions [1] and Gabor filters [2]. AmplitudeModulation FrequencyModulation (AMFM) methods [36] represent an emerging technique that shows great promise in this area.
Multidimensional AMFM models and methods provide us with powerful, image and video decompositions that can effectively describe nonstationary content. They represent an extension to standard Fourier analysis, allowing both the amplitude and the phase functions to vary spatially over the support of the image, following changes in local texture and brightness.
To explain some of the advantages of AMFM methods, we begin with the basic AMFM model. In the 2D model, we expand an input image I(x, y) into a sum of AMFM harmonics using:
where a_{n}(x, y) denote slowlyvarying instantaneous amplitude (IA) functions, φ_{n}(x, y) denote the instantaneous phase (IP) components, and n = 1, 2,..., M indexes the different AMFM harmonics. In (1), the nth AMFM harmonic is represented by a_{n}(x, y) cos φ_{n}(x, y).
With each phase function, the instantaneous frequency (IF) vector field is defined by . Here, the AMFM demodulation problem is defined as one of determining the IA, IP, and IF functions for any given input image.
AMFM decompositions provide physically meaningful texture measurements. Significant texture variations are captured in the frequency components. For single component cases, IF vectors are orthogonal to equiintensity lines of an image, while the IF magnitude provides a measure of local frequency content. In (1), by using AMFM components from different scales, we can produce IF vectors from different scales, at a pixellevel resolution [4,5].
Since AMFM texture features are provided at a pixellevel resolution, AMFM models can be used to segment texture images that are difficult to model with the standard brightnessbased methods [7]. On the other hand, using just histograms of the IF and IA, we can design effective contentbased image retrieval systems using very short image feature vectors [8,9].
In summary, the advantages of AMFM methods include [10]: (i) they provide a large number of physically meaningful texture features, over multiple scales, at a pixellevel resolution, (ii) the image can be reconstructed from the AMFM decompositions, (iii) based on the target application, different AMFM decompositions using different frequency coverage can be designed, and (iv) very robust methods for AMFM demodulation have been recently developed (see some recent examples in [4]).
AMFM decompositions can be extended to represent videos using:
where each AMFM function has been extended to be a function of both space and time. The original phasebased modeling approach was provided in [11] and was recently extended and improved in [12,13].
AMFM decompositions can also be used to reconstruct the input images, allowing us to evaluate their effectiveness on different parts of the image. For continuousspace image decompositions, AMFM reconstruction examples can be found in [14], while [4,5] give several recent, robust multiscale examples for both images and videos. AMFM transform examples were shown in [15], while multidimensional orthogonal FM transforms were demonstrated in [16].
An early example of the use of frequencydomain filtering to target a particular application can be found in the fingerprint examples in [17]. More recently, [18] provided a treegrowth application, where interring spacing was used to design filterbanks that cover a specific part of the spectrum, so as to recover tree ring and tree growth structure from very noisy image inputs. For general images, Gabor filterbank approaches were investigated in [3]. Similarly, for general images and videos, dyadic, multiscale decompositions were introduced in [4].
A summary of AMFM methods is provided in Section 2. Results are presented for several medical imaging applications in Section 3. Finally, conclusions and future work are presented in Section 4.
2. Methods
2.1. Multidimensional AMFM demodulation methods
Recent interest in the development of multidimensional AMFM methods can be traced to early work on speech signal models by Jim Kaiser. In early work, Dr. Kaiser developed on algorithms for estimating the energy of ID signals in [1921]. This led to the development of the ID and 2D energy separation algorithm (ESA) as described in [2224]. Maragos et al. continued with the ID AMFM work previously mentioned for ID application in [22,23,25,26]. Early work on multidimensional energy separation methods appeared in [27]. AMFM demodulation based on Gabor wavelets appeared in [28,29]. Early research on the use of multidimensional energy operators continued in [24,30]. In what follows, we begin with an introduction to multidimensional energy operators.
In 2D, image energy is estimated using the TeagerKaiser operator given by [30]:
where ∇^{2 }= ∂^{2}/∂x^{2 }+ ∂^{2}/∂y^{2 }denotes the Laplacian operator. The IA and IF estimates are obtained using:
where , is estimated by replacing the xderivative by the yderivative in (4). To eliminate sign ambiguities, the TeagerKaiser operator can select the candidate IF for which the image gives the largest projection (e.g., see [31]). Recently, Kokkinos et al. showed a related accurate demodulation method using energy operators in [32] by computing all necessary derivatives by convolving with derivatives of Gabor filters, as opposed to using finite differences. The IA estimates are corrected by dividing by the magnitude response of the Gabor filter at the estimated IF (see [32] for details).
Analytic image methods for AMFM demodulation are based on providing a Hilbertbased extension of the ID Hilbertbased demodulation approach. Here, the basic idea is to simply apply the ID Hilbert operator along the rows (or the columns). The fundamental advantage of this approach is that it preserves the 2D phase and magnitude spectra of the 2D input image. In fact, implementation involves taking the 2D FFT of the input image, removing spectral frequency with a negative rowfrequency component, multiplying the result by 2, and taking the inverse 2D FFT. Given the conjugate symmetry of 2D images, the removal of two frequency quadrants does not result in the loss of any spectral information. Furthermore, it can be shown that for singlecomponent AMFM signals, this can lead to exact demodulation. In practice though, we replace derivatives by finite differences. We will further elaborate on this method in Section 2.2. For early work on this approach we refer to Havlicek's dissertation [3]. Havlicek et al. [3,3341] presented the first results for TVdimensional signals using the quasieigenfunction approximation (QEA) method (see Section 2.2).
Both ESA and Hilbertbased methods share the use of a filterbank prior to AMFM demodulation. The basic idea is to use a filterbank to be able to separate out among different AMFM components. AMFM demodulation is then applied at the output of each channel filter. Here, it is important to note that in the event that two AMFM components fall within the same channel filter, the filterbank approach will not allow us to separate them. Here, new algebraic approaches should be considered.
The filterbank generates AMFM demodulation outputs for each channel. Both ESA and Hilbertbased methods select estimates from a dominant component. For ESA, the dominant component is selected based on an energy criterion. In QEA, the dominant component is often selected based on the maximum IA estimate. Here, please note that a single channel is selected over the entire filterbank. In a multiscale approach, instead of selecting dominant components over the entire filterbank, we select the dominant channel from a collection of channels. The basic idea is to define scales based on the frequency magnitude. The most popular approach is to define low, medium, and high frequency scales (see examples in [4]).
Another robust approach for computing AMFM estimates based on a quasilocal method was developed in [4244] for 1D signals. This methodology was extended to digital images in [5,45]. Furthermore, in [5], we have a comparison of ESA, Hilbertbased, and the quasilocal methods for a variety of 2D AMFM signals. From the comparisons, we note that the choice of the filterbank can have a dramatic effect on the estimation. Flat passbands tend to help with the IA estimation. On the other hand, Gaborbased filterbanks are easy to design and implement and they can perform very well on IF estimation.
Related with the Hilbert based approach, Larkin et al. [46,47] introduced the phase quadrature transform by isotropically extending the Hilbert transform in 2D. Here, instead of applying the Hilbert transform along the rows or the columns, the spiralphase quadrature transform is applied along all directions. It is equivalent to applying the Hilbert transform along each radial direction. Here, a single component AMFM signal is convolved with
with a frequency response given by
Thus, the quadraturephase transform does not alter the magnitude of the AMFM signal. On the other hand, the same cannot be said about the phase. The phase information is not longer preserved. By examining (7), we can see that this is especially problematic for high frequencies (H).
Alternatively, Felsberg and Sommer proposed an nD generalization of the 1D analytic signal based on the Riesz transform, which is used instead of the Hilbert transform [48,49]. This is termed the monogenic signal in [48]. Felsberg and Sommer introduced the following filters in the frequency domain [50]:
with spatial impulse responses given by
The local amplitude of the filter response corresponds to a quantitative measure of a structure (including the contrast) and the local phase corresponds to a qualitative measure of a structure (step, peak, etc.) [50]. The monogenic signal f_{M }is then defined as a 3D vector formed by the signal f(x, y) with its Riesz transform with (h) = (h, h_{2}), using
The phase and magnitude of the monogenic signal f_{M }are then taken as the phase and magnitude of the AMFM signal. This is also extended to multiple scales in the multiresolution framework of [51], and into scalespace in [49]. Other related work on image demodulation based on the Riesz transform extension has been reported by Mellor, Noble, Hahn, Felsberg, Sommer and collaborators in [48,5257]. In (8), it is important to note that while the input signal is 2D, the generated monogenic signal is actually 3D. This was done to extend the ID analytic properties to 2D. On the other hand, it is also clear that these 2D convolutions will also alter the phase of the 2D input AMFM signal. In fact, as can be seen in the fingerprint example of Figure eleven of [51], the estimated amplitudes contains the ridges. In contrast, in the ESA fingerprint examples of [6], the ridges are modeled as a FrequencyModulation process. This is a fundamental difference in the different approaches considered here.
Similar to AMFM methods, we mention the work by Knutsson et al. for representing local structures on phase using tensors [5860]. Furthermore, the complex Wavelet transform provides an extension to the discrete wavelet transform that is related to the 2D Hilbertspace extension [6165].
2.2. Multiscale AMFM methods
In this section, we provide more details on the use of multiplescales in AMFM demodulation. Many of the concepts introduced in this section are shared by the AMFM demodulation methods described in Section 2.1.
We consider multiscale AMFM representations of images given by
where n = 1, 2, ..., M denote different scales [4,5]. In (9), a continuous image I(x, y) is a function of a vector of spatial coordinates (x, y). A collection of M different scales are used to model essential signal modulation structure. The amplitude functions a_{n}(·) are always assumed to be positive.
AMFM models nonstationary image content in terms of its amplitude and phase functions [6]. The aim is to let the frequencymodulated (FM) components cos φ_{n}(·) capture fastchanging spatial variability in the image intensity. For each phase function φ_{n}(·) we define the IF, ∇φ_{n}(·). in terms of the gradient:
The IF vector ∇φ_{n}(·) can vary continuously over the spatial domain of the input signal.
We generalize the concept of scale by considering AMFM demodulation over a collection of bandpass filters. This is clearly depicted in Figure 1. For realvalued images, we only need two quadrants as shown in Figure 2. Note that by applying the ID Hilbert operator along the columns, the upper two frequency quadrants will be set to zero. On the other hand, if the monogenic signal is computed here, we will keep the entire frequency spectrum and generate a 3D signal as given in (8).
Figure 1. Block diagram of the 2D multiscale AMFM demodulation. Dominant AMFM components are selected over different image scales. The bandpass filter selector (upper left) is used to define the bandpass filters that correspond to each scale. The dominant AMFM component is selected using the maximum amplitude at every pixel (see lowest block).
Figure 2. Fourscale filterbank used. (a) Complete frequency spectrum of the filterbank. (b) Zoom on the low frequency bandpass filters.
The effect of either operation is to generate a complexvalued signal estimates of the form: a(x, y) exp (jφ(x, y)), where the IA and phase functions will hopefully approximate the input AMFM components. Then, a collection of bandpass filters is used to isolate the individual AMFM components [3].
The basic assumption here is that different AMFM components will be picked up by different bandpass filters at the same image region. In other words, given any local image region the assumption is that the corresponding AMFM components will be separated by the use of different filters in the interbank.
In discrete terms (see Figure 1), given the digital input image I(k_{1}, k_{2}) (where k_{1 }and k_{2 }are the discrete versions of x and y, respectively), the application of an analytic extension generates an "analytic signal" of the form: I_{AS }(k_{1}, k_{2}). This 2D extension is processed through a collection of bandpass filters (to be described in the Section 2.3) within the desired scale. Each processing block will produce the IA, the IP, and the IF in both x and y directions. At the output of each block, the IA and the IP are estimated by simply taking the magnitude and the phase of the output.
In dominant component analysis (DCA), for each pixel approach, the best AMFM demodulation estimates are selected from the bandpass filter that produces the largest IA estimate. For energybased approaches, the dominant channel is selected based on the maximum energy estimate. Hence, the algorithm adaptively selects the estimates from the bandpass filter with the maximum response. This approach does not assume spatial continuity, and allows the model to quickly adapt to singularities in the input signal (both high and low frequency changes).
A robust method for estimating the IF was produced by Murray et al. [4,5]. Here, for estimating the first IF component for lowfrequencies, four estimates are generated using:
where , and n_{1 }represents a variable displacement that varies from 1 to 4. Variable spacing, local linear phase (VSLLP) produces the most accurate of the four estimates by considering the condition number of the arccos(·) function [4,5]. Similarly, the approach is repeated for the second direction. It turns out that VSLLP reduces to the earlier method of QEA for the special case when n_{1 }= n_{2 }= 1 (see [3] for QEA).
2.3. Filterbank design
For computing the AMFM estimates, the AMFM components need to be isolated from (9) using a multiscale filterbank. Here, the basic idea is to isolate different AMFM components over different frequencies (see [66] for details). Figure 2 (a) depicts the frequency support for a dyadic filterbank decomposition.
In Figure 2, the low pass filter (LPF, with label number 1) has frequency support in [π/16, π/16] radians for both the x and y directions. The filters in the highest frequencies (filters from 2 to 7 in Figure 2), have a bandwidth of π/2 for both x and y directions. The bandwidth is decreased by a factor of 0.5 for each added scale. In Figure 2 (b), a closeup that shows the low frequency filters is provided.
The scales are provided using a collection of bandpass filters as: (i) LPF, (ii) very low frequencies (VL), (iii) low frequencies (L), (iv) medium frequencies (M) and (v) H as given in [4], When the structures are well defined in sizes, the measurements can be related with specific instantaneous wavelengths described in number of pixels [67,68].
For example, for classifying lung opacities in pneumoconiosis (to be fully described in Section 3.3), the correspondence between opacity size and instantaneous frequencies using eleven different combinations of scales (CoS) is shown in Table 1. The pneumoconiosis images that were analyzed in the study have a resolution of 300 dpi. This resolution means that along 0.5 mm we have 6 pixel samples. From the ILO standard [69], the rounded opacities labeled as p, q and r, are in the sizes up to 10 mm. With the knowledge of the sizes, the lengths can be related with the corresponding frequencies. These relationships are summarized in Table 2. Based on the frequency ranges for each rounded opacity type, the cutoff frequencies of the bandpass filters are set up to fit the opacities in them by each scale. Then, a fourscale filterbank is designed such that each scale will be related with the frequency ranges described in Table 2.
Table 1. CoS used for computing the dominant AMFMfeature parameters
Table 2. Relationship among ILO standard grades for rounded opacities and the size in mm, size in pixels and range in frequency
2.4. Postprocessing methods
2.4.1. Histogram features
AMFM features can be summarized using appropriately chosen histograms. Here, we provide a summary of how this can be accomplished.
For each of the 11 CoS (see Table 1), a 96bin feature vector that contains the AMFM histograms is created. Each AMFM histogram has the information of (32bin each): (i) the IA, (ii) the IF magnitude, and (iii) the IF angle. Thus, each image produces 11 different feature vectors, each one corresponding to one of the 11 different CoS.
Since the AMFM estimates of neighboring pixels could be affected by noisy estimates, a 9 × 9 median filter is often applied to the output of the IA and IF estimates (see [4,5,14]). We refer to [70] for methods that can be used to adaptively set the median filter window size.
Due to the ambiguity of the IF vectors (coming from cos φ(k_{1}, k_{2}) = cos(φ(k_{1}, k_{2}))), the IF estimates are mapped to two frequency quadrants:
• Instantaneous frequencies and are kept as they are.
• Instantaneous frequencies and are kept as they are.
• Instantaneous frequencies and are mapped to and .
• Instantaneous frequencies and are mapped to and .
Then, the IF histograms are computed.
For the IF angle, given by arctan , the histogram is centered around the maximum value using:
[h] ← histogram [IF angle using l bins] {Where h is a column vector.}
H ← [hhh]. {Triplicate version of h.}
i ← Location of the maximum value ofh.
h ←H_{l+il/2+l:l+i+l/2}.
As last step, all the computed histograms (IA, IF magnitude, and IF angle) are normalized such that the cumulative sum of each one is equal to 1. Thus, the bins will now correspond to the probability density function (PDF) estimates.
2.5. Partial least squares classification of AMFM features
The histograms of AMFM features can be used in classification applications. Here, we introduce a promising classifier based on Partial least squares.
Partial least squares (PLS, [71]) is a linear regression method. Here, classification is expressed as a regression problem is given by
where y is a n × 1 vector of the classification variables, X is a n × p matrix of the extracted AMFM features, β is a p × 1 vector of regression weights, and ε is a n × 1 vector of residuals. The least squares solution to estimating β is given by the normal equations β = (X^{T}X)^{l }(X^{T}y).
In most of medical imaging applications, there are much more features than images (p <n), and AMFM features in X are highly correlated. Thus, X^{T }X canbe singular or nearly singular and a unique solution to the normal equations will not exist. PLS reduces X to a lower dimensional subspace (k ≪ p). The first step is to factor X as X = TL, where T is an orthogonal n × k matrix of Tscores and L is a k × p matrix of factor loadings. The Tscores matrix are used to find a threshold for classification as outlined in [72],
2.6. Morphological segmentation
AMFM features can also be used for image segmentation. A simple Bayesian classification scheme based on the IA and IF magnitude was demonstrated in [7].
The segmented AMFM regions can be very noisy. Segmentation region denoising can be achieved via the use of an alternating sequential filter (ASF) as given in [73]. The noise in the segmented image is assumed to be present due to small structural components that do not belong to the abnormal regions. Under this assumption, it can be shown that the optimal morphological filter will be an ASF defined by:
where I_{s }is the binary segmented image, B is the structural element; nB is the result of n dilations of B by itself; ○ the standard open operation and • the standard close operation.
The value of n needs to be chosen so that the smaller, noisy structures are removed, while the larger structures are preserved. Once the order of the morphological filter is calculated, the segmented image is denoised using (13).
3. Applications in medical imaging
We provide a selected list of AMFM biomedical imaging applications in Table 3 (see [10] for an earlier version of this table). Here, note that there is a number of medical applications based on energy operators. More recently, medical imaging applications have benefited from the use of multiscale methods summarized in Section 2.2.
Table 3. Select list of AMFM applications in Biomedical Imaging
An interesting application of AMFM models in motion estimation was presented by Fleet and Jepson [11]. In [11], Fleet and Jepson use a 3D Gabor interbank to estimate optical flow motion at pixellevel resolution. Their approach motivated the later application of a harmonic phase (HARP) method in [74].
In one of the earlier applications in electron microcopy, a simple Bayesian method using the IA and the IF was used, in conjunction with morphological filtering to provide segmentations of different abnormalities over 26 images [7]. Elshinawy et al. [75,76] demonstrated the reconstruction of breast cancer images using AMFM components.
Boudraa et al. [77] introduced a new crossenergy operator and used the operator to demonstrate the functional segmentation of dynamic nuclear images. Maragos et al. provide an important application for ultrasound spectroscopy in [78], where AMFM models are used for improving Doppler ultrasound resolution. Vectorvalued based AMFM demodulation is given in [79] by Alexandratou et al. Multiscale AMFM methods were applied to chest radiographs [67,68], ultrasound images of the carotid artery [8083], image retrieval in ophthalmology by Acton et al. [84], retinal image classification [8,85,86], and electron microscopy [7].
We also have some related 1D medical signal applications by different research groups. Relevant 1D AMFM methods appear in [8790]. Medical applications include the classification of surface electromyographic signals in [80], and the analysis of brain rhythms in electroencephalograms [91].
More general (nonmedical) AMFM applications in tracking include the work reported by Prakash et al. [92,93] and Mould et al. [94]. Recent AMFM texture analysis research is also reported by Kokkinos et al. [95] and Tay et al. [96]. Recently, Belaid et al. [97] presented an ultrasound image segmentation method based on the monogenic signal and quadrature filters [49].
In the rest of this section, we provide recent examples from the application of multiscale AMFM decompositions to biomedical imaging. While the applications are focused on the particular models described in Section 2.2, it is important to note that new applications are expected, especially from methods based on multidimensional energy operators and methods based on the monogenic signal.
3.1. Optical flow motion estimation
One of the earliest applications of AMFM methods comes from the application of phasebased method in estimating opticalflow motion. Here, we use the term "opticalflow motion" to differentiate it from actual object motion. We begin with a formulation of optical flow motion and proceed to describe the phasebased approach.
Let the image intensity in a digital video sequence be denoted by I(x, y, t). After a sufficiently small time interval Δt, the intensity at (x, y) will move to a point (x + Δx, y + Δy). Here, we are assuming that the image intensity is preserved: I(x + Δx, y + Δy, t + Δt) = I(x, y, t).
A Taylorseries expansion at starting point (x, y) gives:
where h.o.t refers to higherorder terms. Next, if we ignore the higher order terms, divide by Δt, and then take the limit as Δt → 0, we get the Optical flow constraint equation (OFCE):
To formulate the motion problem, define the motion velocities u and v using: u(x, y, t) = ∂x/∂t and v(x, y, t) = ∂y/∂t. This leads to a new form of the OFCE given by:
In (14), at every pixel, we have a single equation in two unknowns. This forces us to add additional constraints. For example, a continuity constraint requires that the solution keeps the spatial integral of small over the solution small. With this additional constraint we end up with two equations for two unknowns per pixel.
Fleet and Jepson [11] proposed the use of a singlecomponent AMFM model for modeling digital video based on
For estimating the AMFM components, Fleet and Jepson used a 3D Gabor interbank. Then they estimate the IF using
where * denotes the complex conjugate, and ∇ refers to the spatialgradient. The IF is then used to estimate the velocity components that are along the direction of the the IF using
where the first parenthesis gives the velocity magnitude, while the second parenthesis provides the unit vector defined in the direction of the IF. To provide for a method of identifying accurate motion estimation, Fleet and Jepson require that the estimated IF is within the range of the estimating filter. This is expressed as
where (k_{i}, w_{i}) is the peak tuning frequency of the ith filter, σ_{k }is the standard deviation of the filter's amplitude spectrum, and τ is a threshold used to reject unreliable estimates of instantaneous frequencies. Similarly, for the amplitude, they require that the local signal amplitude must be as large as the average local amplitude, and at least 5% of the largest response amplitude across all the filters at that frame. When either one of these two conditions is not met, the method does not provide velocity estimates.
To recover the full velocity components, Fleet and Jepson apply a local linear fit to the projected velocities. This is accomplished by fitting a local linear model over 5 × 5 neighborhoods:
where the local velocity estimate at the central pixel is given by (α_{0}, β_{0}).
This study represents an important early application of an AMFM model to motion estimation. It provided the motivation to the recent work in estimating cardiac motion from MRI (see HARP method in [74]). We expect to see the emergence of new, multiscale AMFM methods that would be applied to motion estimation from medical video.
3.2. Retinal image analysis
Diabetic retinopathy (DR) and agerelated macular degeneration (AMD) are two retinal diseases that present particular characteristics (lesions) on retinal photographs which can be used in early detection and/or classification of the disease. DR and AMD affect individuals in their most productive years of life. Early detection that leads to prevention of vision loss, alone, will lead to significant decrease in risk of early vision loss.
According to the National Eye Institute, DR is one of the leading causes of blindness among workingage Americans, while macular degeneration is a leading cause of blindness among older Americans [98]. It has been shown that regular comprehensive eye exams and timely treatment can lead to improved outcomes and reduced loss of vision. However, to screen the tens of millions at risk for DR would tax the healthcare system beyond capacity. Results of an automatic DR screening system based on AMFM exhibit strong promise in addressing this problem.
AMD is the most common cause of visual loss in the United States and is a growing public health problem. One third of Americans will develop AMD in their lifetimes. To detect AMD, retinal images are graded using the agerelated eye disease study (AREDS) protocol for human grading. AMFM has been used to develop an automated system for characterizing pathological features on these images.
3.2.1. Diabetic retinopathy
Multiscale AMFM decompositions have been used to characterize retinal images with DR. and to screen out patients with different levels of the disease. Figure 3 shows an image of a retina with DR and one of its IA image.
Figure 3. Example of a retinal image with DR. Note the exudates shown with a white square. (a) Original Retinal Image. (b) IA using CoS 4. (c) Thresholded Image of (b) showing large IA values. A valid range for the pixels between [1530%] of the maximum value was used for determining the thresholded image.
The types of lesions that need to be detected by the system are shown in Figure 4. Agurto et al. [9,85,99] presented the extraction of AMFM features for classification of DR images using 376 images from the MESSIDOR database [100].
Figure 4. Lesions typical of a DR retinal image. A) Microaneurysm B) Exudates C) Hemorrahages D) Neovascularization.
These 376 images were processed with AMFM and were divided in two groups for training and testing purposes. Table 4 shows the number of images per risk in each group. The most severe level of DR is grade 3, in which 15 or more microaneurysms or 5 or more hemorrhages are present. For grade 2, the number of microaneurysms in the retina is between 5 and 15 and for grade 1 the retinal images contain between 1 and 5 microaneurysms. The grade 0 corresponds to diabetic patients with no apparent DR in the image. It is important to mention that the images not only contain red lesions such as microaneurysms and hemorrhages, but also contain bright lesions such as exudates and cottonwool spots.
Table 4. Images correctly classified per risk
After all the images were processed with AMFM, they were divided in regions of 100 × 100 pixels avoiding the optic disc. Feature vectors were calculated using the moments in each region: mean, standard deviation, skewness, and kurtosis. Finally, the information per image is the input for the PLS algorithm were the algorithm is trained to classify normal images vs abnormal images.
The results obtained with this approach show an area under the ROC curve of 0.84 with best sensitivity/specificity of 92%/54%. Table 4 (last column) shows the number of images correctly classified per risk level [9]. On this table we can see a large percentage of the highrisk patients are being classified correctly.
3.2.2. Agerelated macular degeneration
Barriga et al. [86] have developed a system based on AMFM to generate multiscale features for classifying pathological structures, such as drusen, on a retinal image. Drusen are subretinal deposits characteristic of AMD. Figure 5 shows an AMD patient with soft and hard drusen present (white spots on the fundus). The standard highresolution AREDS images were downloaded from the Wisconsin Fundus Photo Reading Center's website [101]. A certified ophthalmic technician (reader) selected the retinal features. For these studies, only areas that contain retinal background, vessels, soft drusen, and hard drusen were used. After the reader selected the features, regions of interest (ROI) of 40 × 40 pixels were extracted from the image.
Figure 5. Retinal image from AMD patient from the AREDS standard database.
After calculating the AMFM histograms, principal components analysis (PCA) was applied to the matrix containing histograms of the 120 ROIs (30 for each retinal structure) for three different AMFM features: IA, IF magnitude (IF), and IF angle. Finally, the Mahalanobis distances between the different retinal structures based on PCA decomposition were calculated. These distances are shown in Table 5.
Table 5. Mahalanobis distance between retinal features: Retinal Background (RB), Hard Drusen (DRH), Soft Drusen (DRS), and Vessels
The numbers in Table 5 represent the standard deviations separating the histograms of the retinal structures. The entries under hard drusen (DRH) and soft drusen (DRS) show that they are significantly differentiated from other structures in the retina. A distance of 3 standard deviations represents a classification accuracy of more than 90%. The most interesting entry is the Mahalanobis distance between the two drusen types (2.8 standard deviations). Though still considered a high classification rate (85%), this demonstrates the challenge not only to the algorithm, but for the grader in unequivocally assigning drusen to one class or another.
3.3. Pneumoconiosis in chest radiographs
The chest radiograph is an essential tool used in the screening, surveillance, and diagnosis of dustrelated respiratory illness resulting from silica and coal dust, asbestos, and a variety of other dusts that can lead to disease. The standardized method used by the international labor organization (ILO, [69]) for interpreting the chest radiograph for inorganic dustinduced diseases or pneumoconioses has been widely utilized for the past five decades.
Two limitations of the current ILO system are the intra and interinterpreter variability and the time consuming process of interpreting large numbers of radiographs taken for screening and surveillance programs. Murray et al. [67,68] have developed a technique based on AMFM that detects the level of pneumoconioses based on the nodule formation on the lungs.
Currently, the US remains 6th in the world for pneumoconiosis and interstitial lung disease (ILD) or pulmonary interstitial fibrosis as large numbers of workers continue to be exposed to dust in their work environment. Coal workers' pneumoconiosis (CWP) is identified by a specific pattern of XRay abnormalities and a history of exposure to coal dust. The chest radiograph is the single most useful tool for clinically evaluating both occupationally related and nonoccupational chronic lung diseases. The chest radiograph is an essential tool used in the screening, surveillance, and diagnosis of dustrelated respiratory illness resulting from silica and coal dust, asbestos, and a variety of other dusts that can lead to disease.
An AMFM method for grading chest radiographs according to the ILO standards is presented in [67,68] by Murray et al. Related work, by different research groups, can be found in [102108]. First, a logarithmic transformation to the images is applied to improve the contrast of the XRays [109]. Then, the AMFM estimates of each XRay image is encoded using their histograms as described in the Section 2.4.1 using the tuned bandpass filters and scales showed in Tables 1 and 2. The final classification was computed using the linear regression method PLS (described in Section 2.5).
The results obtained by the system were excellent (area under the ROC curve, AUC = 1.0). However, the significance of the results should be debated based on the small database used. The authors plan as future work the testing of the system using a much bigger database. The authors produced classification results that were significantly better than those presented in [102,103,107].
3.4. Carotid artery ultrasound images
Cardiovascular disease (CVD) is the third leading cause of death and adult disability in the industrial world after heart attack and cancer. Of all the deaths caused by CVD among adults aged 20 and older, an estimated 6 millions are attributed to coronary heart disease and to stroke, with atherosclerosis as the underlying cause [8082,110].
A method based on image analysis of ultrasound images of carotid plaques that can differentiate between the stable plaques that tend to remain asymptomatic and the unstable ones that eventually produce symptoms has the potential to refine the basis for surgery and spare some patients from an unnecessary costly operation which itself carries a risk of stroke.
AMFM methods are being applied for characterizing and analyzing plaques in ultrasound images. Christodoulou et al. present in [81] their investigations for the AMFM characterization of carotid plaques in ultrasound images. In [82], Loizou et al. present how to use AMFM features for describing atherosclerotic plaque features. In what follows, we described the basics of their approaches.
The extraction of features characterizing efficiently the structure of ultrasound carotid plaques is important for the identification of individuals with asymptomatic carotid stenosis at risk of stroke. Christodoloulou et al. present how to characterize the carotid plaques using AMFM features in [81]. They use as basic descriptors the information about the AMFM estimates similar as in the Section 2.4.1: the IA, IF magnitude and IF angle histograms. In addition to the AMFM features, the authors compute Spatial Gray Level Dependence Matrices (SGLDM) computing the following texture measures: (i) angular second moment, (ii) contrast, (iii) correlation, (iv) sum of squares: variance, (v) inverse difference moment, (vi) sum average, (vii) sum variance, (viii) sum entropy, (ix) entropy, (x) difference variance, (xi) difference entropy and (xii) information measures of correlation. Also, the Gray Level Difference Statistics (GLDS) are computed: (i) contrast, (ii) angular second moment, (iii) entropy and (iv) mean. The final step applied by the authors corresponds to the statistical knearest neighbor (KNN) classifier implemented for different values of k (k = 1, 3, 5, 7, 9, 11, 13). The authors applied their method to a database of 274 carotid plaque ultrasound images divided as 137 symptomatic and 137 asymptomatic. Using the AMFM estimates, the authors got a classification success up to 71.5% when the three AMFM estimates and k = 5 were used. This result was better than using the textures features SGLDM and GLDS that gave 68.2%.
Loizou et al. [82] present in an AMFM analysis as an application for investigating the intima media complex (IMC), media layer (ML) and intima layer (IL) of the common carotid artery (CCA). This study represents the first study for the IMC, IL and ML. Clinically, the intimamedia thickness (IMT) is used as a validated measure for the assessment of atherosclerosis, that causes enlargement of the arteries and thickening of the artery walls. The authors use the same AMFM histograms as in Section 2.4.1 but computed using only horizontaloriented bandpass filters. Thus, the authors use similar scales like those described in Table 1 but using only the filters numbered as (see Figure 2a): 1 (LPF), 2, 5 (in the H scale), 8, 11 (in the M scale) and 14 and 17 (in the L scale). The reason for the selection of these horizontaloriented filters is related with the horizontally elongated nature of the atherosclerotic structures. As the last step for the statistical analysis, the authors use the MannWhitney rank sum test in order to identify is there are significant differences (SD) or not (NS) between the extracted AMFM features. The authors use p < 0.05 for significant differences and comparison between different age groups.
Their study, performed on 100 ultrasound images, reveals that the IA of the media layer decreases with age. The authors state that the decreasing of the IA is maybe related to the reduction in calcified, stable plaque components and an increase in stroke risk with age. In terms of the IF, the median of this AMFM estimate of the media layer increases suggesting fragmentation of solid, large plaque components that also increases the risk of stroke.
4. Conclusions and future work
We have provided a summary of recent AMFM applications in medical imaging. In the coming years, we expect that there will be a variety of new applications. Here, it is important to note that medical imaging applications come from a variety of AMFM methods and models.
The advantages of considering AMFM models come from their ability to describe image structures over different frequency scales. For example, in retinal images, structures such as drusen and microaneurysms have a welldefined shape that get captured by the higherfrequency scale AMFM filters. On the other hand, vessels and hemorrhages are better represented in the lower frequencies. In chest radiographs, nodules can also be well represented using the AMFM features. The combination of feature extraction and classification using PLS has produced robust systems for the analysis of chest radiographs, where an automatic grading system for pneumoconiosis has been demonstrated on a limitedsize database. For DR, a screening system has been developed and it is currently being tested on an extensive database for FDA approval.
Future work on AMFM models will undoubtedly yield new methods and new applications. We expect to see new applications from atleast three separate approaches: (i) applications from AMFM methods associated with multidimensional energy operators, (ii) applications from multiscale AMFM decompositions based on Hilbertbased methods, (iii) applications associated with methods associated with the monogenic signal, and (iv) applications from new AMFM methods that are currently being developed.
Competing interests
V. Murray, M.S. Pattichis, and P. Soliz are coinventors on a relevant, pending patent on the use of AMFM methods in medical imaging applications.
V. Murray and M.S. Pattichis declare no competing interests in the research.
E.S. Barriga and P. Soliz work with VisionQuest Biomedical, LLC. who is interested in the commercialization of AMFM technology in medical image analysis applications. Publication of this manuscript represents positive publicity to VisionQuests commercial endeavours.
Authors' contributions
VM worked on the computational examples, worked on all aspects of the writing, and coordinated the writing efforts of the paper. MSP worked on the overall structure and organization of the document, helped with the literature review section, and edited the technical sections of the paper. ESB worked on the literature review of medical imaging applications, in particular for diabetic retinopathy and agerelated macular degeneration. PS provided technical guidance on the AMFM applications reported in the manuscript, including diabetic retinopathy, agerelated macular degeneration, and pneumoconiosis.
References

M Unser, A Aldroubi, A Laine, Guest editorial: Wavelets in medical imaging. IEEE Trans Med Imag 22, 285–288 (2003). Publisher Full Text

A Osareh, B Shadgar, Retinal vessel extraction using gabor filters and support vector machines. Commun Comput Inf Sci 6, 356–363 (2009). Publisher Full Text

JP Havlicek, AMFM image models (The University of Texas at Austin, 1996) (Ph, 1996), . D. thesis

VM Murray Herrera, AMFM methods for image and video processing (University of New Mexico, 2008) (Ph, 2008), . D. thesis

V Murray, P Rodriguez, MS Pattichis, Multiscale AMFM demodulation and image reconstruction methods with improved accuracy. IEEE Trans Image Process 19(5), 1138–1152 (2010). PubMed Abstract  Publisher Full Text

MS Pattichis, AC Bovik, Analyzing image structure by multidimensional frequency modulation. IEEE Trans Pattern Anal Mach Intell 29(5), 753–766 (2007). PubMed Abstract  Publisher Full Text

M Pattichis, C Pattichis, M Avraam, A Bovik, K Kyriacou, AMFM texture segmentation in electron microscopic muscle imaging. IEEE Trans Med Imag 19(12), 1253–1257 (2000). Publisher Full Text

V Murray, M Pattichis, P Soliz, New AMFM analysis methods for retinal image characterization. Asilomar Conference on Signals, Systems and Computers (2008)

C Agurto, V Murray, E Barriga, S Murillo, M Pattichis, H Davis, S Russell, M Abramoff, P Soliz, Multiscale AMFM methods for diabetic retinopathy lesion detection. IEEE Trans Med Imag 29(2), 502–512 (2010)

M Pattichis, Multidimensional AMFM models and methods for biomedical image computing. the 34th IEEE Annual International Conference of the Engineering in Medicine and Biology Society (2009)

DJ Fleet, AD Jepson, Computation of component image velocity from local phase information. Int J Comput Vision 5(1), 77–104 (1990). Publisher Full Text

V Murray, MS Pattichis, AMFM demodulation methods for reconstruction, analysis and motion estimation in video signals. IEEE Southwest Symposium on Image Analysis and Interpretation (IEEE Computer Society, Los Alamitos, CA, USA, 2008) 0, pp. 17–20

V Murray, SE Murillo, MS Pattichis, CP Loizou, CS Pattichis, E Kyriacou, A Nicolaides, An AMFM model for motion estimation in atherosclerotic plaque videos. 41st Asilomar Conference on Signals, Systems and Computers, 746–750 (2007)

J Havlicek, P Tay, A Bovik, Handbook of Image and Video Processing (Elsevier Academic Press, Ch. AMFM Image Models: Fundamental Techniques and Emerging Trends, 2005), pp. 377–395

M Pattichis, A Bovik, AMFM expansions for images. Proc European Signal Processing Conference (1996)

M Pattichis, A Bovik, J Havlicek, N Sidiropoulos, Multidimensional orthogonal fm transforms. IEEE Trans Image Process 10(3), 448–464 (2001). PubMed Abstract  Publisher Full Text

M Pattichis, G Panayi, A Bovik, H ShunPin, Fingerprint Classification using an AMFM model. IEEE Trans Image Process 10(6), 951–954 (2001). Publisher Full Text

J Ramachandran, Image analysis of wood core using instantaneous wavelength and frequency modulation (University of New Mexico, 2008) (Ph, 2008), . D. thesis

J Kaiser, On a simple algorithm to calculate the 'energy' of a signal. International Conference on Acoustics, Speech, and Signal Processing 1, 381–384 (1990)

J Kaiser, On teager's energy algorithm and its generalization to continuous signals. IEEE Digital Signal Processing Workshop (New Paltz, NY, 1990)

J Kaiser, Some useful properties of teager's energy operators. IEEE International Conference on Acoustics, Speech, and Signal Processing 3, 149–152 (1993)

P Maragos, JF Kaiser, TF Quatieri, On amplitude and frequency demodulation using energy operators. IEEE Trans Signal Process 41(4), 1532–1550 (1993). Publisher Full Text

P Maragos, JF Kaiser, TF Quatieri, Energy separation in signal modulations with applications to speech analysis. IEEE Trans Signal Process 41(10), 3024–3051 (1993). Publisher Full Text

P Maragos, AC Bovik, Image demodulation using multidimensional energy separation. J Opt Soc Am A 12(9), 1867–1876 (1995). Publisher Full Text

P Maragos, T Quatieri, J Kaiser, Speech nonlinearities, modulations, and energy operators. International Conference on Acoustics, Speech, and Signal Processing 1, 421–424 (1991)

P Maragos, J Kaiser, T Quatieri, On separating amplitude from frequency modulations using energy operators. IEEE International Conference on Acoustics, Speech, and Signal Processing 2, 1–4 (1992)

P Maragos, A Bovik, T Quatieri, A multidimensional energy operator for image processing. SPIE Symp Visual Commun Image Processing (1992)

A Bovik, N Gopal, T Emmoth, A Restrepo, Localized measurement of emergent image frequencies by gabor wavelets. IEEE Trans Inf Theory 38(2), 691–712 (1992). Publisher Full Text

A Bovik, P Maragos, T Quatieri, Amfm energy detection and separation in noise using multiband energy operators. IEEE Trans Signal Process 41(12), 3245–3265 (1993). Publisher Full Text

P Maragos, A Bovik, Demodulation of images modeled by amplitudefrequency modulations using multidimensional energy separation. IEEE International Conference on Image Processing 3, 421–425 (1994)

MS Pattichis, AMFM transforms with applications (University of Texas at Austin, 1998) (Ph, 1998), . D. thesis

I Kokkinos, G Evangelopoulos, P Maragos, Texture analysis and segmentation using modulation features, generative models, and weighted curve evolution. IEEE Trans Pattern Anal Mach Intell 31(1), 142–157 (2009). PubMed Abstract  Publisher Full Text

J Havlicek, J Havlicek, A Bovik, The analytic image. Image Processing, 1997. Proceedings, International Conference on 2, 446–449 (1997)

J Havlicek, J Havlicek, N Mamuya, A Bovik, Skewed 2d hilbert transforms and computed amfm models. Image Processing, 1998. ICIP 98. Proceedings. 1998 International Conference on 1, 602–606 (1998)

J Havlicek, D Harding, A Bovik, Discrete quasieigenfunction approximation for amfm image analysis. Image Processing, 1996. Proceedings, International Conference on 1, 633–636 (1996)

J Havlicek, D Harding, A Bovik, Extracting essential modulated image structure. Signals, Systems and Computers, 1996. 1996 Conference Record of the Thirtieth Asilomar Conference on 2, 1014–1018 (1996)

J Havlicek, A Bovik, Multicomponent amfm image models and waveletbased demodulation with component tracking. Image Processing, 1994. Proceedings. ICIP94., IEEE International Conference 1, 41–45 (1994)

J Havlicek, A Bovik, P Maragos, Modulation models for image processing and waveletbased image demodulation. Signals, Systems and Computers, 1992. 1992 Conference Record of The TwentySixth Asilomar Conference on 2, 805–810 (1992)

J Havlicek, D Harding, A Bovik, The multicomponent amfm image representation. Image Process. IEEE Trans 5(6), 1094–1100 (1996)

J Havlicek, M Pattichis, D Harding, A Christofides, A Bovik, Amfm image analysis techniques. Image Analysis and Interpretation, 1996., Proceedings of the IEEE Southwest Symposium on, 195–200 (1996)

J Havlicek, D Harding, A Bovik, Reconstruction from the multicomponent amfm image representation. Image Processing, 1995. Proceedings, International Conference on 2, 280–283 (1995)

D Vakman, Signals Oscillations and Waves: A Modern Approach (Artech House, Boston, 1998)

D Vakman, On the Analytic Signal, the TeagerKaiser energy algorithm, and other methods for defining amplitude and frequency. IEEE Trans Signal Process 44(4), 791–797 (1996). Publisher Full Text

G Girolami, D Vakman, Instantaneous frequency estimation and measurement: a QuasiLocal Method. Measurement Sci Tech 13(6), 909–917 (2002). Publisher Full Text

P Rodriguez, M Pattichis, New algorithms for fast and accurate amfm demodulation of digital images. IEEE International Conference on Image Processing 2, II1294–7 (2005)

KG Larkin, DJ Bone, MA Oldfield, Natural demodulation of twodimensional fringe patterns. I. General background of the spiral phase quadrature transform. J Opt Soc Amer A: Opt Image Sci Vision 18(8), 1862–1870 (2001). Publisher Full Text

KG Larkin, Natural demodulation of twodimensional fringe patterns. II. Stationary phase analysis of the spiral phase quadrature transform. J Opt Soc Amer A: Opt Image Sci Vision 18(8), 1871–1881 (2001). Publisher Full Text

M Felsberg, G Sommer, The Monogenic Signal. IEEE Trans Signal Process 49(12), 3136–3144 (2001). Publisher Full Text

M Felsberg, G Sommer, The monogenic scalespace: A unifying approach to phasebased image processing in scalespace. J Math Imag Vision 21, 5–26 (2004)

M Felsberg, G Sommer, A new extension of linear signal processing for estimating local properties and detecting features. In 22. DAGM Symposium Mustererkennung, 195–202 (2000)

M Unser, D Sage, D Van De Ville, Multiresolution monogenic signal analysis using the rieszlaplace wavelet transform. IEEE Trans Image Process 18(11), 2402–2418 (2009). PubMed Abstract  Publisher Full Text

M Mellor, BW Hong, M Brady, Locally rotation, contrast, and scale invariant descriptors for texture analysis. IEEE Trans Pattern Anal Mach Intell 30(1), 52–61 (2008). PubMed Abstract  Publisher Full Text

V Grau, H Becher, J Noble, Registration of multiview realtime 3d echocardiographic sequences. IEEE Trans Med Imag 26(9), 1154–1165 (2007)

M Felsberg, U Köthe, GET: The connection between monogenic scalespace and gaussian derivatives. ScaleSpace 2005 (Springer LNCS, 2005) 3459, pp. 192–203

SL Hahn, The analytic, quaternionic and monogenic 2D complex delta distributions, report 3, Warsaw University Technology, Institute of Radioelectronics, Nowowiejska

T Bülow, G Sommer, Hypercomplex signals  a novel extension of the analytic signal to the multidimensional case. IEEE Trans Signal Process 49(11), 2844–2852 (2001). Publisher Full Text

U Kothe, M Felsberg, RieszTransforms vs. Derivatives: On the relationship between the boundary tensor and the energy tensor. in Scale Space and PDE Methods in Computer Vision, ed. by Kimmel, R, Sochen, N, Weickert, J (Springer, 2005), pp. 179–191 (vol, 2005), . 3459 of LNCS

H Knutsson, CF Westin, M Andersson, Representing local structure using tensors ii. in Image Analysis, ed. by Heyden, A, Kahl, F (Springer Berlin/Heidelberg, 2011), pp. 545–556 (vol, 2011), . 6688 of Lecture Notes in Computer Science

H Knutsson, AM Loglets, Generalized quadrature and phase for local spatiotemporal structure estimation. in Image Analysis, ed. by Bigun, J, Gustavsson, T (Springer Berlin/Heidelberg, 2003), pp. 107–108 (vol, 2003), . 2749 of Lecture Notes in Computer Science

H Knutsson, Representing local structure using tensors, Tech. rep (Computer Vision Laboratory, Linkoping University, 1989)

I Selesnick, R Baraniuk, N Kingsbury, The dualtree complex wavelet transform. Signal Process Mag IEEE 22(6), 123–151 (2005)

N Kingsbury, Complex wavelets for shift invariant analysis and filtering of signals. Appl Comput Harmonic Anal 10(3), 234–253 (2001). Publisher Full Text

N Kingsbury, The dualtree complex wavelet transform: A new technique for shift invariance and directional filters. 8th IEEE DSP Workshop (Bryce Canyon, 1998), pp. 319–322

N Kingsbury, Image processing with complex wavelets. Phil Trans Royal Soc Lond A 357, 2543–2560 (1997)

S Olhede, G Metikas, The monogenic wavelet transform. IEEE Trans Signal Process 57(9), 3426–3441 (2009)

V Murray, VP Rodriguez, MS Pattichis, Robust multiscale AMFM demodulation of digital images. IEEE International Conference on Image Processing 1, 465–468 (2007)

V Murray, MS Pattichis, H Davis, ES Barriga, P Soliz, Multiscale AMFM analysis of pneumoconiosis xray images. IEEE International Conference on Image Processing, 4201–4204 (2009)

V Murray, MS Pattichis, P Soliz, Retrieval of XRay images with different grades of opacities using multiscale amfm methods. Asilomar Conference on Signals, Systems and Computers, 12–16 (2009)

International Labour Office, Guidelines for the use of ILO International Classification of Radiographs of Pneumoconioses (Geneva, Switzerland, 1980)

H Hwang, R Haddad, Adaptive median filters: new algorithms and results. IEEE Trans Image Process 4(4), 499–502 (1995). PubMed Abstract  Publisher Full Text

S de Jong, SIMPLS: An alternative approach to partial least squares regression. Chemometr Intell Lab Syst 18(3), 251–263 (1993). Publisher Full Text

M Barker, W Rayens, Partial least squares for discrimination. J Chemometr 17(3), 166–173 (2003). Publisher Full Text

D Schonfeld, J Goutsias, Optimal morphological pattern restoration from noisy binary images. IEEE Trans Pattern Anal Mach Intell 13(1), 14–29 (1991). Publisher Full Text

NF Osman, ER McVeigh, JL Prince, Imaging heart motion using harmonic phase mri. IEEE Trans Med Imag 19(3), 186–202 (2000). Publisher Full Text

M Elshinawy, M Chouikha, Using amfm image modeling technique in mammograms. MicroNanoMechatronics and Human Science, 2003 IEEE International Symposium on 2, 660–663 (2003)

MY Elshinawy, J Zeng, SCB Lo, MF Chouikha, Breast cancer detection in mammogram with amfm modeling and gabor filtering. Proc 7th International Conference on Signal Processing 3, 2564–2567 (2004)

AO Boudraa, JC Cexus, H Zaidi, Functional segmentation of dynamic nuclear images by crossif_{b}energy operator. Comput Methods Progr Biomed 84(23), 146–152 medical Image Segmentation Special Issue Publisher Full Text

P Maragos, T Loupas, V Pitsikalis, Improving doppler ultrasound spectroscopy with multiband instantaneous energy separation. Proc 14th International Conference on Digital Signal Processing 2, 611–614 (2002)

E Alexandratou, A Sofou, H Papasaika, P Maragos, D Yova, N Kavantzas, Computer vision algorithms in DNA ploidy image analysis. in Proc of the SPIE Imaging, Manipulation, and Analysis ofBiomolecules, Celss, and Tissues, ed. by Farkas DL, Nicolau DV, Leif RC (SPIE, 2006), pp. 180–190 (vol, 2006), . 6088 of IV

C Christodoulou, P Kaplanis, V Murray, M Pattichis, C Pattichis, Classification of surface electromyographic signals using AMFM features. 19th International Conference on Artificial Neural Networks (Limassol, Cyprus, 2009)

CI Christodoulou, C Pattichis, V Murray, M Pattichis, A Nicolaides, AMFM representations for the characterization of carotid plaque ultrasound images. 4th European Conference of the International Federation for Medical and Biological Engineering 22, 546–549 (2009). Publisher Full Text

CP Loizou, V Murray, MS Pattichis, CS Christodoulou, M Pantziaris, A Nicolaides, CS Pattichis, AMFM texture image analysis of the intima and media layers of the carotid artery. 19th International Conference on Artificial Neural Networks (Limassol, Cyprus, 2009)

C Loizou, V Murray, M Pattichis, M Pantziaris, C Pattichis, Multiscale amplitude modulationfrequency modulation (amfm) texture analysis of ultrasound images of the intima and media layers of the carotid artery. IEEE Trans Inf Tech Biomed 15(2), 178–188 (2011)

S Acton, P Soliz, S Russell, M Pattichis, Content based image retrieval: The foundation for future casebased and evidencebased ophthalmology. Proc IEEE International Conference on Multimedia and Expo, 541–544 (2008)

C Agurto, S Murillo, V Murray, M Pattichis, S Russell, M Abramoff, P Soliz, Detection and phenotyping of retinal disease using AMFM processing for feature extraction. Asilomar Conference on Signals, Systems and Computers, 659–663 (2008)

E Barriga, V Murray, C Agurto, M Pattichis, S Russell, M Abramoff, H Davis, P Soliz, Multiscale AMFM for lesion phenotyping on agerelated macular degeneration. IEEE International Symposium on ComputerBasedMedical Systems (Albuquerque, New Mexico, 2009)

T Ezzat, J Bouvrie, T Poggio, AMFM demodulation of spectrograms using localized 2d maxgabor analysis. Proc IEEE International Conference on Acoustics, Speech and Signal Processing ICASSP 4, IV1061–IV1064 (2007)

F Gianfelici, C Turchetti, P Crippa, Multicomponent AMFM demodulation: The state of the art after the development of the iterated hilbert transform. IEEE International Conference on Signal Processing and Communications ICSPC, 1471–1474 (2007)

N Wang, E Ambikairajah, B Celler, N Lovell, Feature extraction using an AMFM model for gait pattern classification. Proc IEEE Biomedical Circuits and Systems Conference BioCAS, 25–28 (2008)

H Li, L Fu, Z Li, Fault detection and diagnosis of gear wear based on teagerhuang transform. Proc InternationalJoint Conference on Artificial Intelligence JCAI '09, 663–666 (2009)

DP Nguyen, R Barbieri, MA Wilson, EN Brown, Instantaneous frequency and amplitude modulation of EEG in the hippocampus reveals state dependent temporal structure. Proc 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society EMBS, 1711–1715 (2008)

R Prakash, R Aravind, Modulationdomain particle filter for template tracking. Proc 19th International Conference on Pattern Recognition ICPR 2008, 1–4 (2008)

PR Senthil, R Aravind, Invariance properties of AMFM image features with application to template tracking. Proc Sixth Indian Conference on Computer Vision, Graphics & Image Processing ICVGIP '08, 614–620 (2008)

N Mould, C Nguyen, J Havlicek, Infrared target tracking with AMFM consistency checks. Proc IEEE Southwest Symposium on Image Analysis and Interpretation SSIAI 2008, 5–8 (2008)

I Kokkinos, G Evangelopoulos, P Maragos, Texture analysis and segmentation using modulation features, generative models, and weighted curve evolution. IEEE Trans Pattern Anal Mach Intell 31(1), 142–157 (2009). PubMed Abstract  Publisher Full Text

P Tay, AMFM image analysis using the hilbert huang transform. Proc IEEE Southwest Symposium on Image Analysis and Interpretation SSIAI 2008, 13–16 (2008)

A Belaid, D Boukerroui, Y Maingourd, JF Lerallut, Phasebased level set segmentation of ultrasound images. IEEE Trans Inf Tech Biomed 15(1), 138–147 (2011)

DC Klonoff, DM Schwartz, An economic analysis of interventions for diabetes. Diabetes care 23, 390–404 (2000). PubMed Abstract  Publisher Full Text

C Agurto, M Pattichis, S Murillo, V Murray, M Abramoff, S Russell, E Barriga, H Davis, P Soliz, Detection of structures in the retina using amfm for diabetic retinopathy classification. 2009 Meeting of the Association for Research in Vision and Ophthalmology (2009)

TECHNOVISION Project:, MESSIDOR: methods to evaluate segmentation and indexing techniques in the field of retinal ophthalmology [http://messidor.crihan.fr/] webcite

AREDS, Areds database, [http://eyephoto.ophth.wisc.edu/ResearchAreas/AREDS/AREDSstdPhotoIndex.htm] webcite

P Soliz, MS Pattichis, J Ramachandran, DS James, Computerassisted diagnosis of chest radiographs for pneumoconioses. in SPIE, vol. 4322, ed. by Sonka, M, Hanson, KM (SPIE, 2001), pp. 667–675

M Pattichis, C Pattichis, C Christodoulou, D James, L Ketai, P Soliz, A screening system for the assessment of opacity profusion in chest radiographs of miners with pneumoconiosis. IEEE Southwest Symposium on Image Analysis and Interpretation, 130–133 (2002)

H Kondo, T Kouda, Detection of pneumoconiosis rounded opacities using neural network. Joint 9th IFSA World Congress and 20th NAFIPS International Conference 3, 1581–1585 (2001)

H Kondo, T Kouda, Computeraided diagnosis for pneumoconiosis using neural network. IEEE Symposium on ComputerBased Medical Systems, 467–472 (2001)

M Pattichis, T Cacoullos, P Soliz, New models for region of interest reader classification analysis in chest radiographs. Pattern Recogn 42(6), 1058–1066 digital Image Processing and Pattern Recognition Techniques for the Detection of Cancer Publisher Full Text

M Pattichis, H Muralldharan, C Pattichis, P Soliz, New image processing models for opacity image analysis in chest radiographs. IEEE Southwest Symposium on Image Analysis and Interpretation, 260–264 (2002)

B Van Ginneken, BM Ter Haar Romeny, M Viergever, Computeraided diagnosis in chest radiography: a survey. IEEE Trans Med Imag 20(12), 1228–1241 (2001). Publisher Full Text

M Pattichis, J Ramachandran, M Wilson, C Pattichis, P Soliz, Optimal scanning, display, and segmentation of the international labor organization (ILO) Xray images set for pneumoconiosis. IEEE Symposium on ComputerBased Medical Systems, 511–515 (2001)

CP Loizou, V Murray, MM Pattichis, I Seimenis, M Pantziaris, CS Pattichis, Multiscale amplitudemodulation frequencymodulation (amfm) texture analysis of multiple sclerosis in brain mri images. IEEE Trans Inf Technol Biomed 15(1), 119–129 (2011)

C Pitris, A Kartakoullis, E Bousi, AMFM techniques in the analysis of optical coherence tomography signals. J Biophotonics 2(67), 364–369 (2009). PubMed Abstract  Publisher Full Text

A Kartakoulis, E Bousi, C Pitris, AMFM analysis of optical coherence tomography signals. in Optical Coherence Tomography and Coherence Domain Optical Methods in Biomedicine XIII, vol. 7168, ed. by Fujimoto, JG, Izatt, JA, Tuchin, VV (SPIE, 2009), p. 71681M

A Kartakoullis, E Bousi, C Pitris, AMFM techniques in optical coherence tomography. in Optical Coherence Tomography and Coherence Techniques IV, vol. 7372, ed. by Andersen, PE, Bouma, BE (SPIE, 2009), p. 73720U

VP Rodriguez, MS Pattichis, Nested random phase sequence sets: a link between AMFM demodulation and increasing operators with application to cardiac image analysis. Proc 6th IEEE Southwest Symposium on Image Analysis and Interpretation, 196–200 (2004)

VP Rodriguez, MS Pattichis, Real time AMFM analysis of ultrasound video. 45th Midwest Symposium on Circuits and Systems 1, I216–19 (2002)

VP Rodriguez, M Pattichis, M Goens, R Abdallah, MHealth: Emerging Mobile Health Systems, Springer US, 491–507 (Ch, 2006), . ObjectBased Ultrasound Video Processing for Wireless Transmission in Cardiology

D Gill, N Gavrieli, N Intrator, Detection and identification of heart sounds using homomorphic envelogram and selforganizing probabilistic model. Proc Computers in Cardiology, 957–960 (2005)