Program, image processing device and image processing method

Abstract

To enable improvement in diagnostic accuracy while reducing the burden on a subject.SOLUTION: Based on a reference pixel value which is a maximum pixel value of an area corresponding to the upper side of the center of a heart in a tomographic image data before smoothing processing is performed, a setting part 3 sets a region where the myocardial distribution of a radiopharmaceutical accumulated in the myocardium is copied to the tomographic image data as a region of interest. An extraction part 4 extracts the myocardial tomographic image data showing the myocardial distribution from the tomographic image data based on the region of interest. An image processing part 5 performs smoothing processing on the myocardial tomographic image data.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] The present disclosure relates to imaging techniques for nuclear medicine. [Background technology]

[0002] Nuclear cardiac medicine tests, which examine the heart using a technique called SPECT (Single Photon Emission Computed Tomography) or PET (Positron Emission Tomography), in which a radiopharmaceutical is administered to a subject and the radiation emitted from the radiopharmaceutical is detected to obtain a three-dimensional nuclear medicine image, are attracting attention.

[0003] In cardiac nuclear medicine examinations, nuclear medicine images corresponding to radiation from radiopharmaceuticals accumulated in the myocardium are acquired, and cardiac examinations are performed based on the nuclear medicine images. However, because radiopharmaceuticals physiologically accumulate in the liver and are excreted through the digestive tract, artifacts due to radiation from the radiopharmaceuticals accumulated in the subdiaphragmatic organs, such as the liver and digestive tract, appear in the nuclear medicine images, affecting images of the myocardium, particularly the inferior wall of the myocardium adjacent to the subdiaphragmatic organs.

[0004] For example, radiation from organs below the diaphragm and radiation from the inferior wall of the myocardium may be affected by the partial volume effect, resulting in an increase in the radiation count value of the inferior wall of the myocardium. In this case, defects in the inferior wall of the myocardium may be obscured, or radiation count values ​​in areas other than the inferior wall of the myocardium may appear to be reduced. Furthermore, when image reconstruction is performed using filtered back-projection (FBP), the radiation count value of the inferior wall of the myocardium may also be reduced (see Non-Patent Document 1).

[0005] Non-Patent Document 2 describes a method for eliminating the influence of radiopharmaceutical accumulation in subdiaphragmatic organs in myocardial nuclear medicine images. Specifically, the method describes repeating imaging until the accumulation of radiopharmaceuticals in the subdiaphragmatic organs decreases or moves, and a method of adding prone position imaging to regular supine position imaging. In prone position imaging, the distance between the heart and the subdiaphragmatic organs increases, making it possible to reduce false images caused by accumulation in the subdiaphragmatic organs. [Prior art documents] [Non-patent literature]

[0006] [Non-Patent Document 1] Steven Burrell, et al. “Artifacts and Pitfalls in Myocardial PerfusionImaging”, Journal of Nuclear Medicine Technology, 2006 Dec;34(4):193-211. [Non-patent document 2] Ryan A Dvorak, et al. “Interpretation of SPECT / CT Myocardial Perfusion Images: Common Artifacts and Quality Control Techniques”, Radiographics, Nov-Dec 2011;31(7):2041-57. Summary of the Invention [Problem to be solved by the invention]

[0007] However, the methods described in Non-Patent Documents 1 and 2 do not necessarily reduce the influence of radiopharmaceuticals on myocardial nuclear medicine images due to accumulation of radiopharmaceuticals in organs below the diaphragm, and may actually reduce diagnostic accuracy. In addition, multiple imaging sessions are required, which increases the examination time and places a heavy burden on the subject.

[0008] The present disclosure has been made in consideration of the above-mentioned problems, and aims to provide a program, an image processing device, and an image processing method that can improve diagnostic accuracy while reducing the burden on the subject. [Means for solving the problem]

[0009] A program according to one aspect of the present disclosure causes a computer to implement a setting unit that sets, as a region of interest in the tomographic image data, a region that depicts the myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium, based on a reference pixel value that indicates the distribution of the radiopharmaceutical that has accumulated in the myocardium and organs below the diaphragm and that is the maximum pixel value of a region corresponding to the upper side of the center of the heart in the tomographic image data before smoothing processing is performed; an extraction unit that extracts myocardial tomographic image data that indicates the myocardial distribution from the tomographic image data based on the region of interest; and an image processing unit that performs smoothing processing on the myocardial tomographic image data.

[0010] In a preferred embodiment, the setting unit measures pixel values ​​of the tomographic image data along each of a plurality of radial lines extending three-dimensionally radially from a predetermined position within the cardiac cavity, and for each radial line, identifies a myocardial outer wall position corresponding to the outer wall of the myocardium in the tomographic image data based on highly concentrated pixels, which are pixels whose pixel values ​​are maximum and equal to or greater than a certain percentage of the reference pixel value, and sets the region of interest based on the myocardial outer wall position.

[0011] In addition, in a preferred embodiment, when there is one highly-accumulated pixel, the setting unit identifies a position that is a certain distance away from the highly-accumulated pixel toward the outside as the position of the myocardial outer wall, and when there are two or more highly-accumulated pixels, the setting unit identifies the position of the pixel having the smallest pixel value between the highly-accumulated pixel closest to the specified position and the highly-accumulated pixel that is second closest to the specified position as the position of the myocardial outer wall.

[0012] In a preferred aspect, the setting unit sets, as the region of interest, an approximately cardiac-shaped region that is a region connecting the positions of the myocardial outer wall on each radial line, or an ellipsoid region that approximates the shape of the approximately cardiac-shaped region to an ellipsoid.

[0013] In a preferred aspect, for each peripheral pixel that is a pixel on the periphery of the region of interest, if there is an outer accumulation pixel that is a pixel whose pixel value is a maximum and is equal to or greater than a certain percentage of the reference pixel value along an outward line extending outward from the peripheral pixel and within a specific distance from the peripheral pixel, the extraction unit extracts the myocardial tomographic image data by subtracting the pixel value of the accumulation below the diaphragm, which is a pixel value corresponding to the radiopharmaceutical that has accumulated in the organ below the diaphragm, from the pixel value of the pixel on the outward line based on the pixel values ​​along the outward line.

[0014] In a preferred aspect, the extraction unit approximates the change in pixel value from the peripheral pixel to the outer accumulation pixel along the outward line by a Gaussian function, and subtracts the pixel value represented by the Gaussian function from the pixel value of the pixel on the outward line as the pixel value of the accumulation below the diaphragm.

[0015] In a preferred aspect, the extraction unit approximates the change in pixel values ​​along the outward line in the region of interest with a Gaussian function, and sets the pixel value represented by the Gaussian function to the pixel value of the pixel on the outward line minus the pixel value of the accumulation below the diaphragm.

[0016] In a preferred aspect, for each of the peripheral pixels, if an outer accumulation pixel is present, the extraction unit sets all pixel values ​​of an area along the outward line that is outside the outer accumulation pixel to zero, and if no outer accumulation pixel is present, sets all pixel values ​​of an area outside a position that is a specific distance outward from the peripheral pixel of the tomographic image data to zero.

[0017] In a preferred embodiment, the image processing unit inversely converts the myocardial tomographic image data into projection data, and performs smoothing processing on the myocardial tomographic image data based on the projection data.

[0018] An image processing device according to one aspect of the present disclosure includes a setting unit that sets, as a region of interest in the tomographic image data, a region that depicts the myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium, based on a reference pixel value that indicates the distribution of the radiopharmaceutical that has accumulated in the myocardium and organs below the diaphragm and that is the maximum pixel value of a region corresponding to the upper side of the center of the heart in the tomographic image data before smoothing processing is performed; an extraction unit that extracts myocardial tomographic image data that indicates the myocardial distribution from the tomographic image data based on the region of interest; and an image processing unit that performs smoothing processing on the myocardial tomographic image data.

[0019] An image processing method according to one aspect of the present disclosure sets, as a region of interest in the tomographic image data, a region that shows the myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium, based on a reference pixel value that indicates the distribution of the radiopharmaceutical that has accumulated in the myocardium and organs below the diaphragm and that is the maximum pixel value in a region corresponding to the upper side of the center of the heart in the tomographic image data before smoothing processing is performed, extracts myocardial tomographic image data that shows the myocardial distribution from the tomographic image data based on the region of interest, and performs smoothing processing on the myocardial tomographic image data. [Effects of the Invention]

[0020] According to the present disclosure, it is possible to improve diagnostic accuracy while reducing the burden on subjects. [Brief explanation of the drawings]

[0021] [Figure 1] 1 is a block diagram illustrating a configuration of an image processing device according to an embodiment of the present disclosure. [Figure 2] 1 is a flowchart illustrating an operation of an image processing device according to an embodiment of the present disclosure. [Figure 3] 10 is a flowchart illustrating an example of a region of interest extraction process. [Figure 4] FIG. 10 is a diagram showing an example of radial lines in a slice image. [Figure 5]10A and 10B are diagrams for explaining a process for identifying the position of the myocardial outer wall. [Figure 6] 10 is a flowchart illustrating an example of myocardial extraction processing. [Figure 7] FIG. 10 is a diagram for explaining a process for identifying pixel values ​​of accumulation in subdiaphragmatic organs. [Figure 8] 1A and 1B are diagrams illustrating an image processing method of the present disclosure and an image processing method of a reference example. [Figure 9] 1A and 1B are diagrams showing myocardial tomographic image data obtained by the image processing method of the present disclosure and myocardial tomographic image data obtained by the image processing method of a reference example. DETAILED DESCRIPTION OF THE INVENTION

[0022] Hereinafter, embodiments of the present disclosure will be described with reference to the drawings.

[0023] Fig. 1 is a block diagram showing the configuration of an image processing device according to an embodiment of the present disclosure. The image processing device 100 shown in Fig. 1 is configured, for example, by a computer system including a processor (computer) and a memory (neither of which is shown). In this case, the components and functions of the image processing device 100 described below are realized, for example, by the processor reading a program and executing the read program. The program can be recorded on a computer-readable recording medium such as a memory.

