WO2013115194A1 - 医用画像診断装置及びその関心領域設定方法 - Google Patents

医用画像診断装置及びその関心領域設定方法 Download PDF

Info

Publication number
WO2013115194A1
WO2013115194A1 PCT/JP2013/051926 JP2013051926W WO2013115194A1 WO 2013115194 A1 WO2013115194 A1 WO 2013115194A1 JP 2013051926 W JP2013051926 W JP 2013051926W WO 2013115194 A1 WO2013115194 A1 WO 2013115194A1
Authority
WO
WIPO (PCT)
Prior art keywords
region
interest
search range
image data
compression
Prior art date
Application number
PCT/JP2013/051926
Other languages
English (en)
French (fr)
Inventor
智章 長野
Original Assignee
日立アロカメディカル株式会社
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 日立アロカメディカル株式会社 filed Critical 日立アロカメディカル株式会社
Priority to CN201380007856.5A priority Critical patent/CN104093363B/zh
Priority to EP13743588.9A priority patent/EP2810602A4/en
Priority to JP2013556419A priority patent/JP6023091B2/ja
Priority to US14/373,837 priority patent/US9514531B2/en
Publication of WO2013115194A1 publication Critical patent/WO2013115194A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • 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/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00—Image coding
    • 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/486—Diagnostic techniques involving arbitrary m-mode
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/56—Details of data transmission or power supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Definitions

  • the present invention relates to a medical image diagnostic apparatus, and more particularly, to a medical image diagnostic apparatus for setting a region of interest including a characteristic part or a measurement position and a method for setting the region of interest.
  • Physical quantity information of the living tissue of the subject is measured by the medical image diagnostic apparatus is an important information for diagnosing the extent of the lesion.
  • the physical quantity information is the length, area, volume, movement, etc. of the living tissue of the subject.
  • an examiner uses an input device (for example, a mouse or a trackball) of a medical image diagnostic apparatus to select a desired position or region (region of interest) to be measured on a medical image displayed on the display unit screen. Manually set the area.
  • an input device for example, a mouse or a trackball
  • the manual setting of the region of interest by all the examiners does not always set the same region of interest under the same conditions.
  • the medical image diagnostic apparatus searches for the outline of the organ region with respect to the entire image data of the medical image, sets the region of interest according to template matching or the like based on the search result, and determines the position of the annulus of the heart
  • Patent Document 1 proposes a method for extracting the contour of the left ventricle based on the above.
  • an object of the present invention is to provide a medical image diagnostic apparatus and a method for setting a region of interest that achieve a reduction in search time for setting the region of interest when the region of interest is set in real time.
  • the present invention captures a medical image of a subject by an image capturing unit, sets image data of the medical image of the subject captured by an image data compression unit as uncompressed image data, and sets a plurality of image data before compression
  • the compressed image data is generated by compression based on the pixels, and the search range of the image data after compression is set by the search range setting unit, and the search range of the image data before compression is set, and the region of interest setting unit setting a region of interest in the medical image based on the search range of the image data after the compression and the search range of the image data before the compression by.
  • the search time for setting the region of interest is shortened.
  • a medical image diagnostic apparatus and a region of interest setting method thereof can be provided.
  • FIG. 1 is a block diagram illustrating an ultrasonic diagnostic apparatus (medical image diagnostic apparatus) according to a first embodiment of the invention.
  • 3 is a flowchart for explaining the operation of the medical image diagnostic apparatus according to the first embodiment of the present invention.
  • the figure which shows that a search range setting part sets a search range in the compressed image data, changing the compression rate (reduction rate) of image data.
  • Image reduction unit the image data compression section
  • FIG. The figure which shows that a template starts from a center point and moves spirally.
  • FIG. 4 is a diagram showing a relationship between a reduction rate of image data and the cumulative search pixel number according to the first embodiment of the present invention. Shows that the coordinate difference based on the plurality of the search range is calculated, the search range as coordinate difference is the smallest is set. The figure which showed setting the other annulus part only in the search range of the left side or the right side rather than the position of the annulus part when one annulus part was set. The figure which showed that a region of interest was set on B mode image data, when using an IMT measurement application. 10 is a flowchart for explaining the operation of the medical image diagnostic apparatus according to the second embodiment of the present invention. The figure which showed the example of 2 direction Doppler measurement.
  • FIG. 1 The figure which showed the example of the left ventricle measurement in the image data of M mode.
  • Figure search range setting unit by the contour extraction process is contour extraction, search range setting section showed that setting the search range.
  • FIG. 1 is a block diagram illustrating an ultrasonic diagnostic apparatus (medical image diagnostic apparatus) 1 as an example of an image capturing unit 20 that captures a medical image of a subject according to the first embodiment of the present invention.
  • the ultrasonic diagnostic apparatus 1 shown in FIG. 1 includes a probe 2, an ultrasonic signal generation unit 3, an ultrasonic image generation unit 4, a storage unit 5, an output / display unit 6, an input unit 7, A region-of-interest setting unit 8, a region-of-interest database 12 including a reference coordinate storage unit, a measurement unit 13, and a control unit 14 are provided.
  • the region-of-interest setting unit 8 includes a search range setting unit 9, an image reduction unit (image data compression unit) 10, and a region-of-interest setting unit 11.
  • the probe 2 is a device that transmits and receives ultrasonic waves from the transducer toward the target tissue, and performs scanning methods such as a linear type, a convex type, and a sector type.
  • the ultrasonic signal generation unit 3 transmits an electrified ultrasonic signal to and from the probe 2.
  • transmission / reception power and timing information from the control unit 14 transmission / reception of ultrasonic waves via the probe 2 is controlled so that a desired ultrasonic signal is generated.
  • the ultrasonic image generation unit 4 generates an ultrasonic image from the ultrasonic signal generated by the ultrasonic signal generation unit 3 based on the imaging setting of the apparatus (for example, the scanning range of the ultrasonic beam and the gain setting). .
  • This ultrasonic image is constantly updated according to the frame rate determined by the imaging setting, and is displayed on the screen of the ultrasonic diagnostic apparatus as an image by the output / display unit 6.
  • the storage unit 5 stores signal data and image data acquired from the subject such as generated ultrasonic signals, ultrasonic images (medical images), and electrocardiograms. Further, the storage unit 5 stores programs for operating various systems constituting the ultrasonic diagnostic apparatus 1.
  • the storage unit 5 is a storage medium such as a semiconductor memory, a hard disk, or an optical disk. Furthermore, the storage unit 5 may be an external storage medium through a network.
  • the output / display unit 6 displays an ultrasonic image (medical image) on the screen of the device, displays the region of interest set by the region of interest setting unit 8 superimposed on the ultrasonic image, and displays a measurement unit The measured value of the living tissue measured in 13 is displayed on the ultrasonic screen.
  • the output / display unit 6 stores an ultrasonic image (medical image) and a measurement value, and outputs the ultrasonic image (medical image) and the measurement value as a measurement report.
  • the input unit 7 is an interface for performing various operations of the ultrasonic diagnostic apparatus 1.
  • the input unit 7 is used to select an item that the examiner wants to measure from the measurement menu.
  • the input unit 7 is an input device such as a keyboard, a trackball, a switch, or a dial.
  • the input unit 7 may include a voice input function.
  • the region-of-interest setting unit 8 sets a region of interest for measurement based on the measurement item set by the input unit 7.
  • Set the search range based on the region-of-interest setting values (reference coordinates that serve as the reference for the location of the region of interest) pre-stored in the region-of-interest database 12, and use it on multiple image data with different reduction rates (compression rates)
  • the region of interest is set recursively and functions to improve the accuracy of setting the region of interest including the characteristic part while increasing the spatial resolution of the setting position of the region of interest.
  • the region of interest setting unit 8 includes a search range setting unit 9, an image reduction unit (image data compression unit) 10, and a region of interest setting unit 11, and includes a search range setting unit 9, an image reduction unit 10, and operation of the region of interest setting unit 11, while changing the reduction ratio of the image data (compression), it is repeatedly performed on image data.
  • the operations of the search range setting unit 9, the image reduction unit 10, and the region-of-interest setting unit 11 are repeatedly performed from small (high) image data with a small reduction rate (compression rate) to large (low) image data.
  • the search range setting unit 9 sets a search range in the image data.
  • the search range setting unit 9 determines a search range for performing template matching in the image data.
  • the initial value of the search range used for determining the first search range (initial search range) (for example, reference coordinates serving as a reference for the position of the region of interest) is read from the region of interest database 12 and used. Further, the search range setting unit 9 changes the search range in accordance with the image reduction rate (compression rate). Furthermore, the search range setting unit 9 performs contour extraction processing on the image data in order to set a region of interest.
  • Image reduction unit (the image data compression unit) 10 performs compression of image data on the basis of a plurality of pixels, to vary the compression ratio of the image data.
  • An image reduction unit (image data compression unit) 10 generates a plurality of reduced image data (compressed image data) having different reduction rates (compression rates).
  • a reduction method compression method
  • an existing method such as a nearest neighbor method or a bicubic method is used.
  • the region-of-interest setting unit 11 sets a region of interest including a characteristic part (such as a heart annulus) in the search range set by the search range setting unit 9.
  • the region-of-interest setting unit 11 uses pattern identification (such as template matching) to set a position that is most likely to be a region of interest within the search range.
  • the region of interest is set in the search range set by the search range setting unit 9 on the reduced image (compressed image) generated by the image reduction unit (image data compression unit) 10.
  • template matching methods include a method based on principal component analysis such as a SAD (Sum of Absolute Difference) method, a subspace method, and a boosting method.
  • Pattern recognition includes k-nearest neighbor discrimination method, subspace method, discrimination function optimization method, decision tree method, and neural network method.
  • the region-of-interest database 12 includes a reference coordinate storage unit that stores reference coordinates serving as a reference for the position of the region of interest.
  • the region of interest database 12 stores historical coordinate data of the region of interest or characteristic site based on a plurality of sample data as the reference coordinates.
  • the region-of-interest database 12 stores region-of-interest template data based on a plurality of past sample data.
  • the region-of-interest setting unit 11 uses the SAD method, the correlation with the template data may be calculated.
  • the region-of-interest setting unit 11 uses the subspace method or Boosting, the learning result is stored in the region-of-interest database. It can be stored in 12 and used for matching calculation.
  • the measurement unit 13 calculates a measurement value for the measurement item set by the input unit 7.
  • description will be mainly made on the measurement of the heart, but it can be used for measuring blood vessels and tumors in addition to the heart. For example, the size and movement of the organ are measured.
  • the measurement value is displayed on the screen of the apparatus by the output / display unit 6 together with the ultrasonic image (medical image).
  • the control unit 14 controls the entire system and controls the synchronization of a series of processes of the region-of-interest setting unit 8, the region-of-interest database 12, the measurement unit 13, and the output / display unit 6 from the setting of measurement items.
  • the control unit 14 uses a control device such as a CPU.
  • the ultrasonic diagnostic apparatus (medical image diagnostic apparatus) 1 will be described.
  • a description will be given of a case where a physical quantity such as a heart volume or a velocity is measured for setting a region of interest including a characteristic part in the heart.
  • FIG. 2 is a flowchart for explaining the operation of the ultrasonic diagnostic apparatus 1 according to the present embodiment.
  • the present embodiment is a method for setting a region of interest according to a measurement item set by an examiner, wherein a template matching search range is set based on reference coordinates stored in the region of interest database 12. Decide and set the region of interest recursively on multiple reduced images (compressed images) to improve the accuracy of setting the region of interest including the featured part while increasing the spatial resolution of the setting position of the region of interest A method will be described.
  • the output / display unit 6 displays an image on which the target biological tissue is depicted on the display screen (step S101). For example, as shown in FIG. 3, when measuring the left ventricle, the output / display unit 6 displays the apex four-chamber image 202 on the display screen 201.
  • the input unit 7 sets measurement items (step S102). For example, as shown in a measurement item (measurement value column) 204 in FIG. 3, the input unit 7 sets LV (Left-Ventricle) Volume (left ventricular volume) as a measurement item.
  • LV Left-Ventricle
  • Volume left ventricular volume
  • the Modified Simpson method is used, extraction of the left ventricular contour 208 is required to measure the left ventricular volume.
  • a region of interest including features such as an annulus (septal side annulus 205, side wall annulus 206) and apex 207 is extracted.
  • the search range setting unit 9 sets an initial search range for performing template matching by setting an initial value of the search range in the image data (step S103). If the initial search range is set to be wide enough to include the region of interest, the region of interest can be set reliably, but the number of search pixels for setting the region of interest becomes excessive, and the search calculation time Becomes longer. Therefore, the initial search range should surely include the region of interest and be as narrow as possible.
  • the reference coordinate for example, based on a plurality of past sample data is stored from the reference coordinate storage unit of the region of interest database 12.
  • the region of interest or the coordinate data of the characteristic part is read out, and the initial search range is set based on the reference coordinates.
  • the initial search range may be set to a range of a predetermined number of pixels with reference coordinates (for example, the average of the coordinate values of sample data) as the center. Further, the range of the predetermined number of pixels may be a standard deviation of the coordinate values of the region of interest or the characteristic part based on a plurality of past sample data.
  • FIG. 4 is a diagram showing a comparative example with the present embodiment.
  • the region of interest is set without reducing (compressing) the image data.
  • the matching calculation is performed while moving the template 403 in the search range 402, and the best matching position is set as the region of interest.
  • the shape of the search range is not limited to a specific shape such as a circle or a rectangle, but in the present embodiment, a rectangle is described as an example.
  • the search range since the number of pixels included in the search range is large, it is necessary to perform a huge number of matching operations. That is, the number of search pixels becomes excessive and the search calculation time becomes long.
  • the region of interest can be reliably set in the search range, but this is not preferable because the calculation time for setting the region of interest becomes long.
  • a first search range for example, an initial search range
  • the region-of-interest setting unit 11 sets a region of interest (including a feature part) in the first search range.
  • the second search range including the region of interest set in the first search range is set
  • the region of interest setting unit 11 is the second A region of interest (including a feature part) is set in the search range.
  • the image reduction unit (image data compression unit) 10 compresses image data based on a plurality of pixels, and changes the compression rate of the image data.
  • a method of simply thinning out the pixels may be used, or the luminance change based on the weighted average of the neighborhood.
  • a smoothing method may be used. In the former case, the calculation time is shortened, and in the latter case, the accuracy of template matching is improved.
  • Other methods for reducing (compressing) image data include the bicubic method (bicubic interpolation), the bilinear method (bilinear interpolation), and the nearest neighbor method (nearest neighbor interpolation).
  • the image reduction unit (image data compression unit) 10 compresses the image data to generate a reduced image (compressed image) (step S104).
  • the image reduction unit (image data compression unit) 10 compresses the image data of the ultrasonic image (medical image) generated by the ultrasonic image generation unit 4 to 1/4 (length and length of horizontal and horizontal sides respectively 1/2).
  • image data 501 is generated.
  • the initial search range created in step S103 is compressed to 1/4 and set as the search range 504. That is, the search range (initial search range) 504 is set in the image data after compression so as to cover the same range as the image data before compression acquired by the ultrasound diagnostic apparatus 1.
  • the coordinates of the center point 505 of the initial search range 504 are reference coordinates stored in the region of interest database 12.
  • the coordinates of the center point 505 are the coordinates of a region of interest or a characteristic part (such as an annulus or apex) based on a plurality of past sample data. That is, the center point 505 of the image data after compression is set at the same position (for example, the position of the same characteristic part) as the image data before compression acquired by the ultrasonic diagnostic apparatus 1.
  • the search range is determined from the reference coordinates of the region of interest stored in the region of interest database 12.
  • a search range is set such that the average (reference coordinates) of the coordinate values of the region of interest or feature portion based on a plurality of past sample data is the center point 505 and the standard deviation of the coordinate values is the search range.
  • a search range that includes all or a part of the coordinate data of the region of interest or the characteristic part based on a plurality of past sample data (for example, 70% of the coordinate data) may be set.
  • a template 506 compressed to 1/4 is set. That is, the template 506 is set in the image data after compression so as to cover the same range as the image data before compression acquired by the ultrasonic diagnostic apparatus 1.
  • the image reduction unit (image data compression unit) 10 reduces (compresses) the image data, instead of reducing (compressing) all of the image data acquired by the ultrasound diagnostic apparatus 1, as shown in FIG.
  • the image reduction unit (image data compression unit) 10 may reduce (compress) the image data or template within the search range.
  • the image reduction unit (image data compression unit) 10 may reduce (compress) the image data by thinning out a search range or a template pixel. In this case, if the image data is reduced by 1/2, one pixel may be thinned out every two pixels, and if the image data is reduced by 1/4, three pixels may be thinned out every four pixels.
  • FIG. 6 is a diagram showing an example of 1/4 reduction of image data.
  • the search range 602 in FIG. 6 is a search range in which 3 pixels out of 4 pixels (2 vertical pixels ⁇ 2 horizontal pixels) are removed (thinned out), and 64 pixels (8 vertical pixels ⁇ 8 horizontal pixels). Image data is compressed to 16 pixels (4 vertical pixels x 4 horizontal pixels).
  • the template 604 in FIG. 6 is a template in which 3 pixels are removed (decimated) out of 4 pixels (2 vertical pixels ⁇ 2 horizontal pixels), and 16 pixels (4 vertical pixels ⁇ 4 horizontal pixels). Image data is compressed to 4 pixels (2 vertical pixels x 2 horizontal pixels).
  • the method is substantially the same as the image reduction by the simple thinning method described above, but this method may be selected in consideration of the ease of mounting.
  • the region-of-interest setting unit 11 of the ultrasonic diagnostic apparatus 1 sets a region of interest including a characteristic part in the search range for the reduced (compressed) image data (step S105). For example, as shown in FIG. 5, the region-of-interest setting unit 11 moves the template 506 in the search range (initial search range) 504 of the image data 501 compressed to 1 ⁇ 4, and sets the most suitable position to the region of interest. Set as the position.
  • a known method such as the SAD method, the subspace method, or Boosting is used, and the learning data used for this may be stored in the region of interest database 12 in advance.
  • the template 506 starts from the upper left of the search range 504, and moves from left to right and from top to bottom.
  • the template 702 may start from the center point 703 of the search range 701 and move in a spiral shape. Since the center point 703 is the coordinates (reference coordinates) of the region of interest or characteristic part based on a plurality of past sample data, there is a high possibility that the region of interest position is set near the center point 703 of the search range 701. Therefore, the region of interest may be searched by moving the template 702 spirally, a position having a correlation with the template 702 equal to or greater than a predetermined threshold may be set as the region of interest position, and the search may be terminated. As a result, the search time for setting the region of interest can be shortened.
  • the region-of-interest setting unit 8 of the ultrasonic diagnostic apparatus 1 determines whether the setting of the region of interest or the characteristic part is finished (step S106).
  • the region-of-interest setting unit 11 sets a region of interest including a characteristic part for image data reduced (compressed) to 1/4, the setting resolution is also 1/4. Accordingly, as the setting resolution is reduced reduction ratio (compression ratio) is small, the set position of the region of interest or characteristic site is roughly set. Whether to end the setting of the region of interest or the characteristic part may be determined by at least one of the reduction ratio (compression ratio), the setting resolution, the setting processing time, and the number of setting processes.
  • the search range setting unit 9 sets the search range so as to include the region of interest (characteristic part) set in step S105.
  • the region of interest (feature part) set in the initial search range (first search range) 504 of the image data 501 compressed at a compression ratio of 1/4 (first compression ratio) are the coordinates of the center point 508 of the search range (second search range) 507 in the image data 502 compressed at the compression rate 1/2 (second compression rate).
  • the search range 507 may be set to a range of a predetermined number of pixels with the center point 508 as the center. Further, the predetermined number of pixels may be smaller than the standard deviation of the coordinate values of the region of interest or the characteristic part based on a plurality of past sample data. In addition, a search range that includes a part of coordinate data of a region of interest or a characteristic part based on a plurality of past sample data (for example, 60% of the coordinate data) may be set. In this case, the range covered by the search range (second search range) 507 may be smaller than the range covered by the search range (first search range) 504.
  • the number of coordinate data covered by the search range (second search range) 507 (for example, 60% of the coordinate data) is the number of coordinate data covered by the search range (first search range) 504 (for example, coordinate data) Less than 70%). This is because the position of the region of interest or feature part is set roughly in the search range (first search range) 504, so the position of the region of interest or feature part is near the center point 508 that is the roughly set position. This is because there is a high possibility that it exists.
  • the image reduction unit (image data compression unit) 10 compresses the image data based on a plurality of pixels, and changes the compression rate of the image data.
  • the image reduction unit (image data compression unit) 10 halves the image data of the ultrasonic image (medical image) generated by the ultrasonic image generation unit 4.
  • the image data 502 is generated by compression. That is, the second compression rate (compression rate 1/2) is larger than the first compression rate (compression rate 1/4).
  • the region-of-interest setting unit 11 sets a region of interest in the second search range. For example, as shown in FIG. 5, the region-of-interest setting unit 11 moves the template 509 in the search range (second search range) 507 of the image data 502 compressed to 1 ⁇ 2, and finds the most suitable position. The position of the region of interest is set (step S105).
  • ROI setting unit 8 of the ultrasonic diagnostic apparatus 1 determines whether to end the setting of the region of interest or characteristic portion (step S106).
  • the reduction rate compression rate
  • the ultrasound image medical image
  • a region of interest or a characteristic part is set by changing.
  • the data volume is 1/4 (compression ratio 1/4) or the data volume 1/2 (compression ratio 1/2) and the position accuracy is sufficient to meet the requirements
  • the reduction ratio in the middle The processing for setting the region of interest or the characteristic part may be finished with (compression rate).
  • the reduction ratio compression ratio
  • the reduction ratio compression ratio
  • the position of the region of interest feature part
  • the coordinates are the coordinates of the center point 511 of the search range (second search range) 510 in the image data 503 compressed at the compression rate 1/1 (second compression rate), and the region of interest setting unit 11
  • the template 512 is moved within the search range (second search range) 511 of the 1/1 image data 503, and the most suitable position is set as the position of the region of interest (step S105).
  • the finally set position is the position of the region of interest including the characteristic parts (the annulus parts 205 and 206 and the apex part 207) in FIG.
  • the reduction ratio (compression ratio) of the image data is set to three stages of 1/4, 1/2, and 1/1 has been described, but an arbitrary reduction ratio (compression ratio) is set. It is possible.
  • FIG. 8 shows the relationship between the reduction rate of the image data (for example, composed of 128 ⁇ 128 pixels and 256 ⁇ 256 pixels) and the cumulative number of search pixels according to the present embodiment.
  • the reduction ratios in FIG. 8 are six levels of 1/50, 1/40, 1/30, 1/20, 1/10, and 1/1. Searching the number of pixels is the number of pixels within the search range is 25 pixels if the search range of the vertical 5 pixels ⁇ 5 horizontal pixels.
  • the cumulative search pixel numbers are 25 (reduction ratio 1/50), 50 (reduction ratio 1/50) from reduction ratio 1/50 to 1/1. 40), 75 (reduction ratio 1/30), 100 (reduction ratio 1/20), 125 (reduction ratio 1/10), and 150 (reduction ratio 1/1). Therefore, when the region of interest or the characteristic part is set from the reduction ratio 1/50 to 1/1, the total number of matching operations is 150.
  • the region-of-interest setting search time can be shortened by recursively setting the region of interest on a plurality of image data having different reduction ratios (compression ratios).
  • the reduction ratio is 6 stages, but it can be an arbitrary number of stages (for example, 4 stages of 1/40, 1/30, 1/20, and 1/10).
  • a search range of a plurality of generated compressed image data can be arbitrarily selected, and a region of interest can be set in the medical image based on the selected search range.
  • the number of matching operations is 50 (reduction ratio 1/40), 75 (reduction ratio 1/30), 100 (reduction ratio 1/20), and 125 (reduction ratio 1/10). Therefore, when the region of interest or the characteristic part is set from the reduction ratio 1/40 to 1/10 of the image data, the total number of matching operations is 125.
  • the number of steps of the reduction rate can be changed according to the required setting accuracy and the allowable calculation time.
  • a plurality of stages in which search ranges of a plurality of generated image data after compression are arranged according to the compression ratio are formed a stage is arbitrarily selected from the formed stages, and the search range of the selected stage
  • the region of interest can be set in the medical image based on the above.
  • the number of stages of the reduction rate (compression rate) is used when the region-of-interest setting unit 8 determines whether or not to finish setting the region of interest or the characteristic part in step S106 of FIG.
  • the region of interest or feature part is set up to the image data with a reduction ratio of 1/1
  • the number of matching operations is 600 times.
  • the number of matching calculations is 2400.
  • step S106 when the region of interest or the characteristic part is set and it is determined that the setting is to be ended, the set region of interest is set as a desired region of interest (step S108). That is, the region of interest is set for measurement by the measurement unit 13.
  • the desired region of interest or characteristic part is a point
  • the point set in step S105 is the position of the region of interest.
  • a predetermined region of interest generation process is performed. For example, when the measuring unit 13 measures the volume of the heart, a region-of-interest generation process for extracting the outline of the heart is performed. As shown in FIG.
  • the measurement unit 13 performs measurement calculation for the measurement item set by the input unit 7 (step S109).
  • the measurement unit 13 calculates the volume value of the left ventricle by applying the Modified Simposon method using the above-described heart contour.
  • the output / display unit 6 displays the measurement result (step S110). For example, as shown in FIG. 3, the left ventricular volume value calculated by the measurement unit 13 is displayed in the measurement value column 204.
  • the set valve annulus parts 205 and 206, the apex part 207, and the contour line 208 may be displayed in a superimposed manner on the ultrasonic image.
  • the search range is set in steps S103 and S107 in FIG.
  • the region-of-interest database 12 including the reference coordinate storage unit stores reference coordinates serving as a reference for the position of the region of interest
  • the search range is set based on the coordinate difference between the set coordinates of the region of interest and the reference coordinates. May be.
  • the search range may be set so that the coordinate difference is minimized.
  • the position of the region of interest or the characteristic part is likely to exist in the vicinity of the region of interest (reference coordinates) in a plurality of past sample data. Therefore, the search range may be determined so that the coordinate difference between the reference coordinates based on a plurality of past sample data and the set region of interest is the smallest.
  • FIG. 9 is a diagram showing an example of determining a search range in image data with a reduction ratio of 1/40.
  • the search range setting unit 9 sets a search range with image data having a reduction ratio of 1/40.
  • the search range setting unit 9 enlarges the coordinate value of the search range set by image data with a reduction ratio 1/48 having a smaller reduction ratio by 48/40 times, and the same range as the reduction ratio 1/48 Set the search range to image data with a reduction ratio of 1/40.
  • the reference coordinates for example, the average of the coordinate values of the sample data stored in the region of interest database 12 are read out.
  • the region-of-interest setting unit 11 may calculate the coordinate difference between the position of the region of interest (reference coordinates) of the sample data and the set position of the region of interest for each sample data using a plurality of sample data. In this case, the average of a plurality of coordinate differences is calculated. Based on the plurality of search ranges, a coordinate difference is calculated, and the search range is set so that the coordinate difference is minimized. The search range is similarly set for other reduction ratios.
  • the region of interest setting unit 11 is based on the set position of the first region of interest.
  • the second region of interest is set limited to a part of the search range.
  • the search range 810 includes the septal valve annulus (first region of interest or first valve annulus) 205 of the heart
  • the search range 811 includes the heart side valve annulus of the heart. (Second region of interest or second annulus part) 206
  • septum side annulus part (first annulus part) 205 is on the left side of the side wall side annulus part (second annulus part) 206.
  • the region-of-interest setting unit 11 is limited to the search range 802 on the right side of the septal valve annulus (first valve annulus) 205, and the side valve annulus (second valve annulus). ) 206 is set.
  • the side wall side valve annulus 206 is set on the right side of the septum side valve annulus 205, when searching for the side wall side valve annulus 206, in the search range 802 on the right side of the septum side valve annulus 205, The setting of the side wall valve annulus 206 may be performed.
  • the region-of-interest setting unit 11 The septum-side valve annulus (second valve annulus) is set only in the search range on the left side of the annulus (first valve annulus) 206.
  • the region-of-interest database 12 stores the region-of-interest coordinates (reference coordinates) set in the past IMT measurement (sample data) and the template of the IMT measurement image. It only has to be set.
  • FIG. 11 is a diagram showing an example in which a region of interest 1005 for measuring the thickness of the intima-media complex 1004 is set on the carotid artery image data 1002.
  • a region of interest suitable for the measurement items can be set.
  • the initial search range based on the reference coordinates stored in the region-of-interest database 12 and changing the reduction rate (compression rate) and performing the search (search) of the region of interest
  • the probability of the region of interest existing A region of interest can be set with a high search range.
  • the search range in each reduced image data compared to image data that is not reduced (reduction ratio 1/1 image data) Since the number of matchings can be reduced, high-speed setting processing can be performed.
  • a medical image of a subject is captured by the image capturing unit 20, the image data of the medical image of the subject captured by the image data compression unit 10 is used as image data before compression, and the uncompressed image data
  • the compressed image data is generated by compressing based on a plurality of pixels
  • the search range setting unit 9 sets the search range of the image data after compression
  • the search range of the image data before compression is set.
  • the region setting unit 8 sets a region of interest in the medical image based on the search range of the image data before compression and the search range of the image data after compression.
  • the medical image diagnostic apparatus of the present invention includes an image imaging unit 20 that captures a medical image of a subject, and image data of the captured medical image of the subject as uncompressed image data.
  • the image data compression unit 10 that generates the compressed image data by compressing based on the plurality of pixels of the image data, and sets the search range of the compressed image data, and sets the search range of the image data before compression
  • the search range setting unit 9 sets the second search range with a size equal to the first search range.
  • the medical image diagnostic apparatus of the present invention includes a reference coordinate storage unit 12 that stores reference coordinates serving as a reference for the position of the region of interest, and the search range setting unit 9 sets the search range based on the quasi-coordinates. To do.
  • the position of the initial search range can be set by using the reference coordinates, and the search range can be set to include the reference coordinates.
  • the search range setting unit 9 is centered on the reference coordinate that is an average of the coordinate values of the region of interest based on sample data, and the search is performed for the standard deviation of the coordinate value of the region of interest. Set as a range.
  • the position of the search range can be set by centering on the reference coordinates, and the search range can be set by using the standard deviation of the coordinate value of the region of interest based on the sample data. it can.
  • the search range setting unit 9 sets, as the center of the search range, a reference coordinate that is a coordinate value of a characteristic part set in the contour extraction process.
  • the search range can be narrowly adjusted while maintaining the positional relationship of the feature part even if the contour extraction process is not highly accurate.
  • the accuracy when setting the area is also improved.
  • the medical image diagnostic apparatus of the present invention includes a reference coordinate storage unit 12 that stores reference coordinates serving as a reference for the position of the region of interest, and the search range setting unit includes the set coordinates of the region of interest and the reference coordinates.
  • the search range is set based on the coordinate difference between and.
  • the search range is set based on the coordinate difference between the set coordinates of the region of interest and the reference coordinates. , Improve accuracy when setting the region of interest.
  • the search range setting unit 9 sets the search range so that the coordinate difference is minimized.
  • the search range is set so that the coordinate difference between the reference coordinates and the set region of interest is the smallest, so that the accuracy in setting the region of interest is improved.
  • the region of interest setting unit 8 is based on the position of the first region of interest, A second region of interest is set limited to a part of the search range.
  • the region of interest setting time can be shortened by setting the region of interest limited to a part of the search range.
  • the region-of-interest setting unit 8 uses the region including the first valve annulus of the heart as the first region of interest, and sets the region including the second valve annulus of the heart as the first region of interest.
  • the second annulus is limited to the search range on the right side of the first annulus part.
  • the first valve ring part is located on the right side of the second valve ring part, the first valve ring part is limited to the search range on the left side of the first valve ring part.
  • the valve annulus since the valve annulus is present on the left and right respectively, when one valve annulus is set, the other valve is limited to the search range on the left or right side of the position of the valve annulus. What is necessary is just to set a ring part. As a result, the search range can be narrowed, and the time for setting the region of interest can be shortened.
  • the medical image diagnostic apparatus of the present invention includes a measuring unit 13 that measures a measurement item in the set region of interest, and the measuring unit includes at least one of the first search range and the second search range. After setting the region of interest, the measurement result of the measurement item is output.
  • the measurement result of the measurement item may be output, and after setting the region of interest in the first search range and the second search range, The measurement result of the measurement item may be output, and the setting of the region of interest can be stopped at a desired timing according to the output result.
  • the region-of-interest setting unit 8 finishes setting the region of interest based on the appropriate value of the measurement result.
  • a medical image of a subject is captured by the image capturing unit 20, and image data of the medical image of the subject captured by the image data compression unit 10 is an image before compression.
  • the data is compressed based on a plurality of pixels of the uncompressed image data to generate compressed image data
  • the search range setting unit 9 sets the search range of the compressed image data
  • An image data search range is set
  • a region of interest setting unit 8 sets a region of interest in the medical image based on the search range of the image data before compression and the search range of the image data after compression. It is characterized by that.
  • a region of interest is set on B-mode image data in order to acquire Doppler or M-mode image data.
  • a region of interest including a predetermined measurement position is set in image data of measurement values related to measurement items in Doppler and M mode.
  • a plurality of regions of interest may be set by combining image data related to the B mode and measurement values.
  • Image reduction unit (image data compression unit) 10 compresses image data based on a plurality of pixels, and changes the compression rate (reduction rate) of the image data.
  • the search range setting unit 9 sets a search range in the compressed image data.
  • the region-of-interest setting unit 11 sets the region of interest including the characteristic part in the search range. Then, the image reduction unit (image data compression unit) 10 compresses the image data (B mode or measurement value image data) at the first compression rate, and the search range setting unit 9 compresses at the first compression rate.
  • a first search range is set, and the region-of-interest setting unit 11 sets a region of interest in the first search range.
  • the image reduction unit (image data compression unit) 10 compresses the image data (B mode or measurement value image data) at the second compression rate
  • the search range setting unit 9 compresses at the second compression rate.
  • a second search range including the region of interest set in the first search range is set, and the region-of-interest setting unit 11 sets the region of interest in the second search range.
  • FIG. 12 is a flowchart for explaining the operation of the medical image diagnostic apparatus according to the present embodiment.
  • FIG. 13 is a diagram illustrating an example of two-way Doppler measurement.
  • FIG. 14 is a diagram illustrating an example of left ventricular measurement in M-mode image data.
  • an image depicting the target biological tissue is displayed on the display screen of the output / display unit 6 (step S201).
  • FIG. 13 shows the image data of the apex four-chamber image
  • FIG. 14 shows the image data of the parasternal long axis image.
  • the measurement items are set by the examiner using the input unit 7 (step S202). For example, in the two-way Doppler of FIG. 13, when the inflow blood flow velocity and the annulus motion velocity are measurement items, it is necessary to set the Doppler gate 1203 on the valve opening and the annulus. Also, as shown in FIG. 14, when the left ventricular wall or the size of the cavity is a measurement item in the M mode, it is necessary to set the M mode beam line 1303 on the B mode image data.
  • the region-of-interest setting unit 8 sets a region of interest including a region of interest (such as an annulus and a valve opening) where the Doppler gate 1203 and the M-mode beam line 1303 are set.
  • the output / display unit 6 displays a Doppler image or an M-mode image.
  • the region-of-interest database 12 including the reference coordinate storage unit stores the positions (reference coordinates) and templates of the annulus, the valve opening, and the M-mode beam line acquired as a plurality of past sample data.
  • a reference coordinate (for example, coordinate data of a region of interest or a feature part based on a plurality of past sample data) is read from the reference coordinate storage unit of the region of interest database 12, and a search range is set based on the reference coordinate.
  • the region of interest is set while changing the reduction rate (compression rate) of the image data, and the region of interest is set on the B-mode image data (step S208).
  • the region of interest including the measurement position of the measurement item is set in the image data related to the measurement item in Doppler or M mode (steps S303 to S308).
  • image data of TDI (measurement value) in the region of interest 1207 set in the B mode image is shown, and in the lower right of FIG. 13, image data of pulse Doppler (measurement value) is shown. Yes. Further, on the right side of FIG. 14, M-mode image data in the M-mode beam line 1303 of the region of interest 1306 set in the B-mode image is displayed.
  • the region-of-interest setting unit 8 sets a region of interest including a predetermined measurement position in the image data of the measurement values related to the Doppler and M-mode measurement items. That is, the region-of-interest setting unit 8 sets a region of interest (measurement position) in the image data.
  • a waveform position 1206 is set as a region of interest (measurement position) of E wave, A wave, e ′ wave, and DT (E wave deceleration time).
  • a wall position 1308 is set as a region of interest (measurement position) of the inner and outer membranes of the myocardium.
  • the region-of-interest database 12 stores the E wave, A wave, e 'wave, and DT positions (reference coordinates) and the template 1208 in the Doppler image, and the inner and outer membrane positions (reference coordinates) and the template in the M mode image. 1307 is stored.
  • a reference coordinate (for example, coordinate data of a region of interest or a feature part based on a plurality of past sample data) is read from the reference coordinate storage unit of the region of interest database 12, and a search range is set based on the reference coordinate.
  • the region of interest is set while changing the reduction rate (compression rate) of the image data, and the region of interest (measurement position) is set on the image data relating to the measurement items of Doppler and M mode (step S308).
  • the region of interest (measurement position) is a point
  • the region of interest set by the region of interest setting unit 11 becomes the measurement position for measurement as it is.
  • a predetermined region of interest generation process is performed.
  • the measurement unit 13 performs measurement calculation in the set region of interest (S409).
  • the measurement unit 13 calculates E wave, A wave, e ′ wave, and DT values from the Doppler waveform at the waveform position (measurement position) 1206. Further, the measurement unit 13 may calculate E / e ′, which is a combination of measurement values.
  • the measurement unit 13 calculates a measurement value related to the size of the myocardium and heart chamber such as LVEdD (left ventricular end diastolic diameter) at the wall position (measurement position) 1308.
  • LVEdD left ventricular end diastolic diameter
  • the output / display unit 6 displays the measurement value calculated by the measurement unit 13 (step S410). They are displayed in the measured value column 1202 in FIG. 13 and the measured value column 1302 in FIG. 14, respectively. Further, the output / display unit 6 may superimpose and display the set regions of interest (Doppler gate 1203, M-mode beam line 1303, waveform position 1206, wall position 1308) on the ultrasound image.
  • the set regions of interest Doppler gate 1203, M-mode beam line 1303, waveform position 1206, wall position 1308
  • a region of interest is set on the B-mode image data.
  • a region of interest including a predetermined measurement position is set in image data of measurement values related to measurement items in Doppler and M mode.
  • the region of interest for measuring the measurement item is set in the Doppler or M mode image data. be able to.
  • a plurality of regions of interest can be set for various image data.
  • Doppler image data or M-mode image data at a position desired by the examiner is generated, and at a predetermined measurement position according to the measurement item input by the examiner.
  • region of interest it is possible to improve the accuracy of region-of-interest setting and shorten the search time. Therefore, it is possible to provide a medical image diagnostic apparatus that reduces the operation burden of the region-of-interest setting of the examiner.
  • the reduction rate (compression rate) is changed to 1/1, but each reduction rate (compression rate) is in the process of gradually increasing the reduction rate (compression rate) of the image data.
  • the measurement results in may be sequentially output, and using the output information of the measurement results, the region of interest setting process may be terminated at a reduction rate (compression rate) that satisfies the measurement accuracy. That is, in the first embodiment, the measurement unit 13 outputs the measurement result of the measurement item after setting the region of interest in the second search range, but in the first search range and the second search range, After setting the region of interest, the measurement result of the measurement item may be output.
  • FIG. 15 is a flowchart for sequentially outputting measurement results and determining whether or not to finish setting the region of interest according to the measurement results.
  • step S106 for determining whether to end the setting is executed after step S110 for displaying the measurement result. That is, the measurement unit 13 outputs the measurement results of the measurement items after setting the region of interest in the first search range and the second search range. Using the region of interest set at a predetermined reduction ratio, the measurement unit 13 performs measurement, and the output / display unit 6 displays the measurement result. If this measurement result is appropriate, the measurement is terminated. If the measurement result is not appropriate, the process returns to step S107, and the region of interest is set at the next reduction rate (compression rate).
  • the determination of whether or not the measurement result is appropriate may be made by the examiner confirming the measurement value, or by comparing it with the appropriate value (threshold value) of the measurement result input by the input unit 7. Alternatively, it may be performed by comparing with an appropriate value (threshold value) of the measurement result stored in the storage unit 5. Further, it may be determined whether or not the value approaches the appropriate value (threshold value), and the setting process may be terminated when the value approaches the appropriate value (threshold value).
  • the region-of-interest setting unit 11 finishes setting the region of interest based on the appropriate value of the measurement result.
  • a method for extracting the outline of a living tissue may be used.
  • the region of interest including the annulus and the apex (feature portion) is set using a method for extracting the intimal contour of the heart.
  • FIG. 16 is a flowchart for explaining setting of a region of interest using a contour extraction method.
  • a step S502 for performing contour extraction processing is added as compared with the first embodiment.
  • the contour extraction process is preferably a contour with a shape restriction such as a dynamic contour model.
  • a model that captures the features of the contour shape of the left ventricle By constraining with a model that captures the features of the contour shape of the left ventricle, the contour unique to the left ventricle can be extracted, and the correct positional relationship between the annulus and the apex can be set.
  • This model may be stored in the storage unit 5 and referred to when the contour is extracted.
  • the search range setting unit 9 performs contour extraction processing (step S502). As shown in FIG. 17, in step S502 for performing contour extraction processing, the search range setting unit 9 extracts a contour 1603.
  • the extracted contour 1603 may be extracted so that the contour line fits the intima in a state in which the positional relationship between the annulus 1601 and the apex 1602 is maintained, even if it is not highly accurate.
  • the search range setting unit 9 sets the search range 1604 (step S107).
  • the search range setting unit 9 sets a search range 1604 centered on the coordinate values (reference coordinates) of the annulus 1601 and the apex 1602 set in the contour extraction process.
  • the search range 1604 may be set using past sample data in the region-of-interest database 12 as in the first embodiment. Since the positions of the annulus 1601 and the apex 1602 are set in the contour extraction process, the search range 1604 may be set to be narrower than the standard deviation of the sample data.
  • the positions of the annulus 1601 and the apex 1602 are set by the contour extraction process, so that the positional relationship between the annulus 1601 and the apex 1602 is maintained even if the outline extraction process is not highly accurate. If the set positions of the valve annulus 1601 and the apex 1602 are used, the search range 1604 can be adjusted narrowly, and the accuracy in setting the annulus 1601 and the apex 1602 is improved. Therefore, since the search range 1604 can be set narrowly, it is possible to improve setting accuracy and reduce setting calculation time.
  • the measurement of the heart in the ultrasonic diagnostic apparatus has been described as an example, but the present invention can be similarly applied to other diagnostic apparatuses and biological tissues.
  • the image data may be at least one of a tomographic image, an elasticity image, and a measurement value image of the subject.
  • Medical images include X-ray images, CT, MRI, ultrasound images (US), angiography (angiography), and the like. That is, the image capturing unit 20 is an X-ray image capturing apparatus, an X-ray CT apparatus, an MRI apparatus, or the like.
  • the medical image diagnostic apparatus sets the search range using the compressed image data and sets the region of interest by template matching or the like, thereby improving the accuracy of the region of interest setting and shortening the search time. Therefore, it is possible to provide a medical image diagnostic apparatus that reduces the burden of setting the region of interest of the examiner, and is useful as a medical image diagnostic device that sets a region of interest including a characteristic part or a measurement position. It is.
  • 1 ultrasonic diagnostic device 2 probe, 3 ultrasonic signal generation unit, 4 ultrasonic image generation unit, 5 storage unit, 6 output / display unit, 7 input unit, 8 region of interest setting unit, 9 search range setting unit , 10 image reduction unit (image data compression unit), 11 region of interest setting unit, 12 region of interest database, 13 measurement unit, 14 control unit, 20 image capturing unit, 203 heart

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Cardiology (AREA)
  • Human Computer Interaction (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Multimedia (AREA)
  • Quality & Reliability (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

　本発明は、画像撮像部によって被検体の医用画像を撮像し、画像データ圧縮部によって撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成し、探索範囲設定部によって圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定し、関心領域設定部によって前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する。

Description

医用画像診断装置及びその関心領域設定方法
　本発明は、医用画像診断装置に関し、特に、特徴部位又は計測位置を含む関心領域を設定する医用画像診断装置及びその関心領域設定方法に関する。
　被検体の生体組織の物理量情報は、医用画像診断装置で計測され、病変の程度を診断する重要な情報である。物理量情報とは、被検体の生体組織の長さ、面積、容積、及び動きなどである。
　従来の技術では、検者が医用画像診断装置の入力装置(例えば、マウスやトラックボール)を用いて、計測したい位置や領域(関心領域)を表示部の画面に表示された医用画像上の所望の領域に手動設定する。
　関心領域の手動設定は、検者の主観でなされるため、すべての検者による関心領域の手動設定が同一の条件下で同じ関心領域を設定するとは限らない。
　そこで、医用画像診断装置が医用画像の全体の画像データに対して臓器領域の輪郭を探索し、その探索結果に基づきテンプレートマッチング等に応じて関心領域を設定して、心臓の弁輪部の位置を基準にした左心室内膜の輪郭を抽出する手法が特許文献1に提案されている。
特開2005-218796号公報
　しかしながら、特許文献1の医用画像の全体の画像データに対するテンプレートマッチングだけでは、リアルタイムで関心領域を設定する場合、関心領域の設定のための探索時間短縮は未解決の問題であると思料する。
　そこで、本発明の目的は、リアルタイムで関心領域を設定する場合、関心領域の設定のための探索時間短縮を達成する医用画像診断装置及びその関心領域設定方法を提供することにある。
　本発明は、画像撮像部によって被検体の医用画像を撮像し、画像データ圧縮部によって撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成し、探索範囲設定部によって圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定し、関心領域設定部によって前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する。
　本発明によれば、設定された圧縮前後の画像データの探索範囲を用いて関心領域を設定することで、リアルタイムで関心領域を設定する場合、関心領域の設定のための探索時間短縮を達成する医用画像診断装置及びその関心領域設定方法を提供できる。
本発明の第1の実施の形態の超音波診断装置(医用画像診断装置)を例示するブロック図。 本発明の第1の実施の形態に係る医用画像診断装置の動作について説明するフローチャート。 出力・表示部が表示する生体組織の画像データの一例を示した図。 第1の実施の形態との比較例を示した図。 画像データの圧縮率(縮小率)を変化させながら、探索範囲設定部が、圧縮された画像データにおいて探索範囲を設定することを示す図。 画像縮小部(画像データ圧縮部)が、探索範囲内の画像データやテンプレートを縮小(圧縮)することを示した図。 テンプレートが、中心点からスタートして渦巻き状に移動することを示す図。 本発明の第1の実施の形態の画像データの縮小率と累積探索画素数の関係を示す図。 複数の探索範囲に基づいて座標差が算出され、座標差が最も小さくなるように探索範囲が設定されることを示す図。 一方の弁輪部が設定された場合は、その弁輪部の位置よりも左側又は右側の探索範囲に限定して他方の弁輪部を設定することを示した図。 IMT計測アプリケーションを用いる場合に、Bモード画像データ上に関心領域が設定されることを示した図。 本発明の第2の実施の形態に係る医用画像診断装置の動作について説明するフローチャート。 2方向ドプラ計測の例を示した図。 Mモードの画像データにおける左室計測の例を示した図。 計測結果を逐次出力して、計測結果に応じて関心領域の設定を終了するか否かを判断するフローチャート。 輪郭抽出方法を用いた関心領域の設定を説明するフローチャート。 輪郭抽出処理によって探索範囲設定部が輪郭抽出し、探索範囲設定部が探索範囲を設定することを示した図。
　(第1の実施の形態)
　以下、本発明の第1の実施の形態の医用画像診断装置について、医用画像診断装置の1例である超音波診断装置を用いて説明する。図1は、本発明の第1の実施の形態の被検体の医用画像を撮像する画像撮像部20の一例として超音波診断装置(医用画像診断装置)1を例示するブロック図である。
　図1に示す超音波診断装置1は、探触子2と、超音波信号生成部3と、超音波画像生成部4と、記憶部5と、出力・表示部6と、入力部7と、関心領域設定部8と、基準座標記憶部を含む関心領域データベース12と、計測部13と、制御部14とを備える。また、関心領域設定部8は、探索範囲設定部9と、画像縮小部(画像データ圧縮部)10と、関心領域設定部11とを含む。
　探触子2は、振動子から対象組織に向かって超音波を送信し受信する装置であり、リニア型、コンベックス型、及びセクタ型などの走査方法を実施する。超音波信号生成部3は、探触子2との間で、電気化された超音波信号を送信する。超音波信号生成部3が制御部14から送受信のパワーやタイミングの情報を受け取ることにより、所望の超音波信号が生成されるように、探触子2を介する超音波の送受信が制御される。
　超音波画像生成部4は、超音波信号生成部3により生成された超音波信号から、装置の撮像設定(例えば、超音波ビームの走査範囲やゲイン設定など)に基づいて超音波画像を生成する。この超音波画像は、撮像設定によって決定されるフレームレートに従って常時更新され、出力・表示部6によって映像として超音波診断装置の画面に表示される。
　記憶部5は、生成された超音波信号、超音波画像(医用画像)、及び心電図などの被検体から取得された信号データ及び画像データが記憶されている。また、記憶部5は、超音波診断装置1を構成する種々のシステムを動作させるためのプログラムを格納する。例えば、記憶部5は、半導体メモリ、ハードディスク、光ディスクなどの記憶媒体である。さらには、記憶部5は、ネットワークを通した外部記憶媒体でもよい。
　出力・表示部6は、装置の画面上に超音波画像(医用画像)を表示したり、関心領域設定部8で設定された関心領域を超音波画像上に重畳して表示したり、計測部13で計測された生体組織の計測値を超音波画面上に表示したりする。また、出力・表示部6は、超音波画像(医用画像)や計測値を記憶したり、超音波画像(医用画像)や計測値を計測レポートとして出力したりする。
　入力部7は、超音波診断装置1の各種操作を行うインターフェイスである。特に、本実施の形態では、検者が計測したい項目を計測メニューから選択するために、入力部7が用いられる。例えば、入力部7は、キーボード、トラックボール、スイッチ、ダイヤルなどの入力機器である。また、入力部7は、音声入力機能を含んでいてもよい。
　関心領域設定部8は、入力部7により設定された計測項目に基づいて、計測のための関心領域を設定する。関心領域データベース12に予め記憶されている関心領域設定値(関心領域の位置の基準となる基準座標など)に基づいて探索範囲を設定し、縮小率(圧縮率)の異なる複数の画像データ上で再帰的に関心領域を設定し、関心領域の設定位置の空間分解能を上げながら、特徴部位を含む関心領域の設定の正確性を向上させるように機能する。この機能を実施するために、関心領域設定部8は、探索範囲設定部9と画像縮小部(画像データ圧縮部)10と関心領域設定部11とを含み、探索範囲設定部9、画像縮小部10、及び関心領域設定部11の動作は、画像データの縮小率(圧縮率)を変化させながら、画像データに繰り返し実施される。探索範囲設定部9、画像縮小部10、及び関心領域設定部11の動作は、縮小率(圧縮率)の小さい(高い)画像データから大きい(低い)画像データに繰り返し実施される。
　探索範囲設定部9は、画像データにおいて探索範囲を設定する。探索範囲設定部9は、画像データにおいてテンプレートマッチングを行うための探索範囲を決定する。最初の探索範囲(初期探索範囲)を決定するために用いられる探索範囲の初期値(例えば、関心領域の位置の基準となる基準座標)は、関心領域データベース12から読み出して使用される。また、探索範囲設定部9は、画像の縮小率(圧縮率)に応じて探索範囲を変化させる。さらには、探索範囲設定部9は、関心領域を設定するために、画像データに輪郭抽出処理を実施する。
　画像縮小部(画像データ圧縮部)10は、複数の画素に基づいて画像データの圧縮を行い、画像データの圧縮率を変化させる。画像縮小部(画像データ圧縮部)10は、縮小率(圧縮率)の異なる複数の縮小画像データ(圧縮画像データ)を生成する。縮小方法(圧縮方法)は、最近傍法やバイキュービック法など既存の方法が利用される。
　関心領域設定部11は、探索範囲設定部9により設定された探索範囲において特徴部位(心臓の弁輪部など)を含む関心領域を設定する。関心領域設定部11は、パターン識別(テンプレートマッチングなど)を用いて、探索範囲の中で関心領域である可能性が最も高い位置を設定する。画像縮小部(画像データ圧縮部)10により生成された縮小画像(圧縮画像)上で探索範囲設定部9により設定された探索範囲において、関心領域の設定が実施される。例えば、テンプレートマッチングの方法として、SAD(Sum of Absolute Difference)法や、部分空間法などの主成分分析を基礎とした方法や、Boostingによる方法などがある。また、パターン認識として、k-最近傍識別法、部分空間法、識別関数の最適化による方法、決定木による方法、及びニューラルネットワークによる方法などがある。
　関心領域データベース12は、関心領域の位置の基準となる基準座標を記憶する基準座標記憶部を含む。例えば、関心領域データベース12は、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データを基準座標として記憶する。また、関心領域データベース12は、過去の複数のサンプルデータに基づく関心領域のテンプレートデータを格納する。関心領域設定部11がSAD法を用いる場合には、このテンプレートデータとの相関を計算すればよいし、関心領域設定部11が部分空間法やBoostingを用いる場合には、学習結果を関心領域データベース12に格納しておき、マッチング演算に利用すればよい。
　計測部13は、入力部7により設定された計測項目について計測値を算出する。本実施の形態では、主に心臓の計測を例に説明するが、心臓の他にも血管や腫瘍の計測にも利用可能である。例えば、臓器のサイズや動きが計測される。計測値は、出力・表示部6によって、超音波画像(医用画像)とともに装置の画面上に表示される。
　制御部14は、システム全体を制御するとともに、計測項目の設定から、関心領域設定部8、関心領域データベース12、計測部13、出力・表示部6の一連の処理の同期を制御する。例えば、制御部14は、CPU等の制御装置を用いたものである。
　次に、本実施の形態に係る超音波診断装置(医用画像診断装置)1の動作について説明する。本実施の形態では、心臓における特徴部位を含む関心領域の設定について、心臓の容積や速度などの物理量を計測する場合を説明する。
　図2は、本実施の形態に係る超音波診断装置1の動作について説明するフローチャートである。図2に示すように、本実施の形態では、検者が設定した計測項目に従って関心領域を設定する方法であって、関心領域データベース12に格納される基準座標に基づいてテンプレートマッチングの探索範囲を決定し、複数の縮小画像(圧縮画像)上で再帰的に関心領域の設定を行って、関心領域の設定位置の空間分解能を上げながら、特徴部位を含む関心領域の設定の正確性を向上させる方法を説明する。
　図2に示すように、出力・表示部6が、対象とする生体組織が描出された画像を、表示画面上に表示する(ステップS101)。例えば、図3に示すように、左室の計測を目的とする場合、出力・表示部6が、心尖部4腔像202を表示画面201に表示する。
　入力部7が計測項目を設定する(ステップS102)。例えば、図3の計測項目(計測値欄)204に示すように、入力部7が、LV(Left Ventricle)Volume(左室容積)を、計測項目として設定する。Modified Simpson法を用いる場合、左室容積を計測するために、左室の輪郭線208の抽出が要求される。超音波画像(医用画像)において左室輪郭線を抽出する場合、弁輪部(中隔側弁輪部205、側壁側弁輪部206)及び心尖部207のような特徴部位を含む関心領域を設定して、設定された関心領域に基づいて輪郭線を抽出することで、左心室内膜の輪郭を抽出する輪郭抽出の正確性が向上する。
　探索範囲設定部9が、画像データにおいて探索範囲の初期値を設定することで、テンプレートマッチングを行うための初期探索範囲の設定を行う(ステップS103)。初期探索範囲は、関心領域を確実に含むように広く設定されると、関心領域を確実に設定することができるが、関心領域を設定するための探索画素数が過剰になり、探索の計算時間が長くなる。したがって、初期探索範囲は、関心領域を確実に含み、且つ可能な限り狭いほうがよい。本実施の形態では、関心領域を確実に含み、且つ可能な限り狭い初期探索範囲を設定するために、関心領域データベース12の基準座標記憶部から基準座標(例えば、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データなど)が読み出され、初期探索範囲が基準座標に基づいて設定される。初期探索範囲は、基準座標(例えば、サンプルデータの座標値の平均)を中心とし、所定の画素数の範囲に設定されればよい。また、所定の画素数の範囲は、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標値の標準偏差であればよい。
　図4は、本実施の形態との比較例を示した図である。図4に示す探索方法では、画像データを縮小(圧縮)せずに関心領域を設定する。図4に示すように、探索範囲402の中でテンプレート403を移動させながらマッチング演算を行い、最も良く適合する位置が関心領域として設定される。なお、探索範囲の形状は円形や矩形など特定の形状に限定されないが、本実施の形態では矩形を例に説明する。比較例では、探索範囲に含まれる画素数が多いので、莫大な回数のマッチング演算を行う必要がある。つまり、探索画素数が過剰になり、探索の計算時間が長くなる。比較例のように、探索範囲を広くすれば確実に探索範囲内に関心領域を設定することができるが、関心領域を設定する計算時間が長くなり好ましくない。
　本実施の形態では、画像縮小部(画像データ圧縮部)10が、複数の画素に基づいて画像データの圧縮を行い、画像データの圧縮率を変化させ、探索範囲設定部9が、縮小(圧縮)された画像データにおいて探索範囲を設定し、関心領域設定部11が、特徴部位を含む関心領域の設定を探索範囲において行う。つまり、画像縮小部(画像データ圧縮部)10が第1の圧縮率(例えば、データ量1/4＝圧縮率1/4)で画像データを圧縮し、探索範囲設定部9が第1の圧縮率で圧縮された画像データにおいて、第1の探索範囲(例えば、初期探索範囲)を設定し、関心領域設定部11が第1の探索範囲において関心領域(特徴部位を含む)を設定する。そして、画像データの圧縮率を変化させて、画像縮小部(画像データ圧縮部)10が第2の圧縮率(例えば、データ量1/2＝圧縮率1/2)で画像データを圧縮し、探索範囲設定部9が第2の圧縮率で圧縮された画像データにおいて、第1の探索範囲で設定された関心領域を含む第2の探索範囲を設定し、関心領域設定部11が第2の探索範囲において関心領域(特徴部位を含む)を設定する。
　画像縮小部(画像データ圧縮部)10は、複数の画素に基づいて画像データの圧縮を行い、画像データの圧縮率を変化させる。画像縮小部(画像データ圧縮部)10が画像データを縮小(圧縮)するとき、単純に画素を間引く方法(単純間引法)を用いてもよいし、近傍の重み付け平均に基づいて輝度変化を滑らかにする方法を用いてもよい。前者であれば、計算時間が短くなり、後者であれば、テンプレートマッチングの正確性が向上する。
　この他、画像データを縮小(圧縮)する手法として、バイキュービック法(双三次補間)、バイリニア法(双一次補間)、ニアレストネイバー法(最近傍補間)などがある。
　図2に示すように、画像縮小部(画像データ圧縮部)10が、画像データを圧縮して縮小画像(圧縮画像)を生成する(ステップS104)。例えば、図5に示すように、画像データ501は、超音波診断装置1により取得された医用画像の画像データを1/4(＝圧縮率1/4)に圧縮した画像データである。つまり、画像縮小部(画像データ圧縮部)10が、超音波画像生成部4により生成される超音波画像(医用画像)の画像データを1/4(縦横辺長をそれぞれ1/2)に圧縮し、画像データ501を生成する。この圧縮画像上に、ステップS103で作成された初期探索範囲が1/4に圧縮され、探索範囲504として設定される。つまり、超音波診断装置1により取得された圧縮前の画像データと同じ範囲をカバーするように、圧縮後の画像データに探索範囲(初期探索範囲)504が設定される。
　初期探索範囲504の中心点505の座標は、関心領域データベース12に格納された基準座標である。中心点505の座標は、過去の複数のサンプルデータに基づく関心領域又は特徴部位(弁輪部や心尖部など)の座標である。つまり、超音波診断装置1により取得された圧縮前の画像データと同じ位置(例えば、同じ特徴部位の位置)に、圧縮後の画像データの中心点505が設定される。本実施の形態では、関心領域データベース12に記憶されている関心領域の基準座標から、探索範囲を決定するようにする。例えば、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標値の平均(基準座標)を中心点505とし、座標値の標準偏差を探索範囲とするような探索範囲が設定される。その他、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データの全て又は一部(例えば、座標データの70％)が含まれるような探索範囲を設定してもよい。
　関心領域データベース12に格納されたテンプレートデータに基づいて、1/4に圧縮されたテンプレート506が設定される。つまり、超音波診断装置1により取得された圧縮前の画像データと同じ範囲をカバーするように、圧縮後の画像データにテンプレート506が設定される。
　なお、画像縮小部(画像データ圧縮部)10が画像データを縮小(圧縮)するとき、図6に示すように、超音波診断装置1により取得された画像データの全てを縮小(圧縮)する代わりに、画像縮小部(画像データ圧縮部)10が、探索範囲内の画像データやテンプレートを縮小(圧縮)すればよい。例えば、画像縮小部(画像データ圧縮部)10が、探索範囲やテンプレートの画素(ピクセル)を間引くことで、画像データを縮小(圧縮)してもよい。この場合、画像データの1/2縮小ならば、2画素ごとに1画素を間引けばよいし、画像データの1/4縮小ならば、4画素ごとに3画素を間引けばよい。図6は、画像データの1/4縮小の例を示した図である。図6の探索範囲602は、4画素(縦2画素×横2画素)のうち3画素が取り除かれた(間引かれた)探索範囲であり、64画素(縦8画素×横8画素)の画像データが16画素(縦4画素×横4画素)に圧縮されている。また、図6のテンプレート604は、4画素(縦2画素×横2画素)のうち3画素が取り除かれた(間引かれた)テンプレートであり、16画素(縦4画素×横4画素)の画像データが4画素(縦2画素×横2画素)に圧縮されている。実質的には、上述の単純間引法による画像縮小と同様であるが、実装の容易さなどを考慮してこの方法が選択されてもよい。
　超音波診断装置1の関心領域設定部11が、縮小(圧縮)された画像データに対して、特徴部位を含む関心領域の設定を探索範囲において行う(ステップS105)。例えば、図5に示すように、関心領域設定部11は、1/4に圧縮された画像データ501の探索範囲(初期探索範囲)504でテンプレート506を移動させて、最も適合する位置を関心領域の位置として設定する。マッチング演算としては、SAD法、部分空間法、Boostingなどの公知の方法が用いられ、これに用いられる学習データは、関心領域データベース12に予め記憶されてもよい。
　図5に示すように、テンプレート506は、探索範囲504の左上からスタートして、左から右へ且つ上から下へ移動する。ただし、図7に示すように、テンプレート702が、探索範囲701の中心点703からスタートして、渦巻き状に移動してもよい。中心点703は、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標(基準座標)であるので、探索範囲701の中心点703の近傍に関心領域位置が設定される可能性が高い。したがって、テンプレート702を渦巻き状に移動させて関心領域を探索して、テンプレート702と所定の閾値以上の相関性を有する位置を関心領域位置と設定し、探索を終了してもよい。この結果、関心領域設定の探索時間を短縮することができる。
　超音波診断装置1の関心領域設定部8が、関心領域又は特徴部位の設定を終了するか判断する(ステップS106)。関心領域設定部11が、1/4に縮小(圧縮)された画像データに対して、特徴部位を含む関心領域の設定を行う場合、設定分解能も1/4となる。したがって、縮小率(圧縮率)が小さいほど設定分解能が小さくなり、関心領域又は特徴部位の設定位置は粗く設定される。関心領域又は特徴部位の設定を終了するか否かは、縮小率(圧縮率)、設定分解能、設定処理時間、及び設定処理回数のうち少なくとも1つによって判断されればよい。
　ステップS106において、関心領域又は特徴部位の設定を継続すると判断された場合は、探索範囲設定部9が、ステップS105において設定された関心領域(特徴部位)を含むように、探索範囲の設定を行う(ステップS107)。この探索範囲は、上述のように、第1の圧縮率(例えば、データ量1/4＝圧縮率1/4)で圧縮された画像データについて、第1の探索範囲(例えば、初期探索範囲)で関心領域(特徴部位を含む)を設定したあと、第2の圧縮率(例えば、データ量1/2＝圧縮率1/2)で圧縮された画像データについて、第1の探索範囲で設定された関心領域の座標を含む第2の探索範囲を設定する条件の下に決定される。例えば、図5に示すように、圧縮率1/4(第1の圧縮率)で圧縮された画像データ501の初期探索範囲(第1の探索範囲)504で設定された関心領域(特徴部位)の座標を、圧縮率1/2(第2の圧縮率)で圧縮された画像データ502における探索範囲(第2の探索範囲)507の中心点508の座標とする。
　探索範囲507は、中心点508を中心とする所定の画素数の範囲に設定されればよい。また、所定の画素数は、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標値の標準偏差よりも小さくてもよい。その他、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データの一部(例えば、座標データの60％)が含まれるような探索範囲を設定してもよい。この場合、探索範囲(第2の探索範囲)507がカバーする範囲は、探索範囲(第1の探索範囲)504がカバーする範囲よりも小さくてもよい。また、探索範囲(第2の探索範囲)507がカバーする座標データ数(例えば、座標データの60％)は、探索範囲(第1の探索範囲)504がカバーする座標データ数(例えば、座標データの70％)よりも小さくてもよい。なぜなら、探索範囲(第1の探索範囲)504で関心領域又は特徴部位の設定位置が粗く設定されているので、粗く設定された位置である中心点508の近傍に関心領域又は特徴部位の位置が存在する可能性が高いからである。
　次に、ステップS104に戻り、画像縮小部(画像データ圧縮部)10は、複数の画素に基づいて画像データの圧縮を行い、画像データの圧縮率を変化させる。本実施の形態では、図5に示すように、画像縮小部(画像データ圧縮部)10が、超音波画像生成部4により生成される超音波画像(医用画像)の画像データを1/2に圧縮し、画像データ502を生成する。つまり、第2の圧縮率(圧縮率1/2)は、第1の圧縮率(圧縮率1/4)より大きい。
　関心領域設定部11は、第2の探索範囲において関心領域を設定する。例えば、図5に示すように、関心領域設定部11は、1/2に圧縮された画像データ502の探索範囲(第2の探索範囲)507でテンプレート509を移動させて、最も適合する位置を関心領域の位置として設定する(ステップS105)。
　超音波診断装置1の関心領域設定部8が、関心領域又は特徴部位の設定を終了するか判断する(ステップS106)。本実施の形態では、超音波画像生成部4により生成される超音波画像(医用画像)の画像データ(データ量1/1＝圧縮率1/1)になるまで、縮小率(圧縮率)を変化させて関心領域又は特徴部位を設定する。ただし、データ量1/4(圧縮率1/4)又はデータ量1/2(圧縮率1/2)で、十分に要求を満たす程度の位置の正確性が得られる場合は、途中の縮小率(圧縮率)で関心領域又は特徴部位の設定の処理を終えてもよい。圧縮率1/1の画像データまで設定を実施する場合は、さらにステップS107、S104、及びS105を実施することにより、縮小率(圧縮率)を変化させて、関心領域(特徴部位)の位置を設定する。
　例えば、図5に示すように、圧縮率1/2(第1の圧縮率)で圧縮された画像データ502の探索範囲(第1の探索範囲)507で設定された関心領域(特徴部位)の座標を、圧縮率1/1(第2の圧縮率)で圧縮された画像データ503における探索範囲(第2の探索範囲)510の中心点511の座標とし、関心領域設定部11は、圧縮率1/1の画像データ503の探索範囲(第2の探索範囲)511でテンプレート512を移動させて、最も適合する位置を関心領域の位置として設定する(ステップS105)。
　最終的に設定された位置が、図3における特徴部位(弁輪部205、206及び心尖部207)を含む関心領域の位置である。本実施の形態では、画像データの縮小率(圧縮率)を1/4、1/2、及び1/1の3段階にした例について説明したが、任意の縮小率(圧縮率)を設定することが可能である。
　次に、縮小率(圧縮率)と探索画素数の関係について説明する。図8は、本実施の形態の画像データ(例えば128×128画素、256×256画素で構成される)の縮小率と累積探索画素数の関係を示している。図8の縮小率は、1/50、1/40、1/30、1/20、1/10、及び1/1の6段階である。探索画素数は、探索範囲内の画素数であり、縦5画素×横5画素の探索範囲であれば25画素である。
　ここでは、いずれの縮小率でも探索画素数が5×5画素＝25画素であるとして説明する。図8に示すように、探索画素数が25画素の場合、累積探索画素数は、縮小率1/50から1/1へ向かって、25(縮小率1/50)、50(縮小率1/40)、75(縮小率1/30)、100(縮小率1/20)、125(縮小率1/10)、150(縮小率1/1)となる。したがって、縮小率1/50から1/1の画像データまで関心領域又は特徴部位の設定を実施する場合は、マッチング演算回数は合計150となる。
　一方、図4に示すような比較例において、縮小率1/50における5×5画素の探索範囲と同じ範囲の探索範囲で、画像データを縮小(圧縮)せずに(縮小率1/1で)関心領域を設定すると、探索範囲画素数が(5画素×50)×(5画素×50)＝62500となり、62500回のマッチング演算回数が必要となる。したがって、本実施の形態によれば、縮小率(圧縮率)の異なる複数の画像データ上で再帰的に関心領域を設定することで、関心領域設定の探索時間を短縮することができる。
　また、図8の例では縮小率を6段階としたが、任意の段階数(例えば、1/40、1/30、1/20、及び1/10の4段階)とすることもできるつまり、複数生成された圧縮後の画像データの探索範囲を任意に選択し、選択された探索範囲に基づき前記医用画像に関心領域を設定することができる。この場合のマッチング演算回数は、50(縮小率1/40)、75(縮小率1/30)、100(縮小率1/20)、125(縮小率1/10)となる。したがって、縮小率1/40から1/10の画像データまで関心領域又は特徴部位の設定を実施する場合は、マッチング演算回数は合計125となる。このように、要求される設定精度や許容される計算時間に応じて、縮小率の段階数を変更することができる。換言すれば、複数生成された圧縮後の画像データの探索範囲を圧縮比に応じて並べた複数の段階を形成し、形成された段階より任意に段階を選択し、選択された段階の探索範囲に基づき前記医用画像に関心領域を設定することができる。縮小率(圧縮率)の段階数は、図2のステップS106において、関心領域設定部8が、関心領域又は特徴部位の設定を終了するか判断するときに用いられる。
　また、図8に示すように、探索画素数が10×10画素＝100画素の場合、縮小率1/1の画像データまで関心領域又は特徴部位の設定を実施すると、マッチング演算回数は600回となる。また、探索画素数が20×20画素＝400画素の場合、マッチング演算回数は2400回となる。このように、探索時間短縮や正確性などの要求に応じて、探索範囲と縮小率(圧縮率)を任意に選択することが可能である。
　ステップS106において、関心領域又は特徴部位が設定され、設定を終了すると判断された場合は、設定された関心領域が所望の関心領域として設定される(ステップS108)。つまり、計測部13による計測のために、関心領域が設定される。所望の関心領域又は特徴部位が点である場合は、ステップS105で設定された点が関心領域の位置となる。また、関心領域が線や領域である場合は、所定の関心領域生成処理が実施される。例えば、計測部13が心臓の容積を計測する場合は、心臓の輪郭を抽出する関心領域生成処理が実施される。図3に示すように、弁輪部205、206と心尖部207が関心領域(特徴部位)として設定された場合、弁輪部205、206と心尖部207を通るような輪郭線(曲線)208が生成される。曲線208を生成する手法として、多項式曲線や動的輪郭モデルなどがある。
　計測部13が、入力部7により設定された計測項目について計測演算を行う(ステップS109)。本実施の形態では、上述の心臓の輪郭を利用して、Modified Simposon法を適用することによって、計測部13が左室の容積値を算出する。
　出力・表示部6が計測結果を表示する(ステップS110)。例えば、図3に示すように、計測部13により算出された左室容積値が、計測値欄204に表示される。設定された弁輪部205、206、心尖部207、及び輪郭線208が超音波画像上に重畳的に表示されてもよい。
　以上、本発明にかかる第1の実施の形態について説明したが、本発明はこれらに限定されるものではなく、請求項に記載された範囲内において変更・変形することが可能である。
　例えば、図2のステップS103及びS107で、探索範囲が設定される。この場合、基準座標記憶部を含む関心領域データベース12が関心領域の位置の基準となる基準座標を記憶し、探索範囲は、設定された関心領域の座標と基準座標との座標差に基づいて設定されてもよい。そして、探索範囲は、座標差が最小になるように設定されてもよい。
　関心領域又は特徴部位の位置は、過去の複数のサンプルデータにおける関心領域位置(基準座標)の近傍に存在する可能性が高い。しがたって、探索範囲は、過去の複数のサンプルデータに基づく基準座標と設定された関心領域との座標差が最も小さくなるように決定されればよい。
　図9は、縮小率1/40の画像データにおいて探索範囲を決定する例を示した図である。図9に示すように、探索範囲設定部9が、縮小率1/40の画像データで探索範囲設定する。例えば、探索範囲設定部9が、縮小率がより小さい縮小率1/48の画像データで設定された探索範囲の座標値を48/40倍に拡大して、縮小率1/48と同範囲の探索範囲を縮小率1/40の画像データに設定する。関心領域データベース12に記憶されている基準座標(例えば、サンプルデータの座標値の平均)が読み出される。基準座標の位置は、縮小率1/1のときに(x，y)＝(120，120)である場合、縮小率1/40の画像データのときには(x，y)＝(3，3)に相当する。関心領域設定部11が、関心領域(特徴部位)の設定位置と基準座標との座標差(誤差)を計算する。例えば、関心領域の設定位置が(x，y)＝(4，4)であって、基準座標の位置が(x，y)＝(3，3)である場合、座標差はユークリッド距離として算出され、となる。
　関心領域設定部11が、複数のサンプルデータを用いて、サンプルデータごとにサンプルデータの関心領域の位置(基準座標)と設定された関心領域の位置との座標差を計算してもよい。この場合、複数の座標差の平均が算出される。複数の探索範囲に基づいて、座標差が算出され、座標差が最も小さくなるように、探索範囲が設定される。他の縮小率でも同様に探索範囲が設定される。
　また、図2のステップS105で、複数の前記探索範囲が重複する範囲において第1の関心領域が設定された場合、関心領域設定部11は、設定された第1の関心領域の位置に基づいて、探索範囲の一部に限定して第2の関心領域を設定する。例えば、図10に示すように、探索範囲810が心臓の中隔側弁輪部(第1の関心領域又は第1の弁輪部)205を含み、探索範囲811が心臓の側壁側弁輪部(第2の関心領域又は第2の弁輪部)206を含み、中隔側弁輪部(第1の弁輪部)205が側壁側弁輪部(第2の弁輪部)206の左側に位置する場合は、関心領域設定部11は、中隔側弁輪部(第1の弁輪部)205の右側の探索範囲802に限定して側壁側弁輪部(第2の弁輪部)206を設定する。
　通常、中隔側弁輪部205より右側に側壁側弁輪部206が設定されるので、側壁側弁輪部206を探索する場合、中隔側弁輪部205の右側の探索範囲802において、側壁側弁輪部206の設定を実施すればよい。
　一方、側壁側弁輪部(第1の弁輪部)206が中隔側弁輪部(第2の弁輪部)205の右側に位置する場合は、関心領域設定部11は、側壁側弁輪部(第1の弁輪部)206の左側の探索範囲に限定して中隔側弁輪部(第2の弁輪部)を設定する。
　また、本実施の形態を他のアプリケーションへ応用する例として、IMTの計測アプリケーションにおける関心領域の設定について示す。図11に示すように、IMT計測アプリケーションを用いる場合にも、Bモード画像データ上に関心領域が設定される。したがって、関心領域データベース12が、過去のIMT計測(サンプルデータ)で設定された関心領域の座標(基準座標)とIMT計測画像のテンプレートとを記憶し、図2のフローチャートと同様に、関心領域の設定が実施されればよい。
　この結果、IMT計測において、検者により選択された計測項目が入力されると、IMT計測位置に関心領域が設定され、計測値が算出されて表示されるアプリケーションを実現することができる。図11は、頸動脈画像データ1002上で、内膜中膜複合体1004の厚さを計測する関心領域1005が設定される例を示した図である。
　本実施の形態によれば、検者が計測項目を設定すると、計測項目に適した関心領域を設定することができる。関心領域データベース12に格納される基準座標に基づいて初期探索範囲を設定し、縮小率(圧縮率)を変更して関心領域の探索(検索)を実施することにより、関心領域が存在する確率の高い探索範囲で関心領域を設定できる。また、縮小率が異なる複数の画像データを用いて再帰的に関心領域の設定を実施することにより、縮小されない画像データ(縮小率1/1画像データ)と比べ、それぞれの縮小画像データにおいて探索範囲を狭くすることができ、マッチング回数が低減することができるので、高速な設定処理を行うことができる。したがって、圧縮された画像データを用いて探索範囲を設定して、テンプレートマッチングなどによって関心領域を設定することで、関心領域設定の正確性向上と探索時間短縮を達成することができるので、検者の関心領域設定の操作負担を低減する医用画像診断装置を提供できる。
　本発明は、画像撮像部20によって被検体の医用画像を撮像し、画像データ圧縮部10によって撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成し、探索範囲設定部9によって圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定し、関心領域設定部8によって前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する。
　具体的には、本発明の医用画像診断装置は、被検体の医用画像を撮像する画像撮像部20と、撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成する画像データ圧縮部10と、圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定する探索範囲設定部9と、前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する関心領域設定部8と、を備える。
　この構成によれば、圧縮された画像データを用いて探索範囲を設定して、テンプレートマッチングなどによって関心領域を設定することで、リアルタイムで関心領域を設定する場合、関心領域の設定のための探索時間短縮を達成することができる。また、この構成によれば、圧縮率が小さい画像データで関心領域を設定し、圧縮率を徐々に大きくすることで、関心領域設定の正確性を向上することができる。
　本発明の医用画像診断装置では、前記探索範囲設定部9は、前記第2の探索範囲を前記第1の探索範囲に等しい大きさで設定する。
　この構成によれば、第2の探索範囲を第1の探索範囲の大きさが等しいので、各探索範囲での処理を複数のプロセッサに配分して行う場合、前記複数のプロセッサの処理開始を一斉にすれば終了時間が同じであるので、複数のプロセッサ間での処理終了までの待機処理を行う必要がないから、待機処理が不要な分について、リアルタイムで関心領域を設定する場合、関心領域の設定のための探索時間を短縮することができる。
　本発明の医用画像診断装置は、前記関心領域の位置の基準となる基準座標を記憶する基準座標記憶部12を備え、前記探索範囲設定部9は、前記探索範囲を前記準座標に基づいて設定する。
　この構成によれば、基準座標を用いることで、初期探索範囲の位置を設定することができ、基準座標を含むように探索範囲を設定することができる。
　本発明の医用画像診断装置では、前記探索範囲設定部9は、サンプルデータに基づく関心領域の座標値の平均である前記基準座標を中心とし、前記関心領域の前記座標値の標準偏差を前記探索範囲として設定する。
　この構成によれば、基準座標を中心とすることで、探索範囲の位置を設定することができ、サンプルデータに基づく関心領域の座標値の標準偏差を用いることで、探索範囲を設定することができる。
　本発明の医用画像診断装置では、前記探索範囲設定部9は、輪郭抽出処理で設定された特徴部位の座標値である基準座標を前記探索範囲の中心として設定する。
　この構成によれば、輪郭抽出処理により特徴部位が設定されるので、輪郭抽出処理が高精度でなくても、特徴部位の位置関係を保持しつつ、探索範囲を狭く調整することができ、関心領域を設定する際の正確性も向上する。
　本発明の医用画像診断装置は、前記関心領域の位置の基準となる基準座標を記憶する基準座標記憶部12を備え、前記探索範囲設定部は、設定された前記関心領域の座標と前記基準座標との座標差に基づいて前記探索範囲を設定する。
　この構成によれば、関心領域の位置は、基準座標の近傍に存在する可能性が高いので、設定された関心領域の座標と基準座標との座標差に基づいて、探索範囲を設定することで、関心領域を設定する際の正確性が向上する。
　本発明の医用画像診断装置では、前記探索範囲設定部9は、前記座標差が最小になるように前記探索範囲を設定する。
　この構成によれば、基準座標と設定された関心領域との座標差が最も小さくなるように探索範囲が設定されることで、関心領域を設定する際の正確性が向上する。
　本発明の医用画像診断装置では、前記関心領域設定部8は、複数の前記探索範囲が重複する範囲において第1の関心領域が設定された場合、前記第1の関心領域の位置に基づいて、前記探索範囲の一部に限定して第2の関心領域を設定する。
　この構成によれば、探索範囲が重複する場合、探索範囲の一部に限定して関心領域を設定することで、関心領域の設定時間を短縮することができる。
　本発明の医用画像診断装置では、前記関心領域設定部8は、心臓の第1の弁輪部を含む領域を前記第1の関心領域とし、心臓の第2の弁輪部を含む領域を前記第2の関心領域とし、前記第1の弁輪部が前記第2の弁輪部の左側に位置する場合は、前記第1の弁輪部の右側の前記探索範囲に限定して前記第2の弁輪部を設定し、前記第1の弁輪部が前記第2の弁輪部の右側に位置する場合は、前記第1の弁輪部の左側の前記探索範囲に限定して前記第2の弁輪部を設定する。
　この構成によれば、弁輪部は左右にそれぞれ存在するので、一方の弁輪部が設定された場合は、その弁輪部の位置よりも左側又は右側の探索範囲に限定して他方の弁輪部を設定すればよい。この結果、探索範囲を狭くすることができ、関心領域の設定時間を短縮することができる。
　本発明の医用画像診断装置は、設定された前記関心領域において計測項目を計測する計測部13を備え、前記計測部は、前記第1の探索範囲及び前記第2の探索範囲の少なくとも1つにおいて前記関心領域を設定した後に、前記計測項目の計測結果を出力する。
　この構成によれば、第2の探索範囲において関心領域を設定した後に、計測項目の計測結果を出力してもよく、第1の探索範囲及び第2の探索範囲において関心領域を設定した後に、計測項目の計測結果をそれぞれ出力してもよくなり、出力結果に応じて、所望のタイミングで関心領域の設定を停止することができる。
　本発明の医用画像診断装置では、前記関心領域設定部8は、前記計測結果の適正値に基づいて、関心領域の設定を終了する。
　この構成によれば、出力結果の適正値に応じて、所望のタイミングで関心領域の設定を停止することができる。
　本発明の医用画像診断装置の関心領域設定方法は、画像撮像部20によって被検体の医用画像を撮像し、画像データ圧縮部10によって撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成し、探索範囲設定部9によって圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定し、関心領域設定部8によって前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する。ことを特徴とする。
　(第2の実施の形態)
　以下、第2の実施の形態にかかる医用画像診断装置について、図面を用いて説明する。
特に言及しない場合は、他の構成は、第1の実施の形態にかかる医用画像診断装置と同様である。
　本実施の形態では、Bモード及び計測値の画像データ上で関心領域を設定して計測を実施するアプリケーションの例を示す。本実施の形態では、ドプラやMモードの画像データを取得するために、Bモードの画像データ上に関心領域が設定される。また、本実施の形態では、ドプラやMモードの計測項目に関する計測値の画像データに、所定の計測位置を含む関心領域が設定される。さらに、Bモード及び計測値に関する画像データを組み合わせて、複数の関心領域が設定されてもよい。
　画像縮小部(画像データ圧縮部)10が、複数の画素に基づいて画像データの圧縮を行い、画像データの圧縮率(縮小率)を変化させる。探索範囲設定部9が、圧縮された画像データにおいて探索範囲を設定する。関心領域設定部11が、特徴部位を含む関心領域の設定を探索範囲において行う。そして、画像縮小部(画像データ圧縮部)10が、第1の圧縮率で画像データ(Bモード又は計測値の画像データ)を圧縮し、探索範囲設定部9が、第1の圧縮率で圧縮された画像データにおいて、第1の探索範囲を設定し、関心領域設定部11が、第1の探索範囲において関心領域を設定する。さらに、画像縮小部(画像データ圧縮部)10は、第2の圧縮率で画像データ(Bモード又は計測値の画像データ)を圧縮し、探索範囲設定部9は、第2の圧縮率で圧縮された画像データにおいて、第1の探索範囲で設定された関心領域を含む第2の探索範囲を設定し、関心領域設定部11は、第2の探索範囲において関心領域を設定する。
　図12は、本実施の形態に係る医用画像診断装置の動作について説明するフローチャートである。図13は、2方向ドプラ計測の例を示した図である。図14は、Mモードの画像データにおける左室計測の例を示した図である。
　図12に示すように、出力・表示部6の表示画面上に、対象である生体組織が描出された画像が表示される(ステップS201)。対象である生体組織として、図13には心尖部4腔像の画像データが示され、図14には傍胸骨長軸像の画像データが示されている。
　入力部7を用いて計測項目が検者によって設定される(ステップS202)。例えば、図13の2方向ドプラにおいて、流入血流速度や弁輪運動速度を計測項目とする場合、弁口と弁輪にドプラゲート1203を設定する必要がある。また、図14に示すように、Mモードで左室壁や腔の大きさを計測項目とする場合、Bモード画像データ上でMモードビームライン1303を設定する必要がある。
　関心領域設定部8が、ドプラゲート1203やMモードビームライン1303が設定される特徴部位(弁輪及び弁口など)を含む関心領域を設定し、関心領域を設定する。出力・表示部6が、ドプラ画像又はMモード画像を表示する。これらの動作は、ステップS203からステップS207を経て実施され、図2に示すステップS103からステップS107と同様である。ただし、基準座標記憶部を含む関心領域データベース12は、過去の複数のサンプルデータとして取得された弁輪、弁口、及びMモードビームラインの位置(基準座標)やテンプレートを記憶する。関心領域データベース12の基準座標記憶部から基準座標(例えば、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データなど)が読み出され、探索範囲が基準座標に基づいて設定される。
　そして、画像データの縮小率(圧縮率)を変化させながら関心領域を設定し、Bモードの画像データ上に関心領域を設定する(ステップS208)。
　また、本実施の形態では、ドプラやMモードの計測項目に関する画像データに、計測項目の計測位置を含む関心領域を設定する(ステップS303～ステップS308)。
　図13の右上には、Bモード画像に設定された関心領域1207におけるTDI(計測値)の画像データが示され、図13の右下には、パルスドプラ(計測値)の画像データが示されている。また、図14の右には、Bモード画像に設定された関心領域1306のMモードビームライン1303におけるMモードの画像データが表示されている。
　関心領域設定部8が、ドプラやMモードの計測項目に関する計測値の画像データに、所定の計測位置を含む関心領域を設定する。つまり、関心領域設定部8が、画像データにおける関心領域(計測位置)を設定する。図13では、E波、A波、e'波、DT(E波減速時間)の関心領域(計測位置)として、波形位置1206が設定される。図14では、心筋の内外膜の関心領域(計測位置)として、壁位置1308が設定される。これらの動作は、ステップS303からステップS307を経て実施され、ステップS103からステップS107又はステップS203からステップS207と同様である。
　ただし、関心領域データベース12は、ドプラ画像におけるE波、A波、e'波、及びDTの位置(基準座標)やテンプレート1208を記憶し、Mモード画像における内外膜の位置(基準座標)やテンプレート1307を記憶する。関心領域データベース12の基準座標記憶部から基準座標(例えば、過去の複数のサンプルデータに基づく関心領域又は特徴部位の座標データなど)が読み出され、探索範囲が基準座標に基づいて設定される。
　そして、画像データの縮小率(圧縮率)を変化させながら関心領域を設定し、ドプラやMモードの計測項目に関する画像データ上に関心領域(計測位置)を設定する(ステップS308)。本実施の形態では、関心領域(計測位置)が点であるので、関心領域設定部11が設定した関心領域がそのまま計測のための計測位置となる。また、関心領域が線や領域である場合は、所定の関心領域生成処理が実施される。
　次に、計測部13が、設定された関心領域で計測演算を行う(S409)。図13のドプラ計測では、計測部13が、波形位置(計測位置)1206において、ドプラ波形よりE波、A波、e'波、及びDTの値を算出する。また、計測部13は、計測値の組み合わせであるE/e'を算出してもよい。図14のMモード計測では、計測部13が、壁位置(計測位置)1308において、LVEdD(左室拡張末期径)などの心筋、心腔の大きさに関する計測値を算出する。
　出力・表示部6が、計測部13により算出された計測値を表示する(ステップS410)。それぞれ、図13の計測値欄1202、図14の計測値欄1302に表示される。また、出力・表示部6は、設定された関心領域(ドプラゲート1203、Mモードビームライン1303、波形位置1206、壁位置1308)を超音波画像上に重畳表示してもよい。
　本実施の形態によれば、ドプラやMモードの画像データを取得するために、Bモードの画像データ上に関心領域が設定される。また、ドプラやMモードの計測項目に関する計測値の画像データに、所定の計測位置を含む関心領域が設定される。つまり、Bモードの画像データに、ドプラやMモードの画像データを取得するための関心領域を設定した後に、ドプラやMモードの画像データに、計測項目の計測を行うための関心領域を設定することができる。この結果、様々な画像データに複数の関心領域の設定が可能になる。したがって、Bモード画像に関心領域を設定して、検者が所望する位置のドプラ画像データ又はMモード画像データが生成されるとともに、検者により入力された計測項目に応じて、所定の計測位置(関心領域)を設定することで、関心領域設定の正確性向上と探索時間短縮を達成することができるので、検者の関心領域設定の操作負担を低減する医用画像診断装置を提供できる。
　以上、本発明にかかる第1及び第2の実施の形態について説明したが、本発明はこれらに限定されるものではなく、請求項に記載された範囲内において変更・変形することが可能である。
　例えば、第1の実施の形態では、縮小率(圧縮率)を1/1まで変化させたが、画像データの縮小率(圧縮率)が徐々に大きくなる過程で、各縮小率(圧縮率)における計測結果を逐次出力し、この計測結果の出力情報を用いて、計測精度を満足する縮小率(圧縮率)で関心領域の設定処理を終了してもよい。つまり、第1の実施の形態では、計測部13は、第2の探索範囲において関心領域を設定した後に、計測項目の計測結果を出力するが、第1の探索範囲及び第2の探索範囲において関心領域を設定した後に、計測項目の計測結果をそれぞれ出力してもよい。
　図15は、計測結果を逐次出力して、計測結果に応じて関心領域の設定を終了するか否かを判断するフローチャートである。図2のフローチャートと比較すると、設定を終了するか否かを判断するステップS106が、計測結果を表示するステップS110の後に実行される。つまり、計測部13は、第1の探索範囲及び第2の探索範囲において関心領域を設定した後に、計測項目の計測結果をそれぞれ出力する。所定の縮小率で設定された関心領域を用いて、計測部13が計測を実行し、出力・表示部6が計測結果を表示する。この計測結果が適切であれば計測を終了し、計測結果が適切でなければステップS107に戻って、次の縮小率(圧縮率)で関心領域の設定が行われる。計測結果が適切か否かの判断は、検者が計測値を確認して行われてもよいし、入力部7により入力された計測結果の適正値(閾値)と比較することで行われてもよいし、記憶部5に記憶された計測結果の適正値(閾値)と比較することで行われてもよい。また、適正値(閾値)に近づいているか否かを判断して、最も適正値(閾値)に近づいたときに設定処理を終了してもよい。関心領域設定部11は、計測結果の適正値に基づいて、関心領域の設定を終了する。
　この結果、縮小率(圧縮率)が1/1に達する前であっても、適正な計測値が計測されると判断された場合は、関心領域の設定処理を終了することができる。したがって、適正な計測精度を確保しつつ、関心領域設定の正確性向上と探索時間短縮を達成することができる。
　また、生体組織の輪郭を抽出する方法を用いてもよい。輪郭抽出方法を用いて関心領域を設定し、その設定位置の周辺を詳細に探索することで、関心領域の設定精度を高めることができる。例えば、心臓の内膜輪郭の抽出方法を用いて、弁輪及び心尖部(特徴部位)を含む関心領域を設定する。
　図16は、輪郭抽出方法を用いた関心領域の設定を説明するフローチャートである。図16に示すように、第1の実施の形態と比較して、輪郭抽出処理を行うステップS502が加わっている。輪郭抽出処理は、動的輪郭モデルのような形状に制約のある輪郭が望ましい。左室の輪郭形状の特徴を捉えたモデルで制約することによって左室特有の輪郭を抽出することができ、弁輪部と心尖部の正しい位置関係を設定することができる。このモデルが、記憶部5に記憶されて、輪郭抽出時に参照されるようにしてもよい。
　探索範囲設定部9が、輪郭抽出処理を行う(ステップS502)。図17に示すように、輪郭抽出処理を行うステップS502において、探索範囲設定部9が輪郭1603を抽出する。抽出された輪郭1603は、高精度ではなくても、弁輪部1601と心尖部1602の位置関係が保持された状態で、輪郭線が内膜にフィットするように抽出されればよい。
　そして、探索範囲設定部9が探索範囲1604を設定する(ステップS107)。探索範囲設定部9は、輪郭抽出処理で設定された弁輪部1601及び心尖部1602の座標値(基準座標)を中心とした探索範囲1604を設定する。探索範囲1604は、第1の実施の形態と同様に、関心領域データベース12の過去のサンプルデータを用いて設定されてもよい。輪郭抽出処理で弁輪部1601と心尖部1602の位置が設定されているので、探索範囲1604の範囲がサンプルデータの標準偏差よりも狭くなるように設定されてもよい。
　関心領域設定部8が、計測部13による計測のための関心領域を設定する(ステップS108)。ここでは、関心領域は輪郭(線)1603である。ステップS108までの処理で、弁輪部1601と心尖部1602の位置が設定されているので、弁輪部2つと心尖部1つの3点を通るように輪郭1603が抽出される。
　この結果、輪郭抽出処理により弁輪部1601及び心尖部1602の位置が設定されるので、輪郭抽出処理が高精度でなくても、弁輪部1601と心尖部1602の位置関係を保持しつつ、設定された弁輪部1601と心尖部1602の位置を用いれば、探索範囲1604を狭く調整することができ、弁輪部1601及び心尖部1602を設定する際の正確性も向上する。したがって、探索範囲1604を狭く設定することができるので、設定精度の向上とともに設定演算時間の短縮を実現することができる。
　以上、本発明にかかる実施の形態について説明したが、本発明はこれらに限定されるものではなく、請求項に記載された範囲内において変更・変形することが可能である。
　第1及び第2の実施の形態では、超音波診断装置における心臓の計測を例に説明したが、本発明は他の診断装置及び生体組織に対しても同様に適用可能である。
　画像データは、被検体の断層画像、弾性画像、及び計測値画像の少なくとも1つであればよい。また、医用画像には、X線画像、CT、MRI、超音波画像(US)、及び血管造影(血管撮影)などがある。つまり、画像撮像部20は、X線画像撮像装置、X線CT装置、MRI装置などである。
　本発明にかかる医用画像診断装置は、圧縮された画像データを用いて探索範囲を設定して、テンプレートマッチングなどによって関心領域を設定することで、関心領域設定の正確性向上と探索時間短縮を達成することができるので、検者の関心領域設定の操作負担を低減する医用画像診断装置を提供できるという効果を有し、特徴部位又は計測位置を含む関心領域を設定する医用画像診断装置などとして有用である。
　1　超音波診断装置、2　探触子、3　超音波信号生成部、4　超音波画像生成部、5　記憶部、6　出力・表示部、7　入力部、8　関心領域設定部、9　探索範囲設定部、10　画像縮小部(画像データ圧縮部)、11　関心領域設定部、12　関心領域データベース、13　計測部、14　制御部、20　画像撮像部、203　心臓

Claims (15)

  1. 　被検体の医用画像を撮像する画像撮像部と、
    　撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成する画像データ圧縮部と、
    　圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定する探索範囲設定部と、
    　前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する関心領域設定部と、
    　を備えたことを特徴とする医用画像診断装置。
  2. 　前記探索範囲設定部は、前記圧縮前の探索範囲を前記圧縮後の探索範囲に等しい大きさで設定することを特徴とする請求項1に記載の医用画像診断装置。
  3. 　前記探索範囲設定部は、輪郭抽出処理で設定された特徴部位の座標値である基準座標を前記探索範囲の中心として設定することを特徴とする請求項2に記載の医用画像診断装置。
  4. 　前記関心領域の位置の基準となる基準座標を記憶する基準座標記憶部を備え、
    　前記探索範囲設定部は、前記探索範囲を前記基準座標に基づいて設定することを特徴とする請求項1に記載の医用画像診断装置。
  5. 　前記探索範囲設定部は、サンプルデータに基づく関心領域の座標値の平均である前記基準座標を中心とし、前記関心領域の前記座標値の標準偏差を前記探索範囲として設定することを特徴とする請求項3に記載の医用画像診断装置。
  6. 　前記関心領域の位置の基準となる基準座標を記憶する基準座標記憶部を備え、
    　前記探索範囲設定部は、設定された前記関心領域の座標と前記基準座標との座標差に基づいて前記探索範囲を設定することを特徴とする請求項1に記載の医用画像診断装置。
  7. 　前記探索範囲設定部は、前記座標差が最小になるように前記探索範囲を設定することを特徴とする請求項6に記載の医用画像診断装置。
  8. 　前記関心領域設定部は、複数の前記探索範囲が重複する範囲において第1の関心領域が設定された場合、前記第1の関心領域の位置に基づいて、前記探索範囲の一部に限定して第2の関心領域を設定することを特徴とする請求項1に記載の医用画像診断装置。
  9. 　前記関心領域設定部は、
    　心臓の第1の弁輪部を含む領域を前記第1の関心領域とし、
    　心臓の第2の弁輪部を含む領域を前記第2の関心領域とし、
    　前記第1の弁輪部が前記第2の弁輪部の左側に位置する場合は、前記第1の弁輪部の右側の前記探索範囲に限定して前記第2の弁輪部を設定し、前記第1の弁輪部が前記第2の弁輪部の右側に位置する場合は、前記第1の弁輪部の左側の前記探索範囲に限定して前記第2の弁輪部を設定することを特徴とする請求項8に記載の医用画像診断装置。
  10. 　設定された前記関心領域において計測項目を計測する計測部を備え、
    　前記計測部は、前記圧縮後の探索範囲及び前記圧縮前の探索範囲の少なくとも1つにおいて前記関心領域を設定した後に、前記計測項目の計測結果を出力することを特徴とする請求項1に記載の医用画像診断装置。
  11. 　前記関心領域設定部は、前記計測結果の適正値に基づいて、関心領域の設定を終了することを特徴とする請求項10に記載の医用画像診断装置。
  12. 　前記画像データ圧縮部は前記圧縮後の画像データを異なる圧縮比で複数生成し、
    前記探索範囲設定部は、複数生成された圧縮後の画像データにそれぞれ探索範囲を設定し、
    前記関心領域設定部は、複数生成された圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定することを特徴とする請求項1に記載の医用画像診断装置。
  13. 　前記関心領域設定部は、複数生成された圧縮後の画像データの探索範囲を任意に選択し、選択された探索範囲に基づき前記医用画像に関心領域を設定することを特徴とする請求項12に記載の医用画像診断装置。
  14. 　前記関心領域設定部は、複数生成された圧縮後の画像データの探索範囲を圧縮比に応じて並べた複数の段階を形成し、形成された段階より任意に段階を選択し、選択された段階の探索範囲に基づき前記医用画像に関心領域を設定することを特徴とする請求項12に記載の医用画像診断装置。
  15. 画像撮像部によって被検体の医用画像を撮像し、
    画像データ圧縮部によって撮像された被検体の医用画像の画像データを圧縮前の画像データとし、当該圧縮前の画像データの複数の画素に基づいて圧縮して圧縮後の画像データを生成し、
    探索範囲設定部によって圧縮後の画像データの探索範囲を設定すると共に、圧縮前の画像データの探索範囲を設定し、
    関心領域設定部によって前記圧縮前の画像データの探索範囲と前記圧縮後の画像データの探索範囲に基づき前記医用画像に関心領域を設定する。
    ことを特徴とする医用画像診断装置の関心領域設定方法。
PCT/JP2013/051926 2012-02-02 2013-01-29 医用画像診断装置及びその関心領域設定方法 WO2013115194A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201380007856.5A CN104093363B (zh) 2012-02-02 2013-01-29 医学图像诊断装置及其设定感兴趣区域的方法
EP13743588.9A EP2810602A4 (en) 2012-02-02 2013-01-29 DIAGNOSTIC APPARATUS FOR MEDICAL IMAGES AND METHOD FOR SETTING A SPECIFIC REGION THEREFOR
JP2013556419A JP6023091B2 (ja) 2012-02-02 2013-01-29 医用画像診断装置及びその関心領域設定方法
US14/373,837 US9514531B2 (en) 2012-02-02 2013-01-29 Medical image diagnostic device and method for setting region of interest therefor

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012021070 2012-02-02
JP2012-021070 2012-02-02

Publications (1)

Publication Number Publication Date
WO2013115194A1 true WO2013115194A1 (ja) 2013-08-08

Family

ID=48905229

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/051926 WO2013115194A1 (ja) 2012-02-02 2013-01-29 医用画像診断装置及びその関心領域設定方法

Country Status (5)

Country Link
US (1) US9514531B2 (ja)
EP (1) EP2810602A4 (ja)
JP (1) JP6023091B2 (ja)
CN (1) CN104093363B (ja)
WO (1) WO2013115194A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016147046A (ja) * 2015-02-12 2016-08-18 東芝メディカルシステムズ株式会社 医用画像処理装置、医用画像処理方法、医用画像処理システム、及び医用画像処理プログラム
JP2020092739A (ja) * 2018-12-10 2020-06-18 キヤノンメディカルシステムズ株式会社 医用画像生成装置、及び医用画像生成プログラム
CN113081030A (zh) * 2020-01-08 2021-07-09 通用电气精准医疗有限责任公司 用于基于m模式分析进行辅助超声扫描平面识别的方法和系统

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5841034B2 (ja) * 2012-09-27 2016-01-06 富士フイルム株式会社 超音波診断装置、超音波画像生成方法およびプログラム
JP2014198144A (ja) * 2013-03-29 2014-10-23 ソニー株式会社 画像処理装置、画像処理方法、情報処理プログラム、蛍光観察システム、および蛍光ナビゲーション・サージェリー・システム
JP5771297B1 (ja) * 2014-02-24 2015-08-26 日立アロカメディカル株式会社 超音波診断装置
JP6197952B2 (ja) * 2014-05-12 2017-09-20 富士通株式会社 商品情報出力方法、商品情報出力プログラムおよび制御装置
JP5842039B1 (ja) * 2014-08-15 2016-01-13 日立アロカメディカル株式会社 超音波画像処理装置
US20170094171A1 (en) * 2015-09-28 2017-03-30 Google Inc. Integrated Solutions For Smart Imaging
JP7099901B2 (ja) * 2018-08-06 2022-07-12 富士フイルムヘルスケア株式会社 超音波画像処理装置およびプログラム
CN111354006A (zh) * 2018-12-21 2020-06-30 深圳迈瑞生物医疗电子股份有限公司 超声图像中目标组织的描迹方法及装置
KR102670514B1 (ko) * 2019-01-29 2024-05-30 삼성메디슨 주식회사 초음파 진단 장치 및 그 동작방법
CN113160116B (zh) * 2021-02-03 2022-12-27 中南民族大学 左心室内外膜自动分割方法、系统及设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005218796A (ja) 2004-02-09 2005-08-18 Matsushita Electric Ind Co Ltd 医用画像処理装置および医用画像処理方法
JP2005287927A (ja) * 2004-04-02 2005-10-20 Konica Minolta Medical & Graphic Inc 画像処理装置、画像処理方法及び医用画像システム
JP2007037864A (ja) * 2005-08-04 2007-02-15 Hitachi Medical Corp 医用画像処理装置
JP2008000586A (ja) * 2006-05-26 2008-01-10 Toshiba Corp 画像処理装置および磁気共鳴イメージング装置
JP2010136779A (ja) * 2008-12-10 2010-06-24 Univ Of Miyazaki 医用画像の処理装置及び方法
JP2011217130A (ja) * 2010-03-31 2011-10-27 Fujifilm Corp 画像符号化装置および方法並びにプログラム
JP2012021904A (ja) * 2010-07-15 2012-02-02 Sony Corp 画像処理装置および方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1420633A (zh) * 2001-10-29 2003-05-28 日本胜利株式会社 图像编码装置以及图像编码方法
TWI226189B (en) 2003-03-28 2005-01-01 Nat Univ Chung Cheng Method for automatically detecting region of interest in the image
DE112006003459B4 (de) * 2005-12-21 2015-02-05 Konica Minolta, Inc. Ultraschalldiagnosevorrichtung
KR20100095833A (ko) 2009-02-23 2010-09-01 주식회사 몬도시스템즈 Ｒｏｉ 의존형 압축 파라미터를 이용하여 영상을 압축하는 장치 및 방법
US9259199B2 (en) 2011-05-27 2016-02-16 Kabushiki Kaisha Toshiba Image processing apparatus and X-ray diagnosis apparatus
CN103607951B (zh) * 2011-09-28 2016-08-24 达卡拉通信系统株式会社 图像处理装置及图像处理方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005218796A (ja) 2004-02-09 2005-08-18 Matsushita Electric Ind Co Ltd 医用画像処理装置および医用画像処理方法
JP2005287927A (ja) * 2004-04-02 2005-10-20 Konica Minolta Medical & Graphic Inc 画像処理装置、画像処理方法及び医用画像システム
JP2007037864A (ja) * 2005-08-04 2007-02-15 Hitachi Medical Corp 医用画像処理装置
JP2008000586A (ja) * 2006-05-26 2008-01-10 Toshiba Corp 画像処理装置および磁気共鳴イメージング装置
JP2010136779A (ja) * 2008-12-10 2010-06-24 Univ Of Miyazaki 医用画像の処理装置及び方法
JP2011217130A (ja) * 2010-03-31 2011-10-27 Fujifilm Corp 画像符号化装置および方法並びにプログラム
JP2012021904A (ja) * 2010-07-15 2012-02-02 Sony Corp 画像処理装置および方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2810602A4

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016147046A (ja) * 2015-02-12 2016-08-18 東芝メディカルシステムズ株式会社 医用画像処理装置、医用画像処理方法、医用画像処理システム、及び医用画像処理プログラム
JP2020092739A (ja) * 2018-12-10 2020-06-18 キヤノンメディカルシステムズ株式会社 医用画像生成装置、及び医用画像生成プログラム
JP7309355B2 (ja) 2018-12-10 2023-07-18 キヤノンメディカルシステムズ株式会社 医用画像生成装置、及び医用画像生成プログラム
CN113081030A (zh) * 2020-01-08 2021-07-09 通用电气精准医疗有限责任公司 用于基于m模式分析进行辅助超声扫描平面识别的方法和系统

Also Published As

Publication number Publication date
JP6023091B2 (ja) 2016-11-09
CN104093363B (zh) 2017-04-05
US20140369564A1 (en) 2014-12-18
EP2810602A1 (en) 2014-12-10
US9514531B2 (en) 2016-12-06
CN104093363A (zh) 2014-10-08
JPWO2013115194A1 (ja) 2015-05-11
EP2810602A4 (en) 2015-06-17

Similar Documents

Publication Publication Date Title
JP6023091B2 (ja) 医用画像診断装置及びその関心領域設定方法
US11986355B2 (en) 3D ultrasound imaging system
WO2017206023A1 (zh) 一种心脏容积识别分析系统和方法
US8913816B2 (en) Medical image dianostic device, region-of-interest setting method, and medical image processing device
US20090149756A1 (en) Method, apparatus and computer program for three-dimensional ultrasound imaging
US20150272546A1 (en) Ultrasound apparatus and method
US9877698B2 (en) Ultrasonic diagnosis apparatus and ultrasonic image processing apparatus
US10922874B2 (en) Medical imaging apparatus and method of displaying medical image
US10646203B2 (en) Ultrasound diagnosis method and apparatus for analyzing contrast enhanced ultrasound image
US9888905B2 (en) Medical diagnosis apparatus, image processing apparatus, and method for image processing
US10004478B2 (en) Method and apparatus for displaying ultrasound image
US20210045717A1 (en) Methods and system for camera-aided ultrasound scan setup and control
US10398411B2 (en) Automatic alignment of ultrasound volumes
KR20170060853A (ko) 의료 영상 장치 및 그 동작방법
CN115426954A (zh) 用于生成路线图图像的双平面和三维超声图像采集以及相关联的系统和设备
CN106028946B (zh) 用于监测病变尺寸趋势的系统及其操作方法
EP3409210B1 (en) Ultrasound diagnosis apparatus and operating method thereof
KR20190085741A (ko) 초음파 영상 장치 및 그의 제어 방법
KR20210114281A (ko) 초음파 영상 장치, 그 제어 방법, 및 컴퓨터 프로그램
US10321893B2 (en) Method and apparatus for generating ultrasound image
US20200281566A1 (en) Ultrasonic imaging apparatus and method of controlling the same

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13743588

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2013556419

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14373837

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2013743588

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE