Distribution control method for anti-overflow administration of isotope application

Abstract

The invention belongs to the technical field of medical automation equipment, and discloses a distribution control method for isotope application anti-overflow administration, which comprises the following steps: acquiring image data of an isotope application drug carrying layer to be administered; based on the image data, extracting a contour line of the medicine carrying layer, and determining a minimum enclosing rectangle of the contour line; according to the length-width ratio of the minimum circumscribed rectangle, performing standardized segmentation on the circumscribed rectangle to form a plurality of primary regions; performing recursive subdivision operation on each primary region, and dividing the current region into a plurality of sub-regions in each time of division until a preset division stop condition is met; screening the finally obtained sub-region, and reserving the sub-region containing the effective region as a drug delivery sub-region; performing boundary corrosion treatment on each drug delivery sub-area, and performing binary search on the optimal corrosion diameter of each drug delivery sub-area; through iterative optimization of the corrosion diameter, the calculation efficiency of optimal corrosion diameter search is significantly improved, and the risk of radioactive liquid medicine overflow is reduced.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of medical automation equipment, and particularly relates to a distribution control method for preventing radioactive liquid overflow of isotopic patch. BACKGROUND

[0002] The core of isotopic patch treatment is to prepare a personalized patch that completely matches the scar shape of a patient. Due to the complex and variable scar shapes commonly seen in clinical practice, such as irregular structures like spiral, centipede or bifurcation, the traditional preparation process needs to accurately cut the filter paper carrier according to the scar profile and ensure that the radioactive isotopic phosphorus-32 liquid can uniformly cover the entire patch area. At present, the ADP drug delivery device is generally used for titration of the liquid, and the liquid is usually slowly dripped from the center of the filter paper, and the coverage is realized by relying on the natural diffusion of the liquid. Although this operation method is barely feasible when dealing with regular convex polygon filter papers, it has many inherent defects when facing complex scar carriers. The process of liquid diffusion from the center to the outside not only takes a long time, significantly reducing the preparation efficiency, but also, due to the physical properties of the surface tension of the liquid, the liquid forms local accumulation in narrow areas, and it is difficult to fully infiltrate in slender branches or corners. This uneven distribution will directly affect the treatment effect, and even may cause radioactive burns due to abnormal local dose; more seriously, when the filter paper center passage is too thin, the liquid may overflow the carrier edge due to the imbalance of surface tension, causing radioactive contamination risk and treatment dose out of control; although some researches have tried to use ant colony algorithm or physical diffusion model to optimize the dripping path, these algorithms often need large computing resources, and are difficult to be practically applied in clinical scenes with strict real-time requirements. Therefore, in the prior art, when titrating the liquid of complex-shaped isotopic patches, it is difficult to accurately control the liquid boundary, radioactive liquid overflow is easy to occur, and contamination risk is caused. SUMMARY

[0003] In view of the deficiencies in the prior art, the purpose of the present application is to provide a distribution control method for preventing radioactive liquid overflow of isotopic patch, which solves the problem that in the prior art, when titrating the liquid of complex-shaped isotopic patches, it is difficult to accurately control the liquid boundary, radioactive liquid overflow is easy to occur, and contamination risk is caused.

[0004] The purpose of the present application can be achieved by the following technical solutions:

[0005] A distribution control method for preventing radioactive liquid overflow of isotopic patch, comprising the following steps:

[0006] Obtaining image data of a drug-loaded layer of an isotopic patch to be administered;

[0007] Based on image data, the contour line of the drug-loaded layer is extracted, and the minimum bounding rectangle of the contour line is determined. The area inside the contour line in the minimum bounding matrix is ​​the effective area, and the rest is the blank area.

[0008] The bounding rectangle is normalized and divided according to its aspect ratio to form multiple primary regions;

[0009] Recursively subdivide each primary region, dividing the current region into several sub-regions each time, until the preset subdivision stopping condition is met;

[0010] Filter the finally obtained sub-regions and retain the sub-regions containing the effective regions as the drug delivery sub-regions;