[0024] The image processing device 100 is also connected to an input / output device 101 and an auxiliary storage device 102. The input / output device 101 includes input devices such as a keyboard, a touch panel, and a pointing device that receive various pieces of information from a user who uses the image processing device 100, and output devices such as a display device and a printer that output various pieces of information to the user. The input / output device 101 may also include a network interface device that transmits and receives various pieces of information via a communication network such as the Internet. The auxiliary storage device 102 is a storage device that stores various pieces of information, such as a large-capacity storage device.

[0025] The image processing device 100 includes an acquisition unit 1, a reconstruction unit 2, a setting unit 3, an extraction unit 4, and an image processing unit 5.

[0026] The acquisition unit 1 acquires projection data based on radiation from a subject administered with a radiopharmaceutical as data to be processed. The projection data is generated by imaging the subject using a radiation detector a certain time after the subject is administered with the radiopharmaceutical. Specifically, the projection data is data including a plurality of pixels arranged two-dimensionally for each predetermined rotation angle of the detector that rotates around the subject, and the pixel value of each pixel. The pixel value is a value determined according to a count value obtained by counting radiation from a part of the subject corresponding to the pixel. The pixel value may be, for example, the count value itself.

[0027] In this embodiment, the projection data is projection data for a myocardial nuclear medicine examination. In this case, the count value usually reflects not only the influence of radiation from the myocardium but also the influence of radiation from radiopharmaceuticals accumulated in organs below the diaphragm, such as the liver. Furthermore, the radiopharmaceutical is not particularly limited as long as it is a radiopharmaceutical used in a myocardial SPECT examination or a myocardial PET examination.

[0028] The reconstruction unit 2 generates tomographic image data based on the projection data acquired by the acquisition unit 1. Specifically, the reconstruction unit 2 performs image reconstruction processing on the projection data to generate horizontal tomographic image data showing horizontal tomographic images (body axis cross-sectional images) of the myocardium, and performs cross-section transformation processing on the horizontal tomographic image data to generate multiple tomographic image data with different cross-sectional orientations. In this embodiment, the tomographic image data includes horizontal long-axis tomographic image data showing horizontal long-axis tomographic images of the myocardium, vertical long-axis tomographic image data showing vertical long-axis tomographic images of the myocardium, and short-axis tomographic image data showing short-axis images of the myocardium.

[0029] Each tomographic image data reflects not only the distribution of radiopharmaceuticals accumulated in the myocardium, but also the distribution of radiopharmaceuticals accumulated in organs adjacent to the heart, such as the liver and digestive tract, below the diaphragm. Each tomographic image data is generated so that the center of the heart (more specifically, the center of the left ventricle of the heart) is positioned at the center of the image. Each tomographic image data is a collection of multiple slice images.

[0030] In this embodiment, the image reconstruction processing is image reconstruction processing using a statistical image reconstruction method including resolution correction, such as the Maximum Likelihood Expectation Maximization (MLEM) method and the Ordered Subset Expectation Maximization (OSEM) method. In this case, it is possible to reduce the influence of radiopharmaceutical accumulation in the subdiaphragmatic organs on the inferior wall of the myocardium. It also facilitates extraction of myocardial tomographic image data by the extraction unit 4, which will be described later. However, the image reconstruction processing may be image processing using a statistical image reconstruction method that does not include resolution correction, or image processing using a conventional filtered back-projection (FBP) method.

[0031] In particular, tomographic image data for myocardial SPECT examinations has a relatively small amount of information and a relatively large amount of statistical noise. Therefore, to improve visibility, a smoothing process is usually performed before or after image reconstruction to smooth out changes in pixel values ​​on the image. However, when smoothing is performed, the effect of radiopharmaceutical accumulation in the subdiaphragmatic organs on the inferior myocardial wall is emphasized. For this reason, in this embodiment, the reconstruction unit 2 does not perform smoothing.

[0032] The acquisition unit 1 may also acquire tomographic image data instead of projection data as data to be processed. In this case, there is no need to perform image reconstruction processing in the reconstruction unit 2. It is also assumed that the tomographic image data acquired by the acquisition unit 1 has not been subjected to smoothing processing.

[0033] The setting unit 3 sets a three-dimensional region that depicts myocardial distribution, which is the distribution of radiopharmaceuticals accumulated in the myocardium, as a region of interest in the tomographic image data. Specifically, the setting unit 3 sets the region that depicts myocardial distribution as a region of interest in the tomographic image data based on a reference pixel value that is the maximum pixel value in an upper region corresponding to the side above the center of the heart in the tomographic image data.

[0034] The extraction unit 4 extracts myocardial tomographic image data indicating myocardial distribution from the tomographic image data based on the region of interest set by the setting unit 3.

[0035] The image processing unit 5 performs various image processing on the myocardial tomographic image data extracted by the extraction unit 4 and outputs the result. Examples of image processing include smoothing and resection transformation. Alternatively, the image processing may involve inversely converting the myocardial tomographic image data into projection data, and then performing image reconstruction, resection transformation, and smoothing on the projection data to generate desired tomographic image data. In this case, it becomes possible to generate tomographic image data with a desired cross-sectional orientation.

[0036] FIG. 2 is a flowchart for explaining the operation of the image processing device 100.

