CN118115047A - 一种城市公共建筑海绵设施运行效果评价方法 - Google Patents

Abstract

本发明涉及图像数据处理技术领域，提出了一种城市公共建筑海绵设施运行效果评价方法，包括：获取下沉式绿地全景灰度图像；根据下沉式绿地全景灰度图像获取下沉式绿地全景二值图像；根据下沉式绿地全景二值图像获取候选倒影区域；根据候选倒影区域计算倒影纹理熵；根据倒影纹理熵计算倒影相似指数；根据倒影相似指数获取倒影相似指数序列；根据倒影指数相似序列获取疑似倒影指数；根据疑似倒影指数计算倒影标记系数；根据倒影标记系数计算积水区域面积；根据积水区域面积对下沉式绿地的运行效果进行评价。本发明根据倒影标记系数计算积水区域面积对下沉式绿地的蓄水能力进行评价，提高城市公共建筑海绵设施运行效果评价的准确性。

Description

一种城市公共建筑海绵设施运行效果评价方法

技术领域

本发明涉及图像数据处理技术领域，具体涉及一种城市公共建筑海绵设施运行效果评价方法。

背景技术

智能绿色建筑理念是一种注重资源效率、环境可持续和居住舒适性的建筑设计和运营理念，通过建筑设计与现代技术、智能系统相结合降低建筑的环境影响，进一步提高了生活质量，促进了资源的可持续利用。其中海绵城市是智能绿色建筑中重要的组成部分，城市公共建筑海绵设施是指利用海绵城市的概念和技术进行创新，在城市公共建筑中应用可持续性雨水管理系统的主要目的是改善水资源的利用效率，降低洪涝风险，提升城市环境质量。

常用的城市公共建筑海绵设施包括绿色屋顶、雨水花园、下沉式绿地、植草沟、透水铺装、雨水湿地、蓄水池、雨水罐等，主要作用是实现对雨水的循环利用。目前城市公共建筑海绵设施的种类纷繁复杂且效果差异较大，对城市公共建筑海绵设施缺乏一个行之有效的评价方法，由于城市公共建筑海绵设施所处环境复杂，传统的评价方法对城市公共建筑海绵设施运行效果的分析误差较大。

发明内容

本发明提供一种城市公共建筑海绵设施运行效果评价方法，以解决城市公共建筑海绵设施中下沉式绿地运行效果分析误差较大的问题，所采用的技术方案具体如下：

本发明一个实施例一种城市公共建筑海绵设施运行效果评价方法，该方法包括以下步骤：

获取下沉式绿地全景灰度图像；

根据图像分割算法获取下沉式绿地全景灰度图像的下沉式绿地二值图像，基于所述下沉式绿地二值图像获取下沉式绿地全景灰度图像的候选倒影区域；根据下沉式绿地全景灰度图像的候选倒影区域获取所述候选倒影区域的倒影判定邻域，基于所述倒影判定邻域获取所述倒影判定邻域的滑动搜索区域；根据下沉式绿地全景灰度图像计算下沉式绿地全景灰度图像中每个像素点的倒影纹理熵；根据候选倒影区域的倒影判定邻域中滑动搜索区域、候选倒影区域中每个像素点的倒影纹理熵计算候选倒影区域的倒影相似指数；

根据候选倒影区域与候选倒影区域的倒影判定邻域内滑动搜索区域之间的倒影相似指数获取候选倒影区域的倒影相似指数序列；根据候选倒影区域的倒影相似指数序列计算候选倒影区域的疑似倒影指数；根据每个候选倒影区域的疑似倒影指数计算每个候选倒影区域的倒影标记系数；根据每个候选倒影区域的倒影标记系数计算下沉式绿地全景灰度图像的积水区域总面积；

根据下沉式绿地全景灰度图像的积水区域面积获取下沉式绿地的积水面积-时间曲线，基于所述积水面积-时间曲线对下沉式绿地的蓄水能力进行评价。

优选的，所述根据图像分割算法获取下沉式绿地全景灰度图像的下沉式绿地二值图像，基于所述下沉式绿地二值图像获取下沉式绿地全景灰度图像的候选倒影区域的方法为：

采用大津阈值分割算法获取下沉式绿地全景灰度图像的二值图像，采用形态学处理对所述二值图像进行腐蚀、膨胀操作，将所述二值图像的腐蚀、膨胀操作处理的结果图像作为所述下沉式绿地全景灰度图像的下沉式绿地二值图像；

采用边缘检测算法和连通域分析算法获取所述下沉式绿地二值图像中的初始积水区域，将所述下沉式绿地二值图像中初始积水区域包含的任意一个其它连通域作为下沉式绿地全景灰度图像的一个候选倒影区域。

优选的，所述根据下沉式绿地全景灰度图像的候选倒影区域获取所述候选倒影区域的倒影判定邻域，基于所述倒影判定邻域获取所述倒影判定邻域的滑动搜索区域的方法为：

获取下沉式绿地全景灰度图像中任意一个候选倒影区域的中心，计算所述任意一个候选倒影区域的中心到所述任意一个候选倒影区域的边界像素点的欧氏距离，将所有所述欧氏距离的计算结果按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的邻域判定序列，将所述邻域判定序列中最大值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的长，将所述邻域序列中最小值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的宽，将所述任意一个候选倒影区域的中心作为所述倒影判定邻域的中心，根据所述长、宽、中心确定所述任意一个候选倒影区域的倒影判定领域；

在所述任意一个候选倒影区域的倒影判定邻域中构建与所述任意一个候选倒影区域相同形状和大小的滑动搜索窗口，将所述滑动搜索窗口作为所述任意一个候选倒影区域的倒影判定邻域的滑动搜索区域。

优选的，所述根据下沉式绿地全景灰度图像计算下沉式绿地全景灰度图像中每个像素点的倒影纹理熵的方法为：

以下沉式绿地全景灰度图像中任意一个像素点为中心构建预设大小的窗口，将所述预设大小的窗口作为所述任意一个像素点的形状特征分析区域，获取所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵，将所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵的熵分别作为所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值，计算所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值的均值作为所述任意一个像素点的倒影纹理熵。

优选的，所述根据候选倒影区域的倒影判定邻域的滑动搜索窗口、候选倒影区域中每个像素点的倒影纹理熵计算候选倒影区域的倒影相似指数的方法为：

；

式中，表示下沉式绿地全景灰度图像中第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域的倒影相似指数；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域中第/>个像素点的倒影纹理熵；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中第/>个像素点的倒影纹理熵；/>分别表示下沉式绿地全景灰度图像中第/>个候选倒影区域、第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中像素点的数量；/>表示调节参数。

优选的，所述根据候选倒影区域与候选倒影区域的倒影判定邻域内滑动搜索区域之间的倒影相似指数获取候选倒影区域的倒影相似指数序列的方法为：

将下沉式绿地全景灰度图像中任意一个候选倒影区域与所述任意一个候选倒影区域的倒影判定邻域内所有滑动搜索区域之间的倒影相似指数按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的倒影相似指数序列。

优选的，所述根据候选倒影区域的倒影相似指数序列计算候选倒影区域的疑似倒影指数的方法为：

采用归一化算法获取下沉式绿地全景灰度图像中任意一个候选倒影区域的倒影相似指数序列的归一化结果，将所述任意一个候选倒影区域的倒影相似指数序列的归一化结果作为所述任意一个候选倒影区域的倒影相似指数归一化序列，将所述任意一个候选倒影区域的倒影相似指数归一化序列中的最大值作为所述任意一个候选倒影区域的疑似倒影指数。

优选的，所述根据每个候选倒影区域的疑似倒影指数计算每个候选倒影区域的倒影标记系数的方法为：

将下沉式绿地全景灰度图像中每个候选倒影区域的疑似倒影指数与预设阈值进行比较，将所述疑似倒影指数大于预设阈值的候选倒影区域的倒影标记系数记为1，将所述疑似倒影指数小于预设阈值的候选倒影区域的倒影标记系数记为0。

