^{1}

^{*}

^{2}

^{*}

Continuous wavelet transform is employed to detect singularities in 2 - D signals by tracking modulus maxima along maxima lines and particularly applied to microcalcification detection in mammograms. The microcalcifications are modeled as smoothed positive impulse functions. Other target property detection can be performed by adjusting its mathematical model. In this application, the general modulus maximum and its scale of each singular point are detected and statistically analyzed locally in its neighborhood. The diagnosed microcalcification cluster results are compared with health tissue results, showing that general modulus maxima can serve as a suspicious spot detection tool with the detection performance no significantly sensitive to the breast tissue background properties. Performed fractal analysis of selected singularities supports the statistical findings. It is important to select the suitable computation parameters- thresholds of magnitude, argument and frequency range - in accordance to mathematical description of the target prop erty as well as spatial and numerical resolution of the analyzed signal. The tests are performed on a set of images with empirically selected parameters for 200 μm/pixel spatial and 8 bits/pixel numerical resolution, appropriate for detection of the suspicious spots in a mammogram. The results show that the magnitude of a singularity general maximum can play a significant role in the detection of microcalcification, while zooming into a cluster in image finer spatial resolu tion both magnitude of general maximum and the spatial distribution of the selected set of singularities may lead to the breast abnormality characterization.