[0037] First, the acquisition unit 1 acquires data to be processed (step S101). For example, the acquisition unit 1 may acquire data from an arbitrary external device, such as an imaging device that images a subject using a radiation detector, via the input / output device 101, or may acquire data stored in the auxiliary storage device 102.

[0038] The acquisition unit 1 determines whether the acquired data is projection data (step S102).

[0039] If the acquired data is projection data, the reconstruction unit 2 generates tomographic image data based on the projection data (step S103). On the other hand, if the acquired data is not projection data, that is, if the acquired data is tomographic image data, the processing of step S103 by the reconstruction unit 2 is skipped. Note that it is assumed that no smoothing processing has been performed on the tomographic image data.

[0040] The setting unit 3 executes a region of interest extraction process (see FIGS. 3 to 5) for setting a region of interest in the tomographic image data acquired by the acquisition unit 1 or generated by the reconstruction unit 2 (step S104).

[0041] The extraction unit 4 executes myocardial extraction processing (see FIGS. 6 and 7) to extract myocardial tomographic image data from the tomographic image data based on the region of interest set by the setting unit 3 (step S105).

[0042] The image processing unit 5 determines whether or not to generate projection data from the myocardial tomographic image data (step S106). For example, whether or not to generate projection data may be specified in advance, or the image processing unit 5 may inquire of the user at this stage whether or not to generate projection data.

[0043] If projection data is not to be generated, the image processing unit 5 performs predetermined image processing including smoothing on the myocardial tomographic image data (step S107). Then, the image processing unit 5 displays the processed myocardial tomographic image data on the input / output device 101 (step S108), and ends the processing. Note that instead of displaying the myocardial accumulation data on the input / output device 101, the image processing unit 5 may transmit the myocardial accumulation data to an external device (not shown) via the input / output device 101, or may store the data in the auxiliary storage device 102.

[0044] On the other hand, when generating projection data, the image processing unit 5 generates the projection data from the myocardial tomographic image data (step S109). Then, the image processing unit 5 generates desired tomographic image data as myocardial tomographic image data from the generated projection data, and performs predetermined image processing including smoothing processing on the myocardial tomographic image data (step S110). Then, the image processing unit 5 proceeds to the processing of step S108.

[0045] FIG. 3 is a flowchart illustrating an example of the region of interest extraction process in step S104.

[0046] In the region of interest extraction process, the setting unit 3 first determines whether or not to perform precision setting for precisely setting the region of interest (step S201). For example, whether or not to perform precision setting may be specified in advance, or the setting unit 3 may inquire of the user at this stage whether or not to perform precision setting.

[0047] If precise setting is not required, the setting unit 3 performs simple setting to simply set the region of interest (step S202), and then ends the process. Methods for simply setting the region of interest include, for example, a method of setting an ellipsoidal region of a predetermined size centered on the center of each tomographic image data as the region of interest, and a method of having the user manually set the region of interest.

[0048] When performing precise setting, the setting unit 3 first obtains a reference pixel value, which is the maximum pixel value in the upper region corresponding to the area above the center of the heart in the tomographic image data (step S203). Since the upper region is an area that is less susceptible to the influence of radiopharmaceutical accumulation in organs below the diaphragm, the reference pixel value can be considered to be the maximum pixel value that is not influenced by the accumulation of radiopharmaceuticals in organs below the diaphragm. In this embodiment, the upper region is assumed to be the upper half of the short-axis tomographic image data. In this case, the reference pixel value can be set more appropriately.

[0049] The upper region may be any region above the center of the heart, and may be limited to, for example, a portion of the upper half of the short-axis tomographic image data. Furthermore, the range of the upper region in the long-axis direction is not limited to the range from the apex to the base of the heart (0% to 100%), assuming that the entire length along the long axis of the heart is 100% and the apex is the starting point, but may more preferably be limited to a range of 30% to 90%.

[0050] Next, the setting unit 3 measures pixel values ​​of the tomographic image data along each of a plurality of radial lines extending three-dimensionally radially from the center of the heart (step S204). In this embodiment, the setting unit 3 measures pixel values ​​from the center of the image outward along each of a plurality of radial lines extending radially from the image center at a fixed angle in the central slice image of each of the tomographic image data, which are the horizontal long-axis tomographic image data, the vertical long-axis tomographic image data, and the short-axis tomographic image data. The fixed angle is, for example, 3 degrees. Note that the starting point of the radial lines is not limited to the center of the heart, and may be any position within the cardiac cavity. Furthermore, the starting point of the three-dimensional radial lines is not limited to the center of the heart, and may be any predetermined position within the cardiac cavity.

[0051] 4 is a diagram showing an example of radial lines in one slice image, in which 40 radial lines X extending from the center O of the heart are shown at a fixed angle of 9 degrees.

[0052] Returning to the explanation of Fig. 3, after measuring the pixel values ​​of each radial line, the setting unit 3 identifies, for each radial line, the pixel position where the pixel value is maximum and equal to or greater than a threshold value that is a certain percentage (for example, 20%) of the reference pixel value as a high accumulation position where the radiopharmaceutical has accumulated (step S205). Note that the high accumulation position is considered to be a position corresponding to the myocardium.