优选的，所述根据每个候选倒影区域的倒影标记系数计算下沉式绿地全景灰度图像的积水区域总面积的方法为：

；

式中，表示下沉式绿地全景灰度图像的积水区域总面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影标记系数；/>表示下沉式绿地全景灰度图像的初始积水面积；/>表示下沉式绿地全景灰度图像中候选倒影区域的数量。

优选的，所述根据下沉式绿地全景灰度图像的积水区域面积获取下沉式绿地的积水面积-时间曲线，基于所述积水面积-时间曲线对下沉式绿地的蓄水能力进行评价的方法为：

在雨后对下沉式绿地公园每隔一段时间采集一幅下沉式绿地全景图像，将每幅下沉式绿地全景图像对应的积水区域总面积作为纵坐标，将每幅下沉式绿地全景图像对应的采集时间作为横坐标，将所述横坐标和纵坐标构成的曲线作为所述下沉式绿地的积水面积-时间曲线，计算积水面积-时间曲线下面积作为下沉式绿地的蓄水能力特征值；

获取与下沉式绿地相同环境条件的非下沉式绿地的全景图像，根据下沉式绿地、非下沉式绿地的全景图像对应的蓄水能力特征值得到下沉式绿地的蓄水能力评价结果。

本发明的有益效果是：通过拍摄下沉式绿地全景图像获取下沉式绿地全景灰度图像，根据下沉式绿地全景灰度图像获取下沉式绿地的积水区域、候选倒影区域，根据下沉式绿地全景灰度图像获取倒影纹理熵，根据候选倒影区域获取倒影判定邻域、滑动搜索区域，根据候选倒影区域、滑动搜索区域、倒影纹理熵计算倒影相似指数，根据倒影相似指数获取疑似倒影指数，根据疑似倒影指数计算倒影标记系数，基于倒影标记系数计算下沉式绿地的积水区域总面积，其有益效果在于避免水中倒影对下沉式绿地积水面积计算的影响，通过积水区域总面积与时间的关系曲线计算下沉式绿地的蓄水指数，根据蓄水指数对下沉式绿地的蓄水能力进行评价，可以进一步提高对城市公共建筑海绵设施中下沉式绿地运行效果评价的准确性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明一个实施例所提供的一种城市公共建筑海绵设施运行效果评价方法的流程示意图；

图2为本发明一个实施例所提供的倒影判定邻域、滑动搜索区域示意图；

图3为本发明一个实施例所提供的积水面积-时间曲线下面积示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1，其示出了本发明一个实施例提供的一种城市公共建筑海绵设施运行效果评价方法的流程图，该方法包括以下步骤：

步骤S001，获取下沉式绿地全景灰度图像。

采用无人机搭载CCD工业相机对下沉式绿地进行图像采集，拍摄过程中无人机和相机的具体型号实施者可以根据实际情况进行选择，将采集的下沉式绿地的图像作为APAP（As-Projective-As-Possible Image Stitching with Moving DLT）算法的输入，根据APAP算法获取下沉式绿地全景RGB图像，将下沉式绿地全景RGB图像转化为灰度图像，输入为所述灰度图像，采用自适应直方图均衡化算法对所述灰度图像进行预处理，将所述预处理的结果图像作为下沉式绿地全景灰度图像，APAP算法、自适应直方图均衡化算法的具体实现过程为公知技术，不再进行赘述。

至此，获取了下沉式绿地全景灰度图像。

步骤S002，根据下沉式绿地全景灰度图像获取候选倒影区域，根据候选倒影区域获取倒影判定邻域、滑动搜索区域，获取下沉式全景图像中每个像素点的倒影纹理熵，根据倒影判定邻域、滑动搜索区域、倒影纹理熵计算倒影相似指数。

输入为下沉式绿地全景灰度图像，采用大津阈值分割算法获取下沉式绿地全景灰度图像的二值图像，将所述二值图像作为下沉式绿地二值图像，所述下沉式绿地二值图像中只包含1、0两种灰度值；输入为下沉式绿地二值图像，采用形态学处理中的膨胀操作、腐蚀操作对下沉式绿地二值图像进行处理，使下沉式绿地二值图像更加完整、连续，形态学处理中的膨胀操作、腐蚀操作的具体实现过程为公知技术，不再进行赘述。

进一步的，由于下沉式绿地中水面反射强度相较于其它物体较大，因此下沉式绿地全景图像中较亮的区域为水面区域，表现在下沉式绿地二值图像中则为像素点灰度值为1的区域，但是积水区域中存在倒影、凸起地面等会影响积水区域大小的判断，因此可通过对积水区域中倒影的特征对连通域进行筛选判断，得到完整的积水区域面积。

具体的，将经过形态学处理后的下沉式绿地二值图像作为输入，采用Canny边缘检测算法和连通域分析算法获取下沉式绿地二值图像中的连通域，将下沉式绿地二值图像中连通域内像素点的灰度值为1的连通域作为下沉式绿地全景灰度图像的初始积水区域，将所述初始积水区域内的其它的连通域作为下沉式绿地全景灰度图像的候选倒影区域，统计所述初始积水区域中像素点的数量，将统计结果作为下沉式绿地全景灰度图像的初始积水面积，Canny边缘检测算法和连通域分析算法的具体实现过程为公知技术，不再进行赘述。

进一步的，根据候选倒影区域所在位置构建候选倒影区域的倒影判定邻域。具体的，如图2所示获取下沉式绿地全景灰度图像中任意一个候选倒影区域的中心，计算所述任意一个候选倒影区域的中心到所述任意一个候选倒影区域的边界像素点的欧氏距离，将所有所述欧氏距离的计算结果按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的邻域判定序列，将所述邻域判定序列中最大值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的长，将所述邻域序列中最小值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的宽，将所述任意一个候选倒影区域的中心作为所述倒影判定邻域的中心，根据所述长、宽、中心确定所述任意一个候选倒影区域的倒影判定领域。

进一步的，在沉式绿地全景灰度图像的任意一个候选倒影区域的倒影判定邻域中构建与所述任意一个候选倒影区域相同形状和大小的滑动搜索窗口，将所述滑动搜索窗口作为所述任意一个候选倒影区域的倒影判定邻域的滑动搜索区域，其中滑动搜索区域的滑动步长为5。

进一步的，根据下沉式绿地全景灰度图像计算每个像素点的倒影纹理熵，通过倒影纹理熵反应下沉式绿地全景灰度图像中每个像素点所在区域的纹理特征，具体的，以下沉式绿地全景灰度图像中任意一个像素点为中心构建预设大小的窗口，将所述预设大小的窗口作为所述任意一个像素点的形状特征分析区域，获取所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵，将所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵的熵分别作为所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值，计算所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值的均值作为所述任意一个像素点的倒影纹理熵，灰度共生矩阵的具体获取过程为公知技术，不再进行赘述。

进一步的，根据候选倒影区域中每个像素点的倒影纹理熵和候选倒影区域的倒影判定邻域中滑动搜索区域内的每个像素点的倒影纹理熵计算倒影相似指数，通过倒影相似指数反应候选倒影区域的倒影判定邻域中是否存在与候选倒影区域相似的滑动搜索区域，即倒影相似指数可以反应候选倒影区域是否为真是倒影区域。具体的倒影相似指数的计算公式如下：

；

式中，表示下沉式绿地全景灰度图像中第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域的倒影相似指数；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域中第/>个像素点的倒影纹理熵；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中第/>个像素点的倒影纹理熵；/>分别表示下沉式绿地全景灰度图像中第/>个候选倒影区域、第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中像素点的数量；/>表示调节参数，大小取经验值0.1。

若下沉式绿地全景灰度图像中第个候选倒影区域可能为真实的倒影，且第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域越相似，则计算得到的/>的值越小，/>的值越大；计算得到的下沉式绿地全景灰度图像中第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域的倒影相似指数/>的值越大，表示下沉式绿地全景灰度图像中第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域相似程度越高，第/>个候选倒影区域真实倒影的可能性越大。

至此，获取了倒影相似指数。

步骤S003，根据倒影相似指数获取倒影相似指数序列，根据倒影相似指数序列计算疑似倒影指数，根据疑似倒影指数计算倒影标记系数，根据倒影标记系数获取下沉式绿地积水区域总面积。

根据下沉式绿地全景灰度图像中每个候选倒影区域与所述每个候选倒影区域的倒影判定领域内滑动搜索区域之间的倒影相似指数获取所述每个候选倒影区域的倒影相似指数序列，根据所述每个候选倒影区域的倒影相似指数计算疑似倒影指数，通过疑似倒影指数反应候选区域为倒影区域的可能性。

具体的，将下沉式绿地全景灰度图像中任意一个候选倒影区域与所述任意一个候选倒影区域的倒影判定邻域内所有滑动搜索区域之间的倒影相似指数按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的倒影相似指数序列。

进一步的，采用最大最小值归一化算法获取下沉式绿地全景灰度图像中任意一个候选倒影区域的倒影相似指数的归一化结果，将所述任意一个候选倒影区域的倒影相似指数的归一化结果作为所述任意一个候选倒影区域的倒影相似指数归一化序列，将所述任意一个候选倒影区域的倒影相似指数归一化序列中的最大值作为所述任意一个候选倒影区域的疑似倒影指数。

进一步的，根据候选倒影区域的疑似倒影指数计算倒影标记系数，根据倒影标记系数计算下沉式绿地全景灰度图像的积水区域总面积。具体的，将下沉式绿地全景灰度图像中每个候选倒影区域的疑似倒影指数与预设阈值进行比较，将疑似倒影指数大于预设阈值的候选倒影区域的倒影标记系数记为1，将疑似倒影指数小于预设阈值的候选倒影区域的倒影标记系数记为0。

进一步的，根据候选倒影区域的倒影标记系数计算下沉式绿地全景灰度图像的积水区总面积，所述积水区域总面积的具体计算公式如下：

；

式中，表示下沉式绿地全景灰度图像的积水区域总面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影标记系数；/>表示下沉式绿地全景灰度图像的初始积水面积；/>表示下沉式绿地全景灰度图像中候选倒影区域的数量。

若下沉式绿地全景灰度图像的积水区域中倒影较多，则计算得到的的值越大，得到的下沉式绿地全景灰度图像的积水区域总面积/>的值越大。

至此，获取了下沉式绿地全景灰度图像的积水区域总面积。

步骤S004，根据下沉式绿地的积水区域总面积获取积水面积-时间曲线，根据积水面积-时间曲线对下沉式绿地的蓄水能力进行评价。

根据下沉式绿地全景灰度图像中积水区域的总面积可以反应在一定时间内下沉式绿地的积水区域面积，因此可在雨后对下沉式绿地公园每隔10分钟采集一幅下沉式绿地全景图像，整体采集时长为24小时。将每幅下沉式绿地全景图像对应的积水区域总面积作为纵坐标，将每幅下沉式绿地全景图像对应的采集时间作为横坐标，将所述横坐标和纵坐标构成的曲线作为所述下沉式绿地的积水面积-时间曲线，如图3所示，计算积水面积-时间曲线下面积作为下沉式绿地的蓄水能力特征值，积水面积-时间曲线下面积的具体计算过程为公知技术，不再进行赘述。

