WO2005025425A1 - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- WO2005025425A1 WO2005025425A1 PCT/JP2004/013140 JP2004013140W WO2005025425A1 WO 2005025425 A1 WO2005025425 A1 WO 2005025425A1 JP 2004013140 W JP2004013140 W JP 2004013140W WO 2005025425 A1 WO2005025425 A1 WO 2005025425A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- area
- image
- elastic
- ultrasonic diagnostic
- scoring
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/5206—Two-dimensional coordinated display of distance and direction; B-scan display
- G01S7/52063—Sector scan display
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
Definitions
- the present invention relates to an ultrasonic diagnostic apparatus capable of displaying a diagnostic site in a subject using an ultrasonic wave as an elastic image representing the hardness or softness of a living tissue.
- a conventional general ultrasonic diagnostic apparatus includes an ultrasonic transmission / reception control unit for controlling ultrasonic transmission / reception, an ultrasonic transmission / reception unit for transmitting / receiving ultrasonic waves to / from a subject, and an ultrasonic transmission / reception unit.
- Tomographic scanning means for repeatedly obtaining tomographic image data in a subject including a moving tissue using a reflected echo signal at predetermined intervals, and image display means for displaying time-series tomographic image data obtained by the tomographic scanning means. It was constituted having. Then, the structure of the living tissue inside the subject was displayed, for example, as a B-mode tomographic image.
- Patent Document 1 JP-A-5-317313
- Patent Document 2 JP-A-2000-060853
- the present invention has been made in view of the above-mentioned points, and is useful for diagnosis of a diseased part. It is an object of the present invention to provide an ultrasonic diagnostic apparatus capable of giving a certain index.
- an ultrasonic diagnostic apparatus of the present invention causes an ultrasonic wave to be transmitted from a probe to a subject, receives a reflected echo signal corresponding to the transmission of the ultrasonic wave, and receives the echo signal.
- a displacement measurement unit configured to measure a displacement of the living tissue of the subject based on a projection echo signal; an elasticity image forming unit configured to determine an amount of distortion or an elastic modulus of the displacement power to form an elasticity image; and display the elasticity image.
- a scoring means for specifying a distortion state or an elastic state based on information of the elasticity image output from the elasticity image forming unit.
- the scoring means specifies a distortion state or an elastic state based on the area or shape of the elastic image.
- a tomographic image forming unit configured to form a tomographic image of the reflected echo signal; and the scoring unit determines a distortion state based on a relationship between a specific region on the tomographic image and a region of a hard part of the elastic image.
- the elastic state is specified.
- FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
- FIG. 2 is a diagram showing an example of a method of attaching a pressure measuring unit (pressure sensor) to an ultrasonic probe and measuring a pressure between a head of the probe and a test object.
- a pressure measuring unit pressure sensor
- FIG. 3 is a diagram showing an example of a scoring process executed by the score input means of FIG. 1.
- FIG. 4 is a diagram showing a display example of an image display when a scoring process is performed using a score input unit.
- FIG. 5 is a diagram showing scoring by a stress signal.
- FIG. 6 is a diagram showing scoring based on a change in a region of an elastic image due to compression.
- FIG. 7 is a diagram showing scoring of an elasticity image based on a hue pattern.
- FIG. 8 is a view showing another embodiment of the scoring process executed by the score input means of FIG. 1.
- FIG. 9 is a diagram showing still another embodiment of the scoring process executed by the score input means of FIG. 1. It is.
- FIG. 10 is a flowchart illustrating an example of an automatic scoring process using an area extraction application.
- FIG. 11 is a diagram schematically showing an example of an operation of the automatic scoring process.
- FIG. 12 is a diagram showing a timing of acquiring an elasticity image for scoring.
- FIG. 13 is a diagram showing one embodiment of an RF signal frame data selection unit in FIG. 1.
- FIG. 14 is a diagram showing another example of the RF signal frame data selection unit in FIG. 1.
- FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
- This ultrasonic diagnostic apparatus obtains a tomographic image of a diagnostic site of the subject 100 using ultrasonic waves, and displays an elastic image representing the hardness or softness of a living tissue.
- this ultrasonic diagnostic apparatus uses an ultrasonic probe 101, a transmitting circuit 102, a receiving circuit 103, a phasing addition circuit 104, and a signal processing unit 105.
- An adder 114, a score input means 115 for inputting a score, and a state of distortion from an elasticity image or an elasticity image and a tomographic image! ⁇ is provided with a scoring unit 116 for specifying the elastic state.
- the ultrasonic probe 101, the transmission circuit 102, the reception circuit 103, the phasing addition circuit 104, and the signal processing unit 105 constitute an ultrasonic transmission / reception unit.
- the ultrasonic transmitting and receiving means obtains one tomographic image by scanning the ultrasonic beam in a fixed direction in the body of the subject 100 using the ultrasonic probe 101.
- the ultrasonic probe 101 is formed by arranging a large number of transducers in a strip shape, and transmits and receives ultrasonic waves to the subject 100 by performing beam scanning mechanically or electronically. Although not shown in the figure, a vibrator that generates a supersonic wave and receives a reflected echo is incorporated therein.
- Each transducer In general, a function of converting an input pulse wave or continuous wave transmission signal into an ultrasonic wave and emitting the ultrasonic wave, and receiving an ultrasonic wave reflected from the inside of the subject 100 into a reception signal of an electric signal. It is formed with the function of converting and outputting.
- the transmission circuit 102 generates a transmission pulse for driving the ultrasonic probe 101 to generate an ultrasonic wave, and transmits an ultrasonic wave transmitted by a built-in transmission phasing addition circuit.
- the convergence point of is set to a certain depth.
- the receiving circuit 103 amplifies the reflected echo signal received by the ultrasonic probe 101 with a predetermined gain.
- the number of received signals corresponding to the number of the amplified vibrators are input to the phasing addition circuit 104 as independent received signals.
- the phasing addition circuit 104 inputs the received signal amplified by the receiving circuit 103, controls the phase thereof, and forms an ultrasonic beam at one or a plurality of convergence points.
- the signal processing unit 105 receives the received signal from the phasing addition circuit 104 and performs various signal processing such as gain correction, log compression, detection, contour enhancement, and filter processing.
- the black-and-white scan converter 106 uses the reflected echo signal output from the signal processing unit 105 of the ultrasonic transmission / reception means to generate RF signal frame data in the subject 100 including the moving tissue using the reflected echo signal.
- the RF signal frame data is acquired at a sound wave cycle, and displayed on the image display 107 via the switching adder 114. Therefore, the black-and-white scan converter 106 is a means for controlling the tomographic scanning means and the system for sequentially reading out the RF signal frame data at the cycle of the television system, for example, the reflected echo signal from the signal processing unit 105 is digitally converted. It includes an A / D converter that converts it into a signal, a plurality of frame memories that store the tomographic image data digitized by the A / D converter in time series, and a controller that controls these operations. Be composed.
- the image display 107 displays the time-series tomographic image data obtained by the black-and-white scan converter 106, that is, the B-mode tomographic image, and is output from the black-and-white scan converter 106 via the switching adder 114. It consists of a D / A converter that converts image data into analog signals, and a color TV monitor that inputs analog video signals from the D / A converter and displays them as images.
- the output-side force of the phasing addition circuit 104 is branched, and an RF signal frame data selection unit 108 and a displacement measurement unit 109 are provided.
- a pressure measuring unit 110 is provided in parallel with this.
- a strain and elasticity calculating unit 111 At the subsequent stage of the displacement measuring unit 109 and the pressure measuring unit 110, there are provided a strain and elasticity calculating unit 111, an elasticity data processing unit 112, and a color scan converter 113, and the output of the color scan converter 113 and the monochrome scan converter 106 is provided.
- a switching adder 114 is provided on the side. Then, the outputs of the black-and-white scan converter 106 and the color scan converter 113 are introduced into a scoring unit 116, where a scoring process is performed. The details of this scoring process will be described later.
- the displacement measuring unit 109 performs one-dimensional or two-dimensional correlation processing based on the set of RF signal frame data selected by the RF signal frame data selecting unit 108, and performs displacement or displacement of each point on the tomographic image. It measures the movement vector (direction and magnitude of displacement).
- a method of detecting the movement vector for example, there are a block matching method and a gradient method as described in Patent Document 1.
- the block matching method divides an image into blocks composed of, for example, NXN pixels, searches for the block closest to the block of interest in the current frame in the previous frame, and performs predictive coding based on this. It is.
- the pressure measurement unit 110 measures or estimates the pressure in the body cavity of the diagnostic site of the subject 100.
- This ultrasonic diagnostic apparatus is installed in the probe head 1011 while transmitting and receiving ultrasonic waves under the control of the control unit 200 using the ultrasonic probe 101 provided in the probe head 1011.
- a pressure is applied or reduced by a pressurizer (not shown) to apply a stress distribution to the body cavity of the diagnostic site of the subject 100.
- the pressure applied to the rod-like member is detected.
- the pressure sensor 1012 is attached to the side of the probe head 1011 and the pressure between the probe head 1011 and the subject 100 is measured at an arbitrary time phase. To be sent to In FIG. 2, a pressurizer that pressurizes and depressurizes the probe head 1011 is omitted.
- the strain and elastic modulus calculating unit 111 calculates the strain or the elastic modulus of each point on the tomographic image based on the displacement (displacement) and the pressure output from the displacement measuring unit 109 and the pressure measuring unit 110, respectively. To generate numerical data of strain or elastic modulus (elastic frame data). This is output to the polarity data processing unit 112.
- the calculation of strain performed by the strain and elastic modulus calculation unit 111 does not require pressure data, for example, and is obtained by calculation by spatially differentiating the displacement.
- the calculation of the Young's modulus, one of the elastic moduli is obtained by calculation by dividing the change in pressure by the change in the amount of movement.
- the color scan converter 113 converts the elasticity frame data output from the elasticity data processing unit 112 with the command output from the control unit 200 of the ultrasonic diagnostic apparatus or the elasticity frame data output from the elasticity data processing unit 112.
- a hue information conversion unit for inputting an upper limit value and a lower limit value as a middle gradation selection range, and providing hue information such as red, green, and blue as elastic image data based on the elastic frame data; You.
- This hue information conversion means converts, for example, a region where large distortion is measured in the elasticity frame data output from the elasticity data processing unit 112 into a red code in the elasticity image data, and conversely, measures a small distortion. It operates to convert the set area into a blue code in the elasticity image data.
- the color scan converter 113 may be constituted by the above-described monochrome scan converter 106.
- the brightness may be increased in the elasticity image data in an area where the distortion is measured to be large, and may be decreased in the elasticity image data in the area where the distortion is measured to be small.
- the switch adder 114 is a means for inputting the black-and-white tomographic image data from the black-and-white scan converter 106 and the color elastic image data from the color scan converter 113 and adding or switching both images. The output is switched so that only black-and-white tomographic image data or only color elastic image data is output, or both image data are added and synthesized. Further, for example, as described in Patent Document 2, a monochrome tomographic image and a monochrome elastic image by a color or monochrome scan converter may be simultaneously displayed on a two-screen display. The image data output from the switching adder 114 is output to the image display 107.
- FIG. 3 is a diagram showing a first embodiment of the scoring process performed by the scoring unit 116 of FIG. 1.
- FIG. 3 (A) shows a monochrome tomographic image output from the monochrome scan converter 106.
- FIG. 4 is a diagram illustrating an example of data (B-mode tomographic image).
- Figure 3 (B)-(F) shows the color scan
- FIG. 9 is a diagram showing an example of color elasticity image data output from the data 113, and schematically shows a case where each symptom is different. These elasticity image data actually perform color display according to the hardness.
- a soft part (a part with a large distortion) is displayed as red
- a hard part (a part with a small distortion) is displayed as blue
- the middle between the two is displayed as a continuous color change.
- Fig. 3 it is not possible to display in such colors, so that the state is indicated by shading, and the soft areas are displayed by rough shading, the hard areas are displayed by fineness, and the areas are displayed by shading.
- the central elliptical portion 61 is a disease site such as a mammary tumor, and is called a low echo area.
- the scoring process is performed by operating the elasticity image data score input means 115 in accordance with how the elasticity image data is displayed in the low echo area 61 according to the state. Each sex image is classified into a score of 1 to 5.
- the scores 115 are based on the following criteria. The observer determines whether or not the power satisfies the criteria, and performs a scoring process. Hereinafter, the criteria for each of the scores 115 will be described.
- Score 2 As shown in Fig. 3 (C), cases where hard regions 63 and 64 based on the elasticity image are partially recognized, avoiding the center of the low echo region 61.
- Score 4 As shown in Fig. 3 (E), a case in which the image is recognized evenly distributed up to the boundary (contour) of the hard region 66 force low echo region 61 based on the elasticity image
- the hard region is a certain threshold value of distortion in a distortion image as an elastic image. For example, if a strain of 0% is set as the threshold value Ts, the hard region is a region in which no compression is achieved even by compression from the body surface. (Measurement point cloud). In the case where the pressing force is captured during scanning of the diagnostic region and the elasticity image is used as the elasticity image, the region where the elasticity is measured is larger than a certain elasticity threshold Ty. do it.
- FIG. 4 is a diagram showing a display example of the image display 107 when performing a scoring process using the score input means 115.
- a translucent color image of elastic image data showing the distribution of tissue hardness is superimposed on a B-mode tomographic image as shown in Fig. 3 (B)-(F).
- a column 1072 for inputting a scoring index is provided on the right side of the column, and a display area 1073 for displaying the index of each score is provided.
- the contents of the respective scores as described above are simply displayed, or the figures such as FIGS. 3 (B) to (F) are also displayed.
- the examiner compares the scoring index 1073 with the scoring index human power column 1072 and inputs a score that is considered appropriate by looking at the display screen 1071. These operations are executed by the score input unit 115.
- a B-mode tomographic image may be displayed in parallel in another window as shown in FIG. This makes it possible to improve the visibility during the scoring process.
- the score entered in the scoring index input field 1072 is also linked to the individual patient reporting function provided in the ultrasound system, so that it can be entered simultaneously in the scoring index input field provided in the reporting function. It may be. Further, a model diagram serving as an index of scoring as shown in FIGS. 3 (B)-(F) may be added to the display area 1073.
- the scoring unit 116 automatically scores using image processing by the force scoring unit 116 described in the case where the examiner performs scoring processing while viewing the image displayed on the display screen 1071 ( Classify.
- the scoring unit 116 performs edge detection using well-known edge detection software on the color elastic image data (distorted image) output from the color scan converter 113, and determines the area B of the hard part area 62-67 as a pixel. Calculate by number.
- scoring is performed using only elasticity image data. For example, an elasticity image serving as a reference image is obtained in advance as reference information, and the number of pixels of the elasticity image area is set as an area B, and a score based on the area of the elasticity image is set as follows.
- scoring is performed by applying the area of the newly acquired elasticity image data to the area B set as described above. If these elastic images are connected in series, set the minimum score to 2 or more. If the images are connected by a smooth circular system, set the minimum score to 3 or more. The score may be determined based on the shape of the image.
- FIG. 5 shows a scoring process using the characteristic that a stress signal is output on an elastic image by pressing a subject.
- FIGS. 5A and 5B are views showing a state where the subject 100 is pressed by the probe 101.
- FIG. 5 (b) shows a malignant tumor.
- FIG. 5 (a) shows a benign tumor.
- scoring section 116 a score corresponding to width 123 of stress signal 122 is set in advance. For example, the score is set as shown in FIG. 4 such that the elastic image is more malignant as the width is larger. Then, by comparing the width set corresponding to the score with the width 123 of the obtained stress signal 122 of the elastic image 121, the elastic image 121 is scored.
- scoring may be performed using the area of the force stress signal 122 in which scoring is performed using the width of the stress signal 122. In this way, scoring can be performed using the stress signal of the elasticity image.
- FIG. 6 shows an example in which a scoring process is performed using an elasticity image that changes when a subject is pressed.
- a malignant tumor there is a characteristic that, when the subject is gradually pressed, for example, the area displayed in blue gradually becomes smaller.
- a benign tumor there is a characteristic that when the subject is gradually pressed, the area displayed in blue gradually increases. So this Using the characteristic, the score for the characteristic corresponding to the color change of the region caused by the compression is set in the front scoring unit 116 as shown in FIG. Then, scoring is performed according to the color of the elastic image that changes with the compression. Thus, scoring can be performed based on a change in the area of the elastic image due to the compression.
- FIG. 7 shows an example of performing scoring processing using a cyst pattern.
- the cyst pattern is a pattern in which a portion 141 of a specific region of the displayed elasticity image with a small depth is displayed in blue, and a region 142 with a large depth is displayed in red in a pair. It should be noted that the area displayed by this cyst pattern is empirically known to be a benign tumor. Therefore, the cyst pattern is separately set as a benign score in the scoring unit 116, and if the cyst pattern is recognized on the elasticity image, it is scored as benign. In this manner, scoring can be performed based on the hue pattern of the elastic image.
- the scoring unit 116 detects the low echo area 61 of the diseased part in the black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106 using known edge detection software.
- the area A of the detected area 61 is obtained by the number of pixels or the like.
- the scoring unit 116 performs similar edge detection on the color elasticity image data (distortion image) output from the color scan converter 113, and performs the same edge detection in the low echo area 61 of the diseased part.
- the area B of the hard part area 62-67 to be measured is obtained by the number of pixels or the like.
- the sum of the areas is defined as the area B.
- the scoring unit 116 performs the following scoring processing based on the obtained relationship between the areas A and B.
- the ratio of the areas A and B is calculated as B / A force 3 ⁇ 4. Also, in the supersonic ultrasound diagnostic apparatus, the threshold of the area ratio is set as
- Score 1 threshold Thl (eg 0.1)
- Score 2 threshold Th2 (for example, 0.3)
- Score 3 threshold Th3 (eg 0.7)
- Score 4 threshold Th4 (for example, 1.0)
- a central portion in the low echo area 61 has a certain ratio smaller than the area of the low echo area (for example, Area AX 0.6)
- a region of interest 61a is provided, and for the region of interest 61a, in the same manner as described above, the hardness of the elasticity image, the extent to which the region occupies, and the evaluation are compared to the scoring method described above. This can be used to evaluate the scoring of the spatial distribution of hard regions inside and outside the low echo region, not just the ratio of the occupied area of the hard region to the low echo region. It can be expected that the result will increase the accuracy of the evaluation.
- a plurality of regions of interest 62b are provided inside or outside the boundary with reference to the boundary of the low echo region as well as the central region of interest.
- the degree of occupancy of the hard region in the elasticity image may be evaluated for the plurality of regions of interest, and the area occupied by each of the plurality of regions of interest 62b may be evaluated.
- Edge detection can be achieved by using a known snake or region growing method. If the edge is hidden by shadows, it can be handled by a method such as setting the ROI or interpolating surrounding edge information.
- scoring is performed using the occupancy of the hard region distributed in the low echo region.
- Statistical processing is performed using a group of gender measurement points as a population, and scoring may be performed based on the statistical features. For example, a method using the average value will be described as follows.
- N the total number of elements in the measurement point group in the low echo area
- the scoring unit 116 calculates the area A of the low-coordinate area 61 of the diseased site in the black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106. It is detected by known region extraction software, for example, the region growing method. The area A of the detected area is obtained by the number of pixels or the like.
- the scoring unit 116 displays the color elastic image data (distorted image) output from the color scan converter 113 as a hard region 62-67 within the area of the area A. Count the number of pixels that can be arbitrarily determined and have a predetermined luminance value or more. The number of pixels is defined as area B. The following scoring process is performed based on the area A and the area B (the number of pixels having a predetermined luminance value or more).
- the ratio of the area (number of pixels) A, B is calculated as B / A force 3 ⁇ 4.
- the area ratio threshold is set in advance.
- Score 1 threshold Thl (eg 0.1)
- Score 2 threshold Th2 (for example, 0.3)
- Score 3 threshold Th3 (eg 0.7)
- Score 4 threshold Th4 (for example, 1.0)
- a hard region 65 ⁇ 3 is obtained by a known region extraction application (eg, a region growing method). An area indicated as 67 and having a predetermined luminance or more is extracted. this At this time, the extraction region is also a region to be extracted that includes the edge that does not remain in the low echo region 61 described above.
- the obtained extraction area is defined as C.
- C / A is calculated as Z 'as the ratio of area A to area C. This value of Z '
- scoring can be performed up to a score of 3 without using a region extraction process for elasticity image data.
- FIG. 10 is a flowchart illustrating an example of an automatic scoring process using the region extraction application.
- FIG. U is a diagram schematically showing an example of the operation of the automatic scoring process.
- the region extraction algorithm (Region Growing) will be described as an example.
- step S91 in order to obtain a strain elasticity image, a low echo portion whose region is to be extracted is selected from a B-mode tomographic image which is drawn after being pressed by an arbitrary means.
- This selection processing is performed by designating an arbitrary position in the low echo area 61 as a source point (denoted as SP in the figure) as shown in FIG.
- a low echo area 61 as shown in FIG. 11B is extracted based on the source point SP specified in step S91 and a threshold given by any means.
- the region growing method is a method of extracting an area adjacent to the designated source point SP, where the luminance value difference falls within a preset value (threshold). Therefore, the area of the extracted low echo area 61 is counted by means such as counting the number of pixels. In FIG. 11B, this area is represented as SBW.
- step S93 in the low echo area, an area having a hardness value and a luminance value representing the region that is harder than the luminance value representing the region arbitrarily determined by the user in the drawn elasticity image.
- this hard region is denoted as Sbl.
- step S94 the area Sbl of the hard region is smaller than the area SBWX 0.7 of the low echo region. It is determined whether or not it is. If the result of the determination is that the area Sbl of the hard area is smaller than the area SBWX 0.7 of the low echo area! / (SBWX 0.7> Sbl), the process proceeds to step S95. Then, in step S95, values up to a score of 13 corresponding to each value are returned.
- step S94 Conversely, as a result of the determination in step S94, the area Sbl of the hard region is changed to the area SBW of the low echo region.
- step S96 If X is 0.7 or more (SBWX 0.7 ⁇ Sbl), the process proceeds to step S96 because there is a possibility that the area Sbl of the hard part is not within the low echo area.
- a region (area Sbl ′) is extracted again from the elasticity image data by the region growing method. In this case, extraction is performed up to the area including the area around the low echo area of the B-mode tomographic image.
- the source point SP may be the one specified in the B-mode tomographic image as it is, or may be newly set separately. Let the area of the region extracted here be Sbl '.
- values of scores 4 and 5 corresponding to the ratio (SBWZSbl ') between the area SBW of the low echo area and the area Sbl' are returned. Specifically, if Sbl 'is SBW, a score of 4 will be returned; otherwise, a score of 5 will be returned.
- the ultrasonic Using the interface the examiner inputs the outline of the diseased part on the displayed B-mode tomogram, obtains the area A of the area 61 based on the outline information, and uses the method described above.
- the scoring may be automatically determined in comparison with the area B according to. Note that the above scoring is an example, and it is important to score using an elastic image including a strain image.
- the numerical values used in the above scoring are merely examples, and actual clinical It is needless to say that it is desirable to use various optimum values in accordance with the above.
- the scoring method particularly for the mammary gland region is used as an example. Needless to say, it defines the coring method.
- the scoring method according to the above-described embodiments may be used not only by itself for scoring but also by combining a plurality of methods and performing one scoring.
- the scoring unit 116 outputs black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106 and color elastic image data (distorted image data) output from the color scan converter 113. ), The same scoring process is performed using the output from the signal processing unit 105, the output from the elasticity data processing unit 112, or the output from the distortion and elastic modulus calculation unit 111. You may do it. Further, in the above-described embodiment, when the force distortion in which the pressure measuring unit is shown in FIG. 1 is obtained and a distorted image is displayed, the pressure measuring unit may not be provided.
- FIG. 12 is a diagram showing the timing at which the distortion or elastic modulus is calculated by the distortion and elastic modulus calculator 111 of the above-described embodiment.
- the state is indicated by 0, and the limit of the pressing is indicated by MAX.
- the timing of acquiring the elasticity image is the timing of the light compression state, and the elasticity image in the compressed state of about 3 to 20% of the compression limit is acquired, and the acquired elasticity image is scored.
- the pressure sensor 1012 detects a pressure value at the limit of compression in advance, and stores the detected pressure value in the scoring unit 116 as a pressure maximum value.
- the pressure value of the elasticity image obtained by pressing is set using the score input means 115.
- the pressure value to be set may be arbitrarily determined as, for example, a pressure value of 10% of the maximum pressure value. Then, the subject 100 is pressed, an elasticity image at the time when the set pressure value is reached is obtained, and scoring is performed on the obtained elasticity image.
- An elasticity image in a state where the subject 100 is strongly pressed may be hard or hard to determine the position and softness of the object 100.
- the scoring accuracy can be improved by ringing.
- FIG. 13 is a diagram showing an embodiment of the RF signal frame data selection unit of FIG. It is.
- the RF signal frame data selection unit 108 arbitrarily selects the number of frames going back in the past (the number of frame intervals with the current frame data) as one RF signal frame data serving as a reference for displacement measurement. That is, the RF signal frame data selection unit 108 includes the RF signal frame data output one after another from the phasing addition circuit 104 one after another at the frame rate of the ultrasonic diagnostic apparatus. In the frame memory 1081.
- the RF signal frame data selection unit 108 sets the data currently secured in the frame memory 1081 as the RF signal frame data N.
- the RF signal frame data selection unit 108 temporally stores the past RF signal frame data N-1, N-2, N-3, N, M in accordance with a control command from the control unit 200 of the ultrasonic diagnostic apparatus. , One RF signal frame data is selected, and it is temporarily stored in the RF signal frame data selection circuit 1082 as RF signal frame data X.
- the RF signal frame data selection unit 108 stores the latest RF signal frame data N stored in the frame memory 1081 and the RF signal frame data X stored in the RF signal frame data selection circuit 1082 in parallel. Then, it outputs to the displacement measuring unit 109.
- FIG. 14 is a diagram illustrating another example of the RF signal frame data selection unit in FIG.
- the RF signal frame data selection unit 108 in FIG. 14 secures the RF signal frame data P acquired in a certain past phase P in the frame memory 1081 according to a control instruction from the control unit 200 of the ultrasonic diagnostic apparatus. .
- the RF signal frame data selection circuit 1082 always refers to the RF signal frame data P secured in the frame memory 1081 as past RF signal frame data at an arbitrary time phase without updating. Therefore, the displacement measurement unit 109 indicates whether the currently secured RF signal frame data N and RF signal frame data P A set of RF signal frame data composed of The settings such as whether to adopt the function shown in Fig. 14 and when to acquire the RF signal frame data P when it is adopted can be arbitrarily switched by the user interface of the ultrasonic diagnostic apparatus. You can change it.
- the interval between past and present RF signal frame data N and P that constitute one set of RF signal frame data is limited to adjacent frames, multiple RF signals acquired during a series of pressurizing or depressurizing operation processes
- the amount of pressurization or depressurization given in the time interval between the RF signal frame data that constitutes the frame data set is sufficient for the amount of pressurization or decompression (generally about 1%) suitable for the rendering of elasticity image data. May not be reached.
- the RF signal frame data selection unit as shown in FIGS. 13 and 14
- the frame interval between the past and present RF signal frame data can be made sufficiently large, and the elasticity can be improved.
- An elastic image based on frame data can be appropriately drawn. This is especially true in ultrasonography, where physical limitations on the physique of the subject do not allow the rate of pressurization or decompression during the series of pressurization or decompression operations to be sufficiently high. Useful.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Acoustics & Sound (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/571,578 US8118746B2 (en) | 2003-09-12 | 2004-09-09 | Ultrasonic diagnostic apparatus |
EP04787803.8A EP1665987B1 (en) | 2003-09-12 | 2004-09-09 | Ultrasonograph |
JP2005513893A JP4685633B2 (ja) | 2003-09-12 | 2004-09-09 | 超音波診断装置 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003-322070 | 2003-09-12 | ||
JP2003322070 | 2003-09-12 | ||
JP2004-067983 | 2004-03-10 | ||
JP2004067983 | 2004-03-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2005025425A1 true WO2005025425A1 (ja) | 2005-03-24 |
Family
ID=34315671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2004/013140 WO2005025425A1 (ja) | 2003-09-12 | 2004-09-09 | 超音波診断装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8118746B2 (ja) |
EP (1) | EP1665987B1 (ja) |
JP (1) | JP4685633B2 (ja) |
WO (1) | WO2005025425A1 (ja) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007312958A (ja) * | 2006-05-25 | 2007-12-06 | Hitachi Medical Corp | 超音波診断装置 |
WO2007138881A1 (ja) | 2006-05-25 | 2007-12-06 | Hitachi Medical Corporation | 超音波診断装置 |
EP1938754A1 (en) * | 2005-10-19 | 2008-07-02 | Hitachi Medical Corporation | Ultrasonograph for creating elastic image |
JP2008212522A (ja) * | 2007-03-07 | 2008-09-18 | Hitachi Medical Corp | 超音波診断装置 |
WO2009063691A1 (ja) * | 2007-11-16 | 2009-05-22 | Hitachi Medical Corporation | 超音波撮像システム |
WO2009123208A1 (ja) * | 2008-03-31 | 2009-10-08 | 株式会社 日立メディコ | 超音波診断装置 |
JP2010193944A (ja) * | 2009-02-23 | 2010-09-09 | Konica Minolta Medical & Graphic Inc | 超音波診断装置 |
JP2011200282A (ja) * | 2010-03-24 | 2011-10-13 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置及びその制御プログラム |
CN1975370B (zh) * | 2005-11-22 | 2012-07-04 | 通用电气公司 | 超声弹性模量计算和成像的方法 |
JP2013146625A (ja) * | 2013-05-08 | 2013-08-01 | Hitachi Medical Corp | 超音波診断装置 |
JP2014030774A (ja) * | 2013-11-22 | 2014-02-20 | Super Sonic Imagine | 関心領域の粘弾性の平均値を測定するための方法および装置 |
JP5438002B2 (ja) * | 2008-06-03 | 2014-03-12 | 株式会社日立メディコ | 医用画像処理装置及び医用画像処理方法 |
US9301732B2 (en) | 2008-03-31 | 2016-04-05 | Hitachi Medical Corporation | Ultrasonic diagnostic arrangements selecting parameters for selecting relevant estimation data for classifying an elasticity image |
JP2016214868A (ja) * | 2015-05-21 | 2016-12-22 | 東芝メディカルシステムズ株式会社 | 超音波診断装置及び画像診断装置 |
US10660614B2 (en) | 2013-06-26 | 2020-05-26 | Sony Corporation | Ultrasonic processing apparatus and method |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1864612A4 (en) * | 2005-03-30 | 2009-10-28 | Hitachi Medical Corp | ULTRASOUND DEVICE |
JP4602825B2 (ja) * | 2005-04-18 | 2010-12-22 | オリンパスメディカルシステムズ株式会社 | 画像表示装置 |
EP1747757B8 (en) * | 2005-07-27 | 2012-02-29 | Samsung Medison Co., Ltd. | Ultrasound system for displaying an elastic image |
FR2899336B1 (fr) | 2006-03-29 | 2008-07-04 | Super Sonic Imagine | Procede et dispositif pour l'imagerie d'un milieu viscoelastique |
JP5479353B2 (ja) * | 2008-10-14 | 2014-04-23 | 株式会社日立メディコ | 超音波診断装置 |
KR101286331B1 (ko) * | 2010-06-09 | 2013-07-15 | 삼성메디슨 주식회사 | 탄성 영상을 제공하는 초음파 시스템 및 방법 |
KR101512291B1 (ko) * | 2013-05-06 | 2015-04-15 | 삼성메디슨 주식회사 | 의료 영상 장치 및 의료 영상 제공 방법 |
JP6594245B2 (ja) * | 2016-03-30 | 2019-10-23 | 富士フイルム株式会社 | 硬さ導出装置、医用撮影システム、硬さ導出方法、及び硬さ導出プログラム |
JP6726744B2 (ja) * | 2016-07-19 | 2020-07-22 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法、及び超音波観測装置の作動プログラム |
US20220148172A1 (en) * | 2020-11-10 | 2022-05-12 | Henry Ford Health System | System and method for detecting medical conditions |
CN114926389B (zh) * | 2022-03-30 | 2024-07-09 | 什维新智医疗科技(上海)有限公司 | 一种超声图像的弹性信号提取装置 |
CN114881923A (zh) * | 2022-03-30 | 2022-08-09 | 什维新智医疗科技(上海)有限公司 | 一种基于超声图像弹性信号的结节回声分析装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4201396B2 (ja) | 1998-08-20 | 2008-12-24 | 株式会社日立メディコ | 超音波診断装置 |
US6120445A (en) * | 1998-10-02 | 2000-09-19 | Scimed Life Systems, Inc. | Method and apparatus for adaptive cross-sectional area computation of IVUS objects using their statistical signatures |
WO2004017815A2 (en) * | 2002-08-23 | 2004-03-04 | The Board Of Regents Of University Of Texas System | Novel risk assessment method based upon coronary calcification distribution pattern imaged by computed tomography |
JP3932482B2 (ja) * | 2002-10-18 | 2007-06-20 | 株式会社日立メディコ | 超音波診断装置 |
JP2004246868A (ja) * | 2003-01-22 | 2004-09-02 | Canon Inc | 画像抽出方法及び該方法をコンピュータに実行させるプログラム |
US7257244B2 (en) * | 2003-02-24 | 2007-08-14 | Vanderbilt University | Elastography imaging modalities for characterizing properties of tissue |
CN1805713B (zh) * | 2003-06-13 | 2010-04-28 | 松下电器产业株式会社 | 超声波诊断装置 |
US8123691B2 (en) * | 2003-08-19 | 2012-02-28 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus for fixedly displaying a puncture probe during 2D imaging |
CA2457376C (en) * | 2003-10-14 | 2015-09-15 | The University Of British Columbia | Method for imaging the mechanical properties of tissue |
-
2004
- 2004-09-09 EP EP04787803.8A patent/EP1665987B1/en not_active Expired - Lifetime
- 2004-09-09 US US10/571,578 patent/US8118746B2/en active Active
- 2004-09-09 JP JP2005513893A patent/JP4685633B2/ja not_active Expired - Fee Related
- 2004-09-09 WO PCT/JP2004/013140 patent/WO2005025425A1/ja active Application Filing
Non-Patent Citations (8)
Title |
---|
HIGASHINO E. ET AL.: "Nyusen shikkan ni taisuru choonpa class bunrui", JOURNAL OF MEDICAL ULTRASONICS, vol. 30, no. SPECIAL EXTRA ISSUE, 15 April 2003 (2003-04-15), pages 125, XP002985734 * |
ITO A. ET AL.: "Nyusen shikkan ni okeru elastography no risho oyo", DAI 11 KAI NIPPON NYUSEN KOJOSEN CHOONPA SHINDAN KAIGI SHOROKUSHU, 11 October 2003 (2003-10-11), pages 50 - 51, XP002985736 * |
KANAI H. ET AL.: "Kekkanheki katasa bunpu no keihiteki hyokaho to denshi senshoku", RINSHOKU KENSA, vol. 51, no. 8, 25 August 2003 (2003-08-25), pages 805 - 812, XP002985733 * |
See also references of EP1665987A4 * |
T. SHIINA ET AL.: "Real time tissue elasticity imaging using the combined autocorrelation method", J. MED. ULTRASONICS, JAPAN, SOCIETY OF ULTRASONICS IN MEDICINE, TOKYO, JP, vol. 29, 21 September 2002 (2002-09-21), pages 119 - 128 |
W.-H. YEH ET AL.: "Elastic modulus measurements of human liver and correlation with pathology", ULTRASOUND IN MEDICINE AND BIOLOGY, NEW YORK, USA, vol. 28, no. 4, 1 April 2002 (2002-04-01), pages 467 - 474, XP004361097, DOI: doi:10.1016/S0301-5629(02)00489-1 |
YAMAKAWA M. ET AL.: "Nyusen shikkan no shoshiki dansei eizoho ni yoru chushutsu", JOURNAL OF MEDICAL ULTRASONICS, vol. 30, no. SPECIAL EXTRA ISSUE, 15 April 2003 (2003-04-15), pages 468, XP002985735 * |
YAMAKAWA M. ET AL.: "Tissue elasticity reconstruction based on 3-dimensional finite-element model", JPN. J. APPL. PHYS., vol. 38, no. 5B, 30 May 1999 (1999-05-30), pages 3393 - 3398, XP002973311 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1938754A1 (en) * | 2005-10-19 | 2008-07-02 | Hitachi Medical Corporation | Ultrasonograph for creating elastic image |
JP5303147B2 (ja) * | 2005-10-19 | 2013-10-02 | 株式会社日立メディコ | 弾性画像を生成する超音波診断装置 |
EP1938754A4 (en) * | 2005-10-19 | 2011-01-05 | Hitachi Medical Corp | ECHOGRAPH FOR CREATING AN ELASTIC IMAGE |
CN1975370B (zh) * | 2005-11-22 | 2012-07-04 | 通用电气公司 | 超声弹性模量计算和成像的方法 |
WO2007138881A1 (ja) | 2006-05-25 | 2007-12-06 | Hitachi Medical Corporation | 超音波診断装置 |
JP2007312958A (ja) * | 2006-05-25 | 2007-12-06 | Hitachi Medical Corp | 超音波診断装置 |
JP2008212522A (ja) * | 2007-03-07 | 2008-09-18 | Hitachi Medical Corp | 超音波診断装置 |
JP4903271B2 (ja) * | 2007-11-16 | 2012-03-28 | 株式会社日立メディコ | 超音波撮像システム |
WO2009063691A1 (ja) * | 2007-11-16 | 2009-05-22 | Hitachi Medical Corporation | 超音波撮像システム |
CN103349557A (zh) * | 2008-03-31 | 2013-10-16 | 株式会社日立医疗器械 | 超声波诊断装置 |
US20110026800A1 (en) * | 2008-03-31 | 2011-02-03 | Akiko Tonomura | Ultrasonic diagnostic apparatus. |
JP2009240464A (ja) * | 2008-03-31 | 2009-10-22 | Hitachi Medical Corp | 超音波診断装置 |
WO2009123208A1 (ja) * | 2008-03-31 | 2009-10-08 | 株式会社 日立メディコ | 超音波診断装置 |
US8718339B2 (en) | 2008-03-31 | 2014-05-06 | Hitachi Medical Corporation | Ultrasonic diagnostic arrangements selecting parameters for selecting relevant estimation data for classifying an elasticity image |
US9301732B2 (en) | 2008-03-31 | 2016-04-05 | Hitachi Medical Corporation | Ultrasonic diagnostic arrangements selecting parameters for selecting relevant estimation data for classifying an elasticity image |
JP5438002B2 (ja) * | 2008-06-03 | 2014-03-12 | 株式会社日立メディコ | 医用画像処理装置及び医用画像処理方法 |
JP2010193944A (ja) * | 2009-02-23 | 2010-09-09 | Konica Minolta Medical & Graphic Inc | 超音波診断装置 |
JP2011200282A (ja) * | 2010-03-24 | 2011-10-13 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置及びその制御プログラム |
JP2013146625A (ja) * | 2013-05-08 | 2013-08-01 | Hitachi Medical Corp | 超音波診断装置 |
US10660614B2 (en) | 2013-06-26 | 2020-05-26 | Sony Corporation | Ultrasonic processing apparatus and method |
JP2014030774A (ja) * | 2013-11-22 | 2014-02-20 | Super Sonic Imagine | 関心領域の粘弾性の平均値を測定するための方法および装置 |
JP2016214868A (ja) * | 2015-05-21 | 2016-12-22 | 東芝メディカルシステムズ株式会社 | 超音波診断装置及び画像診断装置 |
Also Published As
Publication number | Publication date |
---|---|
EP1665987A1 (en) | 2006-06-07 |
US8118746B2 (en) | 2012-02-21 |
US20070112267A1 (en) | 2007-05-17 |
EP1665987B1 (en) | 2016-11-09 |
JPWO2005025425A1 (ja) | 2007-11-08 |
EP1665987A4 (en) | 2009-12-09 |
JP4685633B2 (ja) | 2011-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4685633B2 (ja) | 超音波診断装置 | |
JP4455003B2 (ja) | 超音波診断装置 | |
JP4314035B2 (ja) | 超音波診断装置 | |
JP5264097B2 (ja) | 超音波診断装置 | |
JP4657106B2 (ja) | 超音波診断装置 | |
JP5479353B2 (ja) | 超音波診断装置 | |
EP2263545B1 (en) | Ultrasonograph | |
EP1554982A1 (en) | Ultrasonographic device | |
US8948485B2 (en) | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, ultrasonic image processing program, and ultrasonic image generation method | |
US20110098563A1 (en) | Ultrasonic diagnostic apparatus, ultrasonic image display method and ultrasonic diagnostic program | |
JP2007105400A (ja) | 超音波診断装置及び画像処理装置 | |
JP5016911B2 (ja) | 超音波診断装置 | |
KR101629541B1 (ko) | 초음파 진단 장치 및 그 제어 프로그램 | |
JP5473527B2 (ja) | 超音波診断装置 | |
JP5623609B2 (ja) | 超音波診断装置 | |
JP5802790B2 (ja) | 超音波診断装置 | |
JP4368185B2 (ja) | 超音波診断装置 | |
JP5680703B2 (ja) | 超音波診断装置 | |
JP5128149B2 (ja) | 超音波診断装置 | |
JP4601413B2 (ja) | 超音波診断装置 | |
JP6144990B2 (ja) | 超音波画像撮像装置及び超音波画像撮像方法 | |
JP2012228538A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200480026263.4 Country of ref document: CN |
|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BW BY BZ CA CH CN CO CR CU CZ DK DM DZ EC EE EG ES FI GB GD GE GM HR HU ID IL IN IS JP KE KG KP KZ LC LK LR LS LT LU LV MA MD MK MN MW MX MZ NA NI NO NZ PG PH PL PT RO RU SC SD SE SG SK SY TJ TM TN TR TT TZ UA UG US UZ VN YU ZA ZM |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): BW GH GM KE LS MW MZ NA SD SZ TZ UG ZM ZW AM AZ BY KG MD RU TJ TM AT BE BG CH CY DE DK EE ES FI FR GB GR HU IE IT MC NL PL PT RO SE SI SK TR BF CF CG CI CM GA GN GQ GW ML MR SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2005513893 Country of ref document: JP |
|
REEP | Request for entry into the european phase |
Ref document number: 2004787803 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2004787803 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 2004787803 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2007112267 Country of ref document: US Ref document number: 10571578 Country of ref document: US |
|
WWP | Wipo information: published in national office |
Ref document number: 10571578 Country of ref document: US |