[0053] The setting unit 3 identifies, for each radial line, a myocardial outer wall position, which is a position corresponding to the outer wall of the myocardium, based on the high accumulation positions (step S206). Specifically, the setting unit 3 checks the number of high accumulation positions for each radial line, and identifies the myocardial outer wall position based on the high accumulation positions and the number of high accumulation positions.

[0054] Fig. 5 is a diagram for explaining the process of identifying the position of the myocardial outer wall, showing the change in pixel value along radial lines. In Fig. 5, the horizontal axis represents the position along the radial lines with the center O of the heart set as 0, and the vertical axis represents the pixel value. The threshold value is T.

[0055] 5(a) shows an example in which there is one high accumulation position. In this case, the setting unit 3 determines that there is no accumulation in the subdiaphragmatic organs in the direction along the radial lines, and specifies a position a certain distance L outward from the high accumulation position M as the myocardial outer wall position W. The certain distance L is, for example, 2 cm.

[0056] FIG. 5(b) shows an example where there are two or more high accumulation positions. In this case, the setting unit 3 determines that the high accumulation position M1 closest to the center O of the heart is the position of the myocardium, and the high accumulation position M2 second closest to the center O of the heart is the position of accumulation of the radiopharmaceutical in the organs below the diaphragm, and specifies the position with the smallest pixel value between the high accumulation positions M1 and M2 as the myocardial outer wall position W. In this case, if the distance between the high accumulation positions M1 and M2 is equal to or greater than a predetermined distance (e.g., 2 cm), the setting unit 3 may specify a position that is a predetermined distance outward from the high accumulation position M1 as the myocardial outer wall position W. Note that the predetermined distance may be the same as the constant distance L shown in FIG. 4(a).

[0057] 5(c) shows an example in which the high accumulation position is 0. In this case, the setting unit 3 determines that there is no myocardium in the direction along the radial lines, and does not identify the myocardial outer wall position W.

[0058] Returning to the description of Fig. 3, when the myocardial outer wall position is identified for each radial line, the setting unit 3 sets the region of interest based on the myocardial outer wall position of each radial line (step S207), and ends the process.

[0059] The region of interest is an approximately heart-shaped region formed by connecting the positions of the myocardial outer wall of each radial line of each tomographic image data, or an ellipsoidal region formed by approximating the shape of the approximately heart-shaped region to an ellipsoid. The method for approximating the approximately heart-shaped region to an ellipsoidal region is not particularly limited, and may be, for example, the least squares method. Whether the region of interest is an approximately heart-shaped region or an ellipsoidal region may be set in advance or may be set by the user.

[0060] The region of interest may be adjusted manually. For example, if the region of interest is an ellipsoidal region, adjustment of the region of interest may involve moving the entire region of interest, moving the three axes of the ellipsoid, deforming and rotating the region, etc. Also, if the region of interest is an approximately cardiac region, adjustment of the region of interest may involve mesh deformation, moving the entire region of interest, correcting each pixel, etc.

[0061] FIG. 6 is a flowchart illustrating an example of the myocardium extraction process in step S105.

[0062] In the myocardial extraction process, the extraction unit 4 determines whether to perform estimated removal to remove the influence of the accumulation of the radiopharmaceutical in the subdiaphragmatic organs on the region of interest (step S301). For example, whether to perform estimated removal may be specified in advance, or the extraction unit 4 may inquire of the user at this stage whether to perform estimated removal.

[0063] If estimation and removal is not performed, the extracting unit 4 extracts myocardial tomographic image data by setting all pixel values ​​in regions other than the region of interest of each piece of tomographic image data to zero (step S302), and ends the process.

[0064] When performing estimation and removal, the extraction unit 4 analyzes the tomographic image data and estimates pixel values ​​due to accumulation of the radiopharmaceutical in the subdiaphragmatic organs (step S303).

[0065] Specifically, the extraction unit 4 first measures the pixel value along an outward line extending outward from each periphery pixel, which is a pixel on the periphery of the region of interest in the tomographic image data. Specifically, the outward line is a line passing through the periphery pixel and the center of the region of interest, or a line passing through the periphery pixel and perpendicular to the periphery of the region of interest.

[0066] Next, the extraction unit 4 determines for each peripheral pixel whether or not there is an outer accumulation pixel, which is a pixel whose pixel value is a maximum and is equal to or greater than a threshold value that is a certain percentage (e.g., 20%) of the reference pixel value, within a position a specific distance (e.g., 1 cm) along the outward line.

[0067] If there are no outer accumulation pixels, the extraction unit 4 determines that there are no pixel values ​​due to accumulation of the subdiaphragmatic organs on the outward line, and does not estimate pixel values ​​of accumulation of the subdiaphragmatic organs.

[0068] On the other hand, if outer accumulation pixels are present, the extraction unit 4 estimates pixel values ​​of accumulation in the subdiaphragmatic organs based on pixel values ​​from the peripheral pixels on the outward line to the outer accumulation pixels.

[0069] Figure 7 is a diagram illustrating the process of estimating pixel values ​​of accumulation in subdiaphragmatic organs based on pixel values ​​from the peripheral pixels to the outer accumulation pixels. In Figure 7, the horizontal axis represents the position along the outward line, and the vertical axis represents pixel values. Also, T represents the threshold, R represents the specific distance, C represents the peripheral pixels, and D represents the outer accumulation pixels.

[0070] The extraction unit 4 performs fitting to approximate the change F in pixel values ​​on the outward line from the peripheral pixel C to the outer accumulation pixel D with a Gaussian function, and obtains a Gaussian function G that approximates the change F. Then, the extraction unit 4 identifies the pixel values ​​represented by the Gaussian function G as pixel values ​​of accumulation in the subdiaphragmatic organs. Note that a function other than a Gaussian function may be used as the function that approximates the change F.

[0071] Returning to the explanation of Fig. 6, after estimating the pixel values ​​of the accumulation in the subdiaphragmatic organs, the extraction unit 4 subtracts the pixel values ​​of the accumulation in the subdiaphragmatic organs corresponding to each peripheral pixel from the pixel values ​​of the tomographic image data to extract myocardial image data (step S304), and ends the process. In the example of Fig. 7, the pixel values ​​represented by the curve H obtained by subtracting the pixel values ​​of the accumulation in the subdiaphragmatic organs represented by the Gaussian function G from the pixel values ​​represented by the change in pixel value F become the pixel values ​​of the myocardial image data.

[0072] In step S304, if there is a pixel whose pixel value of the tomographic image data minus the pixel value of the accumulation in the subdiaphragmatic organs is negative, the extraction unit 4 sets the pixel value of the pixel to 0. Furthermore, the extraction unit 4 sets the pixel value of a pixel outside the outer accumulation pixel or outside a position a specific distance away from the outer peripheral pixel along the outward line to 0.

[0073] In step S303, the change F in pixel value from the peripheral pixel C to the outer accumulation pixel D is approximated by the Gaussian function G. However, for example, the change in pixel value from the high accumulation position M in the region of interest to the peripheral pixel C may also be approximated by the Gaussian function. In this case, the myocardial image data is extracted by replacing the pixels outside the high accumulation position M with pixels represented by the Gaussian function.

[0074] As described above, according to this embodiment, the setting unit 3 sets a region of interest in the tomographic image data that depicts the myocardial distribution, which is the distribution of the radiopharmaceutical accumulated in the myocardium, based on a reference pixel value, which is the maximum pixel value in a region above the center of the heart in the tomographic image data before smoothing processing. The extraction unit 4 extracts myocardial tomographic image data that shows the myocardial distribution from the tomographic image data based on the region of interest. The image processing unit 5 performs smoothing processing on the myocardial tomographic image data. Therefore, the smoothing processing is performed after the myocardial tomographic image data that shows the myocardial distribution is extracted from the tomographic image data. Therefore, it is possible to generate image data with high visibility while removing the influence of radiopharmaceutical accumulation in organs below the diaphragm in the tomographic image data, without having to perform multiple imaging of the subject, thereby reducing the burden on the subject and improving diagnostic accuracy.

[0075] In this embodiment, the setting unit 3 measures pixel values ​​of the tomographic image data along each of a plurality of radial lines extending three-dimensionally radially from the center of the heart, and for each radial line, identifies the position of the myocardial outer wall corresponding to the outer wall of the myocardium in the tomographic image data based on the highly-dense pixels, which are pixels whose pixel values ​​are maximum and equal to or greater than a certain percentage of the reference pixel value, and sets a region of interest based on the myocardial outer wall position. This makes it possible to appropriately set the region of interest.

[0076] Furthermore, in this embodiment, when there is one highly accumulating pixel, the setting unit 3 specifies a position a certain distance outward from the highly accumulating pixel as the myocardial outer wall position, and when there are two or more highly accumulating pixels, the setting unit 3 specifies the position of the pixel having the smallest pixel value between the highly accumulating pixel closest to the center and the second-closest to the center as the myocardial outer wall position. This makes it possible to appropriately specify the myocardial outer wall position, thereby enabling more appropriate setting of the region of interest.

[0077] In this embodiment, the region of interest is an approximately heart-shaped region that connects the positions of the myocardial outer wall on the radial lines, or an ellipsoid region that approximates the shape of the approximately heart-shaped region to an ellipsoid, making it possible to set the region of interest to an appropriate region that captures the heart.

[0078] Furthermore, in this embodiment, for each peripheral pixel on the periphery of the region of interest, if there is an outer accumulation pixel, which is a pixel whose pixel value is a maximum and is equal to or greater than a certain percentage of the reference pixel value, along an outward line extending outward from the peripheral pixel and within a specific distance from the peripheral pixel, the extraction unit 4 extracts myocardial tomographic image data by subtracting the pixel value of the accumulation in the subdiaphragmatic organs, which is a pixel value corresponding to the radiopharmaceutical accumulated in the subdiaphragmatic organs, from the pixel values ​​of the pixels on the outward line, based on the pixel values ​​along the outward line. This allows myocardial tomographic image data to be extracted appropriately.

[0079] In this embodiment, the extractor 4 approximates the change in pixel values ​​from the peripheral pixels along the outward line to the outer accumulation pixels using a Gaussian function, and subtracts the pixel values ​​expressed by the Gaussian function from the pixel values ​​of the pixels on the outward line as pixel values ​​of accumulation in the subdiaphragmatic organs. This makes it possible to appropriately determine the pixel values ​​of accumulation in the subdiaphragmatic organs, thereby enabling more appropriate extraction of myocardial tomographic image data.

[0080] In this embodiment, the extraction unit 4 sets all pixel values ​​of the region outside the outer accumulation pixel along the outward line to zero if an outer accumulation pixel is present for each outer peripheral pixel, and sets all pixel values ​​of the region outside a specific distance outward from the outer peripheral pixel of the tomographic image data to zero if no outer accumulation pixel is present. This makes it possible to appropriately extract the myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium.

[0081] In this embodiment, the image processing unit 5 converts the myocardial tomographic image data inversely into projection data and performs smoothing processing on the myocardial tomographic image data based on the projection data, thereby making it possible to generate a desired myocardial tomographic image. [Example]

[0082] Here, in order to evaluate the effect of radiopharmaceutical accumulation in subdiaphragmatic organs on myocardial tomographic image data, myocardial tomographic image data obtained by the image processing method of the present disclosure was compared with myocardial tomographic image data obtained by the image processing method of the reference example.

[0083] 8A and 8B are diagrams showing an image processing method according to the present disclosure and an image processing method according to a reference example. Specifically, Fig. 8A shows the image processing method according to the reference example, and Fig. 8B shows the image processing method according to the present disclosure.

[0084] In the image processing method of the reference example shown in Fig. 8(a), a smoothing process is performed on the tomographic image data generated by the reconstruction process, followed by a resection transformation process and a myocardial extraction process. On the other hand, in the image processing method of the present disclosure shown in Fig. 8(b), as described in the embodiment, a smoothing process is performed on the tomographic image data generated by the reconstruction process, followed by a resection transformation process and a myocardial extraction process. Note that in the myocardial extraction process, all pixel values ​​in regions other than the region of interest in each tomographic image data are set to zero without performing estimation removal to remove the influence of radiopharmaceutical accumulation in organs below the diaphragm on the region of interest.

[0085] FIG. 9 is a diagram showing myocardial tomographic image data obtained by the image processing method of the reference example and myocardial tomographic image data obtained by the image processing method of the present disclosure.

[0086] The projection data was acquired using a Siemens Symbia T6 equipped with a low-energy, high-resolution collimator, and data was processed using a Siemens nuclear medicine image processing workstation, SYNGO Mi Apps (version VA60C), and DRIP and cardioBULL, manufactured by Fujifilm Toyama Chemical Co., Ltd. The dose calibrator used was an Aloka Curiemeter IGC-7.

[0087] The subject used was a Kyoto Chemical Industry Co., Ltd. heart-liver HL-type phantom. A polypropylene cylindrical container (3.5 cm lumen diameter, 10.0 cm length, 100.0 mL volume) simulating a subdiaphragmatic organ collection was attached to the outside of the inferior wall of the left ventricular myocardium. Two types of subjects were prepared: one with no left ventricular myocardial defect and one with a 2.0 cm diameter circular defect chip attached to the center of the inferior wall of the left ventricular myocardium. In the following, tomographic image data acquired when there was no left ventricular myocardial defect will be referred to as defect-free data, and tomographic image data acquired when a circular defect chip was attached will be referred to as defect-containing data.

[0088] Radiopharmaceutical nuclides 99m Tc into the left ventricular myocardium (volume 120 mL) 99mThe amount of Tc enclosed was set to 14.8 MBq (123 kBq / mL), assuming that approximately 2% of the administered dose will accumulate in the myocardium in clinical practice. The cylindrical container was filled with water so that the radioconcentration ratio of the cylindrical container to the left ventricular myocardium was 0 (water only), 0.5, 0.75, 1.0, 1.25, 1.5, 2.0, and 3.0. 99m Tc was sealed in, and Figure 9 shows the myocardial tomographic image data (specifically, the central slice of the short-axis tomographic image data). Only water was sealed in the left and right ventricular cavities, mediastinum, liver, and stomach. The image reconstruction methods used for image reconstruction processing were FBP, OSEM, OSEM-PR, which is OSEM with resolution recovery (RR), and OSEM-ACSCRR, which is OSEM with attenuation correction (AC), scatter correction (SC), and resolution recovery (PR). Figure 9 shows the myocardial tomographic image data obtained by each image reconstruction method.

[0089] Each myocardial tomographic image data was evaluated as follows.

[0090] Specifically, circumferential profile curves (CPCs) were created at 1-degree intervals within the region of interest (ROI) in the central slice of the short-axis tomographic image data, and the maximum pixel value in the upper half of the region was normalized to 100%. Clinically, a %Uptake (the ratio of pixel values ​​to the maximum pixel value) of 70% or less is often considered to indicate significant blood flow abnormalities. Therefore, for data with missing data, the influence of subdiaphragmatic organ accumulation was deemed significant when the minimum value in the 90-degree region on the inferior wall was 70% or greater. For data without missing data, the influence of subdiaphragmatic organ accumulation was deemed significant when the maximum value in the 90-degree region on the inferior wall was 143 (100 / 0.7) or greater, because blood flow abnormalities in walls other than the inferior wall may be apparent.