Wavelet transform modulus maxima method was developed for detection and characterization of signal singularities by Mallat and his collaborators [1-3]. Their method detects signal singularities by tracking the wavelet coefficients magnitude maximum over scale. They proved that, if a wavelet function is derivative of a Gaussian, wavelet transform modulus maxima must propagate towards finer scales. Although the representation by discrete wavelet maxima is not complete since several signals may exhibit the same wavelet maxima [

Arneodo and his team focused on how to recognize a sharp signal transition by tracking its behavior over scale. Specifically algorithms based on continuous wavelet transform modulus maxima method are able to detect singular points in a discrete 2-D and 3-D signals and supported by fractal analysis to give the metrics for the local signal regularity [5-9]. Arneodo’s team developed fractal based algorithm [

Microcalcifications in breast are residual calcium deposits that originate not only from completely normal processes, but also abnormal ones. Shape, morphology, and spatial distribution of individual microcalcifications are some of the features detectable in X-ray mammograms that suggest benign or malignant breast abnormality. Researchers have made few breast tissue classifications and have recognized more than twenty of those conventional features [10-12]. One acceptable simplification is to describe microcalcifications as ellipsoids of diameters between 0.05 mm and 1 mm. For early cancer detection, calcifications with spatial extent less than 0.5 mm are most important for clinical diagnosis. Particularly, this corresponds to calcifications roughly in the order of 0.1 mm in diameter. Also, microcalcifications appearing in clusters may suggest malignancy while individual occurrences are of low clinical significance [

Wavelets as mathematical microscope have been used by many researchers in mammogram analysis. Wang and Karayiannis [

Zhang and coworkers [_{z} = 0.92 under ROC curve, outperforming A_{z} = 0.86 for difference-image technique.

Microcalcifications in mammogram analysis may be modeled as impulse functions smoothed by Gaussian filters. Strickland and Hahn [

Since the Food and Drug Administration approved the technology in 1998, several Computer-Aided Detection (CAD) systems that support microcalcification detection have been deployed for clinical use. Literature reports show mixed results of current CAD systems in practice [17,18]. Results of some of the studies performed by using retrospective analyses or computer modeling suggested that CAD can achieve the main task increasing cancer detection.

Some studies show no practical improvements concluding that the CAD effect on the accuracy of interpretation is unclear [

This leads to a conclusion that different conditions of mammographic characteristics associated with X-ray exposure and breast tissue density should be studied across different institutions [

In this work, we consider microcalcifications as sharp signal transitions modeled as impulse functions smoothed by a Gaussian and characterized by negative Hölder exponent. We detect smoothed impulse functions and identify general modulus maxima and the scales in which they appear. Essentially detected general maxima show the scale in which the singularities are visible as the step functions with Hölder exponent approaching to zero. If in the higher scales the magnitude of the wavelet coefficients starts decreasing the singularity shows its impulse nature.

We detect the wavelet coefficient magnitude in the scale in which a singularity behaves as a step function and in higher scales becomes negative showing impulse nature of the singularity. Notice that normal breast tissue exhibits local signal fluctuation with many smoothed impulse functions in mammograms. This requires that in the next step detected smoothed impulse functions should be classified as significant for further analysis if they have locally dominant magnitudes of general maxima. In short, we take advantage of the capabilities of continuous wavelet transform modulus maxima method to detect signal singularities and differentiate between microcalcifications and normal signal fluctuation keeping track how Hölder exponent is developing over scale.

After detection in coarser spatial resolution, based on the fact that this method can detect step function, a cluster of singularities could be segmented in finer resolutions with appropriately adjusted settings. This method could be used for singularity border transition analysis i.e. its shape analysis for diagnostic purposes and therefore microcalcification estimation malignant vs benign [

In this work we demonstrate use of general maxima metric in mammograms with 200 μm/pixel spatial and 8 bits/pixel numerical resolution to speed up detection not of microcalcification border, but microcalcification itself. The scale of the general maxima detection depends on the size of the singularity and since singularity size in mammograms differs from one to another there is no typical scale but scale distribution that could be identified.

The magnitudes and their spatial distribution are analyzed in 25 regions containing microcalcification clusters from diagnosed 20 mammograms and the results compared to the results of 20 regions with healthy tissue from the same mammograms [

Background on the continuous wavelet transform, modulus maxima method, and fractal dimensions are presented in Section 2. In Section 3 we present the details of the framework and its implementation. The results from both statistical and fractal implementations, including the samples of tested mammograms, are presented in Section 4 while Conclusions are presented in Section 5.

The presence of strong local singularities is characterized by negative Hölder exponent while normal fluctuation is mostly characterized by positive Hölder exponent. This makes an isolated singularity detectable by analyzing the behavior of a modulus maximum along its maxima line.

Mallat and Hwang [

where is magnitude of the wavelet coefficient at point (x, y) in scale a, a is Hölder exponent at (x, y) and A > 0 is a constant.

When Hölder exponent at (x, y) is positive the magnitude of wavelet coefficient will increase with scale while when it is negative, as in case of impulse function, the magnitude of the wavelet coefficient will decrease. Mathematically local microcalcification of the mammogram signal may be modeled as an impulse function convolved by a Gaussian smoothing function.

For further analysis of behavior of magnitude of wavelet coefficients in scale, we use a wavelet coefficient magnitude approximation [

where a targeted microcalcification is modeled as an impulse function smoothed by Gaussian of variance σ^{2} and convolved by a Gaussian wavelet with variance β^{2}. _{}

If a tends to −1, then the magnitude of a wavelet coefficient increases up to the scale corresponding to the variance s^{2} reaching wavelet transform general modulus maximum followed by a decrease as in the case of microcalcification shown in

A Gaussian smoothing function at scale a

has been used as a kernel for wavelet transform discrete filters. The horizontal and vertical wavelet coefficients at scale a are defined as

and

respectively, where the horizontal and vertical wavelets and in scale a are derivatives of smoothing function f_{a} in the same scale a

and

respectively.

Wavelet coefficient at a scale a is represented by its magnitude and argument as follows,

where and

, with horizontal and vertical wavelet coefficients and at scale a, defined in (4) and (5) respectively.

A modulus maximum at (x_{0},y_{0}) at the scale a is represented by its magnitude if it is equal or greater than wavelet coefficient magnitude in its neighborhood. There is no detected modulus maximum at (x_{0},y_{0}) in scale a if any surrounding wavelet coefficient magnitude is grater than magnitude at (x_{0},y_{0}), i.e.

where is 8-neighborhood of (x_{0},y_{0}) in the scale a.

While a maxima line consists of the wavelet transform modulus maxima corresponding to the same local feature but in different scales, the general modulus maximum is a point on the maxima line characterized by its maximal magnitude and the scale in which it is detected. Wavelet transform general modulus maximum will appear with a magnitude and in a scale that depends on the initial variance of the smoothing Gaussian in the initial scale a_{0} = 1.

In the higher scales of corresponding maxima line will exist as long as the smoothed impulse function is the dominant singularity in its neighborhood. In that range the slope of the log-log diagram modulus maxima vs scale will tend to −1 if the maxima line is long enough, i.e. until it is overlapped by a modulus maximum of another dominant singularity. If in higher scales the singularity stops being dominant in its neighborhood, the maxima line tracking process should be interrupted to prevent misleading singularity detection. The magnitude and the scale in which the general modulus maximum is detected as well as the length of maxima line afterward are significant starting points for evaluating the isolated singularities.

A fractal is an object or quantity that displays self-similarity in all scales. A plot of the quantity on a log-log graph versus scale then gives a straight line, which slope is said to be the fractal dimension. The prototypical example for a fractal is the length of a coastline measured with different length ruler. The shorter the ruler, the larger the length measured, a paradox known as the coastline paradox or “the Richardson effect” (L. F. Richardson, 1881-1953) [

Fractal dimension or capacity dimension of a fractal is the exponent D in, where is the minimal number of balls of diameter e needed to cover the fractal set in each scale. Wavelets provide a natural generalization of the classical box-counting techniques to fractal signals, i.e. the wavelets are playing the role of generalized boxes. The wavelet transform as a mathematical microscope can be used to extract microscopic information about scaling properties of fractal objects [

The analogy related to multifractal formalism in thermodynamics allows the following definition of the exponent t(q) for the behavior of partition function along the scale:

where q and t(q) represent inverse temperature and free energy in statistical mechanics. For homogenous fractal functions-monofractals and q = 0, the value of −t(q) is associated to fractal dimension of the set of all selected singularities [

According to (11), fractal dimension of a selected set of singularities is determined as the negative value of the slope of log-log characteristics of partition function Z(q,a) vs scale a with q = 0 as shown in Section 4.

In our experiments the targets are the microcalcifications sized from 0.2 mm to 1 mm in diameters. Accordingly to the mammogram spatial resolution, the scale is changing in three octaves. Notice that algorithm parameters, such as initial smoothing function variance and length of maxima lines necessary for the target detection, depend on the mammogram spatial resolution and the size of expected targets. Also, the general modulus maxima analysis can be used in 1-D and 3-D signals in a similar fashion.

Changing the scale, Gaussian filters are used to compute wavelet transform coefficients of mammographic images in vertical and horizontal direction in three octaves. The framework is based on using wavelet transform modulus maxima to identify the maxima lines in a mammogram and then extract the general modulus maximum on each maxima line. Also in each maxima line a corresponding minimum is detected in higher scales than general maximum and before maxima line fades out. Singularities pointed by maxima lines with no identified corresponding minimums are considered as normal signal fluctuation and rejected from further analysis. Singularities with detected corresponding minimums have negative Hölder exponents. The fact that described general maximum and corresponding minimum are clearly detected means the log-log characteristic modulus maxima vs scale has negative slope and isolated smoothed impulse function is localized at the point where maxima line propagates to scale a = 1. A maxima line might not propagate to a single point at scale a = 1 because of fast oscillations in the cone of influence around terminal point of the maxima line in the lowest detected scale [

An illustration of applied framework is shown in

Can we fix that now? Should we reformat the entire text to avoid empty (blank) space or they would do it?

Variance of smoothed impulse function can be estimated by the scale in which general maximum is detected [

Gaussian kernels and their derivatives are used to enable image processing with no particular bias in the spacescale analysis. Significant fractal properties can change very fast in scale. They often become invisible within less than a 1/2 of an octave and with coarse scale increment the trend of modulus maxima behavior in scale could not be detected.

Each maxima line points to a singularity and each singularity is analyzed based on the modulus maxima behavior along its maxima line.

Detected singularities gradually change in scale both the magnitude M and the argument q of their modulus maxima. Two maxima points in two consecutive scale layers belong to same maxima line if their spatial locations and their arguments are similar. The level of similarity of spatial distance and argument difference are below empirically determined thresholds. The experiments show that the detection process converges if the scale increment between two consecutive scale layers is as low as inc = 2^{1/16}. Then the algorithm will properly form a maxima line with the argument threshold dA = p × log_{2}(inc)/4 radians or Dq < dA = p/64 radians in the simulations presented in Section 4.

In some specific situations when maxima line switches from one point to another and both are the parts of the same feature, the argument difference can be significantly larger than the argument threshold and then magnitude threshold will be checked as an alternative to the argument threshold. The magnitude threshold applied in the simulations presented in this work was adaptive in both space and frequency domain, i.e. the threshold was, where M_{a}_{1} and M_{a}_{2} are the modulus maxima detected in two higher consecutive scales. An example of described situation is presented in

With the parameter setup applied in this work we consider 1/4 of an octave a minimum scale distance between general modulus maximum and corresponding minimum to classify a location as an eligible singularity.

The mammograms from database [^{2} of breast tissue. Numerical resolution is 8 bits/pixel. In this work 12 malignant and 13 benign microcalcification clusters detected in 20 mammograms along with other 20 healthy regions selected from the same mammograms are used.

We find that normal and noise related detected eligible singularities are represented with significantly lower general modulus maxima comparing to those related to diagnosed microcalcification. We have introduced spatially local adaptive threshold in magnitude, to suppress false positive findings. The magnitudes above the local threshold we name (locally) dominant magnitude.

In this work we study magnitude and spatial distribution of general maxima of three different sets of detected singularities in the 64 × 64 pixel sub-images of the diagnosed 25 microcalcification clusters:

1) All detected singularities which show negative

Hölder exponent resembling smoothed impulse function2) The subsets of the set selected in 1, with dominant magnitudes of computed general maxima resembling microcalcifications, and 3) The sets of singularities that correspond to the diagnosed microcalcifications in 25 clusters.

Assuming that the nature of malignant and benign microcalcification clusters is different we additionally study each category described above separately as malignant and benign abnormality.

In addition we study magnitude and spatial distribution of general maxima of three different sets of detected singularities in 20 healthy areas of same mammograms:

1) All detected singularities which show negative Hölder exponent and resembling smoothed impulse function in healthy areas, and 2) The subsets of 1 that have dominant magnitudes of computed general maxima in health areas.

Finally, we compare our results with respect to the background breast tissue type: dense (D), glandular (G), and fatty (F) tissue.

In 25 microcalcification clusters our method detected total of 3115 eligible maxima lines. The histogram of eligible general maxima magnitudes is shown in

For the comparison, the histograms of the detection of eligible general maxima magnitudes in 20 healthy areas

is determined and presented as a dashed black line in

For the purpose of further analysis of the diagrams in

where x is general maxima magnitude,

is cumulative distribution function of general maxima for diagnosed microcalcifications (red line) and is cumulative distribution function of general maxima for total eligible general modulus maxima (blue line) both in [0, x].

While an impulse may be negative, a microcalcification is, in general, a spot brighter than its surroundings. Similarly to the function, here we institute the function, the ratio of detected impulses with the local intensity maximums in their neighborhood and total detected smoothed impulse functions. i.e.

where x is general maxima, magnitude,

is cumulative distribution function of general maxima for smoothed impulses with detected local maximum in each of their neighborhood (red line) and

is cumulative distribution function of general maxima for total eligible general modulus maxima (blue line).

In the same set of 20 mammograms we also choose 20 healthy areas and select all singularities that could be modeled as smoothed impulse function. The histogram of the total number of selected singularities vs gM magnitude is shown as dashed black line in

The diagram of the function in the set of analyzed sub-images with diagnosed microcalcification clusters is the black line in

For each line in ^{−5} and 4 × 10^{−5} in the case of healthy tissue samples, while in the case of diagnosed abnormalities for both lines the linearity (red and full black lines (having positive slope) starts above 4 × 10^{−5}.

To confirm this observation about positive and negative slope we performed same experiment separating the results for dense (D), glandular (G), and fatty (F) tissue. In

In same experiment we notice the difference in behavior of the function in the normal subimages (no microcalcifications-

We conclude that positive vs. negative slope of the linear parts on the right hand side leads to a good differentiation between abnormal and normal tissue sample. We notice that the average magnitude of the general modulus maxima in 25 clusters is 4.18 × 10^{−5}, while in 20 healthy background sample is 3.59 × 10^{−5}. Both averages are in the range of gM magnitude that could be considered in either non-linear or linear part of corresponding lines. This leads to another conclusion that a successful differentiation normal vs. abnormal findings could be expected when an adaptive threshold is set to the local average of gM magnitudes.

For the abnormalities below the local averages the additional microcalcification features should be considered.

In tested 25 clusters of microcalcifications we found 3115 eligible maxima lines that point to smoothed impulses. For each cluster centered in a 64 × 64 pixel subimage we calculated average general modulus maxima magnitude. The results of the calculation of general modulus maxima averages are shown in

The local and total averages, presented in ^{−5} (^{−5} (

(

The minimal computed average magnitude in an abnormal dense area was 3.49 × 10^{−5} (^{−5} (^{−5}) would not result in proper classification. Although the microcalcification related additional criterions would help the proper classification, in this analysis we evaluate the magnitude of general maxima in microcalcification detection leaving the combination of criterions for further research.

As we said, the averages in Column 5 and Column 6 in

In following analysis, a magnitude threshold of the general modulus maxima was determined as the local average of all gM magnitudes of all eligible maxima lines with detected negative Hölder exponent in the area 64 × 64 pixels centered in a diagnosed cluster. This choice of the adaptive threshold has reduced the number of dominant singularities down to 970 (31.1%) of all 3115 eligible singularities (_{z} = 0.77.

In the case of the malignant microcalcifications we find 33.3% as suspicious, while the malignant microcalcifications were 48.7%, meaning that the average threshold is too high. In _{z} = 0.80.

We find a different tendency in the case of benign microcalcification. The average threshold classified 29.1%, while benign microcalcifications were 27.3%, meaning that the average threshold was set a little lower than optimal. False positive and false negative findings were similar (10.7% and 9.0%, respectively). Surprisingly, the other statistical measures were not significantly different, i.e. true positive fraction was 67.1%, while ROC curve area index was A_{z} = 0.76. The details are shown in Tables 2 and 3.

Average gM magnitude of microcalcification was 4.58 × 10^{‑5}, while the average of all eligible maxima lines was 4.18 × 10^{‑5}, leaving additional room for selection improvement based on the magnitude of general maxima. Similar conclusion holds for separately analyzed both malignant and benign clusters where average magnitudes corresponding to microcalcifications were 4.34 × 10^{−5} and 5.02 × 10^{−5} respectively while threshold averages were 4.12 × 10^{−5} and 4.23 × 10^{−5} respectively. The difference of average magnitudes 4.34 × 10^{−5} and 5.02 × 10^{−5} for malignant and benign microcalcifications may lead to microcalcification malignancy probability estimation in further research.

In the following sub-sections, we determine fractal dimension of selected singularities for 3 types of data-all detected singularities that resemble smoothed impulse function, detected smoothed impulse functions with dominant general modulus maxima, and the detected singularities that correspond to singularities in diagnosed microcalcification clusters.

Partition functions Z(q, a), q = 0 and 0.5 £ log_{2}(a) £ 1.5, of all detected singularities in 25 microcalcification clusters are presented in

Specifically the slopes of log_{2}(Z(0,a)) vs. log_{2}(a), 0.5 £ log_{2}(a) £ 1.5, correspond to the fractal dimension of the sets of singularities for all detected singularities (rhombi), as well as for 12 malignant (squares) and 13 benign (triangles) clusters. The trend-lines show the slope is close to 2 in all three analyzed datasets. Notice that in 2-D space, fractal dimension can not be considered higher than 2. The results show that the eligible singularities are uniformly distributed everywhere in the analyzed mammogram sub-images.

Similarly, the slopes of log_{2}(Z(0,a)) vs log_{2}(a), 0.5 £ log_{2}(a) £ 1.5, corresponding to selected dominant singularities, as well as diagnosed microcalcification clusters are determined and shown in

In the case of microcalcification clusters, dominant and diagnosed singularities have similar fractal dimension values (between 1.4 and 1.5), which indirectly supports the assumption that dominant eligible singularities point to microcalcifications. There is no difference between malignant and benign cluster fractal dimensions in 200 mm/pixel spatial resolution. Notice that for microcalcification cluster classification purpose spatial resolution 200 μm/pixel doesn’t allow zooming at 70 μm/pixel where microcalcification edge becomes visible and malignancy analysis by fractal dimension possible.

Dominant singularities in healthy areas have significantly lower fractal dimension, i.e. tending to 1 for the average of all 20 mammograms. Fractal dimension approaching 1 is the sign of linear structures typical for detected edges. Fractal analysis shows that smoother signal transition is uniformly distributed everywhere producing fractal dimension tending to 2 and 1.73 in normal and

abnormal breast tissue respectively. On another hand fractal dimension of detected dominant signal transition in normal tissue is tending to 1 implying that the dominant singularities belong to the curve structures.

Healthy tissue fractal analysis shows that the set of all detected singularities has fractal dimension of 1.73 while the set of locally dominant singularities has fractal dimension tending to 1. This difference approves the significance of general modulus maxima in the detection of microcalcification. Fractal dimension of 1 is typical for the linear structures in the tissue and can be used in mammogram segmentation which we are proposing to be a future extension of this work.

Healthy tissue statistical analysis confirms our expectation that general maxima of normal fluctuation with negative Hölder exponent are lower than the general maxima belonging to the maxima lines representing the abnormalities suggesting that a combination of global and local thresholding need to be employed.

Our proposed algorithm successfully localizes both malignant and benign microcalcification clusters and it is worth to notice that no other features but dominant general modulus maximum with local average threshold is employed in the detection of microcalcification clusters such as those shown in

Wang and Karayiannis [

mammograms. Based on that idea we performed detection of energy blobs of detected dominant general modulus maxima. We consider that a microcalcification is a bright spot with a local intensity maximum in the neighborhood of a singular point related to the microcalcification and eliminate the singular point without a local intensity maximum in its neighborhood. The mammogram mdb223 [

Continuous wavelet transform was employed to detect singularities in mammograms by tracking modulus maxima along maxima lines across the scale and integrating it with fractal properties to differentiate between normal signal fluctuation and irregular ones in breast tissue.

Microcalcifications in mammograms were modeled by smoothed impulse function and detected by their negative Hölder exponent. Proposed procedure for maxima line classification enables selected general maxima to be analyzed by their magnitudes.

This method is capable to make difference between smooth signal transitions inevitable present in mammograms and microcalcification as sharper signal transition. Presented numerical results show that selecting maxima lines with negative Hölder exponent and local thresholding of general modulus maxima make possible the detection of the suspicious spots as local energy peaks of general modulus maxima.

Dominant general modulus maxima highly map into diagnosed microcalcifications in terms of local average magnitude and spatial distribution. Simulation results demonstrated that smoothed impulse function, recognized by their fractal properties of localized wavelet transform modulus maxima, can be considered a reliable feature for microcalcification detection in mammograms. Moreover, comparison of healthy and calcified breast tissue based on both statistical and fractal results show that selected general modulus maximum is a significant feature in microcalcification detection and classification.

Our results show that proposed metric for microcalcification detection is not sensitive on background tissue and therefore particularly convenient for dense parenchyma background in which microcalcification detection can be often a very challenging task.

Although a part of information has been lost by spatial resolution resampling from 50 mm to 200 mm microcalcifications preserve their smoothed impulse function nature, which has been detectable by proposed method. We showed that energy peaks in modulus maxima field correspond to microcalcification suspicious spots.

Coarse resolution makes faster suspicious spots searching procedure in real time during the screening procedure, which can be beneficiary by preventing an additional screening and saving the patients from having to endure the uncomfortable feelings during the procedure and more importantly, emotionally surviving cancer while waiting for additional test results.

Involving local intensity maximum of underlying mammogram in the neighborhood of detected smoothed impulse function significantly reduces false positive findings. Then the algorithm allows smaller threshold value of the general modulus maxima and therefore becomes more sensitive by detecting eligible smoothed impulse function with lower magnitude necessary for early detection of breast cancer.

The number of detected microcalcifications per area unit could be another feature that could reduce false positive findings.

Our findings in the difference of average general maxima magnitudes for malignant and benign microcalcifications leave a space for microcalcification malignancy probability estimation in further research. Namely, the numerical results show that general modulus maxima average of cancerous microcalcifications is lower than the average of benign ones. Lower modulus maximum magnitudes of cancerous findings make the detection more difficult, but open a space for improved classification in terms malignant vs. benign microcalcifications. Additional analysis might show that general maximum magnitude could serve as a parameter for malignancy probability estimation in real time.

The areas with microcalcification clusters tend to have higher fractal dimension than the healthy areas. Practical implementation of how to use fractal dimension in microcalcification detection and/or classification can be another interesting future work.

We showed that initial microcalcification segmentation could be performed based on wavelet transform modulus magnitude and phase and underlying mammogram in coarse resolution not significantly disturbed by noise. If necessary the microcalcification segmentation could be then refined in the finer resolution in which image acquisition was performed.

We found that our implementation of this proposed work including all parameter values is appropriate for the 200 μm/pixel spatial resolution and 8 bits/pixel numerical resolution. This approach could be easily applicable to other mammogram databases with different spatial and numerical resolution if proper adjustment to the estimates of the parameters values were taken into consideration.

Finally, involving additional microcalcification features could decrease false positive findings by enabling lowering modulus maximum threshold, which improves detection sensitivity of applied algorithm and therefore the general modulus maximum of smoothed impulse function can be used not only independently but also integrated with other microcalcification detection algorithms.