CN107423501B - 一种基于光伏电站逻辑图的阵列自动检测和编号方法 - Google Patents

一种基于光伏电站逻辑图的阵列自动检测和编号方法 Download PDF

Info

Publication number
CN107423501B
CN107423501B CN201710581686.1A CN201710581686A CN107423501B CN 107423501 B CN107423501 B CN 107423501B CN 201710581686 A CN201710581686 A CN 201710581686A CN 107423501 B CN107423501 B CN 107423501B
Authority
CN
China
Prior art keywords
array
photovoltaic
pixel
photovoltaic array
rectangular
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
CN201710581686.1A
Other languages
English (en)
Other versions
CN107423501A (zh
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.)
Nanjing University of Posts and Telecommunications
Original Assignee
Nanjing University of Posts and Telecommunications
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 Nanjing University of Posts and Telecommunications filed Critical Nanjing University of Posts and Telecommunications
Priority to CN201710581686.1A priority Critical patent/CN107423501B/zh
Publication of CN107423501A publication Critical patent/CN107423501A/zh
Application granted granted Critical
Publication of CN107423501B publication Critical patent/CN107423501B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • 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/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于光伏电站逻辑图的阵列自动检测和编号方法,该方法首先将光伏电站总逻辑图进行分区得到分区图像,去除分区图像的彩色信息得到灰度图,对灰度图进行二值化阈值分割得到分离背景和前景(即阵列轮廓)的二值图像,然后提取二值图像的矩形轮廓,并通过计算矩形轮廓的众数信息排除假轮廓和噪声干扰,提取准确的阵列矩形,最后对每一个阵列进行行列编号,并计算阵列四个顶点的大地坐标,再转化为GPS坐标,将行列编号和GPS坐标信息保存到数据文件中。本发明能够为无人机拍摄光伏电站红外图像进行巡航时的自动路线规划以及热斑检测与定位等应用提供阵列的位置和编号信息,为提高光伏电站故障检测的自动化水平提供技术支撑。

Description

