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 PDF

Info

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
Application number
CN201380063141.1A
Other languages
Chinese (zh)
Other versions
CN104837410A (en
Inventor
喜屋武弥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Furuno Electric Co Ltd
Original Assignee
Furuno Electric Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Furuno Electric Co Ltd filed Critical Furuno Electric Co Ltd
Publication of CN104837410A publication Critical patent/CN104837410A/en
Application granted granted Critical
Publication of CN104837410B publication Critical patent/CN104837410B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0875Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4514Cartilage
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting 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)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Rheumatology (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

Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means
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.
CN201380063141.1A 2012-12-28 2013-11-05 Soft tissue cartilage boundary face detection method, soft tissue cartilage boundary face detection means Active CN104837410B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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