CN114998271B - A substation high-voltage switchgear secondary warehouse panel image recognition system - Google Patents

A substation high-voltage switchgear secondary warehouse panel image recognition system Download PDF

Info

Publication number
CN114998271B
CN114998271B CN202210651684.6A CN202210651684A CN114998271B CN 114998271 B CN114998271 B CN 114998271B CN 202210651684 A CN202210651684 A CN 202210651684A CN 114998271 B CN114998271 B CN 114998271B
Authority
CN
China
Prior art keywords
image
template
component
picture
voltage switch
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
CN202210651684.6A
Other languages
Chinese (zh)
Other versions
CN114998271A (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.)
Hainan Research Institute Of Zhejiang University
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Hainan Institute of Zhejiang University
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 Zhejiang University ZJU, Hainan Institute of Zhejiang University filed Critical Zhejiang University ZJU
Priority to CN202210651684.6A priority Critical patent/CN114998271B/en
Publication of CN114998271A publication Critical patent/CN114998271A/en
Application granted granted Critical
Publication of CN114998271B publication Critical patent/CN114998271B/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
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/02Affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/28Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/774Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/82Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及变电站高压开关柜图像识别技术领域,尤其是指一种变电站高压开关柜二次仓面板图像识别系统,首先,在每面高压开关柜二次仓面板上张贴一张ID唯一的AprilTag标签,采集二次仓面板图片,并为每面高压开关柜选取1张清晰图片作为处理模板图片;其次,将每张模板图划分为指示灯区域、压板开关区域、表计区域,在切割出的区域图像上用部件矩形框标注每一个部件;然后,检测模板图片和非模板图片中AprilTag标签,获取指示灯和压板开关数据集、数字表数据集;再次训练ResNet50分类网络,获取指示灯和压板开关状态,设计基于模板匹配的自动读数方法,获取数字表读数;最后根据每一个部件预设的正常运行状态下标准阈值,逻辑判断部件是否异常,并及时告警。

Figure 202210651684

The present invention relates to the technical field of substation high-voltage switchgear image recognition technology, in particular to an image recognition system for substation high-voltage switchgear secondary warehouse panels. , collect the panel pictures of the secondary warehouse, and select a clear picture for each high-voltage switchgear as a processing template picture; secondly, divide each template picture into the indicator light area, the pressure plate switch area, and the meter area, and cut out Mark each component with a component rectangle on the region image; then, detect the AprilTag label in the template image and non-template image, obtain the indicator light and pressure plate switch data set, and the digital table data set; train the ResNet50 classification network again to obtain the indicator light and pressure plate For switch status, an automatic reading method based on template matching is designed to obtain digital meter readings; finally, according to the standard threshold value of each component in the normal operating state, it is logically judged whether the component is abnormal, and an alarm is issued in time.

Figure 202210651684

Description

一种变电站高压开关柜二次仓面板图像识别系统An image recognition system for secondary compartment panels of high-voltage switch cabinets in substations

技术领域Technical Field

本发明涉及变电站高压开关柜图像识别技术领域,尤其涉及一种变电站高压开关柜二次仓面板图像识别系统。The present invention relates to the technical field of image recognition of high-voltage switch cabinets in substations, and in particular to an image recognition system for secondary compartment panels of high-voltage switch cabinets in substations.

背景技术Background Art

对于变电站高压开关柜日常巡检,现场工作人员通过查看高压开关柜二次仓面板,记录每一个部件状态,从而获取高压开关柜一次仓内开关设备状态。但目前大多数变电站高压开关柜的巡检主要依赖于人工,通过定期人工巡检抄录每一面高压开关柜上部件状态,一旦开关设备故障无法及时发现,可能带来严重损失。而由于高压开关柜二次仓面板上包含部件种类和数量众多,需要同时识别指示灯状态、压板开关状态和数字表读数,难以通过单一深度学习模型直接实现,且深度学习模型在大量的标注数据基础上训练才有可能得到不错的识别效果。同时,基于深度学习的图像识别算法在实际应用时存在部件漏检和非部件误检为部件的问题,导致识别结果无法和实际位置的部件相对应。因此,仅采用基于深度学习的算法无法实现高压开关柜二次仓面板图像识别的可靠应用,亟需设计一种通用且识别准确率高的图像识别系统,能同时识别不同种类的部件,使其具有推广应用价值。For daily inspection of high-voltage switchgear in substations, on-site staff check the secondary compartment panel of the high-voltage switchgear and record the status of each component, so as to obtain the status of the switchgear in the primary compartment of the high-voltage switchgear. However, at present, the inspection of high-voltage switchgear in most substations mainly relies on manual work. The status of components on each side of the high-voltage switchgear is recorded through regular manual inspection. Once the switchgear failure cannot be discovered in time, it may cause serious losses. However, since the secondary compartment panel of the high-voltage switchgear contains a large number of components, it is necessary to simultaneously identify the indicator light status, the pressure plate switch status and the digital meter reading, which is difficult to achieve directly through a single deep learning model. In addition, the deep learning model can only obtain good recognition effects after training on a large amount of labeled data. At the same time, the image recognition algorithm based on deep learning has the problems of missed detection of components and misdetection of non-components as components in actual application, resulting in the recognition results not corresponding to the components in the actual position. Therefore, it is impossible to achieve reliable application of image recognition of the secondary compartment panel of the high-voltage switchgear by only using algorithms based on deep learning. It is urgent to design a universal image recognition system with high recognition accuracy, which can simultaneously identify different types of components and make it valuable for promotion and application.

发明内容Summary of the invention

本发明的目的是为了解决现有技术中存在的缺点,同时解决高压开关柜依赖人工巡检巡检,解放工人双手,从而提出的一种变电站高压开关柜二次仓面板图像识别系统。The purpose of the present invention is to solve the shortcomings of the prior art and solve the problem that high-voltage switch cabinets rely on manual inspections and free workers' hands, thereby proposing an image recognition system for secondary compartment panels of high-voltage switch cabinets in substations.

为了实现上述目的,本发明采用了如下技术方案:一种变电站高压开关柜二次仓面板图像识别系统,包括以下步骤:In order to achieve the above object, the present invention adopts the following technical solution: a substation high-voltage switch cabinet secondary compartment panel image recognition system, comprising the following steps:

1)在每一面高压开关柜的二次仓面板水平中轴线上,张贴一张ID唯一的AprilTag标签,并记录高压开关柜名称和标签ID一一对应列表;1) On the horizontal center axis of the secondary compartment panel of each high-voltage switch cabinet, post an AprilTag label with a unique ID, and record the one-to-one correspondence between the high-voltage switch cabinet name and the label ID;

2)采集高压开关柜样本图片,构成样本库,每张样本图片包含变电站高压开关柜二次仓面板和AprilTag标签;2) Collect sample images of high-voltage switchgear to form a sample library. Each sample image contains the secondary compartment panel and AprilTag label of the substation high-voltage switchgear;

3)从样本库中,为每面高压开关柜选取1张清晰图片作为模板图;3) Select one clear picture from the sample library as a template picture for each high-voltage switchgear cabinet;

4)调用AprilTag算法检测模板图片中所有标签ID和标签四个角点像素坐标,计算标签中心点水平像素坐标和样本图片水平中轴线像素坐标差值,差值最小的即为当前高压开关柜上对应的标签;4) Call the AprilTag algorithm to detect all tag IDs and the pixel coordinates of the four corner points of the tags in the template image, and calculate the difference between the horizontal pixel coordinates of the center point of the tag and the horizontal central axis pixel coordinates of the sample image. The one with the smallest difference is the corresponding tag on the current high-voltage switchgear;

5)按照步骤4)中获取的对应标签ID,查找步骤1)中的名称和ID列表,获取模板图片所采集的高压开关柜名称,并存储模板图名称、标签ID和四个角点像素坐标到数据库中;5) According to the corresponding tag ID obtained in step 4), search the name and ID list in step 1), obtain the name of the high-voltage switchgear collected by the template image, and store the template image name, tag ID and four corner point pixel coordinates in the database;

6)依据高压开关柜二次仓面板上部件类型和分布特征,将每一张模板图划分指示灯区域、压板开关区域、数字表区域,切割区域图像并命名;6) According to the component types and distribution characteristics on the secondary compartment panel of the high-voltage switch cabinet, each template image is divided into the indicator light area, the pressure plate switch area, and the digital meter area, and the regional images are cut and named;

7)在区域图像中标注部件,命名部件名称,按照部件标注框逐一切割部件图像,并设定每个部件的标准状态;7) Mark the components in the area image, name the components, cut the component images one by one according to the component marking frame, and set the standard state of each component;

8)从样本库中依次读取非模板图片,重复步骤4)获取每张非模板图片中标签ID和四个角点像素坐标;8) Read non-template images from the sample library in sequence, and repeat step 4) to obtain the label ID and four corner pixel coordinates in each non-template image;

9)根据步骤8)中标签角点像素坐标,结合步骤4)中模板图中标签角点像素坐标,计算模板图片和非模板图片之间的仿射变换矩阵;9) Calculate the affine transformation matrix between the template image and the non-template image based on the pixel coordinates of the label corner points in step 8) and the pixel coordinates of the label corner points in the template image in step 4);

10)根据仿射变换矩阵,将模板图片中的区域标注框映射到非模板图片中,切割区域图像并命名保存,再将模板图中部件标注框全部映射到非模板图片中,实现非模板图部件定位;10) According to the affine transformation matrix, the area annotation box in the template image is mapped to the non-template image, the area image is cut and named and saved, and then all the component annotation boxes in the template image are mapped to the non-template image to realize the positioning of the non-template image components;

11)按照定位的部件标注框,切割样本库中每张非模板图片的部件图像,并命名每张部件图像,获得部件图像数据集;11) According to the located component annotation box, cut the component image of each non-template image in the sample library, name each component image, and obtain the component image dataset;

12)将数据集按照部件图像名称分为指示灯数据集、压板开关数据集、数字表数据集;12) Divide the data set into an indicator light data set, a pressure plate switch data set, and a digital table data set according to the component image names;

13)基于步骤12)获取的指示灯和压板开关数据集,进行数据增广,将数据集划分为训练集、验证集和测试集,在算法服务器上导入数据集和ResNet50分类模型,训练基于ResNet50的指示灯和压板开关分类算法模型,获取指示灯“亮”或“灭”状态和压板开关“开”或“合”状态;13) Based on the indicator light and pressure plate switch data set obtained in step 12), data augmentation is performed, the data set is divided into a training set, a validation set, and a test set, the data set and the ResNet50 classification model are imported on the algorithm server, and the indicator light and pressure plate switch classification algorithm model based on ResNet50 is trained to obtain the "on" or "off" state of the indicator light and the "open" or "closed" state of the pressure plate switch;

14)基于步骤12)获取的数字表数据集,设计基于模板匹配的自动读数算法,制作0~9的LED数字二值化模板图,并逐一计算新图片中LED数字表单个数字二值化图像和数字二值化模板图的像素匹配度,对匹配度排序得到匹配度最高的单个数字读数,对数字位置按照水平坐标排序,按照位数计算数字表实际读数;14) Based on the digital table data set obtained in step 12), an automatic reading algorithm based on template matching is designed to produce a binary template image of LED digital numbers from 0 to 9, and the pixel matching degree between the single binary image of the LED digital table and the digital binary template image in the new image is calculated one by one, and the matching degree is sorted to obtain the single digital reading with the highest matching degree, and the digital positions are sorted according to the horizontal coordinates, and the actual reading of the digital table is calculated according to the number of digits;

15)将新采集的高压开关柜图片,输入到高压开关柜图像识别系统中,重复步骤4)和步骤5),获取图片中标签四个角点像素坐标,重复步骤6)获取区域图像;15) Input the newly acquired high-voltage switchgear image into the high-voltage switchgear image recognition system, repeat steps 4) and 5) to obtain the pixel coordinates of the four corner points of the label in the image, and repeat step 6) to obtain the regional image;

16)计算新采集图片和模板图片的仿射变换矩阵,将模板图中所有部件标注框映射到新采集图片上,切割部件图像并命名,按照部件图像名称中的部件类型,调用基于ResNet50的分类算法模型或基于模板匹配的自动读数算法,获取识别结果;16) Calculate the affine transformation matrix of the newly collected image and the template image, map all component annotation boxes in the template image to the newly collected image, cut the component image and name it, and call the classification algorithm model based on ResNet50 or the automatic reading algorithm based on template matching according to the component type in the component image name to obtain the recognition result;

17)读取每一个指示灯、压板开关、数字表预设的标准状态,逻辑判断部件状态是否异常,并及时告警。17) Read the preset standard status of each indicator light, pressure plate switch, and digital meter, logically determine whether the component status is abnormal, and issue an alarm in time.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤1)的标签ID,具体是除0以外的自然数。The tag ID in step 1) is specifically a natural number other than 0.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤6)的二次仓面板上部件类型,主要包含LED指示灯、压板开关、LED数字表,依次命名为LED、SWITCH、METER;所述步骤6)的划分区域,采用区域矩形框划分,每面高压开关柜二次仓面板图像最多划分3个区域,对区域依次命名为“*P-PSx”,其中“*P”表示高压开关柜简称,“PSx”表示区域简称,“x”取值为1、2、3,取值不重复;所述区域矩形框,要求按照部件类型,框到所有的待识别部件,并以区域名称为key值,区域矩形框的左上和右下角点像素坐标为value,存储到Redis数据库中;所述步骤6)的切割区域图像,按照区域矩形框切割图片,切割后按照区域名称对区域图像命名为“*P-PSx-time.jpg”,“time”为计算机当前毫秒数,并将区域图像保存在本地服务器中;所述本地服务器,操作系统为LinuxUbuntu。The component types on the secondary warehouse panel in step 6) mainly include LED indicator lights, pressure plate switches, and LED digital meters, which are named LED, SWITCH, and METER in sequence; the divided areas in step 6) are divided by regional rectangular frames, and each high-voltage switch cabinet secondary warehouse panel image is divided into up to 3 areas, and the areas are named "*P-PSx" in sequence, where "*P" represents the abbreviation of the high-voltage switch cabinet, "PSx" represents the abbreviation of the area, and "x" takes the values of 1, 2, and 3, and the values are not repeated; the regional rectangular frame is required to frame all the components to be identified according to the component type, and use the area name as the key value, and the pixel coordinates of the upper left and lower right corners of the regional rectangular frame as the value, which are stored in the Redis database; the cutting area image in step 6) cuts the picture according to the regional rectangular frame, and after cutting, the regional image is named "*P-PSx-time.jpg" according to the regional name, "time" is the current number of milliseconds on the computer, and the regional image is saved in the local server; the operating system of the local server is LinuxUbuntu.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤7)的标注部件,采用部件矩形框标注,对每一个部件命名,并设定高压开关柜正常运行状态下每一个部件的标准状态,存储部件标注框角点坐标和部件名称到数据库中,所述部件命名,结合区域名称、部件类型和编号,三种部件具体命名规则为“*P-PSx-LEDy”、“*P-PSx-SWITCHy”、“*P-PSx-METERy”,“y”表示非零自然数,在单个区域图像内,y取值唯一,所述存储部件标注框角点坐标和部件名称,以部件名称为key值,区域图像中部件标注框和部件标准状态为value,存储到Redis数据库;所述步骤7)的切割部件图像,对切割出每一个部件图像命名,根据部件名称将部件图像命名为“*P-PSx-LEDy-time.jpg”或“*P-PSx-SWITCHy-time.jpg”或“*P-PSx-METERy-time.jpg”,并将部件图像保存在本地服务器中。The labeled components of step 7) are labeled with component rectangular frames, each component is named, and the standard state of each component under the normal operating state of the high-voltage switch cabinet is set, and the coordinates of the corner points of the component labeling frame and the component name are stored in the database. The component naming is combined with the area name, component type and number. The specific naming rules of the three components are "*P-PSx-LEDy", "*P-PSx-SWITCHy", and "*P-PSx-METERy". "y" represents a non-zero natural number. In a single area image, y has a unique value. The storage component labeling frame corner point coordinates and component name, with the component name as the key value, the component labeling frame and the component standard state in the area image as the value, are stored in the Redis database; the cutting component image of step 7) is named for each cut component image, and the component image is named "*P-PSx-LEDy-time.jpg" or "*P-PSx-SWITCHy-time.jpg" or "*P-PSx-METERy-time.jpg" according to the component name, and the component image is saved in the local server.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤9)的计算仿射变换矩阵,具体为:The calculation of the affine transformation matrix in step 9) is specifically:

Figure BDA0003687934340000051
Figure BDA0003687934340000051