[0011] Each drug delivery region is subjected to boundary erosion treatment, and the optimal erosion diameter of each drug delivery region is searched by binary search. After erosion, a drug delivery safety zone is formed in each drug delivery region.

[0012] Within each safe administration zone, the administration point and dosage are determined, and titration is performed at the corresponding administration point using administration equipment.

[0013] The shape of the drug-loaded layer in an isotope patch can be regular or irregular.

[0014] Regular shapes include circles, ellipses, or convex polygons;

[0015] Irregular shapes include forked, spiral, or irregular structures with recessed edges.

[0016] The bounding rectangle is normalized and divided according to its aspect ratio to form multiple primary regions. The specific steps include:

[0017] Define the longest side of the smallest bounding rectangle as length L and the shortest side as width W;

[0018] Preset segmentation threshold β;

[0019] When aspect ratio When the length of the bounding rectangle is N equal to the length of the bounding rectangle, then N is a number of regions. The integer value rounded to the nearest whole number;

[0020] When aspect ratio When the minimum bounding rectangle is used, it is directly taken as a primary region.

[0021] The segmentation threshold β is equal to 2.

[0022] For each primary region, a recursive subdivision operation is performed, where each subdivision divides the current region into several sub-regions. The specific steps include:

[0023] The recursive subdivision operation is implemented by a quadtree algorithm, which first divides the primary region along the horizontal and vertical directions to generate four equal sub-regions;

[0024] determine whether each sub-region meets a stop condition;

[0025] If a sub-region meets the stop condition, stop the recursive subdivision operation;

[0026] If a sub-region does not meet the stop condition, continue to subdivide the sub-region using the quadtree algorithm;

[0027] Repeat the iterative subdivision until each sub-region meets the stop condition.

[0028] The stop condition includes:

[0029] All of the sub-region is a valid region or a blank region;

[0030] The area of the sub-region is less than the area threshold value one;

[0031] When the sub-region meets any one of the conditions, it is determined that the stop condition is met.

[0032] The optimal erosion diameter of each drug sub-region is determined by binary search, including the following steps:

[0033] Determine the initial erosion diameter search range [D min , D max ] of each drug sub-region;

[0034] D min represents the lower limit of the erosion diameter, and the initial lower limit of the erosion diameter of each drug sub-region is the same preset value;

[0035] D max represents the upper limit of the erosion diameter, and the initial upper limit of the erosion diameter of each drug sub-region is the length of the short side of the corresponding drug sub-region;

[0036] Perform an iterative binary search on each drug sub-region, including the following steps:

[0037] S1, take the middle value D mid of the erosion diameter search range of the drug sub-region as the execution erosion diameter, D mid = (D min + D max ) / 2;

[0038] S2, corrode with the execution erosion diameter;

[0039] If there is a remaining region after corrosion, keep the upper limit D max of the erosion diameter, and take the middle value D mid as the new lower limit D min, update the corrosion diameter search range;

[0040] If there is no remaining area after corrosion, keep the lower limit D min of the corrosion diameter search range as the new upper limit D mid of the corrosion diameter search range; max

[0041] S3, preset a convergence threshold θ, and set the convergence condition as θ≤D max -D min ;

[0042] S4, judge whether the updated corrosion diameter search range meets the convergence condition;

[0043] If the convergence condition is met, the iterative dichotomy search of the current drug sub-zone is terminated;

[0044] If the convergence condition is not met, repeat S1 to S4 until the convergence condition is met;

[0045] S5, take the middle value D mid of the corrosion diameter search range when the convergence condition is met as the optimal corrosion diameter of the current drug sub-zone.

[0046] For each drug safety zone, identify all connected regions inside;

[0047] If there are multiple non-adjacent connected regions, calculate the area of each connected region respectively, and distribute the total drug amount of the drug sub-zone based on the area ratio;

[0048] Preset an area threshold two;

[0049] If the area of the connected region is less than or equal to the area threshold two, randomly select a point near the center point of the connected region as the drug point, and titrate according to the total drug amount corresponding to the connected region;

[0050] If the area of the connected region is greater than the area threshold two, divide the connected region into multiple sub-connected regions with an area less than or equal to the area threshold two, randomly select a point near the center point of each sub-connected region as the drug point, and distribute the drug amount of each sub-connected region according to the area ratio and titrate.

[0051] The lower limit D min of the corrosion diameter search range in each drug sub-zone is one pixel in the image data.

[0052] The convergence threshold θ is one pixel in the image data.

[0053] The beneficial effects of the present application are:

[0054] ​This invention uses image recognition technology to obtain the contour line of the drug-loaded layer, divides the region based on the minimum bounding rectangle of the contour line, and selects each drug delivery region. A binary search algorithm is used to quickly determine the optimal corrosion diameter for each drug delivery region, and corrosion processing is used to accurately locate the safe zone of the drug delivery point. This method employs dynamic boundary corrosion technology, iteratively optimizing the corrosion diameter to significantly improve computational efficiency while ensuring a safe distance between the drug delivery point and the contour boundary. This significantly improves the uniformity of drug distribution and drug delivery efficiency, while effectively reducing the risk of radioactive drug spillage. Attached Figure Description

[0055] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0056] Figure 1 This is a schematic diagram of the drug distribution algorithm of the present invention. Detailed Implementation

[0057] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0058] like Figure 1 As shown, a method for controlling the distribution of isotope patch administration to prevent spillage includes the following steps:

[0059] Acquire image data of the drug-loaded layer of the isotope to be administered;

[0060] Based on image data, the contour line of the drug-loaded layer is extracted, and the minimum bounding rectangle of the contour line is determined. The area inside the contour line in the minimum bounding matrix is ​​the effective area, and the rest is the blank area.

[0061] The bounding rectangle is normalized and divided according to its aspect ratio to form multiple primary regions;

[0062] Recursively subdivide each primary region, dividing the current region into several sub-regions each time, until the preset subdivision stopping condition is met;

[0063] Filter the finally obtained sub-regions and retain the sub-regions containing the effective regions as the drug delivery sub-regions;

[0064] The boundary of each administration sub-area is eroded, and the optimal erosion diameter of each administration sub-area is searched by binary search, and a safe administration area is formed in each administration sub-area after erosion;

[0065] A drug administration point and a drug administration amount are determined in each safe administration area, and titration is performed at the corresponding drug administration point by a drug administration device.

[0066] Preferably, the image data is obtained by a high-speed scanner from a drug-loaded layer of an isotope patch, and the obtained image data is an RGB image.

[0067] Preferably, the contour line of the drug-loaded layer is extracted, and the method specifically comprises the following steps:

[0068] The RGB image is converted into a grayscale image through grayscale processing;

[0069] The grayscale image is preprocessed;

[0070] The contour line of the drug-loaded layer in the preprocessed grayscale image is extracted.

[0071] The preprocessing of the grayscale image specifically comprises the following steps:

[0072] The grayscale image is smoothed and denoised by a Gaussian filter;

[0073] The contrast of the grayscale image is enhanced by histogram equalization.

[0074] The contour line of the drug-loaded layer in the preprocessed grayscale image specifically comprises the following steps:

[0075] The effective edges are extracted from the preprocessed grayscale image by using a Canny edge detection algorithm;

[0076] The isolated points are removed and the disconnected edges are connected by using morphological opening operation;

[0077] The threshold value L of the edge extraction length is set min ;

[0078] The edges with a length greater than or equal to the threshold value L min are extracted by connected region analysis, the edges with a length less than the threshold value L min are removed, and finally the edge point set P = {(x i , y i ) | i = 1, 2,..., N} is obtained;

[0079] The edge point set P is sorted by using a nearest neighbor sorting algorithm, and the sorted point set is optimized to generate an ordered point set;

[0080] The ordered point set is fitted to generate a smooth curve by a cubic spline interpolation method, and the end point of the curve is aligned with the starting point to make the curve closed;

[0081] The curve is converted according to the scaling factor between the pixel size of the RGB image and the actual size of the drug-loaded layer, and finally the contour line with the actual physical size is obtained.

[0082] The shape of the isotope patch drug-loaded layer includes a regular shape or an irregular shape;

[0083] The regular shape includes a circular shape, an elliptical shape or a convex polygon;

