CN106556857B - Response location localization method and device - Google Patents

Abstract

The present invention provides a kind of response location localization method and device.This method comprises: receiving detected by least one photosensor arrays for being coupled with scintillation crystal, gamma photons and scintillation crystal react the energy information of generated optical photon, wherein, scintillation crystal is integrated scintillation crystal, and scintillation crystal has through-hole, through-hole is for accommodating object to be imaged, and gamma photons in the esoteric positron annihilation effect of object to be imaged by generating；And according to received energy information determine response location of the gamma photons in scintillation crystal, wherein response location using cylindrical coordinate system indicate.According to embodiments of the present invention, there is provided a kind of Positron emission tomography equipment for the problem that using integrated scintillation crystal, positioned for the response location to gamma photons by the way of, above-mentioned improved Positron emission tomography equipment can not be suitable for by solving traditional position calculating method.

Description

Response location localization method and device

Technical field

The present invention relates to Positron emission tomography fields, and in particular, to a kind of response location localization method and device.

Background technique

The full name of Positron emission tomography is positron e mission computed tomography (Positron Emission Computed Tomography, abbreviation PET), it is a kind of shown using Radioactive isotope method in human body or animal body The technology of portion's structure is the main means of nuclear medicine studies and clinical diagnosis.

In traditional Positron emission tomography equipment, generally by multiple square detectors (including scintillation crystal, photoelectric transfer The devices such as sensor) it is circularized or the detector system of almost spherical by mechanical structure splicing, for detecting gamma photons.Tool Body, each scintillation crystal (or scintillation crystal array) of traditional Positron emission tomography equipment is coupled with photoelectric sensor Independent detector is formed, and multiple independent detectors are stitched together by complicated mechanical structure and form positive electron hair Penetrate the detector system of imaging device.Since the assembling of detector is spliced, traditional Positron emission tomography equipment is caused to exist Spatial resolution is low, sensitivity is low, edge effect is strong, Machine Design difficulty problem.Therefore, it proposes at present a kind of improved Positron emission tomography equipment uses integrated scintillation crystal to realize, avoid the above-mentioned assembling splicing due to detector and The problem of causing.For this improved Positron emission tomography equipment, traditional position calculating side can not be simply used Method determines response location of the gamma photons in scintillation crystal.

Accordingly, it is desirable to provide a kind of new method for being positioned to the response location of gamma photons, at least portion It solves the above problems with dividing.

Summary of the invention

In order at least be partially solved problems of the prior art, according to an aspect of the present invention, one kind is provided Response location localization method.It is detected this method comprises: receiving at least one photosensor arrays coupled with scintillation crystal To, gamma photons and scintillation crystal react caused by optical photon energy information, wherein scintillation crystal is one The scintillation crystal of change, and scintillation crystal has through-hole, and for accommodating object to be imaged, gamma photons pass through to be imaged through-hole The esoteric positron annihilation effect of object generates；And determine that gamma photons are brilliant in flashing according to the received energy information of institute Intracorporal response location, wherein response location is indicated using cylindrical coordinate system.

According to another aspect of the present invention, a kind of response location positioning device is provided.The device includes: the first reception mould Block, for receiving detected by least one photosensor arrays coupled with scintillation crystal, gamma photons and flashing it is brilliant Body react caused by optical photon energy information, wherein scintillation crystal is integrated scintillation crystal, and is flashed Crystal has through-hole, and for accommodating object to be imaged, gamma photons pass through in object to be imaged esoteric positive electron through-hole Bury in oblivion effect generation；And position determination module, for determining gamma photons in scintillation crystal according to the received energy information of institute Interior response location, wherein response location is indicated using cylindrical coordinate system.

According to the method for the embodiment of the present invention and device, can provide it is a kind of for using integrated scintillation crystal just Electron emission imaging device, the mode that is positioned for the response location to gamma photons, to solve traditional position Calculation method can not be suitable for the problem of above-mentioned improved Positron emission tomography equipment.

A series of concept of simplification is introduced in summary of the invention, these concepts will be in the detailed description section into one Step is described in detail.This part of the disclosure be not meant to attempt to limit technical solution claimed key feature and Essential features do not mean that the protection scope for attempting to determine technical solution claimed more.

Below in conjunction with attached drawing, the advantages of the present invention will be described in detail and feature.

Detailed description of the invention

Following drawings of the invention is incorporated herein as part of the present invention for the purpose of understanding the present invention.Shown in the drawings of this hair Bright embodiment and its description, principle used to explain the present invention.In the accompanying drawings,

Fig. 1 a shows the scintillation crystal and light of improved Positron emission tomography equipment according to an embodiment of the invention The schematic diagram of electric transducer array；

Fig. 1 b show improved Positron emission tomography equipment in accordance with another embodiment of the present invention scintillation crystal and The schematic diagram of photosensor arrays；

Fig. 1 c show according to the present invention the scintillation crystal of the improved Positron emission tomography equipment of another embodiment and The schematic diagram of photosensor arrays；

Fig. 1 d show according to the present invention the scintillation crystal of the improved Positron emission tomography equipment of further embodiment and The schematic diagram of photosensor arrays；

Fig. 2 shows the schematic diagrames of improved Positron emission tomography equipment according to an embodiment of the invention；

Fig. 3 shows the schematic block diagram of response location localization method according to an embodiment of the invention；

Fig. 4 shows the schematic diagram of the cylindrical coordinate system according to an embodiment of the invention applied to hexa-prism crystal；

Fig. 5 shows the flow diagram of the step of angular coordinate of determining response location according to an embodiment of the invention；

Fig. 6 is shown will be shown in Fig. 4, with the photosensor arrays expansion acquisition later of the side surface coupling of scintillation crystal The arrangement schematic diagram of photosensor arrays；

Fig. 7 shows according to an embodiment of the invention for obtaining the structure of the test device of angular error distribution function Schematic diagram；

Fig. 8 shows the process signal of the step of radial coordinate of determining response location according to an embodiment of the invention Figure；

Fig. 9 a to 9i be shown respectively it is according to an embodiment of the invention under different gammaphoton reaction depths to photoelectricity The distribution of photons of optical photon detected by sensor carries out the fitting result of Gauss curve fitting；

Figure 10 shows according to an embodiment of the invention for obtaining the structural schematic diagram of the test device of look-up table；

The process that Figure 11 shows the step of radial coordinate of determining response location in accordance with another embodiment of the present invention is shown It is intended to；

Figure 12 shows the photoelectric sensor battle array that basis according to an embodiment of the invention is coupled with the bottom surface of scintillation crystal The energy information of the detected optical photon of column determines the schematic diagram of radial coordinate；

Figure 13 shows the process signal of the step of axial coordinate of determining response location according to an embodiment of the invention Figure；

Figure 14 shows according to an embodiment of the invention for obtaining the knot of the test device of axial error distribution function Structure schematic diagram；And

Figure 15 shows the schematic block diagram of response location positioning device according to an embodiment of the invention.

Specific embodiment

In the following description, a large amount of details is provided so as to thoroughly understand the present invention.However, this field skill Art personnel will be seen that, only relate to presently preferred embodiments of the present invention described below, and the present invention may not need one or more in this way Details and be carried out.In addition, in order to avoid confusion with the present invention, not for some technical characteristics well known in the art It is described.

To solve the above-mentioned problems, the present invention proposes a kind of response location localization method.This method is suitable for using one The improved Positron emission tomography equipment of the scintillation crystal of change.The improved Positron emission tomography equipment, which uses, has one The scintillation crystal of body structure and the photosensor arrays being adapted to the scintillation crystal form detector, can solve in tradition Positron emission tomography equipment in due to detector splices and the problem of cause.

In order to explain the present invention, the Positron emission tomography equipment that method described herein is applicable in is described first below Structure.

Positron emission tomography equipment may include scintillation crystal, at least one photosensor arrays, reading circuit and Data processing module.

Scintillation crystal has integral structure.Illustratively, scintillation crystal can be in the fabrication process by one at Type technology directly obtains, without any type of splicing and assembling.Illustratively, scintillation crystal can be structured as cylindric Crystal or polygon prism shape crystal.The center of scintillation crystal has a through-hole, can run through two bottom surfaces of scintillation crystal.It is logical Hole for accommodating object to be imaged (toy etc.), i.e., for object to be imaged by the through-hole enter positron emission at As equipment.Illustratively, through-hole can be structured as polygonal through hole or circular through hole.

At least one photosensor arrays is coupled with scintillation crystal, is occurred instead for detecting gamma photons and scintillation crystal Answer generated optical photon.Illustratively, at least one photosensor arrays can be with scintillation crystal direct-coupling or logical Cross optical glue coupling.Gamma photons in the esoteric positron annihilation effect of object to be imaged by generating.Specifically, exist When scanning object to be imaged using Positron emission tomography equipment, the same position of radioactivity can be contained to object internal injection to be imaged The tracer of element.It can bury in oblivion when the positive electron that isotope is released meets with the intracorporal negative electron of object to be imaged, thus Generation a pair of contrary (difference 180 degree), energy are the gamma photons of 511KeV.The contrary gamma of a pair of generation Photon is incident on respectively in two opposite positions in scintillation crystal.The gamma photons being incident in scintillation crystal can be with sudden strain of a muscle Bright crystal reacts, and thus generates a large amount of optical photons.The photosensor arrays coupled with scintillation crystal can detecte this A little optical photons, and when it detects optical photon, it can will be seen that the optical signal of photon is converted to electric signal and will turn Change the electric signal output of acquisition.

Scintillation crystal can be any suitable crystal, and the present invention limits not to this.For example, scintillation crystal can be Bismuth germanium oxide (BGO), yttrium luetcium silicate (LYSO) or lanthanum bromide (LaBr3) etc..Photosensor arrays can be by any suitable The array of photoelectric sensor composition, such as photomultiplier tube (PMT), silicon photomultiplier (SiPM) or avalanche photodide (APD) etc..

As described above, scintillation crystal can be polygon prism shape crystal or cylinder crystal, and through-hole can be circular through hole Or polygonal through hole.Illustratively, in the case where scintillation crystal is polygon prism shape crystal, photosensor arrays can be covered The side of polygon prism shape crystal.In the case where the through-hole of scintillation crystal is polygonal through hole, photosensor arrays can cover The inner wall of lid polygonal through hole.No matter scintillation crystal is polygon prism shape crystal or cylinder crystal, and photosensor arrays are all The upper bottom surface and/or bottom surface of scintillation crystal can be covered.It should be noted that being polygon prism shape crystal or cylindric crystalline substance in scintillation crystal In the case where body, scintillation crystal tool is there are two bottom surface, and wherein which bottom surface is upper bottom surface and which bottom surface is that bottom surface can be with It is set as needed, the present invention limits not to this.

Compared with traditional Positron emission tomography equipment, above-mentioned improved Positron emission tomography equipment has gamma light Sub- positioning accuracy height, the advantages such as high sensitivity, edge effect is weak, Machine Design difficulty is low.

The flashing for describing improved Positron emission tomography equipment according to an embodiment of the present invention below with reference to Fig. 1 a-1d is brilliant Several examples of body and photosensor arrays.

According to an embodiment of the present invention, scintillation crystal is polygon prism shape crystal, and through-hole is circular through hole, and through-hole is through sudden strain of a muscle Two bottom surfaces of bright crystal.Fig. 1 a shows the flashing of improved Positron emission tomography equipment according to an embodiment of the invention The schematic diagram of crystal and photosensor arrays.As shown in Figure 1a, scintillation crystal is polygon prism shape crystal, and through-hole is round logical Hole.In such a case, it is possible to which the inner wall in circular through hole is coated with highly-reflective coating or highly reflecting films.In addition, illustratively, dodging Each side of bright crystal can be coupled with one or more photosensor arrays.Scintillation crystal shown in Fig. 1 a is hexagonal prisms Shape crystal, there are six sides for tool, and each side is coupled with a photosensor arrays.In addition, implementing shown in Fig. 1 a In example, two bottom surfaces of scintillation crystal are coupled with six photosensor arrays respectively.It will be appreciated, however, that shown in Fig. 1 a The configuration mode of photosensor arrays is only exemplary rather than limitation.For example, two bottom surfaces of scintillation crystal can not be with photoelectricity Sensor array coupling, or only a certain bottom surface is coupled with one or more photosensor arrays.

According to a further embodiment of the invention, scintillation crystal is cylinder crystal, and through-hole is polygonal through hole, and through-hole runs through Two bottom surfaces of scintillation crystal.Fig. 1 b shows improved Positron emission tomography equipment in accordance with another embodiment of the present invention The schematic diagram of scintillation crystal and photosensor arrays.As shown in Figure 1 b, scintillation crystal is cylinder crystal, and through-hole is polygon Through-hole.In such a case, it is possible to be coated with highly-reflective coating or highly reflecting films in the side of scintillation crystal.In addition, illustratively, Each inner wall of through-hole can be coupled with one or more photosensor arrays.Scintillation crystal shown in Fig. 1 b is cylindric crystalline substance Body, through-hole are hexagon through-hole, and there are six inner wall, each inner wall couples through-hole tool with a photosensor arrays.This Outside, in the embodiment shown in Fig. 1 b, two bottom surfaces of scintillation crystal are coupled with six photosensor arrays respectively.However, It should be appreciated that the configuration mode of photosensor arrays shown in Fig. 1 b is only exemplary rather than limitation.For example, the two of scintillation crystal A bottom surface can not be coupled with photosensor arrays, or only a certain bottom surface and one or more photosensor arrays couplings It closes.

Another embodiment according to the present invention, scintillation crystal are polygon prism shape crystal, and through-hole is polygonal through hole, and through-hole passes through Wear two bottom surfaces of scintillation crystal.Fig. 1 c shows the improved Positron emission tomography equipment of another embodiment according to the present invention Scintillation crystal and photosensor arrays schematic diagram.As illustrated in figure 1 c, scintillation crystal is polygon prism shape crystal, and through-hole is more Side shape through-hole.In this case, illustratively, each side of scintillation crystal can be with one or more photoelectric sensor battle arrays Column coupling, and/or each inner wall of through-hole can be coupled with one or more photosensor arrays.It is flashed shown in Fig. 1 c brilliant Body is hexa-prism crystal, and there are six sides for tool, and each side is coupled with a photosensor arrays.In addition, Fig. 1 c institute The through-hole of the scintillation crystal shown is hexagon through-hole, and there are six inner wall, each inner wall and a photoelectric sensor battle arrays for through-hole tool Column coupling.In addition, in the embodiment shown in Fig. 1 c, two bottom surfaces of scintillation crystal respectively with six photosensor arrays couplings It closes.It will be appreciated, however, that the configuration mode of photosensor arrays shown in Fig. 1 c is only exemplary rather than limitation.For example, flashing Two bottom surfaces of crystal can not be coupled with photosensor arrays, or only a certain bottom surface and one or more photoelectric sensors Array coupling.

Further embodiment according to the present invention, scintillation crystal are cylinder crystals, and through-hole is circular through hole, and through-hole is through sudden strain of a muscle Two bottom surfaces of bright crystal.Fig. 1 d shows the detection for Positron emission tomography equipment of further embodiment according to the present invention The schematic diagram of device.As shown in Figure 1 d, scintillation crystal is cylinder crystal, and through-hole is circular through hole.In this case, exemplary Ground can be coated with highly-reflective coating or highly reflecting films in the side of scintillation crystal, and the inner wall of through-hole can be with multiple photoelectricity Sensor array coupling.Scintillation crystal shown in Fig. 1 d be cylinder crystal, through-hole be circular through hole, the inner wall of the through-hole with Multiple photosensor arrays couplings.It, illustratively, can will be along row's photoelectricity of the axial alignment of scintillation crystal in Fig. 1 d Sensor is considered as a photosensor arrays.The number of the photosensor arrays coupled with the inner wall of through-hole can be any Suitable number (such as 24), the present invention limits not to this.Comparison it is appreciated that photosensor arrays as far as possible It is covered with the inner wall of through-hole, i.e. the fewer uncovered region on the inner wall of through-hole the better.Certainly, the case where technology allows Under, the inner wall of through-hole can be coupled with a photosensor arrays.It is understood that since the inner wall of through-hole is round , therefore for each photosensor arrays, it can not be bonded the inner wall of through-hole completely, there is sky among the two Gap.Optical glue can be filled in the gap among the inner wall of photosensor arrays and through-hole, so that optical photon can Enter photosensor arrays from scintillation crystal by optical glue.It is understood that the photoelectricity coupled with the inner wall of through-hole The number of sensor array is bigger, and being bonded between each photosensor arrays and the inner wall of through-hole is closer, i.e., in the two Between gap it is smaller.In fact, the gap is generally smaller for improved Positron emission tomography equipment, exist substantially Millimeter magnitude, such as 0.2 millimeter or so, therefore the detection of gamma photons is influenced little.

In addition, in the embodiment shown in Fig. 1 d, two bottom surfaces of scintillation crystal respectively with seven photosensor arrays Coupling.It will be appreciated, however, that the configuration mode of photosensor arrays shown in Fig. 1 d is only exemplary rather than limitation.For example, dodging Two bottom surfaces of bright crystal can not be coupled with photosensor arrays, or only a certain bottom surface and one or more photoelectric sensings The coupling of device array.

According to above description it is found that photosensor arrays can not only be coupled with plane, (ginseng can also be coupled with curved surface Examine Fig. 1 d).Therefore, skilled person will appreciate that, the inner wall of the through-hole of scintillation crystal shown in Fig. 1 a equally can be with one Or multiple photosensor arrays couplings, specific implementation can refer to Fig. 1 d and associated description, repeat no more.

Although scintillation crystal is shown as hexa-prism structure, and in Fig. 1 b and Fig. 1 c, through-hole in Fig. 1 a and Fig. 1 c Be shown as hexagonal structure, it is noted that scintillation crystal be polygon prism shape crystal in the case where scintillation crystal prism number And through-hole is that the number of edges of through-hole in the case where polygonal through hole may each be any suitable number, the present invention is not to this progress Limitation.For example, scintillation crystal can be triangular prism shape crystal, hexa-prism crystal, ten prism-shaped crystal, even 20 quadrangulars Shape crystal, etc..Similarly, through-hole can be tetragonal through hole, hexagon through-hole, 20 tetragonal through hole, etc..It should infuse Meaning, in the case where scintillation crystal is polygon prism shape crystal and through-hole is polygonal through hole, the prism number of scintillation crystal and logical The number of edges in hole may be the same or different.

Reading circuit is connect at least one photosensor arrays, defeated for receiving at least one photosensor arrays Electric signal out, and the energy information of optical photon detected by least one photosensor arrays is exported, electric signal is The optical signal of the optical photon detected by least one photosensor arrays to it, which convert, to be obtained.At data Reason module is connect with reading circuit.Data processing module may include reaction described herein location position device, for realizing Reaction described herein location positioning method.

Fig. 2 shows the schematic diagrames of improved Positron emission tomography equipment according to an embodiment of the invention.Fig. 2 is court The plan view observed to a bottom surface of the scintillation crystal of improved Positron emission tomography equipment.As shown in Fig. 2, positive electron is sent out Penetrating imaging device may include that integrated scintillation crystal 210, multiple photosensor arrays 220, reading circuit (are shown in figure For " signal reading ", corresponding is a reading circuit) 230 and 240 4 parts of data processing module.Scintillation crystal 210 For the scintillation crystal of integration, it is configured to illustrative hexa-prism structure.The center of scintillation crystal 210 has a through-hole, Its two bottom surface for running through scintillation crystal 210.All sides of scintillation crystal 210 and bottom surface can polish, then further according to Need to be coated with highly-reflective coating or highly reflecting films.Photosensor arrays 220 are coupled with scintillation crystal 210.In Fig. 2, flashing There are six photosensor arrays 220 for coupling on the bottom surface of crystal 210.In addition, distinguishing coupling on six sides of scintillation crystal 210 Unify a photosensor arrays 220, i.e. the side of scintillation crystal 210 is coupled with six photosensor arrays altogether.It should be understood that The configuration mode of photosensor arrays shown in Fig. 2 is only exemplary rather than limitation.

As described above, after photosensor arrays 220 acquire optical photon, electric signal is exported.Illustratively, Mei Geguang Electric transducer array 220 can connect a reading circuit 230.It should be understood that the number of reading circuit 230 and its with The connection type of photosensor arrays 220 can be set as needed, and the present invention limits not to this.For example, positive electron All photosensor arrays 220 of transmitting imaging device can be connect with the same reading circuit 230, the reading circuit 230 Received electric signal can be differentiated from which photosensor arrays 220.

Reading circuit 230 can read photosensor arrays 220 output electric signal, and to electric signal measure with Obtain energy information and temporal information.Energy information can indicate the energy for the optical photon that photosensor arrays 220 receive Amount.The energy accumulation for all optical photons that same gamma photons generate is the energy of the gamma photons together.Time letter Breath can indicate the generation time for the optical photon that photosensor arrays 220 receive, and can be considered as generating optical photon Gamma photons reach detector arrival time.Comprehensively consider the time letter for all optical photons that same gamma photons generate Breath, can finally determine the arrival time of gamma photons.Data processing module 240 receives the energy that each reading circuit 230 exports Information and temporal information, and the operation such as data processing, image reconstruction is carried out to these information, to obtain the scanning of object to be imaged Image.

Reading circuit 230 and data processing module 240 can be realized using any suitable hardware, software and/or firmware. Illustratively, data processing module 240 can use field programmable gate array (FPGA), digital signal processor (DSP), answer Miscellaneous programmable logic device (CPLD), micro-control unit (MCU) or central processing unit (CPU) etc. are realized.

The scintillation crystal for the Positron emission tomography equipment that method described herein is applicable in can be, but not limited to be more ribs Column crystal or cylinder crystal.In the following description, the present invention is described mainly in combination with hexa-prism crystal and circular through hole The method of offer, it will be appreciated that being not limited to the sudden strain of a muscle with this particular configuration according to the method for the embodiment of the present invention Bright crystal, the scintillation crystal (such as hexa-prism crystal adds tetragonal through hole) with other constructions are also suitable the present invention.

The positioning accuracy of the response location of gamma photons can directly affect the spatial discrimination of Positron emission tomography equipment Rate and image quality.Energy distribution information (alternatively referred to as " the photon for the optical photon that photoelectric sensor detects can be passed through Distribution "), determine response location of the gamma photons in scintillation crystal.1 to 14 introduction is implemented according to the present invention with reference to the accompanying drawing The response location localization method of example.

Fig. 3 shows the schematic block diagram of response location localization method 300 according to an embodiment of the invention.

As shown in figure 3, response location localization method 300 includes the following steps.

In step S310, receive detected by least one photosensor arrays coupled with scintillation crystal, gamma Photon and scintillation crystal react caused by optical photon energy information, wherein scintillation crystal is integrated flashing Crystal, and scintillation crystal has through-hole, and for accommodating object to be imaged, gamma photons pass through in subject to be imaged through-hole The positron annihilation effect of generation generates.

Gamma photons are hereinbefore described and scintillation crystal reacts and generates optical photon, photosensor arrays inspection It surveys optical photon and exports electric signal, reading circuit measures electric signal and exports the process of the energy information of optical photon, no longer It repeats.In step s310, the energy information of optical photon can be received from reading circuit.

