CN107874773A - Photon detection method, apparatus, equipment and system and storage medium - Google Patents

Photon detection method, apparatus, equipment and system and storage medium Download PDF

Info

Publication number
CN107874773A
CN107874773A CN201710975365.XA CN201710975365A CN107874773A CN 107874773 A CN107874773 A CN 107874773A CN 201710975365 A CN201710975365 A CN 201710975365A CN 107874773 A CN107874773 A CN 107874773A
Authority
CN
China
Prior art keywords
energy
sensor
photon
reading circuit
energy signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710975365.XA
Other languages
Chinese (zh)
Other versions
CN107874773B (en
Inventor
谢思维
杨静梧
黄秋
龚政
应高阳
苏志宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sino Science And Technology (shenzhen) Co Ltd
Original Assignee
Sino Science And Technology (shenzhen) Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sino Science And Technology (shenzhen) Co Ltd filed Critical Sino Science And Technology (shenzhen) Co Ltd
Priority to CN201710975365.XA priority Critical patent/CN107874773B/en
Publication of CN107874773A publication Critical patent/CN107874773A/en
Priority to PCT/CN2018/110255 priority patent/WO2019076268A1/en
Application granted granted Critical
Publication of CN107874773B publication Critical patent/CN107874773B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data

Landscapes

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

Abstract

The present invention provides a kind of photon detection method, apparatus, equipment and system and storage medium.Method includes:The energy signal for the second number that the energy signal of the first number and the independent reading circuit of the second number that the shared reading circuit for the first number that reception is connected with sensor array exports respectively export respectively, wherein, sensor array column average is divided at least two sensor regions, first number is equal to the number of the sensor in each sensor region, and the shared reading circuit of the first number connects all the sensors in each sensor region correspondingly, the single sensor in each individually reading circuit connection sensor array;The regularity of energy distribution of the energy signal of energy signal and the second number based on the first number determines the reaction projected position of high-energy photon.The shared reading circuit and the electric signal of independent reading circuit pickup output that propagation characteristic and distribution character of this method based on scintillation photons are set, to realize that passage reduces.

Description

Photon detection method, apparatus, equipment and system and storage medium
Technical field
The present invention relates to Positron emission tomography field, in particular it relates to a kind of photon detection method, apparatus, equipment and System and storage medium.
Background technology
In recent years, when Positron emission tomography system is by photomultiplier (Photomultiplier Tube, abbreviation PMT) In generation, is increasingly turned to silicon photomultiplier (Silicon Photomultiplier, the abbreviation SiPM) epoch.Positron emission tomography Full name is positron e mission computed tomography (Positron Emission Computed Tomography, abbreviation PET), it is a kind of technology that human body or animal body internal structure are shown using Radioactive isotope method, is that nuclear medicine is ground Study carefully the Main Means with clinical diagnosis.
SiPM sizes are smaller, and detection efficiency is high, so that detector is compact-sized, system sensitivity is high.But exactly because Small for SiPM sizes, the SiPM quantity of the scintillation crystal coupling of same cross-sectional area will be far longer than PMT quantity.Using aperture as Exemplified by 76cm monocyclic whole body PET system, the lighting area of scintillation crystal is about 125cm2, it is necessary to couple 196 shore pines R9800PMT (light sensitive area 25mm), or the SiPM that 3400 sizes of coupling are 6mm.If all SiPM signal is all independent Read, then system channel number will increase by 17 times or so.So technological innovation of the PMT systems to SiPM systems, is PET system In signal read circuit bring certain challenge.
Accordingly, it is desirable to provide a kind of electron channel reduction technology, to solve at least in part on present in prior art State problem.
The content of the invention
In order to solve problems of the prior art at least in part, the present invention provides a kind of photon detection method, dress Put, equipment and system and storage medium.
According to an aspect of the present invention, there is provided a kind of photon detection method, including:Receive what is be connected with sensor array The energy signal for the first number that the shared reading circuit of first number exports respectively and the independent reading circuit of the second number The energy signal of the second number exported respectively, wherein, sensor array column average is divided at least two sensor regions, and first Number is equal to the number of the sensor in each sensor region, and the shared reading circuit of the first number connects correspondingly The all the sensors in each sensor region are connect, the single sensor in each individually reading circuit connection sensor array; And energy signal based on the first number and the regularity of energy distribution of the energy signal of the second number determine the anti-of high-energy photon Projected position is answered, wherein, reaction projected position is reaction position of the high-energy photon in the scintillation crystal coupled with sensor array Put projection on an array of sensors.
Exemplarily, photon detection method also includes:Energy signal based on the first number determines the energy of high-energy photon And/or arrival time.
Exemplarily, the regularity of energy distribution of the energy signal of energy signal and the second number based on the first number determines The reaction projected position of high-energy photon includes:By the energy signal of the first number and the energy signal of the second number input engineering Practise model to be analyzed, to obtain the position data of reaction projected position that machine learning model exports, on high-energy photon.
Exemplarily, photon detection method also includes:Photon reaction event simulation is carried out at sample response location, to obtain The sample energy signal of the first number corresponding with sample response location and the sample energy signal of the second number are obtained, wherein, with Sample projected position is known corresponding to sample response location;And sample energy signal and the second number with the first number Input of the sample energy signal as machine learning model, machine learning is used as using the position data on sample projected position The target output of model, is trained to machine learning model.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
According to another aspect of the present invention, there is provided a kind of photon detection apparatus, including:Sensor array, with scintillation crystal Coupling, reacted caused scintillation photons for detecting high-energy photon and scintillation crystal, wherein, sensor array column average is drawn It is divided at least two sensor regions;Reading circuit, it is connected with sensor array, for receiving the telecommunications of sensor array output Number and export the energy signal related to the energy of high-energy photon, wherein, shared reading of the reading circuit including the first number is electric Road and the independent reading circuit of the second number, the first number are equal to the number of the sensor in each sensor region, and the The shared reading circuit of one number connects all the sensors in each sensor region correspondingly, each individually to read electricity Single sensor in road connection sensor array;Process circuit, the shared reading circuit difference for receiving the first number are defeated The energy signal for the second number that the energy signal of the first number gone out and the independent reading circuit of the second number export respectively, And energy signal based on the first number and the regularity of energy distribution of the energy signal of the second number determine the reaction of high-energy photon Projected position, wherein, reaction projected position is the throwing of response location of the high-energy photon in scintillation crystal on an array of sensors Shadow.
Exemplarily, process circuit be additionally operable to the energy signal based on the first number determine high-energy photon energy and/or Arrival time.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
According to another aspect of the present invention, there is provided a kind of photon detection device, including:Receiving module, for receiving and passing The energy signal and the second number for the first number that the shared reading circuit of the first connected number of sensor array exports respectively The energy signal of the second number that exports respectively of independent reading circuit, wherein, sensor array column average is divided at least two Sensor region, the first number are equal to the number of the sensor in each sensor region, and the shared reading of the first number Circuit connects all the sensors in each sensor region correspondingly, each individually reading circuit connection sensor array In single sensor;And position determination module, believe for the energy signal based on the first number and the energy of the second number Number regularity of energy distribution determine the reaction projected position of high-energy photon, wherein, reaction projected position be high-energy photon with biography The projection of response location on an array of sensors in the scintillation crystal of sensor array coupling.
Exemplarily, photon detection device also includes:Energy or time determining module, for the energy based on the first number Signal determines energy and/or the arrival time of high-energy photon.
Exemplarily, position determination module includes:Input submodule, for the energy signal of the first number and second to be counted Purpose energy signal input machine learning model is analyzed, that machine learning model exports, on high-energy photon to obtain React the position data of projected position.
Exemplarily, photon detection device also includes:Analog module, for carrying out photon reaction at sample response location Event simulation, to obtain the sample energy of the sample energy signal of the first number corresponding with sample response location and the second number Signal, wherein, sample projected position corresponding with sample response location is known;And training module, for the first number The input of purpose sample energy signal and the sample energy signal of the second number as machine learning model, to be projected on sample The position data of position is exported as the target of machine learning model, and machine learning model is trained.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
According to another aspect of the present invention, there is provided a kind of photon detection system, including processor and memory, wherein, deposit Computer program instructions are stored with reservoir, are used to perform following steps when computer program instructions are run by processor:Receive The energy signal and second for the first number that the shared reading circuit for the first number being connected with sensor array exports respectively The energy signal for the second number that the independent reading circuit of number exports respectively, wherein, sensor array column average is divided at least Two sensor regions, the first number is equal to the number of the sensor in each sensor region, and the first number is shared Reading circuit connects all the sensors in each sensor region correspondingly, each individually reading circuit connection sensor Single sensor in array;And the energy signal based on the first number and the Energy distribution rule of the energy signal of the second number Rule determines the reaction projected position of high-energy photon, wherein, reaction projected position is that high-energy photon is coupling with sensor array The projection of response location on an array of sensors in scintillation crystal.
Exemplarily, it is additionally operable to perform following steps when computer program instructions are run by processor:Based on the first number Energy signal determine energy and/or the arrival time of high-energy photon.
Exemplarily, the letter of the energy based on the first number of used execution when computer program instructions are run by processor Number and the regularity of energy distribution of energy signal of the second number the step of determining the reaction projected position of high-energy photon include:By The energy signal of one number and the energy signal of the second number input machine learning model are analyzed, to obtain machine learning mould The position data of reaction projected position that type exports, on high-energy photon.
Exemplarily, it is additionally operable to perform following steps when computer program instructions are run by processor:Position is reacted in sample Put place and carry out photon reaction event simulation, to obtain the sample energy signal of corresponding with sample response location the first number and the The sample energy signal of two numbers, wherein, sample projected position corresponding with sample response location is known;And with first The input of the sample energy signal of number and the sample energy signal of the second number as machine learning model, to be thrown on sample The position data of shadow position is exported as the target of machine learning model, and machine learning model is trained.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
According to another aspect of the present invention, there is provided a kind of storage medium, store programmed instruction, program on a storage medium Instruction is operationally used to perform following steps:Receive the shared reading circuit difference for the first number being connected with sensor array The energy letter for the second number that the energy signal of the first number and the independent reading circuit of the second number of output export respectively Number, wherein, sensor array column average is divided at least two sensor regions, and the first number is equal in each sensor region The number of sensor, and the shared reading circuit of the first number connects all biographies in each sensor region correspondingly Sensor, each individually reading circuit connect the single sensor in sensor array;And the energy signal based on the first number The reaction projected position of high-energy photon is determined with the regularity of energy distribution of the energy signal of the second number, wherein, reaction projection position It is set to the projection of response location of the high-energy photon in the scintillation crystal coupled with sensor array on an array of sensors.
Exemplarily, described program instruction is operationally additionally operable to perform following steps:Energy letter based on the first number Number determine high-energy photon energy and/or arrival time.
Exemplarily, described program instructs the energy signal and second based on the first number of operationally used execution The regularity of energy distribution of the energy signal of number determines that the step of reaction projected position of high-energy photon includes:By the first number Energy signal and the energy signal of the second number input machine learning model are analyzed, to obtain machine learning model output , on high-energy photon reaction projected position position data.
Exemplarily, described program instruction is operationally additionally operable to perform following steps:Carried out at sample response location Photon reacts event simulation, to obtain the sample energy signal of corresponding with sample response location the first number and the second number Sample energy signal, wherein, sample projected position corresponding with sample response location is known;And the sample with the first number The input of this energy signal and the sample energy signal of the second number as machine learning model, with sample projected position Position data is exported as the target of machine learning model, and machine learning model is trained.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
Photon detection method, apparatus, equipment and system according to embodiments of the present invention and storage medium, using based on flicker The shared reading circuit and the telecommunications of independent reading circuit pickup output that the propagation characteristic and distribution character of photon are set Number so that the purpose of passage reduction can be reached while detector is not influenceed, be advantageous to be effectively reduced the work(of PET system Consumption and cost.
A series of concept of simplification is introduced in the content of the invention, these concepts will enter one in specific embodiment part Step describes in detail.Present invention part be not meant to attempt the key feature for limiting technical scheme claimed and Essential features, the protection domain for attempting to determine technical scheme claimed is not meant that more.
Below in conjunction with accompanying drawing, advantages and features of the invention are described in detail.
Brief description of the drawings
The drawings below of the present invention is used to understand the present invention in this as the part of the present invention.Shown in the drawings of this hair Bright embodiment and its description, for explaining the principle of the present invention.In the accompanying drawings,
Fig. 1 shows the schematic diagram of the scintillation photons radiation areas according to an example of the present invention;
Fig. 2 shows the total reflection schematic diagram according to the scintillation photons of an example of the present invention;
Fig. 3 shows single photon reaction event generation according to an embodiment of the invention, using optical software simulation Distribution map of the scintillation photons in SiPM arrays;
Fig. 4 shows the indicative flowchart of photon detection method according to an embodiment of the invention;
Fig. 5 shows sensor region division according to an embodiment of the invention and the individually signal of reading circuit arrangement Figure;
Fig. 6 shows that the reading order of the electric signal of the sensor array shown according to an embodiment of the invention, Fig. 5 is shown It is intended to;
Fig. 7 a-7d be shown respectively it is according to embodiments of the present invention, when the response location of gamma photons is four diverse locations When shared reading circuit output energy signal schematic diagram;
Fig. 8 a show sensor region in accordance with another embodiment of the present invention division and individually reading circuit arrangement Schematic diagram;
Fig. 8 b show sensor region in accordance with another embodiment of the present invention division and individually reading circuit arrangement Schematic diagram;
Fig. 9 shows the schematic block diagram of photon detection apparatus according to an embodiment of the invention;
Figure 10 shows the schematic block diagram of photon detection device according to an embodiment of the invention;And
Figure 11 shows the schematic block diagram of photon detection system according to an embodiment of the invention.
Embodiment
In the following description, there is provided substantial amounts of details is so as to thoroughly understand the present invention.However, this area skill Art personnel will be seen that, described below to only relate to presently preferred embodiments of the present invention, and the present invention can be without one or more so Details and be carried out.In addition, in order to avoid obscuring with the present invention, for some technical characteristics well known in the art not It is described.
In order to solve the above problems, the present invention proposes that a kind of photon detection method, apparatus, equipment and system and storage are situated between Matter.According to embodiments of the present invention, sensor array is divided into different regions, shares some shared readings between different zones Circuit, to reduce port number, (each reading circuit can be considered as a read-out channel, and reduction port number namely reduces constructed Reading circuit number).Which in addition, sensed to assist in identifying the energy signal of shared reading circuit output from Device, mark of the signal of some sensors as identification is individually read using independent reading circuit.With full line or permutation are sensed The passage reduction technology that the electric signal of device output is read with same reading circuit is compared, and passage reduction technology provided by the invention is examined Propagation and distribution situation of the scintillation photons (or optical photon) in sensor array are considered, therefore photosensitive region have been directed to Property is stronger, avoids the limitation due to electronic system, causes detector performance not play completely.Be described below in detail provided herein is Passage reduction technology theoretical foundation.
It is described herein so that high-energy photon is gamma photons as an example.Gamma photons in subject to be imaged by occurring Positron annihilation effect produce.Specifically, can be to treating when using Positron emission tomography device scan object to be imaged Imaging object internal injection contains radioisotopic tracer.When in the positive electron and subject to be imaged that isotope is released Negatron can bury in oblivion when meeting, thus produce a pair (difference 180 degrees) in opposite direction, energy be 511KeV gamma Photon.Caused a pair gamma photons in opposite direction incide in two relative positions in scintillation crystal respectively.Gamma The outer-shell electron that photon incides in scintillation crystal with atom interacts, and outer-shell electron is changed into after absorbing the energy of gamma photons Excitation state, the electronics of excitation state occur energy level transition and produce substantial amounts of scintillation photons.The sensor array coupled with scintillation crystal These scintillation photons can be detected, and when it detects scintillation photons, the optical signal of scintillation photons can be converted to electricity Signal will simultaneously change the electric signal output obtained.
Fig. 1 shows the schematic diagram of the scintillation photons radiation areas according to an example of the present invention.As shown in figure 1, flashing In the case that crystal is discrete crystal, scintillation photons caused by a photon reaction event are transferred to light guide from single small crystals In, the subregion that is subsequently incident in SiPM arrays.It may be noted that photon as described herein reaction event refer to high-energy photon with The event that scintillation crystal reacts.Due to photoconductive layer thinner thickness (be generally less than 5mm), thus scintillation photons in light guide simultaneously Do not scatter and just absorbed by SiPM completely.Therefore, in a photon reaction event, the swept areas of scintillation photons be it is certain, The thickness of photoconductive layer is thinner, and the swept area of scintillation photons is with regard to smaller.
In the case where scintillation crystal is continuous crystal, six faces of crystal polish, the face except coupling SiPM arrays Outside five faces all post high reflection film.For example, if scintillation crystal is yttrium luetcium silicate scintillation crystal (LYSO), its refractive index For 1.82.SiPM surface is glass, and its refractive index is 1.5.Therefore, it is from light that scintillation photons propagate to SiPM from scintillation crystal Close medium is to optically thinner medium, as shown in Figure 2, it may occur that total reflection phenomenon.Fig. 2 shows the passage of scintillation light according to an example of the present invention The total reflection schematic diagram of son.It is appreciated that the incidence angle when scintillation photons are incident to SiPM will be anti-when being more than critical angle Penetrate, so as to be injected into SiPM.That is, most of scintillation photons can be received by the SiPM in a certain region, only There are extremely least a portion of scintillation photons to be received after diffusing reflection by the SiPM outside the region.So no matter scintillation crystal is to connect Continuous crystal or discrete crystal, in a photon reaction event, only part SiPM can receive optical signal in SiPM arrays.
Below using size as 60mm × 60mm × 20mm continuous crystal couple 10 × 10 SiPM arrays (each SiPM is big Small is 6mm) exemplified by the present invention is discussed in detail.Optical software simulated photons reaction event can be utilized, follows the trail of all scintillation photons Track, until scintillation photons are absorbed.Fig. 3 shows list that is according to an embodiment of the invention, being simulated using optical software Distribution map of the scintillation photons in SiPM arrays caused by individual photon reaction event.In figure 3, two coordinates in horizontal direction SiPM sequence number is represented, the scintillation photons number that the coordinate representation SiPM on vertical direction is received.From figure 3, it can be seen that In a photon reaction event, only the SiPM of subregion can receive scintillation photons, wherein, the reaction of gamma photons Projection (react projected position) of the position (position that i.e. gamma photons react in scintillation crystal) on SiPM arrays The quantity for the scintillation photons that the SiPM at place is collected into is most, the energy highest of the optical signal detected, and position is projected apart from the reaction The energy for putting the optical signal that more remote SiPM is detected is lower.It is appreciated that when the response location of gamma photons is along level side To during movement, distributed areas of the scintillation photons in SiPM arrays can also move same distance along horizontal direction, and flash The distributional pattern of photon is basically unchanged.Therefore, it is not necessary to a reading circuit is all distributed for all the sensors, but can root According to distributional pattern of the scintillation photons in SiPM arrays, shared reading circuit is set.
Pass through above-mentioned analysis, it is known that such passage reduction technology be it is possible, i.e., it is big with the radiation scope of scintillation photons Small consistent sensor region shares same reading circuit, to realize the reduction of port number.
Fig. 4 shows the indicative flowchart of photon detection method 400 according to an embodiment of the invention.As shown in figure 4, Photon detection method 400 comprises the following steps.
In step S410, the shared reading circuit for the first number being connected with sensor array exports respectively first is received The energy signal for the second number that the energy signal of number and the independent reading circuit of the second number export respectively, wherein, pass Sensor array is averagely divided at least two sensor regions, and the first number is equal to the number of the sensor in each sensor region Mesh, and the shared reading circuit of the first number connects all the sensors in each sensor region correspondingly, each Single sensor in independent reading circuit connection sensor array.
It is true in step S420, the regularity of energy distribution of the energy signal of energy signal and the second number based on the first number Determine the reaction projected position of high-energy photon, wherein, reaction projected position is high-energy photon in the flicker coupled with sensor array The projection of response location on an array of sensors in crystal.
The example shown in Fig. 3 is continued to use below, i.e., the continuous crystal using size as 60mm × 60mm × 20mm couples 10 × 10 SiPM arrays (each SiPM sizes are 6mm) exemplified by, to describe photon detection method 400.With the response location of gamma photons Centered on, it is assumed that propagation distance of the scintillation photons in continuous crystal or light guide is no more than 2.5 sensors in sensor array The distance of the length of side.That is, in a photon reaction event, the SiPM for receiving optical signal is no more than 25.This In the case of, the row of channels that can be gone forward side by side according to the model split sensor region shown in Fig. 5 reduction.
Fig. 5 shows sensor region division according to an embodiment of the invention and the individually signal of reading circuit arrangement Figure.As shown in figure 5,100 sensors can be divided into tetra- regions of A, B, C, D along center line, each region includes 25 biographies Sensor.In each area, from left to right, 1.1,1.2,1.3 ... 5.5 are followed successively by sensor number from top to bottom.
For each sensor region, the shared reading circuit of the first number connects all sensings therein correspondingly Device so that different sensors region can share the shared reading circuit of the first number.Exemplarily, in different sensors region The consistent sensor of coordinate share same shared reading circuit, the so more convenient later stage is to detecting the sensings of scintillation photons Device is identified.It is of course also possible to the sensor for making the coordinate in different sensors region inconsistent shares same shared reading Circuit, this can set as needed.The coordinate of sensor as described herein is in sensor sensor region where it Position.In the case where sensor is numbered according to mode as shown in Figure 5, the coordinate of sensor can use sensor Numbering represents, and the identical finger numbering of described coordinate is identical.
In the example depicted in fig. 5, the shared reading circuit of the sensor (m, n) in tetra- regions of A, B, C, D (m=1, 2,3,4,5;N=1,2,3,4,5), totally 25 read-out channels.If it is appreciated that only set 25 shared reading circuits, In some cases, such as when on the response location sensor that numbering is 3.3 in tetra- regions of A, B, C, D respectively of gamma photons Fang Shi, the energy signal of 25 shared reading circuits outputs will be same or about, cannot be distinguished by response location be Above the sensor that numbering is 3.3 in which region.Therefore, it is necessary to set some independent reading circuits to be used for assisting in identifying, example 8 independent reading circuits in example as shown in Figure 5.As shown in figure 5, can be respectively 2.2,4.4 by numbering in region A, D, The electric signal list for the sensor (i.e. with the sensor in less darkened boxes circle in Fig. 5) that numbering is 2.4,4.2 in region B, C Solely read, totally 8 read-out channels.
Fig. 6 shows that the reading order of the electric signal of the sensor array shown according to an embodiment of the invention, Fig. 5 is shown It is intended to.According to mode as shown in Figure 6, the electric signal of 33 passage can be sequential read out.As described above, the effect of sensor It is that the optical signals of scintillation photons is subjected to opto-electronic conversion, therefore sensor output is the electric signal obtained after opto-electronic conversion. Reading circuit is used to handle the electric signal from sensor, and exports the optical signal received for representative sensor The energy signal of energy size.It is appreciated that energy size and output representated by the energy signal of each reading circuit output The quantity of scintillation photons received by electric signal to the sensor of the reading circuit is positively related.
A photon reaction event occurs, reading circuit can export the energy signal of 33 passages.In one example, Input that can be using the energy signal of 33 passages as machine learning model, machine learning model can export gamma photons React projected position.Machine learning model can according to first 8 independent read-out channels (A2.2, A4.4, B2.4, B4.2, C2.4, C4.2, D2.2 and D4.4) and rear 25 shared read-out channels output energy signal signal characteristic, to the anti-of gamma photons Projected position is answered to be decoded.
As described above, nearer from the reaction projected position of gamma photons, the energy for the optical signal that SiPM is detected is higher. Assuming that the energy for the optical signal that the SiPM at reaction projected position is detected is E3, the optical signal that adjacent 8 SiPM are detected Energy be E2, the energy for the optical signal that 16 SiPM of outermost are detected is E1.As independent read-out channel A2.2 and A4.4 During with signal, it may appear that such as tetra- kinds of situations of Fig. 7 a-7d.Fig. 7 a-7d be shown respectively it is according to embodiments of the present invention, when gamma light The schematic diagram of the energy signal of shared reading circuit output when the response location of son is four diverse locations.When generation is such as Fig. 7 a institutes During situation about showing, the signal of 33 passages is followed successively by:E2、E2、0、0、0、0、0、0、E1、E1、E1、E1、E1、E1、E2、E2、E2、 E1、E1、E2、E3、E2、E1、E1、E2、E2、E2、E1、E1、E1、E1、E1、E1.In first 8 independent read-out channels only the 1st and 2nd has an energy signal, and the power of energy signal is according to certain aligned transfer in rear 25 shared read-out channels, and shares Energy signal in read-out channel is all from region A.When three kinds of situations as shown in 7b, 7c and 7d occur, first 8 are individually read Going out in passage still only the 1st and the 2nd has an energy signal, and the arrangement of energy is different in rear 25 shared read-out channels. Pay attention to, in the case of shown in Fig. 7 b, the energy signal in 25 shared read-out channels there are four row (20) to come from region A, has One row (5) come from region B.Fig. 7 c are similar with Fig. 7 d, and energy signal comes from more than one region.Machine learning model passes through profit Trained with mass data, there is certain recognition capability to the difference of energy arrangement, gamma can be judged according to the arrangement of energy The reaction projected position of photon.
A kind of set-up mode of independent reading circuit is described above is, i.e., each sensor region sets two individually to read Circuit, 8 independent reading circuits are set altogether.However, above-mentioned example is not limitation of the present invention, independent reading circuit can have There are other rational numbers and set-up mode.For example, in another example, can be by numbering in tetra- regions of A, B, C, D The electric signal of 3.3 sensor is individually read, that is, sets one individually to read respectively at the central sensor in each region Go out circuit, four independent reading circuits are set altogether, so total read-out channel will be 29.It is arranged in each sensor region Independent reading circuit at the heart is enough to distinguish the electric signal from different sensors region.
Photon detection method according to embodiments of the present invention, is set using the propagation characteristic based on scintillation photons and distribution character The shared reading circuit put and the electric signal of independent reading circuit pickup output so that do not influenceing detector performance The purpose of passage reduction can be reached simultaneously, be advantageous to be effectively reduced the power consumption and cost of PET system.
According to embodiments of the present invention, photon detection method 400 also includes:Energy signal based on the first number determines high energy The energy of photon and/or arrival time.Energy corresponding to the energy signal of first number is the energy of high-energy photon.Example Property, the energy signal of the first number can be added, obtain total energy signal.Energy representated by total energy signal is big The small energy size equal to high-energy photon.Exemplarily, impulse level will can earliest occur in the energy signal of the first number Time be considered as arrival time of high-energy photon.Arrival time refers to the time that high-energy photon reaches detector, and it can utilize and pass Sensor array received is weighed to the time of scintillation photons.The determination mode of above-mentioned energy and arrival time are only examples, can be with Energy and/or the arrival time of high-energy photon are determined using other modes.High energy can be obtained by process circuit described below The energy information and temporal information of photon, data processing and image reconstruction are carried out to energy information and temporal information, to be treated The scan image of imaging object.
Because propagation characteristic of the reading circuit based on scintillation photons and distribution character are set, therefore exported based on reading circuit The energy of high-energy photon that obtains of energy signal measurement and arrival time to the with strong points of photosensitive region, accuracy is high.This Outside, due to the reduction of port number so that the measurement efficiency of energy and time are higher.
According to embodiments of the present invention, step S420 can include:By the energy signal of the first number and the energy of the second number Amount signal input machine learning model is analyzed, and is thrown with obtaining reaction that machine learning model exports, on high-energy photon The position data of shadow position.
Machine learning model can be any suitable intelligent algorithm model, and the present invention is not limited its specific category System.Exemplarily, machine learning model can use decision tree, SVMs, neutral net, AdaBoost algorithm models, Bayes classifier etc. is realized.In the following description, machine learning model will be described by taking convolutional neural networks as an example.
Machine learning model can be training in advance or be realized using known models.Machine learning model Reaction projected position can be classified.As described above, machine learning model is trained by using mass data, to energy The difference of arrangement has certain recognition capability, according to the arrangement of energy, can judge the reaction projected position of gamma photons.Cause This, the reaction projected position of gamma photons can simply, be quickly and accurately determined using machine learning algorithm.Machine learning mould Type output be position data, in one example, output be sensor numbering and sensor where sensor regions The numbering in domain.For example, continue to use the sensor array example shown in Fig. 5, Fig. 7 a-7d, it is assumed that the response location of gamma photons exists In the A of region numbering be 2.3 sensor surface, then machine learning model output can be indicated for region A and The data of numbering 2.3.
According to embodiments of the present invention, photon detection method 400 can also include:It is anti-that photon is carried out at sample response location Event simulation is answered, to obtain the sample energy of the sample energy signal of the first number corresponding with sample response location and the second number Signal is measured, wherein, sample projected position corresponding with sample response location is known;And the sample energy with the first number The input of signal and the sample energy signal of the second number as machine learning model, with the positional number on sample projected position According to the target output as machine learning model, machine learning model is trained.
Photon detection method 400 can also include the training step of machine learning model.Training can be a large amount of by collecting The energy signal of photon reaction event is realized known to reaction projected position.Photon reaction event simulation can use optics soft Part is realized.The model of PET system is built using optical software, and sets different sample response locations to be simulated respectively, with Obtain the energy signal that reading circuit exports under different sample response locations.Compare it is appreciated that becoming in sample response location During change, as far as possible sample projected position is also changed.Sample projected position be sample response location on an array of sensors Projection.In the sample projected position difference of gamma photons, the energy signal of reading circuit output is different.Exemplarily, for Each sample projected position, input that can be using energy signal corresponding to the sample projected position as convolutional neural networks, general Sample projected position exports as the target of convolutional neural networks, is trained by back-propagation method.
Machine learning algorithm is a kind of Spontaneous breathing trial method, it is possible to achieve extraordinary classifying quality.To engineering After habit model is trained, the model can be utilized, relatively accurately reaction is thrown according to the energy signal actually measured Shadow position is positioned.
According to embodiments of the present invention, the first number reacts caused flicker not less than high-energy photon and scintillation crystal The number for the sensor that greatest irradiation scope of the photon in sensor array is included.
Example as described above, it is assumed that the sensing that maximum propagation scope of the scintillation photons in sensor array is included The number of device is 5 × 5, then should be no less than by sharing reading circuit by 25.For example, with reference to Fig. 5-7d example showns, share and read The number of circuit can be 25.Certainly, the number for sharing reading circuit can be more, such as 36, that is to say, that can be by each Sensor region is arranged to 6 × 6 sizes.The sensor occupied area of shared reading circuit connection is preferably able to cover scintillation photons Maximum propagation scope in sensor array, to ensure the energy of scintillation photons caused by single photon reaction event substantially Reading circuit reading can be all shared.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.In the first number In the case that mesh is equal to N-1, each sensor region is at most assigned an independent reading circuit.The independent reading of second number Circuit is used to assist in identifying sensor in the case where which sensor region the energy signal that can not differentiate the first number comes from Region.
Assuming that sensor array column average is divided into 4 sensor regions (referring to Fig. 5,7a-7d), then the second number is minimum 3.In the case where independent reading circuit is 3,3 sensor regions are selected from 4 sensor regions, and are this 3 biographies An independent reading circuit is respectively distributed in sensor region.So, the energy signal read using 3 independent reading circuits can be with area Divide which sensor region energy signal comes from.It is appreciated that if individually the number of reading circuit is further reduced, such as only If having two, that may cause in some cases, the energy from unallocated individually two sensor regions of reading circuit Amount signal cannot be distinguished by out.For example, returning to Fig. 5, if setting two independent reading circuits respectively in A, B area, C, D region do not have There is independent reading circuit, that is, 2+25=27 read-out channel is set altogether, then it is 3.3 to be numbered when response location in C regions Sensor directly over when with directly over numbering is 3.3 in D regions sensor when, energy signal that 27 passages are read The regularity of distribution is same or about, will be unable to distinguish.
In addition to the independent read-out channel arrangement shown in Fig. 5-7d, the embodiment of the present invention also proposes other two kinds of sides Case, as shown in figs. 8 a and 8b.In the scheme shown in Fig. 8 a, A3.3, A3.5, B3.3, B5.3, C1.3, C3.3 are individually read, The signal of this 8 passages of D3.1, D3.3, the sensor (m, n) in tetra- regions of remaining A, B, C, D share a passage (m= 1,2,3,4,5;N=1,2,3,4,5), totally 33 passages.In the scheme shown in Fig. 8 b, A3.3, A4.5, B3.3 are individually read, The signal of this 8 passages of B5.2, C1.4, C3.3, D2.1, D3.3, the sensor (m, n) in tetra- regions of remaining A, B, C, D Share a passage (m=1,2,3,4,5;N=1,2,3,4,5), totally 33 passages.
Pay attention to, photon detection method provided by the invention goes for discrete crystal or continuous crystal, and unlimited combinations Body and crystal array size, SiPM and SiPM array sizes.Pay attention to, the present invention does not limit scintillator crystal materials, in text only with Exemplified by LYSO crystal.Scintillation crystal can be any suitable crystal, and the present invention is limited not to this.For example, scintillation crystal Can be bismuth germanium oxide (BGO), yttrium luetcium silicate (LYSO) or lanthanum bromide (LaBr3) etc..The present invention does not limit SiPM arrays and flicker The coupled modes of crystal, can be with direct-coupling, or is coupled by optical glue etc..
Sensor as described herein can be any suitable photoelectric sensor, such as PMT, SiPM or the pole of avalanche optoelectronic two Manage (APD) etc..Although the present invention is mainly described using SiPM as example herein, however, it is not limitation of the present invention, this hair It is bright to can apply to other similar detectors for needing passage reduction technology.
The number of shared reading circuit as described herein depends on the thickness and scintillation photons of crystal and photoconductive layer in crystalline substance The distance propagated in body or light guide.When propagation distance determines, the number of shared reading circuit can also determine therewith.Individually read The selection of circuit has a variety of modes, and the present invention is not intended to limit the selection, arrangement and number of independent reading circuit.The present invention proposes Be a kind of design, i.e., according to photon spread feature reading circuit is carried out regular shared and is separately provided.According to reading The signal characteristic of the energy signal of circuit output calculates the reaction projected position of gamma photons.
According to embodiments of the present invention, sensor array can be divided into 2 × 1,3 × 2,2 × 4,4 × 4 grade arbitrary number groups The sensor region of conjunction.The number and arrangement mode for the sensor region that sensor array includes can be set as needed, this Invention is limited not to this.For example, if scintillation crystal is 60mm × 60mm × 20mm continuous crystal, and scintillation photons exist Propagation distance in sensor is 10mm, then can use be divided into 3 × 3 sensor regions (each sensor region it is big Small is 20mm × 20mm) sensor array carry out photon detection, if the propagation distance of scintillation photons in the sensor is 15mm, then it can use the sensing for being divided into 2 × 2 sensor regions (size of each sensor region is 30mm × 30mm) Device array carries out photon detection.
According to a further aspect of the invention, there is provided a kind of photon detection apparatus.Fig. 9 shows according to an embodiment of the invention The schematic block diagram of photon detection apparatus 900.As shown in figure 9, photon detection apparatus 900 includes sensor array 910, reads electricity Road 920 and process circuit 930.Above with respect in the description of photon detection method 400, it has been described that scintillation crystal, sensing Device array, the circuit structure of reading circuit and operation principle.Process circuit is used for each step for realizing photon detection method 400 Suddenly/function.Those skilled in the art may be referred to understand photon detection apparatus above for the description of photon detection method 400 900 circuit structure and operation principle, is repeated no more.
Sensor array 910 couples with scintillation crystal, for detecting produced by high-energy photon and scintillation crystal react Scintillation photons, wherein, sensor array 910 is averagely divided at least two sensor regions.
Reading circuit 920 is connected with sensor array 910, for receiving the electric signal and defeated of the output of sensor array 910 Go out the energy signal related to the energy of high-energy photon, wherein, the shared reading circuit of reading circuit 920 including the first number and The independent reading circuit of second number, the first number are equal to the number of the sensor in each sensor region, and the first number The shared reading circuit of purpose connects all the sensors in each sensor region correspondingly, and each individually reading circuit connects Connect the single sensor in sensor array.
The energy signal for the first number that the shared reading circuit that process circuit 930 is used for the first number of reception exports respectively And second number the energy signal of the second number that exports respectively of independent reading circuit, and the letter of the energy based on the first number Number and the regularity of energy distribution of energy signal of the second number determine the reaction projected position of high-energy photon, wherein, reaction projection Position is projection of response location of the high-energy photon in scintillation crystal on sensor array 910.Process circuit 930 can be adopted Realized with any suitable hardware, software, and/or firmware.Exemplarily, process circuit 930 can use field programmable gate Array (FPGA), digital signal processor (DSP), CPLD (CPLD), micro-control unit (MCU) or center Processing unit (CPU) etc. is realized.
Exemplarily, process circuit 930 be additionally operable to the energy signal based on the first number determine high-energy photon energy and/ Or arrival time.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array 910 is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
According to a further aspect of the invention, there is provided a kind of photon detection device.Figure 10 is shown according to one implementation of the present invention The schematic block diagram of the photon detection device 1000 of example.
As shown in Figure 10, photon detection device 1000 according to embodiments of the present invention includes receiving module 1010 and position is true Cover half block 1020.The modules can perform each step of the photon detection method above in conjunction with Fig. 1-8b descriptions respectively Suddenly/function.Only the major function of each part of the photon detection device 1000 is described below, and retouched more than omitting The detail content stated.
The shared reading circuit that receiving module 1010 is used to receive the first number being connected with sensor array exports respectively The energy signal of the first number and the energy signal of the second number that exports respectively of the independent reading circuit of the second number, its In, sensor array column average is divided at least two sensor regions, and the first number is equal to the sensing in each sensor region The number of device, and the shared reading circuit of the first number connects all sensings in each sensor region correspondingly Device, each individually reading circuit connect the single sensor in sensor array.
Position determination module 1020 is used for the energy of the energy signal of energy signal and the second number based on the first number The regularity of distribution determines the reaction projected position of high-energy photon, wherein, reaction projected position be high-energy photon with sensor array The projection of response location on an array of sensors in the scintillation crystal of coupling.
Exemplarily, photon detection device 1000 can also include:Energy or time determining module (not shown), for base Energy and/or the arrival time of high-energy photon are determined in the energy signal of the first number.
Exemplarily, position determination module 1020 can include:Input submodule, for by the energy signal of the first number Analyzed with the energy signal of the second number input machine learning model, with obtain machine learning model output, on height The position data of the reaction projected position of energy photon.
Exemplarily, photon detection device 1000 can also include:Analog module, for being carried out at sample response location Photon reacts event simulation, to obtain the sample energy signal of corresponding with sample response location the first number and the second number Sample energy signal, wherein, sample projected position corresponding with sample response location is known;And training module, it is used for Using the sample energy signal of the first number and the sample energy signal of the second number as the input of machine learning model, with The position data of sample projected position is exported as the target of machine learning model, and machine learning model is trained.
Exemplarily, the first number is passing not less than high-energy photon and the scintillation crystal caused scintillation photons that react The number for the sensor that greatest irradiation scope in sensor array is included.
Exemplarily, the consistent sensor of coordinate in different sensors region shares same shared reading circuit.
Exemplarily, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
Those of ordinary skill in the art are it is to be appreciated that the list of each example described with reference to the embodiments described herein Member and algorithm steps, it can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually Performed with hardware or software mode, application-specific and design constraint depending on technical scheme.Professional and technical personnel Described function can be realized using distinct methods to each specific application, but this realization is it is not considered that exceed The scope of the present invention.
Figure 11 shows the schematic block diagram of photon detection system 1100 according to an embodiment of the invention.Photon detection System 1100 includes signal pickup assembly 1110, storage device 1120 and processor 1130.
Signal pickup assembly 1110 is used to gather the energy signal related to the energy of high-energy photon.Signal pickup assembly 1110 be optional, and photon detection system 1100 can not include signal pickup assembly 1110.In such a case, it is possible to utilize Other signal pickup assembly collecting energy signals, and the signal of collection is sent to photon detection system 1100.
The storage device 1120 is stored for realizing the corresponding step in photon detection method according to embodiments of the present invention Rapid computer program instructions.
The processor 1130 is used to run the computer program instructions stored in the storage device 1120, to perform root According to the corresponding steps of the photon detection method of the embodiment of the present invention, and for realizing photon detection according to embodiments of the present invention Receiving module 1010 and position determination module 1020 in device 1000.
In one embodiment, it is used to perform following steps when computer program instructions are run by processor 1130:Receive The energy signal and second for the first number that the shared reading circuit for the first number being connected with sensor array exports respectively The energy signal for the second number that the independent reading circuit of number exports respectively, wherein, sensor array column average is divided at least Two sensor regions, the first number is equal to the number of the sensor in each sensor region, and the first number is shared Reading circuit connects all the sensors in each sensor region correspondingly, each individually reading circuit connection sensor Single sensor in array;And the energy signal based on the first number and the Energy distribution rule of the energy signal of the second number Rule determines the reaction projected position of high-energy photon, wherein, reaction projected position is that high-energy photon is coupling with sensor array The projection of response location on an array of sensors in scintillation crystal.
In one embodiment, it is additionally operable to perform following steps when computer program instructions are run by processor 1130:Base Energy and/or the arrival time of high-energy photon are determined in the energy signal of the first number.
In one embodiment, used execution counts based on first when computer program instructions are run by processor 1130 The regularity of energy distribution of purpose energy signal and the energy signal of the second number determines the step of the reaction projected position of high-energy photon Suddenly include:The energy signal of first number and the energy signal of the second number input machine learning model are analyzed, to obtain Obtain the position data of reaction projected position that machine learning model exports, on high-energy photon.
In one embodiment, it is additionally operable to perform following steps when computer program instructions are run by processor 1130: Photon reaction event simulation is carried out at sample response location, to obtain the sample energy of the first number corresponding with sample response location Signal and the sample energy signal of the second number are measured, wherein, sample projected position corresponding with sample response location is known; And using the sample energy signal of the first number and the sample energy signal of the second number as the input of machine learning model, with Position data on sample projected position is exported as the target of machine learning model, and machine learning model is trained.
In one embodiment, the first number reacts caused passage of scintillation light not less than high-energy photon and scintillation crystal The number for the sensor that greatest irradiation scope of the son in sensor array is included.
In one embodiment, the consistent sensor of coordinate in different sensors region, which is shared, same shared reads electricity Road.
In one embodiment, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
In addition, according to embodiments of the present invention, a kind of storage medium is additionally provided, stores program on said storage Instruction, when described program instruction is run by computer or processor for performing the photon detection method of the embodiment of the present invention Corresponding steps, and for realizing the corresponding module in photon detection device according to embodiments of the present invention.The storage medium Such as can include the storage card of smart phone, the memory unit of tablet personal computer, the hard disk of computer, read-only storage (ROM), Erasable Programmable Read Only Memory EPROM (EPROM), portable compact disc read-only storage (CD-ROM), USB storage or on State any combination of storage medium.
In one embodiment, described program instruction can cause computer or place when by computer or processor operation Reason device realizes each functional module of photon detection device according to embodiments of the present invention, and/or can perform according to this The photon detection method of inventive embodiments.
In one embodiment, described program instruction is operationally used to perform following steps:Reception and sensor array The energy signal for the first number that the shared reading circuit of the first connected number exports respectively and the independent reading of the second number Go out the energy signal for the second number that circuit exports respectively, wherein, sensor array column average is divided at least two sensor regions Domain, the first number is equal to the number of the sensor in each sensor region, and the shared reading circuit of the first number is one by one All the sensors in correspondingly connected each sensor region, each individually reading circuit connect single in sensor array Sensor;And energy signal based on the first number and the regularity of energy distribution of the energy signal of the second number determine high energy light The reaction projected position of son, wherein, reaction projected position is high-energy photon in the scintillation crystal coupled with sensor array The projection of response location on an array of sensors.
In one embodiment, described program instruction is operationally additionally operable to perform following steps:Based on the first number Energy signal determines energy and/or the arrival time of high-energy photon.
In one embodiment, described program instructs the energy signal based on the first number of operationally used execution The step of reaction projected position that high-energy photon is determined with the regularity of energy distribution of the energy signal of the second number, includes:By first The energy signal of number and the energy signal of the second number input machine learning model are analyzed, to obtain machine learning model The position data of reaction projected position exporting, on high-energy photon.
In one embodiment, described program instruction is operationally additionally operable to perform following steps:In sample response location Place carries out photon reaction event simulation, to obtain the sample energy signal and second of the first number corresponding with sample response location The sample energy signal of number, wherein, sample projected position corresponding with sample response location is known;And with the first number The input of purpose sample energy signal and the sample energy signal of the second number as machine learning model, to be projected on sample The position data of position is exported as the target of machine learning model, and machine learning model is trained.
In one embodiment, the first number reacts caused passage of scintillation light not less than high-energy photon and scintillation crystal The number for the sensor that greatest irradiation scope of the son in sensor array is included.
In one embodiment, the consistent sensor of coordinate in different sensors region, which is shared, same shared reads electricity Road.
In one embodiment, the second number is not less than N-1, wherein, N is the number of at least two sensor regions.
Each module in photon detection system according to embodiments of the present invention can pass through reality according to embodiments of the present invention The processor computer program instructions that store in memory of operation of the electronic equipment of photon detection are applied to realize, or can be with The computer instruction stored in the computer-readable recording medium of computer program product according to embodiments of the present invention is counted Calculation machine is realized when running.
Although describe example embodiment by reference to accompanying drawing here, it should be understood that above-mentioned example embodiment is merely exemplary , and be not intended to limit the scope of the invention to this.Those of ordinary skill in the art can carry out various changes wherein And modification, it is made without departing from the scope of the present invention and spiritual.All such changes and modifications are intended to be included in appended claims Within required the scope of the present invention.
Those of ordinary skill in the art are it is to be appreciated that the list of each example described with reference to the embodiments described herein Member and algorithm steps, it can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually Performed with hardware or software mode, application-specific and design constraint depending on technical scheme.Professional and technical personnel Described function can be realized using distinct methods to each specific application, but this realization is it is not considered that exceed The scope of the present invention.
In several embodiments provided herein, it should be understood that disclosed apparatus and method, it can be passed through Its mode is realized.For example, apparatus embodiments described above are only schematical, for example, the division of the unit, only Only a kind of division of logic function, there can be other dividing mode when actually realizing, such as multiple units or component can be tied Another equipment is closed or is desirably integrated into, or some features can be ignored, or do not perform.
In the specification that this place provides, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention Example can be put into practice in the case of these no details.In some instances, known method, structure is not been shown in detail And technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the present invention and help to understand one or more of each inventive aspect, To the present invention exemplary embodiment description in, each feature of the invention be grouped together into sometimes single embodiment, figure, Or in descriptions thereof.However, the method for the invention should be construed to reflect following intention:It is i.e. claimed Application claims features more more than the feature being expressly recited in each claim.More precisely, such as corresponding power As sharp claim reflects, its inventive point is the spy that can use all features less than some disclosed single embodiment Levy to solve corresponding technical problem.Therefore, it then follows thus claims of embodiment are expressly incorporated in this specific Embodiment, wherein each claim is in itself as separate embodiments of the invention.
It will be understood to those skilled in the art that in addition to mutually exclusive between feature, any combinations pair can be used All features and so disclosed any method disclosed in this specification (including adjoint claim, summary and accompanying drawing) Or all processes or unit of equipment are combined.Unless expressly stated otherwise, this specification (including adjoint right will Ask, make a summary and accompanying drawing) disclosed in each feature can be replaced by the alternative features for providing identical, equivalent or similar purpose.
In addition, it will be appreciated by those of skill in the art that although some embodiments described herein include other embodiments In included some features rather than further feature, but the combination of the feature of different embodiments means in of the invention Within the scope of and form different embodiments.For example, in detail in the claims, embodiment claimed it is one of any Mode it can use in any combination.
The all parts embodiment of the present invention can be realized with hardware, or to be run on one or more processor Software module realize, or realized with combinations thereof.It will be understood by those of skill in the art that it can use in practice Microprocessor or digital signal processor (DSP) realize some moulds in photon detection device according to embodiments of the present invention The some or all functions of block.The present invention is also implemented as the part or complete for performing method as described herein The program of device (for example, computer program and computer program product) in portion.Such program for realizing the present invention can store On a computer-readable medium, or can the form with one or more signal.Such signal can be from internet Download and obtain on website, either provide on carrier signal or provided in the form of any other.
It should be noted that the present invention will be described rather than limits the invention for above-described embodiment, and ability Field technique personnel can design alternative embodiment without departing from the scope of the appended claims.In the claims, Any reference symbol between bracket should not be configured to limitations on claims.Word "comprising" does not exclude the presence of not Element or step listed in the claims.Word "a" or "an" before element does not exclude the presence of multiple such Element.The present invention can be by means of including the hardware of some different elements and being come by means of properly programmed computer real It is existing.In if the unit claim of equipment for drying is listed, several in these devices can be by same hardware branch To embody.The use of word first, second, and third does not indicate that any order.These words can be explained and run after fame Claim.
The foregoing is only a specific embodiment of the invention or the explanation to embodiment, protection of the invention Scope is not limited thereto, any one skilled in the art the invention discloses technical scope in, can be easily Expect change or replacement, should all be included within the scope of the present invention.Protection scope of the present invention should be with claim Protection domain is defined.

Claims (15)

1. a kind of photon detection method, including:
The energy signal for the first number that the shared reading circuit for the first number that reception is connected with sensor array exports respectively And second number the energy signal of the second number that exports respectively of independent reading circuit, wherein, the sensor array is put down At least two sensor regions are divided into, first number is equal to the number of the sensor in each sensor region, and And the shared reading circuit of first number connects all the sensors in each sensor region, Mei Gedan correspondingly Only reading circuit connects the single sensor in the sensor array;And
The regularity of energy distribution of the energy signal of energy signal and second number based on first number determines high energy The reaction projected position of photon, wherein, the reaction projected position is that the high-energy photon couples with the sensor array Scintillation crystal in projection of the response location on the sensor array.
2. photon detection method according to claim 1, it is characterised in that the photon detection method also includes:
Energy signal based on first number determines energy and/or the arrival time of the high-energy photon.
3. photon detection method according to claim 1, it is characterised in that the energy letter based on first number Number and the regularity of energy distribution of energy signal of second number determine that the reaction projected position of high-energy photon includes:
The energy signal of first number and the energy signal of second number input machine learning model are analyzed, To obtain the position data of reaction projected position that the machine learning model exports, on the high-energy photon.
4. photon detection method according to claim 3, it is characterised in that the photon detection method also includes:
Photon reaction event simulation is carried out at sample response location, to obtain the first number corresponding with the sample response location The sample energy signal of purpose sample energy signal and the second number, wherein, sample corresponding with the sample response location is thrown Shadow position is known;And
The machine learning is used as using the sample energy signal of first number and the sample energy signal of second number The input of model, exported using the position data on the sample projected position as the target of the machine learning model, it is right The machine learning model is trained.
5. photon detection method according to claim 1, it is characterised in that first number is not less than the high energy light It is sub react with the scintillation crystal caused by greatest irradiation scope of the scintillation photons in the sensor array wrapped The number of the sensor contained.
6. photon detection method according to claim 1, it is characterised in that coordinate in different sensors region is consistent Sensor shares same shared reading circuit.
7. photon detection method according to claim 1, it is characterised in that second number is not less than N-1, wherein, N For the number of at least two sensor region.
8. a kind of photon detection apparatus, including:
Sensor array, coupled with scintillation crystal, for detecting caused by high-energy photon and the scintillation crystal react Scintillation photons, wherein, the sensor array column average is divided at least two sensor regions;
Reading circuit, be connected with the sensor array, for receive the electric signal of sensor array output and export with The related energy signal of the energy of the high-energy photon, wherein, the reading circuit includes the shared reading circuit of the first number With the independent reading circuit of the second number, first number is equal to the number of the sensor in each sensor region, and The shared reading circuit of first number connects all the sensors in each sensor region correspondingly, each independent Reading circuit connects the single sensor in the sensor array;
Process circuit, the energy signal of the first number that the shared reading circuit for receiving first number exports respectively with And the energy signal of the second number that the independent reading circuit of second number exports respectively, and based on first number The regularity of energy distribution of energy signal and the energy signal of second number determines the reaction projected position of high-energy photon, its In, the reaction projected position is response location of the high-energy photon in the scintillation crystal on the sensor array Projection.
9. photon detection apparatus according to claim 8, it is characterised in that the process circuit is additionally operable to based on described the The energy signal of one number determines energy and/or the arrival time of the high-energy photon.
10. photon detection apparatus according to claim 8, it is characterised in that first number is not less than the high energy Photon and the scintillation crystal react caused by greatest irradiation scope institute of the scintillation photons in the sensor array Comprising sensor number.
11. photon detection apparatus according to claim 8, it is characterised in that the coordinate in different sensors region is consistent Sensor share same shared reading circuit.
12. photon detection apparatus according to claim 8, it is characterised in that second number is not less than N-1, wherein, N is the number of at least two sensor region.
13. a kind of photon detection device, including:
Receiving module, the first number that the shared reading circuit for receiving the first number being connected with sensor array exports respectively The energy signal for the second number that purpose energy signal and the independent reading circuit of the second number export respectively, wherein, it is described Sensor array column average is divided at least two sensor regions, and first number is equal to the sensing in each sensor region The number of device, and the shared reading circuit of first number connects all biographies in each sensor region correspondingly Sensor, each individually reading circuit connect the single sensor in the sensor array;And
Position determination module, for the energy signal based on first number and the energy of the energy signal of second number The regularity of distribution determines the reaction projected position of high-energy photon, wherein, it is described reaction projected position for the high-energy photon with institute State projection of the response location in the scintillation crystal of sensor array coupling on the sensor array.
14. a kind of photon detection system, including processor and memory, wherein, it is stored with computer program in the memory Instruction, it is used to perform following steps when the computer program instructions are run by the processor:
The energy signal for the first number that the shared reading circuit for the first number that reception is connected with sensor array exports respectively And second number the energy signal of the second number that exports respectively of independent reading circuit, wherein, the sensor array is put down At least two sensor regions are divided into, first number is equal to the number of the sensor in each sensor region, and And the shared reading circuit of first number connects all the sensors in each sensor region, Mei Gedan correspondingly Only reading circuit connects the single sensor in the sensor array;And
The regularity of energy distribution of the energy signal of energy signal and second number based on first number determines high energy The reaction projected position of photon, wherein, the reaction projected position is that the high-energy photon couples with the sensor array Scintillation crystal in projection of the response location on the sensor array.
15. a kind of storage medium, stores programmed instruction on said storage, described program instruction is operationally used to hold Row following steps:
The energy signal for the first number that the shared reading circuit for the first number that reception is connected with sensor array exports respectively And second number the energy signal of the second number that exports respectively of independent reading circuit, wherein, the sensor array is put down At least two sensor regions are divided into, first number is equal to the number of the sensor in each sensor region, and And the shared reading circuit of first number connects all the sensors in each sensor region, Mei Gedan correspondingly Only reading circuit connects the single sensor in the sensor array;And
The regularity of energy distribution of the energy signal of energy signal and second number based on first number determines high energy The reaction projected position of photon, wherein, the reaction projected position is that the high-energy photon couples with the sensor array Scintillation crystal in projection of the response location on the sensor array.
CN201710975365.XA 2017-10-16 2017-10-16 Photon detection method, device, equipment and system and storage medium Active CN107874773B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201710975365.XA CN107874773B (en) 2017-10-16 2017-10-16 Photon detection method, device, equipment and system and storage medium
PCT/CN2018/110255 WO2019076268A1 (en) 2017-10-16 2018-10-15 Photon detection method, device, apparatus and system, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710975365.XA CN107874773B (en) 2017-10-16 2017-10-16 Photon detection method, device, equipment and system and storage medium

Publications (2)

Publication Number Publication Date
CN107874773A true CN107874773A (en) 2018-04-06
CN107874773B CN107874773B (en) 2020-12-08

Family

ID=61781838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710975365.XA Active CN107874773B (en) 2017-10-16 2017-10-16 Photon detection method, device, equipment and system and storage medium

Country Status (2)

Country Link
CN (1) CN107874773B (en)
WO (1) WO2019076268A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019076268A1 (en) * 2017-10-16 2019-04-25 中派科技(深圳)有限责任公司 Photon detection method, device, apparatus and system, and storage medium
CN110426730A (en) * 2019-07-19 2019-11-08 东软医疗系统股份有限公司 Signal read circuit, signal reading method and device
CN114186166A (en) * 2021-11-12 2022-03-15 苏州瑞派宁科技有限公司 Method and device for calculating reaction depth and computer readable storage medium

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101163988A (en) * 2005-04-22 2008-04-16 皇家飞利浦电子股份有限公司 Digital silicon photomultiplier for tof-pet
WO2012056036A2 (en) * 2010-10-28 2012-05-03 Sensl Technologies Ltd Imaging device and method
CN102944890A (en) * 2012-11-06 2013-02-27 中国科学院高能物理研究所 PS-PMT (position sensitive-photomultiplier tube) based detector signal readout method and system
US20140246594A1 (en) * 2011-05-10 2014-09-04 Eberhard-Karls-Universitaet Tuebingen Universitaetsklinikum Gamma detector based on geigermode avalanche photodiodes
EP2876464A1 (en) * 2012-07-18 2015-05-27 Kabushiki Kaisha Toshiba, Inc. Data acquisition device, pair annihilation gamma ray detector, and pair annihilation gamma ray detection method
US20150285922A1 (en) * 2014-04-03 2015-10-08 Siemens Medical Solutions Usa, Inc. Silicon Photomultiplier Based TOF-PET Detector
CN105093258A (en) * 2015-09-30 2015-11-25 武汉中派科技有限责任公司 Photon measurement front-end circuit
CN105182396A (en) * 2015-06-29 2015-12-23 苏州瑞派宁科技有限公司 Detector signal readout channel multiplexing method
US20160170045A1 (en) * 2014-12-15 2016-06-16 General Electric Company Systems and methods for positron emission tomography signal isolation
CN105816194A (en) * 2015-01-07 2016-08-03 苏州瑞派宁科技有限公司 Detector signal readout channel multiplexing method
CN106461793A (en) * 2015-03-17 2017-02-22 皇家飞利浦有限公司 Scintillation event position determination in a radiation particle detector
CN106443747A (en) * 2016-11-15 2017-02-22 武汉中派科技有限责任公司 Method and device for measuring arrival time of high-energy photons
CN106556857A (en) * 2016-10-19 2017-04-05 武汉中派科技有限责任公司 Response location localization method and device
CN106656390A (en) * 2016-11-15 2017-05-10 武汉中派科技有限责任公司 Apparatus and method for measuring photon time information
CN107076859A (en) * 2014-10-17 2017-08-18 皇家飞利浦有限公司 With shared and depth of interaction estimation the pet detector scintillator arrangement of light
CN107121692A (en) * 2017-06-05 2017-09-01 中派科技(深圳)有限责任公司 Detector and the transmitting imaging device with the detector

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015022354A1 (en) * 2013-08-14 2015-02-19 Koninklijke Philips N.V. Pixel identification for small pitch scintillation crystal arrays
CN107874773B (en) * 2017-10-16 2020-12-08 中派科技(深圳)有限责任公司 Photon detection method, device, equipment and system and storage medium

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101163988A (en) * 2005-04-22 2008-04-16 皇家飞利浦电子股份有限公司 Digital silicon photomultiplier for tof-pet
WO2012056036A2 (en) * 2010-10-28 2012-05-03 Sensl Technologies Ltd Imaging device and method
US20140246594A1 (en) * 2011-05-10 2014-09-04 Eberhard-Karls-Universitaet Tuebingen Universitaetsklinikum Gamma detector based on geigermode avalanche photodiodes
EP2876464A1 (en) * 2012-07-18 2015-05-27 Kabushiki Kaisha Toshiba, Inc. Data acquisition device, pair annihilation gamma ray detector, and pair annihilation gamma ray detection method
CN102944890A (en) * 2012-11-06 2013-02-27 中国科学院高能物理研究所 PS-PMT (position sensitive-photomultiplier tube) based detector signal readout method and system
US20150285922A1 (en) * 2014-04-03 2015-10-08 Siemens Medical Solutions Usa, Inc. Silicon Photomultiplier Based TOF-PET Detector
CN107076859A (en) * 2014-10-17 2017-08-18 皇家飞利浦有限公司 With shared and depth of interaction estimation the pet detector scintillator arrangement of light
US20160170045A1 (en) * 2014-12-15 2016-06-16 General Electric Company Systems and methods for positron emission tomography signal isolation
CN105816194A (en) * 2015-01-07 2016-08-03 苏州瑞派宁科技有限公司 Detector signal readout channel multiplexing method
CN106461793A (en) * 2015-03-17 2017-02-22 皇家飞利浦有限公司 Scintillation event position determination in a radiation particle detector
CN105182396A (en) * 2015-06-29 2015-12-23 苏州瑞派宁科技有限公司 Detector signal readout channel multiplexing method
CN105093258A (en) * 2015-09-30 2015-11-25 武汉中派科技有限责任公司 Photon measurement front-end circuit
CN106556857A (en) * 2016-10-19 2017-04-05 武汉中派科技有限责任公司 Response location localization method and device
CN106443747A (en) * 2016-11-15 2017-02-22 武汉中派科技有限责任公司 Method and device for measuring arrival time of high-energy photons
CN106656390A (en) * 2016-11-15 2017-05-10 武汉中派科技有限责任公司 Apparatus and method for measuring photon time information
CN107121692A (en) * 2017-06-05 2017-09-01 中派科技(深圳)有限责任公司 Detector and the transmitting imaging device with the detector

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019076268A1 (en) * 2017-10-16 2019-04-25 中派科技(深圳)有限责任公司 Photon detection method, device, apparatus and system, and storage medium
CN110426730A (en) * 2019-07-19 2019-11-08 东软医疗系统股份有限公司 Signal read circuit, signal reading method and device
CN110426730B (en) * 2019-07-19 2022-11-01 东软医疗系统股份有限公司 Signal reading circuit, signal reading method and device
CN114186166A (en) * 2021-11-12 2022-03-15 苏州瑞派宁科技有限公司 Method and device for calculating reaction depth and computer readable storage medium
WO2023082307A1 (en) * 2021-11-12 2023-05-19 苏州瑞派宁科技有限公司 Method and apparatus for calculating depth of interaction, and computer-readable storage medium

Also Published As

Publication number Publication date
CN107874773B (en) 2020-12-08
WO2019076268A1 (en) 2019-04-25

Similar Documents

Publication Publication Date Title
CN105452902B (en) For the pixel identification of small spacing scintillation crystal array
CN104076383B (en) Virtual pet detector and the quasi-pixelated readout scheme for PET
US5600144A (en) Three dimensional imaging detector employing wavelength-shifting optical fibers
US7786443B2 (en) Dual-modality imaging
CN102007430B (en) Radiation tomographic equipment
US8466418B2 (en) Gamma ray detector, radiation diagnostic device, tomography device, and method of analyzing tomography device
EP2138867A1 (en) Radioactive ray detecting method and device utilizing energy and position information
ATE128557T1 (en) HIGH RESOLUTION GAMMA RADIATION DETECTOR FOR POSITRONEMISSION TOMOGRAPHY (PET) AND SINGLE PHOTONEMISSION COMPUTER TOMOGRAPHY (SPECT).
CN107874773A (en) Photon detection method, apparatus, equipment and system and storage medium
CN107735694A (en) Apparatus and method for the interaction depth position emissron tomography detector using two points of sensings
CN101925834A (en) Radiation detector and tomographic equipment provided with the same
CN109683188A (en) A kind of Qi Liankefu event and gamma events meet imaging device and method
WO2010074400A2 (en) Multi-planar detector for measuring gamma-ray image and method for detecting 3d position
CN109581460A (en) Composite detecting device
TWI356689B (en) Method for determining location of gamma interacti
JP2000056023A (en) Pet detector capable of providing depth directional information
CN103559391B (en) Judge random coincidence event and judge the method and device of true coincidence event
CN106461792A (en) Scintillation event localization in a radiation particle detector
EP0829022B1 (en) An apparatus for scintigraphic analysis, particularly a mammograph, with sub-millimetric spatial resolution
CN108226988A (en) Radioactive ray method for detecting position, radioactive rays position detector and PET device
JP4843346B2 (en) Mammography equipment
CN209400708U (en) Composite detecting device and system
CN107456235A (en) Method for correcting position and system
Szlávecz et al. A novel model-based PET detector block simulation approach
CN219126405U (en) Crystal array detector and emission imaging device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant