This application is based on and claims priority from Japanese Patent Application No. 2008-311192, filed on Dec. 5, 2008, the entire contents of which are hereby incorporated by reference.
1. Technical Field
The present disclosure relates to a medical image processing apparatus and method. More particularly, the present disclosure relates to a medical image processing apparatus and method for generating a cylindrical projection image of a tubular tissue with no distortion in the aspect ratio.
2. Related Art
In recent years, attention has been focused on a technology of visualizing the inside of a three-dimensional object according to the image processing technology using a computer. Particularly, medical diagnosis using an image of Computed Tomography (CT) or a Magnetic Resonance Imaging (MRI), which makes it possible to visualize the inside of a living body, has been widely conducted in a medical field to find a lesion at an early stage.
Also, volume rendering has been used for obtaining a three-dimensional image of the inside of an object. In the volume rendering (typically raycast method), a virtual ray is projected onto three-dimensional voxels (minute volume elements) constituting volume data. Thus, an image is projected onto a projection plane and volume data are visualized. Particularly, in the ray casting method, sampling is performed at given intervals along the path of the virtual ray, and the voxel value is acquired from the voxel at each sampling point. Then, reflected light at each sampling point is stored, and thus volume data are visualized. Some other volume rendering method, for example, a Maximum Intensity Projection (MIP) method, in which the maximum value of the voxels on the virtual ray is acquired to visualize voxel data, is used.
The voxel is a unit of a three-dimensional region of an object and the voxel value is unique and representing the characteristics of the voxel, such as the density value of the voxel. The voxel value is a scalar value in the CT apparatus, but may be a vector value containing color information. The whole object is represented by the voxel data which are a three-dimensional array of the voxel values. Usually, two-dimensional tomographic image data are acquired by the computed tomography (CT) apparatus. Then, the respective two-dimensional tomographic image data are stacked in a direction perpendicular to the tomographic plane and necessary interpolation is performed. Thus, voxel data of the three-dimensional array are obtained.
In the ray casting method, a virtual ray is applied from a virtual eye point to an object, and a virtual reflected light reflected by the object is produced in response to the opacity value artificially set for the voxel value. To capture a virtual surface, the gradient of voxel data, i.e., a normal vector is found and a shading coefficient is calculated from the cosine of the angle between the virtual ray and the normal vector. The virtual reflected light is calculated by multiplying the strength of the virtual ray applied to the voxel by the opacity value of the voxel and the shading coefficient. Also, artificially-setup color may be added to the voxel value.
In visualizing the tubular tissue in the inside of a living body by volume rendering, a parallel projection method or a perspective projection method can be employed. In the parallel projection method, a virtual ray is projected in parallel from a virtual eye point, and thus it is appropriate for observing the tubular tissue from the outside. On the other hand, in the perspective projection method, a virtual ray is projected radially from a virtual eye point, and thus it is appropriate for observing the tubular tissue from the inside thereof. Thus, in the perspective projection method, the endoscopy of the tubular tissue can be simulated. However, to observe the tubular tissue while moving in the inside of the tubular tissue, it is hard to precisely grasp the position and the size of a polyp in the tube wall.
Meanwhile, in visualizing the tubular tissue in the inside of a living body by the volume rendering, a virtual ray is projected radially from the central path of the tubular tissue, whereby an image can be created as if an cylindrical projection image of the tubular tissue were created using a cylindrical coordinate system. This is so-called a cylindrical projection image. In this cylindrical projection image, the position of a polyp in the tube wall, and the size and the shape of the polyp can be observed with one image. In addition, a curved cylindrical projection image provided by performing cylindrical projection onto a winding tubular tissue with the curved central path is also a kind of cylindrical projection image.
As shown in
As described above, the projection plane 85 defined by the central path 84 is represented by the cylindrical coordinates C (h, α) and cylindrical projection is performed from the central path 84. Thus, a 360-degree panoramic image of the tube wall 83 of the tubular tissue can be created.
By the way, when creating a cylindrical projection image of the tube wall 83 of the tubular tissue, spacing of the virtual rays 92 along the central path 84 is constant. Meanwhile, spacing in the circumferential direction perpendicular to the virtual ray 92 projected from a certain position on the central path 84 is constant with respect to the projection plane 85, but is not constant with respect to the tube wall 83 of the tubular tissue. This is because the diameter of the tubular tissue is not always constant.
Accordingly, the projection spacing of the virtual ray 92 in the circumferential direction perpendicular to the central path 84 and the projection spacing of the virtual ray 92 in the direction along the central path 84 will be now discussed with reference to
As shown in
As described above, with change in the projection spacing of the virtual rays, for example, if an object exists on a tube wall of a tubular tissue, the shape of the object largely changes on the cylindrical projection image of the tube wall of the tubular tissue depending on where the object exists.
The cylindrical projection image of an object existing on the tube wall 83 of the tubular tissue 80 shown in
As shown in
In
Namely, on the cylindrical projection image of the tube wall 83 of the tubular tissue 80 shown in
As described above, if the tubular tissue having different diameters is simply displayed on a cylindrical projection image, the aspect ratio of an object is distorted depending on the position of the object existing on the tubular tissue. Thus, for example, if a polyp (i.e., object) exists on the tube wall of a colon (i.e., on a tubular tissue of a living body), the aspect ratio of the polyp is distorted on the cylindrical projection image. Consequently, for example, it becomes hard to distinguish between the polyp and the tissue of the wall of the colon, which leads to obstruction of image diagnosis.
Exemplary embodiments of the present invention are in relation to the above disadvantages and other disadvantages not described above. However, the present invention is not required to overcome the disadvantages described above, and thus, an exemplary embodiment of the present invention may not overcome any of the problems described above.
It is an illustrative aspect of the present invention to provide a medical image processing apparatus and a method capable of generating a cylindrical projection image of a tubular tissue with no distortion in the aspect ratio even if the diameter of the tubular tissue is changed.
According to one or more illustrative aspects of the present invention, there is provided a medical image processing apparatus for visualizing a tubular tissue contained in volume data. The apparatus comprises: a central path determination section that determines a central path of the tubular tissue; a diameter determination section that determines a diameter of the tubular tissue at a certain point on the central path; a spacing determination section that determines at least two or more different projection spacings for projecting virtual rays along the central path, depending on the diameter of the tubular tissue; a cylindrical projection section that projects the virtual rays along the central path with the projection spacing that depends on the diameter of the tubular tissue; and an image generation section that generates a cylindrical projection image of the tubular tissue, based on information provided by projecting the virtual rays and the volume data.
According to one or more illustrative aspects of the present invention, there is provided a medical image processing method for visualizing a tubular tissue contained in volume data. The method comprises: (a) determining a central path of the tubular tissue based on volume data of voxel space containing the tubular tissue; (b) determining a diameter of the tubular tissue at a certain point on the central path; (c) determining at least two or more different projection spacings for projecting virtual rays along the central path, depending on the diameter of the tubular tissue; (d) projecting the virtual rays along the central path with the projection spacing that depends on the diameter of the tubular tissue; and (e) generating a cylindrical projection image of the tubular tissue, based on information provided by projecting the virtual rays and the volume data.
Other aspects of the invention will be apparent from the following description, the drawings and the claims.
In the accompanying drawings:
Exemplary embodiments of the present invention will be now described with reference to the drawings.
The X-ray source 401 radiates the X-ray beam bundle 402, which is shaped like a pyramid as indicated by the chain line in the figure. The X-ray detector 404 detects the X-ray beam bundle 402 passing through a patient 403 on the table 407. Further, the X-ray detector 404 outputs a signal of the detected X-ray beam bundle 402 to an medical image processing apparatus 100. The X-ray source 401 and the X-ray detector 404 are provided to face each other on the ring-like gantry 405.
The X-ray source 401 and the X-ray detector 404 are configured to rotate around a system axis 406 and move along the system axis 406 (i.e., movable relative to the patient 403). Thus, the X-ray beam bundle 402 is projected onto the patient 403 at various projection angles and various positions with respect to the system axis 406.
The ring-like gantry 405 is supported by a retainer (not shown) and rotatable (see arrow “a”) relative to the system axis 406 passing through the center point of the gantry 405.
The volume data generation section 101 receives a large number of successive tomographic signals in the diagnosis range of the patient 403 from the CT apparatus 400. The volume data generation section 101 generates volume data of voxel space containing a tubular tissue, based on the received signals. The volume data storage section 103 stores the volume data generated by the volume data generation section 101.
The central path determination section 105 determines the area of the tubular tissue existing in the voxel space based on the volume data obtained from the volume data storage section 103 and then determines the central path of the tubular tissue. The central path is a straight line or a curve.
The cross section acquisition section 107 reads the volume data of the cross section of the tubular tissue at a certain point on the central path determined by the central path determination section 105 from the volume data storage section 103. Then, the cross section acquisition section 107 generates a function representing the cross-sectional area of the tubular tissue. The cross-sectional area of the tubular tissue means the area which is surrounded by the tubular tissue.
The diameter determination section 109 determines the diameter of the tubular tissue at a certain point on the central path, based on the function generated by the cross section acquisition section 107. When the diameter determination section 109 determines the diameter of the tubular tissue, the diameter determination section 109 calculates an area S of the cross-sectional area of the tubular tissue at the point on the central path of the tubular tissue from the function generated by the cross section acquisition section 107, and then determines the square root of the area S as a diameter r of the tubular tissue at the point (r=√S).
The diameter determination section 109 may determine the maximum length of the lengths until the virtual rays projected on the cross-sectional area from the points on the central path are attenuated by the tubular tissue as the diameter r of the tubular tissue at each point. The diameter determination section 109 may determine the diameter of a circle inscribing or a circle circumscribing the cross-sectional area of the tubular tissue at each point on the central path as the diameter r of the tubular tissue at the point. The diameter determination section 109 may determine the diameter r of the tubular tissue without calculating a function representing the cross-sectional area of the tubular tissue in the cross section acquisition section 107. In this case, the diameter determination section 109 determines the area of the tubular tissue existing in the voxel space from the volume data stored in the volume data storage section 103, and then determines the diameter of a sphere inscribing the cross section of the area of the tubular tissue at each point on the central path determined by the central path determination section 105 as the diameter r of the tubular tissue at the point. The diameter determination section 109 may make adjustment such as determining the diameter r of the tubular tissue throughout the tubular tissue using the above-described determination methods of the diameter r of the tubular tissue in combination.
The spacing determination section 117 determines spacing on the central path for projecting the virtual rays in accordance with the diameter r of the tubular tissue at each point determined by the diameter determination section 109. The spacing on the central path for projecting the virtual rays is α×r (where α is a constant) ideally, but it is possible to make a correction in view of the characteristic that the diameter of the tubular tissue cannot strictly be defined. The diameters at peripheral points (for example, weighted average) are used in addition to the diameter at each point on the central path, whereby the effect of noise mixed in calculating the diameter r of the tubular tissue can be decreased.
The cylindrical projection section 111 projects the virtual rays along the central path of the tubular tissue with the spacing α×r determined by the spacing determination section 117 according to the cylindrical projection method. The virtual ray projection method of the cylindrical projection section 111 may be “curved cylindrical projection method” or “correction cylindrical projection method” described in U.S. Application Pub. No. 2007/120845 or “umbrella-type projection method” described in U.S. Application Pub. No. 2006/221074. That is, the virtual ray projection method may be a projection method of projecting each virtual ray from the central path used as the reference.
The image generation section 113 generates a cylindrical projection image of the tubular tissue by performing rendering based on the information provided by projecting virtual rays by the cylindrical projection section 111 and the volume data read from the volume data storage section 103. The image generation section 113 displays the generated cylindrical projection image on a display 151.
The image generation section 113 may display the numerical value of the diameter r of the tubular tissue at the attention point on the central path together with the cylindrical projection image. The image generation section 113 may display a ruler for measuring the size of the cylindrical projection image together with the cylindrical projection image. The user of the medical image processing apparatus 100 can visually grasp the size of the lesion or the tube wall of the tubular tissue while seeing the numerical value of the diameter r and the ruler.
The operation section 115 accepts operation of the user of the medical image processing apparatus 100 to set or change the attention point on the central path of the tubular tissue. The operation section 115 may be a keyboard or a mouse, for example.
The operation of the medical image processing apparatus 100 according to the exemplary embodiment will be now described with reference to
First of all, the central path determination section 105 determines a central path C of the tubular tissue based on the volume data obtained from the volume data storage section 103 (step S201). Next, the central path determination section 105 initializes a given position t on the central path C to t=0 (step S203).
Next, the cross section acquisition section 107 sets the position of a point X of the position t on the central path C to C(t) (step S205). Next, the cross section acquisition section 107 reads the volume data of the cross section of the tubular tissue at the point X from the volume data storage section 103 and then generates a function f that represents a cross-sectional area R of the tubular tissue (point X, cross-sectional area R) (step S207).
Next, the diameter determination section 109 determines the diameter r of the tubular tissue at the point X, based on the function f generated at step 5207 (step S209). Next, the cylindrical projection section 111 projects each virtual ray with the point X as the center according to the cylindrical projection method (step S211). That is, the cylindrical projection section 111 projects the virtual ray radially in the circumferential direction perpendicular to the central path C from the point X.
Next, the cylindrical projection section 111 changes the value “t” of the position C(t) of the point X to the value “t+α×r (where a is a constant)” (step S213). According to the step, the position t of projecting the virtual ray by the cylindrical projection section 111 moves by “α×r” along the central path C. That is, the move spacing of the virtual ray projection position t changes in response to the value of the diameter r of the tubular tissue, as shown in
Next, the cylindrical projection section 111 makes a comparison between the value of the position t on the central path C and maximum value t_max (step S215). If the value of the position t on the central path C is less than the maximum value t_max (YES at step S215), the process returns to step S205 and steps 5205 to S215 are repeatedly performed. Therefore, while the point X exists in a given range along the central path of the tubular tissue, the cylindrical projection section 111 performs cylindrical projection of virtual ray at step 5211. On the other hand, if the value of the position t on the central path C is equal to or greater than the maximum value t_max (NO at step S215), the process goes to step 5217.
At step 5217, the image generation section 113 generates a cylindrical projection image of the tubular tissue, based on the volume data read from the volume data storage section 103 and information provided by projecting virtual rays by the cylindrical projection section 111.
Projection of virtual ray and the cylindrical projection image of a tubular tissue will be now described with reference to
As shown in
When the virtual ray 13 is projected with the spacing, the lesions 20A and 20B are displayed in the cylindrical projection image of the tubular tissue 10 such that they have the same aspect ratio independent of the diameters of the tubular tissue where the lesions exist, as shown in
To generate a cylindrical projection image of a tubular tissue, the medical image processing apparatus 100 of the exemplary embodiment unfolds a cylindrical projection plane onto two-dimensional coordinates. The distance between two points and the angle between two line segments on the cylindrical projection image do not change depending on the unfolded position of the tubular tissue (virtual cylinder cut area) as shown in
Thus, as shown in
Next, a display method of a cylindrical projection image of a tubular tissue generated by the medical image processing apparatus 100 will be now described with reference to
A first display example of a cylindrical projection image of a tube wall of a tubular tissue will be now described with reference to
As shown in
If the size of the cylindrical projection image is thus changed in response to the diameter of the tubular tissue, the user can visually grasp the diameter of the tubular tissue and the size of the lesion. For example,
If two or more attention points are set, the image generation section 113 may display a plurality of cylindrical projection images with the sizes which correspond to the respective diameters of the tubular tissue at the attention points.
A second display example of a cylindrical projection image of a tube wall of a tubular tissue will be now described with reference to
As shown in
When the image processing load on the image generation section 113 varies from one method to another, if the method is changed in response to the region of the cylindrical projection image, the load on the image generation section 113 can be reduced. Further, the screen area can be saved.
As described above, in the medical image processing apparatus 100 according to the exemplary embodiment, the spacing for projecting virtual rays along the central path of a tubular tissue varies depending on the diameter of the tubular tissue. Thus, if the tubular tissue has different diameters, a cylindrical projection image of the tubular tissue can be generated with no distortion in the aspect ratio. Consequently, the user of the medical image processing apparatus 100 can precisely grasp the shape of a lesion existing in the tubular tissue.
While the present invention has been shown and described with reference to certain exemplary embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. It is aimed, therefore, to cover in the appended claim all such changes and modifications as fall within the true spirit and scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
P2008-311192 | Dec 2008 | JP | national |