CN108765305B - Out-of-focus QR code image blind restoration method - Google Patents

Out-of-focus QR code image blind restoration method Download PDF

Info

Publication number
CN108765305B
CN108765305B CN201810335921.1A CN201810335921A CN108765305B CN 108765305 B CN108765305 B CN 108765305B CN 201810335921 A CN201810335921 A CN 201810335921A CN 108765305 B CN108765305 B CN 108765305B
Authority
CN
China
Prior art keywords
edge
code image
image
focus
spread function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810335921.1A
Other languages
Chinese (zh)
Other versions
CN108765305A (en
Inventor
谭洪舟
谢志勇
陈荣军
谢舜道
林远鑫
朱雄泳
曾衍瀚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Joint Research Institute
Sun Yat Sen University
SYSU CMU Shunde International Joint Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Joint Research Institute, Sun Yat Sen University, SYSU CMU Shunde International Joint Research Institute filed Critical Joint Research Institute
Priority to CN201810335921.1A priority Critical patent/CN108765305B/en
Publication of CN108765305A publication Critical patent/CN108765305A/en
Application granted granted Critical
Publication of CN108765305B publication Critical patent/CN108765305B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06KGRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K19/00Record carriers for use with machines and with at least a part designed to carry digital markings
    • G06K19/06Record carriers for use with machines and with at least a part designed to carry digital markings characterised by the kind of the digital marking, e.g. shape, nature, code
    • G06K19/06009Record carriers for use with machines and with at least a part designed to carry digital markings characterised by the kind of the digital marking, e.g. shape, nature, code with optically detectable marking
    • G06K19/06037Record carriers for use with machines and with at least a part designed to carry digital markings characterised by the kind of the digital marking, e.g. shape, nature, code with optically detectable marking multi-dimensional coding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

The invention discloses a blind restoration method of an out-of-focus QR code image, which comprises the steps of carrying out graying processing and edge detection processing on the QR code image to obtain an edge image, then carrying out processing on the edge image to obtain an edge straight line and a point with the maximum derivative value change rate, calculating out-of-focus radius according to the point with the maximum edge straight line and derivative value change rate to obtain a point spread function, and finally restoring the blurred QR code image according to the point spread function.

Description

Out-of-focus QR code image blind restoration method
Technical Field
The invention relates to the technical field of image processing, in particular to a blind restoration method for an out-of-focus QR code image.
Background
With the rapid development of the internet of things technology, the application of the QR code image is more and more extensive. Because the QR code image is easily affected by various factors to become fuzzy in the processes of acquisition, compression, storage, transmission and reconstruction, and is not beneficial to the wide application of the bar code, the study on the fuzzy QR code image recovery technology and the increase of the bar code image definition have very important theoretical and practical application values.
The defocused blurred image is obtained by convolving a clear QR code image with a point spread function (or a fuzzy kernel) and adding noise, the image restoration process is an image deconvolution process or a deconvolution process, and the image restoration can be divided into two types according to whether the point spread function is known or not: one is a non-blind deconvolution technology of the blurred image, and the point spread function of the blurred image is assumed to be known, and the deconvolution operation is directly carried out on the blurred image; the other is a blind deconvolution technology of the blurred image, and under the condition that a point spread function of the blurred image is unknown, the clear image needs to be restored by combining the prior knowledge of the blurred image, which is an uncertainty problem.
In the practical application process, a point spread function of a QR code acquisition system is unknown, the point spread function and an original clear bar code image need to be estimated at the same time, an effective blind restoration method is that the point spread function is estimated firstly, the blind deconvolution problem is converted into the non-blind deconvolution problem, and therefore the image restoration difficulty is reduced.
Disclosure of Invention
In order to solve the above problems, an object of the present invention is to provide a blind restoration method for an out-of-focus QR code image, which can effectively reduce the amount of computation and the restoration time in the process of restoring a barcode image while ensuring the quality of the restored image.
The technical scheme adopted by the invention for solving the problems is as follows:
a blind restoration method for an out-of-focus QR code image comprises the following steps:
A. carrying out graying processing on an input QR code image;
B. intercepting the QR code image after the graying treatment to obtain an edge image;
C. carrying out edge detection on the edge image to obtain an edge matrix;
D. scanning the edge matrix column by column to obtain the position of an edge straight line;
E. the edge image is subjected to derivation, and a point with the maximum change rate of the derivative value is obtained through calculation;
F. calculating the distance between the edge straight line and the point with the maximum derivative value change rate to obtain an estimated defocus radius;
G. and calculating according to the defocusing radius to obtain a point spread function, and restoring the defocused and blurred QR code image according to the point spread function.
Further, the step a performs a graying process on the input QR code image, where after the graying process is performed on the QR code image, a pixel matrix of the QR code image is expressed as:
Figure BDA0001629243520000031
and a (i, j) is a pixel value with a position (i, j) in the QR code image, 1< i < M, 1< j < N, N is the width of the QR code image, and M is the height of the QR code image.
Further, in the step B, the grayed QR code image is intercepted, and the specific steps are as follows: and preprocessing the QR code image, and intercepting an 1/4 area at the upper left corner of the QR code image. Due to the position characteristics of the QR code positioning frame and the size of the positioning frame of the QR code of different data amounts, by intercepting the upper left corner 1/4 region as a detection object, it is possible to reduce the amount of calculation without affecting the estimation of the blur radius.
Further, in the step C, in the edge detection of the edge image, a Canny detection operator is adopted to perform edge detection on the edge image. Compared with various edge detection operators, the Canny operator has a better effect of detecting the edge images.
Further, in the step D of scanning the edge matrix column by column, the edge matrix is searched by using a search algorithm, and the edge matrix is compared until a column with a continuous numerical value of 1 is obtained for the first time, that is, the position of the edge straight line is obtained.
Further, the step E of deriving the edge image specifically includes the steps of: and carrying out secondary derivation on the edge image, wherein the formula of the primary derivation is as follows:
dx(i,j)=I(i+1,j)-I(i,j);
dy(i,j)=I(i,j+1)-I(i,j);
wherein I is an edge image, I (I, j) is a value at a position (I, j) in the edge image I, and after a first derivative is obtained, the first derivative is derived again according to the following formula:
G(x,y)=dx(i,j)-dy(i,j)。
further, the step E obtains a point with the maximum change rate of the derivative value through calculation, and obtains a point with the maximum change rate of the derivative value according to the second derivative of the edge image I and the normal direction of the edge straight line. In order to efficiently find the point at which the rate of change of the derivative value is the greatest, it is necessary to further derive the first derivative of the edge image.
Further, the step G obtains a point spread function by calculating according to the defocus radius, and brings the defocus radius into the defocus degradation model to obtain the point spread function, where the defocus degradation model is:
Figure BDA0001629243520000041
where h (x, y) is the point spread function and R is the defocus radius.
Furthermore, in the step G, in restoring the defocused blurred QR code image according to the point spread function, the defocused blurred QR code image is restored according to the point spread function and the RL algorithm.
The invention has the beneficial effects that: according to the blind restoration method for the out-of-focus QR code image, the QR code image is subjected to graying processing and edge detection processing to obtain the edge image, then the edge image is processed to obtain the edge straight line and the point with the maximum derivative value change rate, the out-of-focus radius can be calculated according to the edge straight line and the point with the maximum derivative value change rate, further the point spread function can be obtained, finally the blurred QR code image is restored according to the point spread function, the calculation amount is small, and the restoration speed is high.
Drawings
The invention is further illustrated with reference to the following figures and examples.
FIG. 1 is a flow chart of a blind restoration method of an out-of-focus QR code image according to the present invention;
fig. 2 is a QR code image before restoration;
fig. 3 is a QR code image after edge detection;
FIG. 4 is a graph of the relationship after quadratic derivation of the edge image;
fig. 5 is a restored QR code image.
Detailed Description
Referring to fig. 1, the blind restoration method for the out-of-focus QR code image of the present invention includes the following steps:
A. carrying out graying processing on an input QR code image;
because the input QR code images are different from each other and have certain differences in color, state, and other aspects, in order to make the restoration effect of the final QR code image better, it is necessary to perform the graying process on the QR code image, and after the graying process is performed on the QR code image, the pixel matrix is expressed as:
Figure BDA0001629243520000051
and a (i, j) is a pixel value with a position (i, j) in the QR code image, 1< i < M, 1< j < N, N is the width of the QR code image, and M is the height of the QR code image.
B. Intercepting the QR code image after the graying treatment to obtain an edge image;
firstly, because the QR code image is obtained by shooting or scanning in different ways, the obtained QR code image has certain redundant information, so that the grayed QR code image needs to be preprocessed to eliminate redundant irrelevant information in the QR code image and restore part of real useful information, so as to ensure the restoration effect.
After the QR code image is preprocessed, the QR code image needs to be intercepted, and due to the position characteristics of the QR code positioning frame and the sizes of the positioning frames of the QR codes with different data volumes, the estimation of the fuzzy radius is not influenced while the calculation amount is reduced by intercepting the upper left corner 1/4 area as a detection object.
C. Carrying out edge detection on the edge image to obtain an edge matrix;
after comparing the effects of different edge detection operators, the method selects a Canny detection operator with a better effect to carry out edge detection on the edge image, and the Canny aims at finding an optimal edge detection algorithm, and comprises the following three steps: denoising; finding a brightness gradient in the image; edges are tracked in the image.
The Canny detection operator is used as a multi-stage edge detection algorithm and is suitable for different occasions, and parameters of the Canny detection operator are allowed to be adjusted according to different implementation specific requirements so as to identify different edge characteristics, so that the Canny detection operator has a better detection effect compared with other edge detection operators.
D. Scanning the edge matrix column by column to obtain the position of an edge straight line L;
specifically, in the step D, during the column-by-column scanning of the edge matrix, the edge matrix is searched by using a search algorithm, and the edge matrix is compared until a column with a continuous numerical value of 1 is obtained for the first time, that is, the position of the edge straight line L is obtained.
E. The edge image is subjected to derivation, and a point Q with the maximum change rate of the derivative value is obtained through calculation;
specifically, the derivation of the edge image in step E specifically includes: and carrying out secondary derivation on the edge image, wherein the formula of the primary derivation is as follows:
dx(i,j)=I(i+1,j)-I(i,j);
dy(i,j)=I(i,j+1)-I(i,j);
after the first derivative is obtained, in order to efficiently find the point Q with the maximum change rate of the derivative value, the first derivative of the edge image needs to be subjected to derivation again, and the first derivative is subjected to derivation again according to the following formula:
G(x,y)=dx(i,j)-dy(i,j)。
and after the second derivative of the edge image I is obtained, obtaining a point Q with the maximum derivative value change rate according to the second derivative of the edge image I and the normal direction of the edge straight line.
F. Calculating the distance between the edge straight line L and the point Q with the maximum derivative value change rate, namely the distance between the corresponding columns to obtain an estimated defocus radius R;
G. and calculating to obtain a point spread function according to the defocusing radius R, and restoring the defocused and blurred QR code image according to the point spread function.
Specifically, the step G obtains a point spread function by calculating according to the defocus radius, and brings the defocus radius into the defocus degradation model to obtain the point spread function, where the defocus degradation model is:
Figure BDA0001629243520000071
where h (x, y) is the point spread function.
After the point spread function is obtained, restoring the defocused and blurred QR code image according to the point spread function and an RL algorithm.
In order to verify the restoration effect of the invention on the blurred QR code image, a blurred QR code image is firstly input as shown in FIG. 2, and then is subjected to graying processing to ensure thatThen, after a subsequent restoration effect, performing edge detection on the input QR code image to obtain an edge matrix, as shown in fig. 3, which is the QR code image after edge detection, then using a search algorithm to find edge straight lines, where the columns with the first numerical values all 1 in the edge matrix are the columns where the edge straight lines are located, and performing derivation again on the first derivative of the edge image, as shown in the relationship diagram after second derivation shown in fig. 4, where the vertical coordinate is a gradient value (i.e., the result after second derivation), and the horizontal coordinate is the position of the image column, where the vertical coordinate may be regarded as the position where the edge straight line is located, and the horizontal coordinate may be regarded as the normal of the edge straight line, starting to search the point with the largest derivative change rate of the origin from the normal direction of the edge straight line, that is to determine the position of the Q point, and as shown in fig. 4, the horizontal coordinate of the point with the largest derivative change rate is 10, that is, the coordinate of the position of the point Q where the derivative change rate is the largest can be determined to be 10, and the coordinate of the edge straight line is 0 since the edge straight line is set as the starting point, so that the estimated defocus radius of 10 can be obtained, and then the defocus radius is substituted into the defocus degradation model, and the point spread function is calculated
Figure BDA0001629243520000081
Finally, the RL algorithm is combined with the point spread function to restore the blurred QR code image to obtain a restored image shown in fig. 5, and compared with fig. 2 and 5, the restoration method has the advantages that the restoration effect is better, the image definition is restored, and the image quality is ensured.
Specifically, since the size of the defocus radius is determined by the distance between the edge line and the point Q, the specific coordinate where the edge line is located does not need to be obtained, and only the position of the edge line needs to be known, after the position of the edge line is obtained, the position where the edge line is located is taken as the ordinate, the normal direction of the edge line is taken as the abscissa to establish a coordinate system, and then the point with the maximum derivative change rate is found according to the result of the second derivation, since the intervals between the points in the edge image are consistent, the distance from the point with the maximum derivative change rate to the edge line can be determined as the size of the defocus radius R, that is, the image column position on the abscissa in fig. 4 is not the specific position actually in the edge image, but a relative position.
Meanwhile, the restoration method of the invention is compared with other known methods, an image with the size of 512 x 512 is selected as an input, the time spent in the whole restoration process is compared,
Method document 1 Document 2 Document 3 Document 4 The method of the invention
Recovery time 21.58s 20.87s 24.77s 24.81s 0.45s
The time spent by each method is shown in table 1:
TABLE 1 (recovery time data for the invention and other methods)
As can be seen from the data in Table 1, the recovery time spent by the method of the present invention is considerably shorter, 50-60 times faster than other recovery methods.
Therefore, the invention has good restoration effect, relatively short calculation time in the whole restoration process and small calculation amount.
Reference documents:
document 1: pan J, Hu Z, Su Z, et al, Deblurring Text Images via L0-regulated Intensity and Gradient primer [ C ]// Computer Vision and Pattern recognition. IEEE 2014: 2901-.
Document 2: krishnan D, Tay T, Fergus R.Blind reduction using a normalized specific measure [ C ]// Computer Vision and Pattern recognition. IEEE,2011: 233-.
Document 3: perrone D, Favaro P.A clean Picture of Total Variation Blanket reduction [ J ]. IEEE Transactions on Pattern Analysis & Machine understanding, 2016,38(6): 1041-.
Document 4: gao K, Zhu Z, Dou Z, et al. variable exposure Regulation application for blade Kernel Estimation of Remote Sensing Image blade retrieval [ J ]. IEEE Access,2018,6(99): 4352-.
The above description is only a preferred embodiment of the present invention, and the present invention is not limited to the above embodiment, and the present invention shall fall within the protection scope of the present invention as long as the technical effects of the present invention are achieved by the same means.

Claims (7)

1. A blind restoration method for an out-of-focus QR code image is characterized by comprising the following steps: the method comprises the following steps:
A. carrying out graying processing on an input QR code image;
B. intercepting the QR code image after the graying treatment to obtain an edge image;
C. carrying out edge detection on the edge image to obtain an edge matrix;
D. scanning the edge matrix column by column to obtain the position of an edge straight line;
E. the edge image is subjected to derivation, and a point with the maximum change rate of the derivative value is obtained through calculation;
F. calculating the distance between the edge straight line and the point with the maximum derivative value change rate to obtain an estimated defocus radius;
G. calculating according to the defocusing radius to obtain a point spread function, and restoring the defocused and blurred QR code image according to the point spread function;
the step E of deriving the edge image includes the following specific steps: and carrying out secondary derivation on the edge image, wherein the formula of the primary derivation is as follows:
dx(i,j)=I(i+1,j)-I(i,j);
dy(i,j)=I(i,j+1)-I(i,j);
wherein I is an edge image, I (I, j) is a value at a position (I, j) in the edge image I, and after a first derivative is obtained, the first derivative is derived again according to the following formula:
G(x,y)=dx(i,j)-dy(i,j);
and E, calculating to obtain the point with the maximum change rate of the derivative value, and obtaining the point with the maximum change rate of the derivative value according to the second derivative of the edge image I and the normal direction of the edge straight line.
2. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: the step A is to perform graying processing on the input QR code image, wherein after the graying processing is performed on the QR code image, a pixel matrix of the QR code image is expressed as:
Figure FDA0002812465720000021
and a (i, j) is a pixel value with a position (i, j) in the QR code image, i is more than 1 and less than M, j is more than 1 and less than N, N is the width of the QR code image, and M is the height of the QR code image.
3. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: in the step B, the grayed QR code image is intercepted, and the specific steps are as follows: and preprocessing the QR code image, and intercepting an 1/4 area at the upper left corner of the QR code image.
4. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: and C, in the step of edge detection of the edge image, edge detection is carried out on the edge image by adopting a Canny detection operator.
5. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: and D, in the step of scanning the edge matrix row by row, searching the edge matrix by adopting a search algorithm, and comparing until a row with a continuous numerical value of 1 is obtained for the first time, namely the position of the edge straight line.
6. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: and G, calculating according to the defocusing radius to obtain a point spread function, and bringing the defocusing radius into a defocusing degradation model to obtain the point spread function, wherein the defocusing degradation model is as follows:
Figure FDA0002812465720000031
where h (x, y) is the point spread function and R is the defocus radius.
7. The blind restoration method for the out-of-focus QR code image according to claim 1, characterized in that: and G, in the restoration of the defocused and blurred QR code image according to the point spread function, restoring the defocused and blurred QR code image according to the point spread function and the RL algorithm.
CN201810335921.1A 2018-04-16 2018-04-16 Out-of-focus QR code image blind restoration method Active CN108765305B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810335921.1A CN108765305B (en) 2018-04-16 2018-04-16 Out-of-focus QR code image blind restoration method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810335921.1A CN108765305B (en) 2018-04-16 2018-04-16 Out-of-focus QR code image blind restoration method

Publications (2)

Publication Number Publication Date
CN108765305A CN108765305A (en) 2018-11-06
CN108765305B true CN108765305B (en) 2021-04-06

Family

ID=64010878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810335921.1A Active CN108765305B (en) 2018-04-16 2018-04-16 Out-of-focus QR code image blind restoration method

Country Status (1)

Country Link
CN (1) CN108765305B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112200741B (en) * 2020-10-09 2022-05-17 广东技术师范大学 Out-of-focus two-dimensional code image restoration method based on edge prior
CN113743147B (en) * 2021-07-16 2023-08-29 广东技术师范大学 Anti-fuzzy two-dimensional code image based on fuzzy imaging mechanism and blind restoration method thereof

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236650B (en) * 2008-01-24 2010-11-03 上海交通大学 Real time rotating blurred image restoration method based on gradient loading optimization
JP5505007B2 (en) * 2010-03-18 2014-05-28 富士通株式会社 Image processing apparatus, image processing method, and computer program for image processing
JP5401501B2 (en) * 2011-04-27 2014-01-29 京セラドキュメントソリューションズ株式会社 Image processing device
US10282647B2 (en) * 2015-05-05 2019-05-07 Massachusetts Institute Of Technology Substrate pre-scanning for high throughput microscopy
CN107067373B (en) * 2016-12-13 2019-08-20 南京航空航天大学 A kind of gradient minimisation recovery method of binary image based on 0 norm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A low-complexity pre-processing system for restoring low-quality QR code images;David Munoz-Mejias、Ivan Gonzalez-Diaz、Fernando Diaz-de-Maria;《IEEE Transactions on Consumer Electronics》;20110915;全文 *

Also Published As

Publication number Publication date
CN108765305A (en) 2018-11-06

Similar Documents

Publication Publication Date Title
US8090214B2 (en) Method for automatic detection and correction of halo artifacts in images
CN102156996B (en) Image edge detection method
CN101996406A (en) No-reference structural sharpness image quality evaluation method
CN108765305B (en) Out-of-focus QR code image blind restoration method
CN111489337A (en) Method and system for removing false defects through automatic optical detection
CN111489389A (en) Light spot center detection method
CN105447828A (en) Single-viewpoint image deblurring method for carrying out one-dimensional deconvolution along motion blur path
CN104537381A (en) Blurred image identification method based on blurred invariant feature
CN112200741B (en) Out-of-focus two-dimensional code image restoration method based on edge prior
CN112419369A (en) Anti-occlusion real-time target tracking method
CN107886530A (en) A kind of improved image registration algorithm based on SIFT feature
CN115588109B (en) Image template matching method, device, equipment and application
CN116739943A (en) Image smoothing method and target contour extraction method
CN105825480A (en) Image deblurring method based on sapphire surface defect automatic detection
CN116740501A (en) Training method and application of image blurring region restoration compensation model
US20240005454A1 (en) Method of performing l0 smoothing on the basis of deep gradient prior information to improve sharpness
CN116128826A (en) YOLOv 5-based lithium battery defect detection method, system, equipment and storage medium
CN104966271A (en) Image denoising method based on biological vision receptive field mechanism
Chang et al. A hybrid motion deblurring strategy using patch based edge restoration and bilateral filter
Min et al. Out-of-focus image deblurring for mobile display vision inspection
Vignesh et al. Performance and Analysis of Edge detection using FPGA Implementation
CN104966292A (en) Multi-book object segmentation method based on Gabor filter in vertical direction
JP5855210B2 (en) Information processing apparatus, line noise reduction processing method, and program
CN114581443B (en) Image processing method and device, computer equipment and readable storage medium
Karim et al. Improvement of Harris Algorithm Based on Gaussian Scale Space

Legal Events

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

Effective date of registration: 20240531

Address after: 510000 No. 135 West Xingang Road, Guangdong, Guangzhou

Patentee after: SUN YAT-SEN University

Country or region after: China

Address before: No.9, Nanguo East Road, Yunlu community residents committee, Daliang sub district office, Shunde District, Foshan City, Guangdong Province, 528399

Patentee before: FOSHAN SHUNDE SUN YAT-SEN UNIVERSITY Research Institute

Country or region before: China

Patentee before: SYSU-CMU SHUNDE INTERNATIONAL JOINT Research Institute

Patentee before: SUN YAT-SEN University