In step S320, according to received energy information determine response location of the gamma photons in scintillation crystal, In, response location is indicated using cylindrical coordinate system.

As described above, scintillation crystal is integral structure, can be polygon prism shape crystal or cylinder crystal, through-hole It can be circular through hole or polygonal through hole.Therefore, for such scintillation crystal, position can be carried out using cylindrical coordinate system Set mark.The position of traditional rectangular scintillation crystal calculates the x it needs to be determined that response location, the coordinate in tri- directions y, z, due to The scintillation crystal that the present invention is directed to is prism-shaped or cylindrical-shaped structure, thus may determine that (θ, r, z) coordinate of response location.? Herein, the θ coordinate in cylindrical coordinate system is known as angular coordinate, and r coordinate is known as radial coordinate, and z coordinate is known as axial coordinate.Cylinder Coordinate system defines the space coordinate of object using plane polar coordinates and Z-direction distance, it can be considered that cylindrical coordinate system packet Include polar coordinate system (θ, r).

Fig. 4 shows the schematic diagram of the cylindrical coordinate system according to an embodiment of the invention applied to hexa-prism crystal. In fig. 4 it is shown that the response location (or reflecting point) of illustrative gamma photons.In cylindrical coordinate system shown in Fig. 4, column Areal coordinate system is using the central point of the bottom surface of scintillation crystal as where the polar coordinate system in origin (origin O) and cylindrical coordinate system Plane is parallel with the bottom surface of scintillation crystal.In fact, the polar coordinates in cylindrical coordinate system shown in Fig. 4, in cylindrical coordinate system During the bottom surface of system and scintillation crystal is generally aligned in the same plane.As shown in figure 4, angular coordinate θ indicates response location in scintillation crystal The line and polar axis Ox of projected position in bottom surface (i.e. in polar coordinate system) and origin O angle (range is 0 °~ 360 °), radial coordinate r indicates projected position of the response location in the bottom surface of scintillation crystal (i.e. in polar coordinate system) to original The linear distance of point；Axial coordinate z indicates response location to the vertical range of bottom surface (arriving polar coordinate system).

According to embodiments of the present invention, a kind of Positron emission tomography for using integrated scintillation crystal can be provided Equipment, the mode that is positioned for the response location to gamma photons, so that solving traditional position calculating method can not The problem of suitable for above-mentioned improved Positron emission tomography equipment.

According to embodiments of the present invention, step S320 may include: effective at least one determining photosensor arrays Photosensor arrays；And response location is determined according at least to energy information corresponding with effective photosensor arrays.

At least one photosensor arrays is covered on the different location of scintillation crystal, sits at the angle for calculating response location When mark, radial coordinate and axial coordinate, it may be necessary to based on detected by the photosensor arrays covered at different location The energy information of optical photon calculates the coordinates of different directions.Therefore, the selector from least one photosensor arrays Point photosensor arrays, and determined using selected photosensor arrays response location can be improved location efficiency and Exactness is fixed really to each coordinate of response location.

In one example, scintillation crystal is polygon prism shape crystal and through-hole is circular through hole or polygonal through hole, institute It states and determines that effective photosensor arrays at least one photosensor arrays include: from least one photoelectric sensor battle array It is selected in column with the photosensor arrays of the side surface coupling of scintillation crystal as effective photosensor arrays.Referring back to figure 1a, scintillation crystal is hexa-prism crystal, and through-hole is circular through hole.Distinguish six sides of scintillation crystal shown in Fig. 1 a It is coupled with six photosensor arrays, it can be using this six photosensor arrays as effective photosensor arrays.This Outside, referring back to Fig. 1 c, scintillation crystal is hexa-prism crystal, and through-hole is polygonal through hole.In such a case, it is possible to Using with six photosensor arrays of the six of scintillation crystal side surface couplings as effective photosensor arrays.Certainly, right For the scintillation crystal shown in Fig. 1 c, six photoelectric sensings can also will being coupled with the six of the through-hole of scintillation crystal inner walls Device array is as effective photosensor arrays.The angular coordinate that effective photosensor arrays are used to calculate response location, When radial coordinate and/or axial coordinate, relative to six photoelectric sensings that will be coupled with the six of the through-hole of scintillation crystal inner walls Device array as effective photosensor arrays for, by six photoelectric sensors with the six of scintillation crystal side surface couplings Array as effective photosensor arrays in the case where calculate obtain response location coordinate it is more acurrate.

In another example, scintillation crystal is cylinder crystal and through-hole is circular through hole or polygonal through hole, described Determine that effective photosensor arrays at least one photosensor arrays include: from least one photosensor arrays The middle photosensor arrays for selecting to couple with the inner wall of through-hole are as effective photosensor arrays.Referring back to Fig. 1 b, dodge Bright crystal is cylinder crystal, and through-hole is hexagon through-hole.Six inner walls of the through-hole of scintillation crystal shown in Fig. 1 b point It is not coupled with six photosensor arrays, it can be using this six photosensor arrays as effective photosensor arrays. In addition, scintillation crystal is cylinder crystal, and through-hole is circular through hole referring back to Fig. 1 d.Scintillation crystal shown in Fig. 1 d The inner wall of through-hole coupled with multiple photosensor arrays, all photoelectric sensor battle arrays that the inner wall with through-hole can be coupled Column are used as effective photosensor arrays.

According to embodiments of the present invention, reaction position is determined according at least to energy information corresponding with effective photosensor arrays The step of setting can include determining that the step of angular coordinate of response location.Fig. 5 shows determination according to an embodiment of the invention The flow diagram of the step S500 of the angular coordinate of response location.

As shown in figure 5, step S500 includes the following steps.

In step S510, effective photosensor arrays are unfolded along the direction parallel with the bottom surface of scintillation crystal, so that Effective photosensor arrays arranged in parallel is obtained in same expansion plane.

Fig. 6 is shown will be shown in Fig. 4, with the photosensor arrays expansion acquisition later of the side surface coupling of scintillation crystal The arrangement schematic diagram of photosensor arrays.According to Fig. 4 and Fig. 6 it is found that by six sides of hexa-prism crystal along flashing The side of the bottom surface of crystal is unfolded, i.e., is unfolded along the direction parallel with the bottom surface of scintillation crystal, can obtain a rectangle.Together When, it is unfolded with six photosensor arrays that six sides couple respectively with the expansion of side, so that six photoelectric transfers Sensor array parallel arranges in the same plane, as shown in Figure 6.

In step S520, plane coordinate system, plane are established for the photosensor arrays after expansion in expansion plane The horizontal axis of coordinate system and the bottom surface of scintillation crystal are parallel, and the longitudinal axis of plane coordinate system and the central axis of scintillation crystal are parallel.

As shown in fig. 6, establishing a plane coordinate system, horizontal axis is x-axis, and the longitudinal axis is z-axis.It should be noted that x-axis shown in fig. 6 It is different from polar axis Ox axis shown in Fig. 4.Plane coordinate system shown in fig. 6 is original with a certain vertex on the bottom surface of scintillation crystal Point, using the expansion direction of effective photosensor arrays as X direction, and cylindrical coordinate system shown in Fig. 4 is with scintillation crystal The central point of bottom surface is origin and using the line direction on a certain vertex on the central point and bottom surface as polar axis direction, because Both this is different.In Fig. 6, longitudinal axis z-axis is parallel with the z-axis direction of cylindrical coordinate system.Photoelectric sensing after expansion shown in fig. 6 Device array is made of six photosensor arrays, the corresponding abscissa length of each photoelectric sensing array and scintillation crystal Bottom surface side length it is identical, be shown as L in Fig. 4 and Fig. 6.

In step S530, determined and the photoelectric transfer after expansion according to energy information corresponding with effective photosensor arrays The corresponding specific distribution of photons of sensor array.

When a gamma photons and scintillation crystal react, effective photosensor arrays are (referring to six in Fig. 6 Photosensor arrays) at least partly photosensor arrays can detecte optical photon.According to the anti-of gamma photons Answer the difference of position, the case where optical photon that effective photosensor arrays detect is also different.It can be according to after deployment Photosensor arrays in the arrangement modes of each effectively photosensor arrays will each effective photosensor arrays detection To the energy information of optical photon be integrated into the photosensor arrays corresponding specific photon point obtained together with after expansion Cloth.The specific distribution of photons is the Energy distribution situation on the photosensor arrays of optical photon after deployment, in general Approximate Gaussian distribution is presented in it.

It should be understood that since scintillation crystal is integrated scintillation crystal, the photoelectricity coupled with the inner wall of its side or through-hole Sensor array is arranged around object to be imaged, therefore when a positron annihilation events occur, two generated Contrary gamma photons react on two opposite positions of scintillation crystal, and effective photosensor arrays receive To be optical photon caused by the two gamma photons optical signal.Two gamma photons have respective Energy distribution area Therefore domain when carrying out response location positioning, the Energy distribution region of two gamma photons can be distinguished, be analyzed respectively The response location of each gamma photons.Approach described herein is determined primarily directed to the response location of individual gamma photons Position.

In step S540, the horizontal seat of projected position of the response location in plane coordinate system is calculated according to specific distribution of photons Mark.

Photosensor arrays after expansion are distributed in same expansion plane, and corresponding specific distribution of photons is presented Approximate Gaussian distribution is similarly to the photoelectric sensor coupled in traditional Positron emission tomography with rectangular scintillation crystal The case where optical photon that (or photosensor arrays) detect.Therefore, can based on the photoelectric sensor battle array after expansion Corresponding specific distribution of photons is arranged, the cross of projected position of the response location in plane coordinate system is calculated using traditional location algorithm Coordinate x.Of course, it is possible to the ordinate z of projected position of the response location in plane coordinate system be calculated simultaneously, under can be used for The calculating of the axial coordinate of response location described in text.Traditional location algorithm employed in step S540 can be center of gravity A variety of different methods such as method, artificial neural network, analytic method, the present invention limit not to this.

In step S550, sat according to the angle that the abscissa of the geometrical characteristic of scintillation crystal and projected position calculates response location Mark.

It is calculated after the x coordinate of response location, the geometrical characteristic that can use scintillation crystal calculates and the x coordinate pair The angular coordinate answered.Illustratively, the bottom surface of scintillation crystal can be considered as circle, and the pass based on radian and arc length as a result, It is (i.e. geometrical characteristic), angular coordinate θ can be calculated by the following formula:

Wherein, θ indicates the angular coordinate of response location, and fix indicates to be rounded downwards, and mod expression takes the remainder, and x indicates reaction position The abscissa of the projected position in plane coordinate system is set, a indicates the matrix number of effective photosensor arrays (with Fig. 4's For hexa-prism crystal, the face number of side of matrix number and scintillation crystal of effective photosensor arrays is equal, i.e., a is The face number of the side of scintillation crystal), L indicates each photosensor arrays in effective photoelectric sensor and scintillation crystal The parallel side in bottom surface side length (each photoelectric sensor by taking the hexa-prism crystal of Fig. 4 as an example, in effective photoelectric sensor The side length on array, parallel with the bottom surface of scintillation crystal side is equal with the side length of the bottom surface of scintillation crystal, i.e. L is scintillation crystal Bottom surface side length), R indicates scintillation crystal, the face that couples with effective photosensor arrays on the bottom surface of scintillation crystal (by taking the hexa-prism crystal of Fig. 4 as an example, R is the outer of the bottom surface of scintillation crystal to the radius of the circumscribed circle for the polygon that projection obtains Connect round radius).As described above, cylindrical coordinate system is using the central point of the bottom surface of scintillation crystal as origin and cylindrical coordinates Plane where polar coordinate system in system is parallel with the bottom surface of scintillation crystal.In this case, formula (1) is applicable.

No matter effective photosensor arrays be with the photosensor arrays of the side surface coupling of scintillation crystal or with sudden strain of a muscle The photosensor arrays of the inner wall coupling of the through-hole of bright crystal, formula (1) can be set up.

In scintillation crystal be polygon prism shape crystal and effective photosensor arrays are the side surface couplings with scintillation crystal Photosensor arrays when, the shape of polygon prism more levels off to cylinder, i.e. when the prism number of polygon prism is more, angle calculated Coordinate θ is more accurate.It similarly, is polygonal through hole in the through-hole of scintillation crystal and effective photosensor arrays are and flashing is brilliant When the photosensor arrays of the inner wall coupling of the through-hole of body, the shape of polygon closer to circle, i.e. get over by the number of edges of polygon When more, angular coordinate θ calculated is more accurate.

In above-mentioned steps S500, pass through at least partly photosensor arrays that will be coupled with integrated scintillation crystal Expansion obtains the distribution of photons situation similar with traditional Positron emission tomography, and then calculates the angular coordinate of response location, this Mode is simple and practical, convenience of calculation.

According to embodiments of the present invention, it is determined according at least to energy information corresponding with effective photosensor arrays described anti- It answers position can further include: special angle corresponding with the angular coordinate of response location is determined according to angular error distribution function Error, wherein angular error distribution function is known；And it is carried out according to angular coordinate of the special angle error to response location Calibration.

It is appreciated that when intersection of the response location of gamma photons close to two sides of scintillation crystal, it is seen that light The energy of son may be randomly deposited on the photosensor arrays with the two side surface couplings.It is different from other positions It is that, since two sides have certain angle, the Energy distribution deposited on photosensor arrays also can be with two sides Angle variation and change.Pass through the accuracy of the angular coordinate of the response location of gamma photons determined by above-mentioned localization method Also can be declined.Therefore, the angular coordinate of the response location of gamma photons can be calibrated by certain calibration means.

Angular error distribution function can indicate under different angular coordinates, calculated according to the energy information of optical photon The error of angular coordinate is how many.Angular error distribution function can be obtained by test in advance.It is described below a kind of illustrative Test method, for obtaining angular error distribution function.

According to embodiments of the present invention, spy corresponding with the angular coordinate of response location is being determined according to angular error distribution function Before determining angular error, method 300 be may further include: for each of control angle set to irradiation angle, radiated Angle between source and the polar axis for meeting the polar coordinate system in the line and cylindrical coordinate system of detecting module is equal to this to irradiation angle In the case where, it receives and meets detected by photoelectric sensor in detecting module, the first test gamma photons and meets detection Test crystal in module react caused by optical photon first time information, and receive at least one photoelectric sensing Detected by device array, second test gamma photons and scintillation crystal react caused by optical photon the first energy Information and the second temporal information；For each of control angle set to irradiation angle, when according to first time information and second Between information carry out meeting event detection；For each of control angle set to irradiation angle, refer to when meeting event detection outcome When showing that meeting event occurs, according in the first energy information corresponding with event is met and effective photosensor arrays Corresponding portion of energy information calculates the angular coordinate of test reaction position of the second test gamma photons in scintillation crystal；For Each of angle set is compareed to irradiation angle, calculate test reaction position angular coordinate and this to the difference between irradiation angle, To obtain test angle error；According to the test angle error obtained respectively for all pairs of irradiation angles in control angle set Determine angular error distribution function, wherein radioactive source and meet detecting module synchronously around the center axis rotation of scintillation crystal, To obtain each pair of irradiation angle in control angle set.

Fig. 7 shows according to an embodiment of the invention for obtaining the structure of the test device of angular error distribution function Schematic diagram.As shown in fig. 7, can be arranged one in the side of scintillation crystal meets detecting module.Meet detecting module by one (it is identical as the ingredient of scintillation crystal as described herein, and the size of the test crystal is according to space point for the test crystal of small size Resolution it needs to be determined that, the section of crystal is smaller, and it is higher to meet detection accuracy) and photoelectric sensor composition, for controlling gamma Incoming position range of the ray on scintillation crystal surface.The photoelectric sensor met in detecting module is connect with a reading circuit, Reading circuit of the reading circuit independently of Positron emission tomography equipment.Data acquisition board shown in Fig. 7 can integrate above The data processing module.Illustratively, reaction described herein location position device may include for determining that angle is missed The module of poor distribution function.During determining angular error distribution function, need to carry out signal and meet detection (to meet thing Part detection), as described below.

Meet detecting module together with radioactive source at the uniform velocity around the axis for the Positron emission tomography equipment for including scintillation crystal Rotation, i.e., around the center axis rotation of scintillation crystal.During rotation, radioactive source keeps opposing stationary shape with detecting module is met State, and radioactive source with meet detecting module line and scintillation crystal bottom surface (i.e. with the polar coordinate system in cylindrical coordinate system Place plane) it is parallel.During rotation, a series of pairs of irradiation angles can be obtained, it is final to obtain control angle set.Compare angle Irradiation angle can be distributed in the range of 0 °~360 ° in degree set.In radioactive source and meet the line and column of detecting module In the case that the Mr. Yus such as the angle between the polar axis of the polar coordinate system in areal coordinate system compare angle, θ ' (such as 10 °), record should Angle (namely control angle, θ '), and carry out meeting event detection.When the event that is determined for compliance with occurs, record compares angle herein Spend the Energy distribution of optical photon detected by lower photosensor arrays.

Meet the time that event detection will mainly meet optical photon detected by photoelectric sensor in detecting module The temporal information of optical photon detected by information and at least one photosensor arrays coupled with scintillation crystal carries out Comparison, judge whether two gamma photons are generated by same positron annihilation events, i.e., whether the two is to meet gamma photons It is right.Meeting event detection can be realized using conventional method, be repeated herein not to this.

When the event that is determined for compliance with occurs, can use and effective photosensor arrays (such as side with scintillation crystal The photosensor arrays of face coupling) corresponding part energy information calculates gamma photons that scintillation crystal receives (the Two test gamma photons) test reaction position angular coordinate θ ".The process of calculating angular coordinate θ " and step S500 shown in fig. 5 It is identical, it repeats no more.

Then, by comparing the available control angle, θ of the size of θ ' and θ " ' under test angle error delta θ.For right All pairs of irradiation angles in irradiation angle set are performed both by similar operations, can obtain angular error distribution function.Then, angle is utilized Error distribution function is spent, practical calculated angular coordinate θ can be calibrated.

According to embodiments of the present invention, reaction position is determined according at least to energy information corresponding with effective photosensor arrays The step of setting can include determining that the step of radial coordinate of response location.Fig. 8 shows according to an embodiment of the invention true Determine the flow diagram of the step S800 of the radial coordinate of response location.

As shown in figure 8, step S800 includes the following steps.

In step S810, effective photosensor arrays are unfolded along the direction parallel with the bottom surface of scintillation crystal, so that Effective photosensor arrays arranged in parallel is obtained in same expansion plane.

In step S820, plane coordinate system, plane are established for the photosensor arrays after expansion in expansion plane The horizontal axis of coordinate system and the bottom surface of scintillation crystal are parallel, and the longitudinal axis of plane coordinate system and the central axis of scintillation crystal are parallel.

In step S830, determined and the photoelectric transfer after expansion according to energy information corresponding with effective photosensor arrays The corresponding specific distribution of photons of sensor array.

In step S840, standard deviation and peak value based on specific distribution of photons are to being used to indicate standard deviation, peak value and radius The look-up table of corresponding relationship between coordinate is inquired, to obtain the radial coordinate of response location.

Step S810 to S830 shown in Fig. 8 is identical as step S510 to S530 shown in fig. 5, and those skilled in the art are logical The description that reading is crossed above with respect to each step shown in fig. 5 is understood that step S810 to S830, and details are not described herein.

Reaction depth effect (depth of interaction in traditional Positron emission tomography system Effects, DOI) it is an important factor for influencing the spatial resolution of Positron emission tomography system, it is embodied in gamma photons The inaccuracy of third dimension in three dimension location (crystal decoding can only realize two-dimensional localization).Half proposed in the present invention Diameter coordinate r indicates projected position of the response location in the bottom surface of scintillation crystal (i.e. in polar coordinate system) to the straight line of origin Distance, from the perspective of the photosensor arrays coupled with the inner wall of the side of scintillation crystal or through-hole, radial coordinate r It can reflect reaction depth information of the gamma photons in scintillation crystal.It is therefore possible to use step S800 shown in Fig. 8 is calculated Radial coordinate r, to obtain reaction depth information.

Specifically, photon is established by square shaped scintillation crystal (LYSO crystal, size are 60mm x 60mm x 20mm) Propagation Simulation finds that the distribution of photons for the optical photon being collected on photosensor arrays is approximately Gaussian Profile, and with The change of reaction depth, the standard deviation sigma and energy peak E of Gaussian Profile also change correspondingly.Basis is shown respectively in Fig. 9 a to 9i One embodiment of the invention under different gammaphoton reaction depths to the light of optical photon detected by photoelectric sensor Son distribution carries out the fitting result of Gauss curve fitting.Fig. 9 a to 9i successively indicate the reaction depth of gamma photons be respectively 2mm, 4mm, 6mm ... in the case where 18mm, the distribution of photons of optical photon detected by photoelectric sensor in the x-direction is carried out high Result after this fitting.By Fig. 9 a to 9i as it can be seen that reaction depth of the gamma photons in scintillation crystal is different, on photoelectric sensor The distribution of photons for the optical photon being collected into is also different.Although Fig. 9 a to 9i shows the photon point for rectangular scintillation crystal Cloth situation, but for integrated scintillation crystal of the present invention, the inner wall coupling of the side You Yuqi or through-hole After the expansion, the distribution of photons of optical photon corresponding with the photosensor arrays after expansion has photosensor arrays The regularity of distribution similar with Fig. 9 a to 9i.

Therefore, look-up table relevant to different radii coordinate can be established in advance, according to the photon of actually detected acquisition point Cloth situation (i.e. above-mentioned specific distribution of photons) and the look-up table determine the radial coordinate r of the response location of gamma photons.

One that the method for obtaining the look-up table of the corresponding relationship between standard deviation, peak value and radial coordinate is described below shows Example.

According to embodiments of the present invention, based on specific distribution of photons standard deviation and peak value to being used to indicate standard deviation, peak Before the look-up table of corresponding relationship between value and radial coordinate is inquired, method 300 be may further include: for control Each control radius in radius set, between the origin in radioactive source and the line and cylindrical coordinate system that meet detecting module Distance be equal to the control radius in the case where, receive meet detected by photoelectric sensor in detecting module, third is surveyed Examination gamma photons react the third temporal information of generated optical photon with the test crystal for meeting in detecting module, and Receive detected by least one photosensor arrays, the 4th test gamma photons and scintillation crystal react it is produced Optical photon the second energy information and the 4th temporal information；For control radius set in each control radius, according to Third temporal information and the 4th temporal information carry out meeting event detection；For control radius set in each control radius, When meet event detection outcome instruction meet event occur when, according to it is in the second energy information corresponding with event is met, And the corresponding portion of energy information of effective photosensor arrays determine be unfolded after the corresponding test of photosensor arrays Distribution of photons；According to for the test distribution of photons that obtains respectively of all control radiuses in control radius set standard deviation, Peak value and corresponding control radius determine look-up table, wherein radioactive source and meet detecting module synchronously along scintillation crystal It moves radially, to obtain each control radius in control radius set.

Although the response location of gamma photons can be set in simulation process, obtain what photosensor arrays were collected into Distribution of photons information, to establish look-up table.But the response location due to can not directly set gamma photons in real process, It is therefore proposed that it is a kind of establish in the actual process include distribution of photons information and gamma photons response location radial coordinate it Between corresponding relationship look-up table test method.

The test method is similar with the process of above-mentioned acquisition angular error distribution function, using meeting detecting module come real It is existing.

Figure 10 shows according to an embodiment of the invention for obtaining the structural schematic diagram of the test device of look-up table.Such as Shown in Figure 10, it is arranged one in the top of scintillation crystal and meets detecting module.The structure shown in Fig. 10 for meeting detecting module and Working method is similar with the detecting module shown in Fig. 7 that meets, and repeats no more.

In the embodiment shown in fig. 10, meet detecting module moving radially along scintillation crystal together with radioactive source. In moving process, radioactive source and meet detecting module and keep relative static conditions, and radioactive source and meets detecting module The central axis (i.e. with the z-axis in cylindrical coordinate system) of line and scintillation crystal is parallel.In moving process, it can obtain a series of Radius is compareed, it is final to obtain control radius set.Original in radioactive source and the line and cylindrical coordinate system that meet detecting module The distance between point (in other words with the z-axis of cylindrical coordinate system) is equal in the case that this compares radius, carries out meeting event detection. When the event that is determined for compliance with occurs, record compares the energy of optical photon detected by photosensor arrays under radius herein Information.Furthermore, it is possible to using with effective photosensor arrays (such as photoelectric sensor battle array with the side surface coupling of scintillation crystal Column) corresponding part energy information come determine be unfolded after the corresponding test distribution of photons of photosensor arrays.

Aforesaid operations are performed both by for all control radiuses in control radius set, can obtain and compare radius pair with each The test distribution of photons answered.Then, it can establish between the standard deviation sigma for being used to indicate distribution of photons, peak E and radial coordinate r Corresponding relationship look-up table.During reality positions response location, it can be determined and be reacted based on the look-up table The radial coordinate of position.

It will be understood by those skilled in the art that for the photoelectric sensing coupled with the inner wall of the side of scintillation crystal or through-hole For device array (i.e. effective photosensor arrays), the reaction depth of gamma photons can be logical by subtracting radial coordinate r Pore radius calculates.

Illustratively, the radial coordinate of the response location of gamma photons can also be determined using other modes.Figure 11 is shown The flow diagram of the step S1100 of the radial coordinate of determining response location in accordance with another embodiment of the present invention.

As shown in figure 11, step S1100 includes the following steps.

In step S1110, the specific bottom surface of scintillation crystal is divided into predetermined number destination region.

In step S1120, corresponding region is selected from predetermined number destination region according to the angular coordinate of response location.

In step S1130, the photosensor arrays coupled with corresponding region are selected.

In step S1140, according to the energy of optical photon detected by the photosensor arrays coupled with corresponding region Measure the radial coordinate that information calculates response location.

According to the present embodiment, in the case where determining the angular coordinate of response location of gamma photons, can be sat based on the angle Mark the radial coordinate for determining response location.

Figure 12 shows the photoelectric sensor battle array that basis according to an embodiment of the invention is coupled with the bottom surface of scintillation crystal The energy information of the detected optical photon of column determines the schematic diagram of radial coordinate.

As shown in figure 12, the specific bottom surface of scintillation crystal can be divided by several regions according to angle first.It is specific Bottom surface can be upper bottom surface or bottom surface, and depending on can according to need, the present invention limits not to this.With hexa-prism crystalline substance For body, crystal can be divided into Ith area to VIth area (as shown in figure 12).

Then, corresponding region can be selected from ready-portioned region according to the angular coordinate θ of response location, and selects and is somebody's turn to do The photosensor arrays of corresponding region coupling.For example, can choose the photoelectric sensing coupled with Ith area when working as θ ∈ [0,60 °] Device array can choose the photosensor arrays, etc. coupled with IIth area, be not listed one by one when θ ∈ [60 °, 120 °]. In Figure 12, the response location of gamma photons is corresponding with IIth area, the photosensor arrays of IIth area the Shi Yu coupling of selection.

The distribution of photons situation of optical photon detected by selected photosensor arrays and traditional positive electron The distribution of photons situation for emitting optical photon detected by the photoelectric sensor in imaging system is similar, therefore can be using biography The response location for location algorithm (such as gravity model appoach, artificial neural network, analytic method etc.) calculating gamma photons of uniting is selected The coordinate (x, y) of projected position in the respective coordinates system (being shown as x/y plane coordinate system in Figure 12) of photosensor arrays. Finally radial coordinate r can be solved by the geometrical relationship between the coordinate y of projected position and the radial coordinate r of response location.

According to embodiments of the present invention, reaction position is determined according at least to energy information corresponding with effective photosensor arrays The step of setting can include determining that the step of axial coordinate of response location.Figure 13 shows according to an embodiment of the invention true Determine the flow diagram of the step S1300 of the axial coordinate of response location.

As shown in figure 13, step S1300 includes the following steps.

In step S1310, effective photosensor arrays are unfolded along the direction parallel with the bottom surface of scintillation crystal, so that Effective photosensor arrays arranged in parallel is obtained in same expansion plane.

In step S1320, plane coordinate system, plane are established for the photosensor arrays after expansion in expansion plane The horizontal axis of coordinate system and the bottom surface of scintillation crystal are parallel, and the longitudinal axis of plane coordinate system and the central axis of scintillation crystal are parallel.

In step S1330, determined and the photoelectricity after expansion according to energy information corresponding with effective photosensor arrays The corresponding specific distribution of photons of sensor array.

In step S1340, the vertical of projected position of the response location in plane coordinate system is calculated according to specific distribution of photons Coordinate.

In step S1350, determine that the ordinate of projected position is the axial coordinate of response location.

Step S1310 to S1330 shown in Figure 13 and step S810 to S830 shown in Fig. 8 and step shown in fig. 5 S510 to S530 is identical, and description of the those skilled in the art by reading above with respect to Fig. 5 and each step shown in Fig. 8 can be managed Step S1310 to S1330 is solved, details are not described herein.

Above with respect in the description of step S540, the throwing in response location plane coordinate system shown in Fig. 6 is described The calculation of the abscissa x of shadow position.It similarly, can be based on specific light corresponding with the photosensor arrays after expansion Son distribution, the ordinate z of projected position of the response location in plane coordinate system is calculated using traditional location algorithm.In fact, When projected position in the xz plane coordinate system shown in Fig. 6 to response location positions, it can calculate and obtain its coordinate (x, z), it can while obtaining abscissa x and ordinate z.Ordinate z is that the axial of the response location of gamma photons is sat Mark z.

According to embodiments of the present invention, reaction position is determined according at least to energy information corresponding with effective photosensor arrays Setting may further include: according to axial error distribution function determination specific axial mistake corresponding with the axial coordinate of response location Difference, wherein axial error distribution function is known；And according to specific axial error to the axial coordinate of the response location It is calibrated.

When upper bottom surface or bottom surface of the response location of gamma photons close to scintillation crystal, due to upper bottom surface or bottom surface To the reflex of optical photon, the positioning accuracy that is positioned using response location of traditional location algorithm to gamma photons compared with Difference has certain edge effect.Therefore the present embodiment proposes the side that the axial coordinate z of a kind of pair of response location is calibrated Method.This method can compensate error brought by the edge effect as caused by the bottom surface of scintillation crystal, improve determining for response location Position precision.

Similarly with the calibrating mode of the angular coordinate of response location, it can use known axial error distribution function to anti- The axial coordinate of position is answered to be calibrated.Axial error distribution function can indicate under axially different coordinate, according to visible light The error of the energy information axial coordinate calculated of son is how many.Axial error distribution function can be obtained by test in advance ?.A kind of illustrative test method is described below, for obtaining axial error distribution function.

According to embodiments of the present invention, method 300 may further include: right for each of control axial coordinate set According to axial coordinate, where the polar coordinate system in radioactive source and the line and cylindrical coordinate system that meet detecting module between plane In the case that distance is equal to the control axial coordinate, receive meet detected by photoelectric sensor in detecting module, the 5th Test gamma photons react the 5th temporal information of generated optical photon with the test crystal for meeting in detecting module, And receive detected by least one photosensor arrays, the 6th test gamma photons and scintillation crystal react and are produced The third energy information and the 6th temporal information of raw optical photon；It is axial for each control in control axial coordinate set Coordinate carries out meeting event detection according to the 5th temporal information and the 6th temporal information；For in control axial coordinate set Each control axial coordinate, when meet event detection outcome instruction meet event occur when, according to corresponding with event is met Portion of energy information in third energy information, corresponding with effective photosensor arrays calculates the 6th test gamma photons and exists The axial coordinate of test reaction position in scintillation crystal；For control axial coordinate set in each control axial coordinate, The axial coordinate for calculating test reaction position compares the difference between axial coordinate with this, to obtain test axial error；According to The test axial error obtained respectively for all control axial coordinates in control axial coordinate set determines axial error point Cloth function, wherein radioactive source and meets detecting module and synchronously moved along the direction parallel with the central axis of scintillation crystal, with Obtain each control axial coordinate in control axial coordinate set.

Figure 14 shows according to an embodiment of the invention for obtaining the knot of the test device of axial error distribution function Structure schematic diagram.As shown in figure 14, it is arranged one in the side of scintillation crystal and meets detecting module.Meet detection mould shown in Figure 14 The structure and working method of block meet detecting module and the detecting module shown in Fig. 10 that meets is similar with shown in Fig. 7, no longer superfluous It states.

In the embodiment shown in fig. 14, meet detecting module together with radioactive source along flat with the central axis of scintillation crystal Capable direction is mobile.In moving process, radioactive source keeps relative static conditions, and radioactive source and symbol with detecting module is met It is parallel to close the bottom surface (i.e. with plane where the polar coordinate system in cylindrical coordinate system) of line and the scintillation crystal of detecting module.It is moving During dynamic, a series of control axial coordinates can be obtained, it is final to obtain control axial coordinate set.In radioactive source and meet spy It surveys the line of module and is equal to the feelings for compareing axial coordinate z ' with the distance between plane where the polar coordinate system in cylindrical coordinate system Under condition, the distance (i.e. control axial coordinate z ') is recorded, and carry out meeting event detection.When the event that is determined for compliance with occurs, note Record compares the energy information of optical photon detected by photosensor arrays under axial coordinate herein.

When the event that is determined for compliance with occurs, can use and effective photosensor arrays (such as side with scintillation crystal The photosensor arrays of face coupling) corresponding part energy information calculates gamma photons that scintillation crystal receives (the Six test gamma photons) test reaction position axial coordinate z ".It is walked shown in the process and Figure 13 of calculating axial coordinate z " Rapid S1300 is identical, repeats no more.

Then, by comparing the test axial error Δ z under the available control axial coordinate z ' of the size of z ' and z ".Needle Similar operations are performed both by all control axial coordinates in control axial coordinate set, axial error distribution letter can be obtained Number.Then, using axial error distribution function, practical calculated axial coordinate z can be calibrated.

Illustratively, the axial coordinate of the response location of gamma photons can also be determined using other modes.Show at one In example, the step of determining the radial coordinate of response location may include: according to the photoelectricity coupled with two bottom surfaces of scintillation crystal The energy information of optical photon detected by sensor array determines the axial coordinate of response location.

The optical photon that same gamma photons generate may be distributed on two bottom surfaces of scintillation crystal, by with two bottom surfaces The photosensor arrays of coupling receive simultaneously.It is appreciated that when the response location of gamma photons is located at the half of scintillation crystal When at height, the energy of optical photon detected by the photosensor arrays that are coupled with the upper bottom surface of scintillation crystal and with sudden strain of a muscle The energy of optical photon detected by the photosensor arrays of the bottom surface coupling of bright crystal is of substantially equal, the two ratio Example is 1:1.When the response location of gamma photons is closer from the upper bottom surface of scintillation crystal, coupled with the upper bottom surface of scintillation crystal The energy of optical photon detected by photosensor arrays will be greater than the photoelectric sensing coupled with the bottom surface of scintillation crystal The energy of optical photon detected by device array.That is, the response location of gamma photons is closer apart from a certain bottom surface, with The energy of optical photon detected by the photosensor arrays of bottom surface coupling is relative to another bottom surface with scintillation crystal The ratio of the energy of optical photon detected by the photosensor arrays of coupling will be bigger.Therefore, reason can be passed through in advance By or experiment obtain detected by the photosensor arrays that are coupled under axially different coordinate with two bottom surfaces of scintillation crystal Optical photon energy proportion, then during reality positions response location, can according to scintillation crystal The coupling of two bottom surfaces photosensor arrays detected by the energy information of optical photon determine the axial direction of response location Coordinate.

