CN104586416A - Method and device for generating attenuation diagrams and positron emission cross-sectional images - Google Patents

Method and device for generating attenuation diagrams and positron emission cross-sectional images Download PDF

Info

Publication number
CN104586416A
CN104586416A CN201310533167.XA CN201310533167A CN104586416A CN 104586416 A CN104586416 A CN 104586416A CN 201310533167 A CN201310533167 A CN 201310533167A CN 104586416 A CN104586416 A CN 104586416A
Authority
CN
China
Prior art keywords
pixel
image
die body
pet
profile 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.)
Granted
Application number
CN201310533167.XA
Other languages
Chinese (zh)
Other versions
CN104586416B (en
Inventor
徐天艺
董晓霞
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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201310533167.XA priority Critical patent/CN104586416B/en
Priority to CN201710224826.XA priority patent/CN106963406B/en
Publication of CN104586416A publication Critical patent/CN104586416A/en
Application granted granted Critical
Publication of CN104586416B publication Critical patent/CN104586416B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58Testing, adjusting or calibrating thereof
    • A61B6/582Calibration

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine (AREA)
  • Image Processing (AREA)

Abstract

The invention provides a method and device for generating homogeneous phantom attenuation diagrams. The method comprises the steps of determining pixel points, forming a homogeneous phantom, on PET images of the generated homogeneous phantom, obtaining position information of the pixel points, and generating the attenuation diagrams of the PET images through the position information of the pixel points and attenuation values of the obtained homogeneous phantom. According to the method and device for generating the attenuation diagrams, CT equipment is not needed, and needed attenuation diagrams are directly generated from data in the PET images.

Description

Decay pattern and positron emission fault image generating method and device
Technical field
The present invention relates to field of medical technology, be specifically related to the generation method and apparatus of a kind of decay pattern and positron emission fault image.
Background technology
Positron emission tomography art (Positron Emission Tomography, PET) is the technology that current only employing anatomic form mode carries out function, metabolism and rii receptor, has now been widely used in the diagnoses and treatment of various diseases.
PET make use of the feature of positron emitter and positron annihilation, to determine the function of organ and cell compartment quantitatively.Utilize PET technology, before an examination, to patient, there is radioisotopic suitable radiopharmaceutical label.Along with they decay, radiosiotope radiates positron, and described positron is interacting with the negatron in tissue after a short distance, generation annihilation radiation.In the process of annihilation radiation, produce two γ photons that energy is equal, direction is contrary, described γ photon produces gamma-rays in patient body.
With the gamma-rays distribution situation in different angles detection patient body, the data collected thus are for the formation of the image of radiosiotope distribution in patient body, and be called that PET launches image, the forming process of this image is called image reconstruction.
In PET imaging process, due to the impact by the attenuation factor such as Compton effect, scattering, random signals event, dead time, described in the data that collect inconsistent with practical situation, calculate the PET transmitting image fault of formation thus.Therefore, must take effective measures and correct, just can obtain more real medical image.
At present, usually adopt computerized tomography (Computed Tomography, CT) to launch image to PET and carry out correction for attenuation.Utilize the dampening information of CT technical limit spacing gamma-rays in tissue, according to the dampening information generation of attenuation map got, thus utilize described decay pattern to launch image to PET carry out correction for attenuation, to obtain evenly, higher, the better PET image of picture quality of definition.
But, utilize CT technology to launch image to PET and carry out correction for attenuation, need the equipment of CT technology described in additional application, add corresponding cost and operational complexity.Further, the use of CT technology, causes patient to receive higher radioactive exposure, and the safety of patient is by impact to a certain extent.
In prior art, complete the correction for attenuation that PET launches image, need to utilize CT technology, this just needs on the basis of PET equipment, the equipment of CT technology described in additional application, thus the cost of equipment is increased.Further, operate CT equipment and PET equipment, not only increase operational complexity, the safety of patient also receives impact to a certain degree simultaneously.
For solving the problem, in embodiments of the invention, adopt even die body simulated human tissue, the PET for even die body launches image and provides a kind of decay map generalization method.In described decay map generalization method, by determining the pixel forming described even die body in the PET image of even die body, thus obtain the positional information of the pixel of described even die body, again according to the pad value of described even die body, the decay pattern that even die body PET launches image can be obtained.
Described even die body, is commonly referred to as the object that physical characteristic is relatively uniform, such as: the water of certain volume.
For enabling above-mentioned purpose of the present invention, feature and advantage more become apparent, and are described in detail specific embodiments of the invention below in conjunction with accompanying drawing.
Summary of the invention
The problem that the embodiment of the present invention solves is how when not by being even die body generation of attenuation map when CT technology.
For solving the problem, the embodiment of the present invention provides a kind of evenly die body decay map generalization method, and described method comprises:
Determine that the positron emission fault PET of generated even die body launches pixel image being formed described even die body, obtain the positional information of described pixel;
The pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
Alternatively, the described PET transmitting image determining generated even die body is formed the pixel of described even die body, comprising:
Determine that described PET launches profile region image being formed described even die body, the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
Alternatively, describedly determine that described PET launches profile region image being formed described even die body, comprising:
Obtain the pixel value that described PET launches pixel on image, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
Adjust described first threshold, to obtain on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference.
Alternatively, the described first threshold of described adjustment, to obtain on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image, comprising:
With the initial value that described max pixel value is described first threshold, reduce described first threshold, obtain the first difference and the second difference, until described second difference is greater than described first difference, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image, described second difference for current once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image,
Profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference, comprising: profile region image being formed described even die body is launched as described PET in the first profile region corresponding to described first difference.
Alternatively, the described first threshold of described reduction comprises: each ratio reducing described first threshold is 1% to 5%.
Alternatively, describedly determine that described PET launches pixel image being formed described even die body, comprising:
Obtain the pixel value that described PET launches pixel on image;
Described pixel value is made from big to small to the accumulation calculating of pixel number, obtain added up a pixel numerical value;
When described added up a pixel numerical value equals the pixel number of the position of corresponding described PET transmitting image on described even die body, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
Alternatively, the described accumulation calculating described pixel value being made from big to small to pixel number, comprising: make statistics with histogram to described pixel value, according to described histogrammic statistical result, described pixel value is made from big to small to the accumulation calculating of pixel number.
Alternatively, embodiments of the invention additionally provide a kind of generation method of even die body PET image, and described method comprises: the PET generating even die body launches image; Adopt the decay pattern that above-mentioned decay drawing generating method generates, image is launched to described PET and carries out correction for attenuation, generate the PET image through correction for attenuation.
Embodiments of the invention additionally provide a kind of generating apparatus of even die body decay pattern, and described device comprises:
Acquiring unit, for determining that the PET of generated even die body launches pixel image being formed described even die body, obtains the positional information of described pixel;
Generation unit, for the pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
Alternatively, described acquiring unit comprises: profile region obtains subelement, for determining that described PET launches profile region image being formed described even die body, the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
Alternatively, described profile region acquisition subelement comprises:
First profile region acquisition module, the pixel value of pixel on image is launched for obtaining described PET, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
First threshold adjusting module, for adjusting described first threshold, to obtain on described first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Second profile region acquisition module, for launching as described PET profile region image being formed described even die body using the first profile region corresponding time minimum for described difference.
Alternatively, described first threshold adjusting module is used for the initial value of described max pixel value as described first threshold, reduce described first threshold, obtain the first difference and the second difference, until described first difference is greater than described second difference, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image, described second difference for current once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image,
Described profile region acquisition module is for launching as described PET profile region image being formed described even die body using the first profile region corresponding to described first difference.
Alternatively, described acquiring unit comprises:
First obtains subelement, launches the pixel value of pixel on image for obtaining described PET;
Second obtains subelement, for making the accumulation calculating of pixel number from big to small to described pixel value, obtains added up a pixel numerical value;
3rd obtains subelement, during for equaling the pixel number of the position of corresponding described PET transmitting image on described even die body when described added up a pixel numerical value, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
Alternatively, described second acquisition subelement comprises:
Statistics with histogram module, for making statistics with histogram to described pixel value, obtains statistics with histogram result;
Computing module, for the histogrammic statistical result obtained according to described statistics with histogram module, makes the accumulation calculating of pixel number from big to small to described pixel value.
Embodiments of the invention additionally provide a kind of generating apparatus of even die body PET image, and described device comprises:
First generation unit, described first generation unit launches image for generating PET;
Second generation unit, for the decay pattern adopting above-mentioned decay pattern generating apparatus to generate, launches image to described PET and carries out correction for attenuation, generate the PET image through correction for attenuation.
Compared with prior art, the technical scheme of the embodiment of the present invention has the following advantages:
In embodiments of the invention, by determining the pixel forming described even die body in the PET image of even die body, thus obtain the positional information of the pixel of described even die body.Because the pad value of described even die body can obtain, adopt the pad value of the positional information of described pixel and described even die body can generate the decay pattern of described PET transmitting image, and generate without the need to going by CT technology the decay pattern that described PET launches image.Decay drawing generating method described in the embodiment of the present invention, while the corresponding cost of reduction and operational complexity, decreases the radioactive exposure of patient, improves the safety of patient.
Further, the number of launching the pixel of the position of image due to described PET corresponding on even die body can obtain, using the pixel on the PET transmitting image corresponding to the pixel value being more than or equal to first threshold as the first profile region, by adjustment first threshold, and then to obtain on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image, thus using the first profile region corresponding time minimum for described difference as the profile region forming described even die body, other equipment need not be increased, also simpler in operation.
Further, during due to adjustment first threshold, the initial value that the max pixel value on image is first threshold is launched with PET, reduce described first threshold gradually, when the second difference is greater than the first difference, namely using the first profile region corresponding to the first difference as the profile region being formed described even die body, and without the need to adjusting first threshold again and obtaining the difference corresponding to the first threshold after each adjustment, thus make the deterministic process forming even die body profile region more quick, improve formation efficiency.
In addition, can also by directly making the accumulation calculating of pixel number from big to small to the pixel value of the pixel that PET launches on image, the number that corresponding with on even die body for the number of described pixel described PET launches the pixel of the position of image is compared, thus obtains pixel PET image being formed described even die body.Adopt described decay drawing generating method, reduce the complexity in calculating.
Accompanying drawing explanation
Fig. 1 is the flow chart of drawing generating method of decaying in the embodiment of the present invention;
Fig. 2 determines in the embodiment of the present invention that PET launches flow chart image being formed even die body profile region;
Fig. 3 is that the PET of even die body in embodiments of the invention launches image;
Fig. 4 is a kind of first profile region image in the embodiment of the present invention;
Fig. 5 is the another kind first profile region image in the embodiment of the present invention;
Fig. 6 is the decay pattern image in the embodiment of the present invention;
Fig. 7 is the PET image in the embodiment of the present invention;
Fig. 8 is the flow chart determining PET image is formed even die body pixel in the embodiment of the present invention;
Fig. 9 is to the rectangular histogram that the pixel value of the pixel that PET launches on image is added up in the embodiment of the present invention;
Figure 10 is the structural representation of decay pattern generating apparatus in the embodiment of the present invention;
Figure 11 is the structural representation that in the embodiment of the present invention, profile region obtains subelement;
Figure 12 is the structural representation of acquiring unit in the embodiment of the present invention;
Figure 13 is the structural representation of the second acquisition subelement in the embodiment of the present invention.
Detailed description of the invention
Embodiment one
As shown in Figure 1, present embodiments provide a kind of decay drawing generating method of even die body, described method can adopt following steps:
Step 102: determine that the PET of generated even die body launches pixel image being formed described even die body, obtain the positional information of described pixel;
Step 104: the pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
In the process of PET imaging, owing to being subject to the impact of the factor such as Compton effect, scattering, the PET of generation is caused to launch image fault.Because image is all made up of several pixels, to obtain PET image clearly, can correct by launching to PET each pixel image being formed described even die body, the decay pattern namely utilizing the dampening information of each pixel described to form is launched to described PET each pixel image being formed even die body and is corrected.
Wherein, described decay pattern needs to adopt PET to launch positional information image being formed each pixel of even die body, and described even die body pad value generates.Therefore, as long as determine that PET launches pixel image being formed described even die body, and then obtain the positional information of described pixel, recycle described even die body pad value, just can generate the decay pattern that described PET launches image.
In concrete enforcement, described PET launches the transmitting image that image can be described even die body, also can be the transmitting image in a certain cross section of described even die body, that is, the position that on described even die body, corresponding described PET launches image can be whole even die body, also can be a certain cross section of described even die body.
When described PET transmitting image is the transmitting image of whole even die body, by obtaining the volume of described even die body, then according to the corresponding relation between described volume and pixel number, the pixel number that described even die body physical size is corresponding can be obtained.
When described PET transmitting image is the transmitting image in a certain cross section of described even die body, can by obtaining the physical area in described cross section on described even die body, again according to the corresponding relation of described physical area and pixel number, obtain the pixel number that described in described even die body, physical area is corresponding.The pad value of described even die body can be known when described even die body dispatches from the factory.
In concrete enforcement, determine that the PET of generated even die body launches when image forms the pixel of described even die body, first can determine the profile region described PET image being formed described even die body, then the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
As shown in Figure 2, in concrete enforcement, determine the profile region described PET image being formed described even die body, can following steps be adopted:
Step 202: obtain the pixel value that described PET launches pixel on image, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
Step 204: adjust described first threshold, obtains on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Step 206: profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference.
Below in conjunction with Fig. 3 to Fig. 6, to the process determining the profile region described PET image being formed described even die body in Fig. 2, be described in detail.
As shown in Figure 3, for the PET of a certain cylindrical body launches image.The pixel number that on described cylindrical body, corresponding described PET launches the position of image is 20106, and the max pixel value that described PET launches image is 0.0033.Wherein, the pixel value of described pixel is directly proportional to the brightness of described pixel itself, and it is also the pixel that PET launches that on image, brightness is maximum that PET launches the maximum pixel of image pixel value.
Under normal circumstances, PET launches on image, and the pixel at described even die body edge is the point that whole PET launches that on image, brightness is maximum.But in the process generating described cylindrical body PET transmitting image, due to the impact of partial volume effect, the pixel value causing some pixel at described even die body edge is on the low side, is namely less than described max pixel value.Further, some pixel produced because of scattering may be there is in described PET transmitting image.Therefore, if launch the profile region of the pixel on image as described even die body using the PET corresponding to the max pixel value of described even die body, the pixel of the described even die body of thus obtained formation, because lacking some pixel making pixel value on the low side due to partial volume effect, make the result of the pixel of the described even die body of formation finally determined not accurate enough.
In embodiments of the invention, by arranging first threshold, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, adjust described first threshold, to compare on the first profile region corresponding to the first threshold after adjustment and number that in the first profile region, number and the corresponding described PET on described even die body of pixel launches the pixel of the position of image.
When difference between the number that number and the corresponding described PET on described even die body of pixel on the first profile region corresponding to the first threshold after adjusting and in the first profile region launches the pixel of the position of image is minimum, pixel on first profile region and in the first profile region can comprise those pixels making pixel value on the low side because of partial volume effect, Boundary extracting algorithm or other similar approach is taked to obtain the edge of the first profile region afterwards, and the pixel got rid of because the reasons such as scattering produce outside the first profile region, therefore, the first profile region is now closest to the actual profile district of described even die body, decay pattern according to described first profile region generation is described even die body decay pattern the most accurately.
In concrete enforcement, those skilled in the art can arrange according to practical situation and adjust described first threshold.Such as, as shown in Figure 4, when described first threshold equals 5*10 -4time, the white pixel point part in figure is described first profile region.According to statistics, nearly 20664 pixels on the first profile region now and the first profile region, comprise the pixel produced because of reasons such as scatterings, cause described first profile region larger than the actual profile district of cylindrical body in described pixel.
As shown in Figure 5, when described first threshold equals 8*10 -4time, the white pixel point part in figure is described first profile region.According to statistics, nearly 20161 pixels on the first profile region now and the first profile region, substantially eliminate the pixel because the reasons such as scattering produce outside the first profile region, make the actual profile district of described first profile region and cylindrical body comparatively close.First threshold is adopted to equal 8*10 -4time the pixel positional information that obtains and the pad value of described cylindrical shape die body, generate decay pattern as shown in Figure 6.
After carrying out correction for attenuation according to the decay pattern in Fig. 6 to the PET transmitting image in Fig. 1, the PET through correction for attenuation obtained as shown in Figure 7 launches image, and the described PET through correction for attenuation launches image and is PET image.
As can be seen here, the value of described first threshold is different, and the first corresponding profile region is also constantly change.As long as described first threshold value is suitable, just can get rid of PET and launch edge's pixel loss that image is caused by partial volume effect and the pixel produced in the scattering of contour area exopathogenic factor, and then obtain the actual profile district of described cylindrical body, make thus obtained decay pattern also more accurate.After launching image rectification according to described decay pattern to PET, the PET of acquisition launches the profile region of evenly die body on image also more close to the actual profile district of described even die body.
In concrete enforcement, can adopt and adjust described first threshold with the following method: the initial value taking described max pixel value as described first threshold, reduces described first threshold.Such as: each ratio reducing described first threshold is 1% to 5%, and obtains the first difference and the second difference thus, until described second difference is greater than described first difference.
Wherein, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image; Described second difference for current once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image.
Because described first threshold adjusts according to descending order, described difference first reduces along with the reduction of described first threshold, when described difference is decreased to a certain marginal value, then the difference reducing described first threshold and obtain will increase along with the reduction of described first threshold.That is, when the second difference is greater than the first difference, described first difference is the marginal value of described difference.Therefore, when the second difference is greater than the first difference, can stop adjusting described first threshold, profile region image being formed described even die body is launched as described PET in the first profile region corresponding to described first difference.
As shown in Figure 8, launch on image at PET, when the pixel value of the pixel of described even die body center is greater than the pixel value of the pixel produced because of reasons such as scatterings, determine that described PET launches pixel image being formed described even die body and can also comprise the steps:
Step 802: obtain the pixel value that described PET launches pixel on image;
Step 804: the accumulation calculating described pixel value being made from big to small to pixel number, obtains added up a pixel numerical value;
Step 806: when described added up a pixel numerical value equals the pixel number of the position of corresponding described PET transmitting image on described even die body, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
Because described PET launches on image, the pixel value of pixel reduces to center gradually from the edge of even die body.Therefore, when the pixel value of the pixel of described even die body center is greater than the pixel value of the pixel produced because of reasons such as scatterings, launch on image at described PET, the accumulation calculating of pixel number can be made the pixel value of described pixel from big to small.When described added up a pixel numerical value equals the pixel number of the position of corresponding described PET transmitting image on described even die body, can think that the described PET pixel launched on image corresponding to pixel value of all accumulated calculating is described PET and launches pixel image being formed described even die body.
In concrete enforcement, for convenience of the accumulation calculating described pixel value being made from big to small to pixel number, can add up described pixel value, namely the distribution situation of described pixel value be added up.Histogrammic mode such as can be adopted to add up described pixel value, other modes also can be adopted to add up, as long as visual representation the distribution situation of each pixel value can be gone out.According to described histogrammic statistical result, more described pixel value is made from big to small to the accumulation calculating of pixel number.
Such as: as shown in Figure 9, the statistics with histogram situation that the pixel value launching pixel on image to the PET of a certain die body is done.Wherein, transverse axis represents pixel value, from left to right increases gradually.The longitudinal axis represents the number of pixel, increases gradually from top to bottom.The pixel number that on known described die body, corresponding described PET launches the position of image is a.According to step 802 to step 806, pixel value in Fig. 9 can be more than or equal to 0.65*10 -3pixel launch pixel image being formed described even die body as described PET.
Embodiment two
As shown in FIG. 6 and 7, wherein Fig. 6 is the decay pattern adopting the first profile region determined in Fig. 5 to generate, adopt the decay pattern in Fig. 6, generate even die body PET image as shown in Figure 7, that is: the PET generating even die body launches image, and adopt the decay pattern that the decay drawing generating method described in above-described embodiment generates, image is launched to described PET and carries out correction for attenuation, the PET image through correction for attenuation can be generated.
For making those skilled in the art understand better and realize the present invention, device corresponding to the generation method of above-mentioned decay pattern and PET image is below described in detail.
Embodiment three
As shown in Figure 10, present embodiments provide a kind of generating apparatus of even die body decay pattern, described device can comprise:
Acquiring unit 1002, for determining that the PET of generated even die body launches pixel image being formed described even die body, obtains the positional information of described pixel;
Generation unit 1004, for the pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
Wherein, PET launches image and is made up of several pixels, and the even die body in described PET image also comprises several pixels.Described PET is launched image carries out correction for attenuation, namely pixel image being formed described even die body is launched to PET and carry out correction for attenuation.
Carry out correction for attenuation to the pixel forming described even die body, correspondingly will obtain the dampening information of each pixel of the even die body of formation, described dampening information comprises the positional information and pad value thereof that form described even die body pixel.Because the pad value of even die body can obtain, therefore, after described acquiring unit 1002 obtains the positional information of the pixel described PET image being formed described even die body, the pad value of the described positional information that described generation unit 1004 just can obtain according to described acquiring unit 1002 and described even die body, generates the decay pattern representing and form described even die body dampening information.
In concrete enforcement, in order to obtain the positional information of the pixel forming described even die body, described acquiring unit can comprise: profile region obtains subelement (not shown), for determining that described PET launches profile region image being formed described even die body, the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
As shown in figure 11, in concrete enforcement, described profile region obtains subelement and can comprise:
First profile region acquisition module 1102, the pixel value of pixel on image is launched for obtaining described PET, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
First threshold adjusting module 1104, adjust described first threshold, to obtain on described first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Second profile region acquisition module 1106, launches as described PET profile region image being formed described even die body using the first profile region corresponding time minimum for described difference.
In the present embodiment, by arranging first threshold, the described PET corresponding to the pixel value being more than or equal to first threshold is launched pixel on image as the first profile region by described first profile region acquisition module 1102.Described first threshold adjusting module 1104 adjusts described first threshold, to compare on the first profile region corresponding to the first threshold after adjustment and number that in the first profile region, number and the corresponding described PET on described even die body of pixel launches the pixel of the position of image.When difference between the number that number and the corresponding described PET on described even die body of pixel on the first profile region corresponding to the first threshold after adjusting and in the first profile region launches the pixel of the position of image is minimum, pixel on described first profile region and in the first profile region can comprise those pixels making pixel value on the low side because of partial volume effect, and the pixel got rid of because the reasons such as scattering produce, the first profile region is now closest to the actual profile district of described even die body.Profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference by described profile region acquisition module 1106.
In concrete enforcement, described first threshold adjusting module 1104 may be used for, using the initial value of described max pixel value as described first threshold, reducing described first threshold, such as: each ratio reducing described first threshold is 1% to 5%.Obtain the first difference and the second difference, until described first difference is greater than described second difference.
Wherein, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image on first profile region of described second difference corresponding to the first threshold after current once reduction and in the first profile region.
Now, described profile region acquisition module 1106 is for launching as described PET profile region image being formed described even die body using the first profile region corresponding to described first difference.Adopt the pad value of described first profile region and described even die body can generate decay pattern as described in Figure 6.
Because described first threshold adjusts according to descending order, described difference first reduces along with the reduction of described first threshold, when described difference is decreased to a certain marginal value, then the difference reducing described first threshold and obtain will increase along with the reduction of described first threshold.That is, when the second difference is greater than the first difference, described first difference is the marginal value of described difference.Therefore, when the second difference is greater than the first difference, can stop adjusting described first threshold, profile region image being formed described even die body is launched as described PET in the first profile region corresponding to described first difference.
As shown in figure 12, launch on image at PET, when the pixel value of the pixel of described even die body center is greater than the pixel value of the pixel produced because of reasons such as scatterings, in order to obtain the positional information of the pixel forming described even die body, described acquiring unit can also comprise:
First obtains subelement 1202, launches the pixel value of pixel on image for obtaining described PET;
Second obtains subelement 1204, for making the accumulation calculating of pixel number from big to small to described pixel value, obtains added up a pixel numerical value;
3rd obtains subelement 1206, when described added up a pixel numerical value equals the pixel number of the position of corresponding described PET transmitting image on described even die body, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
In order to determine that described PET launches pixel image being formed described even die body more rapidly, obtain the positional information of described pixel, in concrete enforcement, as shown in figure 13, described second obtains subelement 1204 can comprise:
Statistics with histogram module 1302, for making statistics with histogram to described pixel value, obtains statistics with histogram result;
Computing module 1304, for the histogrammic statistical result obtained according to described statistics with histogram module, makes the accumulation calculating of pixel number from big to small to described pixel value.
It should be noted that, described second obtains subelement when taking statistics to described pixel value, and other forms can also be adopted to add up, as long as visual representation can go out the distribution situation of described pixel value.
Embodiment four
Embodiments of the invention additionally provide the generating apparatus of even die body PET image, the decay pattern that the generating apparatus of described PET image can adopt the generating apparatus of above-mentioned decay pattern to generate, and generate corresponding PET image.Wherein, the image of described decay pattern can with reference to Fig. 6, and described PET image can with reference to Fig. 7.
In concrete enforcement, the generating apparatus of described PET image can comprise the first generation unit and the second generation unit.Wherein, described second generation unit adopts the decay pattern generated in above-described embodiment, carries out correction for attenuation, generate the PET image through correction for attenuation to the PET transmitting image that described first generation unit generates.
One of ordinary skill in the art will appreciate that all or part of step in the various methods of above-described embodiment is that the hardware that can carry out instruction relevant by program has come, this program can be stored in a computer-readable recording medium, and storage medium can comprise: ROM, RAM, disk or CD etc.
Although the present invention discloses as above, the present invention is not defined in this.Any those skilled in the art, without departing from the spirit and scope of the present invention, all can make various changes or modifications, and therefore protection scope of the present invention should be as the criterion with claim limited range.

Claims (15)

1. an even die body decay map generalization method, is characterized in that, comprising:
Determine that the positron emission fault PET of generated even die body launches pixel image being formed described even die body, obtain the positional information of described pixel;
The pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
2. evenly die body decay drawing generating method as claimed in claim 1, is characterized in that, the described PET transmitting image determining generated even die body is formed the pixel of described even die body, comprising:
Determine that described PET launches profile region image being formed described even die body, the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
3. evenly die body decay drawing generating method as claimed in claim 2, is characterized in that, describedly determines that described PET launches profile region image being formed described even die body, comprising:
Obtain the pixel value that described PET launches pixel on image, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
Adjust described first threshold, to obtain on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference.
4. evenly die body decay drawing generating method as claimed in claim 3, it is characterized in that, the described first threshold of described adjustment, to obtain on the first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image, comprise: the initial value taking described max pixel value as described first threshold, reduce described first threshold, obtain the first difference and the second difference, until described second difference is greater than described first difference, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image, described second difference for current once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image,
Profile region image being formed described even die body is launched as described PET in the first profile region corresponding time minimum for described difference, comprising: profile region image being formed described even die body is launched as described PET in the first profile region corresponding to described first difference.
5. evenly die body decay drawing generating method as claimed in claim 4, it is characterized in that, the described first threshold of described reduction comprises: each ratio reducing described first threshold is 1% to 5%.
6. evenly die body decay drawing generating method as claimed in claim 1, is characterized in that, describedly determines that described PET launches pixel image being formed described even die body, comprising:
Obtain the pixel value that described PET launches pixel on image;
Described pixel value is made from big to small to the accumulation calculating of pixel number, obtain added up a pixel numerical value;
When described added up a pixel numerical value equals the pixel number of the position of corresponding described PET transmitting image on described even die body, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
7. evenly die body decay drawing generating method as claimed in claim 6, it is characterized in that, the described accumulation calculating described pixel value being made from big to small to pixel number, comprising:
Statistics with histogram is done to described pixel value, according to described histogrammic statistical result, described pixel value is made from big to small to the accumulation calculating of pixel number.
8. a generation method for even die body PET image, is characterized in that, comprising:
The PET generating even die body launches image;
Adopt the decay pattern that the decay drawing generating method described in any one of claim 1 to 7 generates, image is launched to described PET and carries out correction for attenuation, generate the PET image through correction for attenuation.
9. a generating apparatus for even die body decay pattern, is characterized in that, comprising:
Acquiring unit, for determining that the PET of generated even die body launches pixel image being formed described even die body, obtains the positional information of described pixel;
Generation unit, for the pad value of described even die body adopting the positional information of described pixel and obtain, generates the decay pattern that described PET launches image.
10. the device of evenly die body decay map generalization as claimed in claim 9, it is characterized in that, described acquiring unit comprises:
Profile region obtains subelement, for determining that described PET launches profile region image being formed described even die body, the pixel on described profile region and in profile region is launched as described PET pixel image being formed described even die body.
The device of 11. evenly die body decay map generalizations as claimed in claim 10, it is characterized in that, described profile region obtains subelement and comprises:
First profile region acquisition module, the pixel value of pixel on image is launched for obtaining described PET, described pixel value comprises max pixel value, using the pixel on the described PET transmitting image corresponding to the pixel value being more than or equal to default first threshold as the first profile region, described first threshold is the pixel value being less than or equal to described max pixel value;
First threshold adjusting module, for adjusting described first threshold, to obtain on described first profile region corresponding to first threshold after adjustment and in the first profile region, number and the corresponding described PET on described even die body of pixel launches the difference between the number of the pixel of the position of image;
Second profile region acquisition module, for launching as described PET profile region image being formed described even die body using the first profile region corresponding time minimum for described difference.
The device of 12. evenly die body decay map generalizations as claimed in claim 11, it is characterized in that, described first threshold adjusting module is used for the initial value of described max pixel value as described first threshold, reduce described first threshold, obtain the first difference and the second difference, until described first difference is greater than described second difference, described first difference for before once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image, described second difference for current once reduce after first threshold corresponding to the first profile region on and in the first profile region number and the corresponding described PET on described even die body of pixel launch the difference between the number of the pixel of the position of image,
Described profile region acquisition module is for launching as described PET profile region image being formed described even die body using the first profile region corresponding to described first difference.
The device of 13. evenly die body decay map generalizations as claimed in claim 9, it is characterized in that, described acquiring unit comprises:
First obtains subelement, launches the pixel value of pixel on image for obtaining described PET;
Second obtains subelement, for making the accumulation calculating of pixel number from big to small to described pixel value, obtains added up a pixel numerical value;
3rd obtains subelement, during for equaling the pixel number of the position of corresponding described PET transmitting image on described even die body when described added up a pixel numerical value, the pixel that the described PET corresponding to the pixel value of accumulated calculating launches on image is launched as described PET pixel image being formed described even die body.
The device of 14. evenly die body decay map generalizations as claimed in claim 13, it is characterized in that, described second obtains subelement comprises:
Statistics with histogram module, for making statistics with histogram to described pixel value, obtains statistics with histogram result;
Computing module, for the histogrammic statistical result obtained according to described statistics with histogram module, makes the accumulation calculating of pixel number from big to small to described pixel value.
The generating apparatus of 15. 1 kinds of even die body PET image, is characterized in that, comprising:
First generation unit, described first generation unit launches image for generating PET;
Second generation unit, for the decay pattern adopting the decay pattern generating apparatus described in any one of claim 9 to 14 to generate, launches image to described PET and carries out correction for attenuation, generate the PET image through correction for attenuation.
CN201310533167.XA 2013-10-30 2013-10-30 Decay pattern and positron emission fault image generating method and device Active CN104586416B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310533167.XA CN104586416B (en) 2013-10-30 2013-10-30 Decay pattern and positron emission fault image generating method and device
CN201710224826.XA CN106963406B (en) 2013-10-30 2013-10-30 Attenuation map and positron emission tomography image generation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310533167.XA CN104586416B (en) 2013-10-30 2013-10-30 Decay pattern and positron emission fault image generating method and device

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201710224826.XA Division CN106963406B (en) 2013-10-30 2013-10-30 Attenuation map and positron emission tomography image generation device

Publications (2)

Publication Number Publication Date
CN104586416A true CN104586416A (en) 2015-05-06
CN104586416B CN104586416B (en) 2017-07-04

Family

ID=53112715

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201310533167.XA Active CN104586416B (en) 2013-10-30 2013-10-30 Decay pattern and positron emission fault image generating method and device
CN201710224826.XA Active CN106963406B (en) 2013-10-30 2013-10-30 Attenuation map and positron emission tomography image generation device

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201710224826.XA Active CN106963406B (en) 2013-10-30 2013-10-30 Attenuation map and positron emission tomography image generation device

Country Status (1)

Country Link
CN (2) CN104586416B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106333701A (en) * 2016-09-30 2017-01-18 上海联影医疗科技有限公司 Active die body positioning method utilizing positive ion emission tomography system
US9990718B2 (en) 2016-02-29 2018-06-05 Shanghai United Imaging Healthcare Co., Ltd. System and method for reconstructing ECT image

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111544022B (en) * 2020-05-13 2023-11-28 上海联影医疗科技股份有限公司 PET attenuation correction method and device and computer equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009047602A (en) * 2007-08-21 2009-03-05 Toshiba Corp Positron emission computerd tomograph, attenuation map creating device, and attenuation map creating program
CN101953693A (en) * 2009-07-16 2011-01-26 西门子公司 The estimation based on model of complete or part pet decay pattern
US7923691B2 (en) * 2008-05-08 2011-04-12 Siemens Aktiengesellschaft Method for producing an attenuation map
CN102293662A (en) * 2010-06-17 2011-12-28 西门子公司 Method For Determining Radiation Attenuation In A Positron Emission Tomography Scanner
CN103186882A (en) * 2011-12-30 2013-07-03 沈阳东软派斯通医疗系统有限公司 Image attenuation correction method and image attenuation correction device in position emission computed tomography (PET) system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2210238B1 (en) * 2007-11-09 2013-03-20 Koninklijke Philips Electronics N.V. Apparatus and method for generation of attenuation map
US8611628B2 (en) * 2008-05-15 2013-12-17 Koninklijke Philips N.V. Using non-attenuation corrected PET emission images to compensate for incomplete anatomic images

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009047602A (en) * 2007-08-21 2009-03-05 Toshiba Corp Positron emission computerd tomograph, attenuation map creating device, and attenuation map creating program
US7923691B2 (en) * 2008-05-08 2011-04-12 Siemens Aktiengesellschaft Method for producing an attenuation map
CN101953693A (en) * 2009-07-16 2011-01-26 西门子公司 The estimation based on model of complete or part pet decay pattern
CN102293662A (en) * 2010-06-17 2011-12-28 西门子公司 Method For Determining Radiation Attenuation In A Positron Emission Tomography Scanner
CN103186882A (en) * 2011-12-30 2013-07-03 沈阳东软派斯通医疗系统有限公司 Image attenuation correction method and image attenuation correction device in position emission computed tomography (PET) system

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9990718B2 (en) 2016-02-29 2018-06-05 Shanghai United Imaging Healthcare Co., Ltd. System and method for reconstructing ECT image
US10535166B2 (en) 2016-02-29 2020-01-14 Shanghai United Imaging Healthcare Co., Ltd. System and method for reconstructing ECT image
US11557067B2 (en) 2016-02-29 2023-01-17 Shanghai United Imaging Healthcare Co., Ltd. System and method for reconstructing ECT image
CN106333701A (en) * 2016-09-30 2017-01-18 上海联影医疗科技有限公司 Active die body positioning method utilizing positive ion emission tomography system
CN106333701B (en) * 2016-09-30 2019-05-31 上海联影医疗科技有限公司 The active die body localization method of positron emission tomography

Also Published As

Publication number Publication date
CN104586416B (en) 2017-07-04
CN106963406B (en) 2020-06-02
CN106963406A (en) 2017-07-21

Similar Documents

Publication Publication Date Title
CN102934143B (en) For producing the method for decay pattern in PET-MR
EP3226766B1 (en) System and method for image calibration
US8600139B2 (en) Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information
US11179128B2 (en) Methods and systems for motion detection in positron emission tomography
CN111436958B (en) CT image generation method for PET image attenuation correction
US11918390B2 (en) Methods and systems for motion detection in positron emission tomography
US9569842B2 (en) Medical imaging apparatus
CN103400361A (en) Threshold-based gain correction algorithm in oral computed tomography (CT)
CN104586416A (en) Method and device for generating attenuation diagrams and positron emission cross-sectional images
US20130101193A1 (en) Positron Emission Tomography and Method for Correcting Attenuation of PET Image Using Magnetic Resonance Image
CN105678750A (en) Gray scale mapping curve generation method and apparatus for medical images
KR102283934B1 (en) Method and system for generating virtual CT(Computed Tomography) image and attenuation-corrected PET(Positron Emission Tomography) image based on PET image
US7242004B2 (en) Image correction method, image correction apparatus, and image correction program
Preuhs et al. Over-exposure correction in CT using optimization-based multiple cylinder fitting
US20080073538A1 (en) Application-driven optimization of acquisition and reconstruction of SPECT/PET projection data
CN103405241B (en) Detector afterglow correction method for ray imaging
KR20140043637A (en) Method and apparatus for generating system response of scanner of imaging apparatus and medical image using the same
KR101350496B1 (en) Method to generate a attenuation map of emission tomography and MRI combined imaging system
US20220414832A1 (en) X-ray imaging restoration using deep learning algorithms
Chan et al. Simultaneous partial volume correction and noise regularization for cardiac SPECT/CT
CN110264537B (en) PET image reconstruction method, system, readable storage medium and apparatus
US20190142358A1 (en) Method And System For Dose-Less Attenuation Correction For PET And SPECT
CN109272486A (en) Training method, device, equipment and the storage medium of MR image prediction model
EP3138495B1 (en) Individualized imaging method
CN112634147B (en) PET image noise reduction method, system, device and medium for self-supervision learning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee before: Shanghai Lianying Medical Technology Co.,Ltd.