A method and system for determining the result of regional division in the generation of a position mapping table

Abstract

The present invention provides a method and system for determining the result of region division in the generation of a position mapping table. When the boundary division result after the first determination of the bonding point is 64 partitions, and the partition numbers after region sequence coding are within the range of 0 to 63, and the count of pixels corresponding to the crystal bars in each partition is greater than 0, the region division result is qualified; and when the boundary division result after the first determination of the bonding point is 64 partitions, the boundary attribution and region sequence coding can be directly performed. The present invention mainly proofreads the accuracy of region division and sequence numbering when bonding points appear in the scattered points corresponding to adjacent crystal bars in the process of generating the position mapping table of the interaction position of crystal bars and gamma photons, so as to determine the position mapping table of the system.

Description

Technical Field

[0001] The present invention belongs to the technical field of positron emission tomography, and particularly relates to a method and system for determining the result of regional division in the generation of a position mapping table. Background Art

[0002] Positron emission tomography (PET) is a novel three-dimensional nuclear medicine imaging technique, which belongs to a type of emission tomography. It uses positron radionuclides labeled with human metabolites such as glucose as imaging agents, and reflects the metabolic changes by the uptake of the imaging agent by lesions, providing biological metabolic information of diseases. It can diagnose and guide the treatment of various diseases earlier, more sensitively, and accurately, and has important application values in biology, neuroscience, oncology, pharmacokinetics, clinical diagnosis, and treatment evaluation. Its basic principle is to inject a small amount of positron radionuclide tracer into the human body. The positron annihilates with the negative electron in the human body to generate a pair of γ photons emitted back to back. An external detection device is used to detect this pair of photons to detect the distribution of the radioactive drug in the body. Currently, the commonly used detector system is realized by coupling a crystal array composed of crystal bars with photomultiplier tubes. When the γ photons generated by the annihilation of positive and negative electrons enter the crystal bar and react with the crystal, fluorescence will be generated when the crystal atoms de-excite. The fluorescence propagates along the crystal array and is then converted into an electrical signal by the photomultiplier tube. This electrical signal carries the position information of the incident γ photon. The position information of the crystal bar where the incident γ photon reacts is obtained by using a positioning method, so as to determine the position of the coincidence line of the annihilation event. The accuracy of the position information will directly affect the spatial resolution of the detector, and further affect the imaging quality and image quantitative analysis accuracy of PET. Therefore, the accurate positioning of the crystal bar where the γ photon undergoes energy deposition is crucial.

[0003] However, in fact, due to differences in the arrangement of scintillation crystals, the cutting method of the light-splitting light guide, the signal processing ability of photomultiplier tubes, and various reasons such as the front-end electronic signal acquisition circuit, the position information calculated by the positioning method is prone to distortion, so corresponding position mapping methods must be used for correction. Summary of the Invention

[0004] Aiming at the deficiencies of the prior art, the present invention aims to provide a method for judging the regional division result in the generation of a position mapping table, including: after regionally identifying the boundary position distribution map, making a first judgment on the existence of bonding points according to whether the boundary division result is 64 partitions. Among them, the first judgment on the existence of bonding points is mainly for judging whether there are bonding points for the scattered points corresponding to the adjacent crystal bars in the outer circle; judging whether the partition number is within the range of 0 to 63 after region sequence coding; judging whether the count of pixels corresponding to the crystal bars in each partition is greater than 0; when the boundary division result after the first judgment on the existence of bonding points is 64 partitions, and the partition number after region sequence coding is within the range of 0 to 63, and the count of pixels corresponding to the crystal bars in each partition is greater than 0, the regional division result is qualified; and when the boundary division result after the first judgment on the existence of bonding points is 64 partitions, the boundary attribution and region sequence coding can be directly carried out.

[0005] According to the above method for judging the regional division result in the generation of a position mapping table of the present invention, further, it also includes: when the boundary division result after the first judgment on the existence of bonding points is less than 64 partitions, it is because there are bonding points for the scattered points corresponding to the adjacent crystal bars in the outer circle, and bonding point processing is carried out. The bonding point processing includes: counting the number of regions in the X direction and the Y direction of the outer circle; if the number of regions in a certain direction is less than 8, and the difference between the number of regions in this direction and 8 is n, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into n + 1 in this direction, so as to be divided into 8 in this direction, where n is an integer greater than 0 and less than 8.

