基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法与流程

【技术领域】

本发明属于城市建成风貌环境评价技术交叉领域，特别是涉及一种基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法。

背景技术：

随着我国城市规划过程中对建筑色彩风貌的重视程度不断提升，城市色彩规划管理调研范围大、研究周期长、对审美要求高、无法客观量化等问题不断呈现出来，为城市色彩规划的进一步深化带来了巨大的挑战。同时，互联网技术不断进步，带来了城市数据量的迅速增长，人工智能与图像处理技术的发展也为图像的量化研究提供了新的可能性。

技术实现要素：

本发明的主要目的在于提供一种基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法，利用深度学习技术、孟塞尔颜色系统全面、准确、高效的采集分析城市建筑物色彩信息、并利用地理空间分析技术对孟塞尔色彩信息进行加工并搭建相关地理信息数据库，对后续城市风貌的规划、建设和管理有重要的参考作用。

本发明通过如下技术方案实现上述目的：一种基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法，其包括以下步骤：

步骤s1：从互联网或线下数据库中获得覆盖目标研究范围的城市街道景观影像数据作为原始图像；

步骤s2：采用优化的动态阈值白平衡算法对所述原始图像进行白平衡处理；

步骤s3：利用语义分割算法提取白平衡处理后图像中的建筑物区域作为研究内容；

步骤s4：利用色阶细节劣化方法提取所述研究内容对应图像的色相、饱和度、灰度特征，作为建筑原始色彩数据；

步骤s5：利用孟塞尔色彩空间体系对所述建筑原始色彩数据进行降维处理，形成包含色相、彩度、明度特征的建筑孟塞尔色彩数据；

步骤s6：将步骤s1中所述目标研究范围划分为50m*50m格网，利用空间分析算法分析每个格网的孟塞尔色相、彩度、明度特征；

步骤s7：利用空间分析算法分析所述建筑孟塞尔色彩数据在所述目标研究范围内的分布特征。

进一步的，所述步骤s2采用优化的动态阈值白平衡算法对步骤s1是中所述的原始图像进行白平衡处理，其包括以下步骤：

s21)将宽度、长度分别w、h的图像为从rgb空间转换到ycrcb空间；

s22)将图像切分为宽高比为4∶3的板块，并计算每个板块的cr、cb的平均值为mr、mb，同时利用公式(1)计算每个板块的cr、cb的方差为dr、db：

其中，cb是指蓝色色度分量，cr是指红色色度分量；n为板块的数量；

s23)设定一个大小为w*h的参考白色点亮度矩阵rl，并利用判别式(2)进行判断，若符合判别式(2)，则作为参考白色点，并将该点(i，j)的亮度复制给rl(i，j)；若不符合，则rl(i，j)的值为0：

其中，rl(i，j)为考白色点亮度矩阵rl第i行第j列元素；

s24)选取rl中最大的10％亮度值中的最小值为lu_min，若rl(i，j)＜lu_min，则rl(i，j)＝0，否则rl(i，j)＝1；

s25)将原图r、g、b值分别与rl相乘，得到r1、g1、b1，并计算r1、g1、b1的平均值为rav、gav、bav，并得到调整增益的公式为：

ymax＝double(max(max(y)))/5

rgain＝ymax/rav

ggain＝ymax/gav

bgain＝ymax/bav

其中，ymax、rgain、ggain、bgain为中间计算量，y是指亮度分量；

s26)调整原图像为：

r0＝r*rgain

g0＝g*ggain

b0＝b*bgain。

进一步的，所述步骤s3中的所述语义分割算法为基于深度学习的语义分割算法，采用deeplabv3+模型。

进一步的，所述步骤s4中所述色阶细节劣化方法，包括将所述步骤s3所得图片分割为100*100的图片板块，使用其中一个随机获得的像素值代替其他像素值并提取该像素值hsv色彩空间中的色相、饱和度、明度特征值，其中色相的取值范围为0-360之间的自然数，饱和度、明度的值域为0-1之间的任意数值。

进一步的，所述步骤s5利用孟塞尔色彩空间体系对步骤s4中的原始色彩数据进行降维处理，其降维算法包括将步骤s4中提取的色相、饱和度、明度特征值进行处理，将色相按照红r、红黄yr、黄y、黄绿gy、绿g、绿蓝bg、蓝b、蓝紫pb、紫p、紫红rp分为10个色相带；将明度按深浅从全黑n0至全灰n5到全白n10划分为10个等级；对于表示色调纯度的彩度，其数值从中间0向外随着色调的纯度增加。

进一步的，所述步骤s6包括将步骤s5中所获得的建筑孟塞尔色彩数据根据其空间位置落入50m*50m的格网中，并分别计算提取每个格网中出现频率最高的色相带、明度饱和度的平均等级作为该格网的建筑色彩数据特征。

进一步的，所述步骤s7中所述空间分析算法为优化的热点工具算法，其包括对数据集中的每一个要素计算getis-ordgi*统计得到的gi^*值，分析数据高值或低值要素在空间上发生聚类的位置，从而实现对区域中饱和度、明度的高值中心和低值中心的计算。

进一步的，所述gi^*值的计算公式为：

其中xx是要素j的属性值，wi,j是要素i和j之间的空间权重，n为要素总数，

与现有技术相比，本发明一种基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法的有益效果在于：利用深度学习技术、孟塞尔颜色系统全面、系统、准确、高效的采集城市建筑物色彩信息，并利用地理空间分析技术对建筑色彩信息进行微观、中观、宏观等不同尺度的空间分析，形成城市建筑色彩地理信息数据库，对后续城市风貌的规划、建设和管理有重要的参考作用。具体的，

1)本发明利用城市街道影像大数据实现了对建筑色彩的提取，能够更加全面、高效、量化的提取城市建筑物的配色信息，大幅降低数据收集、整理的工作量，极大的提升了工作效率；

2)本发明采用了基于严格的人类视觉感知的孟塞尔颜色系统对原始数据进行统计分析，使建筑色彩的提取更加精准、更符合实际的色彩视觉体验；

3)本发明加入了基于格网单元及行政区划单元的色彩特征分析，实现了再微观、中观、宏观等不同空间尺度下的建筑色彩分析，并将其按照空间层级录入gis系统，实现了城市建筑色彩特征数据库的搭建；

4)本发明针对不同区域的色彩构成特征提出包括主要不足及主要提升潜力点等的规划管理提升策略，为决策者提供城市建筑色彩的动态管理及可持续发展所需的信息和系统框架。

【附图说明】

图1为本发明实施例的方法流程图；

图2是本发明实施例中对部分街道景观影像提取建筑物部分的结果图；

图3是本发明实施例中对部分街道景观影像提取原始色彩特征数据；

图4、5为本发明实施例中对部分街道景观影响提取孟塞尔色彩特征数据；

图6为本发明实施例中嘉定区50m*50m格网色相分布图；

图7、8是本发明实施例中嘉定区明度、饱和度冷热点区域分布图。

【具体实施方式】

实施例一：

本实施例为一种基于深度学习与孟塞尔色彩的区域建筑空间特征分析方法，其包括以下步骤：

步骤s1：从互联网或线下数据库中获得覆盖目标研究范围的城市街道景观影像数据作为原始图像；

步骤s2：采用优化的动态阈值白平衡算法对所述原始图像进行白平衡处理；

步骤s3：利用语义分割算法提取白平衡处理后图像中的建筑物区域作为研究内容；本实施例对上海嘉定区进行了应用分析，提取的街道景观影像提取建筑物部分图像如图2所示；

步骤s4：利用色阶细节劣化方法提取所述研究内容对应图像的色相、饱和度、灰度特征，作为建筑原始色彩数据；对部分街道景观影像提取的建筑原始色彩特征数据如图3所示；

步骤s5：利用孟塞尔色彩空间体系对所述建筑原始色彩数据进行降维处理，形成包含色相、彩度、明度特征的建筑孟塞尔色彩数据；对部分街道景观影响提取孟塞尔色彩特征数据如图4、图5所示；

步骤s6：将步骤s1中所述目标研究范围划分为50m*50m格网，利用空间分析算法分析每个格网的孟塞尔色相、彩度、明度特征；嘉定区50m*50m格网色相分布图如图6所示；

步骤s7：利用空间分析算法分析所述建筑孟塞尔色彩数据在所述目标研究范围内的分布特征。

嘉定区明度、饱和度冷热点区域分布图如图7、图8所示。饱和度表示色彩的鲜艳程度；明度表示色彩的明暗程度(黑色的明度最低，白色的明度最高)。

图8中显示出嘉定区建筑色彩饱和度主要以中低饱和度为主，少量高饱和度色彩建筑物融入其中。低饱和度的建筑色彩会给人平衡、舒适的感觉。饱和度空间布局呈现以下的特征：

a、中低饱和度建筑主要分布在嘉定区中部，高饱和度建筑物在嘉定区四周的边缘分布；这些高饱和度的建筑物主要构成了嘉定区的饱和度热点建筑。在建筑群中有突出、强调的作用；

b、与城市功能结构耦合程度高，饱和度较低的城市用地主要功能为居住、商业和工业园区。

图7中显示出嘉定区建筑色彩明度主要以中高明度为主，明度最低的建筑物出现在嘉定区的西部和西北地区。高明度的建筑色彩会给人明亮、轻快的感觉，往往会成为该区域的明度热点。一般来说，城市色彩中保持较高的明度比较容易达到高品质的色彩特征。

所述步骤s2采用优化的动态阈值白平衡算法对步骤s1是中所述的原始图像进行白平衡处理，其包括以下步骤：

s21)将宽度、长度分别w、h的图像为从rgb空间转换到ycrcb空间；

s22)将图像切分为宽高比为4∶3的板块，并计算每个板块的cr、cb的平均值为mr、mb，同时利用公式(1)计算每个板块的cr、cb的方差为dr、db：

其中，cb是指蓝色色度分量，cr是指红色色度分量；n为板块的数量；

s23)设定一个大小为w*h的参考白色点亮度矩阵rl，并利用判别式(2)进行判断，若符合判别式(2)，则作为参考白色点，并将该点(i，j)的亮度复制给rl(i，j)；若不符合，则rl(i，j)的值为0：

其中，rl(i，j)为考白色点亮度矩阵rl第i行第j列元素；

s24)选取rl中最大的10％亮度值中的最小值为lu_min，若rl(i，j)＜lu_min，则rl(i，j)＝0，否则rl(i，j)＝1；

s25)将原图r、g、b值分别与rl相乘，得到r1、g1、b1，并计算r1、g1、b1的平均值为rav、gav、bav，并得到调整增益的公式为：

ymax＝double(max(max(y)))/5

rgain＝ymax/rav

ggain＝ymax/gav

bgain＝ymax/bav

其中，ymax、rgain、ggain、bgain为中间计算量，y是指亮度分量；

s26)调整原图像为：

r0＝r*rgain

g0＝g*ggain

b0＝b*bgain。

所述步骤s3中的所述语义分割算法为基于深度学习的语义分割算法，采用deeplabv3+模型。建筑物识别像素级别准确度可达85.32％，平均iou值为0.2518。

所述步骤s4中所述色阶细节劣化方法，包括：将所述步骤s3所得图片分割为100*100的图片板块，使用其中一个随机获得的像素值代替其他像素值并提取该像素值hsv色彩空间中的色相、饱和度、明度特征值，其中色相的取值范围为0-360之间的自然数，饱和度、明度的值域为0-1之间的任意数值。

所述步骤s5利用孟塞尔色彩空间体系对步骤s4中的原始色彩数据进行降维处理，其降维算法包括将步骤s4中提取的色相、饱和度、明度特征值进行处理，将色相按照红r、红黄yr、黄y、黄绿gy、绿g、绿蓝bg、蓝b、蓝紫pb、紫p、紫红rp分为10个色相带；将明度按深浅从全黑n0至全灰n5到全白n10划分为10个等级；对于表示色调纯度的彩度，其数值从中间0向外随着色调的纯度增加，没有理论上的上限(普通的颜色实际上限为10左右，反光、荧光等材料可高达30)。由于人眼对各种颜色的的敏感度不同，彩度不一定与每个色调和明度组合相匹配。

所述步骤s6包括将步骤s5中所获得的建筑孟塞尔色彩数据根据其空间位置落入50m*50m的格网中，并分别计算提取每个格网中出现频率最高的色相带、明度饱和度的平均等级作为该格网的建筑色彩数据特征。

所述步骤s7中所述空间分析算法为优化的热点工具算法，其包括对数据集中的每一个要素计算getis-ordgi*统计得到的gi^*值(即z值)，分析数据高值或低值要素在空间上发生聚类的位置，从而实现对区域中饱和度、明度的高值中心和低值中心的计算。

所述gi^*值的计算公式为：

其中xj是要素j的属性值，wi，j是要素i和j之间的空间权重，n为要素总数，

对于具有显著统计学意义的正的z得分，z得分越高，高值(热点)的聚类就越紧密。对于统计学上的显著性负z得分，z得分越低，低值(冷点)的聚类就越紧密。

以上所述的仅是本发明的一些实施方式。对于本领域的普通技术人员来说，在不脱离本发明创造构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。