一种基于光伏电站逻辑图的阵列自动检测和编号方法
技术领域
本发明涉及一种基于光伏电站逻辑图的阵列自动检测和编号方法,属于光伏电站中的阵 列故障自动分析技术领域。
背景技术
近年来,光伏发电产业在世界范围内进入了突飞猛进的发展阶段,光伏电站的规模也在 不断增长。大型光伏电站的各个阵列在地理位置上通常具有排列规则、数量多、分布面积广 等特点,导致故障检测的工作量和难度大大增加。目前,阵列故障(如热斑)检测的自动化 水平不高,还主要依靠人工通过电气测量的方法对所有阵列组件进行检测,人工成本高且效 率低。而通过人工手持红外热像仪进行检测的方法只能识别小范围局部阵列中的热斑,无法 自动生成故障位置和严重程度的分析报告,且存在一定的危险性。
发明内容
本发明目的是提供一种基于光伏电站逻辑图的阵列自动检测和编号方法,该方法能够对 光伏电站逻辑图中的所有阵列进行自动的检测和识别,并生成光伏阵列的行列编号和四个顶 点的大地坐标,由此为无人机拍摄电站红外图像时的自动规划巡航路线,以及检测到故障(如 热斑)时返回阵列的具体编号和位置信息提供支持,是实现电站故障自动分析和定位功能的 重要技术基础。
为了达到上述目的,本发明采用的技术方案包括如下步骤:
步骤1)定义利用制图工具设计的整个光伏电站总体平面图为总逻辑图I,逻辑图常以浅 色(如接近白色)为背景,以深色(如蓝色或接近黑色)的矩形轮廓表示具体阵列的大小和 位置;总逻辑图I以图像形式存储,具体存储格式不作限制,如为常见的JPEG格式。考虑到 电站具体设计时常根据地形影响分区域地规划阵列的安装位置,首先由人工根据总逻辑图I 中阵列之间不同的行间距进行分区,得到分区图像Ii,其中i=1,2,3,……,K,其中K为分区 图像的数目,分区图像Ii的存储格式亦不作限制,可以与总逻辑图I相同。
步骤2)对各个分区图像Ii依次进行处理,得到分离背景和前景(即光伏阵列)轮廓的二 值图像。具体步骤如下:首先,读入分区图像Ii,去除分区图像Ii的彩色信息,得到灰度图 像并进行归一化;之后,利用最大类间方差法计算灰度图的二值化分割阈值;最后,根据大 量实验统计,取二值化分割阈值为0.9对灰度图像进行二值化,具体操作是把灰度图像中像 素值小于0.9的区域设置为0,用黑色表示分割出的前景(即阵列)轮廓区域;而像素值大 于0.9的区域设置为1,用白色表示分割出的背景区域。
步骤3)对上述步骤2)得到的二值图像进行矩形轮廓提取和存储,主要分为两步。首先, 检测矩形轮廓,其步骤如下:系统性的扫描二值图像直至找到连通轮廓(即由二值图像中黑 色相连像素组成的轮廓形状)的一个像素点,以这个像素点为起始点,跟踪它的轮廓,并标 记边界上的像素;当跟踪轮廓完整闭合,即扫描过程回到此起始点时,完成一个轮廓的跟踪 过程;继续寻找新的连通区域起始点,重复上述跟踪过程,完成所有轮廓的跟踪检测。当没 有新的连通区域起始点时,二值图像的所有轮廓检测完成。
其次,计算包含轮廓的最小矩形,得到二值图像中所有的矩形轮廓。将矩形轮廓存储在 一维结构体数组中,数组中的元素代表矩形轮廓,且每个元素包含矩形轮廓左上角顶点的像 素坐标以及矩形轮廓的宽度和高度,单位为像素;设数组中的元素个数为N,即矩形轮廓的 个数。数组中存储的分区图像Ii的矩形轮廓的行列顺序为由右至左,由下至上,即数组中的 第一个矩形轮廓为分区图像Ii右下角的光伏阵列,数组中的最后一个矩形轮廓为分区图像Ii左上角的光伏阵列。
步骤4)通过计算矩形轮廓的众数信息,排除假轮廓和噪声干扰,提取准确的阵列矩形。 实际光伏电站中的每一个光伏阵列是规则的矩形,且实际宽度和高度为固定数值,但在电站 逻辑图设计和图像处理过程中会产生少许误差并存在噪声干扰,使得检测出来的少数矩形轮 廓的宽度和高度不完全一致。为增加阵列检测的鲁棒性,采取的解决方法如下,遍历上述步 骤3)中得到的数组,对存储在数组中的所有矩形轮廓的宽度和高度进行排序,找出数组中 矩形轮廓的宽度和高度重复出现次数最多的数值,定义该数值为矩形轮廓的宽度和高度的众 数,并记录众数出现的次数;之后利用光伏阵列的宽度和高度众数,计算光伏阵列的面积众 数(PV_Area),通过多次实验统计取面积众数在数值上增减50为阵列允许误差,即认为矩形 面积在[PV_Area-50,PV_Area+50]范围内的为阵列矩形轮廓,排除不在这个范围内的矩形轮 廓。
步骤5)为每一个光伏阵列确定行列编号。首先确定数组中矩形轮廓左上角像素横坐标的 最大值,以及像素横坐标之间的距离众数;再结合上述步骤4)中计算得到的光伏阵列的宽 度和高度众数来计算光伏阵列的行列编号。包括如下子步骤:
a)记光伏阵列的行与列编号为(PV_rows,PV_cols)。初始时置分区图像最右下角阵列 的行列编号为(1,1),之后检测到的阵列编号规则为:同一行中向左排列的阵列列号依次递 增,同一列中向上排列的阵列行号依次递增;
b)遍历数组中矩形轮廓左上角顶点的像素横坐标,将遍历结果按降序排列,则降序排列 中的第一个元素即为矩形轮廓像素横坐标的最大值,用Max_x表示;之后将相邻的两个矩形 轮廓左上角顶点的像素横坐标依次相减,得到像素横坐标之间的距离,取重复出现频数最高 的数值为像素横坐标之间的距离众数,用Distance_x表示;
c)确定每一个光伏阵列的行列编号,也就是对数组中的每个矩形轮廓进行编号,先确定 矩形轮廓的行编号,再确定其列编号。具体方法如下,用{(xi,yi),i=1,2,...,N}表示数组中 每个矩形轮廓左上角顶点的像素坐标,N为数组元素个数,遍历数组,判断当前矩形轮廓左 上角顶点的像素纵坐标yi与前一个矩形轮廓左上角顶点的像素纵坐标yi-1是否相等,若相等, 则该矩形轮廓与前一个矩形轮廓的行编号相同;否则,该矩形轮廓的行编号加1;再利用公 式(1)来确定每个矩形轮廓的列编号,直至数组遍历完成;
Figure BDA0001352458360000036
其中xi表示当前矩形轮廓左上角顶点的像素横坐标,PV_cols表示计算得到的当前矩形 轮廓的列编号,
Figure BDA0001352458360000034
表示向下取整符号,如
Figure BDA0001352458360000035
的值为5。
步骤6)根据人工设定的参考点及其对应的像素坐标来计算每一个光伏阵列四个顶点的 大地坐标。人工设定两个参考点分别为分区图像Ii左上角第一个光伏阵列的左上角大地坐标, 即数组中最后一个矩形轮廓的左上角大地坐标,记为A点;分区图像Ii右下角最后一个光伏 阵列的左上角大地坐标,即数组中第一个矩形轮廓的左上角大地坐标,记为B点。人工输入 A点与B点的大地坐标,以及光伏阵列的实际物理宽度(PV_W)和高度(PV_H),单位为米。 根据输入的A点、B点的大地坐标和光伏阵列的实际物理宽度(PV_W)和高度(PV_H),根据 比例关系来确定每一个光伏阵列四个顶点的大地坐标。包括如下子步骤:
a)输入A点与B点的大地坐标,分别用(PV_xA,PV_yA),(PV_xB,PV_yB)表示,以及光伏 阵列的实际物理宽度(PV_W)和高度(PV_H)。
b)设A点对应的像素坐标为(xN,yN),B点对应的像素坐标为(x1,y1),当 前光伏阵列左上角顶点的像素坐标为(xi,yi),像素坐标均为已知;可根据公式(2)
Figure RE-GDA0001443061100000051
Figure RE-GDA0001443061100000052
c)根据光伏阵列的实际物理宽度(PV_W)和高度(PV_H),结合b)中计 算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000053
利用公式(3) 计算每一个光伏阵列右上角大地坐标
Figure RE-GDA0001443061100000054
Figure RE-GDA0001443061100000055
d)同理,根据光伏阵列的实际宽度(PV_W)和高度(PV_H),结合b)中 计算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000056
利用公式(4) 计算每一个光伏阵列右下角大地坐标
Figure RE-GDA0001443061100000057
Figure RE-GDA0001443061100000058
e)最后,根据光伏阵列的实际宽度(PV_W)和高度(PV_H),结合b)中 计算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000059
利用公式(5) 计算每一个光伏阵列左下角大地坐标
Figure RE-GDA00014430611000000510
Figure RE-GDA00014430611000000511
步骤7)最后根据大地坐标与GPS坐标之间的转换关系,将生成的光伏阵列四个顶点的大 地坐标转换成GPS坐标,并将光伏阵列的GPS坐标和行列编号保存在数据文件中,为无人机 拍摄光伏电站红外图像进行巡航时的自动路线规划和热斑检测与定位等应用提供光伏阵列的 位置和编号信息。
有益效果:
1、本发明基于光伏电站逻辑图的阵列自动检测和编号方法,通过对逻辑图进行图像分析, 可实现自动的阵列编号和准确计算阵列顶点的大地坐标,可辅助无人机自动规划巡航路线, 返回热斑等故障阵列的编号和位置信息,实现阵列故障区域的定位和分析处理。
2、本发明逻辑图的分区处理,简化了实现方法的系统结构和复杂度,适用范围广,同时 也提高了定位的准确度。
3、本发明结合无人机自动巡航和红外图像分析,可以很好地完成光伏电站的自动故障检 测和分析,提高光伏电站的运维效率,并且能够为无人机拍摄光伏电站红外图像进行巡航时 的自动路线规划以及热斑检测与定位等应用提供阵列的位置和编号信息,为提高光伏电站故 障检测的自动化水平提供技术支撑。
附图说明
图1为本发明的光伏电站阵列检测和编号的整体流程图。
图2为本发明的每个分区图像阵列矩形检测和阵列大地坐标计算以及编号的详细流程图。
具体实施方式
下面结合说明书附图对本发明创造作进一步的详细说明。
如图1和图2所示,本发明的基于光伏电站逻辑图的阵列自动检测和编号方法,以光伏 电站的总逻辑图I为待处理图像,其具体实施方案如下:
1)首先由人工根据总逻辑图I中阵列之间不同的行间距进行分区,得到分区图像Ii,其 中i=1,2,3,……,K,其中K为分区图像的数目,分区图像Ii的存储格式亦不作限制,可以与 总逻辑图I相同,如为常见的JPEG格式。
2)对各个分区图像Ii依次进行处理,得到分离背景和前景(即光伏阵列)轮廓的二值图 像。具体步骤如下:首先,读入分区图像Ii,去除分区图像Ii的彩色信息,得到灰度图像并 进行归一化;之后,利用最大类间方差法计算灰度图的二值化分割阈值;最后,根据大量实 验统计,取二值化分割阈值为0.9对灰度图像进行二值化,具体操作是把灰度图像中像素值 小于0.9的区域设置为0,用黑色表示分割出的前景(即阵列)轮廓区域;而像素值大于0.9 的区域设置为1,用白色表示分割出的背景区域。
3)对步骤2)得到的二值图像进行矩形轮廓提取和存储,主要分为两步。首先,检测矩 形轮廓,其步骤如下:系统性的扫描二值图像直至找到连通轮廓的一个像素点,以这个像素 点为起始点,跟踪它的轮廓,并标记边界上的像素;当跟踪轮廓完整闭合,即扫描过程回到 此起始点时,完成一个轮廓的跟踪过程;继续寻找新的连通区域起始点,重复上述跟踪过程, 完成所有轮廓的跟踪检测。当没有新的连通区域起始点时,二值图像的所有轮廓检测完成。
其次,计算包含轮廓的最小矩形,得到二值图像中所有的矩形轮廓。将矩形轮廓存储在 一维结构体数组中,数组中的元素代表矩形轮廓,且每个元素包含矩形轮廓左上角顶点的像 素坐标以及矩形轮廓的宽度和高度,单位为像素;设数组中的元素个数为N,即矩形轮廓的 个数。数组中存储的分区图像Ii的矩形轮廓的行列顺序为由右至左,由下至上,即数组中的 第一个矩形轮廓为分区图像Ii右下角的光伏阵列,数组中的最后一个矩形轮廓为分区图像Ii左上角的光伏阵列。
4)通过计算矩形轮廓的众数信息,排除假轮廓和噪声干扰,提取准确的阵列矩形。为增 加阵列检测的鲁棒性,采取的解决方法如下,遍历步骤3)中得到的数组,对存储在数组中 的所有矩形轮廓的宽度和高度进行排序,找出数组中矩形轮廓的宽度和高度重复出现次数最 多的数值,定义该数值为矩形轮廓的宽度和高度的众数,并记录众数出现的次数;之后利用 光伏阵列的宽度和高度众数,计算光伏阵列的面积众数(PV_Area),通过多次实验统计取面 积众数在数值上增减50为阵列允许误差,即认为矩形面积在[PV_Area-50,PV_Area+50]范围 内的为阵列矩形轮廓,排除不在这个范围内的矩形轮廓。
5)为每一个光伏阵列确定行列编号。首先确定数组中矩形轮廓左上角像素横坐标的最 大值,以及像素横坐标之间的距离众数;再结合步骤4)中计算得到的光伏阵列的宽度和高 度众数来计算光伏阵列的行列编号。包括如下子步骤:
a)记光伏阵列的行与列编号为(PV_rows,PV_cols)。初始时置分区图像最右下角阵列 的行列编号为(1,1),之后检测到的阵列编号规则为:同一行中向左排列的阵列列号依次递 增,同一列中向上排列的阵列行号依次递增;
b)遍历数组中矩形轮廓左上角顶点的像素横坐标,将遍历结果按降序排列,则降序排列 中的第一个元素即为矩形轮廓像素横坐标的最大值,用Max_x表示;之后将相邻的两个矩形 轮廓左上角顶点的像素横坐标依次相减,得到像素横坐标之间的距离,取重复出现频数最高 的数值为像素横坐标之间的距离众数,用Distance_x表示;
c)确定每一个光伏阵列的行列编号,也就是对数组中的每个矩形轮廓进行编号,先确定 矩形轮廓的行编号,再确定其列编号。具体方法如下,用{(xi,yi),i=1,2,...,N}表示数组中 每个矩形轮廓左上角顶点的像素坐标,遍历数组,判断当前矩形轮廓左上角顶点的像素纵坐 标yi与前一个矩形轮廓左上角顶点的像素纵坐标yi-1是否相等,若相等,则该矩形轮廓与前 一个矩形轮廓的行编号相同;否则,该矩形轮廓的行编号加1;再利用公式(1)来确定每个 矩形轮廓的列编号,直至数组遍历完成;
Figure BDA0001352458360000062
其中xi表示当前矩形轮廓左上角顶点的像素横坐标,PV_cols表示计算得到的当前矩形 轮廓的列编号,
Figure BDA0001352458360000061
表示向下取整符号。
6)根据人工设定的参考点及其对应的像素坐标来计算每一个光伏阵列四个顶点的大地坐 标。人工设定两个参考点分别为分区图像Ii左上角第一个光伏阵列的左上角大地坐标,即数 组中最后一个矩形轮廓的左上角大地坐标,记为A点;分区图像Ii右下角最后一个光伏阵列 的左上角大地坐标,即数组中第一个矩形轮廓的左上角大地坐标,记为B点。人工输入A点 与B点的大地坐标,以及光伏阵列的实际物理宽度(PV_W)和高度(PV_H),单位为米。根据 输入的A点、B点的大地坐标和光伏阵列的实际物理宽度(PV_W)和高度(PV_H),根据比例 关系来确定每一个光伏阵列四个顶点的大地坐标。包括如下子步骤:
a)输入A点与B点的大地坐标,分别用(PV_xA,PV_yA),(PV_xB,PV_yB)表示,以及光伏 阵列的实际物理宽度(PV_W)和高度(PV_H)。
b)设A点对应的像素坐标为(xN,yN),B点对应的像素坐标为(x1,y1),当前 光伏阵列左上角顶点的像素坐标为(xi,yi),像素坐标均为已知;可根据公式(2) 计算每一个光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000091
Figure RE-GDA0001443061100000092
c)根据光伏阵列的实际物理宽度(PV_W)和高度(PV_H),结合b)中计 算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000093
利用公式(3) 计算每一个光伏阵列右上角大地坐标
Figure RE-GDA0001443061100000094
Figure RE-GDA0001443061100000095
d)同理,根据光伏阵列的实际宽度(PV_W)和高度(PV_H),结合b)中 计算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000096
利用公式(4) 计算每一个光伏阵列右下角大地坐标
Figure RE-GDA0001443061100000097
Figure RE-GDA0001443061100000098
e)最后,根据光伏阵列的实际宽度(PV_W)和高度(PV_H),结合b)中 计算得到的光伏阵列左上角大地坐标
Figure RE-GDA0001443061100000101
利用公式(5) 计算每一个光伏阵列左下角大地坐标
Figure RE-GDA0001443061100000102
Figure RE-GDA0001443061100000103
7)最后根据大地坐标与GPS坐标之间的转换关系,将生成的光伏阵列四个顶点的大地坐 标转换成GPS坐标,并将光伏阵列的GPS坐标和行列编号保存在数据文件中。

Claims (4)

1.一种基于光伏电站逻辑图的阵列自动检测和编号方法,其特征在于,所述方法包括如下步骤:
步骤1)定义利用制图工具设计的整个光伏电站总体平面图为总逻辑图I,首先由人工根据总逻辑图I中阵列之间不同的行间距进行分区,得到分区图像Ii,其中i=1,2,3,……,K,其中K为分区图像的数目,分区图像Ii的存储格式亦不作限制,能与总逻辑图I相同;
步骤2)对各个分区图像Ii依次进行处理,得到分离背景和前景轮廓的二值图像;
步骤3)对上述步骤2)得到的二值图像进行矩形轮廓提取和存储,分为两步,首先,检测矩形轮廓,其步骤包括:系统性的扫描二值图像直至找到连通轮廓的一个像素点,以这个像素点为起始点,跟踪它的轮廓,并标记边界上的像素,当跟踪轮廓完整闭合,即扫描过程回到此起始点时,完成一个轮廓的跟踪过程;继续寻找新的连通区域起始点,重复上述跟踪过程,完成所有轮廓的跟踪检测,当没有新的连通区域起始点时,二值图像的所有轮廓检测完成;
其次,计算包含轮廓的最小矩形,得到二值图像中所有的矩形轮廓,将矩形轮廓存储在一维结构体数组中,数组中的元素代表矩形轮廓,且每个元素包含矩形轮廓左上角顶点的像素坐标以及矩形轮廓的宽度和高度,单位为像素,数组中存储的分区图像Ii的矩形轮廓的行列顺序为由右至左,由下至上,即数组中的第一个矩形轮廓为分区图像Ii右下角的光伏阵列,数组中的最后一个矩形轮廓为分区图像Ii左上角的光伏阵列;
步骤4)通过计算矩形轮廓的众数信息,排除假轮廓和噪声干扰,提取准确的阵列矩形,为增加阵列检测的鲁棒性,采取的解决方法如下,遍历上述步骤3)中得到的数组,对存储在数组中的所有矩形轮廓的宽度和高度进行排序,找出数组中矩形轮廓的宽度和高度重复出现次数最多的数值,定义该数值为矩形轮廓的宽度和高度的众数,并记录众数出现的次数,之后利用光伏阵列的宽度和高度众数,计算光伏阵列的面积众数PV_Area,通过多次实验统计取面积众数在数值上增减50为阵列允许误差,即认为矩形面积在[PV_Area-50,PV_Area+50]范围内的为阵列矩形轮廓,排除不在这个范围内的矩形轮廓;
步骤5)为每一个光伏阵列确定行列编号,首先确定数组中矩形轮廓左上角像素横坐标的最大值,以及像素横坐标之间的距离众数;再结合上述步骤4)中计算得到的光伏阵列的宽度和高度众数来计算光伏阵列的行列编号;
步骤6)根据人工设定的参考点及其对应的像素坐标来计算每一个光伏阵列四个顶点的大地坐标,人工设定两个参考点分别为分区图像Ii左上角第一个光伏阵列的左上角大地坐标,即数组中最后一个矩形轮廓的左上角大地坐标,记为A点,分区图像Ii右下角最后一个光伏阵列的左上角大地坐标,即数组中第一个矩形轮廓的左上角大地坐标,记为B点,人工输入A点与B点的大地坐标,以及光伏阵列的实际物理宽度PV_W和高度PV_H,单位为米,根据输入的A点、B点的大地坐标和光伏阵列的实际物理宽度PV_W和高度PV_H,根据比例关系来确定每一个光伏阵列四个顶点的大地坐标;
步骤7)最后根据大地坐标与GPS坐标之间的转换关系,将生成的光伏阵列四个顶点的大地坐标转换成GPS坐标,并将光伏阵列的GPS坐标和行列编号保存在数据文件中,为无人机拍摄光伏电站红外图像进行巡航时的自动路线规划和热斑检测与定位应用提供光伏阵列的位置和编号信息。
2.根据权利要求1所述的一种基于光伏电站逻辑图的阵列自动检测和编号方法,其特征在于,所述步骤2)分离背景和前景轮廓的二值图像包括:首先,读入分区图像Ii,去除分区图像Ii的彩色信息,得到灰度图像并进行归一化;之后,利用最大类间方差法计算灰度图的二值化分割阈值;最后,根据大量实验统计,取二值化分割阈值为0.9对灰度图像进行二值化,具体操作是把灰度图像中像素值小于0.9的区域设置为0,用黑色表示分割出的前景轮廓区域;而像素值大于0.9的区域设置为1,用白色表示分割出的背景区域。
3.根据权利要求1所述的一种基于光伏电站逻辑图的阵列自动检测和编号方法,其特征在于,所述步骤5)中确定光伏阵列的行列编号的方法包括:
a)记光伏阵列的行与列编号为(PV_rows,PV_cols),初始时置分区图像最右下角阵列的行列编号为(1,1),之后检测到的阵列编号规则为:同一行中向左排列的阵列列号依次递增,同一列中向上排列的阵列行号依次递增;
b)遍历数组中矩形轮廓左上角顶点的像素横坐标,将遍历结果按降序排列,则降序排列中的第一个元素即为矩形轮廓像素横坐标的最大值,用Max_x表示;之后将相邻的两个矩形轮廓左上角顶点的像素横坐标依次相减,得到像素横坐标之间的距离,取重复出现频数最高的数值为像素横坐标之间的距离众数,用Distance_x表示;
c)确定每一个光伏阵列的行列编号,也就是对数组中的每个矩形轮廓进行编号,先确定矩形轮廓的行编号,再确定其列编号,用{(xi,yi),i=1,2,...,N}表示数组中每个矩形轮廓左上角顶点的像素坐标,N为数组元素个数,遍历数组,判断当前矩形轮廓左上角顶点的像素纵坐标yi与前一个矩形轮廓左上角顶点的像素纵坐标yi-1是否相等,若相等,则该矩形轮廓与前一个矩形轮廓的行编号相同;否则,该矩形轮廓的行编号加1;再利用公式(1)来确定每个矩形轮廓的列编号,直至数组遍历完成;
Figure FDA0002434696700000031
其中xi表示当前矩形轮廓左上角顶点的像素横坐标,PV_cols表示计算得到的当前矩形轮廓的列编号,
Figure FDA00024346967000000410
表示向下取整符号。
4.根据权利要求1所述的一种基于光伏电站逻辑图的阵列自动检测和编号方法,其特征在于,所述步骤6)根据比例关系来确定每一个光伏阵列四个顶点的大地坐标的方法如下:
a)输入A点与B点的大地坐标,分别用(PV_xA,PV_yA),(PV_xB,PV_yB)表示,以及光伏阵列的实际物理宽度PV_W和高度PV_H,大地坐标与物理宽度高度的单位都是米;
b)设A点对应的像素坐标为(xN,yN),B点对应的像素坐标为(x1,y1),当前光伏阵列左上角顶点的像素坐标为(xi,yi),像素坐标均为已知;可根据公式(2)计算每一个光伏阵列左上角大地坐标
Figure FDA0002434696700000041
Figure FDA0002434696700000042
c)根据光伏阵列的实际物理宽度PV_W和高度PV_H,结合b)中计算得到的光伏阵列左上角大地坐标
Figure FDA0002434696700000043
利用公式(3)计算每一个光伏阵列右上角大地坐标
Figure FDA0002434696700000044
Figure FDA0002434696700000045
d)同理,根据光伏阵列的实际宽度PV_W和高度PV_H,结合b)中计算得到的光伏阵列左上角大地坐标
Figure FDA0002434696700000046
利用公式(4)计算每一个光伏阵列右下角大地坐标
Figure FDA0002434696700000047
Figure FDA0002434696700000048
e)最后,根据光伏阵列的实际宽度PV_W和高度PV_H,结合b)中计算得到的光伏阵列左上角大地坐标
Figure FDA0002434696700000049
利用公式(5)计算每一个光伏阵列左下角大地坐标
Figure FDA0002434696700000051
Figure FDA0002434696700000052
至此,所有阵列的编号和四个顶点的大地坐标都已确定。
CN201710581686.1A 2017-07-17 2017-07-17 一种基于光伏电站逻辑图的阵列自动检测和编号方法 Active CN107423501B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710581686.1A CN107423501B (zh) 2017-07-17 2017-07-17 一种基于光伏电站逻辑图的阵列自动检测和编号方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710581686.1A CN107423501B (zh) 2017-07-17 2017-07-17 一种基于光伏电站逻辑图的阵列自动检测和编号方法

Publications (2)

Publication Number Publication Date
CN107423501A CN107423501A (zh) 2017-12-01
CN107423501B true CN107423501B (zh) 2020-05-19

Family

ID=60429969

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710581686.1A Active CN107423501B (zh) 2017-07-17 2017-07-17 一种基于光伏电站逻辑图的阵列自动检测和编号方法

Country Status (1)

Country Link
CN (1) CN107423501B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549309A (zh) * 2018-05-17 2018-09-18 青海黄河上游水电开发有限责任公司光伏产业技术分公司 一种监测系统及监测方法
CN108923749A (zh) * 2018-06-11 2018-11-30 东北电力大学 基于红外视频的光伏组件热斑检测定位方法
CN109299311B (zh) * 2018-09-18 2021-01-12 甘肃启远智能科技有限责任公司 光伏组件经纬度数据配置方法及装置
CN109408916A (zh) * 2018-10-11 2019-03-01 南京国电南自维美德自动化有限公司 一种图形逻辑组态的识别方法
CN109615695B (zh) * 2018-11-13 2023-02-17 远景能源(南京)软件技术有限公司 房屋外部的空间照片到屋顶cad图纸的自动化转换方法
CN111353335A (zh) * 2018-12-21 2020-06-30 核动力运行研究所 一种将控制层逻辑图转换为仿真系统组态文件的方法
CN110084226B (zh) * 2019-03-04 2022-08-09 华为数字能源技术有限公司 应用于光伏阵列的电子布局图获取方法、装置以及系统
CN111104947A (zh) * 2019-11-21 2020-05-05 凌云光技术集团有限责任公司 基于特征一致性的干扰轮廓滤除方法和系统
CN111583176B (zh) * 2020-03-31 2023-09-15 北京全路通信信号研究设计院集团有限公司 一种基于图像的防雷模拟盘元件故障检测方法及系统
CN111740699A (zh) * 2020-05-29 2020-10-02 南京航空航天大学 一种光伏板故障检测识别方法、装置和无人机
CN112132898A (zh) * 2020-09-21 2020-12-25 重庆中电自能科技有限公司 光伏组件定位方法
CN113963370A (zh) * 2021-10-27 2022-01-21 合肥阳光智维科技有限公司 一种光伏电站中光伏设备的编号方法及编号装置
CN116188761B (zh) * 2023-02-02 2024-02-02 杭州图谱光电科技有限公司 一种可用于规则形状阵列图的自动识别定位方法及装置
CN116488573B (zh) * 2023-04-24 2023-11-03 湖北国网华中科技开发有限责任公司 用于光伏站的智能运维方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009086238A3 (en) * 2007-12-21 2010-12-09 E. I. Du Pont De Nemours And Company Photovoltaic array and installation method
CN102308320A (zh) * 2009-02-06 2012-01-04 香港科技大学 从图像生成三维模型
US8410950B2 (en) * 2009-08-17 2013-04-02 Paceco Corp. Photovoltaic panel monitoring apparatus
CN103049035A (zh) * 2012-12-20 2013-04-17 浙江工业大学 光伏电池阵列局部可重构电气系统的重构优化方法
WO2013093153A1 (es) * 2011-12-21 2013-06-27 Abengoa Solar New Technologies, S.A. Método para la inspeccion automatizada de captadores solares fotovoltaicos instalados en plantas
CN103559738A (zh) * 2013-10-17 2014-02-05 河北能源工程设计有限公司 一种山地光伏电站布置方法
CN104268368A (zh) * 2014-08-22 2015-01-07 湖北省电力勘测设计院 一种光伏电站优化施工方法
CN105894487A (zh) * 2014-08-25 2016-08-24 武汉科技大学 钢材图像编号的提取与分割方法
CN106786763A (zh) * 2017-01-13 2017-05-31 东北电力大学 一种风电场增建光伏电站的集电系统网络优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102047432B (zh) * 2007-12-21 2013-07-17 纳幕尔杜邦公司 光伏阵列、框架及其安装和使用方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009086238A3 (en) * 2007-12-21 2010-12-09 E. I. Du Pont De Nemours And Company Photovoltaic array and installation method
CN102308320A (zh) * 2009-02-06 2012-01-04 香港科技大学 从图像生成三维模型
US8410950B2 (en) * 2009-08-17 2013-04-02 Paceco Corp. Photovoltaic panel monitoring apparatus
WO2013093153A1 (es) * 2011-12-21 2013-06-27 Abengoa Solar New Technologies, S.A. Método para la inspeccion automatizada de captadores solares fotovoltaicos instalados en plantas
CN103049035A (zh) * 2012-12-20 2013-04-17 浙江工业大学 光伏电池阵列局部可重构电气系统的重构优化方法
CN103559738A (zh) * 2013-10-17 2014-02-05 河北能源工程设计有限公司 一种山地光伏电站布置方法
CN104268368A (zh) * 2014-08-22 2015-01-07 湖北省电力勘测设计院 一种光伏电站优化施工方法
CN105894487A (zh) * 2014-08-25 2016-08-24 武汉科技大学 钢材图像编号的提取与分割方法
CN106786763A (zh) * 2017-01-13 2017-05-31 东北电力大学 一种风电场增建光伏电站的集电系统网络优化方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Image features for pixel-wise detection of solar photovoltaic arrays in aerial imagery using a random forest classifier;Jordan M. Malof et al;《 2016 IEEE International Conference on Renewable Energy Research and Applications (ICRERA)》;20161123;全文 *
Optimising direct normal insolation of rectangular PV platforms;Ross Edgar et al;《Solar Energy》;20161015;第136卷;全文 *
光伏阵列热斑的红外图像处理的研究;郭宝柱;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20170415;第2017年卷(第4期);全文 *
太阳能光伏阵列故障诊断的研究;杨维翰等;《安徽工业大学学报》;20031031;第20卷(第4期);全文 *

Also Published As

Publication number Publication date
CN107423501A (zh) 2017-12-01

Similar Documents

Publication Publication Date Title
CN107423501B (zh) 一种基于光伏电站逻辑图的阵列自动检测和编号方法
CN107492094A (zh) 一种高压电线绝缘子的无人机视觉检测方法
CN111062991A (zh) 一种光伏场景的巡检目标定位方法、设备及可读存储介质
CN112270320B (zh) 一种基于卫星影像校正的输电线路杆塔坐标校准方法
CN106324581B (zh) 一种基于体元的机载lidar建筑物检测方法
CN113379703A (zh) 基于Yolo-v4网络结构光伏面板暗斑缺陷检测方法
CN112183301B (zh) 一种建筑物楼层智能识别方法及装置
CN110659634A (zh) 一种基于颜色定位和字符分割的集装箱箱号定位方法
CN114863258B (zh) 海天线场景中基于视角转换检测小目标的方法
CN110473174B (zh) 一种基于图像计算铅笔精确数目的方法
CN109190434B (zh) 一种基于亚像素级角点检测的条码识别算法
CN115909059A (zh) 一种自然资源样本库建立方法和装置
CN111598780A (zh) 一种适用于机载LiDAR点云的地形自适应插值滤波方法
CN116168246A (zh) 一种用于铁路工程的弃渣场识别方法、装置、设备及介质
CN116704333B (zh) 一种基于激光点云数据的单株林木检测方法
CN114419430A (zh) 一种基于SE-U-Net++模型的耕地地块提取方法及装置
CN116844147A (zh) 一种基于深度学习的指针式仪表识别和异常警报方法
CN113435266A (zh) 基于极值点特征增强的fcos智能目标检测方法
CN114880730A (zh) 确定目标设备的方法、装置及光伏系统
CN115908370A (zh) 一种基于图像逆透视变换实现水位检测的方法
Carneiro et al. Digital urban morphometrics: automatic extraction and assessment of morphological properties of buildings
CN112446895B (zh) 棋盘格角点自动提取方法、系统、设备及介质
CN109919852B (zh) 一种光学遥感图像的屋顶视角校正方法
CN106909720B (zh) 一种有限元节点坐标快速提取方法
CN116051771A (zh) 一种基于无人机倾斜摄影模型的光伏bim屋顶自动建模方法

Legal Events

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