其中,M2×3为仿射变换矩阵,表示通过线性变化将模板图片中的像素齐次坐标(x0,y0,1)变换到非模板图片中对应位置像素齐次坐标(x1,y1,1),通过模板图和非模板图片中对应三对像素坐标可求解得到M;Among them, M 2×3 is the affine transformation matrix, which means that the pixel homogeneous coordinates (x 0 ,y 0 ,1) in the template image are transformed to the pixel homogeneous coordinates (x 1 ,y 1 ,1) at the corresponding position in the non-template image through linear transformation. M can be solved by the corresponding three pairs of pixel coordinates in the template image and the non-template image;

所述三对像素坐标,选取为AprilTag标签的左上、右上、右下三个角点像素坐标,代入上式即可求解仿射变换矩阵。The three pairs of pixel coordinates are selected as the pixel coordinates of the upper left, upper right and lower right corner points of the AprilTag label, and the affine transformation matrix can be solved by substituting them into the above formula.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤10)的非模板图部件定位,具体为:The non-template image component positioning in step 10) is specifically as follows:

A1:从Redis数据库中读取所有key为“*P”的数据,遍历比较value中标签ID值和当前非模板图片中检测的标签ID值,当ID相等时,得到当前非模板图对应的高压开关柜名称;A1: Read all data with the key "*P" from the Redis database, traverse and compare the tag ID value in the value with the tag ID value detected in the current non-template image, and when the IDs are equal, obtain the name of the high-voltage switchgear corresponding to the current non-template image;

A2:从Redis数据库中读取所有key为“*P-PSx”的数据,“*P”为当前非模板图片对应高压开关柜名称,将模板图片区域矩形框的左上和右下角点的像素坐标代入下式,得到非模板图片的区域矩形框坐标,按照变换后的区域矩形框切割非模板图片;A2: Read all data with the key "*P-PSx" from the Redis database, where "*P" is the name of the high-voltage switchgear corresponding to the current non-template image. Substitute the pixel coordinates of the upper left and lower right corners of the template image's area rectangle into the following formula to obtain the coordinates of the non-template image's area rectangle. Cut the non-template image according to the transformed area rectangle.

Figure BDA0003687934340000052
Figure BDA0003687934340000052

A3:从Redis数据库中读取所有key为“*P-PSx-LEDy”、“*P-PSx-SWITCHy”、“*P-PSx-METERy”的数据,将模板图片部件矩形框的左上和右下角点的像素坐标代入上式,得到非模板图片的部件矩形框坐标,按照变换后的部件矩形框定位非模板图区域图像的部件。A3: Read all data with the key "*P-PSx-LEDy", "*P-PSx-SWITCHy", and "*P-PSx-METERy" from the Redis database, substitute the pixel coordinates of the upper left and lower right corners of the template image component rectangle into the above formula to obtain the coordinates of the component rectangle of the non-template image, and locate the component of the non-template image area image according to the transformed component rectangle.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤14)的制作0~9的LED数字二值化模板图,具体为:The step 14) of making the LED digital binary template image of 0 to 9 is specifically as follows:

B1:收集包含0~9数字的多张数字表部件图像,进行图像灰度处理和二值化处理,二值化阈值设定为150,即像素点的灰度值大于150则设置为255,否则设置为0;B1: Collect multiple digital table component images containing digits 0 to 9, perform image grayscale processing and binarization, and set the binarization threshold to 150, that is, if the grayscale value of the pixel is greater than 150, it is set to 255, otherwise it is set to 0;

B2:通过OpenCV中cv2.findContours方法查找二值化图像的轮廓,得到轮廓点[(x2,y2),...,(xn,yn)],遍历比较x1~xn和y1~yn的最大值和最小值,得到包围数字LED最小矩形的左上角像素坐标(xmin,ymin)和右下角像素坐标(xmax,ymax);B2: Use the cv2.findContours method in OpenCV to find the contour of the binary image and obtain the contour points [(x 2 ,y 2 ),...,(x n ,y n )]. Then, traverse and compare the maximum and minimum values of x 1 to x n and y 1 to yn to obtain the upper left corner pixel coordinates (x min, y min ) and lower right corner pixel coordinates (x max ,y max ) of the minimum rectangle surrounding the digital LED.

B3:对数字表部件图像按照(x,y,w,h)四个参数切割,参数取值为:B3: Cut the digital table component image according to the four parameters (x, y, w, h), the parameter values are:

Figure BDA0003687934340000061
Figure BDA0003687934340000061

B4:重设数字二值化模板图图像尺寸,统一设置为(42,51),保存数字模板图像,并按照数字实际含义命名为“num.jpg”,“num”取值范围是自然数0~9。B4: Reset the image size of the digital binary template image to (42,51), save the digital template image, and name it "num.jpg" according to the actual meaning of the number. The value range of "num" is a natural number from 0 to 9.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤14)的计算像素匹配度,从数字0~9按照图片名称,依次读取单个数字二值化模板图,并计算当前单个数字二值化图像和模板图像的匹配度:The calculation of pixel matching degree in step 14) reads the single digital binary template image in sequence from numbers 0 to 9 according to the picture name, and calculates the matching degree between the current single digital binary image and the template image:

Figure BDA0003687934340000062
Figure BDA0003687934340000062

其中,b0为当前单个数字图像的二值化像素矩阵,b1为读入的模板图的二值化像素矩阵。Among them, b0 is the binary pixel matrix of the current single digital image, and b1 is the binary pixel matrix of the read-in template image.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤14)的匹配度排序,采用冒泡排序法,得到匹配度最大值,对应二值化模板图的数字即为当前单个数字的读数结果,插入到result列表中。The matching degree sorting in step 14) adopts the bubble sorting method to obtain the maximum matching degree. The number corresponding to the binary template image is the reading result of the current single number and is inserted into the result list.

作为上述技术方案的进一步描述:As a further description of the above technical solution:

所述步骤14)的计算数字表实际读数,按照位数大小含义,计算读数:The actual reading of the digital meter in step 14) is calculated according to the meaning of the digit size:

Figure BDA0003687934340000071
Figure BDA0003687934340000071

其中,l是result列表的长度,即检测到的数字的个数,result(k)是result列表的第k个元素。Among them, l is the length of the result list, that is, the number of detected numbers, and result(k) is the kth element of the result list.

本发明具有如下有益效果:The present invention has the following beneficial effects:

与现有技术相比,该一种变电站高压开关柜二次仓面板图像识别系统,通过相比于以往的高压开关柜二次仓面板图像识别算法,本发明方法识别准确率高,对其他类型的高压开关柜具有通用性,仅一个ResNet50分类模型消耗GPU算力,因此还具有算力消耗低的优点;依据本发明方法,对于每面高压开关柜二次仓面板,仅需采集5~10张图像,即可实现算法设计与训练,解决了对大量数据依赖问题;依据本发明方法,每面高压开关柜二次仓面板图像仅需对所有部件标注1次,就可以实现每张图像中部件的精准定位,解决了对繁琐标注数据的问题;根据变电站高压开关柜二次仓面板的特点,本发明选择了传统图像处理和深度学习模型相结合的方式,将整张图像同时识别指示灯、压板开关、数字表问题,简化为指示灯和压板开关分类问题、数字表读数问题,弥补了单一图像处理方式的缺陷。Compared with the prior art, the image recognition system for the secondary compartment panel of the high-voltage switch cabinet of the substation has a high recognition accuracy compared with the previous image recognition algorithm for the secondary compartment panel of the high-voltage switch cabinet, and is universal for other types of high-voltage switch cabinets. Only one ResNet50 classification model consumes GPU computing power, so it also has the advantage of low computing power consumption; according to the method of the present invention, for each secondary compartment panel of the high-voltage switch cabinet, only 5 to 10 images need to be collected to realize algorithm design and training, which solves the problem of dependence on a large amount of data; according to the method of the present invention, for each secondary compartment panel image of the high-voltage switch cabinet, all components only need to be marked once, so that the components in each image can be accurately positioned, which solves the problem of cumbersome marking data; according to the characteristics of the secondary compartment panel of the high-voltage switch cabinet of the substation, the present invention selects a combination of traditional image processing and deep learning models to simultaneously identify the indicator light, pressure plate switch, and digital meter problem in the entire image, which is simplified to the classification problem of the indicator light and the pressure plate switch, and the reading problem of the digital meter, making up for the defects of the single image processing method.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明12面高压开关柜选择的模板图片;FIG1 is a picture of a template selected for a 12-sided high-voltage switch cabinet of the present invention;

图2为本发明对4P和6P高压开关柜图像区域划分的图;FIG2 is a diagram of the image area division of 4P and 6P high-voltage switchgears according to the present invention;

图3为本发明对14P高压开关柜图像区域切割结果的图;FIG3 is a diagram showing the results of cutting the image area of a 14P high-voltage switch cabinet according to the present invention;

图4为本发明在14P高压开关柜区域图像中部件标注结果的图;FIG4 is a diagram showing the component labeling results in the 14P high-voltage switchgear area image according to the present invention;

图5为本发明获取的指示灯数据集与分类的图;FIG5 is a diagram of indicator light data set and classification obtained by the present invention;

图6为本发明获取的压板开关数据集与分类的图;FIG6 is a diagram of a pressure plate switch data set and classification obtained by the present invention;

图7为本发明在ResNet50分类算法模型训练的损失曲线;FIG7 is a loss curve of the present invention in the training of the ResNet50 classification algorithm model;

图8为本发明数字表自动读数算法处理流程框图。FIG8 is a flowchart of the automatic reading algorithm of the digital meter of the present invention.

具体实施方式DETAILED DESCRIPTION

参照图1-8,本发明提供的一种变电站高压开关柜二次仓面板图像识别系统,包括以下步骤:1-8, a substation high-voltage switch cabinet secondary compartment panel image recognition system provided by the present invention includes the following steps:

1)在每一面高压开关柜的二次仓面板水平中轴线上,张贴一张ID唯一的AprilTag标签,并记录高压开关柜名称和标签ID一一对应列表;1) On the horizontal center axis of the secondary compartment panel of each high-voltage switch cabinet, post an AprilTag label with a unique ID, and record the one-to-one correspondence between the high-voltage switch cabinet name and the label ID;

2)采集高压开关柜样本图片,构成样本库,每张样本图片包含变电站高压开关柜二次仓面板和AprilTag标签;2) Collect sample images of high-voltage switchgear to form a sample library. Each sample image contains the secondary compartment panel and AprilTag label of the substation high-voltage switchgear;

3)从样本库中,为每面高压开关柜选取1张清晰图片作为模板图;3) Select one clear picture from the sample library as a template picture for each high-voltage switchgear cabinet;

4)调用AprilTag算法检测模板图片中所有标签ID和标签四个角点像素坐标,计算标签中心点水平像素坐标和样本图片水平中轴线像素坐标差值,差值最小的即为当前高压开关柜上对应的标签;4) Call the AprilTag algorithm to detect all tag IDs and the pixel coordinates of the four corner points of the tags in the template image, and calculate the difference between the horizontal pixel coordinates of the center point of the tag and the horizontal central axis pixel coordinates of the sample image. The one with the smallest difference is the corresponding tag on the current high-voltage switchgear;

5)按照步骤4)中获取的对应标签ID,查找步骤1)中的名称和ID列表,获取模板图片所采集的高压开关柜名称,并存储模板图名称、标签ID和四个角点像素坐标到数据库中;5) According to the corresponding tag ID obtained in step 4), search the name and ID list in step 1), obtain the name of the high-voltage switchgear collected by the template image, and store the template image name, tag ID and four corner point pixel coordinates in the database;

6)依据高压开关柜二次仓面板上部件类型和分布特征,将每一张模板图划分指示灯区域、压板开关区域、数字表区域,切割区域图像并命名;6) According to the component types and distribution characteristics on the secondary compartment panel of the high-voltage switch cabinet, each template image is divided into the indicator light area, the pressure plate switch area, and the digital meter area, and the regional images are cut and named;

7)在区域图像中标注部件,命名部件名称,按照部件标注框逐一切割部件图像,并设定每个部件的标准状态;7) Mark the components in the area image, name the components, cut the component images one by one according to the component marking frame, and set the standard state of each component;

8)从样本库中依次读取非模板图片,重复步骤4)获取每张非模板图片中标签ID和四个角点像素坐标;8) Read non-template images from the sample library in sequence, and repeat step 4) to obtain the label ID and four corner pixel coordinates in each non-template image;

9)根据步骤8)中标签角点像素坐标,结合步骤4)中模板图中标签角点像素坐标,计算模板图片和非模板图片之间的仿射变换矩阵;9) Calculate the affine transformation matrix between the template image and the non-template image based on the pixel coordinates of the label corner points in step 8) and the pixel coordinates of the label corner points in the template image in step 4);

10)根据仿射变换矩阵,将模板图片中的区域标注框映射到非模板图片中,切割区域图像并命名保存,再将模板图中部件标注框全部映射到非模板图片中,实现非模板图部件定位;10) According to the affine transformation matrix, the area annotation box in the template image is mapped to the non-template image, the area image is cut and named and saved, and then all the component annotation boxes in the template image are mapped to the non-template image to realize the positioning of the non-template image components;

11)按照定位的部件标注框,切割样本库中每张非模板图片的部件图像,并命名每张部件图像,获得部件图像数据集;11) According to the located component annotation box, cut the component image of each non-template image in the sample library, name each component image, and obtain the component image dataset;

12)将数据集按照部件图像名称分为指示灯数据集、压板开关数据集、数字表数据集;具体为:指示灯数据可分为两种类别,即指示灯亮、灭,分别记作LED_on、LED_off,数字表数据集即LED数字表显示不同读数,;12) Divide the data set into indicator light data set, pressure plate switch data set, and digital table data set according to the component image name; specifically, the indicator light data can be divided into two categories, namely, indicator light on and off, respectively recorded as LED_on and LED_off, and the digital table data set is the LED digital table showing different readings;

13)基于步骤12)获取的指示灯和压板开关数据集,进行数据增广,将数据集划分为训练集、验证集和测试集,在算法服务器上导入数据集和ResNet50分类模型,训练基于ResNet50的指示灯和压板开关分类算法模型,获取指示灯“亮”或“灭”状态和压板开关“开”或“合”状态;具体为:指示灯和压板开关数据集为指示灯数据集和压板开关数据集混合,分为LED_on、LED_off、SWITCH_on、SWITCH_off共四类部件状态;13) Based on the indicator light and pressure plate switch data set obtained in step 12), data augmentation is performed, and the data set is divided into a training set, a validation set, and a test set. The data set and the ResNet50 classification model are imported on the algorithm server, and the indicator light and pressure plate switch classification algorithm model based on ResNet50 is trained to obtain the "on" or "off" state of the indicator light and the "open" or "closed" state of the pressure plate switch; specifically: the indicator light and pressure plate switch data set is a mixture of the indicator light data set and the pressure plate switch data set, and is divided into four types of component states: LED_on, LED_off, SWITCH_on, and SWITCH_off;

14)基于步骤12)获取的数字表数据集,设计基于模板匹配的自动读数算法,制作0~9的LED数字二值化模板图,并逐一计算新图片中LED数字表单个数字二值化图像和数字二值化模板图的像素匹配度,对匹配度排序得到匹配度最高的单个数字读数,对数字位置按照水平坐标排序,按照位数计算数字表实际读数;14) Based on the digital table data set obtained in step 12), an automatic reading algorithm based on template matching is designed to produce a binary template image of LED digital numbers from 0 to 9, and the pixel matching degree between the single binary image of the LED digital table and the digital binary template image in the new image is calculated one by one, and the matching degree is sorted to obtain the single digital reading with the highest matching degree, and the digital positions are sorted according to the horizontal coordinates, and the actual reading of the digital table is calculated according to the number of digits;

15)将新采集的高压开关柜图片,输入到高压开关柜图像识别系统中,重复步骤4)和步骤5),获取图片中标签四个角点像素坐标,重复步骤6)获取区域图像;15) Input the newly acquired high-voltage switchgear image into the high-voltage switchgear image recognition system, repeat steps 4) and 5) to obtain the pixel coordinates of the four corner points of the label in the image, and repeat step 6) to obtain the regional image;

16)计算新采集图片和模板图片的仿射变换矩阵,将模板图中所有部件标注框映射到新采集图片上,切割部件图像并命名,按照部件图像名称中的部件类型,调用基于ResNet50的分类算法模型或基于模板匹配的自动读数算法,获取识别结果;16) Calculate the affine transformation matrix of the newly collected image and the template image, map all component annotation boxes in the template image to the newly collected image, cut the component image and name it, and call the classification algorithm model based on ResNet50 or the automatic reading algorithm based on template matching according to the component type in the component image name to obtain the recognition result;

17)读取每一个指示灯、压板开关、数字表预设的标准状态,逻辑判断部件状态是否异常,并及时告警;具体为:17) Read the preset standard status of each indicator light, pressure plate switch, and digital meter, logically determine whether the component status is abnormal, and issue an alarm in time; specifically:

对于指示灯和压板开关,当标准状态和识别状态相同时,判断为正常,当标准状态和识别状态不同时,判断为异常;For indicator lights and pressure plate switches, when the standard state and the identification state are the same, it is judged as normal, and when the standard state and the identification state are different, it is judged as abnormal;

对于数字表,当识别读数在预设阈值范围内,判断为正常,当识别读数大于或小于预设阈值范围,判断为异常。For a digital meter, when the identification reading is within the preset threshold range, it is judged as normal, and when the identification reading is greater than or less than the preset threshold range, it is judged as abnormal.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤1)的标签ID,具体是除0以外的自然数。The tag ID of step 1) is specifically a natural number other than 0.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤6)的二次仓面板上部件类型,主要包含LED指示灯、压板开关、LED数字表,依次命名为LED、SWITCH、METER;步骤6)的划分区域,采用区域矩形框划分,每面高压开关柜二次仓面板图像最多划分3个区域,对区域依次命名为“*P-PSx”,其中“*P”表示高压开关柜简称,“PSx”表示区域简称,“x”取值为1、2、3,取值不重复;区域矩形框,要求按照部件类型,框到所有的待识别部件,并以区域名称为key值,区域矩形框的左上和右下角点像素坐标为value,存储到Redis数据库中;步骤6)的切割区域图像,按照区域矩形框切割图片,切割后按照区域名称对区域图像命名为“*P-PSx-time.jpg”,“time”为计算机当前毫秒数,并将区域图像保存在本地服务器中;本地服务器,操作系统为Linux Ubuntu。The component types on the secondary warehouse panel in step 6) mainly include LED indicator lights, pressure plate switches, and LED digital meters, which are named LED, SWITCH, and METER respectively; the divided areas in step 6) are divided by regional rectangular frames, and each high-voltage switch cabinet secondary warehouse panel image is divided into up to 3 areas, and the areas are named "*P-PSx" in sequence, where "*P" represents the abbreviation of the high-voltage switch cabinet, "PSx" represents the abbreviation of the area, and "x" takes the values of 1, 2, and 3, and the values are not repeated; the regional rectangular frame requires that all the components to be identified be framed according to the component type, and the regional name is used as the key value, and the pixel coordinates of the upper left and lower right corners of the regional rectangular frame are used as the value, and are stored in the Redis database; the cutting regional image in step 6) is to cut the picture according to the regional rectangular frame, and after cutting, the regional image is named "*P-PSx-time.jpg" according to the regional name, "time" is the current number of milliseconds on the computer, and the regional image is saved in the local server; the local server has an operating system of Linux Ubuntu.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤7)的标注部件,采用部件矩形框标注,对每一个部件命名,并设定高压开关柜正常运行状态下每一个部件的标准状态,存储部件标注框角点坐标和部件名称到数据库中,部件命名,结合区域名称、部件类型和编号,三种部件具体命名规则为“*P-PSx-LEDy”、“*P-PSx-SWITCHy”、“*P-PSx-METERy”,“y”表示非零自然数,在单个区域图像内,y取值唯一,存储部件标注框角点坐标和部件名称,以部件名称为key值,区域图像中部件标注框和部件标准状态为value,存储到Redis数据库;步骤7)的切割部件图像,对切割出每一个部件图像命名,根据部件名称将部件图像命名为“*P-PSx-LEDy-time.jpg”或“*P-PSx-SWITCHy-time.jpg”或“*P-PSx-METERy-time.jpg”,并将部件图像保存在本地服务器中。In step 7), the marked components are marked with component rectangular frames, each component is named, and the standard state of each component under normal operating conditions of the high-voltage switch cabinet is set. The coordinates of the corner points of the component marking frame and the component name are stored in the database. The component naming is combined with the area name, component type and number. The specific naming rules of the three components are "*P-PSx-LEDy", "*P-PSx-SWITCHy", and "*P-PSx-METERy". "y" represents a non-zero natural number. In a single area image, y has a unique value. The coordinates of the corner points of the component marking frame and the component name are stored, with the component name as the key value, and the component marking frame and the component standard state in the area image as the value, which are stored in the Redis database; in step 7), the cut component image is named for each cut component image, and the component image is named "*P-PSx-LEDy-time.jpg" or "*P-PSx-SWITCHy-time.jpg" or "*P-PSx-METERy-time.jpg" according to the component name, and the component image is saved in the local server.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤9)的计算仿射变换矩阵,具体为:The calculation of the affine transformation matrix in step 9) is specifically:

Figure BDA0003687934340000111
Figure BDA0003687934340000111

其中,M2×3为仿射变换矩阵,表示通过线性变化将模板图片中的像素齐次坐标(x0,y0,1)变换到非模板图片中对应位置像素齐次坐标(x1,y1,1),通过模板图和非模板图片中对应三对像素坐标可求解得到M;Among them, M 2×3 is the affine transformation matrix, which means that the pixel homogeneous coordinates (x 0 ,y 0 ,1) in the template image are transformed to the pixel homogeneous coordinates (x 1 ,y 1 ,1) at the corresponding position in the non-template image through linear transformation. M can be solved by the corresponding three pairs of pixel coordinates in the template image and the non-template image;

三对像素坐标,选取为AprilTag标签的左上、右上、右下三个角点像素坐标,代入上式即可求解仿射变换矩阵。Three pairs of pixel coordinates are selected as the upper left, upper right, and lower right corner pixel coordinates of the AprilTag label. Substituting them into the above formula can solve the affine transformation matrix.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤10)的非模板图部件定位,具体为:The non-template image component positioning in step 10) is specifically as follows:

A1:从Redis数据库中读取所有key为“*P”的数据,遍历比较value中标签ID值和当前非模板图片中检测的标签ID值,当ID相等时,得到当前非模板图对应的高压开关柜名称;A1: Read all data with the key "*P" from the Redis database, traverse and compare the tag ID value in the value with the tag ID value detected in the current non-template image, and when the IDs are equal, obtain the name of the high-voltage switchgear corresponding to the current non-template image;

A2:从Redis数据库中读取所有key为“*P-PSx”的数据,“*P”为当前非模板图片对应高压开关柜名称,将模板图片区域矩形框的左上和右下角点的像素坐标代入下式,得到非模板图片的区域矩形框坐标,按照变换后的区域矩形框切割非模板图片;A2: Read all data with the key "*P-PSx" from the Redis database, where "*P" is the name of the high-voltage switchgear corresponding to the current non-template image. Substitute the pixel coordinates of the upper left and lower right corners of the template image's area rectangle into the following formula to obtain the coordinates of the non-template image's area rectangle. Cut the non-template image according to the transformed area rectangle.

Figure BDA0003687934340000121
Figure BDA0003687934340000121

A3:从Redis数据库中读取所有key为“*P-PSx-LEDy”、“*P-PSx-SWITCHy”、“*P-PSx-METERy”的数据,将模板图片部件矩形框的左上和右下角点的像素坐标代入上式,得到非模板图片的部件矩形框坐标,按照变换后的部件矩形框定位非模板图区域图像的部件。A3: Read all data with the key "*P-PSx-LEDy", "*P-PSx-SWITCHy", and "*P-PSx-METERy" from the Redis database, substitute the pixel coordinates of the upper left and lower right corners of the template image component rectangle into the above formula to obtain the coordinates of the component rectangle of the non-template image, and locate the component of the non-template image area image according to the transformed component rectangle.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤14)的制作0~9的LED数字二值化模板图,具体为:Step 14) produces a binary template image of LED numbers 0 to 9, specifically:

B1:收集包含0~9数字的多张数字表部件图像,进行图像灰度处理和二值化处理,二值化阈值设定为150,即像素点的灰度值大于150则设置为255,否则设置为0;B1: Collect multiple digital table component images containing digits 0 to 9, perform image grayscale processing and binarization, and set the binarization threshold to 150, that is, if the grayscale value of the pixel is greater than 150, it is set to 255, otherwise it is set to 0;

