Bulk phase biological printing scattering correction method and system, electronic equipment and storage medium

Abstract

The invention discloses a face-oriented volume phase biological printing scattering correction method and system, electronic equipment and a storage medium. The face-oriented volume phase biological printing scattering correction method comprises the steps that 100 side window images projected in different light modes are obtained; performing gray binarization on the side window image; extracting a gray value from each column of the side window image and calculating a Fourier transform result; calculating an average frequency spectrum according to a Fourier transform result; calculating to obtain a corrected spectrum according to the average spectrum; performing polar coordinate transformation on each column of the correction frequency spectrum to obtain a two-dimensional mask; performing two-dimensional Fourier transform point multiplication on each two-dimensional mask and the target cube slice to obtain a point multiplication result; performing two-dimensional inverse Fourier transform on the dot product result to obtain a correction result; and performing effective filtering on the correction result, and superposing filtering results of different depths to obtain a scattering correction result of the target cube slice. On the basis of revealing the scattering characteristic of printing ink in bulk phase printing, a scattering correction rotation chromatography algorithm is perfected, and bulk phase biological printing of a high-fidelity phase separation material is achieved.

Description

Technical Field

[0001] The present invention relates to the technical field of software deployment, and in particular to a method, system, electronic device and storage medium for body-oriented bioprinting scatter correction. Background Art

[0002] Three-dimensional bioprinting technology can manufacture complex and high-precision tissue structures and is suitable for biomedicine. 3D bioprinting accurately arranges biomaterials, cells and biological factors to mimic the anatomical characteristics of target tissues and promote tissue regeneration, playing an important role in the research of regenerative medicine and tissue engineering. With the rapid development of engineering and materials, 3D bioprinting technology is also constantly iterating and innovating, mainly including jet-type, extrusion-type bioprinting and digital light-curing (DLP) bioprinting. 3D bioprinting technology represented by inkjet and extrusion takes a long time to build structures, making it difficult to achieve 3D additive manufacturing of complex structures, and the range of compatible biomaterials is limited. Shear stress during printing damages cells. Two-dimensional layer curing printing technology based on digital light curing (DLP) has improved the printing speed to a certain extent, avoided the influence of shear stress on cell viability, and improved the printing resolution, but it is not compatible with high-viscosity bio-inks, and the low-dimensional molding unit may cause the mechanical properties of the printed samples to be anisotropic, which is difficult to predict or control. The printing speed of the above-mentioned traditional 3D bioprinting technology also limits its scalability in the field of bio-additive manufacturing, and the long-term biomanufacturing process will damage cell function. To overcome this problem, bulk bioprinting was developed. Bulk bioprinting changes the layer-by-layer operation of traditional equipment, allowing the creation of geometrically complex centimeter-level structures in tens of seconds. As a new generation of biomanufacturing technology, bulk bioprinting has a higher degree of control over the formed structure, does not require printing support, has a fast forming speed, and can protect cells from shear stress damage. Compared with traditional extrusion bioprinting technology or digital photocuring bioprinting technology, it has many advantages.

[0003] The goal of bioprinting technology is usually to print cell structures or customized and functional 3D cell scaffolds, which can reproduce the physiological microenvironment used for tissue formation and promote cell attachment, proliferation and other activities. The materials currently used for 3D bioprinting are cell-loaded hydrogel systems composed of pure natural polymers. The mechanical properties and printability are generally poor. It is difficult to provide sufficient support for cells in the body, and it is difficult to buffer the complex local tissue tension to maintain the stable three-dimensional structure of the printed body. It is impossible to maintain a long-term stable morphology. Although cells can maintain a certain activity in it, they do not stretch or proliferate, so they cannot form new tissues. For this reason, in many cases, it is necessary to bioprint a hydrogel scaffold with a microporous structure, which will allow nutrient and oxygen exchange and cell diffusion, thereby enhancing cell function. The ATPE system is a two-phase separation phenomenon in a mixture of two hydrophilic components at a critical concentration. Due to its excellent cell compatibility and other advantages, it has been explored and applied to the preparation of porous bio-inks. However, there will be light scattering in the two-phase solution based on body phase bioprinting, which affects the printing effect, so the scattering characteristics of the two-phase solution should be analyzed during printing.