进一步的，获取与下沉式绿地相同环境条件的非下沉式绿地的全景图像，与下沉式绿地的积水面积-时间曲线的获取方式相同，得到非下沉式绿地的蓄水能力特征值，将所述下沉式绿地的蓄水能力特征值与非下沉式绿地的蓄水能力特征值的比值作为所述下沉式绿地的蓄水指数，将所述蓄水指数分别与蓄水能力判断阈值0.8、0.5、0.3进行比较，具体为：当所述蓄水指数小于0.3时，则下沉式绿地的运行效果为优；当所述蓄水指数在区间内时，则下沉式绿地的运行效果为良；当所述蓄水指数大于0.8时，则下沉式绿地的运行效果为差。

至此，完成了对城市公共建筑海绵设施中下沉式绿地的运行效果的评价。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种城市公共建筑海绵设施运行效果评价方法，其特征在于，该方法包括以下步骤：

获取下沉式绿地全景灰度图像；

根据图像分割算法获取下沉式绿地全景灰度图像的下沉式绿地二值图像，基于所述下沉式绿地二值图像获取下沉式绿地全景灰度图像的候选倒影区域；根据下沉式绿地全景灰度图像的候选倒影区域获取所述候选倒影区域的倒影判定邻域，基于所述倒影判定邻域获取所述倒影判定邻域的滑动搜索区域；根据下沉式绿地全景灰度图像计算下沉式绿地全景灰度图像中每个像素点的倒影纹理熵；根据候选倒影区域的倒影判定邻域中滑动搜索区域、候选倒影区域中每个像素点的倒影纹理熵计算候选倒影区域的倒影相似指数；

根据候选倒影区域与候选倒影区域的倒影判定邻域内滑动搜索区域之间的倒影相似指数获取候选倒影区域的倒影相似指数序列；根据候选倒影区域的倒影相似指数序列计算候选倒影区域的疑似倒影指数；根据每个候选倒影区域的疑似倒影指数计算每个候选倒影区域的倒影标记系数；根据每个候选倒影区域的倒影标记系数计算下沉式绿地全景灰度图像的积水区域总面积；

根据下沉式绿地全景灰度图像的积水区域面积获取下沉式绿地的积水面积-时间曲线，基于所述积水面积-时间曲线对下沉式绿地的蓄水能力进行评价。

2.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据图像分割算法获取下沉式绿地全景灰度图像的下沉式绿地二值图像，基于所述下沉式绿地二值图像获取下沉式绿地全景灰度图像的候选倒影区域的方法为：

采用大津阈值分割算法获取下沉式绿地全景灰度图像的二值图像，采用形态学处理对所述二值图像进行腐蚀、膨胀操作，将所述二值图像的腐蚀、膨胀操作处理的结果图像作为所述下沉式绿地全景灰度图像的下沉式绿地二值图像；

采用边缘检测算法和连通域分析算法获取所述下沉式绿地二值图像中的初始积水区域，将所述下沉式绿地二值图像中初始积水区域包含的任意一个其它连通域作为下沉式绿地全景灰度图像的一个候选倒影区域。

3.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据下沉式绿地全景灰度图像的候选倒影区域获取所述候选倒影区域的倒影判定邻域，基于所述倒影判定邻域获取所述倒影判定邻域的滑动搜索区域的方法为：

获取下沉式绿地全景灰度图像中任意一个候选倒影区域的中心，计算所述任意一个候选倒影区域的中心到所述任意一个候选倒影区域的边界像素点的欧氏距离，将所有所述欧氏距离的计算结果按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的邻域判定序列，将所述邻域判定序列中最大值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的长，将所述邻域序列中最小值的8倍作为所述任意一个候选倒影区域的倒影判定邻域的宽，将所述任意一个候选倒影区域的中心作为所述倒影判定邻域的中心，根据所述长、宽、中心确定所述任意一个候选倒影区域的倒影判定领域；

在所述任意一个候选倒影区域的倒影判定邻域中构建与所述任意一个候选倒影区域相同形状和大小的滑动搜索窗口，将所述滑动搜索窗口作为所述任意一个候选倒影区域的倒影判定邻域的滑动搜索区域。

4.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据下沉式绿地全景灰度图像计算下沉式绿地全景灰度图像中每个像素点的倒影纹理熵的方法为：