[0084] The irregular shape includes a bifurcated shape, a spiral shape or a special-shaped structure with edge recess.

[0085] According to the aspect ratio of the minimum circumscribed rectangle, the circumscribed rectangle is normalized and segmented to form a plurality of primary regions, specifically including the following steps:

[0086] The long side of the minimum circumscribed rectangle is defined as the length L, and the short side is defined as the width W;

[0087] The preset segmentation threshold β;

[0088] When the aspect ratio is less than 2, the circumscribed rectangle is evenly divided into N primary regions along the long side, and N is The rounding integer value;

[0089] When the aspect ratio is greater than or equal to 2, the minimum circumscribed rectangle is directly used as a primary region.

[0090] The segmentation threshold β is equal to 2; when the aspect ratio is less than 2, the region shape is relatively regular, and the four-tree process can be directly entered; when the aspect ratio is greater than or equal to 2, the region is too long and narrow, and the phosphorus-32 liquid is easy to diffuse unevenly or escape from the edge due to surface tension; by pre-segmenting into sub-regions with an aspect ratio close to 1:1, the uniform and controllable diffusion of the liquid in each region is ensured, and the generation of abnormally narrow sub-nodes is avoided from the algorithm level, significantly improving the uniformity of the radio dose, the safety of the treatment and the overall robustness of the algorithm.

[0091] Recursive subdivision operation is performed on each primary region, and each time the current region is divided into sub-regions, specifically including the following steps:

[0092] The recursive subdivision operation is realized by the quadtree algorithm, which first divides the primary region into four equal sub-regions along the horizontal and vertical directions;

[0093] Determine whether the generated sub-regions meet the stop condition;

[0094] The sub-regions that meet the stop condition stop the recursive subdivision operation;

[0095] The sub-regions that do not meet the stop condition continue to be subdivided using the quadtree algorithm;

[0096] The iterative subdivision is repeated until each of the resulting sub-regions meets the stop condition.

[0097] The stop condition includes:

[0098] All of the sub-region is a valid region or a blank region;

[0099] The area of the sub-region is less than the area threshold one;

[0100] The area threshold one is determined by the minimum droplet accuracy of the drug delivery device;

[0101] When the sub-region meets any one of the conditions, it is determined that the stop condition is met.

[0102] The optimal etching diameter of each drug delivery sub-region is searched by binary search, including the following steps:

[0103] The initial etching diameter search range [D min , D max ] of each drug delivery sub-region is determined;

[0104] D min represents the lower limit of the etching diameter, and the initial lower limit of the etching diameter of each drug delivery sub-region is the same preset value;

[0105] D max represents the upper limit of the etching diameter, and the initial upper limit of the etching diameter of each drug delivery sub-region is the length of the short side of the corresponding drug delivery sub-region;

[0106] For each drug delivery sub-region, an iterative binary search is performed, including the following steps:

[0107] S1, taking the middle value D mid of the etching diameter search range of the drug delivery sub-region as the execution etching diameter, D mid = (D min + D max ) / 2;

[0108] S2, etching with the execution etching diameter;

[0109] If there is a remaining region after etching, the upper limit of the etching diameter D max is retained, the middle value D mid is taken as the new lower limit of the etching diameter D min , and the etching diameter search range is updated;

[0110] If there is no remaining region after etching, the lower limit of the etching diameter D min is retained, and the middle value D mid is taken as the new upper limit of the etching diameter Dmax , update the corrosion diameter search range;

[0111] S3, preset a convergence threshold θ, and set the convergence condition as θ≤D max -D min ;

[0112] S4, judge whether the updated corrosion diameter search range meets the convergence condition;

[0113] If the convergence condition is met, the iterative dichotomy search of the current drug sub-zone is terminated;

[0114] If the convergence condition is not met, repeat S1 to S4 until the convergence condition is met;

[0115] S5, take the middle value D of the corrosion diameter search range when the convergence condition is met as the optimal corrosion diameter of the current drug sub-zone. mid

[0116] For each drug safety zone, identify all connected regions inside;

[0117] If there are multiple non-adjacent connected regions, calculate the area of each connected region respectively, and distribute the total drug amount of the drug sub-zone based on the area ratio;

[0118] Preset an area threshold two;

[0119] If the area of the connected region is less than or equal to the area threshold two, randomly select a point near the center point of the connected region as the drug point, and titrate according to the total drug amount corresponding to the connected region;

[0120] If the area of the connected region is greater than the area threshold two, divide the connected region into multiple sub-connected regions with an area less than or equal to the area threshold two, randomly select a point near the center point of each sub-connected region as the drug point, distribute the drug amount of each sub-connected region according to the area ratio, and titrate.

[0121] It should be noted that the area threshold one is determined by the minimum single-point titration volume of the drug delivery device and the maximum diffusion area of the drug solution on the isotopic patch drug loading layer;

[0122] The area threshold two is determined by the maximum single-point titration volume of the drug delivery device and the maximum diffusion area of the drug solution on the isotopic patch drug loading layer.

[0123] Preferably, the lower limit D min of the corrosion diameter search range of each drug sub-zone is one pixel in the image data;

[0124] Preferably, the convergence threshold θ is one pixel in the image data;

[0125] ​It should be noted that the lower limit D of the corrosion diameter in the initial corrosion diameter search range of each drug administration sub-area min and the convergence threshold θ can be adjusted as needed, such as being appropriately adjusted according to the drug administration accuracy of the drug administration device.

[0126] The image recognition technology is used to obtain the contour line of the drug loading layer, the region is divided based on the minimum circumscribed rectangle of the contour line, and each drug administration sub-area is obtained through screening; the bisection search algorithm is used to quickly determine the optimal corrosion diameter of each drug administration sub-area, and the precise positioning of the drug administration point safety zone is realized through corrosion processing; the method uses dynamic boundary corrosion technology, and through iterative optimization of the corrosion diameter, the calculation efficiency is significantly improved while ensuring the safety distance between the drug administration point and the contour boundary, thereby significantly improving the uniformity of the drug solution distribution and the drug administration efficiency, and effectively reducing the risk of radioactive drug solution overflow.

[0127] In the description of the present specification, the description of the terms "one embodiment", "example", "specific example" and the like means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0128] The basic principles, main features and advantages of the present application are shown and described above. It should be understood by those skilled in the art that the present application is not limited by the above embodiments, and the above embodiments and descriptions in the specification are only to illustrate the principles of the present application. Without departing from the spirit and scope of the present application, various changes and improvements can be made to the present application, and these changes and improvements all fall within the scope of the claimed present application.

Claims

1. A method for controlling the distribution of isotope patch administration to prevent spillage, characterized in that, Includes the following steps: Acquire image data of the drug-loaded layer of the isotope to be administered; Based on image data, the contour line of the drug-loaded layer is extracted, and the minimum bounding rectangle of the contour line is determined. The area inside the contour line in the minimum bounding matrix is ​​the effective area, and the rest is the blank area. The bounding rectangle is normalized and divided according to its aspect ratio to form multiple primary regions; Recursively subdivide each primary region, dividing the current region into several sub-regions each time, until the preset subdivision stopping condition is met; Filter the finally obtained sub-regions and retain the sub-regions containing the effective regions as the drug delivery sub-regions; Each drug delivery region is subjected to boundary erosion treatment, and the optimal erosion diameter of each drug delivery region is searched by binary search. After erosion, a drug delivery safety zone is formed in each drug delivery region. Determine the administration point and dosage within the safe administration zone, and perform titration after administration.

2. The distribution control method for preventing spillage of isotope patch administration according to claim 1, characterized in that, The shape of the drug-loaded layer in an isotope patch can be regular or irregular. Regular shapes include circles, ellipses, or convex polygons; Irregular shapes include forked, spiral, or irregular structures with recessed edges.

3. The distribution control method for preventing spillage of isotope patch administration according to claim 1, characterized in that, The bounding rectangle is normalized and divided according to its aspect ratio to form multiple primary regions. The specific steps include: Define the longest side of the smallest bounding rectangle as length L and the shortest side as width W; Preset segmentation threshold β; When aspect ratio When the length of the bounding rectangle is N equal to the length of the bounding rectangle, then N is a number of regions. The integer value rounded to the nearest whole number; When aspect ratio When the minimum bounding rectangle is used, it is directly taken as a primary region.