[0004] Volumetric additive manufacturing manufactures centimeter-scale objects in a true 3D manner by irradiating the entire photocurable solution container with one or a group of light projections for photopolymerization. The volume energy dose generated by the accumulated light exposure solidifies the printed material into the desired geometry. Since volumetric additive manufacturing uses three-dimensional voxels as the molding unit, the printing effect is extremely sensitive to the propagation process of structured light in the printing space. The refraction or scattering of light inside the printing ink will significantly affect the molding accuracy. Therefore, the printing ink is generally required to be homogeneous and transparent, with a high scattering mean free path, so that the light mode used for photopolymerization will not be disrupted during its propagation. For aqueous two-phase emulsion system materials with stronger random refraction and scattering effects, when the light modes propagate through complex materials, they will be optically distorted, and the spatial information they carry will also be disrupted. The attenuation and distortion of the light modes that vary with depth will destroy the formation of the required light dose in the material. This directly leads to poor printing fidelity and global resolution loss, and its scattering correction volumetric printing method remains to be established. When using a volumetric printer, it is crucial to achieve scattering correction when printing objects whose size exceeds the scattering mean free path.

[0005] As a new generation of biomanufacturing technology, bulk bioprinting has a higher degree of control over the formed structure, does not require printing support, has a fast forming speed, and can protect cells from shear stress damage. It has many advantages over traditional extrusion bioprinting technology or digital photocuring bioprinting technology. However, it requires high-viscosity ink, and its application is limited to materials with high optical transparency and very little light scattering or absorbance. This technology has not yet been applied to material systems with two-phase solutions that have problems such as refraction and scattering. Therefore, this patent aims to achieve bulk printing in two-phase solutions, and print out products with high fidelity and resolution that can be applied to more fields. Summary of the invention

[0006] In view of the shortcomings existing in the above problems, the present invention provides a method, system, electronic device and storage medium for body-oriented bioprinting scattering correction.

[0007] To achieve the above object, the present invention provides a method for correcting scattering in bulk bioprinting, comprising:

[0008] The projector projects 100 different light patterns onto the square cuvette, obtaining 100 side window images of the projections of different light patterns;

[0009] Binarizing the side window image into grayscale;

[0010] Extract the gray value of each column of each side window image and calculate the one-dimensional Fourier transform result;

[0011] The 100 one-dimensional Fourier transform results are added and summed to obtain an average spectrum;

[0012] Calculating the correction required to keep the spatial frequency uniform at all depths based on the average spectrum, dividing each column of the average matrix by the first column of the average matrix to obtain a corrected spectrum;

[0013] Transforming each column of the modified spectrum into polar coordinates to obtain a two-dimensional mask;

[0014] Each of the two-dimensional masks is dot-multiplied with the two-dimensional Fourier transform of the target cube slice to obtain a dot-multiplication result;

[0015] Performing a two-dimensional inverse Fourier transform on the point product result to obtain correction results of different depths;

[0016] The correction result is effectively filtered, and the filtering results at different depths are superimposed to obtain the scattering correction result of the target cube slice.

[0017] Preferably, the printing test tube is arranged on the rotating platform, the printing test tube is located in the cuvette, printing ink is arranged in the printing test tube, a refractive index matching liquid is arranged in the cuvette, and the projector illuminates the printing test tube.

[0018] Preferably, the grayscale binarization formula is:

[0019] gray_matrix{I n}=im2gray(images{I n})

[0020] n∈[1,100];

[0021] Among them, I represents the input 100 light pattern matrices, I is a q-row and l-column matrix, and the column matrix is ​​denoted as hm(x) (where x = 0, 1, 2...q-1), m∈[1,l].

[0022] Preferably, the grayscale value of each column of each side window image is extracted and the one-dimensional Fourier transform result is calculated, and the formula is:

[0023] fft_result_m=fft(I n {h m (x)});

[0024]

[0025] Preferably, the 100 one-dimensional Fourier transform results are added and summed to obtain an average spectrum, and the formula is:

[0026]

[0027] mean_fft(x)_m=sum_fft(x)_m / n.

[0028] Preferably, the calculation formula of the corrected spectrum is:

[0029] filter(m)=mean_fft_1. / mean_fft_m;

[0030] Filter=filter[1,2,......l].

[0031] Preferably, the two-dimensional mask is recorded as a matrix Bm: q rows and q columns, the target cube slice is recorded as a matrix C(x, y): q rows and q columns, and the calculation formula of the dot product result is:

[0032] result_m(u,v)=FFT2_C.*B m ;

[0033]

[0034] The calculation formula of the correction result is:

[0035]

[0036] The present application also provides a scatter correction system for body-phase bioprinting, comprising:

[0037] A projection module, which is used for a projector to project 100 different light modes onto a square cuvette to obtain 100 side window images of projections of different light modes;

[0038] A grayscale module, used for performing grayscale binarization on the side window image;

[0039] A calculation module, used for extracting grayscale values ​​from each column of each side window image and calculating a one-dimensional Fourier transform result;

[0040] An averaging module, used for adding and summing the 100 one-dimensional Fourier transform results correspondingly to obtain an average spectrum;

[0041] A correction module, used for calculating the correction required to keep the spatial frequency uniform at all depths according to the average spectrum, and dividing each column of the average matrix by the first column of the average matrix to obtain a corrected spectrum;

[0042] A transformation module, used for transforming each column of the modified spectrum into polar coordinates to obtain a one- or two-dimensional mask;

[0043] A point multiplication module, used for performing point multiplication of each of the two-dimensional masks and the two-dimensional Fourier transform of the target cube slice to obtain a point multiplication result;

[0044] A result module, used for performing a two-dimensional inverse Fourier transform on the point product result to obtain correction results of different depths;

[0045] The scattering module is used to effectively filter the correction result and superimpose the filtering results of different depths to obtain the scattering correction result of the target cube slice.

[0046] The present invention also provides an electronic device, comprising at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, the processing unit executes the above method.

[0047] The present invention also provides a storage medium storing a computer program executable by an electronic device, and when the program runs on the electronic device, the electronic device executes the above method.

[0048] Compared with the prior art, the present invention has the following beneficial effects:

[0049] The present invention discloses the scattering characteristics of printing ink in bulk printing, improves the scattering correction rotation tomography algorithm, and realizes bulk biological printing of separation materials with high fidelity. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 It is a flow chart of the scatter correction method for volume bioprinting of the present invention;

[0051] Figure 2 It is a diagram of the hardware device of the body-oriented body-phase bioprinting of the present invention;

[0052] Figure 3 is an uncorrected projection diagram in a certain embodiment of the present invention;

[0053] Figure 4 yes Figure 3 Corresponding corrected projection diagram;

[0054] Figure 5 yes Figure 3 Corresponding to the uncorrected slice image;

[0055] Figure 6 yes Figure 5 Corresponding corrected slice diagram;

[0056] Figure 7 yes Figure 3 Corresponding unedited printed body figure;

[0057] Figure 8 yes Figure 7 Corresponding to the revised printed figure. DETAILED DESCRIPTION

[0058] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0059] Reference Figure 1 The present invention provides a method for correcting scattering of bulk bioprinting, comprising:

[0060] The projector projects 100 different light patterns onto the square cuvette, obtaining 100 side window images of the projections of different light patterns;

[0061] Binarize the side window image into grayscale;

[0062] Specifically, the grayscale binarization formula is:

[0063]

[0064] Among them, I represents the input 100 light pattern matrices, I is a q-row and l-column matrix, and the column matrix is ​​denoted as hm(x) (where x = 0, 1, 2...q-1), m∈[1,l].

[0065] Extract the gray value of each column of each side window image and calculate the one-dimensional Fourier transform result;

[0066] Specifically, the formula is:

[0067] fft_result_m=fft(I n {h m (x)});

[0068]

[0069] Add and sum the 100 one-dimensional Fourier transform results accordingly and obtain the average spectrum;

[0070] Specifically, the formula is:

[0071] The formula for the one-dimensional Fourier summation of m columns is:

[0072] Find the average spectrum of m columns of one-dimensional Fourier transform: mean_fft(x)_m=sum_fft(x)_m / n.

[0073] Calculate the correction required to keep the spatial frequency uniform across all depths based on the average spectrum, and divide each column of the average matrix by the first column of the average matrix to obtain the corrected spectrum;

[0074] Specifically, the formula is:

[0075] Corrected spectrum of column m: fler(m)=mean_fft_1. / mean_fft_m;

[0076] Corrected spectrum matrix l columns q rows: Fliter = flite[1,2,.......l].

[0077] Each column of the corrected spectrum is transformed into polar coordinates to obtain a one- or two-dimensional mask, which is recorded as a matrix Bm: q rows and q columns;

[0078] Each two-dimensional mask is dot-multiplied with the two-dimensional Fourier transform of the target cube slice to obtain a dot-multiplication result;

[0079] Specifically, the cube slice matrix is ​​recorded as matrix C(x,y): q rows and q columns;

[0080] The two-dimensional Fourier transform of the matrix C:

[0081]

[0082] The result of the masked dot product of m columns: result_m(u,v)=FFT2_C.*B m ; result_m is a matrix with q rows and q columns.

[0083] Perform a two-dimensional inverse Fourier transform on the point product result to obtain correction results at different depths;

[0084] Specifically, the formula for the correction result of column m is:

[0085]

[0086] The correction results are effectively filtered, and the filtering results at different depths are superimposed to obtain the scattering correction results of the target cube slice.

[0087] Reference Figure 2 The printing test tube is arranged on a rotating platform, the printing test tube is located in a cuvette, printing ink is arranged in the printing test tube, a refractive index matching liquid is arranged in the cuvette, and a projector illuminates the printing test tube.

[0088] In this embodiment, in terms of hardware equipment selection, a 520nm green light DLP projector is selected as the projection source, and a Wintich PRO4500 series optical engine is used. It is based on the DLPLightCrafter4500 design, uses a 0.45-inch DMD chip, provides HDMI, USB, I2C and synchronous trigger interfaces, and uses a microlens array for light uniformity in the optical path, with an optical uniformity of 95%. The DLP4500 has a central wavelength of 520nm, an output light power of up to 500mW, a lens working distance of 92mm, a projection accuracy of 50μm, a projection surface size of up to 65.6mm×41mm, and a distortion degree of <0.8%. The projection beam is vertically incident on a customized high-transmittance square cuvette, which contains a refractive index matching liquid with a low absorption coefficient for 520nm laser and an adjustable refractive index, and its refractive index is consistent with that of the printing ink. The printing ink (the aqueous two-phase emulsion system material used in this application) is placed in a cylindrical test tube and immersed in the refractive index matching liquid, and is also suspended and fixed under the Zhuoli Hanguan adjustable speed rotating motor. For this type of linear displacement slide, a corresponding series of stepper motor control boxes are matched. The control box can accurately control the speed and acceleration to achieve precise control and adjustment of the printing process. Zhuoli Hanguang adjustable speed rotary motor control box adopts advanced control technology and stable hardware design, has high stability and reliability, and can run stably for a long time. In addition, it also has a rich communication protocol instruction set, supports RS232 and USB communication interfaces, can be connected with other devices or systems, realize data sharing and control interaction, and has good developability.

[0089] Further, in terms of software algorithms, in volume printing, the light projection pattern is in principle determined only according to the 3D shape of the printed object. The target cube STL model to be printed is voxelized into a 3D binary matrix, where "1" indicates the presence of light projection information and "0" indicates the absence of light projection information. The matrix also represents the target light dose required in the printing ink to polymerize it into the desired geometry. A series of dose projections at multiple angles are calculated based on the Radon transform developed from computed tomography. The pattern is obtained using a filtered back-projection algorithm, followed by optimization subject to positivity constraints. This forward model assumes that light is neither attenuated nor distorted during propagation, which is no longer valid in turbid materials. In order to improve the printability of non-transparent materials using the same equipment, it is crucial to optimize the projection pattern set, which is the key to improving the printability of volume printing in non-transparent resins.

[0090] In this embodiment, the pore-forming mechanism of the ATPE system is applied. The material system selected in this application is: a microporous modified bio-ink of a two-phase aqueous emulsion initiated by Ru / Sps with gelatin-methacryloyl (GelMA) as a continuous phase and PEO as a pore-forming agent. The GelMA / PEO two-phase aqueous emulsion ink formula is stabilized by adding rhamnolipid as a biosurfactant.

[0091] Printing process:

[0092] First, generate the projection graphics, input the target cube STL file to be printed, calculate a series of dose projections at multiple angles based on the Radon transform developed by computer tomography technology, and use the filtered back-projection algorithm to obtain the projection graphics. Put the configured double aqueous emulsion ink into a cylindrical test tube, connect the DLP projector and the rotating motor to the computer through the serial port and perform synchronous control. The switching speed of the exposure graphics matches the motor speed. The generated projection pattern is projected onto the double aqueous emulsion ink in the cylindrical test tube through the DLP projector, and the exposure time is controlled to achieve volume printing based on calculated axial exposure.

[0093] In this embodiment, the scatter correction process is:

[0094] Based on the volume printing technology of computational axial exposure, a set of incoherent light patterns are projected through a two-phase solution using the DLP projector of the volume printer. The light pattern selected for scattering characterization is along the x-axis (the x-axis is consistent with the projection direction) and the projection pattern is narrowed as much as possible to improve optical sectioning and increase the contrast of the camera image. The projected pattern falls near the edge of the test tube. The side view image of the side of the square cuvette is captured by a camera with the optical axis orthogonal (z-axis direction). Scattering causes the ballistic light to decrease exponentially with depth, and the light pattern propagating along it becomes blurred. In the frequency domain, the light pattern becomes increasingly blurred with increasing depth, and the scattering material acts as a low-pass filter. Features with high spatial frequencies in the projection decay faster than those with low spatial frequencies. In order to accurately characterize the transmission of all spatial frequencies, the projector projects 100 different light patterns onto the square cuvette. The average transmission of the spatial frequencies of different patterns highlights the strong attenuation of high spatial frequencies when light penetrates the material. 100 side window images with different light modes are taken. In order to alleviate the uneven effect of scattering in the frequency domain space, the amplitude of the high-frequency component should be enhanced. According to the results of the side window image, a correction mask can be calculated to maintain the amplitude of the frequency component on average at all depths z. In the experiment, this correction mask is obtained by dividing the incident average spectrum at z = 0mm by each spectrum taken at different depths. The specific calculation steps are as follows:

[0095] Read side window images: Read 100 side window images and store them in an array.

[0096] Image preprocessing: Convert the image into grayscale and perform mean filtering and smoothing.

[0097] Frequency domain analysis: Extract the grayscale value of each column of each image and calculate its one-dimensional Fourier transform (FFT), move the Fourier transform result to the center of the frequency domain, and calculate the frequency domain features (modulus value).

[0098] Calculate the average: The one-dimensional Fourier results obtained for each column of each image are summed up and the average spectrum is obtained.

[0099] Calculate the corrected spectrum: According to the obtained average spectrum, calculate the correction required to keep the spatial frequency uniform at all depths, and divide each column of the average matrix by the first column of the average matrix to obtain the corrected spectrum. Each column of the corrected spectrum is transformed into a two-dimensional mask by polar coordinates. The row vector of the correction matrix represents the depth, that is, a two-dimensional mask is obtained for different depths. Each two-dimensional mask is dot-multiplied with the two-dimensional Fourier transform of the slice of the target cube, and then the dot-multiplication result is subjected to a two-dimensional inverse Fourier transform to obtain the corrected results at different depths. After that, a filtering that is only effective for a given depth (Z axis) is performed each time, and the filtering results at different depths are superimposed to obtain the scattering correction result of the slice of the target cube. After that, each slice of the target cube is corrected using a loop, and the corrected energy deposition is voxelized and transformed to the projection domain by Radon. The transformation result is high-pass filtered, and the error of the inverse Radon transform result is reduced by combining accurate threshold calculation and high and low threshold erosion, and the optimized axial tomography result is output to achieve high-fidelity bulk bioprinting. Then continue the operation according to the above printing process.

[0100] Overall, the code extracts the frequency characteristics of the image through frequency domain analysis (Fourier transform), especially analyzes the spectrum of the image column, and obtains the overall frequency domain characteristics of the image through averaging method.

[0101] Reference Figure 3-8 Based on the strong scattering effect of the phase separation characteristics of the aqueous two-phase emulsion system material on the projection light in this application, the experimental shooting of the scattering distribution data of the characteristic structured projection in the two-phase printing ink, the analysis of its spatial spectrum characteristics and reverse compensation, and the construction of the scattering-corrected light energy deposition model. The results of the scattering correction algorithm are compared. The projected pattern is a cube stl file. The scattering correction strengthens the intensity around the pattern. After comparison, the shape of the uncorrected print body is obviously distorted, and the outline of the print body after scattering correction is obvious, and the fidelity is enhanced.

[0102] The present application also provides a scatter correction system for body-phase bioprinting, comprising:

[0103] A projection module, which is used for a projector to project 100 different light modes onto a square cuvette to obtain 100 side window images of projections of different light modes;

[0104] Grayscale module, used to convert the side window image into grayscale binarization;

[0105] A calculation module, used for extracting grayscale values ​​from each column of each side window image and calculating a one-dimensional Fourier transform result;

[0106] The averaging module is used to add and sum 100 one-dimensional Fourier transform results and obtain the average spectrum;

[0107] A correction module, for calculating the correction required to keep the spatial frequency uniform at all depths based on the average spectrum, and dividing each column of the average matrix by the first column of the average matrix to obtain a corrected spectrum;

[0108] A transformation module, used for transforming each column of the corrected spectrum into polar coordinates to obtain a one- or two-dimensional mask;

[0109] A point multiplication module is used to perform point multiplication of each two-dimensional mask and the two-dimensional Fourier transform of the target cube slice to obtain a point multiplication result;

[0110] The result module is used to perform a two-dimensional inverse Fourier transform on the point multiplication result to obtain correction results of different depths;

[0111] The scattering module is used to effectively filter the correction results and superimpose the filtering results of different depths to obtain the scattering correction results of the target cube slice.

[0112] The present invention also provides an electronic device, comprising at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, the processing unit executes the above method.

[0113] The present invention also provides a storage medium storing a computer program executable by an electronic device. When the program runs on the electronic device, the electronic device executes the above method.

[0114] The above are only preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for scatter correction of bulk bioprinting, characterized in that: include: The projector projects 100 different light patterns onto the square cuvette, obtaining 100 side window images of the projections of different light patterns; Binarizing the side window image into grayscale; Extract the gray value of each column of each side window image and calculate the one-dimensional Fourier transform result; The 100 one-dimensional Fourier transform results are added and summed to obtain an average spectrum; Calculating the correction required to keep the spatial frequency uniform at all depths based on the average spectrum, dividing each column of the average matrix by the first column of the average matrix to obtain a corrected spectrum; Transforming each column of the modified spectrum into polar coordinates to obtain a two-dimensional mask; Each of the two-dimensional masks is dot-multiplied with the two-dimensional Fourier transform of the target cube slice to obtain a dot-multiplication result; Performing a two-dimensional inverse Fourier transform on the point product result to obtain correction results of different depths; The correction result is effectively filtered, and the filtering results at different depths are superimposed to obtain the scattering correction result of the target cube slice.

2. The method for scatter correction of volume bioprinting according to claim 1, characterized in that: The printing test tube is arranged on the rotating platform, the printing test tube is located in the cuvette, printing ink is arranged in the printing test tube, a refractive index matching liquid is arranged in the cuvette, and the projector illuminates the printing test tube.

3. The method for scatter correction of volume bioprinting according to claim 2, characterized in that: The grayscale binarization formula is: Among them, I represents the input 100 light pattern matrices, I is a q-row and l-column matrix, and the column matrix is ​​denoted as hm(x) (where x = 0, 1, 2...q-1), m∈[1,l].

4. The method for scatter correction of volume bioprinting according to claim 3, characterized in that: The grayscale value of each column of each side window image is extracted and the one-dimensional Fourier transform result is calculated. The formula is: fft_result_m=fft(I n {h m (x)}); 5. The method for scatter correction of volume bioprinting according to claim 4, characterized in that: The 100 one-dimensional Fourier transform results are added and summed to obtain the average spectrum, and the formula is: mean_fft(x)_m=sum_fft(x)_m / n.

6. The method for scatter correction of volume bioprinting according to claim 5, characterized in that: The calculation formula of the corrected spectrum is: filter(m)=mean_fft_1. / mean_fft_m; Filter=filter[1,2,......l].

7. The method for scatter correction of volume bioprinting according to claim 6, characterized in that: The two-dimensional mask is recorded as a matrix Bm: q rows and q columns, and the target cube slice is recorded as a matrix C(x, y): q rows and q columns. The calculation formula of the dot product result is: result_m(u,v)=FFT2_C.*B m ; u=0,1,2...q-1,v=0,1,2...q-1; The calculation formula of the correction result is: x=0,1,2...q-1,y=0,1,2...q-1。 8. A scatter correction system for volume bioprinting, characterized in that: include: A projection module, which is used for a projector to project 100 different light modes onto a square cuvette to obtain 100 side window images of projections of different light modes; A grayscale module, used for performing grayscale binarization on the side window image; A calculation module, used for extracting grayscale values ​​from each column of each side window image and calculating a one-dimensional Fourier transform result; An averaging module, used for adding and summing the 100 one-dimensional Fourier transform results correspondingly to obtain an average spectrum; A correction module, used for calculating the correction required to keep the spatial frequency uniform at all depths according to the average spectrum, and dividing each column of the average matrix by the first column of the average matrix to obtain a corrected spectrum; A transformation module, used for transforming each column of the modified spectrum into polar coordinates to obtain a one- or two-dimensional mask; A point multiplication module, used for performing point multiplication of each of the two-dimensional masks and the two-dimensional Fourier transform of the target cube slice to obtain a point multiplication result; A result module, used for performing a two-dimensional inverse Fourier transform on the point product result to obtain correction results of different depths; The scattering module is used to effectively filter the correction result and superimpose the filtering results of different depths to obtain the scattering correction result of the target cube slice.

9. An electronic device, characterized in that: The method comprises at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the program is executed by the processing unit, the processing unit executes the method according to any one of claims 1 to 7.

10. A storage medium, characterized in that: It stores a computer program executable by an electronic device, and when the program runs on the electronic device, the electronic device executes the method described in any one of claims 1 to 7.

Images