search for




 

Measuring Bone Density Connectivity Using Dual Energy X-Ray Absorptiometry Images
International Journal of Fuzzy Logic and Intelligent Systems 2017;17(4):235-244
Published online December 25, 2017
© 2017 Korean Institute of Intelligent Systems.

Li Chen1, and V. B. Surya Prasath2,3,4

1Department of Computer Science and Information Technology, University of the District of Columbia, Washington, DC, USA, 2Computational Imaging and VisAnalysis (CIVA) Lab, Department of Computer Science, University of Missouri, Columbia, MO, USA, 3Division of Biomedical Informatics, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, USA, 4Department of Biomedical Informatics, College of Medicine, University of Cincinnati, OH, USA
Correspondence to: Li Chen (lchen@udc.edu)
Received August 25, 2017; Revised September 12, 2017; Accepted September 25, 2017.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract

Estimating bone density using mathematical models has direct applications in medical image processing. This paper presents a new measure for bone mineral density analysis based on the dual energy X-ray absorptiometry (DEXA) images. We are proposing an innovative procedure to calculate a scalar value that indicates the connectivity of bone mineral components from DEXA images. This method provides a totally novel measure for bone density study using DEXA scan images that can not only calculate T and Z values and can also determine the quality of the bone in terms of the average intensity of bone density. The new measure proposed is called the λ-measure and it provides new possibilities in terms of finding how well the bone components are connected. We can have λT and λZ corresponding to T and Z scores which augments the traditional values. Combining this method to a λ-connected maximum entropy concept, we obtain good segmentation results. The experiments results and data obtained shows that the new λ-measure should be in the range of [0.962, 0.977].

Keywords : Bone density, DEXA, λ-connectedness, Segmentation
1. Introduction

Bone density is an important feature to humans bone health and measuring it is an important task in medical image processing area. The mathematical models of bone density have direct consequences in possible predictions and making important diagnostics in terms of bone related treatments. Bone density is defined as the amount of bone tissue in a certain volume of bone [1]. The strength of a bone is believed to be related to and can be determined by the density of the bone structure. Moreover, bone density also measures the amount of bone mass present. It predicts fracture risks, observes the usefulness of treatments, and measures the amount of calcium in a specific region of the bone. The dual energy X-ray absorptiometry (DEXA) technique (see, Figure 1) is the most popular method used today to measure bone mineral density (BMD), which is the most important figure of bone density [24]. The bone density measurements allow for the diagnoses of osteoporosis and potential fracture risks at a specific site of a bone. During the DEXA test, a patient lies down on a table as shown in Figure 1(a). The DEXA scanner rapidly directs X-ray energy from two different sources towards the bone being examined in an alternating fashion at a set and fixed frequencies. Since two energy sources are used to measure bone density rather than the traditional X-ray or other radioisotope studies, the precision and accuracy of the measurements are greatly improved [1, 5, 79].

BMD measurements usually decrease with age, and people with a BMD lower than 1.0 g/cm2 need to take action in order to restore their bones. Note that, BMD does depend on the sex and age of the patient, and sometimes on the race of the patient as well. White and black females are at the greatest risk for osteoporosis or any other bone disease. Osteoporosis is known to affect 25 million people in the United States alone and half of all women over the age of 45 and past their menopause stage. It is the most common bone disease in USA today, and causes more than 1.3 million fractures including 250, 000 trendy fractures, 240, 000 wrist fractures, and 500, 000 spinal fractures [10]. In addition, the National Osteoporosis Foundation recommends that BMD testing should be done on all postmenopausal women under the age of 65 who have risk factors such as a low body weight, a history of fractures, smoking , and a family history of fractures. Furthermore, all women over the age of 65 should be tested as well, whether they have risk factors or not [1, 11]. However, normal DEXA measurements only produce two individual numbers: the T and Z scores [12, 13]. The T-score is the standard deviation of the patient’s BMD above or below the mean for the Young-Adult reference population. On the other hand, the Z-score is the standard deviation of the patient’s BMD above or below the mean of the Age-Matched population. The current treatment recommends that all women start the drug therapy if they have the BMD T score of less than −2 without other risk factors and for those with a T score of less than −1.5 if there are other risk factors [6].

