CN104837410B - Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means - Google Patents
Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means Download PDFInfo
- Publication number
- CN104837410B CN104837410B CN201380063141.1A CN201380063141A CN104837410B CN 104837410 B CN104837410 B CN 104837410B CN 201380063141 A CN201380063141 A CN 201380063141A CN 104837410 B CN104837410 B CN 104837410B
- Authority
- CN
- China
- Prior art keywords
- signal
- echo
- state
- subchondral bone
- soft tissue
- 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.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0875—Detecting organic movements or changes, e.g. tumours, cysts, swellings 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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Problem is, detects the boundary face of soft tissue and cartilage exactly in noninvasive mode.Solution is that echo data (S101, S102) is obtained under the 1st state and the 2nd state.According to the echo data of the 1st state and the 2nd state, the echo data (S103) of the 1st state of detection and the subchondral bone (911) under the 2nd state.According to the subchondral bone echo data of the 1st state and the 2nd state, the mobile vector (S104) of subchondral bone echo data is detected.The contraposition (S105) of the sampled data of the comparison other of the echo data of the 1st state and the echo data of the 2nd state is carried out based on mobile vector.Calculate the coefficient correlation (S106) by aligning the echo data of the 1st state as identical comparison other sample position and the echo data of the 2nd state.It is soft tissue by the low region decision of coefficient correlation, and detect their boundary face (S107) using the high region of coefficient correlation as cartilage.
Description
Technical field
The present invention relates to a kind of soft group detected with the ultrasonic wave from outside to the boundary face of cartilage and soft tissue
Knit cartilage boundary detection method.
Background technology
In the past, it is proposed that the various devices for generating the information for diagnosing cartilage state.For example, in the super of patent document 1
In sound wave diagnostic device, the detector for receiving and dispatching ultrasonic wave is connected to the surface of knee, knee is come from what is obtained with the detector
Echo-signal inside lid, diagnoses the state of cartilage.That is, the diagnostic ultrasound equipment of patent document 1 in noninvasive mode to soft
The state of bone is diagnosed.Also, in the diagnostic ultrasound equipment of patent document 1, with the electricity of the echo-signal of depth direction
The difference detection cartilage of flat (intensity).
Citation
Patent document
Patent document 1:JP 2010-305 publications
The content of the invention
Problems to be solved by the invention
But, in the device and method of patent document 1, according to the echo-signal level of soft tissue (muscle or skin)
And the difference detection cartilage surface of the echo-signal level of cartilage.Therefore, if the level of echo-signal does not have between soft tissue and cartilage
There is difference, then can not detect cartilage surface exactly.
Also, generally with conventional ultrasonic signal, boundary face of the echo-signal level in soft tissue and cartilage surface
On not change drastically and exactly, but echo-signal level does not have big difference in the boundary face.Thus, with conventional
Method, it is impossible to detect cartilage surface exactly.
It is an object of the invention to provide a kind of soft tissue and cartilage surface of can exactly being detected in noninvasive mode
The soft tissue cartilage boundary face detection method of boundary face.
The means solved the problems, such as
The present invention is the invention of the soft tissue cartilage boundary face detection method for the boundary face for being related to detection soft tissue and cartilage
And with following characteristics.There is soft tissue cartilage boundary face detection method the 1st echo-signal transmitting-receiving process, the 2nd echo-signal to receive
Send out process, subchondral bone detection process, mobile vector detection process and boundary face detection process.
1st echo-signal receives and dispatches process and sends ultrasonic signal into detected body under the 1st state and obtain the 1st echo
Signal.2nd echo-signal receive and dispatch process under the 2nd state into detected body send ultrasonic signal and obtain the 2nd echo letter
Number.
Subchondral bone detection process detects the subchondral bone echo-signal of the 1st state from the 1st echo-signal, from the 2nd time
The subchondral bone echo-signal of the 2nd state is detected in ripple signal.
Mobile vector detects that process is believed from the subchondral bone echo-signal of the 1st state and the subchondral bone echo of the 2nd state
The mobile vector from the 1st state to the subchondral bone of the 2nd state is detected in number.
Boundary face detects that process is carried out based on mobile vector to the sample position of the 1st state and the echo-signal of the 2nd state
Correction, so as to detect the boundary face of soft tissue and cartilage.
In the method, cartilage make use of to be attached to subchondral bone, soft tissue is non-cohesive on cartilage, and soft tissue can be soft
The situation that bone surface slides laterally.
If make transmission ultrasonic signal detector contact moved in detected body, soft tissue is followed, cartilage with
And subchondral bone is not followed.Therefore, change in location and cartilage position relative to detector of the subchondral bone relative to detector
Change is consistent, and equivalent to mobile vector, but change in location of their change in location with soft tissue relative to detector differs
Cause.In addition, in the state of detector is contacted detected body, making soft tissue with cartilage and having the detected body of cartilage
Lateral bend, similarly, change in location change in location phase with cartilage relative to detector of the subchondral bone relative to detector
Unanimously, subchondral bone is inconsistent relative to the change in location of detector relative to the change in location and soft tissue of detector.
Therefore, if correcting the sample position of the echo-signal of the 1st state and the echo-signal of the 2nd state according to mobile vector
Sample position, and the echo-signal of each sample position is compared, then between soft tissue and cartilage (and subchondral bone)
Comparative result is different.By using the difference, soft tissue and cartilage can be differentiated, the side of soft tissue and cartilage can be also detected
Interface.
In addition, in the soft tissue cartilage boundary face detection method of the present invention, subchondral bone detects process along depth direction
The signal intensity of the 1st echo-signal is obtained successively, will detect that signal more than the subchondral bone detection threshold value is strong
The signal of the scope of degree is detected as the subchondral bone echo-signal of the 1st state.Subchondral bone detection process detection institute
State the subchondral bone echo-signal of the 1st state and the similar degree of the 2nd echo-signal, and by similar degree highest echo-signal
Detected as the subchondral bone echo-signal of the 2nd state.
In addition, in the soft tissue cartilage boundary face detection method of the present invention, subchondral bone detects process along depth direction
Obtain the signal intensity of the 1st echo-signal successively since the side of deep, the subchondral bone detection threshold value will be detected
The signal of the scope of signal intensity above is detected as the subchondral bone echo-signal of the 1st state.Subchondral bone is examined
Survey process and detect the subchondral bone echo-signal of the 1st state and the similar degree of the 2nd echo-signal, and by similar degree most
High echo-signal is detected as the subchondral bone echo-signal of the 2nd state.
The specific detection method of subchondral bone is shown in these methods.
Invention effect
In accordance with the invention it is possible to ultrasonic signal be sent from the outside of the detected bodys such as knee, in the outside of detected body
Receive the echo-signal, and the boundary face of detection soft tissue and cartilage surface exactly.Thereby, it is possible to detect exactly from soft
Bone and the echo come, and can effectively utilize the diagnosis in cartilage.
Brief description of the drawings
Fig. 1 is the block diagram for representing to be related to the structure of the soft tissue cartilage boundary detection means 10 of embodiments of the present invention.
Fig. 2 is the phase of detector 100 for representing to be related to the soft tissue cartilage boundary detection means 10 of embodiments of the present invention
For the figure of the set-up mode of detected body.
Fig. 3 is the figure for illustrating the detection concept for the cartilage surface for being related to embodiments of the present invention.
Fig. 4 is the flow chart for the soft tissue cartilage boundary face detection method for being related to embodiments of the present invention.
Fig. 5 is the figure for the waveform example for representing each echo-signal under the 1st state [T1] and the 2nd state [T2].
Fig. 6 is the oscillogram for illustrating the detection concept for the mobile vector for being related to present embodiment.
Fig. 7 is the figure for representing to be related to the definition of the mobile vector of present embodiment.
Fig. 8 is the figure for representing to be related to the distribution of the mobile vector of present embodiment.
Fig. 9 is for illustrating that detection echo data is to belong to method (the 1st side that cartilage 901 still falls within soft tissue 903
Method) oscillogram.
Figure 10 is for illustrating that detection echo data is to belong to method (the 2nd side that cartilage 901 still falls within soft tissue 903
Method) oscillogram.
Figure 11 is the figure for representing the detection structure based on the mechanical scan for moving oscillator.
Embodiment
Be explained with reference to be related to embodiments of the present invention soft tissue cartilage boundary face detection method and soft group
Knit cartilage boundary face detection means.Fig. 1 represents to be related to the soft tissue cartilage boundary detection means 10 of embodiments of the present invention
The block diagram of structure.Fig. 2 for represent be related to present embodiment soft tissue cartilage boundary detection means 10 detector 100 relative to
The figure of the set-up mode of detected body, Fig. 2 (A) represents the situation of the 1st state (T=T1), and Fig. 2 (B) represents the 2nd state (T=
T2 situation).In addition, in the following description, the example for moving detector 100 is illustrated, but to mobile detected body
Situation, can also be applicable following method or structure.For example, as abutting detector 100 simultaneously on the knee as detected body
Fixed, knee is bent and stretched such situation can also be applicable.As long as that is, soft group between the 1st state [T1] and the 2nd state [T2]
Knitting the method changed with the position relationship of cartilage and structure just can be applicable.
Fig. 3 is the figure for illustrating the detection concept for the cartilage surface for being related to embodiments of the present invention, and Fig. 3 (A) is represented
1st state [T1] (T=T1), Fig. 3 (B) represents the 2nd state [T2] (T=T2).Fig. 3 is the area by have sent ultrasonic signal
The figure that the surface replacement of domain and its near zone is observed for flat plane.
Soft tissue cartilage boundary face detection means 10 possess operating portion 11, transmission control unit 12, echo signal reception portion 13,
Data analysis portion 14 and detector 100.Transmission control unit 12, echo signal reception portion 13 and detector 100 are equivalent to this
" receiving and transmitting part " of invention.
Operating portion 11 accepts the operation input of user.For example, operating portion 11 possesses multiple operator (not shown), according to
Operation of the family to operator, indicates that the processing for detecting cartilage surface starts to perform to transmission control unit 12.
Transmission control unit 12 generates the ultrasonic wave letter that the carrier wave being made up of the frequency of ultrasonic wave is shaped as to pulse type
Number.Transmission control unit 12 generates ultrasonic signal respectively under the 1st state [T1] and the 2nd state [T2].
Transmission control unit 12 exports ultrasonic signal to detector 100.Detector 100 possesses parallel with transmitting-receiving corrugated
Direction on multiple oscillators (reference picture 3) for configuring.The configuration direction of the oscillator is scanning direction.Each oscillator is into detected body
Send the ultrasonic signal being made up of defined transmission field angle.Each oscillator sends ultrasonic signal with predetermined time interval,
And receive its reflection echo signal.
As shown in Figure 2, detector 100 is connected to as detected body detail by the end face of transmitting-receiving corrugated side
Knee soft tissue 903 surface on mode and configure.Herein, as shown in Figure 3, so-called soft tissue 903 is bag
Internal body portion containing skin and muscle, is present in than cartilage 901 closer to the position of the face side of detected body.Cartilage 901
It is attached on subchondral bone 911, subchondral bone 911 is the tissue combined with bone (cancellous bone) 902.
While contacting detector 100 and the surface of soft tissue 903 as shown in Fig. 2 (A), while as shown in Fig. 2 (B)
Detector 100 is moved along surface like that.Thus, as shown in Figure 2, the one side of soft tissue 903 is slided on the surface of cartilage 901
It is dynamic, while following detector 100 and moving.The state for making Fig. 2 (A) before the detector 100 movement is the 1st state (t=T1),
The state for making Fig. 2 (B) after the movement of detector 100 is the 2nd state (t=T2).Now, row of the detector 100 along oscillator is made
Column direction (scanning direction) and move.
Each oscillator sends ultrasonic signal into detected body under the 1st state [T1] and the 2nd state [T1] respectively.This
When, each oscillator of detector 100 turns into the central axis direction for sending wave beam with the direction vertical relative to the surface of soft tissue 903
Mode send ultrasonic signal.
Soft tissue 903 of each oscillator received ultrasonic signal of detector 100 in detected body, cartilage 901 and soft
Echo-signal after being reflected at bone sending down the fishbone 911, and export to echo signal reception portion 13.Detector 100 will be by the 1st state
The 1st echo group SW [T1] that each oscillator is obtained in [T1] echo-signal is constituted with by the 2nd state [T2] each oscillator obtain
Echo-signal constitute the 2nd echo-signal group SW [T2] export respectively to echo signal reception portion 13.
Each echo-signal in 13 pairs of echo signal reception portion carries out defined enhanced processing and exported to data analysis portion 14.Return
Ripple signal receiving part 13 is individually believed the 1st echo group SW [T1] each echo-signal and the 2nd echo crowd SW [T2] each echo
Number it is amplified and handles and export to data analysis portion 14.
Data analysis portion 14 possesses AD transformation components 141, storage part 142 and determination unit 143.AD transformation components 141 are by advise
Fixed time interval makes data discrete to echo-signal sampling.Echo-signal after the data discrete is echo data.Thus,
(sweep) can be scanned by each distance to obtain carrying out the echo data after data sampling at regular intervals in the depth direction.
That is, the echo data being distributed in the 2 dimensional region based on scanning direction and depth direction can be obtained.Hereinafter, by the number of echoes
The echo data of Two dimensional Distribution is referred to as according to group.AD transformation components 141 export each echo data to storage part 142.
Storage part 142 stores each echo data exported by AD transformation components 141.Storage part 142 possesses under the 1st state of storage
The capacity of the echo data of obtained Two dimensional Distribution and the echo data of the Two dimensional Distribution obtained in the 2nd state.
Determination unit 143 detects the boundary face of soft tissue and cartilage, aftermentioned specific processing side according to the flow shown in Fig. 4
Method.Fig. 4 is the flow chart for the soft tissue cartilage boundary face detection method for being related to embodiments of the present invention.
Determination unit 143 reads and obtains the number of echoes for the Two dimensional Distribution obtained under the 1st state for being stored in storage part 142
According to the echo data (S101, S102) with the Two dimensional Distribution obtained under the 2nd state.
Determination unit 143 detects the subchondral bone 911 under the 1st state according to the echo data of the Two dimensional Distribution of the 1st state
Echo data (equivalent to the present invention " the subchondral bone echo-signal of the 1st state ").Determination unit 143 is according to the 2nd state
The echo data of Two dimensional Distribution detects the echo data of the subchondral bone 911 under the 2nd state (equivalent to " the 2nd shape of the present invention
The subchondral bone echo-signal of state ") (S103).
Determination unit 143 is detected according to the subchondral bone echo data of the 1st state and the subchondral bone echo data of the 2nd state
Go out the mobile vector (S104) of subchondral bone echo data.
Determination unit 143 is based on mobile vector for the echo data of the Two dimensional Distribution of the 1st state and the two dimension point of the 2nd state
The echo data of cloth is compared the contraposition (S105) of the sample data of object.
Determination unit 143 calculate by align as identical comparison other sample position the 1st state echo data with
Coefficient correlation between the echo data of 2nd state.Specifically, the setting of determination unit 143 comprising comparison other sample position and
The comparison other region being made up of in the depth direction the region of Rack.Determination unit 143 is calculated in comparison other region
The waveform that constitutes of echo data of the 1st state and the waveform that constitutes of the echo data of the 2nd state in comparison other region
Coefficient correlation (S106).
Determination unit 143 detects that the region that the high comparison other sample position of coefficient correlation gathered and coefficient correlation are low
The region that comparison other sample position is gathered, and detect the border (S107) in this 2 regions.
As described above, the move mode of cartilage 901 is identical with subchondral bone 911, the move mode of soft tissue 903 with it is soft
Bone sending down the fishbone 911 is different.Thus, after the contraposition based on mobile vector, by what is constituted with the identical move mode of subchondral bone 911
The coefficient correlation of the echo data of cartilage 901 increases.On the other hand, after the contraposition based on mobile vector, by with subchondral bone
The coefficient correlation reduction of the echo data for the soft tissue 903 that 911 different move modes are constituted.Therefore, detected in step S107
The boundary face arrived turns into the boundary face of soft tissue 903 and cartilage 901.Like this, by that can detect soft using above-mentioned processing
The boundary face of tissue 903 and cartilage 901.
If in addition, detecting the surface (boundary face of soft tissue 903 and cartilage 901) of cartilage 901, cartilage (not shown)
Partial echo data generation of the diagnosis information generation unit based on cartilage 901 can be used in the information of cartilage degeneration diagnosis.Specifically
For, cartilage diagnosis information generation unit obtains the echo data and subchondral bone near cartilage surface in different multiple periods
Echo data formed by group.Cartilage diagnosis information generation unit according to the composition of the group of these echo datas when during change
Change, detection variable quantity as caused by the denaturation of cartilage surface etc. also may be used.Cartilage diagnosis information generation unit makees the testing result
The information diagnosed to can be used in cartilage degeneration is exported.
Next, referring to the drawings come more specifically describe with data analysis portion 14 perform soft tissue and cartilage border
The detection method in face.In addition, in order that explanation is simple, the detector 100 between the 1st state [T1] and the 2nd state [T2] is (each
Oscillator) displacement Δ x be set to the configuration space of oscillator consistent and illustrate.
First, as the 1st state [T1], for example, in the state of being bent as the knee of detected body with the 1st angle,
Detector 100 is set to be connected on the surface of knee.In other words, detector 100 is made to be connected on the surface of soft tissue 903.This is
Fig. 3 (A) state.
Multiple oscillators on detector 100 are configured respectively in the direction parallel with the surface of soft tissue 903 (with transmitting-receiving ripple
The parallel scanning direction in face) on be configured at regular intervals.Multiple oscillators send ultrasonic wave to the direction vertical with scanning direction
Signal.If Fig. 3 example, five oscillators are equally spaced configured along scanning direction in detector 100, such as Fig. 3 (A) institute
Show like that, each oscillator for being configured in respective position sends ultrasonic signal to the direction vertical with the surface of soft tissue 903.From
The ultrasonic signal that each allocation position is sent reflects in each depth location of soft tissue 903, cartilage 901, subchondral bone 911,
Received, sampled in data analysis portion 14 by each oscillator.The 1st state [T1] receive echo data SWT11, SWT12,
SWT13, SWT14, SWT15 echo data group turn into the 1st echo group SW [T1].
Next, making detector 100 keep being connected in the state of soft tissue 903, parallel with the surface of soft tissue 903
Direction and the direction parallel with scanning direction on by the displacement Δ x amounts of detector 100.The state is the 2nd state [T2], also
For Fig. 3 (B) state.
Now, soft tissue 903 is followed the movement of detector 100 and moved.Therefore, the transmitting-receiving corrugated of detector 100 with it is soft
The relative position relationship of each position of the scanning direction of tissue 903 does not produce change with the movement of detector 100.
On the other hand, because cartilage 901 is fixed on bone 902 by subchondral bone 911, even if there is the shifting of detector 100
Dynamic, cartilage 901 is not also moved.Therefore, the transmitting-receiving corrugated and the scanning direction of cartilage 901 and subchondral bone 911 of detector 100
The relative position relationship of each position change with the movement of detector 100.
After the 2nd state, as shown in Fig. 3 (B), from each oscillator of detector 100 to the surface with soft tissue 903
Vertical direction (direction vertical with transmitting-receiving corrugated (scanning direction)) sends ultrasonic signal.Sent from each allocation position
Ultrasonic signal reflects in each depth location of soft tissue 903, cartilage 901, subchondral bone 911, is received by each oscillator, in data
Analysis unit 14 is sampled.Echo data SWT21, SWT22, SWT23, SWT24, SWT25 for being received in the 2nd state [T2] echo
Data group turns into the 2nd echo group SW [T2].
Like this, detector 100 movement before, obtain by multiple echo data SWT11, SWT12, SWT13, SWT14,
The 1st echo group SW [T1] that SWT15 is constituted.Then, detector 100 movement after, obtain by multiple echo data SWT21,
The 2nd echo group SW [T2] that SWT22, SWT23, SWT24, SWT25 are constituted.
Fig. 5 is expression in the 1st state [T1] and the figure of each echo waveform under the 2nd state [T2].In addition, in order in figure
More understandable earth's surface shows the feature of the present invention in 5, if the interval of distance (amount of movement) Δ x and each oscillator that detector 100 is moved, i.e.
The interval of scan position is equal.In addition, the detection in the soft tissue cartilage boundary face of following explanation under this condition.
Even if (i) cartilage 901 and the detector 100 of subchondral bone 911 are moved, cartilage 901 and subchondral bone 911 are not yet
It is mobile.Therefore, if the displacement Δ x amounts of detector 100, the position (each scan position) of each oscillator of detector 100 and cartilage
901 and subchondral bone 911 each position position relationship along scanning direction deviate amount of movement Δ x amounts.
Now, as Fig. 5 the 1st state [T1] each echo waveform and the 2nd state [T2] each echo waveform shown in that
Sample, the 1st echo group SW [T1] cartilage 901 and subchondral bone 911 region echo data SWT11 not with the 2nd echo group
SW [T2] echo data SWT21 cartilage 901 and the waveform in the region of subchondral bone 911 are consistent, but with the 2nd echo group
SW [T2] echo data SWT22 cartilage 901 and the waveform in the region of subchondral bone 911 are about consistent.
Similarly, in echo data SWT12 cartilage 901 and the region of subchondral bone 911 and in echo data SWT23
Cartilage 901 and subchondral bone 911 region, waveform is about consistent.In echo data SWT13 cartilage 901 and cartilage
The region of sending down the fishbone 911 and the cartilage 901 in echo data SWT24 and the region of subchondral bone 911, waveform are about consistent.
Echo data SWT14 cartilage 901 and the region of subchondral bone 911 and the cartilage 901 and cartilage in echo data SWT25
The region of sending down the fishbone 911, waveform is about consistent.
Therefore, in cartilage 901 and subchondral bone 911, the echo data of each scan position is in the 1st state [T1] and the
Amount of movement Δ x between 2 states [T2], i.e., about one in the state of scan position deviate from an interval by the configuration space of oscillator
Cause.That is, in cartilage 901 and subchondral bone 911, the sample position of detector 100 as each echo data of benchmark is existed
1st state [T1] and change amount of movement Δ x under the 2nd state [T2].Like this, cartilage 901 turns into identical with subchondral bone 911
Move mode.
(ii) soft tissue 903 is as described above, detector 100 is connected on the surface of soft tissue 903, and soft tissue 903 does not have
On the surface for being fixed on cartilage 901.Therefore, if detector 100 moves amount of movement Δ x, soft tissue 903 also follows detector
100 movement, mobile amount of movement Δ x.
Now, as Fig. 5 the 1st state [T1] each echo waveform and the 2nd state [T2] each echo waveform shown in that
Sample, the echo data SWT11 in the region of the 1st echo group SW [T1] soft tissue 903 and the 2nd echo group SW [T2] soft tissue
The echo data SWT21 in 903 region waveform is about consistent.
Similarly, echo data SWT12 soft tissue 903 region and the soft tissue 903 in echo data SWT22
Region, waveform is about consistent.In the region of echo data SWT13 soft tissue 903 and the soft tissue in echo data SWT23
903 region, waveform is about consistent.Echo data SWT14 soft tissue 903 region with the soft of echo data SWT24
The region of tissue 903, waveform is about consistent.In the region of echo data SWT15 soft tissue 903 and in echo data SWT25
Soft tissue 903 region, waveform is about consistent.
Therefore, in soft tissue 903, the echo data of each scan position is in the 1st state [T1] and in the 2nd state [T2]
Under, the position along the scanning direction of correspondence detector 100 is about consistent.That is, in soft tissue 903, using detector 100 as
The sample position of each echo data of benchmark does not produce change in the 1st state [T1] and between the 2nd state [T2].Thus, soft group
Knit the sample position of the move mode and cartilage 901 described above or the such echo data of subchondral bone 911 of 903 echo data
The move mode that change is produced with the movement of detector 100 is different.
Like this, in subchondral bone 911 and cartilage 901 ad-hoc location (the 1st state [T1] with the 2nd state
[T2] is identical) place reflection echo data the sample position using detector 100 as benchmark the 1st state [T1] with the 2nd
Between state [T2] amount of movement Δ x amounts are moved along scanning direction.
On the other hand, ad-hoc location (identical with the 2nd state [T2] in the 1st state [T1]) place in soft tissue 903 is anti-
The sample position using detector 100 as benchmark for the echo data penetrated is not produced in the 1st state [T1] and between the 2nd state [T2]
Changing.
Using the characteristic, with shown in the soft tissue cartilage boundary face detection means step S104 described above of present embodiment
Like that, the mobile vector of the echo data of detection subchondral bone 911.Then, shown in step S105 described above like that, moved with this
The sample position and the sample position of the echo data of the 2nd state [T2] of the echo data of dynamic vector the 1st state [T1] of progress
Contraposition.In other words, the sample position and the 2nd state [T2] of the echo data of the 1st state [T1] of comparison other are decided to become
The combination of the sample position of echo data.
Now, first, shown in step S103 described above like that, the echo data of subchondral bone 911 is detected.Such as Fig. 4 institutes
Show like that, the echo data of subchondral bone 911 turns into high echo level (amplitude).The echo level at cartilage surface of cartilage 901
Although increasing, the region between the cartilage surface and subchondral bone 911, echo level reduction.
Therefore, first, the 1st echo group SW [T1] echo data SWT11-SWT15 time is obtained successively along depth direction
Ripple level.Then, in the range of prescribed depth (suitably being set according to the thickness of cartilage), detection echo level is less than subchondral bone
The situation of detection threshold value, thereafter, echo level more than detection subchondral bone detection threshold value.Thereby, it is possible to by the sample
Detected after position as the echo data of subchondral bone 911.
Next, obtain the echo data of the subchondral bone 911 of the 2nd state [T2], and detect from the 1st state [T1] to the
The mobile vector of the echo data of the subchondral bone 911 of 2 states [T2].Fig. 6 is the movement for illustrating to be related to present embodiment
The oscillogram of the detection concept of vector.In addition, waveform is identical with Fig. 5 oscillogram in itself.
To echo data group's setting region-of-interest Z of the subchondral bone 911 obtained in the 1st state [T1]T1CS.Now, close
Note region ZT1CSProvided with the length of depth direction (time orientation).Specifically, for example, as shown in Figure 6, for by
The 2 dimensional region constituted according to each echo data SWT11-SWT15 scanning directions determined with depth direction, will be along 1
The defined scope of the depth direction (time orientation) of echo data is set to region-of-interest ZT1CS.It is somebody's turn to do next, extracting and being included in
Region-of-interest ZT1CSInterior echo data.This is, in other words, equivalent to region-of-interest ZT1CSThe ripple of the echo data group extracted
Shape.
Next, setting region of search Z to the 2nd echo group SW [T2] echo data SWT21-SWT25R1.Region of search
ZR1With relative to detector 100 by region-of-interest ZT1CSPosition as benchmark and than region-of-interest ZT1CSInclude bigger scope
The mode of echo data determine.
Specifically, according to will be relative to detector 100 and region-of-interest Z under the 2nd state [T2]T1CSIdentical position
As depth direction and scanning direction center and expanded on both depth direction and scanning direction after prescribed limit
Region set.For example, as shown in Figure 6, in the depth direction, in echo-signal SWT12 depth location P1CSSetting
Region-of-interest ZT1CSCenter when, in the echo-signal SWT22 of the 2nd state [T2] depth location P1CSSet region of search
ZR1Center.Then, the scope longer than the length of the depth direction of region-of-interest is set to region of search ZR1Depth direction
Scope.In addition, depth location P1CSAs long as in region of search ZR1It is interior, it is not region of search ZR1Center also may be used.
In addition, setting region-of-interest Z to echo data SWT12 in a scanning directionT1CSWhen, by the center of scanning direction
As echo data SWT22 and comprising echo data SWT21, SWT22, SWT23 mode sets region of search ZR1Scanning side
To scope.
Next, to echo data group's setting sending down the fishbone comparison other of the subchondral bone 911 obtained in the 2nd state [T2]
Region ZT2CS.Now, as shown in Figure 6, by sending down the fishbone comparison other region ZT2CSIt is set to 1 echo-signal of selection and at it
The region of specific length on depth direction.Sending down the fishbone comparison other region ZT2CSThe length of depth direction be set to and pay close attention to
Region ZT1CSDepth direction length it is consistent.
Sending down the fishbone comparison other region ZT2CSIn the region of search Z set as described aboveR1Whole region in be set.
By each sending down the fishbone being set like this comparison other region ZT2CSObtain echo data.
Next, to region-of-interest ZT1CSThe echo data of echo data and sending down the fishbone comparison other region ZT2CS enter
Row relevant treatment, and calculate coefficient correlation.
Change sending down the fishbone comparison other region Z successively on one sideT2CSPosition, while in region of search ZR1Whole region enter
Row region-of-interest ZT1CSEcho data and sending down the fishbone comparison other region ZT2CSThe coefficient correlation of echo data calculate.
Towards region of search ZR1Whole region, perform such a for a region-of-interest ZT1CSCoefficient correlation
Calculate.
Next, by each region-of-interest ZT1CSDetect the sending down the fishbone comparison other region Z that coefficient correlation is maximumT2CS。
This turns into the echo of the subchondral bone 911 under the 2nd state [T2] equivalent to region most like in the region of subchondral bone 911
Data.The echo data of the subchondral bone 911 under the 2nd state [T2] is detected according to the similarity of echo-signal like this.
Next, mobile arrow of the detection from the 1st state [T1] to the echo data of the subchondral bone 911 of the 2nd state [T2]
Measure (moving direction, amount of movement).
Herein, so-called mobile vector can be defined as follows.Fig. 7 is related to the mobile vector of present embodiment for expression
Definition figure.In addition, representing the situation using region-of-interest with sending down the fishbone comparison other region in the figure 7.Mobile vector vm
By be judged as the representative point of mutual most similar region-of-interest and sending down the fishbone with the position of the representative point in comparison other region come
Definition.Mobile vector vmFor the representative of region-of-interest point, as starting point, sending down the fishbone is made with the representative point in comparison other region
For the vector of terminal, defined according to moving direction and amount of movement.
Specifically, for example, in the example of fig. 7, the representative point of region-of-interest is scanning side on echo-signal SWT11
To on the P1 of position, it is present in defined depth location.In addition, sending down the fishbone with the representative point in comparison other region in echo-signal
It is on the P2 of scanning direction position, to be present in defined depth location on SWT22.The depth location of the representative point of region-of-interest is with
Bone is identical with the depth location of the representative point in comparison other region.
Now, mobile vector vmIt is movement as scanning direction is set as moving direction and scan position P1, P2 interval delta x
The vector of amount.
Such a calculating for mobile vector can both be carried out with whole echo datas of subchondral bone 911, can also
Carried out with one or more echo datas of representative.In the case of being carried out with whole echo datas or multiple echo datas, it will divide
The average value for the mobile vector not calculated as subchondral bone 911 mobile vector.
When calculating mobile vector like this, the mobile vector of each echo data of cartilage 901 and soft tissue 903 is Fig. 8
Shown distribution.Fig. 8 is related to the figure of the distribution of the mobile vector of present embodiment for expression.As shown in Figure 8, if soft
The mobile vector v of the echo data in the region of bone sending down the fishbone 911mFor Δ x, the then echo data in the region of cartilage 901 mobile vector vm
For Δ x.But, the mobile vector vm of the echo data of soft tissue 903 is 0, the movement with subchondral bone 911 and cartilage 901
Vector is different.
Using this characteristic, detected with the method shown in next each echo data be belong to cartilage 901 still fall within it is soft
Tissue 903.Fig. 9 is for illustrating that detection echo data is to belong to method (the 1st side that cartilage 901 still falls within soft tissue 903
Method) oscillogram.
To the 1st comparison other region Z of echo data group's setting obtained in the 1st state [T1]T1n(n is by along depth direction
Sample position number determine integer).Now, the 1st comparison other region ZT1nAdvised with the length of depth direction (time orientation)
It is fixed.As specific example, for example, as shown in Figure 9, to the scanning by being determined according to each echo data SWT11-SWT15
The 2 dimensional region that direction and depth direction are constituted, by along the defined scope of the depth direction of 1 echo data (time orientation)
It is set as the 1st comparison other region ZT11.The 1st comparison other region Z is included in next, extractingT11Interior echo data.This
For in other words, equivalent to the 1st comparison other region ZT11The waveform of the echo data group extracted.In addition, can both remove soft
The region of bone sending down the fishbone 911 sets the 1st comparison other region ZT11, can also remove in cartilage 901 and close on subchondral bone 911
The region of low echo level set.
Next, to the 2nd comparison other region Z of echo data group's setting obtained in the 2nd state [T2]T2n(n is by edge
The integer that the sample position number of depth direction is determined).Now, the 2nd comparison other region ZT2nBeing set at makes the 1st comparison other
Region ZT1nOn position after being moved by mobile vector.
For more specifically, by with obtaining the 1st comparison other region Z under the 1st state [T1]T1nEcho data
Oscillator identical oscillator, detection under the 2nd state [T2] with the 1st comparison other region ZT1nIdentical depth areas.Then, will
Depth areas of the region as obtained from being moved mobile vector under the 2nd state [T2] is set to be set as the 2nd comparison other region
ZT2n.That is, the sample position and the echo data of the 2nd state [T2] of the echo data of the 1st state [T1] of comparison other are set to
Sample position is aligned according to mobile vector.
It is used as specific example, when obtaining the mobile vector Δ x only moved on scanning direction as described above, such as Fig. 9
It is shown such, if the 1st comparison other region ZT11It is set at the echo data SWT12 of the 1st state [T1] prescribed depth position
P11, then by the 2nd comparison other region ZT21It is set in from the echo data obtained by identical oscillator under the 2nd state [T2]
The prescribed depth position P21 of echo data SWT23 after SWT22 movement mobile vector Δ x amounts.
Similarly, as shown in Figure 9, if the 1st comparison other region ZT12It is set the echo data of the 1st state [T1]
On SWT12 prescribed depth position P12, then the 2nd comparison other region ZT22Be set at from by identical oscillator in the 2nd state
The prescribed depth position P22 of echo data SWT23 after the echo data SWT22 movement mobile vector Δ x amounts obtained under [T2]
On.
Next, calculating the 1st comparison other region ZT1nEcho data and the 2nd comparison other region ZT2nEcho data
Coefficient correlation.
Herein, as described above, cartilage 901 is moved by with the identical mobile vector of subchondral bone 911, and move mode phase
Together.Therefore, the 1st comparison other region Z after being aligned according to mobile vectorT1nEcho data and the 2nd comparison other region
ZT2nThe coefficient correlation of echo data be about 1.That is, coefficient correlation increases.
On the other hand, because soft tissue 903 does not carry out moving with cartilage 901 and the identical of subchondral bone 911, so
Move mode is different.Therefore, the 1st comparison other region Z after being aligned according to mobile vectorT1nEcho data and the 2nd ratio
Compared with subject area ZT2nEcho data coefficient correlation close to 0.That is, coefficient correlation is reduced.
For example, as shown in Figure 9, the 1st comparison other region Z of setting in cartilage 901T12Echo data and the 2nd
Comparison other region ZT22Echo data be identical waveform, coefficient correlation increases.On the other hand, set in soft tissue 903
The 1st comparison other region ZT11Echo data and the 2nd comparison other region ZT21Echo data be different waveforms, it is related
Coefficient is reduced.
Like this, the echo data that can interpolate that as the high comparison other region of coefficient correlation is the echo in cartilage 901
Data, the echo data in the low comparison other region of coefficient correlation is the echo data in soft tissue 903.
Therefore, such a 1st comparison other area is calculated by each echo data beyond the region to subchondral bone 911
Domain ZT11Echo data and the 2nd comparison other region ZT21Coefficient correlation, it is soft tissue 903 to be capable of detecting when each echo data
In interior or cartilage 901.Then, by detecting the echo data being judged as in soft tissue 903 group with being judged as cartilage 901
The border of interior echo data group, can detect the boundary face of soft tissue 903 and cartilage 901.
As previously discussed, the structure by application present embodiment and processing, can exactly be examined in noninvasive mode
Survey the boundary face of soft tissue and cartilage.In addition, as noninvasive method, whole echo datas being set as into region-of-interest, searched
The high comparison other region of coefficient correlation between rope and detection and respective region-of-interest, also can be by region-of-interest and comparison other
The move mode detection soft tissue and the boundary face of cartilage in region.But, structure and processing by application present embodiment,
Due to being not necessarily to carry out the search process to Zone Full, therefore, it is possible to the border of detection soft tissue and cartilage more at high speed
Face.
In addition, in the above description, the echo data that makes to be obtained under the 2nd state [T2] with mobile vector is illustrated
Sample position is moved, and then sets the example in comparison other region, but makes what is obtained under the 1st state [T1] based on mobile vector
The sample position of echo data, which moves and then sets comparison other region, also may be used.
Figure 10 is for illustrating that detection echo data is to belong to method (the 2nd side that cartilage 901 still falls within soft tissue 903
Method) oscillogram.In the 2nd method, the echo data that will be obtained under the 2nd state [T2] is as benchmark.
To the 2nd comparison other region Z of echo data group's setting obtained under the 2nd state [T2]T2n(n is by along depth side
To sample position number determine integer).As specific example, for example, as shown in Figure 10, to by according to each number of echoes
The 2 dimensional region constituted according to the SWT21-SWT25 scanning directions determined and depth direction, by along the depth side of 1 echo data
It is the 2nd comparison other region Z to the defined range set of (time orientation)T21.Compare pair included in the 2nd next, extracting
As region ZT21Interior echo data.This is, in other words, equivalent to the 2nd comparison other region ZT21The echo data group extracted
Waveform.
Next, to the 1st comparison other region Z of echo data group's setting obtained under the 1st state [T1]T1n(n serves as reasons
The integer determined along the sample position number of depth direction).Now, the 1st comparison other region ZT1nIt is set at based on mobile vector
Make the 2nd comparison other region ZT2nPosition after being moved on the direction opposite with mobile vector and according to the size of mobile vector.
For more specifically, by with obtaining the 2nd comparison other region Z under the 2nd state [T2]T2nEcho data
Oscillator identical oscillator, detection under the 1st state [T1] with the 2nd comparison other region ZT2nIdentical depth areas.Then, will
The depth areas under the 1st state [T1] for obtaining region movement based on mobile vector is set as the 1st comparison other region
ZT1n.That is, relative to the 2nd state [T2] for being set to comparison other echo data sample position, based on mobile vector to the 1st
The sample position of the echo data of state [T1] is aligned.
It is used as specific example, when obtaining the mobile vector Δ x only moved on scanning direction as described above, such as Figure 10
It is shown such, if the 2nd comparison other region ZT21It is set at the echo data SWT22 of the 2nd state [T2] prescribed depth position
P21, then by the 1st comparison other region ZT11It is set in from the echo data obtained by identical oscillator under the 1st state [T1]
SWT12 moved in the opposite direction by mobile vector Δ x amounts after echo data SWT11 prescribed depth position P11.
Similarly, as shown in Figure 10, if the 2nd comparison other region ZT22It is set the echo data of the 2nd state [T2]
On SWT22 prescribed depth position P22, then by the 1st comparison other region ZT12Be set in from by identical oscillator in the 1st state
The echo data SWT12 obtained under [T1] moved in the opposite direction by mobile vector Δ x amounts after echo data SWT11 rule
Determine on depth location P12.
Next, calculating the 2nd comparison other region ZT2nEcho data and the 1st comparison other region ZT1nEcho data
Coefficient correlation.
Even if being handled like this, the border of soft tissue and cartilage can be also detected exactly in noninvasive mode
Face.
In addition, in the above description, illustrating by the 1st comparison other region ZT1nEcho data and the 2nd comparison other
Region ZT2nThe coefficient correlation of echo data differentiate the example of soft tissue and cartilage.But, detect the 1st with other methods
The waveform for the echo-signal that each distance of the waveform for the echo-signal that each distance of state [T1] is scanned and the 2nd state [T2] is scanned
Similar degree, differentiate that soft tissue and cartilage also may be used by the similar degree.
As a specific example, there is the method using norm.Norm Norm can be according to the formula next represented
Definition.
[mathematical expression 1]
Herein, Ddzt1 (i) is the 1st comparison other region ZT1nDepth direction position i point echo data values.
Ddzt2 (i) is the 2nd comparison other region ZT2nDepth direction position i point echo data values.K compares pair equivalent to the 1st
As region ZT1nAnd the 2nd comparison other region ZT2nResolution ratio, represent be included in the 1st comparison other region ZT1nAnd the 2nd ratio
Compared with subject area ZT2nEcho data sum.M is the constant suitably set.
Like this, norm Norm is by the 1st comparison other region ZT1nWith the 2nd comparison other region ZT2nSame position
The absolute value of the residual quantity of echo data takes n powers, and obtained from opening n powers after the n power values are added in region entirety
Value.Therefore, the 1st comparison other region ZT1nWith the 2nd comparison other region ZT2nThe more high then norm Norm of similar degree it is smaller, it is similar
The more low then norm Norm of degree is bigger.
Thus, the 1st small comparison other region Z of normT1nWith the 2nd comparison other region ZT2nIt is high equivalent to coefficient correlation
Situation, can interpolate that as in cartilage 901.Also, the 1st big comparison other region Z of normT1nWith the 2nd comparison other region
ZT2nThe situation low equivalent to coefficient correlation, can interpolate that as in soft tissue 903.
In the case of the mobile vector that detection subchondral bone 911 can also be used using the method for norm like this.Tool
For body, to region-of-interest ZT1CSDetect the minimum sending down the fishbone comparison other region Z of norm NormT2CS.Then, by calculate by
The region-of-interest Z of combination minimum norm NormT1CSRepresentative position be used as starting point and sending down the fishbone comparison other region ZT2CSGeneration
Epitope puts the vector as terminal, can detect mobile vector.
In addition, in the above description, illustrating and realizing each processing detected for cartilage surface by multiple functional blocks
Example.But it is also possible in advance using above-mentioned cartilage surface detection process as program storage, and take the program with calculating is machine-readable
And perform.
In addition, in the above description, illustrating and being abutted on the surface of the soft tissue 903 of the knee as detected body
Detector 100, the example for making the detector 100 move in a scanning direction.But, by fixed detector 100, with fixture etc.
Make knee bends to become the relative position relation of detector 100 and soft tissue 903, cartilage 901 and subchondral bone 911
Change also may be used.
In addition, the detector of described above is by detection of the oscillator only along a direction spread configuration of scanning direction
Device.But, by the 2 dimensional region in scanning direction and the direction vertical with scanning direction and depth direction at predetermined intervals
The mode of configuration oscillator also may be used.
In addition it is also possible to be the structure for setting an oscillator and making the oscillator move in a scanning direction.Figure 11 is expression
According to the figure of the detection structure for the mechanical scan (mechanic scan) for making oscillator movement.
The detector 100A of mechanical scan type possesses the ripple transceiver with oscillator.Ripple transceiver is even along scanning with energy
The mode of axle movement is set.Ripple transceiver is moved along scan axis.The oscillator of ripple transceiver is to the direction vertical with scanning direction
Send ultrasonic signal.While making the ripple transceiver move in a scanning direction, while being parked in multiple scan positions and sending super
Acoustic signals, by receiving its echo-signal, can be obtained and the identical echo-signal of detector 100 with multiple oscillators.
In addition, in the above description, illustrating the various kinds of cartilage 901 and subchondral bone 911 relative to detector 100
The situation that this position is only moved in a scanning direction, to situation about being moved on both scanning direction and depth direction, also can
It is applicable structure as described above and processing.
In addition, the detection method of the subchondral bone under the 1st state [T1] is not limited to open from depth as shallow side as described above
The method for beginning successively to be compared echo level with threshold value, can also be applied since the deep side of depth successively by echo level and threshold value
The method compared.Specifically, because the attenuation than the echo inside subchondral bone bone more in the inner part is bigger, echo level
About 0.Therefore, the setting more slightly larger than 0 is set the threshold to.Then, by echo since deep side, the i.e. bone private side of depth
Level is compared with threshold value, using echo level be the threshold value more than position echo data as subchondral bone echo data and
Obtain.
In addition, in the above description, illustrating and being applied as the detection method of the subchondral bone under the 2nd state [T2]
Can also be with the subchondral bone echo-signal of the example of the similar degree of the echo-signal of the 1st state [T1], but the 2nd state [T2]
The detection of the subchondral bone echo-signal of 1st state [T1] similarly carries out the comparison with threshold value.
Description of symbols:
10:Soft tissue cartilage boundary face detection means, 11:Operating portion, 12:Transmission control unit, 13:Echo signal reception
Portion, 14:Data analysis portion, 100,100A:Detector, 141:AD transformation components, 142:Storage part, 143:Determination unit, 901:Cartilage,
902:Bone, 903:Soft tissue, 911:Subchondral bone.
Claims (17)
1. the boundary face of a kind of soft tissue cartilage boundary face detection method, detection soft tissue and cartilage, it is characterised in that have:
1st echo-signal receive and dispatch process, under the 1st state into detected body send ultrasonic signal and obtain the 1st echo letter
Number;
2nd echo-signal receives and dispatches process, sends ultrasonic signal into the detected body under the 2nd state and obtains the 2nd time
Ripple signal;
Subchondral bone detects process, using the subchondral bone echo-signal of the state of the 1st detection of the backscatter signal the 1st, uses institute
State the subchondral bone echo-signal of the state of the 2nd detection of the backscatter signal the 2nd;
Mobile vector detects process, according to the subchondral bone of the subchondral bone echo-signal of the 1st state and the 2nd state
Echo-signal, detects the mobile vector of the subchondral bone from the 1st state to the 2nd state;And
Boundary face detects process, based on sample of the mobile vector to the 1st state and the echo-signal of the 2nd state
Position is corrected, and detects the boundary face of the soft tissue and the cartilage.
2. soft tissue cartilage boundary face detection method as claimed in claim 1,
In boundary face detection process,
According to the mobile vector to each sample position of the 1st echo-signal and each sample bit of the 2nd echo-signal
Put and aligned,
Based on the echo-signal of the identical sample position after contraposition, the boundary face of the soft tissue and the cartilage is detected.
3. soft tissue cartilage boundary face detection method as claimed in claim 2,
In boundary face detection process,
Detect the similar degree of the 1st echo-signal and the 2nd echo-signal at each sample position after the contraposition, root
The boundary face is detected according to the similar degree.
4. soft tissue cartilage boundary face detection method as claimed in claim 3,
In boundary face detection process,
Calculate the 1st echo-signal comprising the comparison other region including the sample position and the 2nd echo-signal
Coefficient correlation detect the boundary face as the similar degree, and based on the coefficient correlation.
5. soft tissue cartilage boundary face detection method as claimed in claim 3,
In boundary face detection process,
According to including the 1st echo-signal in the comparison other region including the sample position and the 2nd echo-signal
Norm is calculated as the similar degree, and the boundary face is detected based on the norm.
6. the soft tissue cartilage boundary face detection method as described in any one of claim 1 to 5,
In subchondral bone detection process,
Believe using the subchondral bone echo-signal of the 1st detection of the backscatter signal the 1st as the subchondral bone echo of the 1st state
Number, believe using the subchondral bone echo-signal of the 2nd detection of the backscatter signal the 2nd as the subchondral bone echo of the 2nd state
Number;
In mobile vector detection process,
According to the 1st subchondral bone echo-signal and the 2nd subchondral bone echo-signal, detection from the 1st state to
The mobile vector of the subchondral bone of 2nd state.
7. soft tissue cartilage boundary face detection method as claimed in claim 6,
In mobile vector detection process,
Region-of-interest is set to the 1st subchondral bone echo-signal,
Comparison other region is used to the 2nd subchondral bone echo-signal setting sending down the fishbone,
With comparison other region detection it is zone similarity by the high sending down the fishbone of similar degree between the region-of-interest, according to institute
The position relationship for stating region-of-interest and the zone similarity detects the mobile vector.
8. soft tissue cartilage boundary face detection method as claimed in claim 7,
In mobile vector detection process,
The sending down the fishbone region of search bigger than the region-of-interest is set to the 2nd subchondral bone echo-signal, in the sending down the fishbone
With setting the sending down the fishbone comparison other region in region of search.
9. soft tissue cartilage boundary face detection method as claimed in claim 7 or 8,
In mobile vector detection process,
The 1st subchondral bone echo-signal is calculated with the coefficient correlation of the 2nd subchondral bone echo-signal as the class
The zone similarity is detected like degree, and based on the coefficient correlation.
10. soft tissue cartilage boundary face detection method as claimed in claim 7 or 8,
In mobile vector detection process,
The 1st subchondral bone echo-signal is calculated with the norm of the 2nd subchondral bone echo-signal as the similar degree,
And the zone similarity is detected based on the norm.
11. the soft tissue cartilage boundary face detection method as described in any one of claim 1 to 5,
In subchondral bone detection process,
Obtain the signal intensity of the 1st echo-signal successively along depth direction from face side, in the range of prescribed depth, detection
It is less than to the signal intensity after subchondral bone detection threshold value, letter more than the subchondral bone detection threshold value will be detected
The signal detection of the scope of number intensity is the subchondral bone echo-signal of the 1st state,
The subchondral bone echo-signal of the 1st state and the similar degree of the 2nd echo-signal are detected, and by described 2nd time
Among ripple signal is described the with the similar degree highest detection of the backscatter signal of the subchondral bone echo-signal of the 1st state
The subchondral bone echo-signal of 2 states.
12. the soft tissue cartilage boundary face detection method as described in any one of claim 1 to 5,
In subchondral bone detection process,
Obtain the signal intensity of the 1st echo-signal successively along depth direction from deep side, the subchondral bone will be detected
The signal detection of the scope of signal intensity more than detection threshold value is the subchondral bone echo-signal of the 1st state,
The subchondral bone echo-signal of the 1st state and the similar degree of the 2nd echo-signal are detected, and by described 2nd time
Among ripple signal is described the with the similar degree highest detection of the backscatter signal of the subchondral bone echo-signal of the 1st state
The subchondral bone echo-signal of 2 states.
13. the boundary face of a kind of soft tissue cartilage boundary face detection means, detection soft tissue and cartilage, it is characterised in that possess:
Receiving and transmitting part, sends ultrasonic signal into detected body under the 1st state and exports the 1st echo-signal, at described soft group
Knit and institute is sent into the detected body under the 2nd states different relative to the 1st state from the position relationship of the cartilage
State ultrasonic signal and export the 2nd echo-signal;And data analysis portion, use the 1st shape described in the 1st detection of the backscatter signal
The subchondral bone echo-signal of state, uses the subchondral bone echo-signal of the 2nd state described in the 2nd detection of the backscatter signal, root
According to the subchondral bone echo-signal of the 1st state and the subchondral bone detection of the backscatter signal of the 2nd state from the 1st shape
State corrects the 1st echo-signal and described to the mobile vector of the subchondral bone of the 2nd state based on the mobile vector
The sample position of 2nd echo-signal, it is soft according to the 1st echo-signal and the 2nd detection of the backscatter signal after correction
The boundary face of tissue and the cartilage.
14. soft tissue cartilage boundary face detection means as claimed in claim 13,
The receiving and transmitting part possesses the multiple oscillators arranged along scanning direction.
15. soft tissue cartilage boundary face detection means as claimed in claim 13,
The receiving and transmitting part possesses by scanning direction and vertical with the sending direction of the scanning direction and the ultrasonic signal
Direction as defined in region with multiple oscillators of two-dimensional arrangements.
16. soft tissue cartilage boundary face detection means as claimed in claim 13,
The receiving and transmitting part possesses single oscillator and the travel mechanism for making the oscillator move in a scanning direction.
17. the soft tissue cartilage boundary face detection means as described in any one of claim 14 to 16, the scanning direction is
The direction for changing the relative position of the soft tissue and the cartilage.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012286738 | 2012-12-28 | ||
JP2012-286738 | 2012-12-28 | ||
PCT/JP2013/079826 WO2014103512A1 (en) | 2012-12-28 | 2013-11-05 | Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104837410A CN104837410A (en) | 2015-08-12 |
CN104837410B true CN104837410B (en) | 2017-09-22 |
Family
ID=51020611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380063141.1A Active CN104837410B (en) | 2012-12-28 | 2013-11-05 | Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP6266537B2 (en) |
CN (1) | CN104837410B (en) |
WO (1) | WO2014103512A1 (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013066821A2 (en) | 2011-10-28 | 2013-05-10 | Decision Sciences International Corporation | Spread spectrum coded waveforms in ultrasound imaging |
US9844359B2 (en) | 2013-09-13 | 2017-12-19 | Decision Sciences Medical Company, LLC | Coherent spread-spectrum coded waveforms in synthetic aperture image formation |
CN105636520B (en) * | 2013-10-07 | 2018-12-07 | 古野电气株式会社 | Diagnostic ultrasound equipment and characteristic quantity calculating method |
EP3359048B1 (en) * | 2015-10-08 | 2023-07-19 | Decision Sciences Medical Company, LLC | Acoustic orthopedic tracking system and methods |
JP6496042B2 (en) * | 2015-12-04 | 2019-04-03 | 古野電気株式会社 | Ultrasonic analysis apparatus, ultrasonic analysis method, and ultrasonic analysis program |
WO2017110361A1 (en) * | 2015-12-25 | 2017-06-29 | 古野電気株式会社 | Ultrasonic analyzing device, ultrasonic analyzing method and ultrasonic analyzing program |
JP7107522B2 (en) * | 2018-06-07 | 2022-07-27 | 古野電気株式会社 | Ultrasonic analysis device, ultrasonic analysis method and ultrasonic analysis program |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE60136474D1 (en) * | 2000-09-14 | 2008-12-18 | Univ R | ASSESSMENT OF THE CONDITION OF A JOINT AND LOSS OF CARTEL TISSUE |
US7022073B2 (en) * | 2003-04-02 | 2006-04-04 | Siemens Medical Solutions Usa, Inc. | Border detection for medical imaging |
JP4751282B2 (en) * | 2006-09-27 | 2011-08-17 | 株式会社日立製作所 | Ultrasonic diagnostic equipment |
US7862510B2 (en) * | 2007-02-09 | 2011-01-04 | Cyberlogic, Inc. | Ultrasonic bone assessment apparatus and method |
JP5192921B2 (en) * | 2008-06-23 | 2013-05-08 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic equipment |
JP2010094422A (en) * | 2008-10-20 | 2010-04-30 | National Institute Of Advanced Industrial Science & Technology | Apparatus, method and program for discriminating tissue boundary, and recording medium having the program recorded thereon |
JP5377166B2 (en) * | 2009-09-01 | 2013-12-25 | 古野電気株式会社 | Ultrasound bone analyzer |
JP5414581B2 (en) * | 2010-03-12 | 2014-02-12 | 株式会社東芝 | Ultrasonic diagnostic equipment |
JP2011125757A (en) * | 2011-03-30 | 2011-06-30 | Hitachi Aloka Medical Ltd | Ultrasonic image data processor |
-
2013
- 2013-11-05 CN CN201380063141.1A patent/CN104837410B/en active Active
- 2013-11-05 WO PCT/JP2013/079826 patent/WO2014103512A1/en active Application Filing
- 2013-11-05 JP JP2014554214A patent/JP6266537B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP6266537B2 (en) | 2018-01-24 |
JPWO2014103512A1 (en) | 2017-01-12 |
CN104837410A (en) | 2015-08-12 |
WO2014103512A1 (en) | 2014-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104837410B (en) | Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means | |
US20210022704A1 (en) | Ultrasound diagnostic system, ultrasound image generation apparatus, and ultrasound image generation method | |
US8784318B1 (en) | Aberration correction using channel data in ultrasound imaging system | |
Shin et al. | Estimation of average speed of sound using deconvolution of medical ultrasound data | |
EP2047801A1 (en) | Ultrasonographic device | |
US20120116220A1 (en) | Spatially-fine shear wave dispersion ultrasound vibrometry sampling | |
EP2905633B1 (en) | Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method | |
WO2007135450A2 (en) | Image data processing systems | |
CN100531674C (en) | Ultrasonographic device | |
JP5692079B2 (en) | Displacement estimation method and displacement estimation apparatus | |
Grau et al. | Adaptive multiscale ultrasound compounding using phase information | |
US20200397409A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
JP5467922B2 (en) | Ultrasonic diagnostic equipment | |
CN109788944A (en) | The control method of diagnostic ultrasound equipment and diagnostic ultrasound equipment | |
CN112584773A (en) | System and method for performing pulse wave velocity measurements | |
JP6190374B2 (en) | Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program | |
US20190254639A1 (en) | Quantitative Shear Wave Elasticity Imaging Method and System | |
US20040249283A1 (en) | Method and apparatus for breast imaging utilizing ultrasound | |
JP4768315B2 (en) | Ultrasonic signal processing apparatus and ultrasonic signal processing method | |
GB2436655A (en) | Deformation determination in pulse echo or ultrasonic imaging | |
JP6496045B2 (en) | Ultrasonic analysis apparatus, ultrasonic analysis method, and ultrasonic analysis program | |
US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
JP4632517B2 (en) | Transmission / reception wave phasing method and ultrasonic diagnostic apparatus | |
CN108024795A (en) | Frequency spectrum resolver, frequency spectrum analytic method and ultrasonic imaging apparatus | |
WO2013161228A1 (en) | Ultrasonic diagnostic device and ultrasonic diagnostic device control method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |