US9563978B2 - Image generation apparatus, method, and medium with image generation program recorded thereon - Google Patents

Image generation apparatus, method, and medium with image generation program recorded thereon Download PDF

Info

Publication number
US9563978B2
US9563978B2 US14/483,376 US201414483376A US9563978B2 US 9563978 B2 US9563978 B2 US 9563978B2 US 201414483376 A US201414483376 A US 201414483376A US 9563978 B2 US9563978 B2 US 9563978B2
Authority
US
United States
Prior art keywords
pixel
value
opacity curve
pixels
region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active, expires
Application number
US14/483,376
Other languages
English (en)
Other versions
US20140375631A1 (en
Inventor
Jun Masumoto
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
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 Fujifilm Corp filed Critical Fujifilm Corp
Assigned to FUJIFILM CORPORATION reassignment FUJIFILM CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: MASUMOTO, JUN
Publication of US20140375631A1 publication Critical patent/US20140375631A1/en
Application granted granted Critical
Publication of US9563978B2 publication Critical patent/US9563978B2/en
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—Three-dimensional [3D] image rendering
    • G06T15/08—Volume rendering
    • 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/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Clinical applications
    • A61B8/0891—Clinical applications for diagnosis of blood vessels
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461—Displaying means of special interest
    • A61B8/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4887—Locating particular structures in or on the body
    • A61B5/489—Blood vessels
    • 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/483—Diagnostic techniques involving the acquisition of a 3D volume of data
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Definitions

  • the present invention relates to an image generation apparatus, method, and program that generates a pseudo three-dimensional image by performing volume rendering on a three-dimensional image.
  • processing for generating and displaying a pseudo three-dimensional image has been performed in which three-dimensional image data of a subject obtained by CT equipment, MRI equipment, ultrasonic diagnostic equipment, or the like are stereoscopically visualized on a two-dimensional plane using computer graphics technologies in order to facilitate understanding of a three-dimensional structure of the subject and the like.
  • a volume rendering method is known in which opacity and color information of R, G, B are set to each pixel value (voxel value) constituting a three-dimensional image and the visualization is achieved by performing ray-casting on each pixel on a projection plane from the observation side.
  • opacity curve and color map respectively For the opacity or color information set (referred to as opacity curve and color map respectively) which is set to each pixel value, it is common to give one set to one three-dimensional image. But this causes a problem that different tissues having the same signal value cannot be displayed distinguishably.
  • a literature H. Imai, “I see! ! The Bible of the Way of Thinking and Processing—Three-Dimensional Medical Image”, Shujunsha, P. 105, 2003, proposes that different tissues having the same signal value are represented by different color or opacity by applying a different color map and/or opacity curve to each object region.
  • a method that further visualizes and makes that portion extractable by performing imaging after injecting a special liquid called a contrast agent is sometimes used so that pixel values of the blood vessel differ from those of other organs in a three-dimensional image.
  • a variation in pixel value occurs from place to place even in the same blood vessel region due to the imaging timing or the amount of the contrast agent. Therefore, when generating an image that represents a wide range of blood vessel region, even for one target object tissue of blood vessel, the simple application of one color map and/or an opacity curve causes a problem that a three-dimensional morphology of the entire blood vessel cannot be represented accurately.
  • the problem due to the variation in pixel value may possibly occur by a partial volume effect.
  • a blood vessel which should have actually a high pixel value is imaged with a pixel value lower than the real value influenced by the surrounding pixel values at the time of the imaging.
  • This may result in that, in an image generated by applying one color map and opacity curve which is based on the pixel value distribution of a thick blood vessel, a thick blood vessel portion appears clearly but a thin blood vessel portion is not visualized, as shown, for example, in FIG.
  • Japanese Unexamined Patent Publication No 2011-212219 proposes a method for dynamically correcting color map and/or opacity curve with respect to each image, but the method is based on the assumption that the depiction range of each image is so small that the variation in pixel value is negligible and related to how to decide one opacity curve used for generating on image. Therefore, the method does not solve the aforementioned problem when generating an image representing a wider object range.
  • an object of the present invention to provide an image generation apparatus, method, and program capable of generating an image which depicts three-dimensional morphology of a predetermined target object more accurately when generating a pseudo three-dimensional image by performing volume rendering.
  • An image generation apparatus of the present invention includes an image generation section that generates a pseudo three-dimensional image by performing volume rendering on a three-dimensional image using an opacity curve that defines the relationship between pixel value and opacity, a region identification section that identifies a whole region representing a predetermined target object from the three-dimensional image, and an opacity curve setting section that sets a base opacity curve to the identified whole region, obtains, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned, and sets an opacity curve obtained by modifying the base opacity curve using the obtained representative value as the opacity curve to be applied to the pixel concerned in the volume rendering (first image generation apparatus).
  • the opacity curve setting section may be a section that obtains, with respect to each of the at least some pixels, a shift amount of opacity curve based on a subtraction value obtained by subtracting a representative value of pixel values in the whole region from the representative value of pixels in the adjacent region of the pixel concerned and sets an opacity curve obtained by shifting the base opacity curve by the obtained shift amount in a pixel value direction as the opacity curve to be set to the pixel concerned in the volume rendering.
  • pixel value direction refers to, when the opacity curve is set in a coordinate system with one axis (e.g., horizontal axis) representing the pixel value and the other axis (e.g., vertical axis) representing the opacity, a direction of the axis representing the pixel value.
  • one axis e.g., horizontal axis
  • the other axis e.g., vertical axis
  • shift amount When obtaining the shift amount based on the subtraction value, a value having the same positive or negative sign as that of the subtraction value is obtained as the shift amount.
  • the absolute value of the shift amount may be the same as the absolute value of the subtraction value or a value obtained by multiplying the value with a predetermined coefficient other than 1.
  • shifting by the shift amount refers to shifting to plus direction of the axis representing the pixel value if the sign of the shift amount is plus while shifting to minus direction of the axis representing the pixel value if the sign of the shift amount is minus.
  • the opacity curve setting section may be a section that obtains, with respect to each of some pixels obtained by sampling the pixels in the whole region at a predetermined interval, a representative value of pixels in an adjacent region of a pixel concerned, obtains a shift amount of opacity curve based on a difference between the obtained representative value and a representative value of pixel values in the whole region, and sets an opacity curve obtained by shifting the base opacity curve by the obtained shift amount to a pixel value direction as the opacity curve to be applied to the pixel concerned, and obtains, with respect to each of pixels in the whole region other than the sampled some pixels, an estimated value of shift amount at a pixel concerned by interpolation using the shift amount of each of the two or more sampled pixels located adjacent to the pixel concerned and sets an opacity curve obtained by shifting the base opacity curve by the determined estimated value in a pixel value direction as the opacity curve to be applied to the pixel concerned in
  • the representative value of the pixel values in the adjacent region may be a mode value, a median value, or an average value of the pixel values in the adjacent region, an average value of pixel values of all pixel values in the adjacent region which fall in a range of predetermined width of pixel values from the mode value or the median value of the pixels in the adjacent region, or a mode value or a median value of pixel values of all pixel values in the adjacent region which fall in a range of predetermined width of pixel values from the mode value or the median value of the pixels in the adjacent region.
  • the representative value of the pixels in the whole region may be a value of the same type as that of the representative value of the pixels in the adjacent region or a value of different type from that of the representative value of the pixels in the adjacent region.
  • the representative value of the pixels in the adjacent region is the average value
  • the representative value of the pixels in the whole region may be the average value or the mode value of the pixels in the whole region.
  • the image generation section may be a section that generates a pseudo three-dimensional image by performing volume rendering on the three-dimensional image using a color map that defines the relationship between pixel value and display color
  • the apparatus may further includes a color map setting section that sets a base color map to the identified whole region and calculates, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned, and sets a color map obtained by modifying the set base color map using the calculated representative value as the color map to be applied to the pixel concerned in the volume rendering.
  • An image generation apparatus of the present invention includes an image generation section that generates a pseudo three-dimensional image by performing volume rendering on a three-dimensional image using a color map that defines the relationship between pixel value and display color, a region identification section that identifies a whole region representing a predetermined target object from the three-dimensional image, and a color map setting section that sets a base color map to the identified whole region, calculates, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned, and sets a color map obtained by modifying the set base color map using the calculated representative value as the color map to be applied to the pixel concerned in the volume rendering (second image generation apparatus).
  • the predetermined target object may be a blood vessel.
  • First and second image generation methods of the present invention are methods, each performs the processing performed by each section of the first or the second image generation apparatus with at least one computer.
  • First and second image generation programs of the present invention are programs, each causes at least one computer to perform the processing performed by each section of the first or the second image generation apparatus.
  • the programs are recorded on recording media, such as CD-ROM, DVD, and the like, or recorded on an auxiliary storage of a server computer or a network storage in a downloadable manner and supplied to the user.
  • the first image generation apparatus, method, and program of the present invention when generating a pseudo three-dimensional image by performing volume rendering on a three-dimensional image using an opacity curve that defines the relationship between pixel value and opacity, a whole region representing a predetermined target object is identified from the three-dimensional image, a base opacity curve is set to the identified whole region, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned is obtained, and an opacity curve obtained by modifying the base opacity curve using the obtained representative value is set as the opacity curve to be applied to the pixel concerned in the volume rendering.
  • This effect is more significant if the predetermined target object is a blood vessel and further significant if the three-dimensional image is an image obtained by injecting a contrast agent into the blood vessel.
  • the representative value of the pixel values in the adjacent region is the mode value or the median value of the pixel values in the adjacent region, an extreme value due to noise, if present in the adjacent region, is less likely to be reflected to the representative value.
  • the representative value of the pixel values in the adjacent region is an average value, a mode value, or a median value of pixel values of all pixel values in the adjacent region which fall in a range of predetermined width of pixel values from the mode value or the median value of the pixels in the adjacent region, an extreme value that presents on the maximum or the minimum side of the pixel values in the adjacent region is less likely to be reflected to the representative value.
  • a base color map is set to the identified whole region and, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned is calculated, and a color map obtained by modifying the set base color map using the calculated representative value is set as the color map to be applied to the pixel concerned in the volume rendering, a color map appropriate for depicting the pixel may be applied to each of all pixels in the region representing the predetermined target object, whereby an image more accurately depicting three-dimensional morphology of the predetermined target object may be generated.
  • a whole region representing a predetermined target object is identified from the three-dimensional image, a base color map is set to the identified whole region, with respect to each of at least some pixels in the identified whole region, a representative value of pixel values in an adjacent region of a pixel concerned is calculated, and a color map obtained by modifying the set base color map using the calculated representative value is set as the color map to be applied to the pixel concerned in the volume rendering.
  • This allows a color map appropriate for depicting the pixel to be applied to each of all pixels in the region representing the predetermined target object, whereby an image more accurately depicting three-dimensional morphology of the predetermined target object may be generated.
  • FIG. 1 is a schematic block diagram of an image generation apparatus of a first embodiment.
  • FIG. 2 is a flowchart illustrating an operation of the image generation apparatus of the first embodiment.
  • FIG. 3 is a schematic block diagram of an image generation apparatus of a second embodiment.
  • FIG. 4 is a flowchart illustrating an operation of the image generation apparatus of the second embodiment.
  • FIG. 5 is a drawing illustrating an example image generated by the image generation apparatus of the second embodiment.
  • FIG. 6 is a drawing illustrating an example image generated by conventional technology.
  • FIG. 7 is a drawing illustrating an example image generated by conventional technology.
  • FIG. 1 is a block diagram of an image generation apparatus 1 , illustrating a schematic configuration thereof.
  • the configuration of the image generation apparatus 1 illustrated in FIG. 1 is realized by executing an image generation program read into an auxiliary storage unit on a computer.
  • the image generation program is recorded on a recording medium, such as CD-ROM or the like, or distributed via a network, such as the Internet or the like, and installed on the computer.
  • the image generation program defines image obtaining processing, image generation processing, region specification processing, opacity curve setting processing, and display control processing, as the processing to be performed by CPU of the computer.
  • the computer When each processing described above is executed by the CPU according to the definition described above, the computer functions as an image obtaining section 11 , an image generation section 12 , a region identification section 13 , an opacity curve setting section 14 , and a display control section 16 , which will be described later.
  • the image generation apparatus 1 is connected to a storage unit 2 , such as a hard disk drive or the like, and a display unit 3 , such as a display or the like.
  • the storage unit 2 stores three-dimensional image data (three-dimensional images) obtained by imaging predetermined target objects with imaging equipment, such as CT, MRI, PET, SPECT, ultrasonic imaging equipment, and the like.
  • a three-dimensional image is a collection of pixel data in a three-dimensional space and may be obtained by stacking a plurality of tomographic images obtained by imaging equipment.
  • the image obtaining section 11 obtains the three-dimensional image stored in the storage unit 2 and stores the image in a storage device built in the image generation apparatus 1 or a storage unit connected to the image generation apparatus, such as a hard disk drive.
  • the region identification section 13 identifies a whole region representing a predetermined target object from the three-dimensional image obtained by the image obtaining section 11 .
  • a whole region representing a predetermined target object refers to a region that includes all regions representing the predetermined target object and does not include other regions.
  • the predetermined target object is a blood vessel.
  • the region identification section 13 extracts a region representing the blood vessel (blood vessel region) from the three-dimensional image and identifies the extracted region as a whole region.
  • the extraction of the blood vessel region may be obtained by threshold method, Region Growing method, Level Set method, and other various types of image processing.
  • the region identification section 13 searches a linear structure with respect of each local region in the three dimensional image by calculating eigenvalues of a 3 ⁇ 3 Hessian matrix.
  • a region where a linear structure is included one of the three eigenvalues is close to 0 and the other two are relatively large.
  • the eigenvector corresponding to the eigenvalue close to 0 indicates the principal axis direction of the linear structure.
  • the region identification section 13 makes use of this relationship and determines the likelihood of linear structure with respect to each local region based on eigenvalues of the Hessian matrix. For a local region where a linear structure is identified, the center point thereof is detected as a candidate point. Then, the candidate points detected by the search are connected based on a given algorithm.
  • This builds a tree structure formed of the candidate points and blood vessel branches (edges) connecting the candidate points.
  • the coordinate information of the plurality of detected candidate points and vector information indicating the directions of the blood vessel branches are stored in a memory with the identifiers of the candidate points and blood vessel branches.
  • a blood vessel contour blood vessel external wall
  • the shape may be identified by known segmentation methods as represented by the Graph-Cuts.
  • the blood vessel region which is a tubular structure is extracted, and information necessary to identify the extracted blood vessel region is generated and stored in a memory through the aforementioned processing.
  • the opacity curve setting section 14 sets an opacity curve to be applied, in volume rendering, to each pixel in the blood vessel region identified by the region identification section 13 .
  • the opacity curve defines the relationship between pixel value and opacity and may be represented by a function of pixel value. More specifically, a base opacity curve is set to the blood vessel region first. For example, a pixel value distribution of the blood vessel region and a pixel value distribution of the neighboring region are checked with the three-dimensional image, and pixel values serving as the boundary separating the blood vessel region from other regions are obtained based on these distributions.
  • an opacity curve in which opacity “0” changes to a value other than “0” or opacity “1” changes to a value other than “1” is set adjacent to the pixel values.
  • the opacity curve may be a curve in which the opacity changes in stepwise from “0” to “1” or a curve in which the opacity changes at a predetermined slope according to increase or decrease in the pixel value.
  • the opacity curve setting section 14 obtains, with respect to each of some pixels obtained by three-dimensionally sampling the pixels in the blood vessel region at a predetermined interval, a subtraction value obtained by subtracting the average value of the pixel values in the blood vessel region from the average value of the pixel values of an adjacent region (e.g., region in the range of 1 cm in up, down, left, right, front, and back directions) of a pixel concerned as a shift amount of the opacity curve, and sets an opacity curve obtained by shifting the base opacity curve by the determined shift amount in a pixel value direction as the opacity curve to be applied to the pixel concerned.
  • an adjacent region e.g., region in the range of 1 cm in up, down, left, right, front, and back directions
  • the opacity curve setting section 14 obtains, with respect to each of pixels in the blood vessel region other than the sampled some pixels, an estimated value of shift amount at a pixel concerned by interpolation using the shift amount of each of the two or more sampled pixels located adjacent to the pixel concerned, and sets an opacity curve obtained by shifting the base opacity curve by the determined estimated value in a pixel value direction as the opacity curve to be applied to the pixel concerned.
  • the image generation section 12 generates a pseudo three-dimensional image by performing volume rendering on the three-dimensional image using the opacity curve set by the opacity curve setting section 14 .
  • the display control section 16 displays the pseudo three-dimensional image generated by the image generation section 12 on the display unit 3 .
  • FIG. 2 is a flowchart illustrating an operation of the image generation apparatus 1 .
  • the image obtaining section 11 obtains a three-dimensional image stored in the storage unit 2 (S 1 ).
  • the region identification section 13 extracts a region representing a blood vessel from the three-dimensional image obtained by the image obtaining section 11 and identifies the extracted region as a whole region (S 2 ).
  • the opacity curve setting section 14 sets a base opacity curve to the blood vessel region (S 3 ).
  • the opacity curve setting section 14 obtains, with respect to each of some pixels obtained by three-dimensionally sampling the pixels in the blood vessel region at a predetermined interval, a subtraction value obtained by subtracting the average value of the pixel values in the blood vessel region from the average value of the pixel values of the adjacent region of a pixel concerned as a shift amount of the opacity curve and obtains, with respect to each of pixels in the blood vessel region other than the sampled some pixels, an estimated value of the shift amount at a pixel concerned by interpolation using the shift amount of each of the two or more sampled pixels located adjacent to the pixel concerned (S 4 ).
  • the opacity curve setting section 14 sets, with respect to each of the sampled some pixels, an opacity curve obtained by shifting the base opacity curve by the shift amount determined in step S 4 in a pixel value direction as the opacity curve to be applied to the pixel concerned and sets, with respect to each of pixels other than the sampled some pixels, an opacity curve obtained by shifting the base opacity curve by the estimated value determined in step S 4 in a pixel value direction as the opacity curve to be applied to the pixel concerned (S 5 ).
  • the image generation section 12 generate a pseudo three-dimensional image by performing volume rendering on the three-dimensional image using the opacity curve set by the opacity curve setting section 14 (S 6 ). Then, the display control section 16 displays the pseudo three-dimensional image generated by the image generation section 12 on the display unit 3 (S 7 ), and the processing is completed.
  • the region identification section 13 identifies the whole region representing a blood vessel from a three-dimensional image
  • the opacity curve setting section 14 sets a base opacity curve to the identified whole region, and obtains, with respect to each of at least some of the pixels in the whole region, a representative value of pixel values in an adjacent region of a pixel concerned and sets an opacity curve obtained by changing the base opacity curve using the obtained representative value as the opacity curve to be applied to the pixel concerned in volume rendering
  • the image generation section 12 generates a pseudo three-dimensional image by performing volume rendering using the opacity curves set by the opacity curve setting section 14 . This allows an opacity curve appropriate for depicting the pixel to be applied to each of all pixels in the region representing the blood vessel, whereby an image more accurately depicting three-dimensional morphology of the blood vessel may be generated.
  • FIG. 3 illustrates an image generation apparatus according to the second embodiment of the present invention.
  • the image generation apparatus 100 of the present embodiment includes a color map setting section 15 in addition to the configuration of the image generation apparatus 1 according to the first embodiment illustrated in FIG. 1 .
  • Other aspects are identical to those of the first embodiment.
  • the color map setting section 15 sets a color map to be applied to each pixel in the blood vessel region identified by the region identification section 13 in volume rendering.
  • the color map defines the relationship between the pixel value and display color and may be represented by a function of pixel value. More specifically, a base color map is first set to the blood vessel region.
  • the color map setting section 15 obtains, with respect to each of some pixels obtained by three-dimensionally sampling the pixels in the blood vessel region at a predetermined interval, a subtraction value obtained by subtracting the average value of the pixel values in the blood vessel region from the average value of the pixel values of an adjacent region (e.g., region in the range of 1 cm in up, down, left, right, front, and back directions) of a pixel concerned as a shift amount of the color map, and sets a color map obtained by shifting the base color map by the determined shift amount in a pixel value direction to the pixel concerned as the color map to be applied to the pixel concerned.
  • a subtraction value obtained by subtracting the average value of the pixel values in the blood vessel region from the average value of the pixel values of an adjacent region (e.g., region in the range of 1 cm in up, down, left, right, front, and back directions) of a pixel concerned as a shift amount of the color map
  • the color map setting section 15 obtains, with respect to each of pixels other than the sampled some pixels in the blood vessel region, an estimated value of the shift amount at a pixel concerned by interpolation using the shift amount of each of the two or more sampled pixels located adjacent to the pixel concerned, and sets color map obtained by shifting the base color map by the determined estimated value in a pixel value direction as the color map to be applied to the pixel concerned.
  • the color map setting section 15 set a color map to be applied to each pixel in the blood vessel region using the shift amount and the estimated value obtained by the opacity curve setting section 14 without calculating these values. It should, of course, be appreciated that an arrangement may be made in which the opacity curve setting section 14 sets the opacity curve to be applied to each pixel in the blood vessel region using the shift amount and the estimated value obtained by the color map setting section 15 .
  • the image generation section 12 generates a pseudo three-dimensional image by performing volume rendering on the three-dimensional image using the opacity curve set by the opacity curve setting section 14 and the color map set by the color map setting section 15 .
  • FIG. 5 shows an example pseudo three-dimensional image generated by the image generation apparatus 100 of the present embodiment. In FIG. 5 , the three-dimensional morphology of the entire blood vessel is depicted more accurately in comparison with the images generated by performing volume rendering on the same three-dimensional image according to conventional technology shown in FIGS. 6 and 7 .
  • FIG. 4 is a flowchart illustrating an operation of the image generation apparatus 100 .
  • the image obtaining section 11 obtains a three-dimensional image stored in the storage unit 2 (S 11 ).
  • the region identification section 13 extracts a region representing a blood vessel from the three-dimensional image obtained by the image obtaining section 11 and identifies the extracted region as a whole region (S 12 ).
  • the opacity curve setting section 14 sets a base opacity curve to the blood vessel region and the color map setting section 15 sets a base color map to the blood vessel region (S 13 ).
  • the opacity curve setting section 14 or the color map setting section 15 obtains, with respect to each of some pixels obtained by three-dimensionally sampling the pixels in the blood vessel region at a predetermined interval, a subtraction value obtained by subtracting the average value of the pixel values in the blood vessel region from the average value of the pixel values of the adjacent region of a pixel concerned as a shift amount and obtains, with respect to each of pixels other than the sampled some pixels, an estimated value of the shift amount at a pixel concerned by interpolation using the shift amount of each of the two or more sampled pixels located adjacent to the pixel concerned (S 14 ).
  • the opacity curve setting section 14 sets, with respect to each of the sampled some pixels, an opacity curve obtained by shifting the base opacity curve by the shift amount determined in step S 14 in a pixel value direction as the opacity curve to be applied to the pixel concerned and sets, with respect to each of pixels other than the sampled some pixels, an opacity curve obtained by shifting the base opacity curve by the estimated value determined in step S 14 in a pixel value direction as the opacity curve to be applied to the pixel concerned.
  • the color map setting section 15 sets, with respect to each of the sampled some pixels, a color map obtained by shifting the base color map by the shift amount determined in step S 14 in a pixel value direction as the color map to be applied to the pixel concerned and sets, with respect to each of pixels other than the sampled some pixels, a color map obtained by shifting the base color map by the estimated value determined in step S 14 in a pixel value direction as the color map to be applied to the pixel concerned (S 15 ).
  • the image generation section 12 generate a pseudo three-dimensional image by performing volume rendering on the three-dimensional image using the opacity curves set by the opacity curve setting section 14 and the color maps set by the color map setting section 15 (S 16 ). Then, the display control section 16 displays the pseudo three-dimensional image generated by the image generation section 12 on the display unit 3 (S 17 ), and the processing is completed.
  • the region identification section 13 identifies the whole region representing a blood vessel from a three-dimensional image
  • the opacity curve setting section 14 sets a base opacity curve to the identified whole region, and obtains, with respect to each of at least some of the pixels in the whole region, a representative value of pixel values in an adjacent region of a pixel concerned and sets an opacity curve obtained by changing the base opacity curve using the obtained representative value as the opacity curve to be applied to the pixel concerned in volume rendering
  • the color map setting section 15 sets a base color map to the identified whole region, obtains, with respect to each of at least some of the pixels in the whole region, a representative value of pixel values in an adjacent region of a pixel concerned, and sets an color map obtained by changing the base color map using the obtained representative value as the color map to be applied to the pixel concerned in volume rendering
  • the image generation section 12 generates a pseudo three-dimensional image by performing volume rendering using these opacity curve and color map. This allows an o
  • the opacity curve setting section 14 may be a section that obtains, with respect to each of partial regions having a predetermined size and constituting the whole region identified by the region identification section 13 , a representative value of pixel values in partial region concerned and sets an opacity curve obtained by modifying the base opacity curve using the obtained representative value as the opacity curve to be applied to the partial region concerned in volume rendering, or a section that obtains, with respect to each of partial regions having a predetermined size and constituting the whole region identified by the region identification section 13 , a representative value of pixel values in an adjacent region that includes a partial region concerned and sets an opacity curve obtained by modifying the set base opacity curve using the obtained representative value as the opacity curve to be applied to the partial region concerned in the volume rendering.
  • the color map setting section 15 may be a section that obtains, with respect to each of partial regions having a predetermined size and constituting the whole region identified by the region identification section 13 , a representative value of pixel values in partial region concerned and sets a color map obtained by modifying the base color map using the obtained representative value as the color map to be applied to the partial region concerned in volume rendering, or a section that obtains, with respect to each of partial regions having a predetermined size and constituting the whole region identified by the region identification section 13 , a representative value of pixel values in an adjacent region that includes a partial region concerned and sets color map obtained by modifying the set base color map using the obtained representative value as the color map to be applied to the partial region concerned in the volume rendering.
  • the predetermined target object is a blood vessel
  • the predetermined target object may be a structure having tubular structure, such as an intestine, a bronchus, and the like, or various organs, such as a heart, a liver, and the like.
  • each of the representative value of pixels in the adjacent region and the representative value of pixels in the whole region is an average value
  • each of the representative value of pixels in the adjacent region and the representative value of pixels in the whole region may be a mode or median value of the pixels in each region, an average value of pixels of all the pixels in each region which fall in a range of predetermined width of pixel values from the mode or median value of each region, or a mode or median value of pixels of all the pixels in each region which fall in a range of predetermined width of pixel values from the mode or median value of each region.
  • the representative value of the pixels in the whole region may be a value of the same type as that of the representative value of the pixels in the adjacent region or a value of different type from that of the representative value of the pixels in the adjacent region.
  • the surface of the predetermined target object may also be depicted by the volume rendering.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Graphics (AREA)
  • Vascular Medicine (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Human Computer Interaction (AREA)
  • General Engineering & Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Generation (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
US14/483,376 2012-03-14 2014-09-11 Image generation apparatus, method, and medium with image generation program recorded thereon Active 2033-07-26 US9563978B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012056977A JP5814838B2 (ja) 2012-03-14 2012-03-14 画像生成装置、方法及びプログラム
JP2012-056977 2012-03-14
PCT/JP2013/001634 WO2013136783A1 (ja) 2012-03-14 2013-03-13 画像生成装置、方法及びプログラム

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/001634 Continuation WO2013136783A1 (ja) 2012-03-14 2013-03-13 画像生成装置、方法及びプログラム

Publications (2)

Publication Number Publication Date
US20140375631A1 US20140375631A1 (en) 2014-12-25
US9563978B2 true US9563978B2 (en) 2017-02-07

Family

ID=49160711

Family Applications (1)

Application Number Title Priority Date Filing Date
US14/483,376 Active 2033-07-26 US9563978B2 (en) 2012-03-14 2014-09-11 Image generation apparatus, method, and medium with image generation program recorded thereon

Country Status (6)

Country Link
US (1) US9563978B2 (https=)
EP (1) EP2827301B1 (https=)
JP (1) JP5814838B2 (https=)
CN (1) CN104205172A (https=)
IN (1) IN2014DN07725A (https=)
WO (1) WO2013136783A1 (https=)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6543099B2 (ja) * 2015-06-01 2019-07-10 シーフォーシー・メディカル・ソフトウェア，インコーポレーテッド 情報処理装置、情報処理方法、およびプログラム
JP6521250B2 (ja) * 2015-09-25 2019-05-29 京セラドキュメントソリューションズ株式会社 画像形成装置、色変換プログラムおよび色変換方法
JP6671946B2 (ja) * 2015-12-11 2020-03-25 キヤノン株式会社 情報取得装置、撮像装置及び情報取得方法
US10467798B2 (en) * 2016-12-19 2019-11-05 Canon Medical Systems Corporation Rendering a global illumination image from a volumetric medical imaging data set
JP7013849B2 (ja) * 2017-12-22 2022-02-01 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
JP2019114034A (ja) * 2017-12-22 2019-07-11 大日本印刷株式会社 コンピュータプログラム、画像処理装置、画像処理方法及びボクセルデータ
JP7003635B2 (ja) * 2017-12-22 2022-01-20 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
JP6436258B1 (ja) * 2018-03-27 2018-12-12 大日本印刷株式会社 コンピュータプログラム、画像処理装置及び画像処理方法
JP7206617B2 (ja) * 2018-04-12 2023-01-18 大日本印刷株式会社 断層画像表示用カラーマップの最適化装置およびボリュームレンダリング装置
CN110069310B (zh) * 2019-04-23 2022-04-22 北京小米移动软件有限公司 切换桌面壁纸的方法、装置及存储介质
CN110211216B (zh) * 2019-06-14 2020-11-03 北京理工大学 一种基于体绘制不透明度加权的三维图像空域融合方法
US11443476B2 (en) * 2020-05-27 2022-09-13 Canon Medical Systems Corporation Image data processing method and apparatus
JP7653247B2 (ja) * 2020-12-07 2025-03-28 富士フイルム株式会社 磁気共鳴イメージング装置、画像処理装置および画像処理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004283373A (ja) 2003-03-20 2004-10-14 Toshiba Corp 管腔状構造体の解析処理装置
JP2008086658A (ja) 2006-10-04 2008-04-17 Fujifilm Corp 画像表示装置および画像表示プログラム
US20080278488A1 (en) 2007-05-07 2008-11-13 General Electric Company Methods and apparatus for volume rendering
EP2372661A2 (en) 2010-03-31 2011-10-05 Fujifilm Corporation Projection image generation method, apparatus, and program

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5576117B2 (ja) * 2006-07-31 2014-08-20 コーニンクレッカ フィリップス エヌ ヴェ 画像データセットの視覚化のためのプリセットマップを生成する方法、装置及びコンピュータ可読媒体

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004283373A (ja) 2003-03-20 2004-10-14 Toshiba Corp 管腔状構造体の解析処理装置
US20040249270A1 (en) 2003-03-20 2004-12-09 Kabushiki Kaisha Toshiba Processor for analyzing tubular structure such as blood vessels
JP2008086658A (ja) 2006-10-04 2008-04-17 Fujifilm Corp 画像表示装置および画像表示プログラム
US20080278488A1 (en) 2007-05-07 2008-11-13 General Electric Company Methods and apparatus for volume rendering
EP2372661A2 (en) 2010-03-31 2011-10-05 Fujifilm Corporation Projection image generation method, apparatus, and program
CN102208118A (zh) 2010-03-31 2011-10-05 富士胶片株式会社 投影图像生成方法、设备及程序
US20110242097A1 (en) 2010-03-31 2011-10-06 Fujifilm Corporation Projection image generation method, apparatus, and program
JP2011212219A (ja) 2010-03-31 2011-10-27 Fujifilm Corp 投影画像生成装置、方法、及びプログラム

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Chinese Office Action dated May 3, 2016 with an English translation thereof.
Extended European Search Report dated Oct. 22, 2015.
H. Imai, "I see!! The Bible of the Way of Thinking and Processing-Three-Dimensional Medical Image", Shujunsha, P. 105, 2003 and an English translation thereof.
International Search Report (ISR) (PCT Form PCT/ISA/210), in PCT/JP2013/001634, dated Jun. 4, 2013.
Japanese Office Action dated Jun. 2, 2015 with an English translation thereof.

Also Published As

Publication number Publication date
JP5814838B2 (ja) 2015-11-17
US20140375631A1 (en) 2014-12-25
EP2827301A1 (en) 2015-01-21
EP2827301A4 (en) 2015-11-25
JP2013191030A (ja) 2013-09-26
EP2827301B1 (en) 2019-08-21
WO2013136783A1 (ja) 2013-09-19
CN104205172A (zh) 2014-12-10
IN2014DN07725A (https=) 2015-06-26

Similar Documents

Publication Publication Date Title
US9563978B2 (en) Image generation apparatus, method, and medium with image generation program recorded thereon
CN110010249B (zh) 基于视频叠加的增强现实手术导航方法、系统及电子设备
US7990379B2 (en) System and method for coronary segmentation and visualization
CN102208118B (zh) 投影图像生成方法和设备
USRE35798E (en) Three-dimensional image processing apparatus
US8423124B2 (en) Method and system for spine visualization in 3D medical images
JP6059261B2 (ja) マルチモーダル画像統合におけるレジストレーション精度を改善するインテリジェントなランドマーク選択
US8805045B2 (en) Interaction method for regions of-interest in time series images
US8483462B2 (en) Object centric data reformation with application to rib visualization
US20060171585A1 (en) Method for visualization of plaque deposits from 3D image data records of vessel structures
US10198875B2 (en) Mapping image display control device, method, and program
US9619938B2 (en) Virtual endoscopic image display apparatus, method and program
US20170360396A1 (en) Ultrasound imaging apparatus and method for segmenting anatomical objects
JP4212564B2 (ja) 画像処理方法および画像処理プログラム
JP2019010382A (ja) 画像位置合わせ装置、方法およびプログラム
CN116051553B (zh) 一种在三维医学模型内部进行标记的方法和装置
JP3692050B2 (ja) 画像処理装置及びその方法
JP6702902B2 (ja) マッピング画像表示制御装置および方法並びにプログラム
Kiraly et al. 3D human airway segmentation for virtual bronchoscopy
JP2012075806A (ja) 医用画像処理装置、及び医用画像処理プログラム
WO2015111308A1 (ja) ３次元医用画像表示制御装置およびその作動方法並びに３次元医用画像表示制御プログラム
US9536314B2 (en) Image reconstruction
EP3607527B1 (en) Quantitative evaluation of time-varying data
JP6396597B2 (ja) マッピング画像表示制御装置および方法並びにプログラム
WO2013136750A1 (ja) 画像処理装置および方法並びにプログラム

Legal Events

Date Code Title Description
AS Assignment

Owner name: FUJIFILM CORPORATION, JAPAN

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:MASUMOTO, JUN;REEL/FRAME:033722/0126

Effective date: 20140611

STCF Information on status: patent grant

Free format text: PATENTED CASE

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 4

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 8