It is described below another for determining the mode of the axial coordinate of the response location of gamma photons.For brilliant with flashing For the photosensor arrays of upper bottom surface or the bottom surface coupling of body, the axial coordinate z of response location is equal to reaction depth. In the description of the method for determination of the radial coordinate r above with respect to response location, in conjunction with Fig. 8 and Fig. 9 a to 9i be described in detail as What determines radial coordinate r using the look-up table for the corresponding relationship being used to indicate between standard deviation, peak value and radial coordinate.For For axial coordinate z, also there is similar processing mode.That is, can be according to the upper bottom surface or bottom surface with scintillation crystal The energy information of optical photon detected by the photosensor arrays of coupling and be similarly used for instruction standard deviation, peak value and The look-up table of corresponding relationship between radial coordinate determines axial coordinate z.It is noted that for determining axial coordinate z's It is visible detected by the photosensor arrays that are coupled with the upper bottom surface of scintillation crystal or bottom surface involved in look-up table The standard deviation and peak value of the distribution of photons of photon.In this way determine axial coordinate z the step of can with reference to it is above-mentioned about The description of step S840 and Fig. 9 a to 9i, details are not described herein.

It is appreciated that being coupled since the upper bottom surface of scintillation crystal and the area of bottom surface are smaller with upper bottom surface or bottom surface The optical photon that receives of photosensor arrays it is more unlike the photosensor arrays with side surface coupling, therefore compare Preferably between the energy information of the optical photon according to detected by the photosensor arrays coupled with two bottom surfaces Proportionate relationship determines the axial coordinate of response location.

According to a further aspect of the invention, a kind of response location positioning device is provided.Figure 15 shows a reality according to the present invention Apply the schematic block diagram of the response location positioning device 1500 of example.As shown in figure 15, response location positioning device 1500 includes the One receiving module 1510 and position determination module 1520.

First receiving module 1510 is detected for receiving at least one photosensor arrays coupled with scintillation crystal To, gamma photons and scintillation crystal react caused by optical photon energy information, wherein scintillation crystal is one The scintillation crystal of change, and scintillation crystal has through-hole, and for accommodating object to be imaged, gamma photons pass through to be imaged through-hole The esoteric positron annihilation effect of object generates.

Position determination module 1520 be used for according to received energy information determine that gamma photons are anti-in scintillation crystal Answer position, wherein response location is indicated using cylindrical coordinate system.

According to embodiments of the present invention, position determination module 1520 may include: that sensor determines submodule, for determining extremely Effective photosensor arrays in few photosensor arrays；And position determination submodule, for according at least to Effective corresponding energy information of photosensor arrays determines response location.

According to embodiments of the present invention, scintillation crystal is polygon prism shape crystal and through-hole is that circular through hole or polygon are logical Hole, sensor determine that submodule may include: first selecting unit, for from least one photosensor arrays selection with The photosensor arrays of the side surface coupling of scintillation crystal are as effective photosensor arrays.

According to embodiments of the present invention, scintillation crystal is cylinder crystal and through-hole is circular through hole or polygonal through hole, Sensor determines that submodule may include: the second selecting unit, for the selection from least one photosensor arrays and leads to The photosensor arrays of the inner wall coupling in hole are as effective photosensor arrays.

According to embodiments of the present invention, position determination submodule may include: the first expansion unit, be used for effective photoelectric transfer Sensor array is unfolded along the direction parallel with the bottom surface of scintillation crystal, so that effectively photosensor arrays arranged in parallel is same In one expansion plane；First establishment of coordinate system unit, for being built in expansion plane for the photosensor arrays after expansion Vertical plane coordinate system, the horizontal axis of plane coordinate system and the bottom surface of scintillation crystal are parallel, the longitudinal axis and scintillation crystal of plane coordinate system Central axis it is parallel；First distribution of photons determination unit, for according to energy information corresponding with effective photosensor arrays Determine specific distribution of photons corresponding with the photosensor arrays after expansion；Abscissa computing unit, for according to specific light Son distribution calculates the abscissa of projected position of the response location in plane coordinate system；And angular coordinate computing unit, it is used for root The angular coordinate of response location is calculated according to the geometrical characteristic of scintillation crystal and the abscissa of projected position.

According to embodiments of the present invention, angular coordinate computing unit may include computation module, for being calculated based on following formula The angular coordinate of response location:

Wherein, θ indicates the angular coordinate of response location, and fix indicates to be rounded downwards, and mod expression takes the remainder, and x indicates projection position The abscissa set, a indicate that the matrix number of effective photosensor arrays, L indicate each of effective photosensor arrays The side length on photosensor arrays, parallel with the bottom surface of scintillation crystal sides, R indicate scintillation crystal, with effective photoelectric transfer The face of sensor array coupling projects the radius of the circumscribed circle of the polygon of acquisition on the bottom surface of scintillation crystal,

Wherein, cylindrical coordinate system is sat by the pole in origin and cylindrical coordinate system of the central point of the bottom surface of scintillation crystal Plane where mark system is parallel with the bottom surface of scintillation crystal.

According to embodiments of the present invention, position determination submodule may further include: angular error determination unit, be used for root Special angle error corresponding with the angular coordinate of response location is determined according to angular error distribution function, wherein angular error distribution Function is known；And angular coordinate calibration unit, for carrying out school according to angular coordinate of the special angle error to response location It is quasi-.

According to embodiments of the present invention, device 1500 may further include: the second receiving module, for for irradiation angle Each of set is to irradiation angle, the pole of the polar coordinate system in radioactive source and the line and cylindrical coordinate system that meet detecting module In the case that angle between axis is equal to this to irradiation angle, receive detected by the photoelectric sensor met in detecting module, First test gamma photons react first time of generated optical photon with the test crystal for meeting in detecting module Information, and receive detected by least one photosensor arrays, the second test gamma photons and scintillation crystal occur instead The first energy information and the second temporal information of optical photon caused by answering；First meets event checking module, for being directed to Each of angle set is compareed to irradiation angle, carries out meeting event detection according to first time information and the second temporal information； Test angles coordinate calculation module, for, to irradiation angle, referring to when meeting event detection outcome for each of control angle set When showing that meeting event occurs, according in the first energy information corresponding with event is met and effective photosensor arrays Corresponding portion of energy information calculates the angular coordinate of test reaction position of the second test gamma photons in scintillation crystal；Test Angular error computing module, for, to irradiation angle, the angle for calculating test reaction position to be sat for each of control angle set It marks and this is to the difference between irradiation angle, to obtain test angle error；Angular error distribution function determining module is used for basis The test angle error obtained respectively for all pairs of irradiation angles in control angle set determines angular error distribution function, In, radioactive source and meet detecting module synchronously around the center axis rotation of scintillation crystal, to obtain in control angle set Each pair of irradiation angle.

According to embodiments of the present invention, position determination submodule may include: the second expansion unit, be used for effective photoelectric transfer Sensor array is unfolded along the direction parallel with the bottom surface of scintillation crystal, so that effectively photosensor arrays arranged in parallel is same In one expansion plane；Second establishment of coordinate system unit, for being built in expansion plane for the photosensor arrays after expansion Vertical plane coordinate system, the horizontal axis of plane coordinate system and the bottom surface of scintillation crystal are parallel, the longitudinal axis and scintillation crystal of plane coordinate system Central axis it is parallel；Second distribution of photons determination unit, for according to energy information corresponding with effective photosensor arrays Determine specific distribution of photons corresponding with the photosensor arrays after expansion；And query unit, for being based on specific photon The standard deviation and peak value of distribution look into the look-up table for the corresponding relationship being used to indicate between standard deviation, peak value and radial coordinate It askes, to obtain the radial coordinate of response location.

According to embodiments of the present invention, device 1500 may further include: third receiving module, for for control radius Each control radius in set, between the origin in radioactive source and the line and cylindrical coordinate system that meet detecting module away from In the case where being equal to the control radius, receive detected by the photoelectric sensor met in detecting module, third test gal Ma photon reacts the third temporal information of generated optical photon with the test crystal for meeting in detecting module, and receives Detected by least one photosensor arrays, the 4th test gamma photons and scintillation crystal react it is generated can The second energy information and the 4th temporal information of light-exposed son；Second meets event checking module, for for control radius set In each control radius, carry out meeting event detection according to third temporal information and the 4th temporal information；Test distribution of photons Determining module, for meeting event when meeting event detection outcome instruction for each control radius in control radius set When generation, according to and meet in corresponding second energy information of event, corresponding with effective photosensor arrays part Energy information determines test distribution of photons corresponding with the photosensor arrays after expansion；Look-up table determining module is used for root According to for the standard deviation of the test distribution of photons that obtains respectively of all control radiuses in control radius set, peak value and corresponding Control radius determines look-up table, wherein radioactive source and meets detecting module synchronously moving radially along scintillation crystal, to obtain Each control radius in radius set must be compareed.

According to embodiments of the present invention, position determination submodule may include: area division unit, for by scintillation crystal Specific bottom surface is divided into predetermined number destination region；Area selecting unit, for according to the angular coordinate of response location from predetermined number Region in select corresponding region；Sensor selection unit, the photosensor arrays for selecting to couple with corresponding region；With And radial coordinate computing unit, for the optical photon according to detected by the photosensor arrays coupled with corresponding region The radial coordinate of energy information calculating response location.

According to embodiments of the present invention, position determination submodule may include: third expansion unit, be used for effective photoelectric transfer Sensor array is unfolded along the direction parallel with the bottom surface of scintillation crystal, so that effectively photosensor arrays arranged in parallel is same In one expansion plane；Third establishment of coordinate system unit, for being built in expansion plane for the photosensor arrays after expansion Vertical plane coordinate system, the horizontal axis of plane coordinate system and the bottom surface of scintillation crystal are parallel, the longitudinal axis and scintillation crystal of plane coordinate system Central axis it is parallel；Third distribution of photons determination unit, for according to energy information corresponding with effective photosensor arrays Determine specific distribution of photons corresponding with the photosensor arrays after expansion；Ordinate computing unit, for according to specific light Son distribution calculates the ordinate of projected position of the response location in plane coordinate system；And first axis coordinate determination unit, For determining that the ordinate of projected position is the axial coordinate of response location.

According to embodiments of the present invention, position determination submodule may further include: axial error determination unit, be used for root Specific axial error corresponding with the axial coordinate of response location is determined according to axial error distribution function, wherein axial error point Cloth function is known；And axial coordinate calibration unit, for the axial coordinate according to specific axial error to response location It is calibrated.

