CN104837410A - Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program - Google Patents
Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program Download PDFInfo
- Publication number
- CN104837410A CN104837410A CN201380063141.1A CN201380063141A CN104837410A CN 104837410 A CN104837410 A CN 104837410A CN 201380063141 A CN201380063141 A CN 201380063141A CN 104837410 A CN104837410 A CN 104837410A
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Rheumatology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The objective of the invention is to accurate detect boundaries between soft tissue and cartilage non-invasively. Echo data for first and second states are acquired (S101 and S102). On the basis of the echo data for the first and second states, echo data for a subchondral bone (911) in the first state and the second state are detected (S103). On the basis of the subchondral bone echo data for the first state and the second state, movement vectors for the subchondral bone echo data are detected (S104). On the basis of the movement vectors, alignment is performed between sample data for comparison of the echo data for the first state and the echo data for the second state (S105). Correlation coefficients for the echo data for the first state and the echo data for the second state, which are to be identical sample positions for comparison by way of alignment, are calculated (S106). Determinations are made of areas of high coefficients as cartilage and areas of low coefficients as soft tissue, and interfaces therebetween are detected (S107).
Description
Technical field
The present invention relates to a kind of with the soft tissue cartilage boundary detection method detected from the boundary face of ultrasound wave to cartilage and soft tissue of outside.
Background technology
In the past, various generation is proposed for diagnosing the device of the information of cartilage state.Such as, in the diagnostic ultrasound equipment of patent documentation 1, the hyperacoustic detector of transmitting-receiving is connected to the surface of knee, with the echo-signal from knee inside obtained with this detector, the state of diagnosis cartilage.That is, the diagnostic ultrasound equipment of patent documentation 1 is diagnosed with the state of noninvasive mode to cartilage.Further, in the diagnostic ultrasound equipment of patent documentation 1, cartilage is detected by the difference of the level (intensity) of the echo-signal of depth direction.
At first technical literature
Patent documentation
Patent documentation 1: JP 2010-305 publication
Summary of the invention
The problem that invention will solve
But in the device and method of patent documentation 1, the difference according to the echo-signal level of soft tissue (muscle or skin) and the echo-signal level of cartilage detects cartilage surface.Therefore, if the level of echo-signal does not have difference between soft tissue and cartilage, then cartilage surface can not be detected exactly.
Further, use ultrasonic signal in the past substantially, in the boundary face of soft tissue and cartilage surface not sharply and change exactly, but echo-signal level does not have large difference to echo-signal level in this boundary face.Thus, by method in the past, cartilage surface can not be detected exactly.
The object of the present invention is to provide a kind of soft tissue cartilage boundary face detection method that can detect the boundary face of soft tissue and cartilage surface in noninvasive mode exactly.
The means of dealing with problems
The present invention is the invention of the soft tissue cartilage boundary face detection method relating to the boundary face detecting soft tissue and cartilage and has following characteristics.Soft tissue cartilage boundary face detection method has the 1st echo-signal transmitting-receiving operation, the 2nd echo-signal transmitting-receiving operation, subchondral bone detects operation, mobile vector detects operation and boundary face detects operation.
1st echo-signal transmitting-receiving operation sends ultrasonic signal and obtains the 1st echo-signal under the 1st state in detected body.2nd echo-signal transmitting-receiving operation sends ultrasonic signal and obtains the 2nd echo-signal under the 2nd state in detected body.
Subchondral bone detects the subchondral bone echo-signal that operation detects the 1st state from the 1st echo-signal, detects the subchondral bone echo-signal of the 2nd state from the 2nd echo-signal.
The operation that detects mobile vector detects the mobile vector from the 1st state to the subchondral bone of the 2nd state from the subchondral bone echo-signal of the 1st state and the subchondral bone echo-signal of the 2nd state.
Boundary face detects operation and corrects based on the sample position of mobile vector to the echo-signal of the 1st state and the 2nd state, thus detects the boundary face of soft tissue and cartilage.
In the method, make use of cartilage and be attached to subchondral bone, soft tissue is non-cohesive on cartilage, the situation that soft tissue laterally can slide at cartilage surface.
If make the detector of transmission ultrasonic signal contact to move in detected body, then soft tissue is followed, and cartilage and subchondral bone are not followed.Therefore, subchondral bone is consistent relative to the change in location of detector with cartilage relative to the change in location of detector, is equivalent to mobile vector, but their change in location and soft tissue inconsistent relative to the change in location of detector.In addition, even if under the state making detector contact detected body, make soft tissue and cartilage and have the detected body lateral bend of cartilage, similarly, subchondral bone is consistent relative to the change in location of detector with cartilage relative to the change in location of detector, subchondral bone relative to the change in location of detector and soft tissue inconsistent relative to the change in location of detector.
Therefore, if correct the sample position of echo-signal of the 1st state and the sample position of the echo-signal of the 2nd state according to mobile vector, and the echo-signal of each sample position is compared, then between soft tissue from cartilage (and subchondral bone), comparative result is different.By utilizing this difference, soft tissue and cartilage can be differentiated, also can detect the boundary face of soft tissue and cartilage.
In addition, in the detection method of soft tissue cartilage boundary face of the present invention, subchondral bone detects operation obtains described 1st echo-signal successively signal intensity along depth direction, will detect that the signal of the scope of the signal intensity of more than described subchondral bone detection threshold value detects as the subchondral bone echo-signal of described 1st state.Subchondral bone detects operation and detects the subchondral bone echo-signal of described 1st state and the similar degree of described 2nd echo-signal, and echo-signal the highest for similar degree is detected as the subchondral bone echo-signal of described 2nd state.
In addition, in the detection method of soft tissue cartilage boundary face of the present invention, subchondral bone detects operation obtains described 1st echo-signal successively from side, deep signal intensity along depth direction, will detect that the signal of the scope of the signal intensity of more than described subchondral bone detection threshold value detects as the subchondral bone echo-signal of described 1st state.Subchondral bone detects operation and detects the subchondral bone echo-signal of described 1st state and the similar degree of described 2nd echo-signal, and echo-signal the highest for similar degree is detected as the subchondral bone echo-signal of described 2nd state.
Show the concrete detection method of subchondral bone in these methods.
Invention effect
According to the present invention, ultrasonic signal can be sent from the outside of the detected bodys such as knee, in this echo-signal of external reception of detected body, and detect the boundary face of soft tissue and cartilage surface exactly.Thereby, it is possible to detect the echo come from cartilage exactly, and the diagnosis in cartilage can be effectively utilized.
Accompanying drawing explanation
Fig. 1 is the block diagram of the structure representing the soft tissue cartilage boundary checkout gear 10 relating to embodiments of the present invention.
Fig. 2 represents the figure of the detector 100 of the soft tissue cartilage boundary checkout gear 10 relating to embodiments of the present invention relative to the set-up mode of detected body.
Fig. 3 is the figure of the detection concept for illustration of the cartilage surface relating to embodiments of the present invention.
Fig. 4 relates to the flow chart of the soft tissue cartilage boundary face detection method of embodiments of the present invention.
Fig. 5 is the figure of the waveform example of each echo-signal represented under the 1st state [T1] with the 2nd state [T2].
Fig. 6 is the oscillogram of the detection concept for illustration of the mobile vector relating to present embodiment.
Fig. 7 is the figure of the definition representing the mobile vector relating to present embodiment.
Fig. 8 is the figure of the distribution representing the mobile vector relating to present embodiment.
Fig. 9 is the oscillogram of the method (the 1st method) belonging to cartilage 901 for illustration of detection of echoes data or belong to soft tissue 903.
Figure 10 is the oscillogram of the method (the 2nd method) belonging to cartilage 901 for illustration of detection of echoes data or belong to soft tissue 903.
Figure 11 is the figure of the detection architecture represented based on the mechanical scan making oscillator movement.
Detailed description of the invention
With reference to accompanying drawing, the soft tissue cartilage boundary face detection method and soft tissue cartilage boundary face checkout gear that relate to embodiments of the present invention are described.Fig. 1 represents the block diagram of the structure of the soft tissue cartilage boundary checkout gear 10 relating to embodiments of the present invention.Fig. 2 represents that the detector 100 relating to the soft tissue cartilage boundary checkout gear 10 of present embodiment is 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 situation of the 2nd state (T=T2).In addition, in the following description, the example making detector 100 movement is illustrated, but the situation to mobile detected body, also can be suitable for following method or structure.Such as, as the knee of detected body abutting detector 100 and fixing, make knee bend and stretch such situation and also can be suitable for.That is, as long as the method that the position relationship of soft tissue and cartilage changes between the 1st state [T1] and the 2nd state [T2] and structure just can be suitable for.
Fig. 3 is the figure of the detection concept for illustration of the cartilage surface relating to embodiments of the present invention, Fig. 3 (A) represents the 1st state [T1] (T=T1), and Fig. 3 (B) represents the 2nd state [T2] (T=T2).Fig. 3 is be the figure that smooth plane is observed by the surface replacement of the region being sent out ultrasonic signal and its near zone.
Soft tissue cartilage boundary face checkout gear 10 possesses operating portion 11, sends control part 12, echo signal reception portion 13, data parsing portion 14 and detector 100.Send control part 12, echo signal reception portion 13 and detector 100 to be equivalent to " receiving and transmitting part " of the present invention.
The operation input of operating portion 11 accepted user.Such as, operating portion 11 possesses multiple operator (not shown), according to the operation of user to operator, indicates the process detecting cartilage surface to start to perform to sending control part 12.
Send control part 12 and generate the ultrasonic signal carrier wave be made up of hyperacoustic frequency being shaped as pulse type.Send control part 12 and generate ultrasonic signal respectively under the 1st state [T1] with the 2nd state [T2].
Send control part 12 and export ultrasonic signal to detector 100.Detector 100 possesses the multiple oscillators (with reference to Fig. 3) configured on the direction parallel with transmitting-receiving corrugated.The configuration direction of this oscillator is scanning direction.Each oscillator sends the ultrasonic signal be made up of the transmission field angle specified in detected body.Each oscillator sends ultrasonic signal with predetermined time interval, and receives its reflection echo signal.
As shown in Figure 2, detector 100 is connected to as the mode on the surface of the soft tissue 903 of the knee of detected body by the end face of transmitting-receiving side, corrugated and configures detail.Here, as shown in Figure 3, so-called soft tissue 903 be comprise skin and muscle body in part, be present in than the position of cartilage 901 closer to the face side of detected body.Cartilage 901 is attached on subchondral bone 911, and subchondral bone 911 is the tissue be combined with bone (spongy bone) 902.
While make detector 100 and the surface contact of soft tissue 903 as shown in Fig. 2 (A), as shown in Fig. 2 (B), detector 100 is moved surfacewise.Thus, as shown in Figure 2, soft tissue 903 at the surface sliding of cartilage 901, while follow detector 100 and move.The state of the Fig. 2 (A) before this detector 100 is moved is the 1st state (t=T1), and the state of the Fig. 2 (B) after detector 100 is moved is the 2nd state (t=T2).Now, detector 100 is made to move along the orientation (scanning direction) of oscillator.
Each oscillator sends ultrasonic signal respectively under the 1st state [T1] with the 2nd state [T1] in detected body.Now, each oscillator of detector 100 becomes with the direction vertical relative to the surface of soft tissue 903 mode of central axis direction sending wave beam and sends ultrasonic signal.
Echo-signal after the soft tissue 903 of each oscillator received ultrasonic signal in detected body of detector 100, cartilage 901 and the reflection of subchondral bone 911 place, and export echo signal reception portion 13 to.The 2nd echo-signal group SW [T2] that the 1st echo group SW [T1] that the echo-signal obtained by each oscillator in the 1st state [T1] forms by detector 100 is formed with the echo-signal obtained by each oscillator in the 2nd state [T2] exports echo signal reception portion 13 respectively to.
The processing and amplifying that echo signal reception portion 13 specifies each echo-signal also exports data parsing portion 14 to.Echo signal reception portion 13 individually carries out processing and amplifying to each echo-signal of the 1st echo group SW [T1] and each echo-signal of the 2nd echo group SW [T2] and exports data parsing portion 14 to.
Data parsing portion 14 possesses AD transformation component 141, storage part 142 and detection unit 143.AD transformation component 141 is by making data discrete with predetermined time interval to echo-signal sampling.Echo-signal after this data discrete is echo data.Obtain carrying out the echo data after data sampling at regular intervals in the depth direction thereby, it is possible to scan (sweep) by each distance.That is, the echo data distributed in the 2 dimensional region based on scanning direction and depth direction can be obtained.Below, by the echo data of this echo data group referred to as Two dimensional Distribution.AD transformation component 141 exports each echo data to storage part 142.
Storage part 142 stores each echo data exported by AD transformation component 141.The echo data of the Two dimensional Distribution that storage part 142 obtains under possessing storage the 1st state and the capacity of the echo data of the Two dimensional Distribution obtained in the 2nd state.
Detection unit 143 detects the boundary face of soft tissue and cartilage according to the flow process shown in Fig. 4, aftermentioned concrete processing method.Fig. 4 is the flow chart of the soft tissue cartilage boundary face detection method relating to embodiments of the present invention.
Detection unit 143 reads and the echo data obtaining the Two dimensional Distribution obtained under the 1st state being stored in storage part 142 and the echo data (S101, S102) of Two dimensional Distribution obtained under the 2nd state.
Detection unit 143 detects the echo data (being equivalent to " the subchondral bone echo-signal of the 1st state " of the present invention) of the subchondral bone 911 under the 1st state according to the echo data of the Two dimensional Distribution of the 1st state.Detection unit 143 detects the echo data (being equivalent to " the subchondral bone echo-signal of the 2nd state " of the present invention) (S103) of the subchondral bone 911 under the 2nd state according to the echo data of the Two dimensional Distribution of the 2nd state.
Detection unit 143 detects the mobile vector (S104) of subchondral bone echo data according to the subchondral bone echo data of the subchondral bone echo data of the 1st state and the 2nd state.
Detection unit 143 compares the para-position (S105) of the sample data of object for the echo data of the echo data of the Two dimensional Distribution of the 1st state and the Two dimensional Distribution of the 2nd state based on mobile vector.
Detection unit 143 calculates becomes correlation coefficient between the echo data of the 1st state of identical comparison other sample position and the echo data of the 2nd state by para-position.Specifically, detection unit 143 setting comprises comparison other sample position and the comparison other region be made up of the region of Rack in the depth direction.Detection unit 143 calculates the correlation coefficient (S106) of waveform and the waveform be made up of the echo data of the 2nd state in comparison other region be made up of the echo data of the 1st state in comparison other region.
Detection unit 143 detects the region that the comparison other sample position that correlation coefficient is high is gathered and the region that the low comparison other sample position of correlation coefficient is gathered, and detects the border (S107) in these 2 regions.
As described above, the move mode of cartilage 901 is identical with subchondral bone 911, and the move mode of soft element 903 is different from subchondral bone 911.Thus, after based on the para-position of mobile vector, the correlation coefficient of the echo data of the soft tissue 901 be made up of the move mode identical with subchondral bone 911 increases.On the other hand, after based on the para-position of mobile vector, the correlation coefficient of the echo data of the soft tissue 901 be made up of the move mode different from subchondral bone 911 reduces.Therefore, the boundary face detected in step S107 becomes the boundary face of soft tissue 903 and cartilage 901.Like this, the boundary face of soft tissue 903 and cartilage 901 can be detected by applying above-mentioned process.
In addition, if the surface boundary face of cartilage 901 (soft tissue 903 with) of cartilage 901 detected, then not shown cartilage diagnosis information generation unit can be used in the information of cartilage degeneration diagnosis based on the partial echo data genaration of cartilage 901.Specifically, cartilage diagnosis information generation unit obtain echo data near cartilage surface and subchondral bone 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 time during change, detect the variable quantity etc. caused by the degeneration of cartilage surface and also can.This testing result exports as the information that can be used in cartilage degeneration diagnosis by cartilage diagnosis information generation unit.
Next, the detection method of the boundary face of soft tissue and the cartilage performed with data parsing portion 14 is more specifically described with reference to accompanying drawing.In addition, in order to make explanation simple, the displacement Δ x of the detector 100 (each oscillator) between the 1st state [T1] with the 2nd state [T2] is set to illustrate consistent with the configuration space of oscillator.
First, as the 1st state [T1], such as, using as the knee of detected body with under the bending state of the 1st angle, detector 100 is 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 the state of Fig. 3 (A).
The multiple oscillators be configured on detector 100 are configured at regular intervals respectively on the direction (scanning direction parallel with transmitting-receiving corrugated) parallel with the surface of soft tissue 903.Multiple oscillator sends ultrasonic signal to the direction vertical with scanning direction.If the example of Fig. 3, in detector 100, five oscillators are equally spaced configured by along scanning direction, and as shown in Fig. 3 (A), each oscillator being configured in respective position sends ultrasonic signal to the direction vertical with the surface of soft tissue 903.The ultrasonic signal sent from this each allocation position reflects at each depth location of soft tissue 903, cartilage 901, subchondral bone 911, is received, sample in data parsing portion 14 by each oscillator.The echo data group of echo data SWT11, SWT12, SWT13, SWT14, SWT15 of receiving in the 1st state [T1] becomes the 1st echo group SW [T1].
Next, under making detector 100 keep the state being connected to soft tissue 903, in the direction parallel with the surface of soft tissue 903, detector 100 displacement Δ x measures by the direction parallel with scanning direction.This state is the 2nd state [T2], is also the state of Fig. 3 (B).
Now, soft tissue 903 is followed the movement of detector 100 and moves.Therefore, the relative position relationship of the transmitting-receiving corrugated of detector 100 and each position of the scanning direction of soft tissue 903 moving and change not with detector 100.
On the other hand, because cartilage 901 is fixed on bone 902 by subchondral bone 911, even if there is the movement of detector 100, cartilage 901 does not also move.Therefore, the relative position relationship of each position of transmitting-receiving corrugated and the cartilage 901 of detector 100 and the scanning direction of soft tissue 911 changes along with the movement of detector 100.
After becoming the 2nd state, as shown in Fig. 3 (B), send ultrasonic signal from each oscillator of detector 100 to the direction (direction vertical with transmitting-receiving corrugated (scanning direction)) vertical with the surface of soft tissue 903.The ultrasonic signal sent from this each allocation position reflects at each depth location of soft tissue 903, cartilage 901, subchondral bone 911, is received, sample in data parsing portion 14 by each oscillator.The echo data group of echo data SWT21, SWT22, SWT23, SWT24, SWT25 of receiving in the 2nd state [T2] becomes the 2nd echo group SW [T2].
Like this, before detector 100 moves, obtain the 1st echo group SW [T1] be made up of multiple echo data SWT11, SWT12, SWT13, SWT14, SWT15.Then, after detector 100 moves, obtain the 2nd echo group SW [T2] be made up of multiple echo data SWT21, SWT22, SWT23, SWT24, SWT25.
Fig. 5 is for representing the figure of each echo waveform example under the 1st state [T1] with the 2nd state [T2].In addition, feature of the present invention is represented, if the interval at distance (amount of movement) the Δ x of detector 100 movement and the interval of each oscillator, i.e. scanning position is equal in order in Figure 5 more understandable.In addition, the detection in soft tissue cartilage boundary face under this condition is below described.
Even if i () cartilage 901 and subchondral bone 911 detector 100 move, cartilage 901 and subchondral bone 911 also do not move.Therefore, if detector 100 displacement Δ x measures, the position relationship of each position of the position (each scanning position) of each oscillator of detector 100 and cartilage 901 and subchondral bone 911 departs from amount of movement Δ x and measures along scanning direction.
Now, as as shown in each echo waveform of the 1st state [T1] of Fig. 5 and each echo waveform of the 2nd state [T2], the echo data SWT11 in the cartilage 901 of the 1st echo group SW [T1] and the region of subchondral bone 911 is not consistent with the cartilage 901 of echo data SWT21 of the 2nd echo group SW [T2] and the waveform in the region of subchondral bone 911, but with the cartilage 901 of echo data SWT22 of the 2nd echo group SW [T2] and the waveform in the region of subchondral bone 911 approximately consistent.
Similarly, the region of the cartilage 901 of echo data SWT12 and subchondral bone 911 with in the cartilage 901 of echo data SWT23 and the region of subchondral bone 911, waveform is approximately consistent.The region of the cartilage 901 of echo data SWT13 and subchondral bone 911 with in the cartilage 901 of echo data SWT24 and the region of subchondral bone 911, waveform is approximately consistent.The region of the cartilage 901 of echo data SWT14 and subchondral bone 911 with in the cartilage 901 of echo data SWT25 and the region of subchondral bone 911, waveform is approximately consistent.
Therefore, in cartilage 901 and subchondral bone 911, echo data amount of movement Δ x between the 1st state [T1] and the 2nd state [T2] of each scanning position, approximately consistent under the state that namely deviate from an interval by the configuration space of oscillator at scanning position.That is, in cartilage 901 and subchondral bone 911, detector 100 is changed amount of movement Δ x as the sample position of each echo data of benchmark under the 1st state [T1] and the 2nd state [T2].Like this, cartilage 901 becomes identical move mode with subchondral bone 911.
(ii) soft tissue 903 is described above, and detector 100 is connected on the surface of soft tissue 903, and soft tissue 903 is not fixed on the surface of cartilage 901.Therefore, if detector 100 moves amount of movement Δ x, then soft tissue 903 also follows the movement of detector 100, mobile amount of movement Δ x.
Now, as as shown in each echo waveform of the 1st state [T1] of Fig. 5 and each echo waveform of the 2nd state [T2], the echo data SWT11 in the region of the soft tissue 903 of the 1st echo group SW [T1] is approximately consistent with the waveform of the echo data SWT21 in the region of the soft tissue 903 of the 2nd echo group SW [T2].
Similarly, in the region of the soft tissue 903 of echo data SWT12 and in the region of the soft tissue 903 of echo data SWT22, waveform is approximately consistent.In the region of the soft tissue 903 of echo data SWT13 and in the region of the soft tissue 903 of echo data SWT23, waveform is approximately consistent.In the region of the soft tissue 903 of echo data SWT14 and in the region of the soft tissue 903 of echo data SWT24, waveform is approximately consistent.In the region of the soft tissue 903 of echo data SWT15 and in the region of the soft tissue 903 of echo data SWT25, waveform is approximately consistent.
Therefore, in soft tissue 903, the echo data of each scanning position is in the 1st state [T1] with under the 2nd state [T2], and the position along the scanning direction of corresponding detector 100 is approximately consistent.That is, in soft tissue 903, the sample position of detector 100 as each echo data of benchmark is not changed in the 1st state [T1] and between the 2nd state [T2].Thus, move mode and the cartilage 901 described above of the echo data of soft tissue 903 or subchondral bone 911 like that echo data sample position with detector 100 move and the move mode that changes is different.
Like this, the sample position using detector 100 as benchmark of the echo data that ad-hoc location (identical with in the 2nd state [T2] in the 1st state [the T1]) place in subchondral bone 911 and cartilage 901 is reflected is in the 1st state [T1] and between the 2nd state [T2], move amount of movement Δ x along scanning direction measure.
On the other hand, the sample position using detector 100 as benchmark of the echo data that ad-hoc location (identical with in the 2nd state [T2] in the 1st state [the T1]) place in soft tissue 903 is reflected does not change in the 1st state [T1] and between the 2nd state [T2].
Utilize this characteristic, with shown in the soft tissue cartilage boundary face checkout gear step S104 described above of present embodiment, detect the mobile vector of the echo data of subchondral bone 911.Then, carry out the para-position of the sample position of the sample position of the echo data of the 1st state [T1] and the echo data of the 2nd state [T2] with this mobile vector shown in step S105 described above like that.In other words, the combination becoming the sample position of the sample position of the echo data of the 1st state [T1] of comparison other and the echo data of the 2nd state [T2] is determined.
, detect the echo data of subchondral bone 911 now, first, shown in step S103 described above like that.As shown in Figure 4, the echo data of subchondral bone 911 becomes high echo level (amplitude).Although cartilage 901 increases at cartilage surface place echo level, the region between this cartilage surface and subchondral bone 911, echo level reduces.
Therefore, first, the echo level of the echo data SWT11-SWT15 of the 1st echo group SW [T1] is obtained successively along depth direction.Then, in prescribed depth (thickness according to cartilage suitably sets) scope, detection of echoes level is less than the situation of subchondral bone detection threshold value, thereafter, detects the echo level of more than subchondral bone detection threshold value.Thereby, it is possible to detect after this sample 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 mobile vector of the echo data of the subchondral bone 911 of the 2nd state [T2].Fig. 6 is the oscillogram of the detection concept for illustration of the mobile vector relating to present embodiment.In addition, waveform itself is identical with the oscillogram of Fig. 5.
Region-of-interest Z is set to the echo data group of the subchondral bone 911 obtained in the 1st state [T1]
t1CS.Now, region-of-interest Z
t1CSspecify with the length of depth direction (time orientation).Specifically, such as, as shown in Figure 6, for the 2 dimensional region be made up of the scanning direction determined according to each echo data SWT11-SWT15 and depth direction, the scope of the regulation of the depth direction (time orientation) along 1 echo data is set to region-of-interest Z
t1CS.Next, extraction is included in this region-of-interest Z
t1CSinterior echo data.This is in other words, be equivalent to region-of-interest Z
t1CSthe waveform of the echo data group extracted.
Next, to the echo data SWT21-SWT25 setting search region Z of the 2nd echo group SW [T2]
r1.Region of search Z
r1with relative to detector 100 by region-of-interest Z
t1CSposition as benchmark and than region-of-interest Z
t1CSthe mode comprising the echo data of larger scope determines.
Specifically, according to will relative to detector 100 and region-of-interest Z under the 2nd state [T2]
t1CSidentical position as depth direction and the center of scanning direction and the region expand prescribed limit on both depth direction and scanning direction after set.Such as, as shown in Figure 6, in the depth direction, at the depth location P1 of echo-signal SWT12
cSset region-of-interest Z
t1CScenter time, at the depth location P1 of the echo-signal SWT22 of the 2nd state [T2]
cSsetting search region Z
r1center.Then, scope longer for the length of the depth direction than region-of-interest is set to region of search Z
r1the scope of depth direction.In addition, depth location P1
cSas long as at region of search Z
r1in, not region of search Z
r1center also can.
In addition, in a scanning direction region-of-interest Z is set to echo data SWT12
t1CStime, become echo data SWT22 by the center of scanning direction and comprise echo data SWT21, the mode setting search region Z of SWT22, SWT23
r1the scope of scanning direction.
Next, sending down the fishbone comparison other region Z is set to the echo data group of the subchondral bone 911 obtained in the 2nd state [T2]
t2CS.Now, as shown in Figure 6, by sending down the fishbone comparison other region Z
t2CSbe set to selection 1 echo-signal and on its depth direction the region of specific length.Sending down the fishbone comparison other region Z
t2CSthe length of depth direction be set to and region-of-interest Z
t1CSthe length of depth direction consistent.
Sending down the fishbone comparison other region Z
t2CSthe region of search Z set as mentioned above
r1whole region in be set.By each sending down the fishbone comparison other region Z be set like this
t2CSobtain echo data.
Next, to region-of-interest Z
t1CSecho data and the echo data of sending down the fishbone comparison other region ZT2CS carry out relevant treatment, and calculate correlation coefficient.
Change sending down the fishbone comparison other region Z successively on one side
t2CSposition, at region of search Z
r1whole region carry out region-of-interest Z
t1CSecho data and sending down the fishbone comparison other region Z
t2CSthe calculating of correlation coefficient of echo data.
Towards region of search Z
r1whole region, perform such a for a region-of-interest Z
t1CSthe calculating of correlation coefficient.
Next, by each region-of-interest Z
t1CSdetect the sending down the fishbone comparison other region Z that correlation coefficient is maximum
t2CS.This is equivalent to region the most similar in subchondral bone 911 region, and becomes the echo data of the subchondral bone 911 under the 2nd state [T2].The echo data of the subchondral bone 911 under the 2nd state [T2] is detected like this according to the similarity of echo-signal.
Next, the mobile vector (moving direction, amount of movement) to the echo data of the subchondral bone 911 of the 2nd state [T2] from the 1st state [T1] is detected.
Here, so-called mobile vector can define as follows.Fig. 7 represents the figure relating to the definition of the mobile vector of present embodiment.In addition, the situation utilizing region-of-interest and sending down the fishbone comparison other region is represented in the figure 7.Mobile vector v
mdefined by the position of the representative point of the representative point and sending down the fishbone comparison other region that are judged as mutual most similar region-of-interest.Mobile vector v
mfor using the representative point of region-of-interest as starting point, using sending down the fishbone with the representative point in comparison other region as the vector of terminal, define according to moving direction and amount of movement.
Specifically, such as, in the example of fig. 7, the representative point of region-of-interest is, on the P1 of position, scanning direction, be present in the depth location of regulation on echo-signal SWT11.In addition, sending down the fishbone is, on the P2 of position, scanning direction, be present in the depth location of regulation with the representative point in comparison other region on echo-signal SWT22.The depth location of the representative point of region-of-interest is identical with the depth location of the representative point in sending down the fishbone comparison other region.
Now, mobile vector v
mbecome and set scanning direction as moving direction and the interval delta x of scanning position P1, P2 is the vector of amount of movement.
Calculating of such a mobile vector both can be carried out with whole echo data of subchondral bone 911, also can carry out with one or more echo datas of representative.When carrying out with whole echo datas or multiple echo data, using the meansigma methods of mobile vector that calculates the respectively mobile vector as subchondral bone 911.
When calculating mobile vector like this, the mobile vector of each echo data of cartilage 901 and soft tissue 903 is the distribution shown in Fig. 8.Fig. 8 represents the figure relating to the distribution of the mobile vector of present embodiment.As shown in Figure 8, if the mobile vector v of the echo data in subchondral bone 911 region
mfor Δ x, then the mobile vector v of the echo data in cartilage 901 region
malso be Δ x.But the mobile vector vm of the echo data of soft tissue 903 is 0, different from the mobile vector of subchondral bone 911 and cartilage 901.
Utilize this characteristic, detect that each echo data belongs to cartilage 901 or belongs to soft tissue 903 by the method next.Fig. 9 is the oscillogram of the method (the 1st method) belonging to cartilage 901 for illustration of detection of echoes data or belong to soft tissue 903.
1st comparison other region Z is set to the echo data group obtained in the 1st state [T1]
t1n(n is the integer determined by the sample position number along depth direction).Now, the 1st comparison other region Z
t1nspecify with the length of depth direction (time orientation).As concrete example, such as, as shown in Figure 9, to the 2 dimensional region be made up of the scanning direction determined according to each echo data SWT11-SWT15 and depth direction, be the 1st comparison other region Z by the range set of the regulation of the depth direction (time orientation) along 1 echo data
t11.Next, extraction is included in the 1st comparison other region Z
t11interior echo data.This is in other words, be equivalent to the 1st comparison other region Z
t11the waveform of the echo data group extracted.In addition, the region that both can remove subchondral bone 911 sets the 1st comparison other region Z
t11, the region closing on the low echo level of subchondral bone 911 that also can remove in cartilage 901 sets.
Next, the 2nd comparison other region Z is set to the echo data group obtained in the 2nd state [T2]
t2n(n is the integer determined by the sample position number along depth direction).Now, the 2nd comparison other region Z
t2nbe set at and make the 1st comparison other region Z
t1non position after moving by mobile vector.
More specifically, the 1st comparison other region Z is obtained by with under the 1st state [T1]
t1nthe identical oscillator of the oscillator of echo data, detect under the 2nd state [T2] with the 1st comparison other region Z
t1nidentical depth areas.Then, the depth areas under the 2nd state [T2] obtained making this region move by mobile vector is set as the 2nd comparison other region Z
t2n.That is, the sample position being set to the sample position of the echo data of the 1st state [T1] of comparison other and the echo data of the 2nd state [T2] carries out para-position according to mobile vector.
As concrete example, when obtaining the mobile vector Δ x of only movement on scanning direction as above, as shown in Figure 9, if the 1st comparison other region Z
t11be set at the prescribed depth position P11 of the echo data SWT12 of the 1st state [T1], then by the 2nd comparison other region Z
t21be set in the prescribed depth position P21 moving the echo data SWT23 after mobile vector Δ x amount from the echo data SWT22 obtained under the 2nd state [T2] by identical oscillator.
Similarly, as shown in Figure 9, if the 1st comparison other region Z
t12be set on the prescribed depth position P12 of the echo data SWT12 of the 1st state [T1], then the 2nd comparison other region Z
t22be set at and move on the prescribed depth position P22 of the echo data SWT23 after mobile vector Δ x amount from the echo data SWT22 obtained under the 2nd state [T2] by identical oscillator.
Next, the 1st comparison other region Z is calculated
t1necho data and the 2nd comparison other region Z
t2nthe correlation coefficient of echo data.
Here, as mentioned above, cartilage 901 moves by the mobile vector identical with subchondral bone 911, and move mode is identical.Therefore, the 1st comparison other region Z after para-position is carried out according to mobile vector
t1necho data and the 2nd comparison other region Z
t2nthe correlation coefficient of echo data be approximately 1.That is, correlation coefficient increases.
On the other hand, because soft tissue 903 does not carry out the movement identical with cartilage 901 and subchondral bone 911, institute is different in a mobile fashion.Therefore, the 1st comparison other region Z after para-position is carried out according to mobile vector
t1necho data and the 2nd comparison other region Z
t2nthe correlation coefficient of echo data close to 0.That is, correlation coefficient reduces.
Such as, as shown in Figure 9, the 1st comparison other region Z of setting in cartilage 901
t12echo data and the 2nd comparison other region Z
t22echo data be identical waveform, correlation coefficient increases.On the other hand, the 1st comparison other region Z of setting in soft tissue 903
t11echo data and the 2nd comparison other region Z
t21echo data be different waveforms, correlation coefficient reduce.
Like this, can be judged as that the echo data in the comparison other region that correlation coefficient is high is the echo data in cartilage 901, the echo data in the comparison other region that correlation coefficient is low is the echo data in soft tissue 903.
Therefore, such a 1st comparison other region Z is calculated by each echo data beyond the region to subchondral bone 911
t11echo data and the 2nd comparison other region Z
t21correlation coefficient, can detect that each echo data is in soft tissue 903 or in cartilage 901.Then, by detecting the echo data group be judged as in soft tissue 903 and the border being judged as the echo data group in cartilage 901, the boundary face of soft tissue 903 and cartilage 901 can be detected.
As previously discussed, by applying structure and the process of present embodiment, the boundary face of soft tissue and cartilage can be detected exactly in noninvasive mode.In addition, as noninvasive method, whole echo datas is set as region-of-interest, searches for and detect and comparison other region that between respective region-of-interest, correlation coefficient is high, also can be detected the boundary face of soft tissue and cartilage by the move mode in region-of-interest and comparison other region.But, by applying structure and the process of present embodiment, owing to there is no need to carry out the search process to Zone Full, therefore, it is possible to detect the boundary face of soft tissue and cartilage more at high speed.
In addition, in the above description, illustrate and with mobile vector, the sample position of the echo data obtained under the 2nd state [T2] is moved, and then the example in setting comparison other region, but make the sample position of the echo data obtained under the 1st state [T1] move and then set comparison other region also can based on mobile vector.
Figure 10 is the oscillogram of the method (the 2nd method) belonging to cartilage 901 for illustration of detection of echoes data or belong to soft tissue 903.In the 2nd method, the echo data that will obtain under the 2nd state [T2] is as benchmark.
2nd comparison other region Z is set to the echo data group obtained under the 2nd state [T2]
t2n(n is the integer determined by the sample position number along depth direction).As concrete example, such as, as shown in Figure 10, to the 2 dimensional region be made up of the scanning direction determined according to each echo data SWT21-SWT25 and depth direction, be the 2nd comparison other region Z by the range set of the regulation of the depth direction (time orientation) along 1 echo data
t21.Next, extraction is included in the 2nd comparison other region Z
t21interior echo data.This is in other words, be equivalent to the 2nd comparison other region Z
t21the waveform of the echo data group extracted.
Next, the 1st comparison other region Z is set to the echo data group obtained under the 1st state [T1]
t1n(n is the integer determined by the sample position number along depth direction).Now, the 1st comparison other region Z
t1nbe set at and make the 2nd comparison other region Z based on mobile vector
t2nposition after moving on the direction contrary with mobile vector and according to the size of mobile vector.
More specifically, the 1st comparison other region Z is obtained by with under the 2nd state [T2]
t2nthe identical oscillator of the oscillator of echo data, detect under the 1st state [T1] with the 2nd comparison other region Z
t2nidentical depth areas.Then, the depth areas under the 1st state [T1] making this region move to obtain based on mobile vector is set as the 1st comparison other region Z
t1n.That is, relative to the sample position of echo data of the 2nd state [T2] being set to comparison other, para-position is carried out based on the sample position of mobile vector to the echo data of the 1st state [T1].
As concrete example, when obtaining the mobile vector Δ x of only movement on scanning direction as above, as shown in Figure 10, if the 2nd comparison other region Z
t21be set at the prescribed depth position P21 of the echo data SWT22 of the 2nd state [T2], then by the 1st comparison other region Z
t11be set in the prescribed depth position P11 measuring the echo data SWT11 after moving in the opposite direction from the echo data SWT12 obtained under the 1st state [T1] by identical oscillator by mobile vector Δ x.
Similarly, as shown in Figure 10, if the 2nd comparison other region Z
t22be set on the prescribed depth position P22 of the echo data SWT22 of the 2nd state [T2], then by the 1st comparison other region Z
t12be set in and measure on the prescribed depth position P12 of the echo data SWT11 after moving in the opposite direction by mobile vector Δ x from the echo data SWT12 obtained under the 1st state [T1] by identical oscillator.
Next, the 2nd comparison other region Z is calculated
t2necho data and the 1st comparison other region Z
t1nthe correlation coefficient of echo data.
Even if process like this, the boundary face of soft tissue and cartilage also can be detected exactly in noninvasive mode.
In addition, in the above description, illustrate by the 1st comparison other region Z
t1necho data and the 2nd comparison other region
t1nthe correlation coefficient of echo data differentiate the example of soft tissue and cartilage.But by this similar degree, the similar degree of the waveform of the echo-signal that the waveform of echo-signal scanned by each distance that other method detects the 1st state [T1] and each distance of the 2nd state [T2] are scanned, differentiates that soft tissue and cartilage also can.
As concrete example, there is the method for application norm.Norm Norm can according to the formula definition next represented.
[mathematical expression 1]
Here, Ddzt1 (i) is the 1st comparison other region Z
t1nthe echo data values of point of depth direction position i.Ddzt2 (i) is the 2nd comparison other region Z
t2nthe echo data values of point of depth direction position i.K is equivalent to the 1st comparison other region Z
t1nand the 2nd comparison other region Z
t2nresolution, represent be included in the 1st comparison other region Z
t1nand the 2nd comparison other region Z
t2nthe sum of echo data.M is the constant of suitably setting.
Like this, norm Norm is by the 1st comparison other region Z
t1nwith the 2nd comparison other region Z
t2nthe absolute value of residual quantity of echo data of same position get n power, and the value opened n power after this n power value being added in regional integration and obtain.Therefore, the 1st comparison other region Z
t1nwith the 2nd comparison other region Z
t2nthe higher then norm Norm of similar degree less, lower then norm Norm is larger for similar degree.
Thus, the 1st comparison other region Z that norm is little
t1nwith the 2nd comparison other region Z
t2nbe equivalent to the situation that correlation coefficient is high, can be judged as in cartilage 901.Further, the 1st comparison other region Z that norm is large
t1nwith the 2nd comparison other region Z
t2nbe equivalent to the situation that correlation coefficient is low, can be judged as in soft tissue 903.
The method applying norm like this also can use in the situation of the mobile vector detecting subchondral bone 911.Specifically, to region-of-interest Z
t1CSdetect the sending down the fishbone comparison other region Z that norm Norm is minimum
t2CS.Then, by calculating the region-of-interest Z of combination minimum for norm Norm
t1CSrepresentative position as starting point and sending down the fishbone comparison other region Z
t2CSrepresentative position as the vector of terminal, can mobile vector be detected.
In addition, in the above description, the example being realized each process detected for cartilage surface by multiple functional device is illustrated.But, also in advance above-mentioned cartilage surface check processing can be stored as program, and get this program with calculating is machine-readable and performs.
In addition, in the above description, illustrate and abut detector 100 on the surface of the soft tissue 903 of the knee as detected body, make the example of this detector 100 movement in a scanning direction.But, by fixed detector 100, make knee bends that the relative position relation of detector 100 and soft tissue 903, cartilage 901 and subchondral bone 911 is changed with fixture etc. and also can.
In addition, the detector of above-mentioned explanation is only along the detector of a direction spread configuration of scanning direction by oscillator.But, also can by the mode configuring oscillator in the 2 dimensional region of scanning direction and the direction vertical with scanning direction and depth direction at predetermined intervals.
In addition, also can be a structure oscillator being set and making the movement in a scanning direction of this oscillator.Figure 11 is the figure of the detection architecture represented according to the mechanical scan (mechanic scan) making oscillator movement.
The detector 100A of mechanical scan type possesses the ripple transceiver with oscillator.Ripple transceiver is can the even mode along scan axis movement arrange.Ripple transceiver moves along scan axis.The oscillator of ripple transceiver sends ultrasonic signal to the direction vertical with scanning direction.While make this ripple transceiver move in a scanning direction, being parked in multiple scanning position and sending ultrasonic signal, by receiving its echo-signal, the echo-signal identical with the detector 100 with multiple oscillator can be obtained.
In addition, in the above description, illustrate cartilage 901 and subchondral bone 911 situation relative to each sample position only movement in a scanning direction of detector 100, to the situation of movement on scanning direction and depth direction, also can be suitable for structure as described above and process.
In addition, the detection method of the subchondral bone under the 1st state [T1] be not limited to as described above from depth as shallow side successively by the method that echo level compares with threshold value, also can to apply from the dark side of the degree of depth successively by method that echo level compares with threshold value.Specifically, because the attenuation of the echo than subchondral bone bone inside is more in the inner part larger, echo level is approximately 0.Therefore, the setting slightly larger than 0 is set the threshold to.Then, being compared with threshold value by echo level from the dark side of the degree of depth, i.e. bone private side, is that the echo data of the position of more than this threshold value obtains as the echo data of subchondral bone using echo level.
In addition, in the above description, the detection method illustrated as the subchondral bone under the 2nd state [T2] applies the example with the similar degree of the echo-signal of the 1st state [T1], but the subchondral bone echo-signal of the 2nd state [T2] also can be carried out and the comparing of threshold value in the same manner as the detection of the subchondral bone echo-signal of the 1st state [T1].
Description of symbols:
10: soft tissue cartilage boundary face checkout gear, 11: operating portion, 12: send control part, 13: echo signal reception portion, 14: data parsing portion, 100,100A: detector, 141:AD transformation component, 142: storage part, 143: detection unit, 901: cartilage, 902: bone, 903: soft tissue, 911: subchondral bone.
Claims (17)
1. a soft tissue cartilage boundary face detection method, detects the boundary face of soft tissue and cartilage, it is characterized in that having:
1st echo-signal transmitting-receiving operation, sends ultrasonic signal and obtains the 1st echo-signal under described 1st state in described detected body;
2nd echo-signal transmitting-receiving operation, sends ultrasonic signal and obtains the 2nd echo-signal under described 2nd state in described detected body;
Subchondral bone detects operation, uses the subchondral bone echo-signal of described 1st detection of the backscatter signal the 1st state, uses the subchondral bone echo-signal of described 2nd detection of the backscatter signal the 2nd state;
Mobile vector detects operation, according to the subchondral bone echo-signal of described 1st state and the subchondral bone echo-signal of described 2nd state, detects the mobile vector to the subchondral bone of described 2nd state from described 1st state; And
Boundary face detects operation, corrects, detect the boundary face of described soft tissue and described cartilage based on the sample position of described mobile vector to the echo-signal of described 1st state and described 2nd state.
2. soft tissue cartilage boundary face as claimed in claim 1 detection method,
Detect in operation in described boundary face,
Para-position is carried out according to each sample position of described mobile vector to each sample position of described 1st echo-signal and described 2nd echo-signal,
Based on the echo-signal of the same sample position after para-position, detect the boundary face of described soft tissue and described cartilage.
3. soft tissue cartilage boundary face as claimed in claim 2 detection method,
Detect in operation in described boundary face,
Detect described 1st echo-signal at each sample position place after described para-position and the similar degree of described 2nd echo-signal, detect described boundary face according to described similar degree.
4. soft tissue cartilage boundary face as claimed in claim 3 detection method,
Detect in operation in described boundary face,
The correlation coefficient of described 1st echo-signal and described 2nd echo-signal that calculate the comparison other region comprising described sample position as described similar degree, and detects described boundary face based on this correlation coefficient.
5. soft tissue cartilage boundary face as claimed in claim 3 detection method,
Detect in operation in described boundary face,
Described 1st echo-signal in comparison other region and described 2nd echo-signal according to comprising described sample position calculate norm as described similar degree, and detect described boundary face based on this norm.
6. the soft tissue cartilage boundary face detection method as described in any one of claim 1 to 5,
Detect in operation at described mobile vector,
Use described 1st detection of the backscatter signal the 1st subchondral bone echo-signal, use described 2nd detection of the backscatter signal the 2nd subchondral bone echo-signal;
According to described 1st subchondral bone echo-signal and described 2nd subchondral bone echo-signal, detect the mobile vector to the subchondral bone of described 2nd state from described 1st state.
7. soft tissue cartilage boundary face as claimed in claim 6 detection method,
Detect in operation at described mobile vector,
To described 1st subchondral bone echo-signal setting region-of-interest,
To described 2nd subchondral bone echo-signal setting sending down the fishbone comparison other region,
Be zone similarity by described sending down the fishbone comparison other region detection high for similar degree between described region-of-interest, the position relationship according to described region-of-interest and described zone similarity detects described mobile vector.
8. soft tissue cartilage boundary face as claimed in claim 7 detection method,
Detect in operation at described mobile vector,
The sending down the fishbone region of search larger than described region-of-interest to described 2nd subchondral bone echo-signal setting, sets described sending down the fishbone comparison other region in this sending down the fishbone region of search.
9. soft tissue cartilage boundary face as claimed in claim 7 or 8 detection method,
Detect in operation at described mobile vector,
Calculate the correlation coefficient of described 1st subchondral bone echo-signal and described 2nd subchondral bone echo-signal as described similar degree, and detect described zone similarity based on this correlation coefficient.
10. soft tissue cartilage boundary face as claimed in claim 7 or 8 detection method,
Detect in operation at described mobile vector,
Calculate the norm of described 1st subchondral bone echo-signal and described 2nd subchondral bone echo-signal as described similar degree, and detect described zone similarity based on this norm.
11. soft tissue cartilage boundary face detection methods as described in any one of claim 1 to 10,
Detect in operation at described subchondral bone,
The signal intensity of described 1st echo-signal is obtained successively along depth direction from face side, within the scope of prescribed depth, after detecting that described signal intensity is less than subchondral bone detection threshold value, to detect that the signal detection of the scope of the signal intensity of more than described subchondral bone detection threshold value is the subchondral bone echo-signal of described 1st state
Detect the subchondral bone echo-signal of described 1st state and the similar degree of described 2nd echo-signal, and be the subchondral bone echo-signal of described 2nd state by detection of the backscatter signal the highest for similar degree.
12. soft tissue cartilage boundary face detection methods as described in any one of claim 1 to 10,
Detect in operation at described subchondral bone,
Obtain the signal intensity of described 1st echo-signal along depth direction from side, deep successively, will detect that the signal detection of the scope of the signal intensity of more than described subchondral bone detection threshold value is the subchondral bone echo-signal of described 1st state,
Detect the subchondral bone echo-signal of described 1st state and the similar degree of described 2nd echo-signal, and be the subchondral bone echo-signal of described 2nd state by detection of the backscatter signal the highest for similar degree.
13. 1 kinds of soft tissue cartilage boundary face checkout gears, detect the boundary face of soft tissue and cartilage, it is characterized in that possessing:
Receiving and transmitting part, under the 1st state, in detected body, send ultrasonic signal and export the 1st echo-signal, under the 2nd state that described soft tissue is different relative to described 1st state from the position relationship of described cartilage, in described detected body, send described ultrasonic signal and export the 2nd echo-signal; And
Data parsing portion, use the subchondral bone echo-signal of the 1st state described in described 1st detection of the backscatter signal, use the subchondral bone echo-signal of the 2nd state described in described 2nd detection of the backscatter signal, according to the subchondral bone echo-signal of described 1st state and the subchondral bone detection of the backscatter signal of described 2nd state from described 1st state to the mobile vector of the subchondral bone of described 2nd state, the sample position of described 1st echo-signal and described 2nd echo-signal is corrected based on described mobile vector, the boundary face of soft tissue and described cartilage according to described 1st echo-signal after correction and described 2nd detection of the backscatter signal.
14. soft tissue cartilage boundary face as claimed in claim 13 checkout gears,
Described receiving and transmitting part possesses the multiple oscillators along scanning direction arrangement.
15. soft tissue cartilage boundary face as claimed in claim 13 checkout gears,
Described receiving and transmitting part possesses in the region specified by described scanning direction and the direction vertical with the sending direction of described ultrasonic signal with this scanning direction with multiple oscillators of two-dimensional arrangements.
16. soft tissue cartilage boundary face as claimed in claim 13 checkout gears,
Described receiving and transmitting part possesses single oscillator and makes the travel mechanism of this oscillator movement in a scanning direction.
17. soft tissue cartilage boundary face checkout gears as described in any one of claim 13 to 16,
Described scanning direction is the direction that the relative position of described soft tissue and described cartilage is changed.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-286738 | 2012-12-28 | ||
JP2012286738 | 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 true CN104837410A (en) | 2015-08-12 |
CN104837410B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108601584A (en) * | 2015-12-04 | 2018-09-28 | 古野电气株式会社 | Ultrasonic wave resolver, ultrasonic wave analytic method and ultrasonic wave analysis program |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8939909B2 (en) | 2011-10-28 | 2015-01-27 | 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 |
JP6038338B2 (en) * | 2013-10-07 | 2016-12-07 | 古野電気株式会社 | Ultrasonic diagnostic apparatus, ultrasonic diagnostic method, and ultrasonic diagnostic program |
KR20180095796A (en) * | 2015-10-08 | 2018-08-28 | 디시전 사이선씨즈 메디컬 컴패니, 엘엘씨 | Acoustic Orthopedic Tracking Systems and Methods |
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 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040197015A1 (en) * | 2003-04-02 | 2004-10-07 | Siemens Medical Solutions Usa, Inc. | Border detection for medical imaging |
JP2008079792A (en) * | 2006-09-27 | 2008-04-10 | Hitachi Ltd | Ultrasonic diagnostic apparatus |
CN101243981A (en) * | 2007-02-09 | 2008-08-20 | 计算机逻辑公司 | Ultrasonic bone assessment apparatus and method |
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 |
US20110054313A1 (en) * | 2009-09-01 | 2011-03-03 | Furuno Electric Co., Ltd. | Ultrasonic bone analyzing device |
JP2011125757A (en) * | 2011-03-30 | 2011-06-30 | Hitachi Aloka Medical Ltd | Ultrasonic image data processor |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7184814B2 (en) * | 1998-09-14 | 2007-02-27 | The Board Of Trustees Of The Leland Stanford Junior University | Assessing the condition of a joint and assessing cartilage loss |
JP5192921B2 (en) * | 2008-06-23 | 2013-05-08 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic equipment |
JP5414581B2 (en) * | 2010-03-12 | 2014-02-12 | 株式会社東芝 | Ultrasonic diagnostic equipment |
-
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
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040197015A1 (en) * | 2003-04-02 | 2004-10-07 | Siemens Medical Solutions Usa, Inc. | Border detection for medical imaging |
JP2008079792A (en) * | 2006-09-27 | 2008-04-10 | Hitachi Ltd | Ultrasonic diagnostic apparatus |
CN101243981A (en) * | 2007-02-09 | 2008-08-20 | 计算机逻辑公司 | Ultrasonic bone assessment apparatus and method |
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 |
US20110054313A1 (en) * | 2009-09-01 | 2011-03-03 | Furuno Electric Co., Ltd. | Ultrasonic bone analyzing device |
JP2011050555A (en) * | 2009-09-01 | 2011-03-17 | Furuno Electric Co Ltd | Ultrasonic bone analyzing device |
JP2011125757A (en) * | 2011-03-30 | 2011-06-30 | Hitachi Aloka Medical Ltd | Ultrasonic image data processor |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108601584A (en) * | 2015-12-04 | 2018-09-28 | 古野电气株式会社 | Ultrasonic wave resolver, ultrasonic wave analytic method and ultrasonic wave analysis program |
CN108601584B (en) * | 2015-12-04 | 2021-05-28 | 古野电气株式会社 | Ultrasonic analysis device, ultrasonic analysis method, and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN104837410B (en) | 2017-09-22 |
WO2014103512A1 (en) | 2014-07-03 |
JPWO2014103512A1 (en) | 2017-01-12 |
JP6266537B2 (en) | 2018-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104837410A (en) | Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program | |
EP1591068B1 (en) | Ultrasonographic device | |
Shin et al. | Estimation of average speed of sound using deconvolution of medical ultrasound data | |
CN101505664B (en) | Ultrasonograph | |
JP5766303B2 (en) | Radar equipment | |
CN103476344B (en) | Subject information acquisition device | |
US20130123630A1 (en) | Adaptive Image Optimization in Induced Wave Ultrasound Imaging | |
CN103269639A (en) | Shear wave velocity estimation using center of mass | |
CN102460356A (en) | Method for determining the locations of one or more impacts or touches on a surface of an object including two or more transducers | |
EP3359047A1 (en) | Systems and methods for ultrasound elastography with continuous transducer vibration | |
US11357480B2 (en) | Quantitative shear wave elasticity imaging method and system | |
CN102695457A (en) | Ultrasonic diagnostic device, and method for measuring intima-media thickness | |
JP5692079B2 (en) | Displacement estimation method and displacement estimation apparatus | |
CN106199582A (en) | The method being applied to the two-frequency CW radar human body locating and tracking of indoor | |
US20160135694A1 (en) | Medical radar method and system | |
JP6190374B2 (en) | Soft tissue cartilage interface detection method, soft tissue cartilage interface detection device, and soft tissue cartilage interface detection program | |
CN112584773A (en) | System and method for performing pulse wave velocity measurements | |
CN102525549B (en) | Ultrasound wave cartilage resolver, ultrasound wave cartilage analytic method | |
CN105748100B (en) | Quasi-static ultrasonic elastograph imaging displacement calculates method and apparatus | |
CN105636520A (en) | Ultrasound diagnosis device, ultrasound diagnosis method, and ultrasound diagnosis program | |
CN109452954A (en) | Ultrasonic imaging method and device | |
JPWO2014045925A1 (en) | Ultrasonic analysis apparatus, ultrasonic analysis method and program | |
JP4815607B2 (en) | Non-destructive measuring device and non-destructive measuring method of reinforcing bar diameter in reinforced concrete structure by electromagnetic wave radar | |
KR101117544B1 (en) | Ultrasound Diagnostic System and Method For Forming Elastic Image | |
US20170065257A1 (en) | Ultrasound diagnostic device and ultrasound image processing 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 |