In order to understand more about human bone density, researchers have developed several methods for bone modeling and the connectivity of bone mass for the detailed analysis of the bone components, especially for the trabecular bones [1417]. These methods include fractal analysis, finite elements, and the Monte-Carlo simulation. The finite elements method is for 3D bone modeling and is not currently directly used for DEXA images. Even though it is proven that fractal analysis is an effective measure for BMD, other methods shall be investigated intentionally to obtain a better or improved industry standard for the detailed bone density other than the classical T and Z scores.

In this paper, we present a new measure for BMD analysis based on DEXA images. We propose a determinative procedure to calculate a scalar value that indicates the connectivity of bone mineral components. This method provides a totally novel measure for bone density study using DEXA scan images that originally only calculates T and Z value that only provide the quality of the bone in terms of average intensity of bone density. The new measure proposed called λ-Measure that shows how well the bone components are connected. We could have λT and λZ corresponding to T and Z scores. Combining with this method to λ-connected maximum entropy, we can see the similar result is obtained. The data experiments showed that λ-measure should be in the range [0.962, 0.977]. This opens up new research directions in terms of bone mineral connectivity as well as density analysis from DEXA images. A preliminary research can be found in [18].

2. Background

Bone density is the amount of bone tissue in a certain volume of bone (g/cm3). (See Figure 2 [1]) It is often very difficult to establish a measure without damaging the bone; therefore, the BMD (g/cm2) is used instead [10, 13]. The strength of a bone can be determined based on BMD. In addition, BMD measures the amount of bone mass, predicts fracture risks, observes the usefulness of treatment, and measures the amount of calcium in a specific region of the bone. BMD can be measured using DEXA, the most accurate and less expensive technology available today. The normal average BMD is about 3.88 g/cm2 in males and 2.90 g/cm2 in females. The range of BMD in the forearm is from 700 to 800 g/cm2 and in the spine region is 1, 000 to 1, 200 g/cm2. The density of hard bone is a constant number of 1, 900 kg/m3 [10]. Considerable amount of research are conducted in this area [10, 15, 1924].

Three main companies that manufacture DEXA detector densitometers are: Hologic, Norlan, and Lunar. Usually, the detector scans two or more areas, the hip, spine and/or total body. Figure 2(a) shows the division of the four lumbar and Figure 2(b) shows a picture of a hipbone that was scanned by a Lunar-machine, and Figure 2(c) shows a report made by [13].

After the bone image has been scanned onto the computer screen, an algorithmic program in the QDR software automatically computes the results based on the bone scan using the following formula: The T-score is the standard deviation (SD) of the patient’s BMD above or below the mean for the Young-Adult normal reference population. The Z-score is the SD of the patient’s BMD above or below the mean of the Age-Matched population:

TScore=(BMD of patient-Average BMD ofYoung-Adult)SD of Young-Adult,ZScore=(BMD of patient-Average BMD ofAge-Matched)SD of Age-Matched.
3. Connected Component Segmentation

Image segmentation is a task which requires a scheme to partition an image into connected components with similar characteristics [25]. Each component is a targeted object. There are many kinds of segmentation, e.g. measurement of space clustering, region growing, split-and-merge segmentation, edge detection, active contours, etc. Different types of data need different segmentation methods. For gray-scale images such as velocity profiles, we can use region growing or edge detection. For texture images such as waveform profiles or seismic sections, a rule-based segmentation scheme may be needed.