According to embodiments of the present invention, device 1500 may further include: the 4th receiving module, for axial for control Each control axial coordinate in coordinate set, the pole in radioactive source and the line and cylindrical coordinate system that meet detecting module are sat In the case that the distance between plane is equal to the control axial coordinate where mark system, the photoelectric sensing met in detecting module is received Detected by device, the 5th test gamma photons react generated visible light with the test crystal for meeting in detecting module 5th temporal information of son, and receive detected by least one photosensor arrays, the 6th test gamma photons and dodge Bright crystal react caused by optical photon third energy information and the 6th temporal information；Third meets event detection mould Block, for being believed according to the 5th temporal information and the 6th time for each control axial coordinate in control axial coordinate set Breath carries out meeting event detection；Axial coordinate computing module is tested, for for each control in control axial coordinate set Axial coordinate, when meet event detection outcome instruction meet event occur when, according to third energy corresponding with event is met Portion of energy information in information, corresponding with effective photosensor arrays calculates the 6th test gamma photons in scintillation crystal The axial coordinate of interior test reaction position；Axial error computing module is tested, for in control axial coordinate set Each control axial coordinate, the axial coordinate for calculating test reaction position compare the difference between axial coordinate with this, to obtain Test axial error；Axial error distribution function determining module, for according to for all right in control axial coordinate set Axial error distribution function is determined according to the test axial error that axial coordinate obtains respectively, wherein radioactive source and meets detection mould Block is synchronously moved along the direction parallel with the central axis of scintillation crystal, to obtain each control in control axial coordinate set Axial coordinate.

According to embodiments of the present invention, position determination submodule may include: the second axial coordinate determination unit, be used for basis The energy information determination of optical photon detected by the photosensor arrays coupled with two bottom surfaces of scintillation crystal is reacted The axial coordinate of position.

Those skilled in the art according to the description of each step above with respect to response location localization method and attached drawing 1 to 14, it is to be understood that embodiment and its advantage of response location positioning device 1500 disclosed herein etc., for sake of simplicity, herein It is repeated not to this.

The present invention has been explained by the above embodiments, but it is to be understood that, above-described embodiment is only intended to The purpose of citing and explanation, is not intended to limit the invention to the scope of the described embodiments.Furthermore those skilled in the art It is understood that the present invention is not limited to the above embodiments, introduction according to the present invention can also be made more kinds of member Variants and modifications, all fall within the scope of the claimed invention for these variants and modifications.Protection scope of the present invention by The appended claims and its equivalent scope are defined.

Claims (26)

1. a kind of response location localization method, comprising:

It receives detected by least one photosensor arrays for couple with scintillation crystal, gamma photons and flashing crystalline substance Body react caused by optical photon energy information, wherein the scintillation crystal is integrated scintillation crystal, and The scintillation crystal has through-hole, and for accommodating object to be imaged, the gamma photons pass through described to be imaged the through-hole The esoteric positron annihilation effect of object generates；And

According to received energy information determine response location of the gamma photons in the scintillation crystal, wherein it is described Response location is indicated using cylindrical coordinate system；

Wherein, it is described according to received energy information determine response location packet of the gamma photons in the scintillation crystal It includes:

Determine effective photosensor arrays at least one described photosensor arrays；And

The response location is determined according at least to energy information corresponding with the effective photosensor arrays；

Wherein, described to determine the response location packet according at least to energy information corresponding with the effective photosensor arrays It includes:

Effective photosensor arrays are unfolded along the direction parallel with the bottom surface of the scintillation crystal, so that described have Photosensor arrays arranged in parallel is imitated in same expansion plane；

Plane coordinate system is established for the photosensor arrays after expansion in the expansion plane, the plane coordinate system Horizontal axis is parallel with the bottom surface of the scintillation crystal, and the longitudinal axis of the plane coordinate system is parallel with the central axis of the scintillation crystal；

According to the determining photoelectric sensor battle array with after the expansion of energy information corresponding with the effective photosensor arrays Arrange corresponding specific distribution of photons；

The abscissa of projected position of the response location in the plane coordinate system is calculated according to the specific distribution of photons； And

The angular coordinate of the response location is calculated according to the abscissa of the geometrical characteristic of the scintillation crystal and the projected position.

2. the method according to claim 1, wherein the scintillation crystal is polygon prism shape crystal and described logical Hole is circular through hole or polygonal through hole, effective photoelectric sensor at least one photosensor arrays described in the determination Array includes:

The photoelectric sensor battle array with the side surface coupling of the scintillation crystal is selected from least one described photosensor arrays Column are used as effective photosensor arrays.

3. the method according to claim 1, wherein the scintillation crystal is cylinder crystal and the through-hole It is circular through hole or polygonal through hole, effective photoelectric sensor battle array at least one photosensor arrays described in the determination Column include:

The photosensor arrays coupled with the inner wall of the through-hole are selected to make from least one described photosensor arrays For effective photosensor arrays.

4. the method according to claim 1, wherein the geometrical characteristic according to the scintillation crystal and described The angular coordinate that the abscissa of projected position calculates the response location is implemented based on following formula:

Wherein, θ indicates the angular coordinate of the response location, and fix indicates to be rounded downwards, and mod expression takes the remainder, and x indicates the throwing The abscissa of shadow position, a indicate that the matrix number of effective photosensor arrays, L indicate effective photoelectric sensor The side length on each photosensor arrays in array, parallel with the bottom surface of scintillation crystal sides, R indicate the flashing Crystal, the face that is coupled with effective photosensor arrays the polygon of acquisition is projected on the bottom surface of the scintillation crystal Circumscribed circle radius,

Wherein, the cylindrical coordinate system is using the central point of the bottom surface of the scintillation crystal as origin and the cylindrical coordinate system In polar coordinate system where plane it is parallel with the bottom surface of the scintillation crystal.

5. the method according to claim 1, wherein it is described according at least to effective photosensor arrays Corresponding energy information determines that the response location further comprises:

Special angle error corresponding with the angular coordinate of the response location is determined according to angular error distribution function, wherein institute It is known for stating angular error distribution function；And

The angular coordinate of the response location is calibrated according to the special angle error.

6. according to the method described in claim 5, it is characterized in that, it is described according to angular error distribution function it is determining with it is described Before the corresponding special angle error of the angular coordinate of response location, the method further includes:

For control each of angle set to irradiation angle,

In radioactive source and meet the angle between the line of detecting module and the polar axis of the polar coordinate system in the cylindrical coordinate system Equal to this to irradiation angle in the case where, meet detected by the photoelectric sensor in detecting module described in reception, the first test Gamma photons and the test crystal met in detecting module react the first time information of generated optical photon, And receive detected by least one described photosensor arrays, the second test gamma photons and the scintillation crystal occur The first energy information and the second temporal information of optical photon caused by reacting；

It carries out meeting event detection according to the first time information and second temporal information；

When meeting event detection outcome instruction and meeting event and occur, believe according to corresponding first energy of event that meets Portion of energy information in breath, corresponding with the effective photosensor arrays calculates the second test gamma photons in institute State the angular coordinate of the test reaction position in scintillation crystal；

Calculate the test reaction position angular coordinate and this to the difference between irradiation angle, to obtain test angle error；

The angle is determined according to the test angle error obtained respectively for all pairs of irradiation angles in the control angle set Error distribution function is spent,

Wherein, the radioactive source and the center axis rotation for meeting detecting module and synchronously surrounding the scintillation crystal, to obtain Obtain each pair of irradiation angle in the control angle set.

7. method according to any one of claims 1 to 6, which is characterized in that it is described according at least to effective photoelectricity The corresponding energy information of sensor array determines that the response location includes:

Effective photosensor arrays are unfolded along the direction parallel with the bottom surface of the scintillation crystal, so that described have Photosensor arrays arranged in parallel is imitated in same expansion plane；

Plane coordinate system is established for the photosensor arrays after expansion in the expansion plane, the plane coordinate system Horizontal axis is parallel with the bottom surface of the scintillation crystal, and the longitudinal axis of the plane coordinate system is parallel with the central axis of the scintillation crystal；

According to the determining photoelectric sensor battle array with after the expansion of energy information corresponding with the effective photosensor arrays Arrange corresponding specific distribution of photons；And

Standard deviation and peak value based on the specific distribution of photons are to pair being used to indicate between standard deviation, peak value and radial coordinate The look-up table that should be related to is inquired, to obtain the radial coordinate of the response location.

8. the method according to the description of claim 7 is characterized in that in the standard deviation based on the specific distribution of photons and Before peak value inquires the look-up table for the corresponding relationship being used to indicate between standard deviation, peak value and radial coordinate, the side Method further comprises:

For control radius set in each control radius,

Radioactive source and meet the line of detecting module between the origin in the cylindrical coordinate system at a distance from be equal to the control Meet detected by the photoelectric sensor in detecting module in the case where radius, described in reception, third test gamma photons with The test crystal met in detecting module react caused by optical photon third temporal information, and receive described in Detected by least one photosensor arrays, the 4th test gamma photons and produced by the scintillation crystal reacts Optical photon the second energy information and the 4th temporal information；

It carries out meeting event detection according to the third temporal information and the 4th temporal information；

When meeting event detection outcome instruction and meeting event and occur, believe according to corresponding second energy of event that meets The determining photoelectric sensor with after the expansion of portion of energy information in breath, corresponding with the effective photosensor arrays The corresponding test distribution of photons of array；

According to the standard deviation of the test distribution of photons obtained respectively for all control radiuses in the control radius set, peak Value and corresponding control radius determine the look-up table,

Wherein, the radioactive source and it is described meet detecting module synchronously moving radially along the scintillation crystal, with obtain Each control radius in the control radius set.

9. method according to any one of claims 1 to 6, which is characterized in that it is described according at least to effective photoelectricity The corresponding energy information of sensor array determines that the response location includes:

The specific bottom surface of the scintillation crystal is divided into predetermined number destination region；

Corresponding region is selected from the predetermined number destination region according to the angular coordinate of the response location；

Select the photosensor arrays coupled with the corresponding region；And

According to the energy information meter of optical photon detected by the photosensor arrays coupled with the corresponding region Calculate the radial coordinate of the response location.

10. method according to any one of claims 1 to 6, which is characterized in that it is described according at least to effective photoelectricity The corresponding energy information of sensor array determines that the response location includes:

Effective photosensor arrays are unfolded along the direction parallel with the bottom surface of the scintillation crystal, so that described have Photosensor arrays arranged in parallel is imitated in same expansion plane；

Plane coordinate system is established for the photosensor arrays after expansion in the expansion plane, the plane coordinate system Horizontal axis is parallel with the bottom surface of the scintillation crystal, and the longitudinal axis of the plane coordinate system is parallel with the central axis of the scintillation crystal；

According to the determining photoelectric sensor battle array with after the expansion of energy information corresponding with the effective photosensor arrays Arrange corresponding specific distribution of photons；

The ordinate of projected position of the response location in the plane coordinate system is calculated according to the specific distribution of photons； And

The ordinate for determining the projected position is the axial coordinate of the response location.

11. according to the method described in claim 10, it is characterized in that, it is described according at least to effective photoelectric sensor battle array It arranges corresponding energy information and determines that the response location further comprises:

Specific axial error corresponding with the axial coordinate of the response location is determined according to axial error distribution function, wherein The axial error distribution function is known；And

The axial coordinate of the response location is calibrated according to the specific axial error.

12. according to the method for claim 11, which is characterized in that described according to the determination of axial error distribution function and institute Before the corresponding specific axial error of axial coordinate for stating response location, the method further includes:

For control axial coordinate set in each control axial coordinate,

In radioactive source and where meeting polar coordinate system in the line and the cylindrical coordinate system of detecting module between plane away from In the case where being equal to the control axial coordinate, meet detected by the photoelectric sensor in detecting module described in reception, the Five test gamma photons and the test crystal met in detecting module react generated optical photon the 5th when Between information, and receive detected by least one described photosensor arrays, the 6th test gamma photons and the flashing Crystal react caused by optical photon third energy information and the 6th temporal information；

It carries out meeting event detection according to the 5th temporal information and the 6th temporal information；

When meeting event detection outcome instruction and meeting event and occur, believe according to the corresponding third energy of event that meets Portion of energy information in breath, corresponding with the effective photosensor arrays calculates the 6th test gamma photons in institute State the axial coordinate of the test reaction position in scintillation crystal；

