OpenCV(四)—— 车牌号识别

本节是车牌识别的最后一部分 —— 车牌字符识别,从一个完整的车牌图片到识别出车牌上的字符大致需要如下几步:

  • 预处理:将车牌图片灰度化、二值化,并去除识别时的干扰因素,比如车牌铆钉
  • 字符分割:将整个车牌图片按照每个字符分割成 7 个单独的字符图片保存到集合中
  • 字符识别:使用经过训练的数字、英文字符、中文字符的特征集合对字符图片进行识别,得到最终的车牌号

下面详解以上步骤。

1、预处理

图像识别的预处理工作一般都是灰度化、二值化这些图像“降噪”处理,当然我们这里还有一个特殊的处理,就是要祛除车牌图像上的铆钉,它是对识别准确度影响较大的一个因素。

1.1 AnnPredictor 预处理

新建一个类 AnnPredictor 用于进行字符识别:

#ifndef ANNPREDICTOR_H
#define ANNPREDICTOR_H#define ANNPREDICTOR_DEBUG#include <opencv2/opencv.hpp>
#include <string>
#include <opencv2/ml.hpp>using namespace std;
using namespace cv;
using namespace ml;class AnnPredictor {
public:AnnPredictor(const char* ann_model, const char* ann_zh_model);~AnnPredictor();// 字符识别string predict(Mat plate);private:// 用于数字和英文字符识别Ptr<ANN_MLP> ann;// 用于中文字符识别Ptr<ANN_MLP> ann_zh;// ANN 的 HOG 特征HOGDescriptor* annHog = nullptr;// 与 SvmPredictor 中的函数相同void getHOGFeatures(HOGDescriptor* svmHog, Mat src, Mat& dst);// 去除铆钉bool clearRivet(Mat &plate);// 验证单个字符尺寸bool verifyCharSize(Mat src);// 获取城市汉字字符在排序后的矩形集合中的索引int getCityIndex(vector<Rect> rects);// 获取中文字符所在的矩形void getChineseRect(Rect cityRect, Rect& chineseRect);// 识别车牌字符保存到 str_plate 中void predict(vector<Mat> plateCharMats, string& str_plate);// 汉字字符集合static string ZHCHARS[];// 数字与英文字符集合static char CHARS[];
};#endif // !ANNPREDICTOR_H

这样在 LicensePlateRecognizer 中调用 predict() 即可获取到车牌字符串:

LicensePlateRecognizer::LicensePlateRecognizer(const char* svm_model, const char* ann_model, const char* ann_zh_model)
{...annPredictor = new AnnPredictor(ann_model, ann_zh_model);
}LicensePlateRecognizer::~LicensePlateRecognizer()
{...if (annPredictor){delete annPredictor;annPredictor = nullptr;}
}string LicensePlateRecognizer::recognize(Mat src)
{// 1.车牌定位,使用 Sobel 算法定位// 2.精选车牌定位得到的候选车牌图,找到最有可能是车牌的图...// 3.对车牌图进行字符识别string str_plate = annPredictor->predict(plate);plate.release();return str_plate;
}

predict() 内,先进行预处理,像灰度化和二值化这些操作前面已出现过多次就不再赘述了:

string AnnPredictor::predict(Mat plate)
{// 1.预处理 // 1.1 灰度化Mat gray;cvtColor(plate, gray, COLOR_BGR2GRAY);// 1.2 二值化(非黑即白,对比更强烈)Mat shold;threshold(gray, shold, 0, 255, THRESH_OTSU + THRESH_BINARY);// 1.3 去铆钉if (!clearRivet(shold)){return string("未识别到车牌");}...
}

主要说一下如何去掉车牌图片上的铆钉。

1.2 去掉车牌上的铆钉

二值化后的图片仍能看到车牌上的铆钉,这是影响识别准确程度的一个干扰因素:

去柳钉前

因此我们要去掉它,去掉后的效果:

去柳钉后

去除铆钉的思路是,对车牌图像进行逐行扫描,如果这一行是铆钉,那么颜色跳变的次数应该为 4 次,远远少于正常字符的颜色跳变次数:

2024-4-4.去柳钉前颜色跳变示意图1

上面红线是扫描到铆钉的行,先是黑色,扫描到铆钉变为白色,离开铆钉再变为黑色,第二颗铆钉重复上述过程,因此有 4 次黑白之间的颜色跳变。而第二条红线扫描到正常字符,跳变次数远远大于 4,我们就用这个思路去除铆钉:

/**
* 通过一行的颜色跳变次数判断是否扫描到了铆钉,
* 一行最小跳变次数为 12,最大为 12 + 8 * 6 = 60。
* 如果该行是铆钉行,则将该行所有像素都涂成黑色(像素值为 0)
*/
bool AnnPredictor::clearRivet(Mat &plate)
{// 1.逐行扫描统计颜色跳变次数保存到集合中int minChangeCount = 12;vector<int> changeCounts;int changeCount;for (int i = 0; i < plate.rows; i++){for (int j = 0; j < plate.cols - 1; j++){int pixel_front = plate.at<char>(i, j);int pixel_back = plate.at<char>(i, j + 1);if (pixel_front != pixel_back){changeCount++;}}changeCounts.push_back(changeCount);changeCount = 0;}// 2.计算字符高度,即满足像素跳变次数的行数int charHeight = 0;for (int i = 0; i < plate.rows; i++){if (changeCounts[i] >= 12 && changeCounts[i] <= 60){charHeight++;}}// 3.判断字符高度 & 面积占整个车牌的高度 & 面积的百分比,排除不符合条件的情况// 3.1 高度占比小于 0.4 则认为无法识别float heightPercent = float(charHeight) / plate.rows;if (heightPercent <= 0.4){return false;}// 3.2 面积占比小于 0.15 或大于 0.5 则认为无法识别float plate_area = plate.rows * plate.cols;// countNonZero 返回非 0 像素点(即白色)个数,或者自己遍历找像素点为 255 的个数也可float areaPercent = countNonZero(plate) * 1.0 / plate_area;// 小于 0.15 就是蓝背景白字车牌确实达不到识别标准,大于 0.5 是因为// 黄背景黑子二值化会把背景转化为白色,由于前面的处理逻辑只能处理// 蓝背景车牌,所以黄色车牌的情况也直接认为不可识别if (areaPercent <= 0.15 || areaPercent >= 0.5){return false;}// 4.将小于最小颜色跳变次数的行全部涂成黑色for (int i = 0; i < plate.rows; i++){if (changeCounts[i] < minChangeCount){for (int j = 0; j < plate.cols; j++){plate.at<char>(i, j) = 0;}}}return true;
}

2、字符分割

接下来开始分割字符,主要可以分为三部分:

  1. 找到各个字符的轮廓,并生成轮廓对应的图片
  2. 对汉字字符的轮廓图片进行特殊处理
  3. 将所有 7 个字符的轮廓图片保存到集合中为字符识别做准备

下面详解。

2.1 生成字符轮廓图片

先通过寻找轮廓的函数 findContours() 找到轮廓,生成轮廓矩形集合 vec_ann_rects:

string AnnPredictor::predict(Mat plate)
{// 1.预处理 ...// 2.字符分割// 2.1 找轮廓// vector<Point>是点的集合,可以连成线,线的集合就是轮廓了vector<vector<Point>> contours;findContours(shold, // 输入的图像contours, // 轮廓,接收结果RETR_EXTERNAL, // 轮廓检索模式:外轮廓CHAIN_APPROX_NONE // 轮廓近似算法模式:不进行轮廓近似,保留所有的轮廓点);vector<Rect> vec_ann_rects;// 在原图上克隆一个用来画矩形Mat src_clone = plate.clone();for each (vector<Point> points in contours) {Rect rect = boundingRect(points);Mat rectMat = shold(rect);// rectangle(src_clone, rect, Scalar(0, 0, 255));// 尺寸判断,符合规格的放入 vec_sobel_rects 集合中if (verifyCharSize(rectMat)) {vec_ann_rects.push_back(rect);}}...
}

遍历 vec_ann_rects 生成轮廓矩形并生成与矩形对应的图片,这就是分割的字符图片。当然,在将它们存入集合前你需要过滤一下,因为不是所有轮廓都刚好是一个完整的字符。比如“渝”字,由于汉字相比于英文字符和数字,结构复杂,因此无法识别为整个字,而是识别出“渝”字中点或者某个局部部分:

2024-4-8.字符找轮廓效果图2

因为我们导出的车牌宽度只有 136 像素,所以放大后并不清晰,但是能看出来,英文字母和数字的轮廓只有一个,而“渝”的轮廓有多个。面对这种情况,我们先用 verifyCharSize() 进行尺寸校验,将不符合规格的字符图片过滤掉:

bool AnnPredictor::verifyCharSize(Mat src)
{// 最理想情况 车牌字符的标准宽高比float aspect = 45.0f / 90.0f;// 当前获得矩形的真实宽高比float realAspect = (float)src.cols / (float)src.rows;// 最小的字符高float minHeight = 10.0f;// 最大的字符高float maxHeight = 35.0f;// 1、判断高符合范围  2、宽、高比符合范围// 最大宽、高比 最小宽高比float error = 0.7f;float maxAspect = aspect + aspect * error;//0.85float minAspect = 0.05f;int plate_area = src.cols * src.rows;float areaPercent = countNonZero(src) * 1.0 / plate_area;if (areaPercent <= 0.8 && realAspect >= minAspect && realAspect <= maxAspect&& src.rows >= minHeight &&src.rows <= maxHeight) {return true;}return false;
}

汉字被识别为多个部分,无法通过 verifyCharSize() 的校验,会被过滤掉。也就是说,这一步中,我们只获取了英文和数字图片,汉字图片暂时还未获取。

2.2 获取汉字轮廓

思路是,先定位到汉字后面表示城市的那一位英文字符,再向左侧推导获取汉字轮廓。

首先,根据图片的横坐标对 vec_ann_rects 集合内的字符图片进行从左至右的排序:

string AnnPredictor::predict(Mat plate)
{...// 2.2 对矩形轮廓从左至右排序sort(vec_ann_rects.begin(), vec_ann_rects.end(), [](const Rect& rect1, const Rect& rect2) {return rect1.x < rect2.x;});...
}

然后获取到表示城市字符的轮廓索引:

string AnnPredictor::predict(Mat plate)
{...// 2.3 获取城市字符轮廓的索引int cityIndex = getCityIndex(vec_ann_rects);...
}

getCityIndex() 获取城市字符索引的思路是:车牌上共有 7 个字符,那么城市字符位于第 2 位,该字符中间横坐标一定在车牌水平方向的 1/7 ~ 2/7 之间:

