CN106778435A - Feature extracting method based on Image neighborhood structure tensor equation - Google Patents

Feature extracting method based on Image neighborhood structure tensor equation Download PDF

Info

Publication number
CN106778435A
CN106778435A CN201611083692.6A CN201611083692A CN106778435A CN 106778435 A CN106778435 A CN 106778435A CN 201611083692 A CN201611083692 A CN 201611083692A CN 106778435 A CN106778435 A CN 106778435A
Authority
CN
China
Prior art keywords
image
neighborhood
gradient
pixel
structure tensor
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.)
Granted
Application number
CN201611083692.6A
Other languages
Chinese (zh)
Other versions
CN106778435B (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.)
Shanghai Zhengya Dental Technology Co Ltd
Original Assignee
SHANGHAI SMARTEE DENTAL TECHNOLOGY Co Ltd
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 SHANGHAI SMARTEE DENTAL TECHNOLOGY Co Ltd filed Critical SHANGHAI SMARTEE DENTAL TECHNOLOGY Co Ltd
Priority to CN201611083692.6A priority Critical patent/CN106778435B/en
Publication of CN106778435A publication Critical patent/CN106778435A/en
Application granted granted Critical
Publication of CN106778435B publication Critical patent/CN106778435B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06KGRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K7/00Methods or arrangements for sensing record carriers, e.g. for reading patterns
    • G06K7/10Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation
    • G06K7/14Methods or arrangements for sensing record carriers, e.g. for reading patterns by electromagnetic radiation, e.g. optical sensing; by corpuscular radiation using light without selection of wavelength, e.g. sensing reflected white light
    • G06K7/1404Methods for optical code recognition
    • G06K7/1408Methods for optical code recognition the method being specifically adapted for the type of code
    • G06K7/14172D bar codes

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Electromagnetism (AREA)
  • General Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

The invention discloses a kind of feature extracting method based on Image neighborhood structure tensor, carried out according to following steps order:1)Determine the lower boundary of Quick Response Code;2)The extracted region carried out using region-growing method;3)Characteristic value based on structure tensor is calculated;4)Using least square method, three borders of other discrete to obtaining are fitted, and obtain final standard datamatix two-dimension code areas.The present invention calculates the neighborhood border for determining Quick Response Code using the characteristic value of region-growing method and structure tensor, simpler compared with the method for determining each borderline region respectively, and the border for determining is more accurate.Extraction of the present invention suitable for two-dimension code area.

Description

