CN104080404B - Bone density measurement device - Google Patents
Bone density measurement device Download PDFInfo
- Publication number
- CN104080404B CN104080404B CN201280068207.1A CN201280068207A CN104080404B CN 104080404 B CN104080404 B CN 104080404B CN 201280068207 A CN201280068207 A CN 201280068207A CN 104080404 B CN104080404 B CN 104080404B
- Authority
- CN
- China
- Prior art keywords
- bone density
- dimensional
- dimension
- image
- labelling
- 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.)
- Expired - Fee Related
Links
- 230000037182 bone density Effects 0.000 title claims abstract description 160
- 238000001739 density measurement Methods 0.000 title claims abstract description 24
- 238000002372 labelling Methods 0.000 claims abstract description 68
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 210000000988 bone and bone Anatomy 0.000 claims description 13
- 239000000284 extract Substances 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 239000003550 marker Substances 0.000 abstract 1
- 238000000034 method Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 8
- 238000005259 measurement Methods 0.000 description 6
- 239000002184 metal Substances 0.000 description 6
- 238000013211 curve analysis Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000002308 calcification Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 206010010214 Compression fracture Diseases 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 208000036626 Mental retardation Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/505—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Human Computer Interaction (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention relates to a kind of bone density measurement device, in bone density measurement device, when setting lines on bone density image, distribution curve corresponding thereto is shown as the distribution of one-dimensional bone density.By the analysis to distribution curve, and the position of characteristic element (characteristic point or characteristic interval) is specific.This position is shown as one-dimensional position labelling on distribution curve, is shown as two-dimensional position labelling on bone density image.When making above-mentioned lines move, distribution curve is real time updated, and the position of the most one-dimensional labelling and two-dimensional marker is real time updated.
Description
Technical area
The present invention relates to a kind of bone density measurement device, particularly relate to one and together show two dimension bone density
The bone density measurement device that image is distributed with one-dimensional bone density.
Background technology
Bone density measurement device is, based on by the predetermined measurement object being in measured body is irradiated X
Ray and obtain through data, and the device measuring the bone density measuring object is (with reference to patent
Document 1).Bone density measurement device is also sometimes referred to as bone mineral measuring device.
When bone density measurement, generally, high-energy X-rays and mental retardation are irradiated to the 2 dimensional region of measured body
Amount X-ray, and according to the X-ray transparent data of two kinds obtained by irradiating these X-ray
Two dimensional Distribution (the Dual Energy X-ray Absorptiometry method: dual intensity X penetrates of computing bone density
Line absorption method).This Two dimensional Distribution is shown as two dimension bone density image.Existing bone density measurement device has
Standby following function, i.e. show divide corresponding with the horizontal line specified by user on bone density image
The function of cloth curve (one-dimensional bone density distribution).Specifically, on display picture, in bone density figure
In the white space that the surrounding of picture is existing, distribution curve is shown as smaller image.Distribution curve
Transverse axis be the trunnion axis (horizontal line) in the axle representing distance, the scale of this axle and bone density image
Scale is the most consistent.
Citation
Patent documentation
Patent documentation 1: Japanese Unexamined Patent Publication 2011-167388 publication
Summary of the invention
Invent problem to be solved
In above-mentioned existing display packing, it is difficult to identify between bone density image and distribution curve is corresponding
Relation.Such as, it is impossible to promptly identify that the specific peak value on distribution curve is corresponding to bone density picture
Which position, additionally, cannot promptly identify that specifically corresponding partly on bone density image is distributed song
Which part on line.
It is an object of the invention to, in the situation that two dimension bone density image is shown with the distribution of one-dimensional bone density
Under diagnostic imaging is assisted.Or, it is an object of the invention to, it is possible to quickly and easily grasp
Corresponding relation between two dimension bone density image and the distribution of one-dimensional bone density.Or, the purpose of the present invention exists
In, using the teaching of the invention it is possible to provide the specific of the analysis subject area on two dimension bone density image or correction are carried out auxiliary
The information helped.
For the method solving problem
Bone density measurement device involved in the present invention is characterised by, comprises: bone density image forming part,
It forms two dimension bone density image according to by the detection data obtaining measured's X-ray irradiation；
Line of reference configuration part, it is for setting line of reference on described two dimension bone density image；Extraction unit, its from
Described two dimension bone density image extracts the one-dimensional bone density distribution corresponding with described line of reference；Analysis portion,
The distribution of described one-dimensional bone density is analyzed by it, and using characteristic point or characteristic interval as characteristic element
And carry out specific；Display part, its display described two dimension bone density image and the distribution of described one-dimensional bone density,
And on described two dimension bone density image, show the two-dimensional position of the two-dimensional position representing described characteristic element
Labelling, and in the distribution of described one-dimensional bone density, show one-dimensional position one-dimensional representing described characteristic element
Position mark.
According to said structure, when setting line of reference on two dimension bone density image, corresponding thereto
Dimension bone density distribution is extracted from two dimension bone density image.Display part shows two dimension bone density image
It is distributed with one-dimensional bone density.The result of analysis is performed, one-dimensional bone as the distribution of one-dimensional bone density
Characteristic element in Density Distribution is specific.Characteristic element is characterized a little or characteristic interval.Characteristic point is
Represent the point of such as maximum, minima, maximum value minimum, boundary position or position of centre of gravity.Special
Levying interval is, specific mountain part, specific valley are graded.Can also be distributed according to one-dimensional bone density
Analyze, carry out specific multiple characteristic elements.On two dimension bone density image, display represents the two of characteristic element
Dimension position mark, and on one-dimensional bone density is distributed, display represents the one of the one-dimensional position of this feature key element
Dimension position mark.Set two images by correspondence and observe such that it is able to be rapid and accurate
Ground identify one-dimensional bone density distribution on characteristic element corresponding to two dimension bone density image on which position,
Characteristic element on two dimension bone density image is corresponding to which position in the distribution of one-dimensional bone density on the contrary.
Thereby, it is possible to easily whether judging characteristic key element is equivalent to the object position of Analysis of bone mineral density, such as,
The setting of analysis object range relative to two dimension bone density image or region-of-interest can be implemented exactly
Setting.
Being preferably, multiple characteristic elements are carried out specific in described one-dimensional bone density is distributed by described analysis portion,
Described display part shows the two-dimensional position representing multiple described characteristic elements on described two dimension bone density image
Multiple two-dimensional position labellings, and described one-dimensional bone density distribution on show expression multiple described features will
Multiple one-dimensional position labellings of the one-dimensional position of element.It is preferably, for each characteristic element, makes Two-dimensional Position
Tagging and one-dimensional position are labeled as one group and distribute specific tone etc., it is possible to by respective group from other
Group difference.That is, it is preferably configured as, for each characteristic element, and two labellings can be identified moment
Between corresponding relation.In order to be identified in units of group, it is possible to so that shape beyond tone,
The not phase sample such as the kind of lines, brightness.Can also be identified by the mark of word or labelling.
It is preferably, on described two dimension bone density image during mobile described line of reference, on described line of reference
The position of described two-dimensional position labelling is real time updated, and, it is described that described one-dimensional bone density is distributed
The position of one-dimensional position labelling is real time updated.It is preferably, mobile institute on described two dimension bone density image
When stating line of reference, show the motion track of described two-dimensional position labelling.Can be according to the shape of motion track
Identify axis of bone, boundary edge etc..
Being preferably, described characteristic interval (i.e. characteristic area) is the part with one-dimensional extension, described
Two-dimensional position is labeled as the labelling of position and the size representing described characteristic interval, described one-dimensional position labelling
For representing position and the labelling of size of described characteristic interval.In that case, two-dimensional position labelling phase
When in first region, below labelling, one-dimensional position labelling is equivalent to following second area labelling.
Bone density measurement device involved in the present invention is characterised by, comprises；Bone density image forming part, its
Two dimension bone density image is formed according to the detection data by measured's X-ray irradiation is obtained；Ginseng
According to line configuration part, it is for setting line of reference on described two dimension bone density image；Extraction unit, it is from institute
State and two dimension bone density image extracts the one-dimensional bone density distribution corresponding with described line of reference；Analysis portion, its
The distribution of described one-dimensional bone density is analyzed, and characteristic interval is carried out specific using as characteristic element；
Display part, its described bone density image of display and the distribution of described one-dimensional bone density, and close at described bone
Show the first area labelling representing described characteristic interval on degree image, and be distributed in described one-dimensional bone density
Upper display represents the second area labelling of described characteristic interval.Can be at two dimension bone density figure according to this structure
As upper promptly specific interval is carried out specific and one-dimensional bone density distribution on the most promptly to specific interval
Carry out specific.And both corresponding relations also are able to easily identify.First area labelling can also be
Representing interval line portion, second area labelling can also be for all filling part.Characteristic interval is such as
Be equivalent to metal part, calcification part, compression fracture part etc..Such as, by bone density is distributed
Carry out threshold process thus extract characteristic.
It is preferably, on described two dimension bone density image during mobile described line of reference, on described line of reference
The position of described two-dimensional position labelling is real time updated, and, it is described that described one-dimensional bone density is distributed
The position of one-dimensional position labelling is real time updated.It is preferably, mobile institute on described two dimension bone density image
When stating line of reference, the two dimensional image of the motion track being equivalent to described two-dimensional position labelling is shown.
Structure accordingly, for example, it is possible to only give the metal colored parts of high brightness.
Accompanying drawing explanation
Fig. 1 is the block diagram of the preferred implementation representing bone density measurement device involved in the present invention.
Fig. 2 is the figure representing the first display example.
Fig. 3 is the figure representing the second display example.
Fig. 4 is the figure representing the 3rd display example.
Fig. 5 is the figure of the display representing rendered image.
Fig. 6 is the figure for illustrating motion track.
Fig. 7 is the figure for arbitrary line of reference is set for explanation.
Fig. 8 is the flow chart of the action example representing the device shown in Fig. 1.
Detailed description of the invention
Hereinafter, preferred embodiment illustrating the present invention with reference to the accompanying drawings.
In fig. 1 it is illustrated that the preferred implementation of bone density measurement device involved in the present invention, figure
1 for representing the block diagram of its total.This bone density measurement device is, is arranged in medical institutions,
And the bone density device that i.e. bone mineral amount measures to measured.
In FIG, bone density measurement device generally has measurement portion 10 and operational part 12.First,
Measurement portion 10 is illustrated.Measurement portion 10 has X-ray generation unit 18.This X-ray produces single
Unit 18 has the X-ray tube 20 producing X-ray.For this X-ray generation unit 18 by high voltage source
(HV) 21 supply electric power.In the present embodiment, it is formed with tool by X-ray generation unit 18
There is the X-ray beam 23 of shape i.e. fan beam (fan-beam) shape that end dissipates.Can also be formed
For pencil beam (pencil beam) or cone beam (cone beam).
Platform 14 is the mounting measured i.e. device of measured body 16.X-ray beam 23 is through measured 16
And arrive X-ray detector unit 22.In the present embodiment, X-ray detector unit 22 is by one-dimensional row
Row multiple X ray sensors and constitute.In the present embodiment, X-ray generation unit 18 and X
Ray detector units 22, is carried on body direction of principal axis i.e. paper vertical direction by conveying mechanism.By
This forms two dimension irradiation area.In order to bone density is carried out computing, i.e. in order to by skeleton portion from soft tissue
Part picks out, in the present embodiment, alternately irradiates high-energy X-rays and low-energy X-ray.
That is, each measurement position all obtains two detection data.
Bone density image forming part 24 is, forms two dimension bone density figure based on two this detection data
The module of picture.Although being formed with bone density image as the image representing skeleton portion in the present embodiment,
But soft tissue can also be reflected on this image.The view data of the bone density image formed via
Display process portion 26 and be sent to display part 28, thus show that bone is close on the display picture of display part 28
Degree image.
In the present embodiment, it is possible to utilize input unit 30, and set on bone density image by examiner
Line of reference.Thus, the distribution curve as the one-dimensional bone density distribution corresponding with this line of reference shows automatically
On picture.About this content described below.
Distribution curve preparing department 32 extracts from two dimension bone density image and is set for this image
The one-dimensional bone density distribution (the data row formed by multiple bone density data) that line of reference is corresponding.This is one-dimensional
Bone density distribution constitutes distribution curve.Represent that the data of this distribution curve are sent via display process portion 26
To display part 28.Additionally, this distribution curve is sent to distribution curve analysis portion 34.Distribution curve analysis portion
34 are, according to the condition specified by examiner, and implement for characteristic element (characteristic point or feature
Interval) parts specifically analyzed.In the case of this kind, it is also possible to specific multiple characteristic elements.
Here, as above-mentioned characteristic point, it is possible to list the point of expression peak, represent minimum
Point, represent maximum point, represent minimizing point, boundary position etc., as feature in addition
Point, it is possible to list position of centre of gravity etc..It is interval as features described above, it is possible to list specific mountain part,
Specific valley is graded.Waveform analysis is carried out such that it is able to by right by being distributed in one-dimensional bone density
Such as peak point carries out specific next specific peak.Same method can be utilized and automatically to each feature
Click on row operation.Such as, about boundary point, it is possible to by the rising point of waveform and drop point are carried out spy
Fixed, thus automatically identify this boundary point.About maximum and minimum, it is possible to according to carrying out waveform
The change of the differential value that differential processes and obtains, and this maximum and minimum are carried out specific.Center of gravity
Point is asked for by carrying out center of gravity computing for the part more than threshold value in waveform.Characteristic interval can lead to
Cross such as threshold process and carry out specific.
The analysis result of distribution curve analysis portion 34 is output to the first graph image forming portion 36 and
Two graph image forming portions 38.First graph image forming portion 36 is, is formed and is shown by overlap and i.e. to overlap
The module of the display the first graph image on two dimension bone density image.First graph image comprises expression ginseng
According to the figure of line, the figure of the one or more labellings represented on line of reference and other figure.Second figure
Shape image is, forms display second in the distribution of one-dimensional bone density i.e. wheel Guo that shown by overlap or overlap
The module of graph image.Second graph image comprise as one or more labellings one or more vertically
Line.
Display process portion 26 has image complex functionality, and to two dimension bone density image and the first graph image
Synthesize, and by defeated to display part 28 for the view data of the first composograph by being synthesized into
Go out.Additionally, distribution curve is synthesized by display process portion 26 with second graph image, and will be by closing
The view data of the second composograph become and obtain exports to display part 28.Subsidiary proposition, with each figure
The function that block represents, can be realized by software function.
Illustrate the first display example in fig. 2.In fig. 2, display picture 40 is set with first
Viewing area 42 and the second viewing area 48.Bone density image 44 is shown in the first viewing area 42,
And to show in the way of overlapping with this bone density image 44 as with reference to image or assistant images
First graph image 46.Bone density image 44 is for showing the bone density value of each position as brightness
The image gone out.Although the most only showing skeleton portion, but the part being equivalent to soft tissue is worked as
So can also be by image conversion.First graph image 46 has the lines 47 that distribution curve makes, and also has
There is multiple labelling 54,56,58,60.Here, labelling 54 is the position representing the peak on lines 47
The labelling put, labelling 56 is the labelling of the position representing the minimum on lines 47, labelling 58,60
For representing the labelling of two ends that is two boundary point in the skeleton portion on lines 47.Subsidiary proposition, symbol 62
Illustrate the region-of-interest (ROI) delimiting the region that mean bone density value is carried out computing.At bone density figure
In 44, transverse axis is X-direction, and the longitudinal axis is Y-direction.
In the second viewing area 48, distribution curve 50 is shown as waveform.Additionally, to overlap with it
Mode shows second graph image 52.The transverse axis of distribution curve 50 is equivalent to X-direction, distribution curve
The longitudinal axis of 50 corresponds to bone density value D.The two ends correspondence of the transverse axis on distribution curve 50 is set in lines
The two ends of 47, distribution curve 50 illustrates the one-dimensional bone density distribution on lines 47.Sometimes can also be by
Referred to herein as scatter chart.
As it has been described above, when being specified the position of lines 47 on bone density image 44 by examiner, from bone
Density image 44 extracts the bone density value data row being equivalent to lines 47, distribution curve 50 conduct accordingly
Waveform and be formed.And then, by performing the waveform analysis relative to this distribution curve 50 such that it is able to
The position of one or more characteristic point the most specific.In the example shown in this Fig. 2, had by specific
High level, minima, two skeleton portion boundary points.Labelling 54A represents the position of peak, labelling 56A
Representing the position of minimum, labelling 58A, 60A represent two boundary points.
Subsidiary proposition, in the illustrated example shown in fig. 2, forms each mark of the key element of the first graph image 46
Note has a little or the form of circle, on the other hand, forms the labelling of the key element of second graph image 52
54A, 56A, 58A, 60A are configured as vertical line.It is assigned different for each characteristic point
Tone, such as labelling 54 are shown by identical tone each other with labelling 54A.Similarly, labelling 56
It is shown by identical tone with labelling 56A, labelling 58 and labelling 58A and labelling 60 and labelling
60A is also shown by identical tone in the way of respective one group.Thus, user can identify moment
Between two graph images 46,52, the group of labelling is to relation, i.e. can be promptly to about each feature
Point position on respective image carries out specific.
When being changed the lines 47 position in vertical direction i.e. Y direction by examiner, second graph
The content of image 52 is real time updated, i.e. the generation of distribution curve is carried out in real time with analyzing.Meanwhile,
The content of the first graph image 46 is also real time updated.By such process, thus user will appreciate that
About the detailed content of the various piece of bone density image 44, therefore, it is possible to implement region-of-interest exactly
The setting of 62.Additionally, take from analyzing rejecting metal part object, from analyzing rejecting compressing object
Property fracture portions or, object is rejected the various measures such as calcification part will become easy from analyzing.
Additionally, apply the additional part etc. not being automatically recognized as skeleton portion owing to bone content is low to skeleton portion
Measure also will become easy.
Although in the above-described embodiment, utilize tone to implement the identification of each characteristic point, but also may be used
To utilize other display mode of the change of lines kind, carry out the identification of each characteristic point.
Illustrate the second display example in figure 3.It addition, for identical with the structure shown in Fig. 1
Also the description thereof will be omitted for structure tag same symbol.Other is shown that example is also same.
In the first viewing area 42, show bone density image 44A and the first graph image 46A.
On the other hand, the second viewing area 48 shows distribution curve 50 and second graph image 52A.
Distribution curve 50A represents the one-dimensional bone density distribution on lines 47, by being analyzed it and specific
Multiple very big positions and multiple minimum position, each position is expressed as labelling 64A～76A.?
This, labelling 70A, 72A, 74A, 76A are the labelling representing very big position, labelling 64A, 66A, 68A
For representing the labelling of each minimum position.
Corresponding to these first graph images 46A have on lines 47 shown multiple labellings 64～
76.Here, labelling 70,72,74,76 is the labelling representing very big position, labelling 64,66,68
For representing the labelling of minimum position respectively.Can be according to the arrangement of multiple labellings, interval, closeness, week
Phase property, identifies the circumstance or state that the bone density on lines 47 is distributed intuitively.With above-mentioned first
Show example similarly, even if in this second display example, when the position being changed lines 47 by examiner
Time, the content of second graph image 52A and the first graph image 46A is real time updated corresponding to which.
Illustrate the 3rd display example in the diagram.Bone density image is shown in the first viewing area 42
44B, also shows the first graph image 46B.Distribution curve 50B is shown in the second viewing area 48,
Also show second graph image 52B.In the 3rd display example, in distribution curve 50B, only
Automatically extracting the part beyond predetermined threshold value, this part 82 is shown by colour as coloured part.
Such as, such part is equivalent to metal part.
Corresponding to this, in the first graph image 46B, lines 47 are equivalent to above-mentioned part 82
Interval 80 be shown as coloring interval.The Show Color of part 82 is identical with the Show Color of interval 80,
Both corresponding relations are shown in the way of having same hue.Subsidiary proposing, symbol 78 represents and is buried
Enter the metalwork in skeleton.In the case of containing such metalwork 78, only this part is observed
Extra high bone density value, and such part in distribution curve 50B as the part all filled
Part 82 and specific, by specifically for interval 80 in this external bone density image 44B.
So, user is by observing such display such that it is able to take rapidly from average bone close
The object of degree computing is removed the measure of this part etc..
Illustrate the 4th display example in Figure 5.In this display example, on bone density image 44C
Lines 47 are scanned with the form shown in symbol 84, former as after image between colour attaching area at this moment
Sample retains, and its result is, produces extension rendered image 86 two-dimensionally.That is, the part of metalwork is equivalent to
Shown by colour as the image all filled.By in this way using interval motion track as
Color part and former state retain such that it is able to easily generate rendered image 86.
Illustrate the 5th display example in figure 6.In the 5th display example, picture shows
Specific characteristic point on lines 47 represents the labelling 88 of the position of peak the most in this example, with this
The motion track 90 of this labelling 88 is shown as curve simultaneously.That is, when with the side shown in symbol 84
When formula makes lines 47 move from the bottom to top, the extreme higher position on lines 47 is by real-time monitored, and its result is,
The motion track 90 that the position with peak can be made to be connected shows.According to this motion track 90
And be capable of identify that axis of bone, or can understand that along axis of bone how part that bone density value is higher intuitively
Change.The motion track of the labelling representing minimum can certainly be shown, if making expression skeleton in addition
Two motion tracks at the two ends in portion show, can show the profile in skeleton portion clearly.Can also show
Center of gravity line.
Illustrate the 6th display example in the figure 7.In this display example, on bone density image 44E,
Taken arbitrary direction and arbitrary length by examiner and be set with lines 92.One point of lines 92 is with A
Representing, another point represents with B.Corresponding and generate have distribution curve 50, on this transverse axis
End is corresponding to A, and its other end is corresponding to B.Its result is, specifies arbitrarily on two dimension bone density image
Path, and can make on this path one-dimensional bone density distribution show.Although lines are at Fig. 7
Shown example is made up of straight line, but these lines can also be made up of curve.
The action example of the device shown in Fig. 1 with flow chart illustration in fig. 8.
In S10, measured is placed on platform, and performs the measurement relative to measured.S12
In on picture, show bone density image.In S14, examiner on bone density image, set conduct
Mean bone density etc. is carried out the region-of-interest of the object range of computing.Such as can also be for multiple lumbar vertebras
In each lumbar vertebra and each set region-of-interest.The setting automatization of these region-of-interests can also be made.
In S16, carry out and add in region-of-interest or delete the operand becoming mean bone density value
Part.Such as, the rejecting of suitable metal part, the rejecting of compression fracture part, the picking of calcification part
Remove or the various operations such as addition of Low BMD part.S14 Yu S16 may be reversed order and carry out
Setting and when the adding or delete of analysis object of region-of-interest as above, when being judged to
Break as being highly difficult so needing to carry out S18 owing to only implementing various judgements according to bone density image
Investigation in the case of, in the case of being i.e. judged as needing the display of distribution curve and analyzing, perform S20
Each later operation.
In S20, bone density image is specified lines by examiner.In S22, from bone density figure
The bone density data row corresponding with lines are extracted in Xiang, and based on these bone density data row as waveform
Distribution curve is produced.This distribution curve is displayed on image.In S24, distribution curve is carried out
Analyze.That is, according to condition selected by the user, and one or more characteristic point is specific.And,
In S26, generate the first graph image and the second graph image representing one or more characteristic points,
These graph images are together displayed on picture with bone density image and distribution curve.In S28,
Judge whether the operation proceeded as above.When returning to master operation, perform to start from S14 each
Operation.That is, by referring to above distribution curve analysis result, and if needing to revise to so far
The region-of-interest being set, then perform so to revise operation, similarly, if needing to chase after in S14
Add or eliminate analysis object, then in S16, perform so operation.Can also be by S14 and S16
Operation be programmed in the operation of S20～S28.Example shown in Fig. 8 is only an example.Also
Can substitute for characteristic point and characteristic interval is specific.
At the end of being judged as investigation in S18, carry out according to the predetermined input carried out by user
S30, carries out the analysis of bone density image in S30.Specifically, known with reference in being region-of-interest
Not by being carried out the bone density value being distributed in the region part in skeleton portion, and arithmetic average bone based on this
Density.In S32, this analysis result is displayed on picture by numerical value etc..
According to above action example, region-of-interest setting, add or delete analyze object time, because of
On two dimensional image and one dimensional image, the analysis result of distribution curve can be identified, therefore obtain energy
Enough advantages quickly and accurately implementing the operations such as above-mentioned setting operation.At this moment, if application is such as schemed
Method shown in 4 grades, then owing to metal part can be identified as all filling chrominance section, therefore can
Easily this part is rejected.Furthermore, it is also possible to be easily implemented with the separation etc. of each vertebra.According to upper
State embodiment, because of can intuitively and easily identify position between bone density image and distribution curve and close
System, therefore can obtain the advantage implementing the auxiliary to diagnostic imaging etc.
Claims (8)
1. a bone density measurement device, it is characterised in that including:
Bone density image forming part, it is according to by the detection data obtaining measured's X-ray irradiation
Form two dimension bone density image；
Line of reference configuration part, it is for setting line of reference on described two dimension bone density image；
Extraction unit, it is close that it extracts the one-dimensional bone corresponding with described line of reference from described two dimension bone density image
Degree distribution；
Analysis portion, the distribution of described one-dimensional bone density is analyzed, and by characteristic point or characteristic interval by it
Carry out specific as characteristic element；
Display part, its display described two dimension bone density image and the distribution of described one-dimensional bone density, and
The two-dimensional position labelling of the two-dimensional position representing described characteristic element is shown on described two dimension bone density image,
And in the distribution of described one-dimensional bone density, show the one-dimensional position mark of the one-dimensional position representing described characteristic element
Note.
2. bone density measurement device as claimed in claim 1, it is characterised in that
Multiple characteristic elements are carried out specific in described one-dimensional bone density is distributed by described analysis portion,
Described display part shows the two dimension representing multiple described characteristic elements on described two dimension bone density image
Multiple two-dimensional position labellings of position, and show the multiple described spies of expression in the distribution of described one-dimensional bone density
Levy multiple one-dimensional position labellings of the one-dimensional position of key element.
3. bone density measurement device as claimed in claim 1, it is characterised in that
Described characteristic element is described characteristic point,
Described Based on Feature Points maximum, minima, maximum value minimum, boundary position or center of gravity position
Put.
4. bone density measurement device as claimed in claim 1, it is characterised in that
On described two dimension bone density image during mobile described line of reference, the described two dimension on described line of reference
The position of position mark is real time updated, and, the described one-dimensional position in the distribution of described one-dimensional bone density
The position of labelling is real time updated.
5. bone density measurement device as claimed in claim 4, it is characterised in that
On described two dimension bone density image during mobile described line of reference, the movement of described two-dimensional position labelling
Track is shown.
6. bone density measurement device as claimed in claim 1, it is characterised in that
Described characteristic element is described characteristic interval,
Described characteristic interval is distributed the specific mountain part or valley being had by described one-dimensional bone density
Point.
7. bone density measurement device as claimed in claim 6, it is characterised in that
On described two dimension bone density image during mobile described line of reference, the described two dimension on described line of reference
Position and the size of position mark are real time updated, and, it is described that described one-dimensional bone density is distributed
Position and the size of one-dimensional position labelling are real time updated.
8. bone density measurement device as claimed in claim 7, it is characterised in that
On described two dimension bone density image during mobile described line of reference, be equivalent to described two-dimensional position labelling
Motion track two dimensional image be shown.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012-017531 | 2012-01-31 | ||
| JP2012017531A JP5345715B2 (en) | 2012-01-31 | 2012-01-31 | Bone density measuring device |
| PCT/JP2012/079946 WO2013114708A1 (en) | 2012-01-31 | 2012-11-19 | Bone densitometer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104080404A CN104080404A (en) | 2014-10-01 |
| CN104080404B true CN104080404B (en) | 2016-10-12 |
Family
ID=48904770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201280068207.1A Expired - Fee Related CN104080404B (en) | 2012-01-31 | 2012-11-19 | Bone density measurement device |
Country Status (4)
| Country | Link |
|---|---|
| JP (1) | JP5345715B2 (en) |
| CN (1) | CN104080404B (en) |
| TW (1) | TWI508706B (en) |
| WO (1) | WO2013114708A1 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105451655B (en) * | 2013-09-02 | 2018-05-04 | 株式会社日立制作所 | X-ray imaging device and radioscopy method for displaying image |
| JP6267938B2 (en) * | 2013-11-11 | 2018-01-24 | 株式会社日立製作所 | Medical X-ray measurement apparatus and boundary determination method |
| JP6129145B2 (en) * | 2014-12-03 | 2017-05-17 | 株式会社日立製作所 | Medical X-ray measuring device |
| JP6758249B2 (en) * | 2017-05-18 | 2020-09-23 | 富士フイルム株式会社 | Image processing equipment, radiation imaging system, image processing method, and image processing program |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09511418A (en) * | 1991-02-13 | 1997-11-18 | ルナー コーポレイション | Method for measuring bone mineral density around a prosthesis |
| JP2005192657A (en) * | 2003-12-26 | 2005-07-21 | Aloka System Engineering Co Ltd | X-ray ct equipment |
| DE102004026524A1 (en) * | 2004-05-25 | 2005-12-22 | Aesculap Ag & Co. Kg | Bone based coordinate system determination procedure for implant operations uses tomography to find centre of gravity and load paths |
| CN101422371A (en) * | 2007-10-19 | 2009-05-06 | 通用电气公司 | Image guided acquisition of quantitative dual energy data |
| CN101835425A (en) * | 2007-10-23 | 2010-09-15 | 阿洛卡株式会社 | Bone measurement device and bone image processing method |
| CN102217944A (en) * | 2010-04-14 | 2011-10-19 | 深圳市艾克瑞电气有限公司 | Method for measuring bone density by using digital X-ray images |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5306306A (en) * | 1991-02-13 | 1994-04-26 | Lunar Corporation | Method for periprosthetic bone mineral density measurement |
| AU772012B2 (en) * | 1998-09-14 | 2004-04-08 | Board Of Trustees Of The Leland Stanford Junior University | Assessing the condition of a joint and preventing damage |
| JP5324779B2 (en) * | 2005-02-09 | 2013-10-23 | 株式会社日立メディコ | Image diagnosis support system, image diagnosis support program |
| JP4801567B2 (en) * | 2006-11-27 | 2011-10-26 | 株式会社リガク | Volume data processing apparatus and volume data processing program |
-
2012
- 2012-01-31 JP JP2012017531A patent/JP5345715B2/en not_active Expired - Fee Related
- 2012-11-19 WO PCT/JP2012/079946 patent/WO2013114708A1/en active Application Filing
- 2012-11-19 CN CN201280068207.1A patent/CN104080404B/en not_active Expired - Fee Related
- 2012-12-28 TW TW101150949A patent/TWI508706B/en not_active IP Right Cessation
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09511418A (en) * | 1991-02-13 | 1997-11-18 | ルナー コーポレイション | Method for measuring bone mineral density around a prosthesis |
| JP2005192657A (en) * | 2003-12-26 | 2005-07-21 | Aloka System Engineering Co Ltd | X-ray ct equipment |
| DE102004026524A1 (en) * | 2004-05-25 | 2005-12-22 | Aesculap Ag & Co. Kg | Bone based coordinate system determination procedure for implant operations uses tomography to find centre of gravity and load paths |
| CN101422371A (en) * | 2007-10-19 | 2009-05-06 | 通用电气公司 | Image guided acquisition of quantitative dual energy data |
| CN101835425A (en) * | 2007-10-23 | 2010-09-15 | 阿洛卡株式会社 | Bone measurement device and bone image processing method |
| CN102217944A (en) * | 2010-04-14 | 2011-10-19 | 深圳市艾克瑞电气有限公司 | Method for measuring bone density by using digital X-ray images |
Also Published As
| Publication number | Publication date |
|---|---|
| TWI508706B (en) | 2015-11-21 |
| JP2013153999A (en) | 2013-08-15 |
| JP5345715B2 (en) | 2013-11-20 |
| WO2013114708A1 (en) | 2013-08-08 |
| CN104080404A (en) | 2014-10-01 |
| TW201332519A (en) | 2013-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3627145B1 (en) | Data processing device and data processing method for x-ray inspection, and x-ray inspecting device with said device mounted therein | |
| Miranda et al. | Static and dynamic error of a biplanar videoradiography system using marker-based and markerless tracking techniques | |
| US20140050378A1 (en) | System and method for detecting materials or disease states using multi-energy computed tomography | |
| US20180020999A1 (en) | Bone analyzing device | |
| CN103442638B (en) | Bone density measurement device | |
| CN104080404B (en) | Bone density measurement device | |
| JP6798587B2 (en) | Diagnostic imaging system | |
| JP6144915B2 (en) | Biological tissue image reconstruction method, acquisition method and apparatus | |
| Yang et al. | Improving bone strength prediction in human proximal femur specimens through geometrical characterization of trabecular bone microarchitecture and support vector regression | |
| CN103002809A (en) | Adaptive road mapping | |
| IL275987B1 (en) | Methods for x-ray imaging of a subject using multipleenergy decomposition | |
| US7801350B2 (en) | Method for absorptiometry radiographic imaging using a three-dimensional generic model | |
| Sarve et al. | Extracting 3D information on bone remodeling in the proximity of titanium implants in SRμCT image volumes | |
| EP2995255B1 (en) | Image processing device | |
| KR20210028559A (en) | Image analyzing method, image processing apparatus, bone mineral density measuring apparatus and learning model creation method | |
| CN104939848A (en) | Generation of monochromatic images | |
| CN107007294A (en) | X-ray imaging apparatus and bone density measurement method | |
| CN109567837A (en) | The system and method for the finite element analysis of DXA TOMO based on bone | |
| CN104039225A (en) | Medical device | |
| Garson et al. | Characterization of speckle in lung images acquired with a benchtop in-line x-ray phase-contrast system | |
| Subramanian et al. | Quantitative evaluation of bone microstructure using high-resolution extremity cone-beam CT with a CMOS detector | |
| Holcombe et al. | The accuracy of local rib bone geometry measurement using full body CT images | |
| JPH09294740A (en) | Osteometry | |
| Rubinacci et al. | Application of high resolution pQCT analysis for the assessment of a bone lesion: A technical note | |
| Lima et al. | Micro imaging analysis for osteoporosis assessment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20161110 Address after: Tokyo, Japan, Japan Patentee after: Hitachi Ltd. Address before: Tokyo, Japan Patentee before: Hitachi Aloka Medical Ltd. |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20161114 Address after: Tokyo, Japan, Japan Patentee after: Hitachi Ltd. Address before: Tokyo, Japan Patentee before: Hitachi Aloka Medical Ltd. |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161012 Termination date: 20171119 |