/**
* 寻找城市字符(7 位字符中的第 2 位)轮廓索引
*/
int AnnPredictor::getCityIndex(vector<Rect> rects)
{int cityIndex = 0;for (int i = 0; i < rects.size(); i++){Rect rect = rects[i];int midX = rect.x + rect.width / 2;// 如果字符水平方向中点坐标在整个车牌水平坐标的// 1/7 ~ 2/7 之间,就认为是目标索引。136 是我们// 训练车牌使用的素材的车牌宽度if (midX < 136 / 7 * 2 && midX > 136 / 7){cityIndex = i;break;}}return cityIndex;
}

获取城市字符索引后,可以根据其横坐标推断出汉字字符的轮廓:

string AnnPredictor::predict(Mat plate)
{...// 2.4 推导汉字字符的轮廓Rect chineseRect;getChineseRect(vec_ann_rects[cityIndex], chineseRect);...
}

getChineseRect() 会将推断出的矩形保存到 chineseRect 中:

/**
* 通过城市字符的矩形,确定汉字字符的矩形
*/
void AnnPredictor::getChineseRect(Rect cityRect, Rect& chineseRect)
{// 把宽度稍微扩大一点以包含完整的汉字字符// 还有一层理解,就是汉字与城市字符之间的空隙也要计算进去float width = cityRect.width * 1.15;// 城市轮廓矩形的横坐标int x = cityRect.x;// 用城市矩形的横坐标减去汉字宽度得到汉字矩形的横坐标int newX = x - width;chineseRect.x = newX > 0 ? newX : 0;chineseRect.y = cityRect.y;chineseRect.width = width;chineseRect.height = cityRect.height;
}

2.3 保存所有字符图片

最后将 7 个字符的图片保存到 plateCharMats 集合中等待字符识别:

string AnnPredictor::predict(Mat plate)
{...// 2.5 将字符图像保存到集合中// 先保存汉字字符图像vector<Mat> plateCharMats;plateCharMats.push_back(shold(chineseRect));// 再获取汉字之后的 6 个字符并保存int count = 6;if (vec_ann_rects.size() < 6){return string("未识别到车牌");}for (int i = cityIndex; i < vec_ann_rects.size() && count; i++, count--){plateCharMats.push_back(shold(vec_ann_rects[i]));}...
}

3、字符识别

3.1 识别过程

调用 predict() 传入字符图片集合 plateCharMats,识别的字符结果保存在 str_plate 中:

string AnnPredictor::predict(Mat plate)
{...// 3.字符识别string str_plate;predict(plateCharMats, str_plate);for (Mat m : plateCharMats) {m.release();}// 4.释放 Matgray.release();shold.release();src_clone.release();return str_plate;
}

predict() 内遍历 plateCharMats,提取 HOG 特征后对汉字和数字英文分开识别,注意 ZHCHARS 与 CHARS 内的字符顺序要和训练样本存放顺序相同:

string AnnPredictor::ZHCHARS[] = { "川", "鄂", "赣", "甘", "贵", "桂", "黑", "沪", "冀", "津", "京", "吉", "辽", "鲁", "蒙", "闽", "宁", "青", "琼", "陕", "苏", "晋", "皖", "湘", "新", "豫", "渝", "粤", "云", "藏", "浙" };
char AnnPredictor::CHARS[] = { '0', '1', '2', '3', '4', '5', '6', '7', '8', '9', 'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'J', 'K', 'L', 'M', 'N', 'P', 'Q', 'R', 'S', 'T', 'U', 'V', 'W', 'X', 'Y', 'Z' };void AnnPredictor::predict(vector<Mat> plateCharMats, string& result)
{for (int i = 0; i < plateCharMats.size(); i++){Mat mat_plate_char = plateCharMats[i];// 提取 HOG 特征Mat features;getHOGFeatures(annHog, mat_plate_char, features);Mat sample = features.reshape(1, 1);Mat response;Point maxLoc;Point minLoc;if (i){// 字母和数字ann->predict(sample, response);minMaxLoc(response, 0, 0, &minLoc, &maxLoc);int index = maxLoc.x;result += CHARS[index];}else{// 汉字ann_zh->predict(sample, response);minMaxLoc(response, 0, 0, &minLoc, &maxLoc);int index = maxLoc.x;result += ZHCHARS[index];}}
}

至此代码结束,先来看一下效果:

2024-4-8.最终识别结果1

对于【渝G 83666】的识别结果为【渝G 80666】,错了一位,这与识别的算法,还有训练样本的数量都有关系。总的来说,Demo 提供了一种车牌识别的思路,但是准确度还是有限的。

3.2 样本制作

最后来说说样本是如何制作的。

识别车牌字符,需要所有字符训练的特征集合,即首位的汉字共 31 个字符、第二位英文字符 24 个(刨除 I 和 O 两个容易被误识别为 1 和 0)、以及后续位数中需要用到的 10 个数字。我们将训练样本分为两类:英文字符和数字的训练样本放入 ann 文件夹,汉字字符放入 ann_zh 文件夹。两个文件夹内都需要对字符进行编号:

  • ann 中 0 号文件夹是数字 0 的训练素材,1 号文件夹是数字 1 的训练素材,以此类推,总共是 10 + 24 = 24 个文件夹
  • ann_zh 中 1 号文件夹是“川”字的训练素材,2 号文件夹是“鄂”字的训练素材,共计 31 个文件夹

这些字符的编号顺序需要记录在一个额外的文档中,内容如下:

"川", "鄂", "赣", "甘", "贵", "桂", "黑", "沪", "冀", "津", "京", "吉", "辽", "鲁", "蒙", "闽", "宁", "青", "琼", "陕", "苏", "晋", "皖", "湘", "新", "豫", "渝", "粤", "云", "藏", "浙"
'0', '1', '2', '3', '4', '5', '6', '7', '8', '9', 'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'J', 'K', 'L', 'M', 'N', 'P', 'Q', 'R', 'S', 'T', 'U', 'V', 'W', 'X', 'Y', 'Z'

训练后可以得到两个特征集合文件 ann.xml(数字和英文字符)和 ann_zh.xml(汉字字符),正是初始化 AnnPredictor 的 ann 和 ann_zh 所加载的文件。目录结构如下:

2024-4-8.车牌字符样本训练目录

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/pingmian/6316.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux——命名管道

管道特点 只能用于具有具体祖先的进程之间的通信&#xff0c;通常&#xff0c;一个管道由一个进程创建&#xff0c;然后该进程调用fork&#xff0c;创建子进程&#xff0c;关闭相应的读写端&#xff0c;然后父子进程就可以通信了管道提供流式服务一般而言&#xff0c;进程退出…

一个肉夹馍思考的零耦合设计

刷抖音听说知识付费是普通人的一个收入增长点&#xff0c;写了三十几篇文章一毛钱没赚&#xff0c;感觉有点沮丧。天上下着小雨雨&#xff0c;稀稀嗦嗦的&#xff0c;由于了很久还是买了一个&#x1f928;。 忽然觉得生活有点悲催&#xff0c;现在已经变得斤斤计较&#xff0c;…

Modelsim自动仿真平台的搭建

Modelsim自动仿真平台的搭建 如果要搭建自动仿真平台脚本那就需要更改下面3个文件。run_simulation.bat、complie.do和wave.do文件。注&#xff1a;前提是安装了modulsim并且配置好了环境变量&#xff0c;这里不过多介绍。 一、下面是run_simulation.bat文件的内容 : 注释的…

四足机器人摆线规划程序

一、标准摆线公式 { x r ∗ ( θ − sin ⁡ ( θ ) ) y r ∗ ( 1 − cos ⁡ ( θ ) ) \left\{\begin{array}{l} xr *(\theta-\sin (\theta)) \\ yr *(1-\cos (\theta)) \end{array}\right. {xr∗(θ−sin(θ))yr∗(1−cos(θ))​ 这里的r表示摆线的圆的半径&#xff0c; θ \…

3:容器之分类和各种测试

文章目录 array测试array容器bsearch之前用qsort 可以极大提高搜索效率 测试vector容器list容器forward_list容器测试slistdeque容器stack容器quque为什么stack 和 queue这样的容器 没有迭代器的相关内容 multiset代码里面有使用 stl自带的find 和 容器自带的find 一般容器自带…

计算机视觉科普到实践

第一部分&#xff1a;计算机视觉基础 引言&#xff1a; 计算机视觉作为人工智能领域的一个重要分支&#xff0c;近年来取得了显著的进展。本文将带领读者深入了解计算机视觉的基础知识&#xff0c;并通过实践案例展示其应用。让我们一同探索这个令人着迷的领域吧&#xff01;…

docker compose mysql主从复制及orchestrator高可用使用

1.orchestrator 功能演示&#xff1a; 1.1 多级级联&#xff1a; 1.2 主从切换&#xff1a; 切换成功后&#xff0c;原来的主库是红色的&#xff0c;需要在主库的配置页面点击“start replication ”&#xff0c;重新连接上新的主库。 1.3 主从故障&#xff0c;从库自动切换新…

高精度(加减乘除)

1.加法 我们第一位存低位&#xff08;倒着存方便&#xff09; 下面是AC代码&#xff1a; #include<bits/stdc.h> using namespace std; const int N1e610; vector<int> add(vector<int> &A,vector<int> &B) {vector<int> c;int t0;for…

【JVM】从硬件层面和应用层面的有序性和可见性,到Java的volatile和synchronized

Java的关键字volatile保证了有序性和可见性&#xff0c;这里我试着从底层开始讲一下有序性和可见性。 一&#xff0c;一致性 数据如果同时被两个cpu读取了&#xff0c;如何保证数据的一致性&#xff1f;或者换句话说&#xff0c;cpu1改了数据&#xff0c;cpu2的数据就成了无效…

基于AT89C52单片机的智能热水器控制系统

点击链接获取Keil源码与Project Backups仿真图&#xff1a; https://download.csdn.net/download/qq_64505944/89242443?spm1001.2014.3001.5503 C 源码仿真图毕业设计实物制作步骤05 题 目 基于单片机的智能热水器系统 学 院 专 业 班 级 学 号 学生姓名 指导教师 完成日期…

242 基于matlab的3D路径规划

基于matlab的3D路径规划&#xff0c;蚁群算法&#xff08;ACO&#xff09;和天牛须&#xff08;BAS&#xff09;以及两种结合的三种优化方式&#xff0c;对3D路径规划的最短路径进行寻优。程序已调通&#xff0c;可直接运行。 242 3D路径规划 蚁群算法和天牛须 - 小红书 (xiaoh…

ant-design中的穿梭框提示文字修改

ant-design中的穿梭框提示文字修改 1.ant-design中的穿梭框提示文字修改 <a-transferv-model:target-keys"targetKeys":data-source"transform.list":filter-option"filterOption":list-style"{width: 100%,height: 500px,}":rowK…

unity入门学习笔记

文章目录 unity学习笔记熟悉界面窗口页面快捷键视图特点移动、旋转、缩放快捷键聚焦和隐藏 一些基本概念模型模型的导入一些补充 资源文件资源包的导出资源包的导入 轴心物体的父子关系空物体Global与localpivot与center 组件脚本基础我的第一个脚本 获取脚本组件本地坐标播放模…

顶顶顶顶顶顶顶顶顶顶顶顶

欢迎关注博主 Mindtechnist 或加入【智能科技社区】一起学习和分享Linux、C、C、Python、Matlab&#xff0c;机器人运动控制、多机器人协作&#xff0c;智能优化算法&#xff0c;滤波估计、多传感器信息融合&#xff0c;机器学习&#xff0c;人工智能等相关领域的知识和技术。关…

【电子通识】为什么用双绞线?双绞线抗干扰的原理是什么?

使用双绞线最大的理由是抗干扰。不仅可以防止别人干扰,也可以防止自己干扰别人。这与EMC中的EMS和EMI相对应(参考【EMC专题】电磁兼容--基本概念)。 双绞线是由一对带有绝缘层的铜线(绝缘层使两根线中的金属导体不会因为互碰而导致短路)以螺旋的方式缠绕在一起所构成的。通…

ElementUI从unpkg.com完整下载到本地的方法 - 解决unpkg.com不稳定的问题 - 自建镜像站 - 不想打包只想cdn一下

方法 方法1&#xff09;随便弄个文件夹&#xff0c;根据官网npm方法下载包&#xff0c;提取即可 npm i element-ui -S cd /node_modules/element-ui/ ls src 安装npm方法&#xff1a;https://nodejs.org/en 方法2&#xff09;不推荐 - 在github中搜索对应的库zip包&#xff0…

【网络原理】HTTP 协议的基本格式和 fiddler 抓包工具的用法

系列文章目录 【网络通信基础】网络中的常见基本概念 【网络编程】网络编程中的基本概念及Java实现UDP、TCP客户端服务器程序&#xff08;万字博文&#xff09; 【网络原理】UDP协议的报文结构 及 校验和字段的错误检测机制&#xff08;CRC算法、MD5算法&#xff09; 【网络…

PG实例连接访问控制

实例访问控制可以控制来自于不同主机&#xff0c;不同用户是否允许访问指定的数据库&#xff0c;以及验证方式。 与oracle中的连接管理器的功能相同&#xff0c;之前有写过一篇oracleCMAN连接管理器的配置实操&#xff1a; 配置oracle连接管理器&#xff08;cman&#xff09;…

深度学习每周学习总结P7(咖啡豆识别)

&#x1f368; 本文为&#x1f517;365天深度学习训练营 中的学习记录博客&#x1f356; 原作者&#xff1a;K同学啊 | 接辅导、项目定制 –来自百度网盘超级会员V5的分享 数据链接 提取码&#xff1a;7zt2 –来自百度网盘超级会员V5的分享 目录 0. 总结1. 数据导入及处理部分…

文章解读与仿真程序复现思路——电力自动化设备EI\CSCD\北大核心《电-氢-混氢天然气耦合的城市综合能源系统低碳优化调度》

本专栏栏目提供文章与程序复现思路&#xff0c;具体已有的论文与论文源程序可翻阅本博主免费的专栏栏目《论文与完整程序》 论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html 电网论文源程序-CSDN博客电网论文源…