Feature extracting method based on Image neighborhood structure tensor equation
Technical field
The invention belongs to algorithm field, it is related to a kind of Region Feature Extraction method, and in particular to one kind is based on Image neighborhood The feature extracting method of structure tensor equation.
Background technology
The widely used Quick Response Code in Data Matrix Shi world manufacture fields.Data Matrix are the one of Quick Response Code Individual member, was invented with 1989 by U.S. world data company, was widely used in false proof, the pool mark of commodity.It is that one kind can be with The coding of the surface of solids is marked directly on, the coding can automatically be read as usual bar code by corresponding scanning means, Favored by manufacturing industry very much.Current Data Matrix are widely used in product identification, false proof, quality tracing, automatic stored, logistics The system such as management and control.Data Matrix employ the error-correcting code technique of complexity so that the coding has superpower antipollution Ability.Even if coded portion is damaged, reading full detail is not interfered with equally.The print characteristics of Data Matrix cause it into For currently the only support can be marked directly(Print, scribe, photoetching, the mode such as burn into punching press)In product or component surface Coding.Its efficient fault freedom allows it to bear the pollution identified to component surface in manufacture or the process of circulation, because This receives manufacturing welcome very much.For a variety of applications, the Data of diversified forms has been promulgated in the world Matrix symbol standards systems.The minimum dimension of Data Matrix is minimum in current all bar codes, especially particularly suitable In the mark of finding, and it is directly printed on physically.
Data Matrix can be divided into ECC000-140 and ECC200 two types again, and ECC000-140 has various differences The error correcting function of grade, and ECC200 then produces polynomial computation to make mistake correction through Reed-Solomon algorithms Code, its size can on demand be printed as different size, but the error correcting code for using should with dimensional fits, due to its algorithm compared with For easy, and size is more resilient, therefore general more universal with ECC200.Data Matrix code densities are high, and size is small, information Amount is big, domestic also less to DM yards of research to this identification there is provided possible.Data Matrix yards is a kind of matrix type two-dimension Bar code, its maximum feature is exactly density high, and its minimum dimension is the code of minimum in current all bar codes.DM yards can be only 25mm2Area on encode 30 numerals.DM employs the error-correcting code technique of complexity so that the coding has superpower antipollution Ability.Data Matrix can still deposit rational data content because providing minimum and highdensity label, therefore particularly suitable Identified in finding, commodity counterfeit prevention, circuit identifier etc..Due to its outstanding error correcting capability, DM yards oneself turn into South Korea's mobile phone two dimension The mainstream technology of bar code.For QR, DM yards few due to information capacity difference, using simple, is referred to as in the industry " simple Code ", not high to demanding terminal, the mobile phone of 300,000 pixels just can recognize that, it is more the increment based on WAP.Quick Response Code is to hand Machine online brings new entrance, and by scanning all kinds of bar codes, user will soon enter WAP site, carry out fast browsing. Data Matrix symbols look like a chessboard being made up of two kinds of colors of the depth, each formed objects Black or white boxes are referred to as a data unit, and Data Matrix symbols are exactly to be made up of many such data units. Seek border area outer layer to have width is a dead zone for data unit.The border that border area is " chessboard " is sought, is served only for positioning and defining number According to unit-sized, and any coding information is not contained.Coding information is included by the data field for seeking border area encirclement.
Mostly printed using Data Matrix in the prior art, scribed, photoetching, the mode, these modes such as burn into punching press The RM of the Quick Response Code of generation is simpler, because the resolution at its edge is higher, but by the way of 3D printing, directly Generation Quick Response Code, its identification edge is relatively obscured, it is more difficult to realize positioning, therefore, research one kind can be generated for 3D printing method The Quick Response Code method that carries out zone location, have great importance.
The content of the invention
The technical problem to be solved in the present invention, is to provide a kind of feature extraction side based on Image neighborhood structure tensor equation Method, the neighborhood border for determining Quick Response Code is calculated using the characteristic value of region-growing method and structure tensor, and each side is determined more respectively The method in battery limit (BL) domain is simpler, and the border for determining is more accurate.
In order to solve the above technical problems, the technical solution used in the present invention is:
A kind of feature extracting method based on Image neighborhood structure tensor, is carried out according to following steps order:
1)Determine the lower boundary of Quick Response Code;
2)The extracted region carried out using region-growing method;
3)Characteristic value based on structure tensor is calculated;
4)Using least square method, three borders of other discrete to obtaining are fitted, and obtain final standard Datamatix two-dimension code areas.
As a kind of restriction of the invention, described step 1)Carry out in the following order:
a)Calculate image centroid;
To the digital picture of a width two-dimensional discrete, its rank geometric moment is defined as:
;Wherein, p, q=0,1,2..., M, N are the width and height of image;
Then the computing formula of barycenter is:
;Wherein x0It is the abscissa of barycenter, y0It is the ordinate of barycenter;
b)The region of barycenter surrounding pixel is intercepted, and its gray histogram curve is calculated by the row of image array;
c)Obtain estimation region of the row where finding curve global minimum after histogram curve as Quick Response Code lower boundary;
d)Calculate the gradient of minimum region;
Image pixel horizontal directionGradient be:gx
Image pixel vertical directionGradient be:gy =
The gradient magnitude of image is:
e)Gradient magnitude to being calculated carries out non-maxima suppression computing, obtains the local maximum of gradient;
f)The concavo-convex curve matching that will be made up of maximum of gradients using least square method is straight line, is calculated Quick Response Code Lower boundary;Wherein, the object function J of least square method is:J()=min()。
As another restriction of the invention, described step 2)Carried out according to following steps order:
a)Growth starting point of the seed point as region is set;
b)Build similitude and judge criterion, by the region where there is the potting gum to sub-pixel of same nature with seed point In;
c)The pixel that will be newly merged into is used as new seed point;
d)Repeat a)、b)And c)The iteration of algorithm is carried out, until not meeting step b)In conditional pixel be put together.
As the third restriction, described step 3 of the invention)Carried out according to following steps order:
a)Covariance matrix is built using the structure tensor of image;
b)The k gradient in direction of each pixel in image is calculated with gradient operator respectively, k gradient vector square is obtained Battle array, the value of each point is corresponding pixel points Grad in this direction in original image in each matrix;
c)There is the Grad G in its corresponding k direction using each pixel in original image, then can obtain each picture in image The neighborhood of element and the gradient in k direction of each neighborhood, if the matrix that the gradient that L is certain neighborhood of pixels is constituted:
L=, wherein i is the scope of neighborhood of pixels, and k is the gradient direction of neighborhood;
The gradient of the neighborhood of each pixel and the different directions of neighborhood is recycled, the structure tensor based on neighborhood region is constructed Matrix J, then the structure tensor matrix J of each pixel be:
J=LT*L=;Wherein, LTIt is the transposition square of matrix L Battle array;J is symmetrical matrix, using the property of covariance matrix, calculates the covariance matrix, can obtain k and not have correlation Characteristic value, if the pixel is smooth region, k characteristic value is equal0。
d)The all pixels of traversing graph picture, calculate the structure tensor matrix of all pixels point, using structure tensor square respectively The result of calculation of battle array characteristic value can extract the smooth region of image as similarity criterion, as Quick Response Code in itself where Region.
The present invention also has a kind of restriction, and the object function of described least square method is:J()=min[]。
As a result of above-mentioned technical scheme, compared with prior art, acquired technological progress is the present invention:
The present invention calculates the neighborhood border for determining Quick Response Code using the characteristic value of region-growing method and structure tensor, determines more respectively The method of each borderline region is simpler, and the border for determining is more accurate.
Extraction of the present invention suitable for two-dimension code area.
The present invention is described in further detail below in conjunction with Figure of description with specific embodiment.
Brief description of the drawings
Fig. 1 is step 1 in the embodiment of the present invention 1)Extract the schematic diagram of result;
Fig. 2 is step 3 in the embodiment of the present invention 1)Extract the schematic diagram of result;
Fig. 3 is step 4 in the embodiment of the present invention 1)Extract the schematic diagram of result.
Specific embodiment
A kind of feature extracting method based on Image neighborhood structure tensor of embodiment 1
A kind of feature extracting method based on Image neighborhood structure tensor, is carried out according to following steps order:
1)Determine the lower boundary of Quick Response Code;
a)Calculate image centroid;
To the digital picture of a width two-dimensional discrete, its rank geometric moment is defined as:
;Wherein, p, q=0,1,2..., M, N are the width and height of image;
Then the computing formula of barycenter is:
;Wherein x0It is the abscissa of barycenter, y0It is the ordinate of barycenter;
b)The region of barycenter surrounding pixel is intercepted, and its gray histogram curve is calculated by the row of image array;
c)Obtain estimation region of the row where finding curve global minimum after histogram curve as Quick Response Code lower boundary;
d)Calculate the gradient of minimum region;
Image pixel horizontal directionGradient be:gx
Image pixel vertical directionGradient be:gy =
The gradient magnitude of image is:
e)Gradient magnitude to being calculated carries out non-maxima suppression computing, obtains the local maximum of gradient;
f)The concavo-convex curve matching that will be made up of maximum of gradients using least square method is straight line, is calculated Quick Response Code Lower boundary, as shown in Figure 1;Wherein, the object function J of least square method is:J()=min()。
2)The extracted region carried out using region-growing method;
a)Growth starting point of the seed point as region is set;
b)Build similitude and judge criterion, by the region where there is the potting gum to sub-pixel of same nature with seed point In;
c)The pixel that will be newly merged into is used as new seed point;
d)Repeat a)、b)And c)The iteration of algorithm is carried out, until not meeting step b)In conditional pixel be put together.
3)Characteristic value based on structure tensor is calculated;
a)Covariance matrix is built using the structure tensor of image;
b)The k gradient in direction of each pixel in image is calculated with gradient operator respectively, k gradient vector square is obtained Battle array, the value of each point is corresponding pixel points Grad in this direction in original image in each matrix;
c)There is the Grad G in its corresponding k direction using each pixel in original image, then can obtain each picture in image The neighborhood of element and the gradient in k direction of each neighborhood, if the matrix that the gradient that L is certain neighborhood of pixels is constituted:
L=, wherein i is the scope of neighborhood of pixels, and k is the gradient direction of neighborhood;
The gradient of the neighborhood of each pixel and the different directions of neighborhood is recycled, the structure tensor based on neighborhood region is constructed Matrix J, then the structure tensor matrix J of each pixel be:
J=LT*L=;Wherein, LTIt is the transposition square of matrix L Battle array;J is symmetrical matrix, using the property of covariance matrix, calculates the covariance matrix, can obtain k and not have correlation Characteristic value, if the pixel is smooth region, k characteristic value is equal0。
d)The all pixels of traversing graph picture, calculate the structure tensor matrix of all pixels point, using structure tensor square respectively The result of calculation of battle array characteristic value can extract the smooth region of image as similarity criterion, as Quick Response Code in itself where Region, as shown in Figure 2.
4)Using least square method, three borders of other discrete to obtaining are fitted, and obtain final standard Datamatix two-dimension code areas, as shown in Figure 3;Wherein, the object function of described least square method is:J()=min[]。
The present embodiment is applied to the extraction of two-dimension code area, is particularly suited for the Quick Response Code of 3D printing, when specifically used, Quick Response Code identification sign is printed on dentognathic model, but due to the more common printing type Quick Response Code of the two-dimension code area of resin printing Identification is more difficult, therefore in 3D printing dentognathic model, more accurate RM, the present embodiment is needed exist for for identification The extracting method for being provided can effectively determine the identification border of Quick Response Code, facilitate follow-up Quick Response Code to recognize.
The above, is only presently preferred embodiments of the present invention, is not the restriction for making other forms to the present invention, is appointed What those skilled in the art is changed as enlightenment possibly also with above-mentioned technology contents or is modified as equivalent variations Equivalent embodiments.But, it is every without departing from the technology of the present invention design, above example is made according to technical spirit of the invention The simple modification for going out, equivalent variations and remodeling, still fall within the protection domain of the claims in the present invention.

Claims (5)

1. a kind of feature extracting method based on Image neighborhood structure tensor, it is characterised in that it enters according to following steps order OK:
1)Determine the lower boundary of Quick Response Code;
2)The extracted region carried out using region-growing method;
3)Characteristic value based on structure tensor is calculated;
4)Using least square method, three borders of other discrete to obtaining are fitted, and obtain final standard Datamatix two-dimension code areas.
2. the feature extracting method based on Image neighborhood structure tensor according to claim 1, it is characterised in that:Described Step 1)Carry out in the following order:
a)Calculate image centroid;
To the digital picture of a width two-dimensional discrete, its rank geometric moment is defined as:
Wherein, p, q=0,1,2..., M, N are the width and height of image;
Then the computing formula of barycenter is:
Wherein x0It is the abscissa of barycenter, y0It is the ordinate of barycenter;
b)The region of barycenter surrounding pixel is intercepted, and its gray histogram curve is calculated by the row of image array;
c)Obtain estimation region of the row where finding curve global minimum after histogram curve as Quick Response Code lower boundary;
d)Calculate the gradient of minimum region;
Image pixel horizontal directionGradient be:gx
Image pixel vertical directionGradient be:gy =
The gradient magnitude of image is:
e)Gradient magnitude to being calculated carries out non-maxima suppression computing, obtains the local maximum of gradient;
f)The concavo-convex curve matching that will be made up of maximum of gradients using least square method is straight line, is calculated Quick Response Code Lower boundary;Wherein, the object function J of least square method is:J()=min()。
3. the feature extracting method based on Image neighborhood structure tensor according to claim 1, it is characterised in that:Described Step 2)Carried out according to following steps order:
a)Growth starting point of the seed point as region is set;
b)Build similitude and judge criterion, by the region where there is the potting gum to sub-pixel of same nature with seed point In;
c)The pixel that will be newly merged into is used as new seed point;
d)Repeat a)、b)And c)The iteration of algorithm is carried out, until not meeting step b)In conditional pixel be put together.
4. the feature extracting method based on Image neighborhood structure tensor according to claim 1, it is characterised in that:Described Step 3)Carried out according to following steps order:
a)Covariance matrix is built using the structure tensor of image;
b)The k gradient in direction of each pixel in image is calculated with gradient operator respectively, k gradient vector square is obtained Battle array, the value of each point is corresponding pixel points Grad in this direction in original image in each matrix;
c)There is the Grad G in its corresponding k direction using each pixel in original image, then can obtain each picture in image The neighborhood of element and the gradient in k direction of each neighborhood, if the matrix that the gradient that L is certain neighborhood of pixels is constituted:
L=, wherein i is the scope of neighborhood of pixels, and k is the gradient direction of neighborhood;
The gradient of the neighborhood of each pixel and the different directions of neighborhood is recycled, the structure tensor based on neighborhood region is constructed Matrix J, then the structure tensor matrix J of each pixel be:
J=LT*L=;Wherein, LTIt is the transposed matrix of matrix L;J It is symmetrical matrix, using the property of covariance matrix, calculates the covariance matrix, the k feature without correlation can be obtained Value, if the pixel is smooth region, k characteristic value is equal0;
d)The all pixels of traversing graph picture, calculate the structure tensor matrix of all pixels point respectively, special using structure tensor matrix The result of calculation of value indicative can extract the smooth region of image as similarity criterion, as Quick Response Code in itself where area Domain.
5. the feature extracting method based on Image neighborhood structure tensor equation according to claim 1, it is characterised in that:Institute The object function of the least square method stated is:J()=min()。
CN201611083692.6A 2016-11-30 2016-11-30 Feature extraction method based on image neighborhood structure tensor equation Active CN106778435B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611083692.6A CN106778435B (en) 2016-11-30 2016-11-30 Feature extraction method based on image neighborhood structure tensor equation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611083692.6A CN106778435B (en) 2016-11-30 2016-11-30 Feature extraction method based on image neighborhood structure tensor equation

Publications (2)

Publication Number Publication Date
CN106778435A true CN106778435A (en) 2017-05-31
CN106778435B CN106778435B (en) 2020-04-10

Family

ID=58913586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611083692.6A Active CN106778435B (en) 2016-11-30 2016-11-30 Feature extraction method based on image neighborhood structure tensor equation

Country Status (1)

Country Link
CN (1) CN106778435B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549878A (en) * 2018-04-27 2018-09-18 北京华捷艾米科技有限公司 Hand detection method based on depth information and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102346850A (en) * 2011-10-13 2012-02-08 西北工业大学 DataMatrix bar code area positioning method under complex metal background
KR20130012623A (en) * 2011-07-26 2013-02-05 정지환 Teeth attachment for identifying information and its manufacturing method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130012623A (en) * 2011-07-26 2013-02-05 정지환 Teeth attachment for identifying information and its manufacturing method
CN102346850A (en) * 2011-10-13 2012-02-08 西北工业大学 DataMatrix bar code area positioning method under complex metal background

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
周慧: "视屏图像的超分辨率重建技术研究", 《中国优秀硕士学位论文全文数据库》 *
姚林昌: "嵌入式二维条码识别技术的研究与开发", 《中国优秀硕士学位论文全文数据库》 *
杨家红: "结合分水岭与自动种子区域生长的彩色图像分割算法", 《中国图像图形学报》 *
邵云龙: "基于稀疏测度PSF估计的天文图像复原改进算法", 《桂林电子科技大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549878A (en) * 2018-04-27 2018-09-18 北京华捷艾米科技有限公司 Hand detection method based on depth information and system
CN108549878B (en) * 2018-04-27 2020-03-24 北京华捷艾米科技有限公司 Depth information-based hand detection method and system

Also Published As

Publication number Publication date
CN106778435B (en) 2020-04-10

Similar Documents

Publication Publication Date Title
CN106529368B (en) The method of two dimensional code identification region positioning
Paul et al. Remote sensing optical image registration using modified uniform robust SIFT
CN107031033B (en) It is a kind of can 3D printing hollow out two dimensional code model generating method and system
CN109101981B (en) Loop detection method based on global image stripe code in streetscape scene
CN105894069B (en) A kind of generation method and recognition methods of the CRC two dimensional code of vision guided navigation
CN107247985B (en) Coding, positioning and identifying method of two-dimensional code
CN110114781B (en) Method for detecting and identifying remote high density visual indicia
BR112013011943B1 (en) method for identifying a two-dimensional bar code in digital image data of the bar code, non-transitory computer readable medium and apparatus configured to identify a two-dimensional bar code in digital image data of the bar code
CN115293311A (en) Color watermark anti-counterfeiting method and device based on micro-point code
CN106778435A (en) Feature extracting method based on Image neighborhood structure tensor equation
CN106875357A (en) Image in 2 D code processing method
CN112507751A (en) QR code positioning method and system
Ye et al. A local derivative pattern based image forensic framework for seam carving detection
CN107358244A (en) A kind of quick local invariant feature extraction and description method
CN106056575A (en) Image matching method based on object similarity recommended algorithm
CN109711223A (en) A kind of promotion QR code decoding rate method and apparatus
Elwan et al. SAR image matching based on local feature detection and description using convolutional neural network
US20230196707A1 (en) Fiducial patterns
Martynov et al. Aztec core symbol detection method based on connected components extraction and contour signature analysis
US20230267642A1 (en) Fiducial location
Ye et al. A hybrid feature model for seam carving detection
Yi et al. Efficient Localization of Multitype Barcodes in High‐Resolution Images
CN112733748A (en) Voronoi constraint image uniform matching method considering textures
CN111597853B (en) Concrete mark extraction method
CN111191759A (en) Two-dimensional code generation method and positioning and decoding method based on GPU

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201210 Shanghai city Pudong New Area Zu Road No. 2305 A block 122 Building 2 North Gate God's favored one

Applicant after: Shanghai Ya Fang dental Polytron Technologies Inc

Address before: 201210 Shanghai city Pudong New Area Zu Road No. 2305 A block 122 Building 2 North Gate God's favored one

Applicant before: Shanghai Smartee Dental Technology Co., Ltd.

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 201210 Shanghai Pudong New Area Zuchong Road 2305 Tianzhijianzi A Block 122 North Gate 2 Floor

Patentee after: Zhengya Dental Technology (Shanghai) Co.,Ltd.

Address before: 201210 Shanghai Pudong New Area Zuchong Road 2305 Tianzhijianzi A Block 122 North Gate 2 Floor

Patentee before: SHANGHAI SMARTEE DENTI-TECHNOLOGY Co.,Ltd.