CN117495990A - PET sinogram data compression storage method, system and equipment based on monolithic FPGA - Google Patents
PET sinogram data compression storage method, system and equipment based on monolithic FPGA Download PDFInfo
- Publication number
- CN117495990A CN117495990A CN202410001693.XA CN202410001693A CN117495990A CN 117495990 A CN117495990 A CN 117495990A CN 202410001693 A CN202410001693 A CN 202410001693A CN 117495990 A CN117495990 A CN 117495990A
- Authority
- CN
- China
- Prior art keywords
- data
- response line
- line data
- sinogram
- fpga
- 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.)
- Pending
Links
- 238000003860 storage Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000013144 data compression Methods 0.000 title claims description 22
- 230000004044 response Effects 0.000 claims abstract description 87
- 238000012636 positron electron tomography Methods 0.000 claims abstract description 43
- 238000005259 measurement Methods 0.000 claims abstract description 23
- 238000012879 PET imaging Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000007906 compression Methods 0.000 claims abstract description 10
- 230000006835 compression Effects 0.000 claims abstract description 10
- 238000005520 cutting process Methods 0.000 claims abstract description 6
- 238000000691 measurement method Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 abstract description 11
- 230000001133 acceleration Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000001959 radiotherapy Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- -1 carrier Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012831 peritoneal equilibrium test Methods 0.000 description 1
- 238000012877 positron emission topography Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007858 starting material Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Abstract
The invention discloses a PET sinogram compression storage method, system and equipment based on a single-chip FPGA, belonging to the field of digital data processing of medical images. The method can compress, store and output the sinogram data acquired in the PET imaging process in real time in the FPGA, and comprises the following steps: after the PET detector and the electronic system measure the binary data of the distance and the angle of the corresponding response line of each coincidence event in real time, the corresponding compression address and the corresponding counting label are respectively obtained by a bit cutting and splicing mode, the compression address and the corresponding multiple counts are counted and stored in the FPGA, new data obtained by subsequent measurement are inserted or accumulated on the existing data according to the compression address and by a traversing comparison mode, and sine graph data is finally output according to the coordinate sequence after the number of the required instance of statistical imaging. The invention can obviously reduce the storage space required by storing the sinogram in the FPGA, realize the real-time statistics and storage of the sinogram in the FPGA and improve the reconstruction speed of the PET image.
Description
Technical Field
The invention relates to the field of digital data processing of medical images, in particular to a PET sinogram data compression storage method, system and equipment based on a single-chip FPGA.
Background
Positron emission tomography (Positron Emission Tomography, PET) is an advanced nuclear medicine imaging modality whose basic imaging principle is to inject a radionuclide-labeled tracer into a living body, decay the radionuclide to produce positrons, annihilate the negative electrons of the organism tissue, and emit a pair of annihilation photons in opposite directions, and a detector, in combination with an electronics system, detects a pair of photons generated by the same annihilation event by conforming to a measurement technique, and localizes the location of annihilation by a Line of the photon locations. Sinograms (sinograms) are two-dimensional statistical histograms with projection coordinate data of a large number of response lines as storage addresses and counts as storage data, and final PET images are obtained by data correction and image reconstruction of the sinograms.
The electronics system in a conventional PET system typically outputs raw position information corresponding to the event and stores and calculates the position information in a computer system to obtain sinogram data. In an application scenario requiring fast PET imaging, for example, in a scenario of conducting real-time guidance on radiotherapy by using PET imaging information in image-guided radiotherapy, a process of acquiring a sinogram by using an FPGA in an electronic system needs to be accelerated, which includes a storage and output process of the sinogram.
For the current high-precision PET system, the addresses of the sinograms are finely divided, which can make the storage space required for completely storing the sinograms greater than or equal to 1Mb, for example, the size of the sinograms is 512×256, and the length of each count storage is 8 bits. While current PET systems typically use more detector rings to construct longer axial detector systems in order to achieve higher sensitivity, this can result in PET systems having tens or even hundreds of imaging slices, i.e., the sinograms corresponding to these imaging slices need to be stored in the FPGA at the same time. The RAM storage resources in the FPGA are very limited, for example, the maximum RAM storage resources that can be provided by the Xilinx Spartan-7, artix-7 and kinex-7 series FPGAs are 4.32, 13.14 and 34.38Mb respectively, so that it is difficult to implement the complete storage of the sinogram in a monolithic FPGA in a manner of taking the projection coordinate data as the storage address and the count as the storage data, and further, the acceleration of PET imaging cannot be implemented based on the monolithic FPGA, or the system integration is poor and the cost is high due to the use of multiple FPGAs.
In view of this, the present invention has been made.
Disclosure of Invention
The invention aims to provide a PET sinogram data compression storage method, a system and equipment of a monolithic FPGA, which can compress and store sinogram data in a PET rapid imaging application scene, greatly reduce RAM storage resources required by sinogram storage, store all sinograms of the system in the monolithic FPGA, realize acceleration of PET imaging process, improve system integration level and reduce cost, and further solve the technical problems in the prior art.
The invention aims at realizing the following technical scheme:
a PET sinogram data compression storage method of a monolithic FPGA comprises the following steps:
step 1: real-time measurement system composed of PET detector and electronic system is used to obtain the corresponding response line data of each annihilation event by real-time measurement according to the measurement method, and each response line data is represented by its projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
step 2: obtaining corresponding compressed addresses of each response line data obtained in the step 1 in a bit cutting mode, and marking the corresponding compressed addresses asMeanwhile, a counting label corresponding to each response line data is obtained in a bit splicing mode and is marked as N;
step 3: for each response line data processed in step 2Is stored in the data format of (1) and is compressed according to the compression address corresponding to each response line data>After the sizes are sequenced, the data are inserted into the RAM of the monolithic FPGA as new response line data or accumulated on stored response line data, so that statistics and storage of the response line data are realized; in the data format of the response line data store +.>A count index representing corresponding data in the sinogram;
step 4: judging whether the number of cases of response line data stored in the RAM of the monolithic FPGA meets the PET imaging requirement, if not, repeating the steps 1 to 3, and if so, executing the step 5;
step 5: and outputting the response line data stored in the RAM of the FPGA according to the coordinate sequence of the sinogram.
A PET sinogram data compression storage system of a monolithic FPGA is used for realizing the method of the invention, which comprises the following steps:
the system comprises a real-time measurement system and a sinogram data processing module comprising a single-chip FPGA; wherein,
the real-time measurement system consists of a PET detector and an electronic system, and can obtain response line data corresponding to each annihilation event by real-time measurement according to a measurement method, wherein each response line data is represented by projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
the sinogram data processing module is used for receiving response line data output by the real-time measurement system, and obtaining a corresponding compressed address in an FPGA (field programmable gate array) in a bit truncation manner for each received response line dataMeanwhile, the counting label N of the response line data is obtained in a bit splicing mode; each response line is then dataData format of->Count marks representing corresponding data in the sinogram, and compressed addresses corresponding to the respective response line data +.>After the sizes are sequenced, the data are inserted into the RAM of the monolithic FPGA as new response line data or accumulated on stored response line data, so that statistics and storage of the response line data are realized; and repeating the processing until the number of the cases of the response line data stored in the RAM of the single-chip FPGA meets the requirement of PET imaging, and outputting the response line data stored in the RAM of the single-chip FPGA according to the coordinate sequence of the sinogram.
A processing apparatus, comprising:
at least one memory for storing one or more programs;
at least one processor capable of executing one or more programs stored in the memory, which when executed by the processor, enable the processor to implement the methods of the present invention.
Compared with the prior art, the PET sinogram data compression storage method, system and equipment for the monolithic FPGA have the beneficial effects that:
by utilizing the continuity and sparsity of the sinogram in the PET rapid imaging application scene, address data of a plurality of adjacent counts are compressed and stored in a mode of combining the plurality of counts according to compressed addresses, so that the storage of invalid count data is greatly avoided, compared with the mode of completely storing the sinogram, the storage resource of a RAM is greatly saved, the fact that all sinograms of a system are stored in a single-chip FPGA is realized, the acceleration of the PET imaging process by the single-chip FPGA is further realized, and meanwhile, the system integration level is improved and the cost is reduced.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic diagram of a small animal PET system imaging principle according to an embodiment of the present invention.
Fig. 2 is a flowchart of a method for storing PET sinogram data based on a monolithic FPGA according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a sinogram data compression storage method according to an embodiment of the present invention.
FIG. 4 is a schematic diagram of test results provided by an embodiment of the present invention, wherein (a) is the maximum number of instances received by each of the 39 sinograms; (b) The maximum number of compression cases stored for each sinogram RAM after data compression.
Detailed Description
The technical scheme in the embodiment of the invention is clearly and completely described below in combination with the specific content of the invention; it will be apparent that the described embodiments are only some embodiments of the invention, but not all embodiments, which do not constitute limitations of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to fall within the scope of the invention.
The terms that may be used herein will first be described as follows:
the term "and/or" is intended to mean that either or both may be implemented, e.g., X and/or Y are intended to include both the cases of "X" or "Y" and the cases of "X and Y".
The terms "comprises," "comprising," "includes," "including," "has," "having" or other similar referents are to be construed to cover a non-exclusive inclusion. For example: including a particular feature (e.g., a starting material, component, ingredient, carrier, formulation, material, dimension, part, means, mechanism, apparatus, step, procedure, method, reaction condition, processing condition, parameter, algorithm, signal, data, product or article of manufacture, etc.), should be construed as including not only a particular feature but also other features known in the art that are not explicitly recited.
The term "consisting of … …" is meant to exclude any technical feature element not explicitly listed. If such term is used in a claim, the term will cause the claim to be closed, such that it does not include technical features other than those specifically listed, except for conventional impurities associated therewith. If the term is intended to appear in only a clause of a claim, it is intended to limit only the elements explicitly recited in that clause, and the elements recited in other clauses are not excluded from the overall claim.
Unless specifically stated or limited otherwise, the terms "mounted," "connected," "secured," and the like should be construed broadly to include, for example: the connecting device can be fixedly connected, detachably connected or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the terms herein above will be understood by those of ordinary skill in the art as the case may be.
The terms "center," "longitudinal," "transverse," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counterclockwise," etc. refer to an orientation or positional relationship based on that shown in the drawings, merely for ease of description and to simplify the description, and do not explicitly or implicitly indicate that the apparatus or element in question must have a particular orientation, be constructed and operated in a particular orientation, and therefore should not be construed as limiting the present disclosure.
The PET sinogram data compression storage method and system based on the monolithic FPGA provided by the invention are described in detail below. What is not described in detail in the embodiments of the present invention belongs to the prior art known to those skilled in the art. The specific conditions are not noted in the examples of the present invention and are carried out according to the conditions conventional in the art or suggested by the manufacturer. The reagents or apparatus used in the examples of the present invention were conventional products commercially available without the manufacturer's knowledge.
As shown in fig. 2, an embodiment of the present invention provides a PET sinogram data compression storage method based on a monolithic FPGA, including the following steps:
step 1: real-time measurement system composed of PET detector and electronic system is used to obtain the corresponding response line data of each annihilation event by real-time measurement according to the measurement method, and each response line data is represented by its projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
step 2: step by step1 to obtain corresponding compressed address by bit cutting mode, which is marked asMeanwhile, a counting label corresponding to each response line data is obtained in a bit splicing mode and is marked as N;
for a sinogram of size A B, the projection space coordinates of the original LOR data are inThe corresponding value range is 0-A-1, and the bit number required by binary storage in the FPGA is +.>，/>The corresponding value range is 0-B-1, and the number of bits required by binary number storage in the FPGA is +.>Adjacent 2 in the sinogram m ×2 n Compressing and storing data of the counting addresses, wherein m and n are non-negative integers and are not 0 at the same time, and compressing address +_ corresponding to each response line data obtained in a bit cutting mode under a binary data format>The method comprises the following steps:
；
；
meanwhile, the counting index N can be calculated by the following formula:
；
the counting label N of each response line data can be obtained directly by a bit-wise splicing mode under a binary data format, and is expressed as follows:
。
the data format stored in the RAM of the FPGA can be expressed asWhereinCount numbers for corresponding data in the sinogram, a total of 2 m ×2 n A plurality of;
step 3: counting and storing LOR data in RAM of FPGA, and storing the dataOrdering according to the size. For any LOR data obtained by real-time measurement, performing traversal comparison between a compressed address serving as new LOR data and a compressed address of stored LOR data in RAM, determining the RAM address where the new LOR data is stored, if the compressed address is stored in RAM, according to the counting label obtained by calculation in the step 2, and comparing the counting label corresponding to the stored LOR data with the counting label corresponding to the stored LOR data>Adding 1; if the compressed address is not stored in RAM, the new LOR data is compressed according to the corresponding compressed address +.>Is inserted into RAM and the calculation label corresponding to the LOR data is +.>Marked 1, others->And is noted as 0.
Step 4: judging whether the instance number of LOR data stored in the RAM of the monolithic FPGA meets the PET imaging requirement, if not, repeating the steps 1 to 3, and if so, executing the step 5;
step 5: the LOR data stored in the RAM of the monolithic FPGA are output according to the coordinate sequence of the sinogram, and the following 3 modes are adopted according to different specific application occasions:
mode one: directly according to storageThe data format is read out sequentially from front to back according to the stored RAM address;
mode two: restoring the compressed address and the count label into original LOR coordinate data, combining the count obtained by statistics under the coordinate, and according to the following stepsIs read out according to the data format of the data;
mode three: and outputting the complete sinogram count according to the coordinate sequence without outputting LOR coordinate data, namely outputting A multiplied by B count data in sequence, wherein the count output which is not counted is 0.
The embodiment of the invention also provides a PET sinogram data compression storage system of the monolithic FPGA, which is used for realizing the method, and comprises the following steps:
the system comprises a real-time measurement system and a sinogram data processing module comprising a single-chip FPGA; wherein,
the real-time measurement system consists of a PET detector and an electronic system, and can obtain response line data corresponding to each annihilation event by real-time measurement according to a measurement method, wherein each response line data is represented by projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
the sinogram data processing module is used forReceiving response line data output by a real-time measurement system, and obtaining a corresponding compressed address in an FPGA (field programmable gate array) in a bit truncation manner for each received response line dataMeanwhile, the counting label N of the response line data is obtained in a bit splicing mode; each response line is then dataData format of->Count marks representing corresponding data in the sinogram, and compressed addresses corresponding to the respective response line data +.>After the sizes are sequenced, the data are inserted into the RAM of the monolithic FPGA as new response line data or accumulated on stored response line data, so that statistics and storage of the response line data are realized; and repeating the processing until the number of the cases of the response line data stored in the RAM of the single-chip FPGA meets the requirement of PET imaging, and outputting the response line data stored in the RAM of the single-chip FPGA according to the coordinate sequence of the sinogram.
The embodiment of the invention further provides a processing device, which comprises:
at least one memory for storing one or more programs;
at least one processor capable of executing one or more programs stored in the memory, which when executed by the processor, enable the processor to implement the methods described above.
In order to clearly show the technical scheme and the technical effects provided by the invention, the PET sinogram data compression storage method and system based on the monolithic FPGA provided by the embodiment of the invention are described in detail in the following by using a specific embodiment.
Example 1
The complete exemplary implementation of the invention is accomplished based on a small animal PET prototype electronics system autonomously designed by a nuclear detection and nuclear power national emphasis laboratory, which constitutes a complete small animal PET system with a PET detector and an image reconstruction system. In the imaging process of the PET system, as shown in fig. 1, the PET detector is composed of 20 detector sub-rings, 20×20 sinograms are generated by adopting a 3D data acquisition mode, in order to realize rapid imaging, the system stores acquired data in 39 2D sinograms by adopting a data recombination mode, and finally, image reconstruction is realized through a 2D filtered back projection (Filtered Back Projection, FBP) algorithm. In the PET imaging process, the FPGA is used for accelerating the sinogram acquisition process, and the storage resource limitation of the FPGA and the subsequent expandability of the system are considered, so that the complete storage of all sinograms is difficult to realize in a single-chip FPGA.
The embodiment provides a single-chip FPGA-based PET sinogram data compression and storage method, which is used for compressing and storing the sinogram data, and the flow of the method is shown in a figure 2, and the method comprises the following steps:
step 1: real-time measurement system composed of PET detector and electronic system, and each LOR data corresponding to each annihilation event is obtained by real-time measurement according to the measurement method, and each LOR data is represented by its projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
in this embodiment, the distance coordinate of the sinogram is divided into 512 steps, the angle coordinate is divided into 256 steps, the total address number of each sinogram is 131072, the depth of the storage count under each address is 8 bits, and the total required storage space is the total required storage space when all 39 sinograms are completely stored in the FPGAThe method comprises the following steps:
。
step 2: the compressed address is recorded asFor a sinogram of size A B +.>The corresponding value range is 0-A-1, and the bit number required by binary storage in the FPGA is +.>，/>The corresponding value range is 0-B-1, and the number of bits required by binary number storage in the FPGA is +.>Adjacent 2 in the sinogram m ×2 n The counted address data are compressed and stored, wherein m and n are non-negative integers and are not 0 at the same time, and the compressed address stored under the binary data format obtained in a bit truncation mode is:
；
；
meanwhile, the counting label can be directly calculated by the following formula:
，
and can be directly obtained by a bit-wise splicing mode under a binary data format, which is expressed as follows:
。
the data format stored in the FPGA can be expressed asWherein->Count numbers for corresponding data in the sinogram, a total of 2 m ×2 n A plurality of;
in this embodiment, the size of the sinogram is 512×256, address data of adjacent 1×4 counts in the sinogram is compressed and stored, that is, m=0 and n=2 are taken, and the calculation of the compressed address can be expressed as:
；
；
the calculation of the count index may be expressed as:
；
the data format stored in the RAM of the FPGA can be expressed asA schematic diagram of the sinogram data compression process is shown in fig. 3.
Step 3: in RAM of FPGA, LOR data is counted and stored, and the stored data is according to the following stepsOrdering the sizes of (3); for any LOR data obtained by real-time measurement, performing traversal comparison on the compressed address of the new data and the compressed address of the data stored in the RAM, determining the RAM address where the new data is stored, and if the compressed address is stored in the RAM, calculating the counting mark according to the step 2Number +.>Adding 1; if the compressed address is not stored in RAM, new data is written as +.>Is inserted into RAM and the LOR data therein is corresponding +.>Marked 1, others->And is noted as 0.
Step 4: judging whether the LOR instance number stored in the RAM of the monolithic FPGA meets the PET imaging requirement, if not, repeating the steps 1 to 3, and if so, executing the step 5.
Step 5: the stored data are output according to the coordinate sequence of the sinogram, and the following 3 types are different according to specific application occasions:
mode one: directly according to storageThe data format is read out sequentially from front to back according to the stored RAM address;
mode two: restoring the compressed address and the count label into original LOR coordinate data, combining the count obtained by statistics under the coordinate, and according to the following stepsIs read out according to the data format of the data;
mode three: and outputting the complete sinogram count according to the coordinate sequence without outputting LOR coordinate data, namely outputting A multiplied by B count data in sequence, wherein the count output which is not counted is 0.
In this embodiment, in order to reduce the pressure stored and calculated by the back-end computer, the data is outputted by using the output method 2.
Simulation data generated using a GATE Monte Carlo simulation platform is used to advance the data compression method used in the embodimentsThe line test, simulation time was set to 50ms. Fig. 4 (a) is a maximum number of cases received by each of the 39 sinograms, the total number of cases is 157541, wherein the maximum number of cases of a single sinogram is 8049, which is far smaller than the total number of addresses of the sinograms, so that the data of the sinograms has sparseness; in FIG. 4 (b), the maximum number of compression cases is stored in each sinogram RAM after data compression, the total number of compression cases is 41955, and the maximum count is 15, so that the data is compressed in binaryThe storage depth of the compressed data is set to be 4 bits, the total data length of the compressed data is 31 bits, and the size of the storage space required after compression can be calculated out +.>The method comprises the following steps:
；
by adopting the sinogram data compression method in the embodiment, all sinogram data can be stored in the monolithic FPGA, and the data compression rate R can be calculated as follows:
。
in summary, the embodiment of the invention utilizes the continuity and sparsity of the sinogram in the PET rapid imaging application scene, compresses address data of a plurality of adjacent counts and stores the address data in a mode of combining the compressed addresses with the plurality of counts, thereby greatly avoiding the storage of invalid count data, greatly saving RAM storage resources compared with a mode of completely storing the sinogram, realizing the storage of all sinograms of the system in a single-chip FPGA, further realizing the acceleration of the PET imaging process by the single-chip FPGA, and improving the integration level of the system and reducing the cost.
Those of ordinary skill in the art will appreciate that: all or part of the flow of the method implementing the above embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and the program may include the flow of the embodiment of each method as described above when executed. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), or the like.
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions easily contemplated by those skilled in the art within the scope of the present invention should be included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims. The information disclosed in the background section herein is only for enhancement of understanding of the general background of the invention and is not to be taken as an admission or any form of suggestion that this information forms the prior art already known to those of ordinary skill in the art.
Claims (6)
1. A PET sinogram data compression storage method of a monolithic FPGA is characterized by comprising the following steps:
step 1: real-time measurement system composed of PET detector and electronic system is used to obtain the corresponding response line data of each annihilation event by real-time measurement according to the measurement method, and each response line data is represented by its projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
step 2: obtaining corresponding compressed addresses of each response line data obtained in the step 1 in a bit cutting mode, and marking the corresponding compressed addresses asMeanwhile, a counting label corresponding to each response line data is obtained in a bit splicing mode and is marked as N;
step 3: for each response line data processed in step 2Is stored in the data format of (1) and is compressed according to the compression address corresponding to each response line data>After the sizes are sequenced, the data are inserted into the RAM of the monolithic FPGA as new response line data or accumulated on stored response line data, so that statistics and storage of the response line data are realized; in the data format of the response line data store +.>A count index representing corresponding data in the sinogram;
step 4: judging whether the number of cases of response line data stored in the RAM of the monolithic FPGA meets the PET imaging requirement, if not, repeating the steps 1 to 3, and if so, executing the step 5;
step 5: and outputting the response line data stored in the RAM of the FPGA according to the coordinate sequence of the sinogram.
2. The method for storing PET sinogram data compression of a monolithic FPGA of claim 1, wherein, in the method, if the sinogram to be PET imaged is a x B in size, in the projection space coordinates of the original response line data,the corresponding value range is 0-A-1, and the bit number required by binary storage in the FPGA is +.>，/>The corresponding value range is 0-B-1, and the number of bits required by binary number storage in the FPGA is +.>Adjacent 2 in the sinogram m ×2 n Compressing and storing data of the counting addresses, wherein m and n are non-negative integers and are not 0 at the same time, and each response line data obtained in a bit cutting mode corresponds to a compressed address under a binary data formatThe method comprises the following steps:
；
；
meanwhile, the counting label N of each response line data obtained in a bit splicing mode under a binary data format is as follows:
。
3. the method for storing PET sinogram data of a monolithic FPGA according to claim 1 or 2, wherein in the step 3, any response line data measured in real time is inserted as new response line data in a RAM of the monolithic FPGA or accumulated on stored response line data in the following manner, comprising:
performing traversal comparison on the compressed address of the new response line data and the compressed address of the stored response line data in the RAM of the monolithic FPGA, determining the RAM address where the new response line data is stored, and if the compressed address is stored in the RAM, sounding the stored RAM according to the counting label calculated in the step 2Counting label of line-applied dataAdding 1; if the compressed address is not stored in RAM, the new response line data is compressed according to the corresponding compressed address +.>Is inserted into the RAM, and the corresponding count index +.>Marked 1, others->And is noted as 0.
4. The method for storing PET sinogram data of a monolithic FPGA according to claim 1 or 2, wherein in step 5, the outputting the response line data stored in the RAM of the FPGA according to the coordinate sequence of the sinogram includes:
mode one: directly stored asThe data format is read out sequentially from front to back according to the stored RAM address;
mode two: restoring the compressed address and the count label into original response line data coordinates, combining the count counted under the coordinates, and according to the following stepsIs read out according to the data format of the data;
mode three: and outputting complete sinogram count according to the coordinate sequence without outputting the coordinates of the response line data, namely outputting A multiplied by B count data in sequence, wherein the count output without statistics is 0.
5. A monolithic FPGA-based PET sinogram data compression storage system for implementing the method of any one of claims 1-4, comprising:
the system comprises a real-time measurement system and a sinogram data processing module comprising a single-chip FPGA; wherein,
the real-time measurement system consists of a PET detector and an electronic system, and can obtain response line data corresponding to each annihilation event by real-time measurement according to a measurement method, wherein each response line data is represented by projection space coordinates and is recorded asWherein->For the length of the perpendicular from the origin of coordinates to the line of response, +.>An included angle between the vertical line and the x axis;
the sinogram data processing module is used for receiving response line data output by the real-time measurement system, and obtaining a corresponding compressed address in an FPGA (field programmable gate array) in a bit truncation manner for each received response line dataMeanwhile, the counting label N of the response line data is obtained in a bit splicing mode; each response line is then dataData format of->Count marks representing corresponding data in the sinogram, and compressed addresses corresponding to the respective response line data +.>After the sizes are sequenced, the data are inserted into the RAM of the monolithic FPGA as new response line data or accumulated on stored response line data, so that statistics and storage of the response line data are realized; repeating the above stepsAnd outputting the response line data stored in the RAM of the monolithic FPGA according to the coordinate sequence of the sinogram until the number of instances of the response line data stored in the RAM of the monolithic FPGA meets the requirement of PET imaging.
6. A processing apparatus, comprising:
at least one memory for storing one or more programs;
at least one processor capable of executing one or more programs stored in the memory, which when executed by the processor, cause the processor to implement the method of any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410001693.XA CN117495990A (en) | 2024-01-02 | 2024-01-02 | PET sinogram data compression storage method, system and equipment based on monolithic FPGA |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410001693.XA CN117495990A (en) | 2024-01-02 | 2024-01-02 | PET sinogram data compression storage method, system and equipment based on monolithic FPGA |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117495990A true CN117495990A (en) | 2024-02-02 |
Family
ID=89669450
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410001693.XA Pending CN117495990A (en) | 2024-01-02 | 2024-01-02 | PET sinogram data compression storage method, system and equipment based on monolithic FPGA |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117495990A (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050253075A1 (en) * | 2004-05-14 | 2005-11-17 | Cti Pet Systems, Inc. | Nearest-neighbor rebinning in clinical PET using on-line three dimensional LOR-to-bin mapping |
| US20070269093A1 (en) * | 2006-05-18 | 2007-11-22 | Siemens Medical Solutions Usa, Inc. | On-Line Time-Of-Flight Mashing: DMA Rebinning Applied to a Long-Axis PET-TOF System for Reduced Transaxial Angular Sampling With 3-D Nearest-Neighbor Projection-Space Rebinning in Clinical PET/CT |
| US20070278409A1 (en) * | 2006-06-06 | 2007-12-06 | General Electric Company | Methods and apparatus for pet time of flight |
| US20100074500A1 (en) * | 2008-09-22 | 2010-03-25 | Siemens Medical Solutions Usa, Inc. | System and method for 3d time of flight pet forward projection based on an exact axial inverse rebinning relation in fourier space |
| CN109829869A (en) * | 2019-03-04 | 2019-05-31 | 中山市明峰医疗器械有限公司 | Numerical value domain interpolation method for reducing PET image noise |
| CN111080734A (en) * | 2019-11-25 | 2020-04-28 | 中国科学院深圳先进技术研究院 | Method and terminal for processing Positron Emission Tomography (PET) data |
| US20210074034A1 (en) * | 2019-09-05 | 2021-03-11 | Siemens Medical Solutions Usa, Inc. | Methods and apparatus for neural network based image reconstruction |
| CN112907664A (en) * | 2021-02-08 | 2021-06-04 | 明峰医疗系统股份有限公司 | Image pixel position searching method and system for PET detector based on FPGA module and computer readable storage medium |
| WO2021114306A1 (en) * | 2019-12-14 | 2021-06-17 | 深圳先进技术研究院 | High-precision pet reconstruction method and apparatus based on gpu acceleration |
| US20230386036A1 (en) * | 2022-05-24 | 2023-11-30 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for medical imaging |
-
2024
- 2024-01-02 CN CN202410001693.XA patent/CN117495990A/en active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050253075A1 (en) * | 2004-05-14 | 2005-11-17 | Cti Pet Systems, Inc. | Nearest-neighbor rebinning in clinical PET using on-line three dimensional LOR-to-bin mapping |
| US20070269093A1 (en) * | 2006-05-18 | 2007-11-22 | Siemens Medical Solutions Usa, Inc. | On-Line Time-Of-Flight Mashing: DMA Rebinning Applied to a Long-Axis PET-TOF System for Reduced Transaxial Angular Sampling With 3-D Nearest-Neighbor Projection-Space Rebinning in Clinical PET/CT |
| US20070278409A1 (en) * | 2006-06-06 | 2007-12-06 | General Electric Company | Methods and apparatus for pet time of flight |
| US20100074500A1 (en) * | 2008-09-22 | 2010-03-25 | Siemens Medical Solutions Usa, Inc. | System and method for 3d time of flight pet forward projection based on an exact axial inverse rebinning relation in fourier space |
| CN109829869A (en) * | 2019-03-04 | 2019-05-31 | 中山市明峰医疗器械有限公司 | Numerical value domain interpolation method for reducing PET image noise |
| US20210074034A1 (en) * | 2019-09-05 | 2021-03-11 | Siemens Medical Solutions Usa, Inc. | Methods and apparatus for neural network based image reconstruction |
| CN111080734A (en) * | 2019-11-25 | 2020-04-28 | 中国科学院深圳先进技术研究院 | Method and terminal for processing Positron Emission Tomography (PET) data |
| WO2021114306A1 (en) * | 2019-12-14 | 2021-06-17 | 深圳先进技术研究院 | High-precision pet reconstruction method and apparatus based on gpu acceleration |
| CN112907664A (en) * | 2021-02-08 | 2021-06-04 | 明峰医疗系统股份有限公司 | Image pixel position searching method and system for PET detector based on FPGA module and computer readable storage medium |
| US20230386036A1 (en) * | 2022-05-24 | 2023-11-30 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for medical imaging |
Non-Patent Citations (5)
| Title |
|---|
| ANTONIO J 等: "A PET Design Based on SiPM and Monolithic LYSO Crystals: Performance Evaluation", IEEE TRANSACTIONS ON NUCLEAR SCIENCE, 31 December 2016 (2016-12-31), pages 2471 - 2477 * |
| JIAMING LU等: "Design of the Time Measurement Circuit, Single Data Frame Sorting, and Packaging Logic of a TOF PET System", 2020 6TH INTERNATION COFERENCE ON EVENT-BASED CONTROL, 31 December 2020 (2020-12-31), pages 1 - 7 * |
| ZHONGHUA KUANG等: "Design and performance of SIAT aPET: a uniform high-resolution small animal PET scanner using dual-ended readout detectors", PHYSICS IN MEDICINE & BIOLOGY, 30 March 2020 (2020-03-30), pages 1 - 13 * |
| 林琦琛: "基于FPGA的PET图像快速重建与预处理系统研究", 中国优秀硕士学位论文全文数据库 信息科技辑, 15 July 2021 (2021-07-15), pages 138 - 561 * |
| 郭勇;严伟;黄仰锐;屈军乐;: "FPGA技术在生物医学成像中的研究进展", 生物化学与生物物理进展, no. 06, 28 April 2020 (2020-04-28), pages 483 - 496 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8265365B2 (en) | Time of flight scatter distribution estimation in positron emission tomography | |
| CN109887048B (en) | PET scattering correction method, image reconstruction device and electronic equipment | |
| US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
| Farquhar et al. | ROC and localization ROC analyses of lesion detection in whole-body FDG PET: effects of acquisition mode, attenuation correction and reconstruction algorithm | |
| JP6185262B2 (en) | Nuclear medicine bone image analysis technology | |
| US20100072375A1 (en) | Iterative algorithms for crystal efficiencies estimations from TOF compressed normalization data | |
| US7409077B2 (en) | Nearest-neighbor rebinning in clinical PET using on-line three dimensional LOR-to-bin mapping | |
| US8359345B2 (en) | Iterative algorithms for variance reduction on compressed sinogram random coincidences in PET | |
| EP1631844B1 (en) | Generating detector efficiency estimates for a oet scanner | |
| CN117495990A (en) | PET sinogram data compression storage method, system and equipment based on monolithic FPGA | |
| Yamakawa et al. | Attenuation correction using level set method for application specific PET scanners | |
| Nowak et al. | A multiscale MAP estimation method for Poisson inverse problems | |
| Pal et al. | 2D linear and iterative reconstruction algorithms for a PET-insert scanner | |
| CN113712583B (en) | Dead time correction method, device, computer equipment and storage medium | |
| CN112652029A (en) | PET imaging method, device and equipment | |
| CN112052885B (en) | Image processing method, device, equipment and PET-CT system | |
| KR102283454B1 (en) | image reconstruction method to reconstruct the image by correcting the response depth information included in the observed data using the flight time information of the positron emission tomography | |
| Liu et al. | System response matrix calculation using symmetries for dual‐head PET scanners | |
| Vidal et al. | Artificial evolution for 3D PET reconstruction | |
| US7890282B2 (en) | Estimation of crystal efficiency with axially compressed sinogram | |
| CN108648807B (en) | Image reconstruction method and device | |
| WO2022116143A1 (en) | Pet imaging method, apparatus, and device | |
| Goertzen et al. | A method for measuring the energy spectrum of coincidence events in positron emission tomography | |
| Vidal et al. | PET reconstruction using a cooperative coevolution strategy in LOR space | |
| McLennan et al. | SORTEO: Monte Carlo-based simulator with list-mode capabilities |
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 |