[0006] According to the above method for judging the regional division result in the generation of a position mapping table of the present invention, further, n is 1, that is, the number of regions in a certain direction is 7, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into 2 in this direction, so as to be divided into 8 in this direction; or n is 2, that is, the number of regions in a certain direction is 6, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into 3 in this direction, so as to be divided into 8 in this direction.

[0007] According to the above method for judging the regional division result in the generation of a position mapping table of the present invention, further, it also includes: after the bonding point processing, making a second judgment on the existence of bonding points, judging whether the boundary division result is 64 partitions. Among them, the second judgment on the existence of bonding points is mainly for judging whether there are bonding points for the scattered points corresponding to the adjacent crystal bars inside; and after the boundary division result after the second judgment on the existence of bonding points is 64 partitions, then carry out the boundary attribution and region sequence coding.

[0008] According to the method for determining the result of region division in the generation of the position mapping table described above in the present invention, further, it further includes: when the boundary division result after the determination of the existence of the second bonding point is less than 64 partitions, re-partition it manually according to the pixel size of the scatter plot obtained from the statistics of the interaction between the crystal bar and gamma photons to obtain the position mapping table.

[0009] According to the method for determining the result of region division in the generation of the position mapping table described above in the present invention, further, it further includes: when the partition number is not within the range of 0 to 63 after the region sequence coding or the count of the pixels corresponding to the crystal bar in one or more partitions is 0, re-partition it manually according to the pixel size of the scatter plot obtained from the statistics of the interaction between the crystal bar and gamma photons to obtain the position mapping table.

[0010] The present invention also provides a system for determining the result of region division in the generation of the position mapping table, including: a bonding point existence judgment module, a bonding point processing module, and a partition coding judgment module. The bonding point existence judgment module includes a first bonding point existence judgment module and a second bonding point existence judgment module. The first bonding point existence judgment module is directly connected between the region identification module and the boundary attribution module. The bonding point processing module and the second bonding point existence judgment module are sequentially connected between the first bonding point existence judgment module and the boundary attribution module. The partition coding judgment module includes a first partition coding judgment module and a second partition coding judgment module. The first partition coding judgment module and the second partition coding judgment module are sequentially connected between the region sequence encoder and the position mapping table generation module.

[0011] According to the system for determining the result of region division in the generation of the position mapping table described above in the present invention, further, the bonding point processing module integrates an X-direction region number statistics unit, a Y-direction region number statistics unit, and a processing unit. The input ports and output ports of the X-direction region number statistics unit and the Y-direction region number statistics unit are respectively connected to the output port of the first bonding point existence judgment module and the input port of the processing unit. The output port of the processing unit is connected to the second bonding point existence judgment module.

[0012] According to the system for determining the result of region division in the generation of the position mapping table described above in the present invention, further, it further includes: a scatter plot manual boundary division module. The input ports of the scatter plot manual boundary division module are respectively connected to the output ports of the second bonding point existence judgment module, the first partition coding judgment module, and the second partition coding judgment module. The output port of the scatter plot manual boundary division module is connected to the position mapping table generation module. Description of the Drawings

[0013] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented in accordance with the contents of the specification. In order to make the above technical solution and other features and advantages of the present invention more obvious and easy to understand, the following specifically cites preferred embodiments and describes them in detail with reference to the accompanying drawings.

[0014] Figure 1 The present invention is a schematic structural diagram of a system for generating a method for mapping a crystal bar to a gamma photon interaction position;

[0015] Figure 2 The figure is a schematic diagram of the structure of the crystal bar center position estimation module of the system of the method for generating the crystal bar and gamma photon action position mapping table. In addition, the figure shows the connection relationship between the median filter processor, the mean filter processor, the crystal bar center position estimation module and the distance calculation module. The rest of the system is similar to Figure 1 Similar but omitted;

[0016] Figure 3 The figure is a schematic diagram of the structure of the crystal strip boundary determination module of the system for generating a method for mapping the position of crystal strips and gamma photons. In addition, the figure shows the connection relationship between the distance calculation module, the crystal strip boundary determination module and the area identification module. The rest of the system is similar to Figure 1 Similar but omitted;

[0017] Figure 4 The present invention is a flow chart of a method for generating a mapping table of crystal strips and gamma photon interaction positions;

[0018] Figure 5 It is a two-dimensional position spectrum obtained by statistics of the method for generating a mapping table of the interaction positions of crystal strips and gamma photons, i.e., a scatter plot;

[0019] Figure 6 The peak search result along the X direction of the scatter plot of a method for generating a mapping table of crystal strips and gamma photon interaction positions;

[0020] Figure 7 The peak search result along the Y direction of the scatter plot of a method for generating a mapping table of crystal strips and gamma photon interaction positions;

[0021] Figure 8 A map of the center position of a crystal strip estimated by a method for generating a mapping table of the interaction positions of a crystal strip and gamma photons, i.e., a binary position spectrum;

[0022] Figure 9 The distance distribution diagram of a method for generating a mapping table of the interaction position between a crystal bar and a gamma photon is shown. Figure 8 The local distance distribution map corresponding to area A shown in FIG;

[0023] Figure 10 It is a display diagram of the boundary division result of a method for generating a mapping table of the interaction position between a crystal bar and gamma photons;

[0024] Figure 11 For Figure 10 It is a display diagram of the crystal efficiency of each partition after boundary division in

[0025] Figure 12 It is a position mapping table of a method for generating a mapping table of the interaction position between a crystal bar and gamma photons, only showing the partial position mapping table corresponding to the B region shown in Figure 10 ;

[0026] Figure 13 It is a structural diagram of a system for determining and processing the region division result in the generation of a position mapping table disclosed by the present invention;

[0027] Figure 14 It is a structural diagram of the bonding point processing module of a system for determining and processing the region division result in the generation of a position mapping table disclosed by the present invention;

[0028] Figure 15 It is a flowchart of a method for determining and processing the region division result in the generation of a position mapping table of the system disclosed by the present invention. Detailed implementation manners

[0029] Next, in combination with the accompanying drawings and specific implementation manners, the present invention will be further described:

[0030] As Figure 1 shown, the generation of the position mapping table mentioned in the present invention is executed by a system for generating a mapping table of the interaction position between a crystal bar and gamma photons, which includes an array crystal 01, a data statistical detector 02, a median filtering processor 03, mean filtering processors 04 and 05, a central position estimation module 06 of the crystal bar, a distance calculation module 07, a crystal bar boundary determination module 08, a region identification module 09, a boundary attribution module 10, a region sequence encoder 11, a position mapping table generation module 12, and a position mapping table output module 13. The array crystal 01 has a plurality of crystal bars (not shown in the figure). The data statistical detector 02 is coupled with the array crystal 01. The data statistical detector 02, the median filtering processor 03, the mean filtering processors 04 and 05, the central position estimation module 06 of the crystal bar, the distance calculation module 07, the crystal bar boundary determination module 08, the region identification module 09, the boundary attribution module 10, the region sequence encoder 11, the position mapping table generation module 12, and the position mapping table output module 13 are connected in sequence.

[0031] According to the system of the method for generating a mapping table of crystal bars and gamma photon interaction positions of the present invention, further, the data statistics detector 02 is a PET detector. In one embodiment of the present invention, an 8*8 BGO array crystal and a Hamamatsu PMT coupled PET detector are used for data acquisition, and the time position distribution of the incident photons detected by the detector 02 is statistically analyzed, and the position distribution diagram is displayed in a 64*64 pixel image. Figure 5 As shown, the total number of incident photon times is 736124.

[0032] In addition, according to the system of the method for generating a mapping table of crystal bar and gamma photon interaction positions of the present invention, further, the median filter processor 03 and the mean filter processors 04 and 05 are all Gaussian filter processors, the median filter processor 03 is a median filter processor of a 3*3 two-dimensional template, and the mean filter processors 04 and 05 include an X-direction mean filter processor 04 and a Y-direction mean filter processor 05, and the X-direction mean filter processor 04 and the Y-direction mean filter processor 05 are both connected to the median filter processor 03 and the crystal bar center position estimation module 06, the X-direction mean filter processor 04 is a mean filter processor of a 1*3 template, and the Y-direction mean filter processor 05 is a mean filter processor of a 3*1 template.

[0033] like Figure 2 As shown, according to the system of the method for generating a crystal bar and gamma photon interaction position mapping table of the present invention, further, the crystal bar center position estimation module 06 has an "AND operation" processor 14 and a binarization processor 15, and the "AND operation" processor 14 and the binarization processor 15 are integrated in series into the crystal bar center position estimation module 06.

[0034] According to the system of the method for generating a mapping table of crystal bars and gamma photons interaction positions of the present invention, further, the distance calculation module 07 is a Euclidean distance calculator.

[0035] like Figure 3 As shown, according to the system of the method for generating a mapping table of crystal bar and gamma photon action positions of the present invention, further, the crystal bar boundary determination module 08 includes an X-direction one-dimensional peak-seeking unit 16, a Y-direction one-dimensional peak-seeking unit 17 and an "OR operation" processor 18. The X-direction one-dimensional peak-seeking unit 16, the Y-direction one-dimensional peak-seeking unit 17 and the "OR operation" processor 18 are integrated into the crystal bar boundary determination module 08. Both ends of the X-direction one-dimensional peak-seeking unit 16 and the Y-direction one-dimensional peak-seeking unit 17 are respectively connected to the distance calculation module 07 and the "OR operation" processor 18. The other end of the "OR operation" processor 18 is connected to the area identification module 09. The encoder of the area identification module 09 is a random encoder.

[0036] likeFigure 4 As shown in the figure, a method for generating a position mapping table related to the interaction position of a crystal bar and gamma photons mentioned in the present invention includes the following steps:

[0037] Step 1: The data statistical detector 02 performs two-dimensional statistics on the position information of incident photon events to obtain a two-dimensional position spectrum, that is, a scatter plot, as Figure 5 shown. The brightness of each pixel point in the two-dimensional position spectrum represents the number of events at that position.

[0038] Step 2: The median filter processor 03 performs median filtering on the obtained two-dimensional position spectrum using a 3×3 two-dimensional template.

[0039] Step 3: The X-direction mean filter processor 04 and the Y-direction mean filter processor 05 respectively perform mean filtering on the two-dimensional position spectrum along the X direction and the Y direction using 1×3 and 3×1 templates.

[0040] Among them, in the above technical solution, steps 2 and 3 are to remove the singular points in the two-dimensional position spectrum, reduce the influence of statistical fluctuations on the original data, and avoid the situation of over-division in the subsequent boundary division method. In step 3, mean filtering with 1×3 and 3×1 templates is performed on the two-dimensional position spectrum along the X direction and the Y direction respectively to obtain the two-dimensional position spectrum 2X and the two-dimensional position spectrum 2Y, and then the average of the same pixel positions of these two position spectra is taken to obtain the processed two-dimensional position spectrum, that is, the value of the pixel (X, Y) in the processed two-dimensional position spectrum = (two-dimensional position spectrum 2X(X, Y) + two-dimensional position spectrum 2Y(X, Y)) / 2.

[0041] Step 4: For the processed two-dimensional position spectrum, the central position estimation module 06 estimates the central position of each crystal bar using the local maximum method for the two-dimensional position spectrum, and the binarization processor 15 performs binarization processing on the data. The pixel value of the brightest point is set to 1, and the values of the remaining pixels are set to 0 to obtain a binarized position map, as Figure 8 shown.

[0042] Among them, when the central position estimation module 06 performs one-dimensional peak searching on the processed two-dimensional position spectrum, we use the derivative peak searching method. One-dimensional peak searching is performed on the processed two-dimensional position spectrum along the X direction and the Y direction respectively. The pixel point value at the peak position is set to 1, and the pixel points at non-peak positions are set to 0. See Figure 6 and Figure 7 , Figure 6 is the peak searching result along the X direction; Figure 7This is the peak search result along the Y direction. Then, the peak search results obtained along the X direction and the Y direction are subjected to an "AND operation" by the "AND operation" processor 14 at the same pixel position. If this pixel point is a peak point along both the X direction and the Y direction, then the binarization processor 15 sets the value of this pixel point to 1; otherwise, it sets it to 0. The binarization process is completed. Finally, the center position of each crystal bar is estimated, as shown in Figure 8 shown.

[0043] Step Five: The distance calculation module 07 calculates the distance from each pixel point to the nearest bright point (the point with a pixel value of 1), obtaining a distance distribution map, as shown in Figure 9 shown.

[0044] Among them, the Euclidean distance from each pixel point to the nearest bright point is calculated to obtain a distance distribution map. The calculated results are shown in Figure 9 shown. Only a part of it is shown. In the distance distribution map, the crystal bar boundary determination module 08 determines that the local maximum value is the boundary position corresponding to the crystal bar we need.

[0045] Step Six: The one-dimensional peak search unit 16 in the X direction and the one-dimensional peak search unit 17 in the Y direction respectively perform local maximum value search on the distance distribution map along the X direction and the Y direction to determine the center position of the crystal bar and divide the boundaries of each crystal bar.

[0046] Among them, the local maximum value method is also used in this Step Six. The specific steps are as follows:

[0047] The one-dimensional peak search unit 16 in the X direction and the one-dimensional peak search unit 17 in the Y direction perform one-dimensional peak search on the distance distribution map along the X direction and the Y direction respectively using the derivative peak search method;

[0048] The "OR operation" processor 18 performs an "OR operation" on the peak search results in the X direction and the Y direction at the same pixel point. The value with the largest number corresponds to the boundary position of each crystal bar, obtaining a boundary position distribution map, as shown in Figure 10 shown; at the same time, the crystal efficiency of each partition is calculated, as shown in Figure 11 shown.

[0049] Step Seven: Using the region identification method, each region is corresponding to a crystal bar number one by one to complete the mapping relationship between the scatter plot and the crystal bar.

[0050] The specific steps are as follows:

[0051] The region identification module 09 performs region identification on the boundary position distribution map. At this time, the formed region identification is not encoded in sequence. It only calculates which pixel points belong to the same region. The pixel points belonging to the same region use the same encoding, and the encodings between each region are ensured not to repeat;

[0052] The boundary position obtained in Step 6 above occupies one pixel, and it is necessary to determine which region this pixel point belongs to. According to the characteristics of the one-dimensional peak-seeking method adopted in Step 6, the boundary attribution module 10 makes boundary attribution judgments according to the following rules: The boundary points in the X direction belong to the region on the right, and the boundary points in the Y direction belong to the region below;

[0053] The region sequence encoder 11 adopts a sequence coding method to make each region correspond to the crystal table coding sequence, and completes the mapping relationship between the scatter plot and the crystal bar.

[0054] And Step 8: The position mapping table generation module 12 adopts a boundary extraction position mapping method to generate a position mapping table. As Figure 12 shown, only the partial position mapping table corresponding to Region B shown in Figure 10 is shown. Among them, it is required that the generated position mapping table should match the subsequent software interface, and the generated position mapping table is output from the system through the position mapping table output module 13.

[0055] As Figure 13 shown, a system for determining the result of region division in the generation of a position mapping table disclosed in the present invention includes: an adhesive point existence judgment module, an adhesive point processing module 20, and a partition coding judgment module. The adhesive point existence judgment module includes a first adhesive point existence judgment module 19 and a second adhesive point existence judgment module 21. The first adhesive point existence judgment module 19 is directly connected between the region identification module 09 and the boundary attribution module 10. The adhesive point processing module 20 and the second adhesive point existence judgment module 21 are sequentially connected between the first adhesive point existence judgment module 19 and the boundary attribution module 10. The partition coding judgment module includes a first partition coding judgment module 22 and a second partition coding judgment module 23. The first partition coding judgment module 22 and the second partition coding judgment module 23 are sequentially connected between the region sequence encoder 11 and the position mapping table generation module 12.

[0056] As Figure 14 shown, according to the above-mentioned system for determining the result of region division in the generation of a position mapping table of the present invention, further, the adhesive point processing module 20 integrates an X-direction region number statistics unit 25, a Y-direction region number statistics unit 26, and a processing unit 27. The input ports and output ports of the X-direction region number statistics unit 25 and the Y-direction region number statistics unit 26 are respectively connected to the output port of the first adhesive point existence judgment module 19 and the input port of the processing unit 27. The output port of the processing unit 27 is connected to the second adhesive point existence judgment module 21.

[0057] As Figure 13As shown in the figure, a system for determining the result of regional division in the generation of a position mapping table according to the above-mentioned present invention further includes: a scatter plot manual boundary division module 24. The input ports of the scatter plot manual boundary division module 24 are respectively connected to the output ports of the second bonding point existence judgment module 21, the first partition code judgment module 22, and the second partition code judgment module 23. The output port of the scatter plot manual boundary division module 24 is connected to the position mapping table generation module 12. The structure of a system for a method of generating a position mapping table between a crystal bar and a gamma photon interaction, which is outside the system for determining and processing the result of regional division in the generation of a position mapping table disclosed in the present invention, is as Figures 1 to 3 the same as the above relevant description, which is omitted here.

[0058] As Figure 15 shown, a method for determining the result of regional division in the generation of a position mapping table of the system disclosed in the present invention includes:

[0059] After the region identification module 09 performs region identification on the boundary position distribution map, the first bonding point existence judgment module 19 makes a first bonding point existence judgment according to whether the boundary division result is 64 partitions. Among them, the first bonding point existence judgment is mainly for judging whether there is a bonding point for the scatter points corresponding to the adjacent crystal bars in the outer circle;

[0060] When the boundary division result after the first bonding point existence judgment is less than 64 partitions, it is because there is a bonding point for the scatter points corresponding to the adjacent crystal bars in the outer circle, and the bonding point processing module 20 performs bonding point processing. The bonding point processing includes: the X-direction region number statistical unit 25 and the Y-direction region number statistical unit 26 perform statistical counting on the regions in the X direction and the Y direction of the outer circle; the processing unit 27 proceeds as follows: if the number of regions in a certain direction is less than 8, and the difference between the number of regions in this direction and 8 is n, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count in this direction into n + 1, so as to be divided into 8 in this direction, where n is an integer greater than 0 and less than 8;

[0061] That is: when n is 1, that is, the number of regions in a certain direction is 7, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count in this direction into 2, so as to be divided into 8 in this direction; or when n is 2, that is, the number of regions in a certain direction is 6, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count in this direction into 3, so as to be divided into 8 in this direction;

[0062] After the bonding point processing, the second bonding point existence judgment module 21 makes a second judgment on the existence of the bonding point, and judges whether the boundary division result is 64 partitions. Among them, the second judgment on the existence of the bonding point is mainly for judging whether there is a bonding point for the scattered points corresponding to the adjacent internal crystal bars; and after the boundary division result after the second judgment on the existence of the bonding point is 64 partitions, the boundary attribution module 10 and the region sequence encoder 11 perform boundary attribution and region sequence encoding again;

[0063] If the boundary division result after the second judgment on the existence of the bonding point is less than 64 partitions, the scattered point diagram manual boundary division module 24 re-partitions it manually according to the pixel size of the scattered point diagram obtained from the statistics of the interaction between the crystal bar and the gamma photon, and the position mapping table generation module 12 obtains the position mapping table;

[0064] After the region sequence encoding, the first partition encoding judgment module 22 judges whether the partition number is within the range of 0 to 63;

[0065] The second partition encoding judgment module 23 judges whether the count of the pixels corresponding to the crystal bar in each partition is greater than 0;

[0066] When the boundary division result after the first judgment on the existence of the bonding point is 64 partitions, and the partition number after the region sequence encoding is within the range of 0 to 63, and the count of the pixels corresponding to the crystal bar in each partition is greater than 0, the region division result is qualified. Because during the region sequence encoding process, due to the existence of the bonding point, even after the bonding point processing, some regions will still be missed, and these regions are not recognized in the region sequence encoding, so the count of the pixels corresponding to the crystal bar in them will be displayed as 0; and

[0067] When the boundary division result after the first judgment on the existence of the bonding point is 64 partitions, the boundary attribution and region sequence encoding can be directly performed;

[0068] If the partition number is not within the range of 0 to 63 after the region sequence encoding or the count of the pixels corresponding to the crystal bar in one or more partitions is 0, the scattered point diagram manual boundary division module 24 re-partitions it manually according to the pixel size of the scattered point diagram obtained from the statistics of the interaction between the crystal bar and the gamma photon, and the position mapping table generation module 12 obtains the position mapping table.

[0069] The present invention provides a method and system for judging the region division result in the generation of the position mapping table. During the process of generating the position mapping table of the interaction position between the crystal bar and the gamma photon, when there is a bonding point for the scattered points corresponding to adjacent crystal bars during the region division, the bonding point will directly affect the accuracy of the region division and the sequence numbering. The present invention mainly proofreads this to determine the position mapping table of the system.

[0070] Those skilled in the art can make various corresponding changes and deformations according to the technical solutions and concepts described above, and all such changes and deformations should fall within the protection scope of the claims of the present invention.

Claims

1. A method for determining the result of region division in the generation of a position mapping table, characterized in that: The method includes: After region identification of the boundary position distribution map, a first determination of the existence of bonding points is made according to whether the boundary division result is 64 partitions. Among them, the first determination of the existence of bonding points is mainly for determining whether there are bonding points in the scatter points corresponding to adjacent crystal bars in the outer ring; when the boundary division result after the first determination of the existence of bonding points is less than 64 partitions, it is because there are bonding points in the scatter points corresponding to adjacent crystal bars in the outer ring, and bonding point processing is performed. The bonding point processing includes: Counting the number of regions in the X direction and the Y direction for the outer ring; If the number of regions in a certain direction is less than 8, and the difference between the number of regions in this direction and 8 is n, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into n + 1 in this direction, so as to be divided into 8 in this direction, where n is an integer greater than 0 and less than 8; when n is 1, that is, the number of regions in a certain direction is 7, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into 2 in this direction, so as to be divided into 8 in this direction; or when n is 2, that is, the number of regions in a certain direction is 6, then count the total count of each region in this direction, and evenly divide the pixel points of the region with the largest count into 3 in this direction, so as to be divided into 8 in this direction; After the bonding point processing, a second determination of the existence of bonding points is made to determine whether the boundary division result is 64 partitions. Among them, the second determination of the existence of bonding points is mainly for determining whether there are bonding points in the scatter points corresponding to adjacent crystal bars inside; and after the boundary division result after the second determination of the existence of bonding points is 64 partitions, boundary attribution and region sequence coding are performed; when the boundary division result after the second determination of the existence of bonding points is less than 64 partitions, a position mapping table is obtained by manually re-partitioning it according to the pixel size of the scatter point diagram obtained by counting the interaction between crystal bars and gamma photons; After the region sequence coding, it is determined whether the partition number is within the range of 0 to 63; when the partition number after the region sequence coding is not within the range of 0 to 63 or the count of pixels corresponding to crystal bars in one or more partitions is 0, a position mapping table is obtained by manually re-partitioning it according to the pixel size of the scatter point diagram obtained by counting the interaction between crystal bars and gamma photons; Determine whether the count of pixels corresponding to crystal bars in each partition is greater than 0; When the boundary division result after the first determination of the existence of bonding points is 64 partitions, and the partition number after the region sequence coding is within the range of 0 to 63, and the count of pixels corresponding to crystal bars in each partition is greater than 0, the region division result is qualified; and When the boundary division result after the first adhesive point existence judgment is 64 partitions, the boundary attribution and area sequence coding can be directly carried out; if the partition number is not within the range of 0 to 63 after the area sequence coding or the count of pixels corresponding to the crystal bars in one or more partitions is 0, the pixels of the scatter plot obtained according to the statistics of the interaction between the crystal bars and gamma photons are manually repartitioned to obtain a position mapping table.