The axial coordinate for calculating the test reaction position compares difference between axial coordinate with this, with obtain test it is axial accidentally Difference；

It is true according to the test axial error obtained respectively for all control axial coordinates in the control axial coordinate set The fixed axial error distribution function,

Wherein, the radioactive source and described meet detecting module synchronously along the side parallel with the central axis of the scintillation crystal To movement, to obtain each control axial coordinate in the control axial coordinate set.

13. method according to any one of claims 1 to 6, which is characterized in that it is described according at least to effective photoelectricity The corresponding energy information of sensor array determines that the response location includes:

The energy of optical photon detected by the photosensor arrays coupled according to two bottom surfaces with the scintillation crystal Information determines the axial coordinate of the response location.

14. a kind of response location positioning device, comprising:

First receiving module, for receiving detected by least one photosensor arrays coupled with scintillation crystal, gal Ma photon and the scintillation crystal react caused by optical photon energy information, wherein the scintillation crystal is one The scintillation crystal of body, and the scintillation crystal has through-hole, and the through-hole is for accommodating object to be imaged, the gamma light Son in the esoteric positron annihilation effect of object to be imaged by generating；And

Position determination module, for according to received energy information determine that the gamma photons are anti-in the scintillation crystal Answer position, wherein the response location is indicated using cylindrical coordinate system；

Wherein, the position determination module includes:

Sensor determines submodule, for determining effective photoelectric sensor battle array at least one described photosensor arrays Column；And

Position determination submodule, described in being determined according at least to energy information corresponding with the effective photosensor arrays Response location；

Wherein, the position determination submodule includes:

First expansion unit, for by the effective photosensor arrays along the direction parallel with the bottom surface of the scintillation crystal Expansion, so that effective photosensor arrays arranged in parallel is in same expansion plane；

First establishment of coordinate system unit, for establishing plane for the photosensor arrays after expansion in the expansion plane Coordinate system, the horizontal axis of the plane coordinate system is parallel with the bottom surface of the scintillation crystal, the longitudinal axis of the plane coordinate system and institute The central axis for stating scintillation crystal is parallel；

First distribution of photons determination unit, for according to energy information corresponding with the effective photosensor arrays determine with The corresponding specific distribution of photons of photosensor arrays after the expansion；

Abscissa computing unit, for calculating the response location in the plane coordinate system according to the specific distribution of photons Projected position abscissa；And

Angular coordinate computing unit, for calculating institute according to the geometrical characteristic of the scintillation crystal and the abscissa of the projected position State the angular coordinate of response location.

15. device according to claim 14, which is characterized in that the scintillation crystal is polygon prism shape crystal and described Through-hole is circular through hole or polygonal through hole, and the sensor determines that submodule includes:

First selecting unit, for selecting the side coupling with the scintillation crystal from least one described photosensor arrays The photosensor arrays of conjunction are as effective photosensor arrays.

16. device according to claim 14, which is characterized in that the scintillation crystal is cylinder crystal and described logical Hole is circular through hole or polygonal through hole, and the sensor determines that submodule includes:

Second selecting unit, for from least one described photosensor arrays selection couple with the inner wall of the through-hole Photosensor arrays are as effective photosensor arrays.

17. device according to claim 14, which is characterized in that the angular coordinate computing unit includes computation module, is used In the angular coordinate for calculating the response location based on following formula:

Wherein, θ indicates the angular coordinate of the response location, and fix indicates to be rounded downwards, and mod expression takes the remainder, and x indicates the throwing The abscissa of shadow position, a indicate that the matrix number of effective photosensor arrays, L indicate effective photoelectric sensor The side length on each photosensor arrays in array, parallel with the bottom surface of scintillation crystal sides, R indicate the flashing Crystal, the face that is coupled with effective photosensor arrays the polygon of acquisition is projected on the bottom surface of the scintillation crystal Circumscribed circle radius,

Wherein, the cylindrical coordinate system is using the central point of the bottom surface of the scintillation crystal as origin and the cylindrical coordinate system In polar coordinate system where plane it is parallel with the bottom surface of the scintillation crystal.

18. device according to claim 14, which is characterized in that the position determination submodule further comprises:

Angular error determination unit, for corresponding with the angular coordinate of the response location according to the determination of angular error distribution function Special angle error, wherein the angular error distribution function is known；And

Angular coordinate calibration unit, for being calibrated according to the special angle error to the angular coordinate of the response location.

19. device according to claim 18, which is characterized in that described device further comprises:

Second receiving module, for, to irradiation angle, in radioactive source and meeting detecting module for each of control angle set Line and the polar coordinate system in the cylindrical coordinate system polar axis between angle be equal to this to irradiation angle in the case where, receive Detected by the photoelectric sensor met in detecting module, the first test gamma photons and described meet in detecting module Test crystal react caused by optical photon first time information, and receive at least one described photoelectric sensor Detected by array, second test gamma photons and the scintillation crystal react caused by optical photon the first energy Measure information and the second temporal information；

First meets event checking module, for being directed to each of control angle set to irradiation angle, when according to described first Between information and second temporal information carry out meeting event detection；

Test angles coordinate calculation module, for for control each of angle set to irradiation angle, when meeting event detection knot When fruit instruction meets event and occurs, according to it is described meet it is in corresponding first energy information of event, with effective light The corresponding portion of energy information of electric transducer array calculates test of the second test gamma photons in the scintillation crystal The angular coordinate of response location；

Test angle error calculating module, for, to irradiation angle, it is anti-to calculate the test for each of control angle set Answer position angular coordinate and this to the difference between irradiation angle, to obtain test angle error；

Angular error distribution function determining module, for according to for all pairs of irradiation angles difference in the control angle set The test angle error of acquisition determines the angular error distribution function,

Wherein, the radioactive source and the center axis rotation for meeting detecting module and synchronously surrounding the scintillation crystal, to obtain Obtain each pair of irradiation angle in the control angle set.

20. 4 to 19 described in any item devices according to claim 1, which is characterized in that the position determination submodule includes:

Second expansion unit, for by the effective photosensor arrays along the direction parallel with the bottom surface of the scintillation crystal Expansion, so that effective photosensor arrays arranged in parallel is in same expansion plane；

Second establishment of coordinate system unit, for establishing plane for the photosensor arrays after expansion in the expansion plane Coordinate system, the horizontal axis of the plane coordinate system is parallel with the bottom surface of the scintillation crystal, the longitudinal axis of the plane coordinate system and institute The central axis for stating scintillation crystal is parallel；

Second distribution of photons determination unit, for according to energy information corresponding with the effective photosensor arrays determine with The corresponding specific distribution of photons of photosensor arrays after the expansion；And query unit, for being based on the specific light The standard deviation and peak value of son distribution carry out the look-up table for the corresponding relationship being used to indicate between standard deviation, peak value and radial coordinate Inquiry, to obtain the radial coordinate of the response location.

21. device according to claim 20, which is characterized in that described device further comprises:

Third receiving module, in radioactive source and meeting detecting module for each control radius in control radius set The distance between line and origin in the cylindrical coordinate system be equal in the case that this compares radius, spy is met described in reception It surveys detected by the photoelectric sensor in module, third tests gamma photons and the test crystal met in detecting module The third temporal information of optical photon caused by reacting, and receive at least one described photosensor arrays and detected To, the 4th test gamma photons and the scintillation crystal react caused by optical photon the second energy information and the Four temporal informations；

Second meets event checking module, each control radius for being directed in control radius set, when according to the third Between information and the 4th temporal information carry out meeting event detection；

Distribution of photons determining module is tested, for for each control radius in control radius set, when meeting event detection When resulting indicator closes event and occurs, according to it is described meet it is in corresponding second energy information of event, with it is described effectively The corresponding portion of energy information of photosensor arrays determines test light corresponding with the photosensor arrays after the expansion Son distribution；

Look-up table determining module, for according to the test obtained respectively for all control radiuses in the control radius set Standard deviation, peak value and the corresponding control radius of distribution of photons determine the look-up table,

Wherein, the radioactive source and it is described meet detecting module synchronously moving radially along the scintillation crystal, with obtain Each control radius in the control radius set.

22. 4 to 19 described in any item devices according to claim 1, which is characterized in that the position determination submodule includes:

Area division unit, for the specific bottom surface of the scintillation crystal to be divided into predetermined number destination region；

Area selecting unit selects corresponding area for the angular coordinate according to the response location from the predetermined number destination region Domain；

Sensor selection unit, the photosensor arrays for selecting to couple with the corresponding region；And

Radial coordinate computing unit, for according to detected by the photosensor arrays coupled with the corresponding region The energy information of optical photon calculates the radial coordinate of the response location.

23. 4 to 19 described in any item devices according to claim 1, which is characterized in that the position determination submodule includes:

Unit is unfolded in third, for by the effective photosensor arrays along the direction parallel with the bottom surface of the scintillation crystal Expansion, so that effective photosensor arrays arranged in parallel is in same expansion plane；

Third establishment of coordinate system unit, for establishing plane for the photosensor arrays after expansion in the expansion plane Coordinate system, the horizontal axis of the plane coordinate system is parallel with the bottom surface of the scintillation crystal, the longitudinal axis of the plane coordinate system and institute The central axis for stating scintillation crystal is parallel；

Third distribution of photons determination unit, for according to energy information corresponding with the effective photosensor arrays determine with The corresponding specific distribution of photons of photosensor arrays after the expansion；

Ordinate computing unit, for calculating the response location in the plane coordinate system according to the specific distribution of photons Projected position ordinate；And

First axis coordinate determination unit, for determining that the ordinate of the projected position is that the axial of the response location is sat Mark.

24. device according to claim 23, which is characterized in that the position determination submodule further comprises:

Axial error determination unit, for corresponding with the axial coordinate of the response location according to the determination of axial error distribution function Specific axial error, wherein the axial error distribution function is known；And

Axial coordinate calibration unit, for carrying out school according to axial coordinate of the specific axial error to the response location It is quasi-.

25. device according to claim 24, which is characterized in that described device further comprises:

4th receiving module, in radioactive source and meeting for each control axial coordinate in control axial coordinate set The line of detecting module is equal to this with the distance between plane where the polar coordinate system in the cylindrical coordinate system and compares axial seat Meet detected by the photoelectric sensor in detecting module in the case of target, described in reception, the 5th test gamma photons and institute State meet the crystal of the test in detecting module react caused by optical photon the 5th temporal information, and receive it is described extremely Detected by few photosensor arrays, caused by the 6th test gamma photons and the scintillation crystal react The third energy information and the 6th temporal information of optical photon；

Third meets event checking module, each control axial coordinate for being directed in control axial coordinate set, according to institute It states the 5th temporal information and the 6th temporal information carries out meeting event detection；

Axial coordinate computing module is tested, for for each control axial coordinate in control axial coordinate set, when meeting When event detection outcome instruction meets event and occurs, according to it is described meet it is in the corresponding third energy information of event, with The corresponding portion of energy information of effective photosensor arrays calculates the 6th test gamma photons in the flashing crystalline substance The axial coordinate of intracorporal test reaction position；

Axial error computing module is tested, for calculating institute for each control axial coordinate in control axial coordinate set The axial coordinate for stating test reaction position compares difference between axial coordinate with this, to obtain test axial error；

Axial error distribution function determining module, for according to axial for all controls in the control axial coordinate set The test axial error that coordinate obtains respectively determines the axial error distribution function,

Wherein, the radioactive source and described meet detecting module synchronously along the side parallel with the central axis of the scintillation crystal To movement, to obtain each control axial coordinate in the control axial coordinate set.

26. 4 to 19 described in any item devices according to claim 1, which is characterized in that the position determination submodule includes:

Second axial coordinate determination unit, for according to the photosensor arrays coupled with two bottom surfaces of the scintillation crystal The energy information of detected optical photon determines the axial coordinate of the response location.