[0091] For missing data, in the image processing method of the reference example, when the image reconstruction method was FBP or OSEM, the radioactivity concentration ratio of the cylindrical container was 0.75 or more and the %Uptake was 70% or more, and when the image reconstruction method was OSEM-PR or OSEM-ACSCPR using PR, the radioactivity concentration ratio of the cylindrical container was 1.0 or more and the %Uptake was 70% or more.

[0092] On the other hand, in the image processing method of the present disclosure, when the image reconstruction method was FBP, the %Uptake was 70% or more when the radioactivity concentration ratio of the cylindrical container was 1.75 or more, when the image reconstruction method was OSEM, the %Uptake was 70% or more when the radioactivity concentration ratio of the cylindrical container was 1.25 or more, and when the image reconstruction method was OSEM-PR, the %Uptake was 70% or more when the radioactivity concentration ratio of the cylindrical container was 3.0 or more. Furthermore, when the image reconstruction method was OSEM-ACSCPR, the %Uptake was less than 70% at all radioactivity concentration ratios.

[0093] Furthermore, for data without defects, in the image processing method of the reference example, when the image reconstruction method was FBP, the %Uptake was 143% or more when the radioactivity concentration ratio of the cylindrical container was 2.0 or more, and in the other methods, the %Uptake was 143% or more when the radioactivity concentration ratio of the cylindrical container was 1.5 or more. On the other hand, in the image processing method of the present disclosure, the Uptake was less than 143% for all image reconstruction methods.

[0094] Therefore, it was demonstrated that the image processing method of the present disclosure can reduce the influence of accumulation of subdiaphragmatic organs, thereby improving diagnostic accuracy.

[0095] The above-described embodiments and examples of the present disclosure are merely illustrative examples of the present disclosure, and are not intended to limit the scope of the present disclosure to these embodiments alone. Those skilled in the art may implement the present disclosure in various other forms without departing from the scope of the present disclosure. [Explanation of symbols]

[0096] 1: Acquisition unit, 2: Reconstruction unit, 3: Setting unit, 4: Extraction unit, 5: Image processing unit, 100: Image processing device, 101: Input / output device, 102: Auxiliary storage device

Claims

1. An image processing method performed by an image processing device, comprising: the image processing device sets a region that shows a myocardial distribution, which is the distribution of the radiopharmaceutical accumulated in the myocardium, on a tomographic image that shows the distribution of the radiopharmaceutical used in the myocardial SPECT examination or the myocardial PET examination; the image processing device removes the influence of accumulation of organs under the diaphragm in a region other than the region imaging the myocardial distribution from the tomographic image, thereby extracting a myocardial tomographic image consisting of the region imaging the myocardial distribution and the region other than the region imaging the myocardial distribution from which the influence has been removed; The image processing method includes performing a smoothing process on the myocardial tomographic image by the image processing device.

2. The image processing method according to claim 1 , wherein in extracting the myocardial tomographic image, the image processing device removes the influence by subtracting pixel values ​​indicating accumulation of the organs below the diaphragm from the tomographic image.

3. The image processing method according to claim 1 , wherein, in extracting the myocardial tomographic image, the image processing device removes the influence by setting pixel values ​​of regions other than the region showing the myocardial distribution to zero.

4. The image processing method according to claim 1 , wherein in setting the region showing the myocardial distribution, the image processing device sets a region specified by a user as the region showing the myocardial distribution.

5. The image processing method according to claim 1 , wherein the tomographic image is a tomographic image before smoothing processing is performed.

6. The image processing method according to claim 1 , wherein the smoothing process is a process for smoothing changes in pixel values ​​on the image.

7. a setting unit that sets a region that shows a myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium, in a tomographic image that shows the distribution of the radiopharmaceutical used in the myocardial SPECT examination or the myocardial PET examination; an extracting unit that removes the influence of accumulation of organs under the diaphragm in a region other than the region where the myocardial distribution is imaged from the tomographic image, thereby extracting a myocardial tomographic image consisting of the region where the myocardial distribution is imaged and the region other than the region where the myocardial distribution is imaged from which the influence has been removed; a smoothing unit that performs smoothing processing on the myocardial tomographic image; and a program for causing a computer to realize the above-mentioned smoothing unit.

8. a setting unit that sets a region that shows a myocardial distribution, which is the distribution of the radiopharmaceutical that has accumulated in the myocardium, in a tomographic image that shows the distribution of the radiopharmaceutical used in the myocardial SPECT examination or the myocardial PET examination; an extracting unit that removes the influence of accumulation of organs under the diaphragm in a region other than the region where the myocardial distribution is imaged from the tomographic image, thereby extracting a myocardial tomographic image consisting of the region where the myocardial distribution is imaged and the region other than the region where the myocardial distribution is imaged from which the influence has been removed; a smoothing unit that performs a smoothing process on the myocardial tomographic image.

Images