2. A region division result determination processing system in position mapping table generation, characterized in that: The system includes: an adhesive point existence judgment module, an adhesive point processing module, and a partition coding judgment module. The adhesive point existence judgment module includes a first adhesive point existence judgment module and a second adhesive point existence judgment module. The first adhesive point existence judgment module is directly connected between the area identification module and the boundary attribution module. The adhesive point processing module and the second adhesive point existence judgment module are sequentially connected between the first adhesive point existence judgment module and the boundary attribution module. The partition coding judgment module includes a first partition coding judgment module and a second partition coding judgment module. The first partition coding judgment module and the second partition coding judgment module are sequentially connected between the area sequence encoder and the position mapping table generation module; The system further includes a scatter plot manual boundary division module. The input ports of the scatter plot manual boundary division module are respectively connected to the output ports of the second adhesive point existence judgment module, the first partition coding judgment module, and the second partition coding judgment module. The output port of the scatter plot manual boundary division module is connected to the position mapping table generation module; Among them, the first adhesive point existence judgment module makes the first adhesive point existence judgment according to whether the boundary division result is 64 partitions. Among them, the first adhesive point existence judgment is mainly for judging whether there is an adhesive point in the scatter points corresponding to the outer ring adjacent crystal bars; When the boundary division result after the first adhesive point existence judgment is less than 64 partitions, it is because there is an adhesive point in the scatter points corresponding to the outer ring adjacent crystal bars, and the adhesive point processing module performs adhesive point processing. The adhesive point processing includes: the X-direction area number statistical unit and the Y-direction area number statistical unit perform the X-direction and Y-direction area number statistics on the outer ring; the processing unit proceeds as follows: if the area number in a certain direction is less than 8, the difference between the area number in this direction and 8 is n, then the total count of each area in this direction is statistically calculated, and the pixel points of the area with the largest count are evenly divided into n + 1 in this direction, so that it is divided into 8 in this direction, where n is an integer greater than 0 and less than 8; when n is 1, that is, the area number in a certain direction is 7, then the total count of each area in this direction is statistically calculated, and the pixel points of the area with the largest count are evenly divided into 2 in this direction, so that it is divided into 8 in this direction; or when n is 2, that is, the area number in a certain direction is 6, then the total count of each area in this direction is statistically calculated, and the pixel points of the area with the largest count are evenly divided into 3 in this direction, so that it is divided into 8 in this direction; After the bonding point processing, the second bonding point presence judgment module makes a second judgment on the presence of bonding points, determining whether the boundary division result is 64 partitions. Among them, the second judgment on the presence of bonding points mainly focuses on whether there are bonding points for the scattered points corresponding to adjacent internal crystal bars; and after the boundary division result after the second judgment on the presence of bonding points is 64 partitions, the boundary attribution module and the region sequence encoder perform boundary attribution and region sequence encoding; if the boundary division result after the second judgment on the presence of bonding points is less than 64 partitions, the scattered point diagram manual boundary division module re-partitions it manually according to the pixel size of the scattered point diagram obtained from the statistics of the interaction between the crystal bar and gamma photons, and the position mapping table generation module obtains the position mapping table. After the region sequence encoding, the first partition encoding judgment module determines whether the partition number is within the range of 0 to 63. After the region sequence encoding, when the partition number is not within the range of 0 to 63 or the count of pixels corresponding to the crystal bar in one or more partitions is 0, it is re-partitioned manually according to the pixel size of the scattered point diagram obtained from the statistics of the interaction between the crystal bar and gamma photons to obtain the position mapping table. The second partition encoding judgment module determines whether the count of pixels corresponding to the crystal bar in each partition is greater than 0; when the boundary division result after the first judgment on the presence of bonding points is 64 partitions, and the partition number after the region sequence encoding is within the range of 0 to 63, and the count of pixels corresponding to the crystal bar in each partition is greater than 0, the region division result is qualified. When the boundary division result after the first judgment on the presence of bonding points is 64 partitions, boundary attribution and region sequence encoding can be directly performed; if the partition number is not within the range of 0 to 63 or the count of pixels corresponding to the crystal bar in one or more partitions is 0 after the region sequence encoding, the scattered point diagram manual boundary division module re-partitions it manually according to the pixel size of the scattered point diagram obtained from the statistics of the interaction between the crystal bar and gamma photons, and the position mapping table generation module obtains the position mapping table.

Images