B2:通过OpenCV中cv2.findContours方法查找二值化图像的轮廓,得到轮廓点[(x2,y2),...,(xn,yn)],遍历比较x1~xn和y1~yn的最大值和最小值,得到包围数字LED最小矩形的左上角像素坐标(xmin,ymin)和右下角像素坐标(xmax,ymax);B2: Use the cv2.findContours method in OpenCV to find the contour of the binary image and obtain the contour points [(x 2 ,y 2 ),...,(x n ,y n )]. Then, traverse and compare the maximum and minimum values of x 1 to x n and y 1 to yn to obtain the upper left corner pixel coordinates (x min, y min ) and lower right corner pixel coordinates (x max ,y max ) of the minimum rectangle surrounding the digital LED.

B3:对数字表部件图像按照(x,y,w,h)四个参数切割,参数取值为:B3: Cut the digital table component image according to the four parameters (x, y, w, h), the parameter values are:

Figure BDA0003687934340000131
Figure BDA0003687934340000131

B4:重设数字二值化模板图图像尺寸,统一设置为(42,51),保存数字模板图像,并按照数字实际含义命名为“num.jpg”,“num”取值范围是自然数0~9。B4: Reset the image size of the digital binary template image to (42,51), save the digital template image, and name it "num.jpg" according to the actual meaning of the number. The value range of "num" is a natural number from 0 to 9.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤14)的计算像素匹配度,从数字0~9按照图片名称,依次读取单个数字二值化模板图,并计算当前单个数字二值化图像和模板图像的匹配度:Step 14) calculates the pixel matching degree, reads the single digital binary template image in sequence from numbers 0 to 9 according to the picture name, and calculates the matching degree between the current single digital binary image and the template image:

Figure BDA0003687934340000132
Figure BDA0003687934340000132

其中,b0为当前单个数字图像的二值化像素矩阵,b1为读入的模板图的二值化像素矩阵。Among them, b0 is the binary pixel matrix of the current single digital image, and b1 is the binary pixel matrix of the read-in template image.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤14)的匹配度排序,采用冒泡排序法,得到匹配度最大值,对应二值化模板图的数字即为当前单个数字的读数结果,插入到result列表中。The matching degree sorting in step 14) adopts the bubble sorting method to obtain the maximum matching degree. The number corresponding to the binary template image is the reading result of the current single number and is inserted into the result list.

作为上述技术方案的进一步的实施方式:As a further implementation of the above technical solution:

步骤14)的计算数字表实际读数,按照位数大小含义,计算读数:Step 14) calculates the actual reading of the digital meter, and calculates the reading according to the meaning of the digit size:

Figure BDA0003687934340000141
Figure BDA0003687934340000141

其中,l是result列表的长度,即检测到的数字的个数,result(k)是result列表的第k个元素。Among them, l is the length of the result list, that is, the number of detected numbers, and result(k) is the kth element of the result list.

工作原理:Working principle:

通过无人机巡检采集高压开关柜二次仓面板图片,共计12面高压开关柜,单次巡检任务中,无人机对每面高压开关柜各采集1张图片,在不同时刻共执行10次巡检任务,即每面高压开关柜采集图像10张,构建样本库。单张图片像素为3000×4000,为每面高压开关柜各选取模板图片1张。The images of the secondary compartment panels of the high-voltage switch cabinets were collected through drone inspections. There were 12 high-voltage switch cabinets in total. In a single inspection task, the drone collected 1 image for each high-voltage switch cabinet. A total of 10 inspection tasks were performed at different times, that is, 10 images were collected for each high-voltage switch cabinet to build a sample library. The pixel of a single image is 3000×4000, and 1 template image was selected for each high-voltage switch cabinet.

不同批次采集的图片偏差较大,统计其余9次采图与模板图的像素偏差,平均偏差达到132.46个像素点,因此需要进行部件定位。There are large deviations between the images collected in different batches. The pixel deviations between the remaining 9 images and the template image are statistically analyzed, and the average deviation reaches 132.46 pixels. Therefore, component positioning is required.

根据高压开关柜部件类型和分布特征,对典型4P和6P高压开关柜二次仓面板图片划分区域。对所有12张模板图片划分区域,如表1所示。According to the type and distribution characteristics of high-voltage switchgear components, the secondary compartment panel images of typical 4P and 6P high-voltage switchgear are divided into regions. All 12 template images are divided into regions, as shown in Table 1.

表1 12面高压开关柜全称及区域划分Table 1 Full name and area division of 12 high-voltage switchgear

Figure BDA0003687934340000142
Figure BDA0003687934340000142

Figure BDA0003687934340000151
Figure BDA0003687934340000151

按照表1划分区域,对14P高压开关柜二次仓面板模板图片进行区域切割,切割出14P-PS1指示灯区域、14P-PS2指示灯区域、14P-PS3压板开关区域。Divide the areas according to Table 1, cut the area of the 14P high-voltage switch cabinet secondary compartment panel template image, and cut out the 14P-PS1 indicator light area, 14P-PS2 indicator light area, and 14P-PS3 pressure plate switch area.

对14P模板图片的区域图像14P-PS1、14P-PS2、14P-PS3进行部件标注,产生部件类型、部件标注框和部件标准状态。Component annotation is performed on the regional images 14P-PS1, 14P-PS2, and 14P-PS3 of the 14P template image to generate component types, component annotation frames, and component standard states.

读入14P高压开关柜的非模板图片,检测非模板图片中标签角点C1,读取Redis数据库存储的模板图片中标签角点C0:Read the non-template image of the 14P high-voltage switchgear, detect the label corner point C1 in the non-template image, and read the label corner point C0 in the template image stored in the Redis database:

Figure BDA0003687934340000152
Figure BDA0003687934340000152

Figure BDA0003687934340000153
Figure BDA0003687934340000153

选取前三对角点像素坐标,计算仿射变化矩阵M:Select the pixel coordinates of the first three diagonal points and calculate the affine transformation matrix M:

Figure BDA0003687934340000154
Figure BDA0003687934340000154

将14P模板图片中区域矩形框通过仿射变换映射到14P非模板图片中,并切割区域图像。The regional rectangular frame in the 14P template image is mapped to the 14P non-template image through affine transformation, and the regional image is cut.

将14P模板图片的区域图像中部件矩形框通过仿射变换映射到14P非模板图片的区域图像中,实现精准部件定位。The component rectangular frame in the regional image of the 14P template image is mapped to the regional image of the 14P non-template image through affine transformation to achieve accurate component positioning.

对于样本库中其余高压开关柜二次仓面板图片,重复上述过程,并按照仿射变换后部件标注框切割出部件图像,得到指示灯数据集、压板开关数据集、数字表数据集。For the remaining high-voltage switchgear secondary compartment panel images in the sample library, the above process is repeated, and the component images are cut out according to the component annotation box after affine transformation to obtain the indicator light data set, pressure plate switch data set, and digital table data set.

通过旋转、翻转等数据增广处理后,得到指示灯和压板开关数据集各类数据数量,如表2所示,并将数据集中每个类别部件图像分为训练集、验证集、测试集。After data augmentation processing such as rotation and flipping, the number of each type of data in the indicator light and pressure plate switch dataset is obtained, as shown in Table 2, and each category of component images in the dataset is divided into training set, verification set, and test set.

表2指示灯和压板开关数据集Table 2 Indicator light and pressure plate switch data set

Figure BDA0003687934340000161
Figure BDA0003687934340000161

基于ResNet50的LED指示灯和压板开关分类模型训练在算法服务器中进行,服务器操作系统为64位的Ubuntu16.04,基本配置如表3所示,分类模型是基于Pytorch1.6.0深度学习框架实现的,编程语言使用的是Python2.7。The training of the LED indicator and pressure plate switch classification model based on ResNet50 is carried out in the algorithm server. The server operating system is 64-bit Ubuntu16.04. The basic configuration is shown in Table 3. The classification model is implemented based on the Pytorch1.6.0 deep learning framework, and the programming language used is Python2.7.

表3算法服务器配置参数Table 3 Algorithm server configuration parameters

Figure BDA0003687934340000162
Figure BDA0003687934340000162

Figure BDA0003687934340000171
Figure BDA0003687934340000171

在ResNet50模型训练过程中,首先使用在ImageNet上预训练的ResNet50模型来初始化网络权重,并设置网络的主要训练参数如表4所示。During the ResNet50 model training process, the ResNet50 model pre-trained on ImageNet is first used to initialize the network weights, and the main training parameters of the network are set as shown in Table 4.