以下沉式绿地全景灰度图像中任意一个像素点为中心构建预设大小的窗口，将所述预设大小的窗口作为所述任意一个像素点的形状特征分析区域，获取所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵，将所述任意一个像素点的形状特征分析区域的0°、45°、90°、135°方向上的灰度共生矩阵的熵分别作为所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值，计算所述每个像素点的第一特征值、第二特征值、第三特征值、第四特征值的均值作为所述任意一个像素点的倒影纹理熵。

5.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据候选倒影区域的倒影判定邻域的滑动搜索窗口、候选倒影区域中每个像素点的倒影纹理熵计算候选倒影区域的倒影相似指数的方法为：

；

式中，表示下沉式绿地全景灰度图像中第/>个候选倒影区域与第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域的倒影相似指数；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域中第/>个像素点的倒影纹理熵；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中第/>个像素点的倒影纹理熵；/>分别表示下沉式绿地全景灰度图像中第/>个候选倒影区域、第/>个候选倒影区域的倒影判定邻域内第/>个滑动搜索区域中像素点的数量；/>表示调节参数。

6.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据候选倒影区域与候选倒影区域的倒影判定邻域内滑动搜索区域之间的倒影相似指数获取候选倒影区域的倒影相似指数序列的方法为：

将下沉式绿地全景灰度图像中任意一个候选倒影区域与所述任意一个候选倒影区域的倒影判定邻域内所有滑动搜索区域之间的倒影相似指数按照由小到大的顺序排序组成的序列作为所述任意一个候选倒影区域的倒影相似指数序列。

7.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据候选倒影区域的倒影相似指数序列计算候选倒影区域的疑似倒影指数的方法为：

采用归一化算法获取下沉式绿地全景灰度图像中任意一个候选倒影区域的倒影相似指数序列的归一化结果，将所述任意一个候选倒影区域的倒影相似指数序列的归一化结果作为所述任意一个候选倒影区域的倒影相似指数归一化序列，将所述任意一个候选倒影区域的倒影相似指数归一化序列中的最大值作为所述任意一个候选倒影区域的疑似倒影指数。

8.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据每个候选倒影区域的疑似倒影指数计算每个候选倒影区域的倒影标记系数的方法为：

将下沉式绿地全景灰度图像中每个候选倒影区域的疑似倒影指数与预设阈值进行比较，将所述疑似倒影指数大于预设阈值的候选倒影区域的倒影标记系数记为1，将所述疑似倒影指数小于预设阈值的候选倒影区域的倒影标记系数记为0。

9.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据每个候选倒影区域的倒影标记系数计算下沉式绿地全景灰度图像的积水区域总面积的方法为：

；

式中，表示下沉式绿地全景灰度图像的积水区域总面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的面积；/>表示下沉式绿地全景灰度图像中第/>个候选倒影区域的倒影标记系数；/>表示下沉式绿地全景灰度图像的初始积水面积；/>表示下沉式绿地全景灰度图像中候选倒影区域的数量。

10.根据权利要求1所述的一种城市公共建筑海绵设施运行效果评价方法，其特征在于，所述根据下沉式绿地全景灰度图像的积水区域面积获取下沉式绿地的积水面积-时间曲线，基于所述积水面积-时间曲线对下沉式绿地的蓄水能力进行评价的方法为：

在雨后对下沉式绿地公园每隔一段时间采集一幅下沉式绿地全景图像，将每幅下沉式绿地全景图像对应的积水区域总面积作为纵坐标，将每幅下沉式绿地全景图像对应的采集时间作为横坐标，将所述横坐标和纵坐标构成的曲线作为所述下沉式绿地的积水面积-时间曲线，计算积水面积-时间曲线下面积作为下沉式绿地的蓄水能力特征值；

获取与下沉式绿地相同环境条件的非下沉式绿地的全景图像，根据下沉式绿地、非下沉式绿地的全景图像对应的蓄水能力特征值得到下沉式绿地的蓄水能力评价结果。