A digital image F is a mapping from a grid-space ∑ to the real set ℝ (or ℝn in general). Two points p, q in ∑ are adjacent if the distance between these two points is equal to 1. Two points p, q in S ⊂ ∑ are connected if there is a path which consists of adjacent points that links p and q. The set S is connected if each pair of points in S is connected. We say S is uniform if F on S (a sub-image) has properties of uniformity. For instance, a popular uniformity measure is defined in terms of the maximum difference between any pixel value and the mean value of S. Segmentation can be defined as a partition of an image into connected subsets, each of which is uniform, but no union of adjacent subsets is uniform. The formal definition of segmentation is then: In a digital image F, if there exists a non-empty segmentation F1, F2, . . . , Fm , satisfying

  1. FiFj = ∅︀, if i = 1, . . . , m, and j = 1, ..,m, ij,

  2. i=1,...,mFi = F,

  3. Each Fi is connected,

  4. Each Fi is uniform,

  5. If Fi and Fj are adjacent, then FiFj is not uniform.

Then, {F1, F2, . . . , Fm} is called a segmentation of the given image F and each Fi represent a connected region.

We next briefly describe a discrete connected component based segmentation known as λ-connected segmentation. The λ-connected segmentation uses a graph-theoretic method to find a connected component in an image when only small changes appear in adjacent pixels [2633]. It is partially an inversion of gradually varied fitting [34]. The technique is defined on a graph as G = (V,E) with an associated potential function αε. A measure C(x, y) is defined in λ-connectedness on two vertices x, yG. For λ ∈ [0, 1], x and y are said to be λ-connected if C(x, y) ≥ λ. Figure 3(a) shows an application on seismic image segmentation with λ = 0.98. The left frame is the original image, and the right frame shows the five layers of segmented image. Each layer represents a connected component present in the image. This program can segment bone density images (See Figure 3(b), the original picture shown at top is from a DEXA scan image). Note the similarities in both the images, the objective is to segment distinct regions according to their grayscale value distribution. We will explain how the bone image analysis can be done using the generic λ-connectedness in the following section in more detail. We recall the definitions and results associated with λ-connectedness and refer to [35, Chapter 12] for more details.

Suppose that F is a digital image. Then, F is a mapping: F : ∑2 → ℝ, where ∑2 is a 2D array and ℝ is a set of real numbers. Thus, a pixel is a point in ∑2 × ℝ. Let d(x, y) be any standard distance function, in this work we use the Euclidean distance.

Definition 1

Let ε be a small real number. If p = (x, F(x)) and q = (y,F(y)) satisfy d(x, y) ≤ ε , then p and q are called ε-neighbors.

Definition 2

If a function αε : P × P → [0, 1], and p = (x, F(x)) and q = (y,F(y)) satisfy

  1. αε(p, p) = 1,

  2. αε(p, q) = αε(q, p),

  3. αε(p, q) = 0, if d(x, y) > ε,

  4. αε(p, q) ≥ αε(p, q′), if d(f(x), f(y)) ≤ d(f(x), f(y′)),

then αε is called the neighbor-connectivity with respect to ε.

A finite sequence Ω(p(x1), p(xn))={p(x1), p(x2), . . . , p(xn)} is called an ε-path, if d(xi, yj) ≤ ε, j = i + 1 and xixj , ij except that x1 = xn. The neighbor-connectivity of an ε-path Ω(p(x1), p(xn)) is defined as

αɛ(Ω(p(x),p(xn)))=min{αɛ(p(xi),p(xi+1)):i=1,,n-1},

or

αɛ(Ω(p(x),p(xn))=Παɛ(p(xi),p(xi+1)),i=1,,n-1.

Finally,

C(p(x),p(y))=max{αɛ(Ω(p(x),p(y)):Ω(p(x),p(y))is an ɛ-path}

is called the connectivity of p(x) and p(y) with respect to the ε-paths.

Definition 3

p(x) and p(y) are called λ-connected if

C(p(x),p(y))λ,         λ[0,1].

The following theorem shows the characteristics of the λ-connectedness according the neighborhood connectivity rule used.

Theorem 1

The λ-connectedness is an equivalence relation if (3) is used, i.e., it can be used for classification. The λ-connectedness is a similarity relation if (4) is used.

4. Maximum Entropy Criterion for λ Estimation

According to a study [36], we use the maximum entropy method to determine the λ value. It is λ-connected maximum entropy. The maximum entropy method was first proposed by Kapur et al. [37]. It is based on the maximization of inner entropy in both the foreground and background. The purpose of finding the best threshold is to make both objects in the foreground and background, respectively, as smooth as possible. Such an idea can be used for λ-connected segmentation. However, the total inner entropy for the image is to calculate the entropy for each segment (λ-connected component), not for the thresholding clipped foreground/background. This is because in λ-connected segmentation there is no specific background. Each λ-connected segment can be viewed as foreground, and the rest may be viewed as the background. It is different from the original maximum entropy where the range of pixel values determines the inclusion of pixels. Therefore, we need to summarize all inner entropies in all segments.

H(λ)=(inner entropy of each λ-connect.comp.).

We will select the λ such that H(λ) will be maximized. We call this λ-value the maximum entropy connectedness. This unique value is a new measure for images. Since the maximum entropy means the minimum amount of information or minimum variation, we want the minimum change inside each segment. This matches the philosophy of the original maximum entropy method. In other words, the λ-connected maximum entropy has a better meaning in some applications. We use the λe such that

H(λe)=max{H(x):x[0,1]}.

Assume there are m number of λ-components, define inner entropy of each λ-component Si:

Hi(λ)=k=0255(-Histogram[k]/n)(log(Histogram[k]/n)),

where n is the number of points in the component Si. Histogram [k] is the number of pixels whose values are k in the segment. Thus,

H(λ)=i=1mHi(λ).

Figure 4 shows the result for the DEXA image shown in Figure 3(b), using the maximum entropy connectedness presented in this section, we got λe = 0.95. Finally, Figure 5 shows an example of a four part lumbar data and its corresponding λ-connected segmentation results for various λ values. The original image shows a good strength bone. Figure 5(b), λ = 0.98, contains a some separated segmentation in each lumbar especially in the leftmost lumbar. It means λ = 0.98 is too large. Figure 5(c), λ = 0.97, and Figure 5(d), λ = 0.96, showed the good connectivity. Figure 5(e), λ = 0.95, and Figure 5(f), λ = 0.93, have connected more data points that may be noisy pixels. This example also proved that the previous result of [0.962, 0.977] range for the λ-measure is reasonable. A large data set for running the methods described above are necessary for obtaining a more accurate estimation and more research shall be done in the future.

5. A New Measure for Bone Connectivity

To determine the value of λ in segmentation is a critical issue. How to determine the value of λ that will yield a clear segmentation? The question is even though we can find such a λ by test a numbers of different value in segmentation. The value is not unique. We need to find a definitive procedure to determine the value of λ.

The selection of λ depends on how many segments we want to separate. In the following section, we will discuss the maximum-connectivity spanning tree method, which can be used to find the value of λ. Since the process of finding the maximum-connectivity spanning tree is time consuming, a histogram or an experimental test may be initiated to solve the problem. A fast algorithm to find the value of λ is the binary search. It is used to test the value of λ. First, test for λ= 0.5. If λ is too big then λ= 0.25 is tested. If λ is too small then λ= 0.75 is tested and so on. Similarly, for λ-connected fitting, the same method may also be used. The speed of the binary search method compared to the minimal spanning tree method for different sized bone images will also be studied. Finding the corrected λ value for different bones or for a different age group of patients is more important than finding a value for λ in terms of good segmentation. The value that can lead to segment a DEXA-image in a good understandable component may indicate the quality of the bone, just like the fractal dimensions. An appropriate λ value means to collect each bone component correctly. A large amount of bones from the young adults or a different age group may be needed to determine the best λ value. The best λ value means that it produces the satisfied bone components and the λ value is not too big or too small. The relationship between the λ value and the value of fractal dimensions is another important issue. If there is a clear relationship, the result will be very impressive at the science level. It could lead to new developments in human bone density analysis research.

When λ is small, the segmentation is rough. When λ is big, the segmentation is fine. We define

λmin={the largest value that segments the bone into only one component},λmax={the smallest value that segments the bone is as the same as when λ=1}.

The measure λm = 0.5(λmin + λmax). Or we can choose Golden-Cut-Ratio λ.618 = (.618λmin+.382λmax) and λ.618 = (.618λmax + .382λmin). With this value if it produces more segments that means the bone connectivity is not good. If it produces less segments, it means that the bone connectivity is good.

λTScore=(λmof patient-Average λmof Young-Adult)SD of Young-Adult,λZScore=(λmof patient-Average λmof Age-Matched)SD of Age-Matched.

Compared to the original measures given in Eqs. (1) and (2), we see that our new measure is based on the λ-connectedness thereby discerned directly from the DEXA images. The new measures are therefore robust as it depends on the data itself and can be computed automatically. Further, the the SD of Young-Adult and SD of Age-Matched are kept the same between the original scores and the new measure considered here. For example, in Figure 6 we check the λ-measure with the example showed in Figure 3(b). In this image, we do not display any segment that is smaller than 30 pixels. So, λmin = 0.940, λmax = 1, and λm = 0.970. λ.618 = (.618λmin +.382λmax) = 0.962 and λ.618 = (.618λmax+.382λmin) = 0.977. The measure of connectivity should be place in the range of [0.962,0.977]. The Table 1 lists the calculation results based on different λ-values and Figure 7 shows the corresponding graph. Based on this curve, we suggest a modeling function for the distribution.

P(λ)=p0d(λ-λmin)/(δλ)+c,

where λminλλmax and d is the dimension of the image. Again, λmin is the maximum λ-value making the image to be only one segment, and λmax is the minimum λ-value making the most of segments (It also can count the small percentage cut at the each end).

Remark 1

Note that the λ-connected segmentation is to segment an image if a value of λ is provided. λ-measure is to find a λ value that is the best to describe the connectivity or connectedness of an image. λ-connected segmentation is basically a dual segmentation technique comparing to the (multiple-) threshold segmentation method that translates a gray-scale image into a binary image based on a clip-level.

One can think about the thresholding as a vertical processing technique, but λ-connected search is a horizontal search. These two methods share the following properties: (1) they are simple to implement, and (2) The segmentation results generate equivalence relations. They can be cooperated by (a) using λ-connected segmentation first (to generate major segments) then multiple-thresholding (to generate detailed segments), or (b) using multiple-thresholding first (to generate major segments) then λ-connected segmentation (to generate detailed segments). If we use the same λ-value and thresholds, then we would guess that we will get the same detailed segments. This statement is not true since we have the following counter-example:

Mathematically, let Ta,b(P) be the threshold result for all pixel in image P that is greater than or equal to a and smaller than b, and let Sλ(P) be the results of λ-connected segmentation respect to value λ. Thus, apply λ-connected segmentation then use thresholding is denoted by Ta,b(Sλ(P)). Two pixels in Ta,b(P) is λ-connected then they must be λ-connected in P. Therefore, Sλ(Ta,b(P)) is the refinement of Ta,b(Sλ(P)), see Table 2 for illustration.

Theorem 2

Sλ(Ta,b(P)) ≤ref Ta,b(Sλ(P)).

Remark 2

Applying λ-connectedness based segmentation result in finding the appropriate thresholds quickly. Whereas to apply thresholding first will help in stopping the extreme region linking in the segmentation.

6. Conclusion

We considered a new measure for bone density based on λ-connectedness in DEXA images. Our proposed method can not only calculate BMD scores for the patients, but also determines the quality of the bon in terms of the average intensities. Our main aim here was to derive a new measure for BMD analysis based on the λ-connectedness analysis and comparisons with the original T- and Z-scores created from the DEXA exams is required on a larger dataset for understanding the upgraded guidance obtained by our automatic BMD calculations. For a general image, the best description of the image connectivity, may still be a λ-measure. However, the λ range might only be determined in a small range, how to find this range? Can we use information theory or statistical methods to this problem? Note that it is related to the research on the image pyramid and scale space based segmentation schemes [3840]. Some new investigation on bone density research can be found in [41]. These research directions are worthwhile to pursue in the future, along with the new developments on graph-based segmentation are still active, and it can be found in [42].

Acknowledgement

Many thanks to Drs. Richard Squier and Chin-Shoou Lin for providing helps in this research. The testing image in Figure 5(a) was from NIH.

Conflict of Interest

No potential conflict of interest relevant to this article was reported.


Figures
Fig. 1.

(a) Dual-energy X-ray absorptiometry (DEXA) procedure [5] (Dr. Caroline Lebreton, CHU Raymond Poincare, Garches, France). (b) Assessment of bone mineral density of the femoral neck (A) and the lumbar spine (B): T scores of −4.2 and −4.3 were found at the hip (A) and lumbar spine (B), respectively in a 53-year-old male patient affected with Fabry disease [6].


Fig. 2.

Typical DEXA scan images and report with annotations [13].


Fig. 3.

The λ-connectedness based image segmentations provides better region separation using a discrete connected measure as illustrated by (a) seismic image with three distinct regions separated by interfacial boundaries and (b) DEXA image which is also segmented accurately with different gray regions representing the density levels along with boundaries given in blue color. For color referred in the figure, the reader is referred to the PDF version of the article.


Fig. 4.

Maximum entropy criteria, λe = 0.95 based result of λ-connected segmentation.


Fig. 5.

Four part lumbar data and its λ-connected component segmentations show in different color codings. The segmentation results in (c), and (d) show good connectivity, and represent good results whereas the results in (b) is fragmented, and (e)–(f) include noisy pixels in the segments. For color referred in the figure, the reader is referred to the PDF version of the article.


Fig. 6.

Different λ values and their corresponding segmentation results. Note that the blue color indicates the border and background of the segments. For color referred in the figure, the reader is referred to the PDF version of the article.


Fig. 7.

Different λ ∈ [0, 1] values and the number of segments based on Table 1. The index of x and y axis represent the λ value and the corresponding number of segments obtained respectively. Note that we restricted our range of the λ between 0.9600 to 0.9950, and the number of segmentations above 22.


TABLES

Table 1

Partitioned segments according to the λ values

λPPlog P
0.9600224.6904161.342423
0.9650275.1961521.431364
0.97006481.80618
0.975021514.662882.332438
0.980052622.934692.720986
0.9850123735.171013.09237
0.9900256550.645833.409087
0.9915312355.883813.494572
0.9925354759.55673.549861
0.9935384762.024193.585122
0.9950402063.403473.604226

Table 2

Simple illustration of λ-connectedness and thresholding (example of the λ-connected component of P has two λ-connected components after thresholding)

Orignal PVariation 1T1,6
56556555
4014141
321321321

References
  1. Chawla, NV (2003). A preliminary paper on dual-energy X-ray absorptiometry (DEXA): Independent Research Tutorial, Georgetown University
  2. Lorente-Ramos, R, Azpeitia Arman, J, Munoz Hernandez, A, Garcia Gomez, JM, Diez-Martinez, P, and Grande-Barez, M (2011). Dual-energy X-ray absorptiometry in the diagnosis of osteoporosis: a practical guide. American Journal of Roentgenology. 196, 897904.
    CrossRef
  3. Lorente Ramos, RM, Azpeitia Arman, J, Arevalo Galeano, N, Munoz Hernandez, A, Garcia Gomez, JM, and Gredilla Molinero, J (2012). Absorciometria con rayos X de doble energia. Fundamentos, metodologia y aplicaciones clnicas. Radiologia. 54, 410-423.
    Pubmed CrossRef
  4. Gluer, CC (2017). 30 years of DXA technology innovations. Bone. 104, 7-12.
    CrossRef
  5. ,. What is DEXA Scanning. Available http://www.gorhams.dk/html/what_is_dexa_scanning.html
  6. Germain, DP (2010). Fabry disease. Orphanet Journal of Rare Diseases. 5.
    Pubmed KoreaMed CrossRef
  7. Lu, RS, Taylor, M, and Bottema, MJ (2016). Texture analysis improves the estimate of bone fracture risk from DXA images. Osteoarthritis and Cartilage. 24, S319-S320.
    CrossRef
  8. Staron, RB, Greenspan, R, Millar, TT, Bilezikian, JP, Shane, E, and Haramati, N (1999). Computerized bone densitometric analysis: operator-dependent errors. Radiology. 211, 467470.
    CrossRef
  9. Spence, AP (1990). Basic Human Anatomy. Redwood City, CA: Benjamin/Cummings
  10. Yarusskaya, A, and Donina, K. (2002) . Density of bone. Available https://hypertextbook.com/facts/2002/AnnaYarusskaya.shtml
  11. ,. Results for bone density testing beyond the basics. Available http://www.uptodate.com/contents/bone-density-testingbeyond-the-basics
  12. ,. (2013) . Bone density. Available http://courses.washington.edu/bonephys/opbmd.html
  13. ,. (2002) . DEXA reports. Available http://courses.washington.edu/bonephys/opDEXA1.html
  14. Richardson, ML, and Gillespy, T. Fractal analysis of trabecular bone. Available http://www.fractal.org/Life-Science-Technology/Publications/Fractal-Analysis-of-Trabecular-Bone.htm
  15. Ruimerman, R, van Rietbergen, B, Hilbers, PAJ, and Huiskes, R (2003). A 3-dimensional computer model to simulate trabecular bone metabolism. Biorheology. 40, 315320.
  16. Thomsen, JS, Barlach, J, and Mosekilde, L (1996). Determination of connectivity density in human iliac crest bone biopsies assessed by a computerized method. Bone. 18, 459465.
    CrossRef
  17. Michael, GJ, Sim, LH, and van Doorn, T (1997). A Monte Carlo model for bone mineral measurement using dual energy X-ray absorptiometry. Australasian Physical & Engineering Sciences in Medicine. 20, 8491.
  18. Chen, L 2006. -measure for bone density connectivity., Proceedings of IEEE International Symposium on Industrial Electronics, Montreal, Canada, Array, pp.489494.
  19. Yeh, OC, and Keaveny, TM (1999). Biomechanical effects three-dimensional finite element study. Bone. 25, 223228.
    CrossRef
  20. Verhulp, E, van Rietbergen, B, and Huiskes, R (2004). A three-dimensional digital image correlation technique for strain measurements in microstructures. Journal of Biomechanics. 37, 13131320.
    CrossRef
  21. Samarabandu, J, Acharya, R, Hausmann, E, and Allen, K (1993). Analysis of bone X-rays using morphological fractals. IEEE Transactions on Medical Imaging. 12, 466470.
    CrossRef
  22. Majumdar, S, Lin, J, Link, T, Millard, J, Augat, P, and Ouyang, X (1999). Fractal analysis of radiographs: assessment of trabecular bone structure and prediction of elastic modulus and strength. Medical Physics. 26, 13301340.
    CrossRef
  23. Michael, GJ, and Henderson, CJ (1998). Monte Carlo modeling of an extended DXA technique. Physics in Medicine and Biology. 43, 25832596.
    CrossRef
  24. Kinney, JH, and Ladd, AJ (1998). The relationship between three-dimensional connectivity and the elastic properties of trabecular bone. Journal of Bone and Mineral Research. 13, 839845.
    CrossRef
  25. Gonzalez, RC, and Wood, R (1993). Digital Image Processing. Reading, MA: Addison-Wesley
  26. Chen, L, Adjei, O, and Cooley, DH 2000. λ-Connectedness: method and application., Proceedings of IEEE Conference on System, Man, and Cybernetics, Nashville, TN, Array, pp.11571562.
  27. Chen, L (1994). Gradually varied surface and its optimal uniform approximation. Image and Video Processing II. 2182, 300307.
  28. Chen, L (1991). The -connected segmentation and the optimal algorithm for split-and-merge segmentation. Chinese Journal of Computers. 14, 321331.
  29. Chen, L (2002). λ-connected approximations for rough sets. Rough Sets and Current Trends in Computing. Heidelberg: Springer, pp. 572577
  30. Chen, L, Berkey, FT, and Johnson, SA (1994). Application of a fuzzy object search technique to geophysical data processing. Nonlinear Image Processing. 2180, 300309.
  31. Chen, L, Cheng, HD, and Zhang, J (1994). Fuzzy subfiber and its application to seismic lithology classification. Information Science: Applications. 1, 7795.
  32. Chen, L, Cooley, DH, and Zhang, L (1998). Intelligent data fitting technique for 3D velocity reconstruction. Application and Science of Computational Intelligence. 3390, 103112.
  33. Chen, L, and Adjei, O 2004. -connected segmentation and fitting., Proceedings of IEEE Conference on System, Man and Cybernetics, Hague, The Netherlands, Array, pp.35003506.
  34. Chen, L (1990). The necessary and sufficient condition and the efficient algorithms for gradually varied fill. Chinese Science Bulletin. 35, 870-873.
  35. Chen, L (2013). Digital Functions and Data Reconstruction: Digital-Discrete Methods. New York: Springer
    CrossRef
  36. Chen, L . Lambda-connectedness determination for image segmentation., Proceedings of 36th IEEE Applied Image Pattern Recognition Workshop, Washington, DC, Array, pp.7179.
  37. Kapur, JN, Sahoo, PK, and Wong, AKC (1985). A new method of gray level picture thresholding using the entropy of the histogram. Computer Vision, Graphics, and Image Processing. 29, 273285.
    CrossRef
  38. Prasath, VBS (2009). Color image segmentation based on vectorial multiscale diffusion with inter-scale linking. Pattern Recognition and Machine Intelligence. 5909, 339344.
  39. Meer, P, Sher, CA, and Rosenfeld, A (1990). The chain pyramid: Hierarchical contour processing. IEEE Transactions on Pattern Analysis and Machine Intelligence. 12, 363376.
    CrossRef
  40. Mathieu, C, Magnin, IE, and Baldy-Porcher, C (1995). Optimal stochastic pyramid: segmentation of MRI data. Medical Imaging VI: Image Processing. 1652, 1422.
  41. Bansal, SC, Khandelwal, N, Rai, DV, Sen, R, Bhadada, SK, Sharma, KA, and Goswami, N (2011). Comparison between the QCT and the DEXA scanners in the evaluation of BMD in the lumbar spine. Journal of Clinical and Diagnostic Research. 5, 694699.
  42. Mikesell, D, and Hicks, IV (2017). Image segmentation via weighted carving decompositions. Combinatorial Image Analysis. 10256, 268-279.
    CrossRef
Biographies

Li Chen is currently an Associate Professor of computer science at the University of the District of Columbia. He received his BS, MS, and PhD in CS from Wuhan University (1982), Utah State University (1995), and the University of Bedfordshire (2001), respectively. Chen’s research interests are broad in computer science and applied mathematics and include applied algorithm design, digital and discrete geometry, image processing, and applications to data science. He has made contributions to several research areas of computer science and its applications including : (1) lambda-connected image segmentation methods , (2) fuzzy relation equations , (3) the digital form of the Gauss-Bonnet theorem, (4) the polynomial time algorithm for finite Abelian group decomposition, (5) the definition of digital manifolds and classification of 3D digital surfaces, and (6) the optimum algorithm for the check matrix of the optimal SEC-DED code (optimal Hamming code). In 2014, Chen chaired the Satellite Conference on Data Science of International Congress of Mathematicians (ICM14). He is an ACM Distinguished Speaker since 2015.


Surya Prasath is an assistant professor (research) in the Computer Science Department at the University of Missouri, USA. He received his PhD in Mathematics from the Indian Institute of Technology Madras, India in 2009 (defended in March 2010). He has been a postdoctoral fellow at the Department of Mathematics, University of Coimbra, Portugal, for two years. Since 2012 he is with the Computational Imaging and VisAnalysis (CiVA) Lab at the University of Missouri, USA working on various mathematical image processing and computer vision problems. He had summer fellowships/visits at Kitware Inc. NY, USA, The Fields Institute, Canada, and IPAM, University of California Los Angeles (UCLA), USA. His main research interests include nonlinear PDEs, regularization methods, inverse & ill-posed problems, variational, PDE based image processing, and computer vision with applications in remote sensing, medical imaging domains.