表4网络训练参数Table 4 Network training parameters

参数名称Parameter name 参数值Parameter Value 迭代次数/EpochEpoch 120120 初始学习率/LearningrateInitial learning rate/Learningrate 1.01.0 权重衰减因子/WeightdecayWeight decay factor/Weightdecay 0.00010.0001 动量衰减因子/MomentumMomentum 0.90.9 批量大小/BatchSizeBatch Size 3232

训练结束时网络训练产生的训练损失(training loss)和验证损失(validationloss)曲线。由图可知,网络的训练损失和验证损失同步收敛,ResNet50分类模型训练效果较好。The training loss and validation loss curves generated by network training at the end of training. As can be seen from the figure, the training loss and validation loss of the network converge synchronously, and the training effect of the ResNet50 classification model is good.

使用训练好的分类算法模型,对测试集中的部件图像进行分类测试,测试结果统计如表4-8所示,平均识别准确率达到99%,具有非常好的分类效果。Using the trained classification algorithm model, the component images in the test set were classified and tested. The test results are shown in Table 4-8. The average recognition accuracy reached 99%, which has a very good classification effect.

表4-8算法测试结果Table 4-8 Algorithm test results

Figure BDA0003687934340000172
Figure BDA0003687934340000172

对于数字表数据集,共计10张。按照本发明所提方法制作0~9单个数字模板图片。再通过灰度处理、二值化、查找轮廓等步骤,对数字表部件图像自动读数。在数字表数据集上测试,数字表自动读数准确率达到100%。For the digital table data set, there are 10 in total. According to the method proposed in the present invention, single digital template images of 0 to 9 are produced. Then, the digital table component images are automatically read through grayscale processing, binarization, contour search and other steps. Tested on the digital table data set, the automatic reading accuracy of the digital table reaches 100%.

由此可见,本发明能够实现变电站高压开关柜二次仓面板图像的自动识别,具有较高的准确率,并且具有采图数据量极少、通用性高等优点,能够应用于变电站智能巡检系统。It can be seen that the present invention can realize the automatic recognition of the secondary compartment panel image of the substation high-voltage switch cabinet with high accuracy, and has the advantages of extremely small amount of image collection data and high versatility, and can be applied to the intelligent inspection system of the substation.

最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。Finally, it should be noted that the above is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the aforementioned embodiments, it is still possible for those skilled in the art to modify the technical solutions described in the aforementioned embodiments, or to make equivalent substitutions for some of the technical features therein. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. A secondary bin panel image recognition system of a transformer substation high-voltage switch cabinet comprises the following steps:
1) Pasting an AprilTag label with unique ID on the horizontal central axis of a secondary bin panel of each high-voltage switch cabinet, and recording a list of the names and the label IDs of the high-voltage switch cabinets in one-to-one correspondence;
2) Collecting sample pictures of the high-voltage switch cabinet to form a sample library, wherein each sample picture comprises a secondary bin panel and an april tag of the high-voltage switch cabinet of the transformer substation;
3) 1 clear picture is selected from a sample library for each high-voltage switch cabinet to serve as a template picture;
4) Invoking an april tag algorithm to detect all tag IDs and tag four corner pixel coordinates in a template picture, calculating the difference value between the tag center point horizontal pixel coordinates and the sample picture horizontal central axis pixel coordinates, and obtaining the tag corresponding to the current high-voltage switch cabinet as the smallest difference value;
5) Searching the name and ID list in the step 1) according to the corresponding tag ID obtained in the step 4), obtaining the name of the high-voltage switch cabinet acquired by the template picture, and storing the name of the template picture, the tag ID and four corner pixel coordinates into a database;
6) Dividing each template map into an indicator light area, a pressing plate switch area and a digital table area according to the types and distribution characteristics of the components on the secondary bin panel of the high-voltage switch cabinet, cutting the images of the areas and naming the images;
7) Marking parts in the region image, naming the names of the parts, cutting the part images one by one according to the part marking frames, and setting the standard state of each part;
8) Sequentially reading non-template pictures from a sample library, and repeating the step 4) to obtain a tag ID and four corner pixel coordinates in each non-template picture;
9) Calculating an affine transformation matrix between the template picture and the non-template picture according to the coordinates of the marked corner pixels in the step 8) and the coordinates of the marked corner pixels in the template picture in the step 4);
10 Mapping the region labeling frame in the template picture into the non-template picture according to the affine transformation matrix, cutting the region image, naming and storing, and mapping the part labeling frame in the template picture into the non-template picture to realize the positioning of the non-template picture component;
11 Cutting part images of each non-template picture in the sample library according to the positioned part labeling frame, and naming each part image to obtain a part image data set;
12 Dividing the data set into an indicator lamp data set, a pressing plate switch data set and a digital table data set according to the names of the component images;
13 Based on the pilot lamp and the pressing plate switch data set obtained in the step 12), carrying out data augmentation, dividing the data set into a training set, a verification set and a test set, importing a data set and a ResNet50 classification model on an algorithm server, training the pilot lamp and the pressing plate switch classification algorithm model based on the ResNet50, and obtaining a pilot lamp on state or an off state and a pressing plate switch on state or an on state;
14 Based on the digital table data set obtained in the step 12), an automatic reading algorithm based on template matching is designed, LED digital binarization template diagrams of 0-9 are manufactured, pixel matching degrees of single digital binarization images of the LED digital table and the digital binarization template diagrams in the new pictures are calculated one by one, single digital reading with the highest matching degree is obtained by sequencing the matching degrees, the digital positions are sequenced according to horizontal coordinates, and actual reading of the digital table is calculated according to the number of bits;
15 Inputting the newly acquired high-voltage switch cabinet picture into a high-voltage switch cabinet image recognition system, repeating the step 4) and the step 5), acquiring the coordinates of four corner pixels of the label in the picture, and repeating the step 6) to acquire an area image;
16 Calculating affine transformation matrixes of the newly acquired picture and the template picture, mapping all component labeling frames in the template picture onto the newly acquired picture, cutting and naming component images, and calling a classification algorithm model based on ResNet50 or an automatic reading algorithm based on template matching according to the component types in the component image names to acquire an identification result;
17 Reading the preset standard state of each indicator light, the pressure plate switch and the digital table, logically judging whether the state of the component is abnormal or not, and giving an alarm in time.
2. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the tag ID of step 1) is specifically a natural number other than 0.
3. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the types of the components on the secondary bin panel in the step 6) mainly comprise LED indicator lamps, a pressing plate switch and an LED digital meter, which are named LED, SWITCH, METER in sequence; the dividing area of the step 6) is divided by adopting an area rectangular frame, the maximum of 3 areas are divided by the secondary bin panel image of each high-voltage switch cabinet, the areas are sequentially named as 'P-PSx', wherein 'P' represents the short name of the high-voltage switch cabinet, 'PSx' represents the short name of the area, and 'x' takes values of 1, 2 and 3 and the values are not repeated; the regional rectangular frame is required to be framed to all the parts to be identified according to the type of the parts, and pixel coordinates of upper left corner points and lower right corner points of the regional rectangular frame are value values according to the regional name as key values and are stored in a Redis database; the step 6) of cutting the area image, cutting the picture according to the area rectangular frame, naming the area image as 'P-PSx-time. Jpg' according to the area name after cutting, wherein 'time' is the current millisecond number of the computer, and storing the area image in a local server; and the local server is provided with an operating system of Linux Ubuntu.
4. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the labeling component in the step 7) adopts component rectangular frame labeling, names each component, sets the standard state of each component in the normal operation state of the high-voltage switch cabinet, stores the corner coordinates of the component labeling frame and the component names in a database, and stores the component names, the component standard states of the component labeling frame and the component in the region image as value in terms of key values by combining the region names, the component types and the numbers, wherein the specific naming rules of the three components are "xP-PSx-LEDy", "x-P-PSx-SWITCHy", "x-P-PSx-METERy", "y" represent nonzero natural numbers; the cutting part image of the step 7) is named after each part image is cut, the part image is named as 'P-PSx-LEDy-time. Jpg' or 'P-PSx-SWITCHy-time. Jpg' or 'P-PSx-METERy-time. Jpg' according to the part name, and the part image is stored in a local server.
5. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the calculated affine transformation matrix in the step 9) is specifically:
Figure FDA0003687934330000041
wherein M is 2×3 For affine transformation matrix, representing the pixel homogeneous coordinates (x 0 ,y 0 1) transforming to corresponding position pixel homogeneous coordinates (x) 1 ,y 1 1), obtaining M through solving corresponding three pairs of pixel coordinates in a template picture and a non-template picture;
and selecting three pairs of pixel coordinates, namely, the pixel coordinates of the upper left corner, the upper right corner and the lower right corner of the april tag label, and substituting the pixel coordinates into the pixel coordinates to solve the affine transformation matrix.
6. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the positioning of the non-template diagram component in the step 10) is specifically as follows:
a1: reading data with all keys of P from a Redis database, traversing and comparing a tag ID value in a value with a tag ID value detected in a current non-template picture, and obtaining a high-voltage switch cabinet name corresponding to the current non-template picture when the IDs are equal;
a2: reading data with all keys of P-PSx from a Redis database, wherein P is the name of a high-voltage switch cabinet corresponding to the current non-template picture, substituting pixel coordinates of upper left corner points and lower right corner points of a template picture region rectangular frame into the following formula to obtain region rectangular frame coordinates of the non-template picture, and cutting the non-template picture according to the transformed region rectangular frame;
Figure FDA0003687934330000042
a3: reading data of which all keys are "+" -P-PSx-LEDy "," -P-PSx-SWITCHy "," -P-PSx-METERy ", substituting the pixel coordinates of the upper left corner and the lower right corner of the rectangular frame of the template picture component into the upper formula, and obtaining the coordinates of the rectangular frame of the component of the non-template picture, and positioning the component of the non-template picture region image according to the transformed rectangular frame of the component.
7. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the step 14) of manufacturing the LED digital binarization template chart with 0 to 9 comprises the following specific steps:
b1: collecting a plurality of digital table part images containing 0-9 digits, performing image gray scale processing and binarization processing, wherein a binarization threshold value is set to 150, namely, if the gray scale value of a pixel point is greater than 150, 255 is set, otherwise, 0 is set;
b2: searching the outline of the binary image by using a cv2.findContours method in OpenCV to obtain outline points [ (x) 2 ,y 2 ),...,(x n ,y n )]Traversal comparison x 1 ~x n And y 1 ~y n And the minimum and maximum of (2) to obtain the upper left corner pixel coordinate (x) min, y min ) And lower right corner pixel coordinates (x max ,y max );
B3: cutting the digital table part image according to four parameters of (x, y, w and h), wherein the parameters are as follows:
Figure FDA0003687934330000051
b4: the image size of the digital binarization template image is reset and uniformly set as (42, 51), the digital template image is saved, and the numerical template image is named as 'num. Jpg' according to the actual meaning of the digital, and the value range of 'num' is natural numbers 0-9.
8. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the step 14) of calculating the pixel matching degree, sequentially reading single digital binarization template images from the numbers 0 to 9 according to the picture names, and calculating the matching degree of the current single digital binarization image and the template image:
Figure FDA0003687934330000061
wherein b 0 Binarized pixel matrix for current single digital image, b 1 The pixel matrix is binarized for the read template map.
9. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: and 14) sorting the matching degree, namely obtaining the maximum value of the matching degree by adopting an bubbling sorting method, wherein the number corresponding to the binarization template diagram is the reading result of the current single number, and inserting the reading result into a result list.
10. The substation high-voltage switch cabinet secondary bin panel image recognition system according to claim 1, wherein: the actual reading of the digital table is calculated in the step 14), and the reading is calculated according to the meaning of the number of bits:
Figure FDA0003687934330000062
where l is the length of the result list, i.e. the number of digits detected, and result (k) is the kth element of the result list.
CN202210651684.6A 2022-06-10 2022-06-10 A substation high-voltage switchgear secondary warehouse panel image recognition system Active CN114998271B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210651684.6A CN114998271B (en) 2022-06-10 2022-06-10 A substation high-voltage switchgear secondary warehouse panel image recognition system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210651684.6A CN114998271B (en) 2022-06-10 2022-06-10 A substation high-voltage switchgear secondary warehouse panel image recognition system

