US20030065264A1 - Bone age assessment using ultrasound - Google Patents
Bone age assessment using ultrasound Download PDFInfo
- Publication number
- US20030065264A1 US20030065264A1 US10/042,735 US4273501A US2003065264A1 US 20030065264 A1 US20030065264 A1 US 20030065264A1 US 4273501 A US4273501 A US 4273501A US 2003065264 A1 US2003065264 A1 US 2003065264A1
- Authority
- US
- United States
- Prior art keywords
- ossification
- acoustic
- actuated
- skeletal structure
- bone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4504—Bones
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
- A61B8/0858—Clinical applications involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
- A61B8/0875—Clinical applications for diagnosis of bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4514—Cartilage
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4528—Joints
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/58—Testing, adjusting or calibrating the diagnostic device
Definitions
- the present invention relates to assessing bone age using ultrasound.
- Bone age assessment in growing subjects is a well-known diagnostic tool that is especially useful in predicting stature and/or growth problems in children, teenagers and adults.
- Ultrasonic methods of bone density measurement are known. For instance, according to prior art ultrasound measurement systems, (for example, WO 00/28316 and U.S. Pat. No. 5,564,423, the disclosure of which are incorporated herein by reference), ultrasound is used to determine the density of a non-cartilage osseous structure.
- An aspect of some embodiments of the invention relates to estimating bone age from acoustic signals of ossifying structures.
- bone age is estimated by measuring an acoustic velocity in cartilage structures that are in the process of ossifying wherein the velocity is expected to increase as a function of the ossification during the human maturation process.
- acoustic signal velocity for example, is measured in primary ossification centers such as the bones of the wrist or secondary ossification centers such as the distal regions of the ulna and radius.
- two or more acoustic signals of one or more ossification centers are measured to determine bone age. Additionally or alternatively, the ratio of acoustic signals between two or more ossification centers is used to measure bone age. Additionally or alternatively, acoustic signals from structures associated with ossifying structures are used to measure bone age, for example, the fibrocartilage of the pubic symphysis, skull suture ligaments and tooth and mandibular changes.
- Bone age may be used to predict adult stature or other aspects of the maturation process. Such predictions are based on bone age derived by other methods known in the art. Alternatively or additionally, tracking of ossification in bone is used to detect and/or track the progress of various disease states and/or disorders, with, for instance, a more accurate profile than X-ray evaluation due to its non-ionizing nature, allowing frequent monitoring without harm.
- parameters other than acoustic velocity such as broadband ultrasound attenuation (BUA) and dispersion of ultrasound signal
- BOA broadband ultrasound attenuation
- dispersion of ultrasound signal are used to estimate bone age, for example by correlating these parameters with the known BA assessment of a group of children.
- signals reflected from bone are used to measure bone age, for example, by measuring attenuation of backscatter intensity of the ultrasound signal.
- a scanning or multi-beam measurement system may be used.
- the acoustic signal provides a spatial measure, for example, indicating a profile of velocity along a bone axis or a radial profile of an ossification center.
- Different bones, different measures and/or different measurement systems may be used for different situations and/or for analyzing different bone ages or disease states.
- An aspect of some embodiments of the invention relates to using an existing osteoporosis measurement device, possibly with minimal changes, to assess bone age.
- a device designed for measuring osteoporosis in a finger is reprogrammed with a table that associates acoustic velocities with bone ages, rather than with states of osteoporosis.
- the transducers for measuring acoustic signal are modified specifically for application to growth centers. It is noted that velocity limits used in osteoporosis measurement-devices are designed to obtain measurements from non-growth center areas.
- a method for measuring bone age comprising:
- said ossification-actuated skeletal structure comprises one or more areas undergoing ossification.
- said ossification-actuated skeletal structure comprises one or more bones.
- said ossification-actuated skeletal structure comprises one or more regions of cartilage.
- said ossification-actuated skeletal structure comprises one or more regions of non-cartilage soft tissue.
- said ossification-actuated skeletal structure comprises one or more regions of fibrocartilage.
- said ossification-actuated skeletal structure comprises a region with one or more primary ossification centers.
- said ossification-actuated skeletal structure comprises one or more of: the bones of the wrist, the bones of the palm, the bones of the tarsus, the mandible.
- said ossification-actuated skeletal structure comprises a region with one or more secondary ossification centers.
- said ossification-actuated skeletal structure contains an epiphysis.
- said ossification-actuated skeletal structure comprises a region of one or more of: an ulna, a radius, a femur, a bone of a ray of an extremity.
- said receiving comprises utilizing two or more different acoustic signals to provide a measure of bone age.
- said two or more acoustic signals are associated with the same bone.
- said two or more acoustic signals are associated with paths in different bones.
- said two or more acoustic signals are received from the same direction.
- said two or more acoustic signals are received from the different directions.
- said signal passes through said one or more structures including an ossification-actuated skeletal structure.
- said signal echoes from said one or more structures including an ossification-actuated skeletal structure.
- said analysis of said signal is responsive to speed of sound from said one or more structures including an ossification-actuated skeletal structure.
- said analysis of said signal is responsive to broadband ultrasound attenuation from said one or more structures including an ossification-actuated skeletal structure.
- said analysis of said signal is responsive to dispersion of ultrasound from said one or more structures including an ossification-actuated skeletal structure.
- said analysis of said signal is performed, at least in part, in the frequency domain.
- said analysis of said signal is performed, at least in part, in the time domain.
- said analysis of said signal is responsive to attenuation of an ultrasound signal in said one or more structures including an ossification-actuated skeletal structure.
- said analysis is used to predict adult stature.
- said analysis is compared to a database having correlation with one or more of: conventional radiographs, CT images, MRI images and Nuclear Medicine scans.
- said transmitting is by a scanning acoustic signal transmitter.
- said transmitting is by a multi-beam acoustic signal transmitter.
- said receiving provides two or more acoustic signal measures along an axis of said one or more structures including an ossification-actuated skeletal structure.
- said receiving provides two or more acoustic signal measures radially around said one or more structures including an ossification-actuated skeletal structure.
- said analysis is correlated with a known bone age measurement system.
- said analysis is responsive to a formula providing a correlation with a known bone age measurement system.
- an estimate of bone age is responsive to time of flight of an acoustic signal between two transducers, with said ossification activated skeletal structure being situated intermediate to said transducers.
- said acoustic information is constructed into a database of bone age measurements.
- said database is arranged according to one or more of: sex, ethnic group, geographic location, nutrition and general inheritance.
- said database includes two or more measurements of one or more of said one or more structures including an ossification-actuated skeletal structure.
- said database includes one or more measurements of two or more growth stages from said one or more structures including an ossification-actuated skeletal structure
- said database includes one or more measurements of said one or more structures including an ossification-actuated skeletal structure in two or more populations.
- said received signals are compared to similar signals in a database to predict one or more of predict one or more of adult bone length, density, thickness and resilience and adult stature.
- said received signals are compared to similar signals in a database to indicate a bone-growth related disorder.
- said received signals are compared to similar signals in a database to track the progress of a bonegrowth related disorder.
- said received signals are compared to similar signals in a database to track hormone therapy in a growth stature disorder.
- said received signals are compared to similar signals in a database to indicate one or more growth-plate related disease states, including osteogenic sarcoma, slipped growth plate, premature arrest of growth plate growth and inflammation of growth plate.
- two or more acoustic measurements are made on a single subject and entered into said database.
- said two or more acoustic measurements are compared to track one or more growth-related disorders, including precocious puberty, delayed puberty, rickets, kwashiorkor, hypoparathyroidism, pituitary dwarfism and diabetes.
- said two or more acoustic measurements are compared to track treatment of one or more growth-related disorders, including precocious puberty, delayed puberty, rickets, kwashiorkor, hypoparathyroidism, pituitary dwarfism and diabetes.
- an apparatus for estimating bone age comprising:
- an acoustic transmitter and an acoustic receiver positioned on either side of one or more structures including an ossification-actuated skeletal structure;
- an electronic moveable gantry that adjusts the position of said acoustic transmitter and said acoustic receiver in relation to said ossification-actuated structure
- said apparatus transmits and receives one or more acoustic signals linearly along an axis through said ossification-actuated structure.
- said apparatus transmits and receives one or more acoustic signals radially around an axis through said ossification-actuated structure.
- said computer system controls said acoustic signal transmitter to provide an acoustic signal appropriate for said ossification-actuated structure.
- said computer system estimates said bone age responsive to one of more of broadband ultrasound attenuation, acoustic backscatter, dispersion of acoustic signal and speed of sound in said ossification-actuated structure.
- said computer system uses an imager to control the position of said acoustic signal receiver and said acoustic signal transmitter.
- said computer system contains a visual display to provide information on said bone age.
- said visual display comprises a graph.
- said computer system is comprised in a computer network.
- said computer system comprises a neural network.
- said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures including an ossification-actuated skeletal structure to provide an estimate of bone age.
- said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures, including an ossification-actuated skeletal structure to predict stature.
- said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures including an ossification-actuated skeletal structure to indicate, track or follow treatment of one or more of: a bone-growth related disorder, a growth plate disorder, and a growth related disorder.
- FIG. 1 schematically illustrates ultrasound measurements of ossifying wrist and arm bones of a 6.5-year-old male, in accordance with an embodiment of the present invention
- FIG. 2A schematically illustrates ultrasound measurements of a portion of ossifying arm, wrist and hand bones of an 11 year old female in accordance with an embodiment of the present invention
- FIG. 2B schematically illustrates a typical graph of ultrasound measurements associated with a path along the bones of FIG. 2A, in accordance with an embodiment of the present invention
- FIG. 3A schematically illustrates ultrasound measurements of ossifying femora and pelvis of a male, age 14, shown in a radiographic representation in accordance with an embodiment of the present invention.
- FIG. 3B schematically illustrates a pubic symphysis, shown in cross-section and a typical graph of its ultrasound measurement in accordance with an embodiment of the present invention.
- FIG. 4 is a graph of Bone Age of the hand assessed by ultrasonic multi parameter model as compared to standard X-ray.
- FIG. 1 illustrates ultrasound measurements of a portion of an arm 150 and a wrist 100 of a 6.5-year-old male in accordance with an exemplary embodiment of the present invention.
- Each bone of wrist 100 exhibits an ossification center, represented by a crosshatched area, that begins from a central portion and progresses toward the periphery and is referred to as a “primary ossification center”.
- Bones with only primary ossification centers continue to show changes in acoustic velocity due to increased ossification, increased bone size and/or changes in bone shape through the 15th-19th years.
- Some bones have additional ossification centers that appear as a strip of cartilage between ossified bone sections, referred to as a “secondary ossification center.”
- a radius 232 contains a secondary ossification center consisting of a cartilaginous growth plate 242 that is situated between an epiphysis 232 b and a metaphysis 232 a , both of which are ossified.
- Growth plate 242 first appears in the first postnatal year and ossifies in the seventeenth year in females and in the nineteenth year in males.
- a secondary ossification center ossifies, at a later age, epiphysis 232 b and metaphysis 232 a fuse together with growth plate 242 that ossifies as well. After the stage when fusion occurs, bone 232 is fully ossified lacking the bands of cartilage and hence is radio-opaque throughout.
- Ulna 244 similarly contains a growth plate 244 a that appears in the fifth year in females and in the sixth year in males. Ulnar growth plate 244 a ossifies in the seventeenth year in females and eighteenth year in males.
- two transducers such as a receiver 122 and a transmitter 120
- a receiver 122 and a transmitter 120 are held in a gantry 174 and placed so that they are in contact with the skin over ulna 244 and radius 232 .
- the thickness of soft tissue covering the bone is measured or the effect of soft tissue is optionally ignored.
- Transducers 122 and 120 are positioned so that they transmit along a region 276 , to obtain acoustic signals from ulna 244 and radius 232 and their respective ossification centers. This method of ultrasound transmission is referred to herein as a “through” method.
- the frequency of transducers 122 and 120 is between 300 kHz to 2 MHz.
- the distance between transducers 122 and 120 is measured using a ruler in which the distance is indicated by a digital readout 182 .
- the transducers are controlled and the received signals processed by a controller 142 (e.g., a computer), which optionally includes a display 140 and/or a printer 144 .
- a controller 142 e.g., a computer
- the transducers are controlled and the received signals processed by a controller 142 (e.g., a computer), which optionally includes a display 140 and/or a printer 144 .
- the acoustic signal from a growth plate, such as ulna growth plate 244 a is markedly different than that of bone, such as bone 244 .
- One parameter that can be measured from the obtained acoustic signal is the speed of sound.
- the distance between transducers is evaluated by digital ruler 182 that is connected to transducers 120 and 122 and the average speed of sound is calculated by dividing the distance between the transducers by the transmission time. Alternatively or additionally the distance between the transducers can be measured using optic or acoustic distance measurement devices.
- the speed of sound in ulnar cartilage 244 a is about 1700 meters per second (m/s) while speed of sound in bone 244 is about 2000-4500 m/s, depending on probe localization and age.
- the transmission time of a signal varies based upon the composition of the structures in the path, which changes with age.
- An ultrasound signal passing along region 276 will pass through ossified ulna 244 , ulnar growth plate cartilage 244 a , and ulnar epiphysis 248 that is also ossified. It continues through ossified radial epiphysis 232 b , radial growth plate cartilage 242 and ossified radius 232 .
- the acoustic path in the measurement of speed of sound is based upon the time of flight of the first acoustic signal to arrive at the receiver. As bone provides the fastest acoustic signal transmission time, the first received ultrasound signal indicates the path with the greatest ossification in a given area.
- the first ultrasound signal passing through region 276 will be markedly faster than the first ultrasound signal transmitted in an area along axis 152 , due to the greater amount of ossification in area 276 .
- the acoustic signal from transducers 122 and 120 , taken along region 276 is sent to a controlling unit 142 , and can be analyzed either locally, such as using a PC, or remotely through a computer network. Such analysis can be immediate or the information can be stored and analyzed at a later time. Further, the analysis may be based upon, for instance, spectral analysis (frequency domain) and/or temporal analysis (time domain), which are known to characterize biological tissue.
- spectral analysis frequency domain
- time domain time domain
- Dispersion is the slope at which the speed of sound changes with frequency.
- Dispersion is calculated from the dispersion graph, which is the curve of speed of sound versus frequency, calculated using the phase of each Fourier component of the signal and the distance. In a linear graph, the dispersion value is the slope of the graph and is given in units of m ⁇ s ⁇ 1 ⁇ MHz ⁇ 1 .
- Broadband ultrasound attenuation is the value of the slope at which attenuation changes with frequency. Broadband ultrasound attenuation is calculated from an attenuation graph, which is the curve of the logarithm of amplitude of the Ultrasound signal versus frequency, calculated using the amplitude of each Fourier component of the signal and the distance. If the attenuation graph is linear, the broadband ultrasound attenuation is the slope of the graph in units of db ⁇ cms ⁇ 1 ⁇ MHz ⁇ 1 .
- a database is established of children of different ages that contains the children's ultrasound parameters and their bone age as assessed by conventional radiograph.
- Statistical algorithms for instance, can be used to correlate the ultrasound parameters and the bone age.
- such an algorithm correlates speed of sound measured from the radius-ulna site, including the metaphysis and epiphysis, to radiographic bone age of the wrist and palm as read by the Greulich and Pyle (see Greulich and Pyle reference noted above) whereby acoustic signals provide the Bone Age (BA) utilizing equations such as:
- SOS speed of sound
- BA bone age
- two or more of the following parameters are obtained from the ultrasound signal are used to provide an estimation of bone age:
- the linear fit is calculated on several data points around the first half-amplitude level of the first half-wave of the ultrasound pulse. This type of processing requires that the first half wave of the ultrasound signal will not be saturated.
- all three parameters are formulated to calculate two new parameters.
- a1, a2, a3, a4 and b are constants in the model.
- this model uses different sets of constants for males and females.
- FIG. 4 is an exemplary embodiment of a graph 410 of bone age of hand 200 , assessed by ultrasonic multiparameter model as compared to standard X-ray for a population of girls represented by square marks and a population of boys, represented by triangular marks.
- a y axis 420 represents the bone age calculated by a multiparameter model on graph 440 and an x axis 422 represents the bone age calculated using X-ray on graph 440 , with line 424 demonstrating the best linear fit to the results.
- Graph 410 demonstrates an R 2 obtained for this correlation of 0.92 where R 2 is the square of the Pearson Correlation Coefficient.
- Wrist bones 100 comprise a Lunate 102 , a Scaphoid 104 , a Capitate 106 , a Trapezium 108 , a Trapezoid 110 , a Hamate 112 , a Triquetral 114 and a Pisiform 116 .
- Wrist bones 100 are cartilaginous at birth. Ossification, the process by which cartilage is replaced with bone, begins soon after birth in Capitate 106 . Ossification centers appear in bones on a regular pattern and predictably increase in size until skeletal maturation in the late teens. At age 6 . 5 , most of the bones of wrist 100 contain a region of ossification, represented by the shaded areas, while the clear areas represent regions that are still made up of cartilage. Ossified regions contain calcium that is relatively opaque to X-ray, whereas cartilage does not contain calcium and is not well differentiated over other soft tissue in the body on X-ray images. Thus the shaded areas would appear less dense on an X-ray image on a fairly uniform dense background.
- the ossified areas of wrist 100 of a 6.5-year-old male are: the Lunate 102 a , the Scaphoid 104 a , the Capitate 106 a , the Trapezium 108 a , the Trapezoid 110 a , the Hamate 112 a and the Triquetral 114 a .
- Pisiform 116 will not begin to ossify until the 9 th to 12 th year.
- velocity in wrist 100 is measured by transmitting an ultrasonic signal from transducer 120 to transducer 122 with the transducers 120 and 122 moved so that they transmit and receive along axis 152 .
- the distance between transducers 120 and 122 is fixed or known. For instance, the distance is measured by ruler 182 with a digital readout that is connected to ultrasound transducers 120 and 122 so that the average velocity through wrist 100 can be determined.
- the measured average velocity is assumed to be a function of the ratio between soft tissue (including cartilage) and ossified or partly ossified tissue. Therefore, the velocity of the acoustic signal is an indicator of bone age. By varying the location of measurement, a more exact estimate may be provided, for example, indicating which particular bones have ossified and which have not.
- a statistic of the velocity may be used for bone age assessment, for example, maximal, average or minimal velocity or any percentile or the velocity along a certain line in wrist 100 .
- these should be compared to the bone age A formula using one or more of these statistics, or a combination of these statistics is correlated with the bone age as assessed by standard age determination methods such as using the Greulich and Pyle method.
- the best correlative formula is utilized to provide a standard reference upon which bone age assessment is based.
- bone age can be measured in premature babies, from birth through infancy, childhood, puberty, and teenage years and beyond. For instance, in the second month following birth, Capitate 106 ossification center appears within area marked 106 a and an acoustic signal traveling through Capitate 106 will travel faster than a signal transmitted immediately following birth when Capitate 106 and all wrist 100 bones are cartilaginous. Even prior to the second month when Capitate ossification center 106 a appears, the density of parts of Capitate 106 may increase, indicating the amount of time that has passed since birth. Such increase may occur in other bones as well. Thus, by determining the acoustic velocity associated with ossification of infant wrist 100 , bone age can be estimated.
- Capitate ossification center 106 a continues to increase in size. At the end of the third month, the ossification center of Hamate 112 a appears and bone age can be measured based upon the amount of ossification in either bone or both bones simultaneously. Hamate 112 and Capitate 106 continue to ossify, so that when transducers 120 and 122 are put in a position so that the acoustic signal travels along dashed line 162 , it has a higher velocity in the ninth month and a still higher velocity in the second year, providing a measure of bone age.
- Such measurements soon after birth have application, for instance, in judging crucial maturation progress of premature babies who are born with fetal bone structure according to the month of gestation at which they were born. For example, a premature baby born after 4.5 months of gestation will not have bone ossification of Capitate 106 a until the equivalent of 9 months gestation plus one month.
- appearance of ossification in Capitate 106 a signifies a real age of one-month though the baby is 5.5 months postpartum.
- Ultrasound being non-ionizing as is a radiograph beam, is much more suited for making these measurements on pediatric subjects, particularly when multiple measurements are made.
- Triquetral ossification center 114 a appears in the third year, while Lunate ossification center 102 a appears during the fourth year, providing increased ossification density and faster speed of sound for measuring bone age along paths that include these structures.
- the ossification centers in the Scaphoid ( 104 a ), the trapezium ( 108 a ), and the trapezoid ( 110 a ), appear in the fourth year in females and the fifth year in males and continue to increase in size, providing further acceleration of acoustic signals for measuring bone age.
- Pisiform 116 ossification center 116 a appears in the ninth or tenth year in females and the twelfth year in males and continues to increase in size, providing further acceleration of acoustic signals for measuring bone age.
- ultrasound transmitter 120 and ultrasound receiver 122 remain stationary during the transmission and reception of ultrasound signals that pass through two or more bones of wrist 100 .
- Ultrasound transducers 120 and 122 are positioned in a specific relationship to wrist 100 , such as with both transducers in contact with the skin over one of more bones.
- ultrasound transmitter 120 is set in position near base of metacarpal 132
- Pisiform bone 116 and ultrasound receiver 122 is set in position near metacarpal 124 , over Scaphoid bone 104 .
- This position allows a signal to travel along axis 152 , through non-ossified Pisiform 116 , ossified 114 a and non-ossified portions of Triquetral 114 , Capitate non-ossified portion 106 , and ossified portion 104 a and non-ossified portions of Scaphoid 104 .
- the amount of skin and muscle tissue over these structures is small, so the acoustic signal transmitted along axis 152 has a velocity that reflects the amount of ossification along its path.
- ultrasound transmitter 120 is positioned on the dorsal surface of the hand, directly over Trapezoid 108 .
- Ultrasound receiver 122 is positioned on the ventral surface of the hand directly below Trapezoid 108 . The amount of ossification in the bone results in a specific acoustic velocity that provides a measure of bone age.
- Ultrasound transmitter 120 and ultrasound receiver 122 can be positioned over a variety of individual wrist structures to measure age. For instance, in the second month, over the dorsal and ventral surfaces of Capitate 106 and in the third month over the dorsal and ventral surfaces of Hamate 112 .
- both individual bones and multi-bone paths can be measured for determining a bone age single estimate using statistical analysis.
- the statistical analysis can be performed using discriminant analysis, binaric logistic regression, multinomial logistic regression or neural networks methods.
- a ratio between the acoustic signal velocity of two ossifying bones (or paths) is used to measure bone age. For instance, at the end of the fourth month, one month after Hamate 112 begins ossification; ultrasound transducers 120 and 122 are positioned over the dorsal and ventral surfaces over Hamate 112 to provide one measure of acoustic velocity. Transducers 120 and 122 are then placed over the dorsal and ventral surfaces of Capitate 106 to provide a second measure of acoustic velocity. These two acoustic signal velocities provide a ratio as a measure of bone age. Alternatively, two or more sets of transducers are used.
- transducers 120 and 122 are modified specifically for obtaining measurements from the wrist, with the transducer shape being, for instance, of a size and shape that matches the size and shape of wrist bones 100 .
- the gantry is designed to place the transducers against the skin automatically or manually. For instance one transducer remains fixed in position and the other transducer is moved in relation to the hand to place a standard pressure against the skin.
- Such a transducer may incorporate a calibrated pressure-sensitive sensor within the head to provide optimal accuracy of pressure.
- the hand is held in a bed that can move up or down in relation to the position of transducers 120 and 122 , so that a standard acoustic path is traced through wrist 100 .
- An imager may be provided to better locate the path 276 of acoustic signal, either signaling the operator to make positional adjustments of wrist 100 in relation to transducers 120 and 122 , or, optionally, making such adjustments automatically.
- measurement using ultrasound transducers 120 and 122 is repeated at a plurality of locations, for example, locations along a reference line 146 or in a radial fashion according to a reference 148 .
- Such measurement provides spatial acoustic information of one or more of a plurality of ossification centers from one or more of a plurality of reception and transmission points.
- Such acoustic velocities provide an acoustic map of wrist 100 for measuring bone age. Additionally or alternatively, two or more such acoustic velocities can be used to provide a ratio as a measure of bone age. Additionally or alternatively, two or more such acoustic velocities are averaged to provide a measure of bone age.
- spatial information is provided as a profile of velocities.
- a profile may be used, for example, to locate an area of minimum or maximum velocity, to locate an ossified or partially ossified area and/or to detect and/or identify abnormal ossification patterns.
- Such a profile can be generated, for example, using a computer system that maintains a record of the coordinates of the probes along with acoustic signal measurements.
- This data is displayed by the computer system, for example, as a 3D graph with 2 spatial axes and one speed of sound axis. Additionally or alternatively, the data provides a basis for directives on the further placement, movement or changes in signal frequency or signal analysis from probes 122 and 120 .
- computer system 142 automatically controls gantry 174 and transducers 120 and 122 .
- a first acoustic velocity is measured along one axis such as along line 146 and a second acoustic velocity is measured along another axis such as along line 154 and these values are averaged or placed in a ratio to provide a measure of bone age.
- transducer 120 and/or transducer 122 comprises an acoustic transmission grid that transmits multiple signals at specific time/space intervals.
- the multiple input and output points provide spatial measurements of the bone velocity and hence the extent and location of ossification and ossification centers.
- transducers 120 and 122 is a phased array beam scanning transducer.
- FIG. 2A illustrates ultrasound measurements of a portion of an arm 250 , a wrist 252 and digit structures 254 of an 11-year-old female in accordance with an exemplary embodiment of the present invention.
- the ossified portions of wrist 252 comprise a Scaphoid 104 a , Capitate 106 a , Trapezium 108 a , Trapezoid 110 a , Hamate 112 a and Triquetral 114 a , which show greater ossification than in 6.5 year-old male wrist 100 of FIG. 1.
- Pisiform 116 has begun to ossify in region 116 a .
- the acoustic signal from each of these bones will show increase in velocity over those of 6.5-year-old wrist 100 illustrated in FIG. 1.
- a ray 270 contains a distal phalanx 226 with a secondary growth plate 216 and epiphysis 226 a ; a middle phalanx 224 with a secondary growth plate 214 and an epiphysis 224 a ; a proximal phalanx 222 with a secondary growth plate 212 and epiphysis 222 a ; and a metacarpal 220 with a secondary growth plate 210 and epiphysis 220 a . Shown for reference, without ossification patterns, are metacarpals and associated rays of thumb, 124 , middle finger 128 , ring finger 130 and pinkie finger 132 .
- ultrasound transmitter 256 and ultrasound receiver 258 are moved along an axis 270 .
- the transducers 256 and 258 are moved in position so that the signal between them passes from the dorsal to ventral surfaces of portions of radius 232 , Scaphoid 104 , Trapezoid 110 , and through a second ray 218 of digit structures 254 .
- FIG. 2B shows an exemplary graph 260 of such ultrasound measurements taken along axis 270 . Where the y axis represents the speed of sound and the x axis represents the distance along ray 218 in centimeters.
- a graph area 266 demonstrates a decrease in velocity, reflecting a decrease in ossification in base of Scaphoid 104 .
- Transducers 256 and 258 pass through second ray 218 of digit structures 254 , passing a metacarpal 220 . As they pass a growth plate 210 , there is a corresponding decrease in velocity indicated by a graph region 210 a .
- Transducers 256 and 258 pass a proximal phalanx 222 , showing an increased velocity.
- Transducers 256 and 258 pass a middle phalanx 224 with an increase in velocity due to ossification. As they pass a growth plate 214 , there is a corresponding decrease in velocity indicated by a graph region 214 a . Transducers 256 and 258 pass a distal phalanx 226 with an increase in acoustic velocity. As they pass a growth plate 216 , there is a corresponding slow down in acoustic signal indicated by a graph region 216 a.
- Reduced velocity corresponding to the growth plate position will change until fusion when the growth plate becomes ossified.
- the length and shape of the bones will continue to remodel and grow in size until fusion is completed.
- the bone age is estimated by comparing the profile of velocity changes along multiple ossifying structures, such as those values contained in graph 260 to a database of standardized profiles. Alternatively or additionally, the bone age is estimated by comparing the velocity at one or more particular anatomical points (e.g., metacarpal 220 ) to a database of standardized velocities.
- final bone age is an average of a plurality of such bone ages. Alternatively, the bone age is based on an averaging of the velocity at the plurality of points. Additionally or alternatively, the bone age is based on a ratio between average acoustic velocity of a bone, such as metacarpal 220 , and a second bone, such as proximal phalanx 222 .
- the average acoustic velocity across the width of a bone, such as metacarpal 220 , and the acoustic velocity through the length of the same bone 220 are divided to provide a measure of bone age.
- the finger is bent and measurements are taken along the length of the proximal, middle and distal phalanges either separately or through two or more at along the same axis.
- a database is constructed of bone measurements.
- the database is used to analyze acoustic signals. Extraction of information from the database can be, for example, via matching received signals to similar signal and values contained in the database, or by interpolation.
- the information extracted from the database is used as an aid in diagnosis, e.g. estimating a bone age or a disease state. Additionally or alternatively, the extracted information is used for adult height prediction, e.g. final stature.
- height prediction can be made utilizing any of the known formulae and tables stated above. Such a derivation can be arrived at, for instance, by substituting the ultrasound bone age obtained by quantitative ultrasound in place of the bone age obtained by an X-ray film.
- the database associates one of more of the following items of information with acoustic velocity or acoustic signal parameters: sex, ethnic group, geographic location, nutrition, genetic inheritance.
- the database is created by acquiring acoustic signal measurement in a plurality of subjects, in a plurality of bones, at a plurality of growth stages, and in a plurality of populations.
- the bone age is estimated by comparing acoustic signals to known bone age measurement systems such as those based on X-ray, CT, MRI and Nuclear Medicine.
- the measurements may be used to train a neural network.
- the database is used as an aid in medical diagnosis and for monitoring. For instance, in conjunction with repeated bone age measurements the database monitors growth hormone therapy that is often administered to children who have short stature. Alternatively, the database can be used to diagnose such conditions as precocious puberty or delayed puberty so that appropriate hormone therapy can be administered.
- multiple measurements of a single subject are stored in a personal database.
- a personal database may be used, for example, for tracking changes in bone age over time or in diagnosing bone-age related disorders. Such disorders, for instance, can be related to malnutrition syndromes such as rickets and kwashiorkor. Additionally or alternatively, such a personal database may be used in tracking treatment of metabolic disease states such as diabetes. Additionally or alternatively, such a personal database may be used in diagnosis and/or treatment of endocrine-related growth disorders such as primary hyperparathyroidism or pituitary dwarfism.
- the measurement or measurements may be compared to data in a data-base that compares such measurements with one or more of adult bone length, density, thickness and resilience to predict one or more of adult bone length, density, thickness and resilience.
- transducer frequency is between 300 kHz to 5 MHz.
- a frequency closer to 5 MHz can be used as there is less attenuation of the ultrasound signal and a higher signal to noise ratio that improves the quality of signal.
- a frequency near 200 kHz is used, as a higher frequency may be attenuated excessively while going through several centimeters of tissue.
- a higher ultrasound frequency such as 10 MHz or even 50 MHz or higher may be used.
- higher or lower frequencies can be utilized.
- an existing osteoporosis measurement device such as the DBM Sonic 1200 ultrasound system by IGEA of Carpi (Mo), Italy, is modified to provide a bone age estimate.
- a database is reprogrammed with a table that associates acoustic velocities with bone ages, rather than with states of osteoporosis.
- a gantry 274 and ultrasound transducers 256 and 258 are electronically controlled by computer system 142 .
- parameters such as transducer positional changes and acoustic signal, are automatically effected by the computer system's analysis of the acoustic signal against an ultrasound bone age database for the purpose of obtaining optimal bone age data from ray 218 .
- Growth plates such as growth plate 212 , exhibit patterns of interdigitation with adjacent areas of bone 222 and this interdigitation may exhibit changes that can be used to forecast stature, metabolic bone health and age.
- multiple measurements with transducers 256 and 258 are taken radially while rotating the transducers around a growth plate 212 of bone 222 to provide a spatial map of growth plate 222 interdigitation.
- acoustic signals of interdigitation can be used for detecting signs of disease associated with growing bones.
- Osteogenic sarcoma for example, is a cancer that often begins from a growth plate such as in the knee. Ossification patterns in a growth plate affected by osteogenic sarcoma generally exhibit differences as compared to a normal growth plate, which differences can support early detection using ultrasound.
- ultrasonic measurements are used for screening for various conditions.
- Acoustic signals from the radial movement—a movement along line 270 —of transducers 256 and 258 are entered into computing means 142 and can be viewed on visual display 140 , for instance as a graph 260 . Additionally or alternatively, display 140 provides an image of the bone and accompanying data, or a database of information based on sex, race and/or other factors. Such information can be printed on printer 144 .
- FIG. 3A Illustrates ultrasound measurements of femora 310 and a pelvis 300 of a male age 14 using the echo method of transducer measurement in accordance with an embodiment of the present invention. Shown in FIG. 3A are an illium 320 , a sacrum 322 and an ishium 324 that are shaded to demonstrate ossification.
- An ultrasound transducer 344 contains one or more receivers and transmitters of acoustic signals arranged, for example, in a grid pattern. Probe 344 is shown sending and receiving acoustic signals to femur 310 to scan a femoral head 312 and a femoral trochanter 314 .
- Femoral head 312 develops a secondary ossification center 312 a at six months after birth and fuses with femur 310 at age 14 in females and 17 in males.
- ultrasound signal reflecting off a bone produces backscatter, picked up by the receiver in the echo method.
- the backscatter attenuates according to the density of material from which it is received, so that the bone content of a particular tissue is assessed.
- transducer 344 may be used to measure backscatter attenuation to determine the ossification level of the location of growth plates in femora 310 .
- greater trochanter 314 develops a secondary ossification center 314 a in the fourth year and fuses with femura 310 at age 14 in females and 17 in males.
- Ultrasound probe 344 measures bone age based on the backscatter attenuation of an acoustic signal in one or more of these structures to determine age.
- probe 344 is rotated along an arc 340 to provide an ossification ratio of femur 310 and pelvis 300 using backscatter attenuation, for instance, as a measure of bone age.
- backscatter attenuation for instance, as a measure of bone age.
- such measurement includes acoustic signal information of a pubic symphysis 328 .
- probe 344 containing both a receiver and transmitter, is rotated fully around pelvis 300 using the echo method of measurement, obtaining information from backscatter attenuation of ultrasound signal.
- a receiver probe is placed along line III, on one side of pubic symphysis 328 and a transmitter probe is placed at the opposite end of line III so that Broadband ultrasound attenuation can be measured from the pubic symphysis, 328 .
- a ratio between broadband ultrasound attenuation velocities of these structures, using two or more broadband ultrasound attenuation measurements and/or averages of broadband ultrasound attenuation within each structure can be used to measure bone age.
- probe 344 indicates and/or tracks a disease state such as slipped capitus femoris, where femoral head 312 is displaced on femur 310 . Additionally or alternatively, ultrasound probe 344 is positioned over pubic symphysis 328 to give acoustic velocity information. To determine SOS, two or more transducers should be contained in the probe or two or more separate transducers should be used.
- FIG. 3B illustrates a cross-sectional diagram of a pubic symphysis 328 along line III-III in a mature male and an associated graph 370 of acoustic backscatter attenuation.
- Pubic symphysis 328 is made up of an interpubic disc 352 that is fibrocartilaginous hyaline cartilage. One either side is hyaline cartilage 350 and 350 a and a portion of a pubic tubercle 356 and 356 a .
- Pubic tubercle 356 is ossified at age 14 .
- the signal backscatter increases representing the greater density of hyaline cartilage as shown in an area 378 in graph 370 .
- the signal backscatter increases even further representing the greater density of bone as represented by an area 380 of graph 370 .
- bone age can continue to be measured.
- Such measurements can also be used to indication the presence of a disease such as acromegaly, a diagnosis that can be confirmed by measuring the amount of somatotropin hormone production in the blood.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Rheumatology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
- This application is a continuation of PCT Application PCT/IL01/00683, filed Jul. 24, 2001, which designates the United States.
- The present invention relates to assessing bone age using ultrasound.
- Bone age assessment in growing subjects is a well-known diagnostic tool that is especially useful in predicting stature and/or growth problems in children, teenagers and adults.
- Many methods of assessing bone age are based upon radiographic analysis, such as the Greulich and Pyle (GP) method (Greulich WW, Pyle SI, Radiographic atlas of skeletal development of the hand and wrist, 2nd ed. Stanford Calif., Stanford University Press, 1959.) and the Tanner and Whitehouse (TW2) method [Tanner, J M, Whitehouse, R H, Marshall, W A, et al. “Assessment of skeletal maturity and prediction of adult height” (TW2 method). NY, Academic Press], both of which assess bone age by the radiographic presentation of the bones of the wrist and hand. In the Greulich and Pyle method a comparison is made between the child's radiograph and the corresponding standard in the Greulich and Pyle atlas. In the
TW2 method 20 bones in hand and wrist are scored according to their stage, thus producing a total score for which a skeletal age may be read directly from the tables. - These methods present problems of accessibility as X-ray units are often available only in secondary care centers such as hospitals. Additionally, ionizing radiation is undesirable in elective procedures, especially to children.
- Allessandro Castriota-Scanderbeg et al., in “Skeletal age assessment in children and young adults: Comparison between a newly developed sonographic method and conventional methods,” Skeletal Radiology 1998 27:271-277, propose a method for assessing skeletal age using ultrasound imaging measurements of the thickness of femoral head articular cartilage. In this method, a non-cartilage structure in an ultrasound image, an epiphysis, is linearly measured along a cross sectional plane without regard to structural aspects such as bone density. Further, this method is inaccurate and likely requires an imaging specialist for its administration, detracting from its cost-effectiveness.
- Chalana et al. U.S. Pat. No. 5,605,155 uses ultrasound images for the measurement of fetal head structure to predict fetal head size. Hechard Patrick, FR 2768322 uses X-rays or ultrasound to measure the thickness of the epiphysis and metaphysis of a bone to establish an index used for assessing a “Ratio of Residual Growth.” Holmberg, U.S. Pat. No. 6,135,960 proposes using ultrasound transducers placed in a Cartesian coordinate system for characterizing objects within the body. The above publications assess object boundaries utilizing dimensional ultrasound imaging techniques without regard to structural aspects such as bone density.
- Ultrasonic methods of bone density measurement are known. For instance, according to prior art ultrasound measurement systems, (for example,
WO 00/28316 and U.S. Pat. No. 5,564,423, the disclosure of which are incorporated herein by reference), ultrasound is used to determine the density of a non-cartilage osseous structure. - Thus use of backscatter attenuation to determine bone density is described by Wear, K A and Garra, B S, “Assessment of bone density using ultrasonic backscatter”,Ultrasound Med Biol, 1998 Jun; 24(5):689-95.
- Langton, et al. in “Quantitative Ultrasound”
Chapter 17, p. 311-312, measure bone's speed of sound and broadband ultrasound attenuation in children using the center of the posterior portion of the Calcaneal, specifically so the ultrasound signal does not pass through the Calcaneal growth plate. This measurement area begins ossification prenatally so that this method estimates bone's speed of sound and broadband ultrasound attenuation based upon non-cartilage properties. - An aspect of some embodiments of the invention relates to estimating bone age from acoustic signals of ossifying structures.
- In an exemplary embodiment of the present invention, bone age is estimated by measuring an acoustic velocity in cartilage structures that are in the process of ossifying wherein the velocity is expected to increase as a function of the ossification during the human maturation process. Such acoustic signal velocity, for example, is measured in primary ossification centers such as the bones of the wrist or secondary ossification centers such as the distal regions of the ulna and radius.
- Optionally, two or more acoustic signals of one or more ossification centers are measured to determine bone age. Additionally or alternatively, the ratio of acoustic signals between two or more ossification centers is used to measure bone age. Additionally or alternatively, acoustic signals from structures associated with ossifying structures are used to measure bone age, for example, the fibrocartilage of the pubic symphysis, skull suture ligaments and tooth and mandibular changes.
- Bone age may be used to predict adult stature or other aspects of the maturation process. Such predictions are based on bone age derived by other methods known in the art. Alternatively or additionally, tracking of ossification in bone is used to detect and/or track the progress of various disease states and/or disorders, with, for instance, a more accurate profile than X-ray evaluation due to its non-ionizing nature, allowing frequent monitoring without harm.
- In an exemplary embodiment, parameters other than acoustic velocity, such as broadband ultrasound attenuation (BUA) and dispersion of ultrasound signal, are used to estimate bone age, for example by correlating these parameters with the known BA assessment of a group of children. Additionally or alternatively, signals reflected from bone are used to measure bone age, for example, by measuring attenuation of backscatter intensity of the ultrasound signal.
- Optionally, especially for some types of measures and/or where object boundaries are unclear, a scanning or multi-beam measurement system may be used. Optionally, the acoustic signal provides a spatial measure, for example, indicating a profile of velocity along a bone axis or a radial profile of an ossification center.
- Different bones, different measures and/or different measurement systems may be used for different situations and/or for analyzing different bone ages or disease states.
- An aspect of some embodiments of the invention relates to using an existing osteoporosis measurement device, possibly with minimal changes, to assess bone age. In one example, a device designed for measuring osteoporosis in a finger is reprogrammed with a table that associates acoustic velocities with bone ages, rather than with states of osteoporosis.
- Further, the transducers for measuring acoustic signal are modified specifically for application to growth centers. It is noted that velocity limits used in osteoporosis measurement-devices are designed to obtain measurements from non-growth center areas.
- There is thus provided, in accordance with an exemplary embodiment of the invention a method for measuring bone age comprising:
- transmitting an acoustic energy into the body of a subject;
- receiving an acoustic signal from one or more structures including an ossification-actuated skeletal structure or a cranial structure that changes with age, responsive to said transmitted acoustic energy;
- analyzing the acoustic signal to determine at least one effect of said structure on said signal; and
- estimating the age of the structure from said determined effect.
- Optionally, said ossification-actuated skeletal structure comprises one or more areas undergoing ossification. Optionally, said ossification-actuated skeletal structure comprises one or more bones. Optionally, said ossification-actuated skeletal structure comprises one or more regions of cartilage. Optionally, said ossification-actuated skeletal structure comprises one or more regions of non-cartilage soft tissue. Optionally, said ossification-actuated skeletal structure comprises one or more regions of fibrocartilage.
- In an embodiment of the invention, said ossification-actuated skeletal structure comprises a region with one or more primary ossification centers. Optionally, said ossification-actuated skeletal structure comprises one or more of: the bones of the wrist, the bones of the palm, the bones of the tarsus, the mandible.
- In an embodiment of the invention, said ossification-actuated skeletal structure comprises a region with one or more secondary ossification centers. Optionally, said ossification-actuated skeletal structure contains an epiphysis. Optionally said ossification-actuated skeletal structure comprises a region of one or more of: an ulna, a radius, a femur, a bone of a ray of an extremity.
- Optionally, said receiving comprises utilizing two or more different acoustic signals to provide a measure of bone age. Optionally, said two or more acoustic signals are associated with the same bone. Optionally, said two or more acoustic signals are associated with paths in different bones. Optionally, said two or more acoustic signals are received from the same direction. Optionally, said two or more acoustic signals are received from the different directions.
- In an embodiment of the invention, said signal passes through said one or more structures including an ossification-actuated skeletal structure.
- In an embodiment of the invention, said signal echoes from said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said analysis of said signal is responsive to speed of sound from said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said analysis of said signal is responsive to broadband ultrasound attenuation from said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said analysis of said signal is responsive to dispersion of ultrasound from said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said analysis of said signal is performed, at least in part, in the frequency domain. Optionally, said analysis of said signal is performed, at least in part, in the time domain.
- In an embodiment of the invention, said analysis of said signal is responsive to attenuation of an ultrasound signal in said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said analysis is used to predict adult stature.
- In an embodiment of the invention, to provide an estimate of bone age, said analysis is compared to a database having correlation with one or more of: conventional radiographs, CT images, MRI images and Nuclear Medicine scans.
- In an embodiment of the invention, said transmitting is by a scanning acoustic signal transmitter.
- In an embodiment of the invention, said transmitting is by a multi-beam acoustic signal transmitter.
- Optionally, said receiving provides two or more acoustic signal measures along an axis of said one or more structures including an ossification-actuated skeletal structure.
- Optionally, said receiving provides two or more acoustic signal measures radially around said one or more structures including an ossification-actuated skeletal structure.
- In an embodiment of the invention, said analysis is correlated with a known bone age measurement system.
- In an embodiment of the invention, said analysis is responsive to a formula providing a correlation with a known bone age measurement system. Optionally, is responsive to at least one of speed of sound, broadband ultrasound attenuation, scattering and dispersion of acoustic signal through or from an ossification activated skeletal structure. Optionally, an estimate of bone age is responsive to time of flight of an acoustic signal between two transducers, with said ossification activated skeletal structure being situated intermediate to said transducers.
- Optionally, separate formulas are used to correlate known bone age data with acoustic signals from males and females.
- In an embodiment of the invention, said acoustic information is constructed into a database of bone age measurements. Optionally, said database is arranged according to one or more of: sex, ethnic group, geographic location, nutrition and general inheritance. Optionally, said database includes two or more measurements of one or more of said one or more structures including an ossification-actuated skeletal structure. Optionally, said database includes one or more measurements of two or more growth stages from said one or more structures including an ossification-actuated skeletal structure Optionally, said database includes one or more measurements of said one or more structures including an ossification-actuated skeletal structure in two or more populations.
- Optionally, said received signals are compared to similar signals in a database to predict one or more of predict one or more of adult bone length, density, thickness and resilience and adult stature. Optionally, said received signals are compared to similar signals in a database to indicate a bone-growth related disorder. Optionally, said received signals are compared to similar signals in a database to track the progress of a bonegrowth related disorder. Optionally, said received signals are compared to similar signals in a database to track hormone therapy in a growth stature disorder. Optionally, said received signals are compared to similar signals in a database to indicate one or more growth-plate related disease states, including osteogenic sarcoma, slipped growth plate, premature arrest of growth plate growth and inflammation of growth plate.
- Optionally, two or more acoustic measurements are made on a single subject and entered into said database. Optionally, said two or more acoustic measurements are compared to track one or more growth-related disorders, including precocious puberty, delayed puberty, rickets, kwashiorkor, hypoparathyroidism, pituitary dwarfism and diabetes.
- Optionally, said two or more acoustic measurements are compared to track treatment of one or more growth-related disorders, including precocious puberty, delayed puberty, rickets, kwashiorkor, hypoparathyroidism, pituitary dwarfism and diabetes.
- There is further provided, in accordance with an exemplary embodiment of the invention, an apparatus for estimating bone age comprising:
- an acoustic transmitter and an acoustic receiver positioned on either side of one or more structures including an ossification-actuated skeletal structure;
- an electronic moveable gantry that adjusts the position of said acoustic transmitter and said acoustic receiver in relation to said ossification-actuated structure;
- a computer system that performs one or more functions of:
- positioning of said moveable gantry;
- controlling acoustic signals transmitted by said acoustic transmitter;
- receiving acoustic signals from said receiver responsive to said transmitted signals; and
- estimating said bone age responsive to said received signals.
- In an embodiment of the invention, said apparatus transmits and receives one or more acoustic signals linearly along an axis through said ossification-actuated structure.
- In an embodiment of the invention, said apparatus transmits and receives one or more acoustic signals radially around an axis through said ossification-actuated structure.
- Optionally, said computer system controls said acoustic signal transmitter to provide an acoustic signal appropriate for said ossification-actuated structure.
- Optionally, said computer system estimates said bone age responsive to one of more of broadband ultrasound attenuation, acoustic backscatter, dispersion of acoustic signal and speed of sound in said ossification-actuated structure.
- Optionally, said computer system uses an imager to control the position of said acoustic signal receiver and said acoustic signal transmitter.
- Optionally, said computer system contains a visual display to provide information on said bone age. Optionally, said visual display comprises a graph.
- Optionally, said computer system is comprised in a computer network.
- Optionally, said computer system comprises a neural network.
- Optionally, said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures including an ossification-actuated skeletal structure to provide an estimate of bone age.
- Optionally, said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures, including an ossification-actuated skeletal structure to predict stature. Optionally, said computer system compares said received acoustic signal to a database containing information of one or more acoustic signals from one or more structures including an ossification-actuated skeletal structure to indicate, track or follow treatment of one or more of: a bone-growth related disorder, a growth plate disorder, and a growth related disorder.
- Exemplary non-limiting embodiments of the present invention described in the following description, read with reference to the figures attached hereto. In the figures, identical and similar structures, elements or parts thereof that appear in more than one figure are generally labeled with the same or similar references in the figures in which they appear. Dimensions of components and features shown in the figures are chosen primarily for convenience and clarity of presentation and are not necessarily to scale. The attached figures are:
- FIG. 1 schematically illustrates ultrasound measurements of ossifying wrist and arm bones of a 6.5-year-old male, in accordance with an embodiment of the present invention;
- FIG. 2A schematically illustrates ultrasound measurements of a portion of ossifying arm, wrist and hand bones of an 11 year old female in accordance with an embodiment of the present invention;
- FIG. 2B schematically illustrates a typical graph of ultrasound measurements associated with a path along the bones of FIG. 2A, in accordance with an embodiment of the present invention;
- FIG. 3A schematically illustrates ultrasound measurements of ossifying femora and pelvis of a male,
age 14, shown in a radiographic representation in accordance with an embodiment of the present invention; and - FIG. 3B schematically illustrates a pubic symphysis, shown in cross-section and a typical graph of its ultrasound measurement in accordance with an embodiment of the present invention.
- FIG. 4 is a graph of Bone Age of the hand assessed by ultrasonic multi parameter model as compared to standard X-ray.
- FIG. 1 illustrates ultrasound measurements of a portion of an
arm 150 and awrist 100 of a 6.5-year-old male in accordance with an exemplary embodiment of the present invention. Each bone ofwrist 100 exhibits an ossification center, represented by a crosshatched area, that begins from a central portion and progresses toward the periphery and is referred to as a “primary ossification center”. Bones with only primary ossification centers continue to show changes in acoustic velocity due to increased ossification, increased bone size and/or changes in bone shape through the 15th-19th years. Some bones have additional ossification centers that appear as a strip of cartilage between ossified bone sections, referred to as a “secondary ossification center.” - A
radius 232, for instance, contains a secondary ossification center consisting of acartilaginous growth plate 242 that is situated between anepiphysis 232 b and a metaphysis 232 a, both of which are ossified.Growth plate 242 first appears in the first postnatal year and ossifies in the seventeenth year in females and in the nineteenth year in males. When a secondary ossification center ossifies, at a later age,epiphysis 232 b and metaphysis 232 a fuse together withgrowth plate 242 that ossifies as well. After the stage when fusion occurs,bone 232 is fully ossified lacking the bands of cartilage and hence is radio-opaque throughout. -
Ulna 244 similarly contains agrowth plate 244 a that appears in the fifth year in females and in the sixth year in males.Ulnar growth plate 244 a ossifies in the seventeenth year in females and eighteenth year in males. - In an exemplary embodiment, two transducers, such as a
receiver 122 and atransmitter 120, are held in agantry 174 and placed so that they are in contact with the skin overulna 244 andradius 232. The thickness of soft tissue covering the bone is measured or the effect of soft tissue is optionally ignored.Transducers region 276, to obtain acoustic signals fromulna 244 andradius 232 and their respective ossification centers. This method of ultrasound transmission is referred to herein as a “through” method. In an exemplary embodiment, the frequency oftransducers transducers digital readout 182. - Generally, the transducers are controlled and the received signals processed by a controller142 (e.g., a computer), which optionally includes a
display 140 and/or aprinter 144. - The acoustic signal from a growth plate, such as
ulna growth plate 244 a, is markedly different than that of bone, such asbone 244. One parameter that can be measured from the obtained acoustic signal is the speed of sound. The distance between transducers is evaluated bydigital ruler 182 that is connected totransducers ulnar cartilage 244 a is about 1700 meters per second (m/s) while speed of sound inbone 244 is about 2000-4500 m/s, depending on probe localization and age. The transmission time of a signal varies based upon the composition of the structures in the path, which changes with age. - An ultrasound signal passing along
region 276 will pass throughossified ulna 244, ulnargrowth plate cartilage 244 a, andulnar epiphysis 248 that is also ossified. It continues through ossifiedradial epiphysis 232 b, radialgrowth plate cartilage 242 andossified radius 232. - The acoustic path in the measurement of speed of sound, for example, is based upon the time of flight of the first acoustic signal to arrive at the receiver. As bone provides the fastest acoustic signal transmission time, the first received ultrasound signal indicates the path with the greatest ossification in a given area. The first ultrasound signal passing through
region 276, for example, will be markedly faster than the first ultrasound signal transmitted in an area alongaxis 152, due to the greater amount of ossification inarea 276. - The acoustic signal from
transducers region 276 is sent to a controllingunit 142, and can be analyzed either locally, such as using a PC, or remotely through a computer network. Such analysis can be immediate or the information can be stored and analyzed at a later time. Further, the analysis may be based upon, for instance, spectral analysis (frequency domain) and/or temporal analysis (time domain), which are known to characterize biological tissue. - Other parameters that can be obtained from the acoustic signal include, dispersion and Broadband Ultrasound Attenuation. Dispersion is the slope at which the speed of sound changes with frequency. Dispersion is calculated from the dispersion graph, which is the curve of speed of sound versus frequency, calculated using the phase of each Fourier component of the signal and the distance. In a linear graph, the dispersion value is the slope of the graph and is given in units of m·s−1 ·MHz−1.
- Broadband ultrasound attenuation is the value of the slope at which attenuation changes with frequency. Broadband ultrasound attenuation is calculated from an attenuation graph, which is the curve of the logarithm of amplitude of the Ultrasound signal versus frequency, calculated using the amplitude of each Fourier component of the signal and the distance. If the attenuation graph is linear, the broadband ultrasound attenuation is the slope of the graph in units of db·cms−1·MHz−1.
- In an exemplary embodiment of the invention, a database is established of children of different ages that contains the children's ultrasound parameters and their bone age as assessed by conventional radiograph. Statistical algorithms, for instance, can be used to correlate the ultrasound parameters and the bone age.
- In an exemplary embodiment, such an algorithm correlates speed of sound measured from the radius-ulna site, including the metaphysis and epiphysis, to radiographic bone age of the wrist and palm as read by the Greulich and Pyle (see Greulich and Pyle reference noted above) whereby acoustic signals provide the Bone Age (BA) utilizing equations such as:
- BA=(SOS−1566)/27.9 for females [Eq. 1]
- BA=(SOS−1655)/15.6 for males, [Eq. 2]
- Where speed of sound (SOS) is in units of m/s and bone age (BA) is given in years. These formulae, for instance, are valid for an age range of 4-18 years.
- Multi Parameter Bone Age Assessment
- In an exemplary embodiment of the invention two or more of the following parameters are obtained from the ultrasound signal are used to provide an estimation of bone age: The distance (D), in mm, between the transmitter and receiver probes, the time (T1) of the arrival of the pulse, calculated by extrapolating a linear fit, as explained below, to the zero-signal level and the time (T2) of receipt of the pulse front, calculated by interpolating the same linear fit as in TI to the signal half-amplitude level.
- The linear fit is calculated on several data points around the first half-amplitude level of the first half-wave of the ultrasound pulse. This type of processing requires that the first half wave of the ultrasound signal will not be saturated.
- In an exemplary embodiment, all three parameters are formulated to calculate two new parameters. C1=D/T1 and C2=D/T2, and uses a linear combination of T1, T2, C1 and C2 as an estimator to the bone age:
- BA=a1·C1+a2·C2+a3·T1+a4·T2+b, [Eq. 3]
- Where a1, a2, a3, a4 and b are constants in the model. Optionally, this model uses different sets of constants for males and females.
- For males: a1=0.35; a2=−0.32; a3=22.7; a4=−22.5; b=−56.8;
- For females: a1=0.14; a2=−0.13; a3=10.1; a4=−9.8; b=−31.0;
- FIG. 4 is an exemplary embodiment of a graph410 of bone age of
hand 200, assessed by ultrasonic multiparameter model as compared to standard X-ray for a population of girls represented by square marks and a population of boys, represented by triangular marks. Ay axis 420 represents the bone age calculated by a multiparameter model ongraph 440 and anx axis 422 represents the bone age calculated using X-ray ongraph 440, withline 424 demonstrating the best linear fit to the results. Graph 410 demonstrates an R2 obtained for this correlation of 0.92 where R2 is the square of the Pearson Correlation Coefficient. - Bone Age Estimated From Other Sites
-
Wrist bones 100 comprise aLunate 102, aScaphoid 104, aCapitate 106, aTrapezium 108, aTrapezoid 110, aHamate 112, aTriquetral 114 and aPisiform 116. -
Wrist bones 100 are cartilaginous at birth. Ossification, the process by which cartilage is replaced with bone, begins soon after birth inCapitate 106. Ossification centers appear in bones on a regular pattern and predictably increase in size until skeletal maturation in the late teens. At age 6.5, most of the bones ofwrist 100 contain a region of ossification, represented by the shaded areas, while the clear areas represent regions that are still made up of cartilage. Ossified regions contain calcium that is relatively opaque to X-ray, whereas cartilage does not contain calcium and is not well differentiated over other soft tissue in the body on X-ray images. Thus the shaded areas would appear less dense on an X-ray image on a fairly uniform dense background. - The ossified areas of
wrist 100 of a 6.5-year-old male are: theLunate 102 a, theScaphoid 104 a, theCapitate 106 a, theTrapezium 108 a, theTrapezoid 110 a, theHamate 112 a and theTriquetral 114 a.Pisiform 116 will not begin to ossify until the 9th to 12th year. - In an exemplary embodiment of the invention, velocity in
wrist 100 is measured by transmitting an ultrasonic signal fromtransducer 120 totransducer 122 with thetransducers axis 152. In an exemplary embodiment of the invention, the distance betweentransducers ruler 182 with a digital readout that is connected toultrasound transducers wrist 100 can be determined. - As indicated, the measured average velocity is assumed to be a function of the ratio between soft tissue (including cartilage) and ossified or partly ossified tissue. Therefore, the velocity of the acoustic signal is an indicator of bone age. By varying the location of measurement, a more exact estimate may be provided, for example, indicating which particular bones have ossified and which have not.
- Alternatively or additionally, a statistic of the velocity may be used for bone age assessment, for example, maximal, average or minimal velocity or any percentile or the velocity along a certain line in
wrist 100. In an exemplary embodiment these should be compared to the bone age A formula using one or more of these statistics, or a combination of these statistics is correlated with the bone age as assessed by standard age determination methods such as using the Greulich and Pyle method. Ultimately, the best correlative formula is utilized to provide a standard reference upon which bone age assessment is based. - The above-defined methods are in distinct contrast to those of Langton, et al. (“Quantitative Ultrasound”
Chapter 17, p. 311-312) who use a portion of non-cartilage Calcaneus as a basis for measuring bone's speed of sound and broadband ultrasound attenuation as noted above. There is no mention of bone age measurement in this reference. - Assessing Bone Age Using Ultrasound Along A Defined Path
- Using various methods, as described below, bone age can be measured in premature babies, from birth through infancy, childhood, puberty, and teenage years and beyond. For instance, in the second month following birth,
Capitate 106 ossification center appears within area marked 106 a and an acoustic signal traveling throughCapitate 106 will travel faster than a signal transmitted immediately following birth whenCapitate 106 and allwrist 100 bones are cartilaginous. Even prior to the second month whenCapitate ossification center 106 a appears, the density of parts ofCapitate 106 may increase, indicating the amount of time that has passed since birth. Such increase may occur in other bones as well. Thus, by determining the acoustic velocity associated with ossification ofinfant wrist 100, bone age can be estimated. -
Capitate ossification center 106 a continues to increase in size. At the end of the third month, the ossification center ofHamate 112 a appears and bone age can be measured based upon the amount of ossification in either bone or both bones simultaneously.Hamate 112 andCapitate 106 continue to ossify, so that whentransducers line 162, it has a higher velocity in the ninth month and a still higher velocity in the second year, providing a measure of bone age. - Such measurements soon after birth have application, for instance, in judging crucial maturation progress of premature babies who are born with fetal bone structure according to the month of gestation at which they were born. For example, a premature baby born after 4.5 months of gestation will not have bone ossification of
Capitate 106 a until the equivalent of 9 months gestation plus one month. By measuring premature baby bone age, appearance of ossification inCapitate 106 a, signifies a real age of one-month though the baby is 5.5 months postpartum. Ultrasound, being non-ionizing as is a radiograph beam, is much more suited for making these measurements on pediatric subjects, particularly when multiple measurements are made. - With time, other bones develop ossification centers. For instance,
Triquetral ossification center 114 a appears in the third year, whileLunate ossification center 102 a appears during the fourth year, providing increased ossification density and faster speed of sound for measuring bone age along paths that include these structures. - The ossification centers in the Scaphoid (104 a), the trapezium (108 a), and the trapezoid (110 a), appear in the fourth year in females and the fifth year in males and continue to increase in size, providing further acceleration of acoustic signals for measuring bone age.
Pisiform 116 ossification center 116 a appears in the ninth or tenth year in females and the twelfth year in males and continues to increase in size, providing further acceleration of acoustic signals for measuring bone age. - Shown for reference only, without ossification, are the bases of metacarpals, associated with Thumb,124,
index finger 126,middle finger 128,ring finger 130 andpinkie finger 132. - In one embodiment of the invention,
ultrasound transmitter 120 andultrasound receiver 122 remain stationary during the transmission and reception of ultrasound signals that pass through two or more bones ofwrist 100.Ultrasound transducers wrist 100, such as with both transducers in contact with the skin over one of more bones. In the present set-up, for instance,ultrasound transmitter 120 is set in position near base ofmetacarpal 132, atPisiform bone 116 andultrasound receiver 122 is set in position nearmetacarpal 124, overScaphoid bone 104. This position allows a signal to travel alongaxis 152, throughnon-ossified Pisiform 116, ossified 114 a and non-ossified portions ofTriquetral 114, Capitatenon-ossified portion 106, andossified portion 104 a and non-ossified portions ofScaphoid 104. The amount of skin and muscle tissue over these structures is small, so the acoustic signal transmitted alongaxis 152 has a velocity that reflects the amount of ossification along its path. - Additionally or alternatively, measurement is made of a single wrist bone. For instance,
ultrasound transmitter 120 is positioned on the dorsal surface of the hand, directly overTrapezoid 108.Ultrasound receiver 122 is positioned on the ventral surface of the hand directly belowTrapezoid 108. The amount of ossification in the bone results in a specific acoustic velocity that provides a measure of bone age. -
Ultrasound transmitter 120 andultrasound receiver 122 can be positioned over a variety of individual wrist structures to measure age. For instance, in the second month, over the dorsal and ventral surfaces ofCapitate 106 and in the third month over the dorsal and ventral surfaces ofHamate 112. In addition, both individual bones and multi-bone paths can be measured for determining a bone age single estimate using statistical analysis. The statistical analysis can be performed using discriminant analysis, binaric logistic regression, multinomial logistic regression or neural networks methods. - Alternatively, a ratio between the acoustic signal velocity of two ossifying bones (or paths) is used to measure bone age. For instance, at the end of the fourth month, one month after
Hamate 112 begins ossification;ultrasound transducers Hamate 112 to provide one measure of acoustic velocity.Transducers Capitate 106 to provide a second measure of acoustic velocity. These two acoustic signal velocities provide a ratio as a measure of bone age. Alternatively, two or more sets of transducers are used. - In an
exemplary embodiment transducers wrist bones 100. Optionally, the gantry is designed to place the transducers against the skin automatically or manually. For instance one transducer remains fixed in position and the other transducer is moved in relation to the hand to place a standard pressure against the skin. Such a transducer may incorporate a calibrated pressure-sensitive sensor within the head to provide optimal accuracy of pressure. Optionally, the hand is held in a bed that can move up or down in relation to the position oftransducers wrist 100. An imager may be provided to better locate thepath 276 of acoustic signal, either signaling the operator to make positional adjustments ofwrist 100 in relation totransducers - Optionally, measurement using
ultrasound transducers reference line 146 or in a radial fashion according to areference 148. Such measurement provides spatial acoustic information of one or more of a plurality of ossification centers from one or more of a plurality of reception and transmission points. Such acoustic velocities provide an acoustic map ofwrist 100 for measuring bone age. Additionally or alternatively, two or more such acoustic velocities can be used to provide a ratio as a measure of bone age. Additionally or alternatively, two or more such acoustic velocities are averaged to provide a measure of bone age. - In an exemplary embodiment of the invention, spatial information is provided as a profile of velocities. Such a profile may be used, for example, to locate an area of minimum or maximum velocity, to locate an ossified or partially ossified area and/or to detect and/or identify abnormal ossification patterns. Such a profile can be generated, for example, using a computer system that maintains a record of the coordinates of the probes along with acoustic signal measurements. This data is displayed by the computer system, for example, as a 3D graph with 2 spatial axes and one speed of sound axis. Additionally or alternatively, the data provides a basis for directives on the further placement, movement or changes in signal frequency or signal analysis from
probes computer system 142 automatically controlsgantry 174 andtransducers - Additionally or alternatively, a first acoustic velocity is measured along one axis such as along
line 146 and a second acoustic velocity is measured along another axis such as alongline 154 and these values are averaged or placed in a ratio to provide a measure of bone age. - Additionally or alternatively to moving transducers,
transducer 120 and/or transducer 122 (e.g., receiver and/or transmitter) comprises an acoustic transmission grid that transmits multiple signals at specific time/space intervals. The multiple input and output points provide spatial measurements of the bone velocity and hence the extent and location of ossification and ossification centers. Alternatively or additionally, one or both oftransducers - FIG. 2A illustrates ultrasound measurements of a portion of an
arm 250, awrist 252 anddigit structures 254 of an 11-year-old female in accordance with an exemplary embodiment of the present invention. The ossified portions ofwrist 252, comprise a Scaphoid 104 a,Capitate 106 a,Trapezium 108 a,Trapezoid 110 a,Hamate 112 a andTriquetral 114 a, which show greater ossification than in 6.5 year-oldmale wrist 100 of FIG. 1. Additionally,Pisiform 116 has begun to ossify in region 116 a. Corresponding to these changes, the acoustic signal from each of these bones will show increase in velocity over those of 6.5-year-old wrist 100 illustrated in FIG. 1. - Secondary growth plates are also seen in the bones of the rays of each digit. A
ray 270, for instance, contains adistal phalanx 226 with asecondary growth plate 216 and epiphysis 226 a; a middle phalanx 224 with asecondary growth plate 214 and anepiphysis 224 a; aproximal phalanx 222 with a secondary growth plate 212 and epiphysis 222 a; and a metacarpal 220 with a secondary growth plate 210 andepiphysis 220 a. Shown for reference, without ossification patterns, are metacarpals and associated rays of thumb, 124,middle finger 128,ring finger 130 andpinkie finger 132. - In an exemplary embodiment,
ultrasound transmitter 256 andultrasound receiver 258 are moved along anaxis 270. Thetransducers radius 232,Scaphoid 104,Trapezoid 110, and through asecond ray 218 ofdigit structures 254. FIG. 2B shows anexemplary graph 260 of such ultrasound measurements taken alongaxis 270. Where the y axis represents the speed of sound and the x axis represents the distance alongray 218 in centimeters. - As
ultrasound transducers radius 232, the velocity of the acoustic signal, shown ingraph portion 262, decreases. Astransducers pass epiphysis 232 b, the acoustic velocity increases, as shown in agraph section 264, reflecting the presence of increased ossification within the cartilage. However, as the ossification is not complete, signal 264 does not reflect as high a velocity as prior tosignal area 262. - As
transducers Trapezoid 110 alongaxis 270, agraph area 266 demonstrates a decrease in velocity, reflecting a decrease in ossification in base ofScaphoid 104.Transducers second ray 218 ofdigit structures 254, passing ametacarpal 220. As they pass a growth plate 210, there is a corresponding decrease in velocity indicated by agraph region 210 a.Transducers proximal phalanx 222, showing an increased velocity. As they pass a growth plate 212, there is a corresponding decrease in velocity indicated by a graph region 212 a.Transducers growth plate 214, there is a corresponding decrease in velocity indicated by agraph region 214 a.Transducers distal phalanx 226 with an increase in acoustic velocity. As they pass agrowth plate 216, there is a corresponding slow down in acoustic signal indicated by agraph region 216 a. - Reduced velocity corresponding to the growth plate position will change until fusion when the growth plate becomes ossified. Typically, the length and shape of the bones will continue to remodel and grow in size until fusion is completed.
- In an exemplary embodiment of the invention, the bone age is estimated by comparing the profile of velocity changes along multiple ossifying structures, such as those values contained in
graph 260 to a database of standardized profiles. Alternatively or additionally, the bone age is estimated by comparing the velocity at one or more particular anatomical points (e.g., metacarpal 220) to a database of standardized velocities. Optionally, final bone age is an average of a plurality of such bone ages. Alternatively, the bone age is based on an averaging of the velocity at the plurality of points. Additionally or alternatively, the bone age is based on a ratio between average acoustic velocity of a bone, such asmetacarpal 220, and a second bone, such asproximal phalanx 222. - In another example, the average acoustic velocity across the width of a bone, such as
metacarpal 220, and the acoustic velocity through the length of thesame bone 220, are divided to provide a measure of bone age. To provide a measure along the bone length, for example, the finger is bent and measurements are taken along the length of the proximal, middle and distal phalanges either separately or through two or more at along the same axis. - In an exemplary embodiment a database is constructed of bone measurements. The database is used to analyze acoustic signals. Extraction of information from the database can be, for example, via matching received signals to similar signal and values contained in the database, or by interpolation.
- In an exemplary embodiment, the information extracted from the database is used as an aid in diagnosis, e.g. estimating a bone age or a disease state. Additionally or alternatively, the extracted information is used for adult height prediction, e.g. final stature. In an exemplary embodiment, height prediction can be made utilizing any of the known formulae and tables stated above. Such a derivation can be arrived at, for instance, by substituting the ultrasound bone age obtained by quantitative ultrasound in place of the bone age obtained by an X-ray film.
- In an exemplary embodiment, the database associates one of more of the following items of information with acoustic velocity or acoustic signal parameters: sex, ethnic group, geographic location, nutrition, genetic inheritance. In an exemplary embodiment, the database is created by acquiring acoustic signal measurement in a plurality of subjects, in a plurality of bones, at a plurality of growth stages, and in a plurality of populations. Optionally, the bone age is estimated by comparing acoustic signals to known bone age measurement systems such as those based on X-ray, CT, MRI and Nuclear Medicine. Alternatively to creating a database, the measurements may be used to train a neural network.
- In an exemplary embodiment, the database is used as an aid in medical diagnosis and for monitoring. For instance, in conjunction with repeated bone age measurements the database monitors growth hormone therapy that is often administered to children who have short stature. Alternatively, the database can be used to diagnose such conditions as precocious puberty or delayed puberty so that appropriate hormone therapy can be administered. In an exemplary embodiment, multiple measurements of a single subject are stored in a personal database. Such a personal database may be used, for example, for tracking changes in bone age over time or in diagnosing bone-age related disorders. Such disorders, for instance, can be related to malnutrition syndromes such as rickets and kwashiorkor. Additionally or alternatively, such a personal database may be used in tracking treatment of metabolic disease states such as diabetes. Additionally or alternatively, such a personal database may be used in diagnosis and/or treatment of endocrine-related growth disorders such as primary hyperparathyroidism or pituitary dwarfism.
- Alternatively or additionally, the measurement or measurements may be compared to data in a data-base that compares such measurements with one or more of adult bone length, density, thickness and resilience to predict one or more of adult bone length, density, thickness and resilience.
- In an exemplary embodiment, transducer frequency is between 300 kHz to 5 MHz. When the structures being measured are more superficial, as in
digit structures 254, a frequency closer to 5 MHz can be used as there is less attenuation of the ultrasound signal and a higher signal to noise ratio that improves the quality of signal. Along the mid shaft of the femur, for instance, a frequency near 200 kHz is used, as a higher frequency may be attenuated excessively while going through several centimeters of tissue. - When greater clarity of microstructure is desired, such as in spatial radial measurement of the bone growth center, a higher ultrasound frequency, such as 10 MHz or even 50 MHz or higher may be used. With alternative transducer design, higher or lower frequencies can be utilized.
- In an exemplary embodiment, an existing osteoporosis measurement device, such as the DBM Sonic 1200 ultrasound system by IGEA of Carpi (Mo), Italy, is modified to provide a bone age estimate. In an exemplary embodiment, a database is reprogrammed with a table that associates acoustic velocities with bone ages, rather than with states of osteoporosis. Additionally or alternatively, a gantry274 and
ultrasound transducers computer system 142. In an exemplary embodiment, parameters such as transducer positional changes and acoustic signal, are automatically effected by the computer system's analysis of the acoustic signal against an ultrasound bone age database for the purpose of obtaining optimal bone age data fromray 218. - Growth plates, such as growth plate212, exhibit patterns of interdigitation with adjacent areas of
bone 222 and this interdigitation may exhibit changes that can be used to forecast stature, metabolic bone health and age. In an exemplary embodiment, multiple measurements withtransducers bone 222 to provide a spatial map ofgrowth plate 222 interdigitation. - Additionally or alternatively, acoustic signals of interdigitation can be used for detecting signs of disease associated with growing bones. Osteogenic sarcoma, for example, is a cancer that often begins from a growth plate such as in the knee. Ossification patterns in a growth plate affected by osteogenic sarcoma generally exhibit differences as compared to a normal growth plate, which differences can support early detection using ultrasound. Optionally, ultrasonic measurements are used for screening for various conditions.
- Acoustic signals from the radial movement—a movement along
line 270—oftransducers visual display 140, for instance as agraph 260. Additionally or alternatively,display 140 provides an image of the bone and accompanying data, or a database of information based on sex, race and/or other factors. Such information can be printed onprinter 144. - FIG. 3A Illustrates ultrasound measurements of
femora 310 and apelvis 300 of amale age 14 using the echo method of transducer measurement in accordance with an embodiment of the present invention. Shown in FIG. 3A are anillium 320, asacrum 322 and anishium 324 that are shaded to demonstrate ossification. Anultrasound transducer 344 contains one or more receivers and transmitters of acoustic signals arranged, for example, in a grid pattern.Probe 344 is shown sending and receiving acoustic signals tofemur 310 to scan afemoral head 312 and afemoral trochanter 314.Femoral head 312 develops asecondary ossification center 312 a at six months after birth and fuses withfemur 310 atage 14 in females and 17 in males. In backscatter attenuation measurements, ultrasound signal reflecting off a bone produces backscatter, picked up by the receiver in the echo method. The backscatter attenuates according to the density of material from which it is received, so that the bone content of a particular tissue is assessed. Hence,transducer 344 may be used to measure backscatter attenuation to determine the ossification level of the location of growth plates infemora 310. For example,greater trochanter 314 develops asecondary ossification center 314 a in the fourth year and fuses withfemura 310 atage 14 in females and 17 in males.Ultrasound probe 344 measures bone age based on the backscatter attenuation of an acoustic signal in one or more of these structures to determine age. - Additionally or alternatively,
probe 344 is rotated along anarc 340 to provide an ossification ratio offemur 310 andpelvis 300 using backscatter attenuation, for instance, as a measure of bone age. Optionally such measurement includes acoustic signal information of apubic symphysis 328. Additionally or alternatively,probe 344, containing both a receiver and transmitter, is rotated fully aroundpelvis 300 using the echo method of measurement, obtaining information from backscatter attenuation of ultrasound signal. - Additionally or alternatively, a receiver probe is placed along line III, on one side of
pubic symphysis 328 and a transmitter probe is placed at the opposite end of line III so that Broadband ultrasound attenuation can be measured from the pubic symphysis, 328. As in other measurements, a ratio between broadband ultrasound attenuation velocities of these structures, using two or more broadband ultrasound attenuation measurements and/or averages of broadband ultrasound attenuation within each structure can be used to measure bone age. - In an exemplary embodiment,
probe 344 indicates and/or tracks a disease state such as slipped capitus femoris, wherefemoral head 312 is displaced onfemur 310. Additionally or alternatively,ultrasound probe 344 is positioned overpubic symphysis 328 to give acoustic velocity information. To determine SOS, two or more transducers should be contained in the probe or two or more separate transducers should be used. - FIG. 3B illustrates a cross-sectional diagram of a
pubic symphysis 328 along line III-III in a mature male and an associatedgraph 370 of acoustic backscatter attenuation.Pubic symphysis 328 is made up of aninterpubic disc 352 that is fibrocartilaginous hyaline cartilage. One either side ishyaline cartilage pubic tubercle Pubic tubercle 356 is ossified atage 14. - In
graph 370, as the acoustic signal measurespubic tubercle 356, the acoustic backscatter is high due to the ossification content of the bone, corresponding to aregion 372 ofgraph 370. Athyaline cartilage 350, the acoustic backscatter decreases as seen in aregion 374 ofgraph 370, due to the lack of ossification of this structure. Atinterpubic disc 352 the signal backscatter decreases even further, as seen in aregion 376 ofgraph 370, as this structure is fibrocartilaginous, with less dense fibers interspersed in the cartilage. - At
hyaline cartilage 350 a the signal backscatter increases representing the greater density of hyaline cartilage as shown in anarea 378 ingraph 370. Atpubic tubercle 356 a, the signal backscatter increases even further representing the greater density of bone as represented by anarea 380 ofgraph 370. - According to T. W. Todd, in “Age changes in the pubic bones,”Am. J. Phys. Anthropology, 3:285-334 and T. W. McKern and T. D. Stewart, “Skeletal age changes in young American males” Tech Rep. EP 45. Environmental Protection Research Div. Natick, Mass., changes within the pubic symphysis, such as changes in the bony surfaces, occur from the late teens to the fifties and beyond. Similarly, skeletal age is measured by the appearance and changes in the cranial sutures above 25 years of age. (T. W. Todd and D. W. Lyon, “Endocranial suture closure, its progress and age relationship. Am. J. Phys. Anthropol. 7:325-384). Thus, using acoustic signals in one or more structures that exhibit changes after bone growth has ceased in the long bones of the body, bone age can continue to be measured. Such measurements can also be used to indication the presence of a disease such as acromegaly, a diagnosis that can be confirmed by measuring the amount of somatotropin hormone production in the blood.
- While the invention has been described with respect to limited number of embodiments, it will be appreciated that many variations, modifications and other applications of the invention may be made. Also, combination of elements from variations may be combined and single elements may be used. Any and all such variations and modifications, as well as others that may become apparent to those skilled in the art are intended to be included within the scope of the invention, as defined by the appended claims.
- The terms “include”, “comprise” and “have” and their conjugates as used herein mean “including but not necessarily limited to.”
- It will be appreciated by a person skilled in the art that the present invention is not limited by what has thus far been described. Rather, the scope of the present invention is limited only by the following claims.
Claims (61)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IL2001/000683 WO2003009758A1 (en) | 2001-07-24 | 2001-07-24 | Bone age assessment using ultrasound |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IL2001/000683 Continuation WO2003009758A1 (en) | 2001-07-24 | 2001-07-24 | Bone age assessment using ultrasound |
Publications (2)
Publication Number | Publication Date |
---|---|
US20030065264A1 true US20030065264A1 (en) | 2003-04-03 |
US7678049B2 US7678049B2 (en) | 2010-03-16 |
Family
ID=11043073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/042,735 Expired - Fee Related US7678049B2 (en) | 2001-07-24 | 2001-10-25 | Bone age assessment using ultrasound |
Country Status (3)
Country | Link |
---|---|
US (1) | US7678049B2 (en) |
JP (1) | JP2004535882A (en) |
WO (1) | WO2003009758A1 (en) |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020162031A1 (en) * | 2001-03-08 | 2002-10-31 | Shmuel Levin | Method and apparatus for automatic control of access |
FR2889803A1 (en) * | 2005-08-22 | 2007-02-23 | Patrick Hechard | Child`s final size estimating method, involves determining size ratio based on epiphyseal ratio value with preset table, and calculating estimated final size by multiplying size of child from representation of skeleton part and size ratio |
US20070237365A1 (en) * | 2006-04-07 | 2007-10-11 | Monro Donald M | Biometric identification |
US20080139912A1 (en) * | 2006-12-07 | 2008-06-12 | Electronics And Telecommunications Research Institute | Method of acquiring phalange edge image and system therefor |
US20080161674A1 (en) * | 2006-12-29 | 2008-07-03 | Donald Martin Monro | Active in vivo spectroscopy |
CN100434039C (en) * | 2006-11-22 | 2008-11-19 | 山东建筑大学 | Measurement method of fetal irregularly curved bones |
US20090048514A1 (en) * | 2006-03-09 | 2009-02-19 | Slender Medical Ltd. | Device for ultrasound monitored tissue treatment |
US20100056898A1 (en) * | 2005-06-13 | 2010-03-04 | Mckenna David R | Objective Methods of Estimating Age of Animals and Carcasses |
US20100063394A1 (en) * | 2006-11-14 | 2010-03-11 | Kagoshima University | Osteoporosis progress stage measuring device and osteoporosis progress stage measuring method |
US20100274161A1 (en) * | 2007-10-15 | 2010-10-28 | Slender Medical, Ltd. | Implosion techniques for ultrasound |
US7828734B2 (en) | 2006-03-09 | 2010-11-09 | Slender Medical Ltd. | Device for ultrasound monitored tissue treatment |
US20110060247A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | Methods and apparatus for measuring bone lengths |
US20110058726A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | System and method for performing bone densitometer measurements |
US20110058725A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | Methods and apparatus for measuring body circumference |
US7983461B2 (en) | 2007-12-10 | 2011-07-19 | Electronics And Telecommunications Research Institute | Method and system for extracting distal radius metaphysis |
US20110178541A1 (en) * | 2008-09-12 | 2011-07-21 | Slender Medical, Ltd. | Virtual ultrasonic scissors |
EP2454997A1 (en) * | 2010-11-18 | 2012-05-23 | Furuno Electric Company Limited | Method, device and program for analyzing cartilage using ultrasonic wave |
US8300911B1 (en) | 2009-09-10 | 2012-10-30 | General Electric Company | Methods and apparatus for measuring visceral fat mass |
RU2472441C1 (en) * | 2011-10-24 | 2013-01-20 | Государственное бюджетное образовательное учреждение высшего профессионального образования "Смоленская государственная медицинская академия" Министерства здравоохранения и социального развития Российской Федерации | Method of early diagnostics of bone formation impairment in infants and young children |
US20150196275A1 (en) * | 2014-01-14 | 2015-07-16 | Cyberlogic, Inc. | Ultrasound Measurement Device and Method for Ultrasonic Measurement |
US9107798B2 (en) | 2006-03-09 | 2015-08-18 | Slender Medical Ltd. | Method and system for lipolysis and body contouring |
JP2016538043A (en) * | 2013-11-30 | 2016-12-08 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Automated detection of patient age by ECG |
WO2017100322A1 (en) * | 2015-12-08 | 2017-06-15 | Kneevoice, Inc. | Assessing joint condition using acoustic sensors |
WO2018057714A1 (en) * | 2016-09-21 | 2018-03-29 | The General Hospital Corporation | Systems, methods and media for automatically generating a bone age assessment from a radiograph |
CN109255786A (en) * | 2018-09-30 | 2019-01-22 | 杭州依图医疗技术有限公司 | A kind of method and device detecting the stone age |
CN110009605A (en) * | 2019-03-21 | 2019-07-12 | 浙江工业大学 | A method and system for bone age prediction based on deep learning |
CN112545493A (en) * | 2020-12-04 | 2021-03-26 | 石家庄喜高科技有限责任公司 | Height evaluation method and terminal equipment |
CN113100827A (en) * | 2021-04-10 | 2021-07-13 | 汕头市超声仪器研究所股份有限公司 | A kind of ultrasonic bone age detection method |
US11367181B2 (en) * | 2018-12-29 | 2022-06-21 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for ossification center detection and bone age assessment |
US11963790B2 (en) | 2020-11-19 | 2024-04-23 | Merative Us L.P. | Estimating spinal age |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006014979A (en) * | 2004-07-02 | 2006-01-19 | Oyo Denki Kk | Bone strength measuring device |
JP2007065892A (en) * | 2005-08-30 | 2007-03-15 | Ginga Kobo:Kk | Antiaging management system |
JP4690204B2 (en) * | 2006-01-16 | 2011-06-01 | 富士フイルム株式会社 | Image reproduction apparatus and program thereof |
FR2946753B1 (en) * | 2009-06-11 | 2011-07-22 | Centre Nat Rech Scient | ULTRASONIC METHOD AND DEVICE FOR CHARACTERIZING A MEDIUM |
KR20130137203A (en) * | 2011-01-06 | 2013-12-16 | 더루우브리졸코오포레이션 | Improvements in or relating to ultrasound generating apparatus, and methods for generating ultrasound |
AU2012332817B2 (en) | 2011-10-28 | 2017-03-30 | Decision Sciences International Corporation | Spread spectrum coded waveforms in ultrasound imaging |
JP6303905B2 (en) * | 2014-08-06 | 2018-04-04 | コニカミノルタ株式会社 | Ultrasonic diagnostic apparatus, image processing method of ultrasonic diagnostic apparatus, and image processing program |
US9844359B2 (en) | 2013-09-13 | 2017-12-19 | Decision Sciences Medical Company, LLC | Coherent spread-spectrum coded waveforms in synthetic aperture image formation |
JP6071101B1 (en) | 2014-01-17 | 2017-02-01 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Ultrasound apparatus and method for evaluating a subject's bone |
JP2016123733A (en) * | 2015-01-05 | 2016-07-11 | パイオニア株式会社 | Portable device |
CA2977975A1 (en) | 2015-02-25 | 2016-09-01 | Decision Sciences Medical Company, LLC | Acoustic signal transmission couplants and coupling mediums |
JP6751550B2 (en) * | 2015-05-15 | 2020-09-09 | 森永乳業株式会社 | Intestinal age estimation method and intestinal age estimation device |
AU2016334258B2 (en) | 2015-10-08 | 2021-07-01 | Decision Sciences Medical Company, LLC | Acoustic orthopedic tracking system and methods |
CN107515250B (en) * | 2017-09-04 | 2020-09-22 | 中国航空工业集团公司基础技术研究院 | Acoustic microscopic imaging device for tooth body |
CN108334899A (en) * | 2018-01-28 | 2018-07-27 | 浙江大学 | Quantify the bone age assessment method of information integration based on hand bone X-ray bone and joint |
CN108606811A (en) * | 2018-04-12 | 2018-10-02 | 上海交通大学医学院附属上海儿童医学中心 | A kind of ultrasound stone age detecting system and its method |
CN109036560B (en) * | 2018-07-09 | 2021-10-01 | 广西壮族自治区妇幼保健院 | Bone Development Analysis System |
JP6999812B2 (en) * | 2018-08-01 | 2022-01-19 | 中國醫藥大學附設醫院 | Bone age evaluation and height prediction model establishment method, its system and its prediction method |
KR102355327B1 (en) * | 2018-09-10 | 2022-02-21 | 주식회사 소노엠 | A system for measuring bone age |
WO2020181213A1 (en) | 2019-03-06 | 2020-09-10 | Decision Sciences Medical Company, LLC | Methods for manufacturing and distributing semi-rigid acoustic coupling articles and packaging for ultrasound imaging |
KR102128997B1 (en) * | 2019-03-25 | 2020-07-02 | 본와이즈 주식회사 | Apparatus, method and recording medium storing command for determining dental age |
WO2020219705A1 (en) | 2019-04-23 | 2020-10-29 | Allan Wegner | Semi-rigid acoustic coupling articles for ultrasound diagnostic and treatment applications |
KR102243969B1 (en) * | 2019-05-07 | 2021-04-22 | 서강대학교산학협력단 | Bladder volume measuring device and operating method thereof |
DE112020007252T8 (en) * | 2020-05-29 | 2023-05-04 | Mitsubishi Electric Corporation | Occupant detection apparatus, left alone child alarm system, method of detecting an occupant and method of issuing an alarm about a left alone child |
WO2021240777A1 (en) * | 2020-05-29 | 2021-12-02 | 三菱電機株式会社 | Occupant detection device and occupant detection method |
RU2747163C1 (en) * | 2020-07-21 | 2021-04-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Смоленский государственный медицинский университет" министерства здравоохранения Российской Федерации | A method for predicting the risk of a decrease in bone strength in newborns from mothers with excessive gestational weight gain who entered pregnancy with a normal body weight |
KR102436358B1 (en) * | 2020-07-29 | 2022-08-26 | 전남대학교산학협력단 | Age estimation method and system from PMCT scan image |
CN116685847A (en) | 2020-11-13 | 2023-09-01 | 决策科学医疗有限责任公司 | System and method for synthetic aperture ultrasound imaging of objects |
CN112447302B (en) * | 2020-12-04 | 2023-01-06 | 石家庄喜高科技有限责任公司 | Height growth speed evaluation method and terminal equipment |
CN114795288A (en) * | 2021-01-22 | 2022-07-29 | 复旦大学 | Human body osteoporosis detection method based on ultrasonic guided wave frequency dispersion curve |
KR102776035B1 (en) * | 2022-03-21 | 2025-03-06 | 주식회사 바디맵 | Apparatus and method for predicting bone growth based on finger joint ultrasound image |
WO2024162741A1 (en) * | 2023-01-31 | 2024-08-08 | 고려대학교 산학협력단 | Method and apparatus for assessing pubertal elbow bone age using lateral elbow x-ray images and artificial intelligence model |
Citations (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4490640A (en) * | 1983-09-22 | 1984-12-25 | Keisuke Honda | Multi-frequency ultrasonic transducer |
US4752918A (en) * | 1983-06-23 | 1988-06-21 | Etat Francais | Electrio-acoustic transducers |
US4774959A (en) * | 1986-01-10 | 1988-10-04 | Walker Sonix A/S | Narrow band ultrasonic frequency attentuation bone measurement system |
US4870972A (en) * | 1987-03-19 | 1989-10-03 | Thomson-Csf | Multiple-frequency acoustic transducer, especially for medical imaging |
US4913157A (en) * | 1986-06-03 | 1990-04-03 | Analog Devices, Inc. | Ultrasound method and apparatus for evaluating, in vivo, bone conditions |
US4926870A (en) * | 1988-08-30 | 1990-05-22 | Osteo-Technology, Inc. | Method and apparatus for ultrasonic analysis of bone strength in vivo |
US4963782A (en) * | 1988-10-03 | 1990-10-16 | Ausonics Pty. Ltd. | Multifrequency composite ultrasonic transducer system |
US5197475A (en) * | 1988-08-10 | 1993-03-30 | The Board Of Regents, The University Of Texas System | Method and apparatus for analyzing material properties using ultrasound |
US5452722A (en) * | 1992-06-22 | 1995-09-26 | Mccue Ultrasonics, Ltd. | Ultrasound bone analyser |
US5458130A (en) * | 1992-07-30 | 1995-10-17 | Orthologic Corporation | Ultrasonic therapy and assessment apparatus and method |
US5483965A (en) * | 1988-05-11 | 1996-01-16 | Lunar Corporation | Ultrasonic densitometer device and method |
US5533519A (en) * | 1993-06-24 | 1996-07-09 | Radke; John C. | Method and apparatus for diagnosing joints |
US5615681A (en) * | 1994-12-22 | 1997-04-01 | Aloka Co., Ltd | Method for measuring speed of sound in tissue and tissue assessment apparatus |
US5761156A (en) * | 1995-04-03 | 1998-06-02 | Marco Systemanalyse Und | Piezoelectric ultrasonic transducer |
US5806520A (en) * | 1994-03-25 | 1998-09-15 | Centre National De La Recherche Scientifique (C.N.R.S.) | Method and device for evaluating and characterizing the properties of bone |
US5824085A (en) * | 1996-09-30 | 1998-10-20 | Integrated Surgical Systems, Inc. | System and method for cavity generation for surgical planning and initial placement of a bone prosthesis |
US5895364A (en) * | 1996-06-20 | 1999-04-20 | The Trustees Of Stevens Institute Of Technology | Non-invasive technique for bone mass measurement |
US6086538A (en) * | 1995-11-22 | 2000-07-11 | Osteometer Meditech A/S | Methods and apparatus for evaluation of bone condition |
US6095979A (en) * | 1996-01-29 | 2000-08-01 | Aloka Co., Ltd. | Bone assessment apparatus |
US6244101B1 (en) * | 1999-06-01 | 2001-06-12 | Battelle Memorial Institute | Photoacoustic method for measuring concentration of chemical species |
US6344024B1 (en) * | 1997-05-12 | 2002-02-05 | Dwl Elektronische Systeme Gmbh | Multifrequency ultrasound probe |
US6413215B1 (en) * | 2000-04-05 | 2002-07-02 | The University Of Vermont | Implant wear debris detection apparatus and method |
US6463167B1 (en) * | 1996-09-19 | 2002-10-08 | Philips Medical Systems Technologies Ltd. | Adaptive filtering |
US20020162031A1 (en) * | 2001-03-08 | 2002-10-31 | Shmuel Levin | Method and apparatus for automatic control of access |
US20030035773A1 (en) * | 2001-07-27 | 2003-02-20 | Virtualscopics Llc | System and method for quantitative assessment of joint diseases and the change over time of joint diseases |
US6585666B2 (en) * | 1998-10-13 | 2003-07-01 | The Administrators Of The Tulane Educational Fund | Arthroscopic diagnostic probe to measure mechanical properties of articular cartilage |
US6607478B2 (en) * | 1999-03-02 | 2003-08-19 | Medtronic Ave, Inc. | Brachytherapy device assembly and method of use |
US20030229299A1 (en) * | 2002-06-05 | 2003-12-11 | Fuji Photo Film Co., Ltd. | Method of and system for evaluating joint and computer program for carrying out the method |
US6799066B2 (en) * | 2000-09-14 | 2004-09-28 | The Board Of Trustees Of The Leland Stanford Junior University | Technique for manipulating medical images |
US20040193048A1 (en) * | 2001-07-24 | 2004-09-30 | Liat Tsoref | Joint analysis using ultrasound |
US6836557B2 (en) * | 2002-10-02 | 2004-12-28 | VirtualS{tilde over (c)}opics, LLC | Method and system for assessment of biomarkers by measurement of response to stimulus |
Family Cites Families (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2827423C2 (en) * | 1978-06-22 | 1987-04-16 | Philips Patentverwaltung Gmbh, 2000 Hamburg | Device for determining the internal structure of a body using sound rays |
JPS5757538A (en) * | 1980-09-25 | 1982-04-06 | Tokyo Shibaura Electric Co | Ultrasonic imaging apparatus |
JPS6098799A (en) | 1983-11-02 | 1985-06-01 | Olympus Optical Co Ltd | Layer-built ultrasonic transducer |
US4973666A (en) | 1987-11-06 | 1990-11-27 | Washington Research Foundation | Peptide fragments containing HP and LP cross-links |
US4930511A (en) | 1988-05-11 | 1990-06-05 | Lunar Radiation, Inc. | Ultrasonic densitometer device and method |
US5343863A (en) * | 1988-05-11 | 1994-09-06 | Lunar Corporation | Ultrasonic densitometer device and method |
DE3925014A1 (en) | 1989-07-27 | 1991-01-31 | Joachim Dr Med Grifka | Testing stability of human knee joint - applying pressure to bendable holder for exact parallel sliding of shin bone w.r.t. thigh bone |
US5509042A (en) | 1991-02-13 | 1996-04-16 | Lunar Corporation | Automated determination and analysis of bone morphology |
US6160866A (en) | 1991-02-13 | 2000-12-12 | Lunar Corporation | Apparatus for bilateral femur measurement |
US5577089A (en) | 1991-02-13 | 1996-11-19 | Lunar Corporation | Device and method for analysis of bone morphology |
GB2257253B (en) | 1991-06-17 | 1995-01-11 | Christian Mcdonald Langton | Ultrasound bone analyser |
US5309898A (en) * | 1992-07-30 | 1994-05-10 | Kaufman Jonathan J | Ultrasonic bone-therapy and assessment apparatus and method |
GB2278436A (en) * | 1993-05-28 | 1994-11-30 | Kevin Hill | Image processing system and method for automatic feature extraction |
JPH0716224A (en) * | 1993-06-30 | 1995-01-20 | Shimadzu Corp | Ultrasonic transmission inspection device |
JP2664628B2 (en) * | 1993-11-08 | 1997-10-15 | アロカ株式会社 | Bone evaluation device for calcaneus |
US5439001A (en) | 1993-11-17 | 1995-08-08 | Ivac Corporation | Flexible diaphragm tonometer |
IT1268599B1 (en) | 1994-01-14 | 1997-03-06 | Igea Srl | ULTRASONIC MEASURING SYSTEM FOR DETECTION OF DENSITY AND BONE STRUCTURE. |
JP3207318B2 (en) | 1994-05-11 | 2001-09-10 | アロカ株式会社 | Ultrasonic bone evaluation device |
US5487388A (en) | 1994-11-01 | 1996-01-30 | Interspec. Inc. | Three dimensional ultrasonic scanning devices and techniques |
DE69634714T2 (en) | 1995-03-31 | 2006-01-19 | Kabushiki Kaisha Toshiba, Kawasaki | Therapeutic ultrasound device |
US5651363A (en) * | 1996-02-16 | 1997-07-29 | Orthologic Corporation | Ultrasonic bone assessment method and apparatus |
JPH09224935A (en) * | 1996-02-27 | 1997-09-02 | Furuno Electric Co Ltd | Ultrasonic bone measuring method and device therefor |
US5785041A (en) * | 1996-03-26 | 1998-07-28 | Hologic Inc. | System for assessing bone characteristics |
US5605155A (en) | 1996-03-29 | 1997-02-25 | University Of Washington | Ultrasound system for automatically measuring fetal head size |
US5704356A (en) | 1996-05-06 | 1998-01-06 | Raycont Ltd. | System and method of diagnosis of malformed hips in babies or small animals |
JP3119825B2 (en) | 1996-11-08 | 2000-12-25 | アロカ株式会社 | Ultrasonic bone evaluation device |
GB2321704A (en) | 1997-02-01 | 1998-08-05 | Huntleigh Technology Plc | Ultrasound device for measuring bone density |
FR2768322B1 (en) * | 1997-09-18 | 1999-12-10 | Patrick Hechard | SIZE FORECASTING METHOD BY DETERMINING THE GROWTH INDEX |
US6074352A (en) | 1998-03-26 | 2000-06-13 | Brigham And Women's Hospital | Method for the treatment of joint diseases characterized by unwanted pannus |
US6135960A (en) | 1998-08-31 | 2000-10-24 | Holmberg; Linda Jean | High-resolution, three-dimensional whole body ultrasound imaging system |
CN1335941A (en) | 1998-11-06 | 2002-02-13 | 康宁股份有限公司 | Athermal optical waveguide grating device |
WO2000028316A1 (en) | 1998-11-08 | 2000-05-18 | Sunlight Medical, Ltd. | Bone velocity determination |
EP1173097B1 (en) | 1999-04-16 | 2008-07-23 | Texon Technologies Ltd. | Ultrasonic scanning apparatus |
JP4702968B2 (en) * | 1999-08-25 | 2011-06-15 | 株式会社東芝 | Ultrasonic diagnostic equipment |
US6306089B1 (en) * | 1999-09-24 | 2001-10-23 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with customized measurements and calculations |
US6607487B2 (en) | 2001-01-23 | 2003-08-19 | The Regents Of The University Of California | Ultrasound image guided acetabular implant orientation during total hip replacement |
JP2005504563A (en) | 2001-05-25 | 2005-02-17 | イメージング セラピューティクス,インコーポレーテッド | Methods and compositions for resurfacing joints |
US6468215B1 (en) * | 2001-07-16 | 2002-10-22 | Artann Laboratories | Method and device for multi-parametric ultrasonic assessment of bone conditions |
-
2001
- 2001-07-24 JP JP2003515156A patent/JP2004535882A/en active Pending
- 2001-07-24 WO PCT/IL2001/000683 patent/WO2003009758A1/en active Application Filing
- 2001-10-25 US US10/042,735 patent/US7678049B2/en not_active Expired - Fee Related
Patent Citations (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4752918A (en) * | 1983-06-23 | 1988-06-21 | Etat Francais | Electrio-acoustic transducers |
US4490640A (en) * | 1983-09-22 | 1984-12-25 | Keisuke Honda | Multi-frequency ultrasonic transducer |
US4774959A (en) * | 1986-01-10 | 1988-10-04 | Walker Sonix A/S | Narrow band ultrasonic frequency attentuation bone measurement system |
US4913157A (en) * | 1986-06-03 | 1990-04-03 | Analog Devices, Inc. | Ultrasound method and apparatus for evaluating, in vivo, bone conditions |
US4870972A (en) * | 1987-03-19 | 1989-10-03 | Thomson-Csf | Multiple-frequency acoustic transducer, especially for medical imaging |
US5483965A (en) * | 1988-05-11 | 1996-01-16 | Lunar Corporation | Ultrasonic densitometer device and method |
US5197475A (en) * | 1988-08-10 | 1993-03-30 | The Board Of Regents, The University Of Texas System | Method and apparatus for analyzing material properties using ultrasound |
US4926870A (en) * | 1988-08-30 | 1990-05-22 | Osteo-Technology, Inc. | Method and apparatus for ultrasonic analysis of bone strength in vivo |
US4963782A (en) * | 1988-10-03 | 1990-10-16 | Ausonics Pty. Ltd. | Multifrequency composite ultrasonic transducer system |
US5452722A (en) * | 1992-06-22 | 1995-09-26 | Mccue Ultrasonics, Ltd. | Ultrasound bone analyser |
US5458130A (en) * | 1992-07-30 | 1995-10-17 | Orthologic Corporation | Ultrasonic therapy and assessment apparatus and method |
US5533519A (en) * | 1993-06-24 | 1996-07-09 | Radke; John C. | Method and apparatus for diagnosing joints |
US5806520A (en) * | 1994-03-25 | 1998-09-15 | Centre National De La Recherche Scientifique (C.N.R.S.) | Method and device for evaluating and characterizing the properties of bone |
US5615681A (en) * | 1994-12-22 | 1997-04-01 | Aloka Co., Ltd | Method for measuring speed of sound in tissue and tissue assessment apparatus |
US5761156A (en) * | 1995-04-03 | 1998-06-02 | Marco Systemanalyse Und | Piezoelectric ultrasonic transducer |
US6086538A (en) * | 1995-11-22 | 2000-07-11 | Osteometer Meditech A/S | Methods and apparatus for evaluation of bone condition |
US6095979A (en) * | 1996-01-29 | 2000-08-01 | Aloka Co., Ltd. | Bone assessment apparatus |
US5895364A (en) * | 1996-06-20 | 1999-04-20 | The Trustees Of Stevens Institute Of Technology | Non-invasive technique for bone mass measurement |
US6463167B1 (en) * | 1996-09-19 | 2002-10-08 | Philips Medical Systems Technologies Ltd. | Adaptive filtering |
US5824085A (en) * | 1996-09-30 | 1998-10-20 | Integrated Surgical Systems, Inc. | System and method for cavity generation for surgical planning and initial placement of a bone prosthesis |
US6344024B1 (en) * | 1997-05-12 | 2002-02-05 | Dwl Elektronische Systeme Gmbh | Multifrequency ultrasound probe |
US6585666B2 (en) * | 1998-10-13 | 2003-07-01 | The Administrators Of The Tulane Educational Fund | Arthroscopic diagnostic probe to measure mechanical properties of articular cartilage |
US6607478B2 (en) * | 1999-03-02 | 2003-08-19 | Medtronic Ave, Inc. | Brachytherapy device assembly and method of use |
US6244101B1 (en) * | 1999-06-01 | 2001-06-12 | Battelle Memorial Institute | Photoacoustic method for measuring concentration of chemical species |
US6413215B1 (en) * | 2000-04-05 | 2002-07-02 | The University Of Vermont | Implant wear debris detection apparatus and method |
US6799066B2 (en) * | 2000-09-14 | 2004-09-28 | The Board Of Trustees Of The Leland Stanford Junior University | Technique for manipulating medical images |
US20020162031A1 (en) * | 2001-03-08 | 2002-10-31 | Shmuel Levin | Method and apparatus for automatic control of access |
US20040193048A1 (en) * | 2001-07-24 | 2004-09-30 | Liat Tsoref | Joint analysis using ultrasound |
US7204807B2 (en) * | 2001-07-24 | 2007-04-17 | Sunlight Medical Ltd. | Joint analysis using ultrasound |
US20030035773A1 (en) * | 2001-07-27 | 2003-02-20 | Virtualscopics Llc | System and method for quantitative assessment of joint diseases and the change over time of joint diseases |
US20030229299A1 (en) * | 2002-06-05 | 2003-12-11 | Fuji Photo Film Co., Ltd. | Method of and system for evaluating joint and computer program for carrying out the method |
US6836557B2 (en) * | 2002-10-02 | 2004-12-28 | VirtualS{tilde over (c)}opics, LLC | Method and system for assessment of biomarkers by measurement of response to stimulus |
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040088553A1 (en) * | 2001-03-08 | 2004-05-06 | Shmuel Levin | Method and apparatus for automatic control of access |
US8267866B2 (en) | 2001-03-08 | 2012-09-18 | Shmuel Levin | Method and apparatus for assessment of bone age using ultrasound |
US7370208B2 (en) | 2001-03-08 | 2008-05-06 | Shmuel Levin | Method and apparatus for automatic control of access |
US20020162031A1 (en) * | 2001-03-08 | 2002-10-31 | Shmuel Levin | Method and apparatus for automatic control of access |
US20100056898A1 (en) * | 2005-06-13 | 2010-03-04 | Mckenna David R | Objective Methods of Estimating Age of Animals and Carcasses |
FR2889803A1 (en) * | 2005-08-22 | 2007-02-23 | Patrick Hechard | Child`s final size estimating method, involves determining size ratio based on epiphyseal ratio value with preset table, and calculating estimated final size by multiplying size of child from representation of skeleton part and size ratio |
US20090048514A1 (en) * | 2006-03-09 | 2009-02-19 | Slender Medical Ltd. | Device for ultrasound monitored tissue treatment |
US7828734B2 (en) | 2006-03-09 | 2010-11-09 | Slender Medical Ltd. | Device for ultrasound monitored tissue treatment |
US9107798B2 (en) | 2006-03-09 | 2015-08-18 | Slender Medical Ltd. | Method and system for lipolysis and body contouring |
US20070237365A1 (en) * | 2006-04-07 | 2007-10-11 | Monro Donald M | Biometric identification |
US20100063394A1 (en) * | 2006-11-14 | 2010-03-11 | Kagoshima University | Osteoporosis progress stage measuring device and osteoporosis progress stage measuring method |
CN100434039C (en) * | 2006-11-22 | 2008-11-19 | 山东建筑大学 | Measurement method of fetal irregularly curved bones |
US8050472B2 (en) | 2006-12-07 | 2011-11-01 | Electronics And Telecommunications Research Institute | Method of acquiring phalange edge image and system therefor |
US20080139912A1 (en) * | 2006-12-07 | 2008-06-12 | Electronics And Telecommunications Research Institute | Method of acquiring phalange edge image and system therefor |
US20080161674A1 (en) * | 2006-12-29 | 2008-07-03 | Donald Martin Monro | Active in vivo spectroscopy |
US20100274161A1 (en) * | 2007-10-15 | 2010-10-28 | Slender Medical, Ltd. | Implosion techniques for ultrasound |
US7983461B2 (en) | 2007-12-10 | 2011-07-19 | Electronics And Telecommunications Research Institute | Method and system for extracting distal radius metaphysis |
US20110178541A1 (en) * | 2008-09-12 | 2011-07-21 | Slender Medical, Ltd. | Virtual ultrasonic scissors |
US8280138B2 (en) * | 2009-09-10 | 2012-10-02 | General Electric Company | System and method for performing bone densitometer measurements |
US20110060247A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | Methods and apparatus for measuring bone lengths |
US8300911B1 (en) | 2009-09-10 | 2012-10-30 | General Electric Company | Methods and apparatus for measuring visceral fat mass |
US20110058725A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | Methods and apparatus for measuring body circumference |
US20110058726A1 (en) * | 2009-09-10 | 2011-03-10 | General Electric Company | System and method for performing bone densitometer measurements |
US8295570B2 (en) | 2009-09-10 | 2012-10-23 | General Electric Company | Methods and apparatus for measuring body circumference |
CN102525549A (en) * | 2010-11-18 | 2012-07-04 | 古野电气株式会社 | Method, device and program for analyzing cartilage using ultrasonic wave |
EP2454997A1 (en) * | 2010-11-18 | 2012-05-23 | Furuno Electric Company Limited | Method, device and program for analyzing cartilage using ultrasonic wave |
RU2472441C1 (en) * | 2011-10-24 | 2013-01-20 | Государственное бюджетное образовательное учреждение высшего профессионального образования "Смоленская государственная медицинская академия" Министерства здравоохранения и социального развития Российской Федерации | Method of early diagnostics of bone formation impairment in infants and young children |
US11744502B2 (en) | 2013-11-30 | 2023-09-05 | Koninklijke Philips N.V. | Automated detection of patient age by ECG |
JP2016538043A (en) * | 2013-11-30 | 2016-12-08 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Automated detection of patient age by ECG |
US20150196275A1 (en) * | 2014-01-14 | 2015-07-16 | Cyberlogic, Inc. | Ultrasound Measurement Device and Method for Ultrasonic Measurement |
WO2015108843A1 (en) * | 2014-01-14 | 2015-07-23 | Cyberlogic, Inc. | Ultrasound measurement device and method for ultrasonic measurement |
CN106102563A (en) * | 2014-01-14 | 2016-11-09 | 计算机逻辑公司 | Ultrasonic measurement devices and method for ultrasonic measurement |
WO2017100322A1 (en) * | 2015-12-08 | 2017-06-15 | Kneevoice, Inc. | Assessing joint condition using acoustic sensors |
US10448919B2 (en) | 2015-12-08 | 2019-10-22 | Kneevoice, Inc. | Assessing joint condition using acoustic sensors |
US11844643B2 (en) | 2015-12-08 | 2023-12-19 | Kneevoice, Inc. | Assessing joint condition using acoustic sensors |
US10991093B2 (en) | 2016-09-21 | 2021-04-27 | The General Hospital Corporation | Systems, methods and media for automatically generating a bone age assessment from a radiograph |
WO2018057714A1 (en) * | 2016-09-21 | 2018-03-29 | The General Hospital Corporation | Systems, methods and media for automatically generating a bone age assessment from a radiograph |
CN109255786A (en) * | 2018-09-30 | 2019-01-22 | 杭州依图医疗技术有限公司 | A kind of method and device detecting the stone age |
US11367181B2 (en) * | 2018-12-29 | 2022-06-21 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for ossification center detection and bone age assessment |
US20220309663A1 (en) * | 2018-12-29 | 2022-09-29 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for ossification center detection and bone age assessment |
US11735322B2 (en) * | 2018-12-29 | 2023-08-22 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for ossification center detection and bone age assessment |
CN110009605A (en) * | 2019-03-21 | 2019-07-12 | 浙江工业大学 | A method and system for bone age prediction based on deep learning |
US11963790B2 (en) | 2020-11-19 | 2024-04-23 | Merative Us L.P. | Estimating spinal age |
CN112545493A (en) * | 2020-12-04 | 2021-03-26 | 石家庄喜高科技有限责任公司 | Height evaluation method and terminal equipment |
CN113100827A (en) * | 2021-04-10 | 2021-07-13 | 汕头市超声仪器研究所股份有限公司 | A kind of ultrasonic bone age detection method |
Also Published As
Publication number | Publication date |
---|---|
US7678049B2 (en) | 2010-03-16 |
WO2003009758A1 (en) | 2003-02-06 |
JP2004535882A (en) | 2004-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7678049B2 (en) | Bone age assessment using ultrasound | |
Laugier et al. | Ultrasound parametric imaging of the calcaneus: in vivo results with a new device | |
US6899680B2 (en) | Ultrasound measurement techniques for bone analysis | |
US7204807B2 (en) | Joint analysis using ultrasound | |
US8419643B2 (en) | Ultrasonic method and apparatus for assessment of bone | |
Blake et al. | Applications of bone densitometry for osteoporosis | |
Moilanen et al. | Assessment of the tibia using ultrasonic guided waves in pubertal girls | |
WO1990001903A1 (en) | Method and apparatus for ultrasonic bone analysis | |
Barkmann et al. | Femur ultrasound (FemUS)—first clinical results on hip fracture discrimination and estimation of femoral BMD | |
Castriota-Scanderbeg et al. | Precision of sonographic measurement of articular cartilage: inter-and intraobserver analysis | |
CN101199429B (en) | Ultrasonic bone assessment apparatus and method | |
Laugier | An overview of bone sonometry | |
Damilakis et al. | Ultrasound velocity through the cortex of phalanges, radius, and tibia in normal and osteoporotic postmenopausal women using a new multisite quantitative ultrasound device | |
US20070016061A1 (en) | System for measuring and tracking human body fat | |
CN102186422B (en) | Ultrasonic bone assessment apparatus and method | |
Wan et al. | Summation of ossification ratios of radius, ulna and femur: a new parameter to evaluate bone age by ultrasound | |
RU2598642C2 (en) | Method and device for estimating mineral density of a bone | |
JP2596701B2 (en) | Bone evaluation device | |
McLean et al. | Elderly cohort study subjects unable to return for follow-up have lower bone mass than those who can return | |
Padilla et al. | In vitro ultrasound measurement at the human femur | |
Laugier | Quantitative ultrasound instrumentation for bone in vivo characterization | |
Bi et al. | Ultrasonic backscatter measurements of human cortical and trabecular bone densities in a head-down bed-rest study | |
Bozkurt et al. | The diagnostic accuracy of ultrasonography in determining the reduction success of distal radius fractures | |
Shimura et al. | Assessment of measurement of children's bone age ultrasonically with Sunlight BonAge | |
Määttä et al. | Correlation of tibial low-frequency ultrasound velocity with femoral radiographic measurements and BMD in elderly women |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SUNLIGHT MEDICAL LTD., ISRAEL Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TSOREF, LIAT;FURMAN, DOV;REEL/FRAME:013898/0195 Effective date: 20011024 Owner name: SUNLIGHT MEDICAL LTD.,ISRAEL Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TSOREF, LIAT;FURMAN, DOV;REEL/FRAME:013898/0195 Effective date: 20011024 |
|
AS | Assignment |
Owner name: BEAM-MED LTD., ISRAEL Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SUNLIGHT MEDICAL LTD.;REEL/FRAME:020419/0858 Effective date: 20070621 Owner name: BEAM-MED LTD.,ISRAEL Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SUNLIGHT MEDICAL LTD.;REEL/FRAME:020419/0858 Effective date: 20070621 |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.) |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.) |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Expired due to failure to pay maintenance fee |
Effective date: 20180316 |