4. The distribution control method for preventing spillage of isotope patch administration according to claim 3, characterized in that, The segmentation threshold β is equal to 2.

5. The distribution control method for preventing spillage of isotope patch administration according to claim 3, characterized in that, For each primary region, a recursive subdivision operation is performed, where each subdivision divides the current region into several sub-regions. The specific steps include: The recursive subdivision operation is implemented using the quadtree algorithm, which first divides the primary region into four equal subregions by dividing it into two equal parts horizontally and vertically. Determine whether each generated sub-region meets the stopping condition; For subregions that meet the stopping condition, stop the recursive subdivision operation; Subregions that do not meet the stopping condition are further subdivided using the four-part tree algorithm; Repeat the iterative subdivision process until all the resulting sub-regions meet the stopping condition.

6. The distribution control method for preventing spillage of isotope patch administration according to claim 5, characterized in that, Stopping conditions include: The entire sub-region is either a valid region or a blank region; The area of ​​the sub-region is less than the area threshold. If a subregion satisfies any one of the conditions, it is determined that the stopping condition has been met.

7. The distribution control method for preventing spillage of isotope patch administration according to claim 6, characterized in that, The optimal corrosion diameter for each dosing region is determined by binary search, specifically including the following steps: Determine the initial corrosion diameter search range for each dosing region [D] min D max ]; D min This indicates the lower limit of the corrosion diameter. The initial lower limit of the corrosion diameter for each dosing zone is the same preset value. D max This indicates the upper limit of the corrosion diameter. The initial upper limit of the corrosion diameter for each dosing region is the length of the shorter side of the corresponding dosing region. Perform an iterative binary search for each dosing region, including the following steps: S1. Take the median value D of the corrosion diameter search range of the dosing area. mid To achieve the corrosion diameter, D mid =(D min +D max ) / 2; S2. Perform corrosion according to the corrosion diameter; If a residual area exists after corrosion, the upper limit of the corrosion diameter D is retained. max With the median value D mid The new lower limit of corrosion diameter D min Update the search range for corrosion diameter; If no remaining area exists after corrosion, then retain the lower limit of the corrosion diameter D. min With the median value D mid The new upper limit for corrosion diameter D max Update the search range for corrosion diameter; S3. Preset convergence threshold θ, and set the convergence condition as θ≤D. max -D min ; S4. Determine whether the updated corrosion diameter search range meets the convergence condition. If the convergence condition is met, the iterative binary search of the current dosing region terminates. If the convergence condition is not met, repeat S1 to S4 until the convergence condition is met. S5. Take the median value D of the search range for the corrosion diameter that satisfies the convergence condition. mid This represents the optimal corrosion diameter for the current dosing region.

8. The distribution control method for preventing spillage of isotope patch administration according to claim 7, wherein determining the administration point and dosage within the safe administration zone and performing titration of the administration specifically includes the following steps: For each safe drug administration zone, identify all connected regions within it; If there are multiple non-adjacent connected regions, calculate the area of ​​each connected region separately, and allocate the total dosage to the herb region based on the area ratio. Preset area threshold two; If the area of ​​the connected region is less than or equal to the area threshold two, a point is randomly selected near the center point of the connected region as the drug administration point, and titration is performed according to the total drug administration amount corresponding to the connected region. If the area of ​​the connected region is greater than the area threshold two, the connected region is divided into multiple sub-connected regions with areas less than or equal to the area threshold two. A point is randomly selected near the center point of each sub-connected region as the drug administration point, and the drug dosage of each sub-connected region is allocated according to the area ratio and titrated.

9. The distribution control method for preventing spillage of isotope patch administration according to claim 7, characterized in that, The lower limit D of the corrosion diameter in the initial corrosion diameter search range for each dosing zone min Each of these represents a single pixel in the image data.

10. The distribution control method for preventing spillage of isotope patch administration according to claim 7, characterized in that, The convergence threshold θ is a pixel in the image data.

Images