Publications (2)

Publication Number Publication Date
CN114998271A CN114998271A (en) 2022-09-02
CN114998271B true CN114998271B (en) 2023-06-13

Family

ID=83032379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210651684.6A Active CN114998271B (en) 2022-06-10 2022-06-10 A substation high-voltage switchgear secondary warehouse panel image recognition system

Country Status (1)

Country Link
CN (1) CN114998271B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118038199B (en) * 2023-12-25 2024-07-26 浙江大学 A one-shot semi-automatic annotation method and system for switchgear images based on homography transformation

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111386529A (en) * 2017-06-22 2020-07-07 宾通机器人公司 System and method for quickly identifying and processing image regions of interest
CN112115807A (en) * 2020-08-28 2020-12-22 国网福建省电力有限公司厦门供电公司 An error-proof simulation method and system for switching off operation of secondary equipment in substation based on image recognition

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107563396B (en) * 2017-08-10 2019-08-30 南京大学 A method for constructing an intelligent identification system for protection screens in power inspections
WO2022037281A1 (en) * 2020-08-19 2022-02-24 广西电网有限责任公司贺州供电局 Real-time substation field operation monitoring and alarm system based on machine vision
CN112183369A (en) * 2020-09-29 2021-01-05 国网上海市电力公司 A pointer meter reading identification method for unmanned inspection of substations
CN112150366A (en) * 2020-09-30 2020-12-29 广西亚像科技有限责任公司 Method for identifying states of upper pressure plate and indicator lamp of transformer substation control cabinet

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111386529A (en) * 2017-06-22 2020-07-07 宾通机器人公司 System and method for quickly identifying and processing image regions of interest
CN112115807A (en) * 2020-08-28 2020-12-22 国网福建省电力有限公司厦门供电公司 An error-proof simulation method and system for switching off operation of secondary equipment in substation based on image recognition

Also Published As

Publication number Publication date
CN114998271A (en) 2022-09-02

Similar Documents

Publication Publication Date Title
CN112070135B (en) Power equipment image detection method and device, power equipment and storage medium
CN111401419A (en) Improved RetinaNet-based employee dressing specification detection method
CN108921166A (en) Medical bill class text detection recognition method and system based on deep neural network
CN105574161B (en) A kind of brand logo key element recognition methods, device and system
CN112308054B (en) Automatic reading method of multifunctional digital meter based on target detection algorithm
CN112163553B (en) Material price accounting method, device, storage medium and computer equipment
CN112613339B (en) Automatic identification and examination method and device for electrical drawing
CN109446689A (en) DC converter station electrical secondary system drawing recognition methods and system
CN108764134A (en) A kind of automatic positioning of polymorphic type instrument and recognition methods suitable for crusing robot
CN111144749A (en) Intelligent labeling crowdsourcing platform for power grid images and working method
CN113065556A (en) Digital instrument positioning and digital identification method and device and computer equipment
CN114998271B (en) A substation high-voltage switchgear secondary warehouse panel image recognition system
CN109858428A (en) ANA flourescent sheet automatic identifying method based on machine learning and deep learning
CN114511718B (en) Intelligent management method and system for materials for building construction
CN110309843A (en) An Automatic Recognition Method for Multiple Types of Components in Power Equipment Images
CN113076730A (en) Test report intelligent management method, device, system and storage medium
CN109033322A (en) A kind of test method and device of multidimensional data
CN116129459A (en) Automatic generation method of plant station model
CN110941910A (en) A kind of intelligent auxiliary method and system for three-dimensional design review of power grid
CN110688900A (en) A method of withdrawing meter management based on image recognition
CN114240924A (en) Power grid equipment quality evaluation method based on digitization technology
CN114299523A (en) Auxiliary operation identification and correction analysis method and analysis system
CN109598293B (en) Unmanned aerial vehicle inspection aerial photo batch processing training method based on class balance sampling
CN117475221A (en) A method and system for identifying the operation status of AC and DC distribution panels.
CN114720097B (en) Optical detection system for detecting TGG wafer

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 572024 area C310, third floor, phase II standard plant, yazhouwan science and Technology City, Yazhou District, Sanya City, Hainan Province

Patentee after: Hainan Research Institute of Zhejiang University

Patentee after: ZHEJIANG University

Address before: 572024 area C310, third floor, phase II standard plant, yazhouwan science and Technology City, Yazhou District, Sanya City, Hainan Province

Patentee before: Research Institute of Hainan Zhejiang University

Patentee before: ZHEJIANG University