本发明提供基于图像脑敏感数据的建成环境主导色测度方法和系统，涉及城市品质测度领域。该基于图像脑敏感数据的建成环境主导色测度方法，包括获取建成环境图像样本对应脑电数据；基于脑电数据计算环境主导色的敏感度；根据建成环境图像样本提取主导色特征参数；搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练；将待分析的环境图像输入到训练好的模型中，得到预测的主导色敏感度结果，解决了对融合图像色彩特征与环境品质的非线性模型预测的效果和传统环境图像评价过程繁琐，整体流程周期过长，导致无法同步地对主导色信息进行反馈、影响环境主导色品质预测的精度和效率的问题。
1.一种基于图像脑敏感数据的建成环境主导色测度方法，其特征在于，包括：获取建成环境图像样本对应脑电数据；
搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练；
将待分析的环境图像输入到训练好的模型中，得到预测的主导色敏感度结果；
所述根据建成环境图像样本提取主导色特征参数的步骤具体包括：将图像样本 进行数据维度转换，令图像缩放后的尺寸设为1024×600像素；
对图像进行色彩识别与分割，输出颜色簇划 ，具体为：其中， 表示各个颜色簇畸变程度之和， 表示该像素的颜色值， 表示颜色簇的像素数目， 表示环境图像像素点的坐标， 表示第 类颜色的质心， 表示颜色簇数，为二分量，用于判断 是否属于第 类颜色， 表示第 个颜色簇的像素数目；
根据图像样本色彩提取结果，获取与图像颜色簇关联的所有颜色名称，计算图像样本第 个颜色类别的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；颜色簇色块的边界由像素颜色的平均值计算得到，并经过适度平滑，避免因简化边界造成的测算误差；
构造环境主导色特征参数，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
对环境主导色特征进行min‑max归一化处理，具体为：其中， 表示归一化前的特征值， 表示归一化后所述特征值的结果。
2.根据权利要求1所述的一种基于图像脑敏感数据的建成环境主导色测度方法，其特征在于：所述获取建成环境图像样本对应脑电数据包括对 位被试者在相同的实验室环境下进行 张建成环境图像的脑电数据采集，得到 组脑电数据，每一组数据的数据量均为 ，其中， 是每组数据的主导色特征维度， 是单次采集的脑电数据样本数量。
3.根据权利要求2所述的一种基于图像脑敏感数据的建成环境主导色测度方法，其特征在于：所述基于脑电数据计算环境主导色敏感度的具体步骤包括：选取枕叶区域O1、OZ、O2、POZ、PO3、PO4、PO7、PO8八个导联对于建成环境图像样本刺激的前后3秒脑电信号；
通过原始脑电数据获取图像样本视觉刺激数据的差异波；
对每一导脑电信号做短时傅里叶变换，分别提取预处理后脑电数据的 频带8 13Hz、~频带14 41Hz、 频带4 8Hz功率谱密度；
根据 频带、 频带、 频带的平均相对功率谱计算脑电敏感度指数，从而得到图像样本建成环境主导色敏感度，计算过程如下：其中， 表示脑电敏感度指数，1≤ ≤8，表示所述的八个导联； 、 、分别是导联 、 、 频带的平均相对功率谱；
针对被试者个体差异的影响，对脑敏感度指数进行无量纲处理，具体为：其中， 表示第 位被试者对第 张图像样本无量纲处理后的脑敏感度， 表示图像样本的数量；
4.根据权利要求3所述的一种基于图像脑敏感数据的建成环境主导色测度方法，其特征在于：所述搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练的具体步骤包括：将建成环境图像及其脑敏感数据转化为若干建成环境序列样本，使用XGBoost决策树算法搭建建成环境主导色测度模型，将75%序列样本数据进行训练，其余作为测试集；
利用串联拼接方法融合所述8个维度的环境主导色特征，得到总体环境主导色特征；
将敏感度数据和主导色特征参数输入建成环境主导色测度模型，具体为：其中， 表示第 张图像样本的总体环境主导色特征， 表示该图像样本的主导色敏感度数值， 表示图像样本的数量；
对输入特征参数进行Kaiser‑Meyer‑Olkin检验和Bartlett球面判断。
5.根据权利要求4所述的一种基于图像脑敏感数据的建成环境主导色测度方法，其特征在于：所述将待分析的环境图像输入到训练好的模型中，得到预测的主导色敏感度结果的步骤具体包括：将色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度作为影响指标得到关于建成环境主导色预测的非线性回归模型为：其中， 表示预测的主导色敏感度数据， ，HS表示色相比例，BS表示饱和度，VS表示明度，MCA表示最大颜色簇面积，DPC表示颜色簇分割度，NPC表示颜色簇多样性，IPS表示相似颜色簇蔓延度，CDS表示颜色簇形状复杂度；
针对达到最终训练的权重更加平滑，以避免过拟合现象，采用的损失函数为：其中， 表示所述模型回归树所有预测参数与真实参数之差的集合， 表示测量预测参数与目标参数之差， 表示正则项优化函数，以避免过拟合， 表示所述回归树的叶子结点数， 表示每个叶子结点的得分，和 表示需要调参的系数；
其中， 表示所述建成环境序列样本第 个主导色特征值的平均数， 、 分别表示所有阳性样本和阴性样本特征值的平均数，表示第 个环境主导色特征对应的实例；
计算所述模型的主导色特征权重，并根据特征权重进行环境主导色品质评估，具体处理过程如下：其中， 表示所述建成环境序列样本第 个环境主导色特征的权重值， 表示所述模型回归树叶片的所有样本的梯度统计的总和； 表示所述模型回归树叶片所有样本的二阶统计的总和；
其中， 表示颜色簇的色相总数， 表示主导色特征的权重值， 表示主导色特征的参变量，将 归一化得到最终的环境主导色品质分数。
6.一种基于图像脑敏感数据的建成环境主导色测度系统，其特征在于，所述系统包括：数据采集处理模块，用于获取若干建成环境图像及其对应的脑电数据，转化为若干建成环境序列样本；
脑敏感度提取模块，用于从所述脑电数据中提取脑敏感度指数，得到建成环境主导色的敏感度数值；
主导色特征提取模块，用于从图像样本中进行图像色彩识别与分割，得到图像颜色簇和主导色特征参数；
环境主导色测度模型训练模块，用于搭建建成环境主导色测度模型，将敏感度数据和主导色特征参数输入，利用XGBoost决策树算法进行训练；
特征重要性识别模块，用于识别重要主导色特征，根据环境主导色特征选择表，建立综合环境主导色测度体系；
品质量化评估模块，用于建成环境测度方法应用，根据主导色特征权重进行环境主导色品质评估；
其中，所述用于从图像样本中进行图像色彩识别与分割，得到图像颜色簇和主导色特征参数的步骤具体包括：将图像样本 进行数据维度转换，令图像缩放后的尺寸设为1024×600像素；
对图像进行色彩识别与分割，输出颜色簇划 ，具体为：其中， 表示各个颜色簇畸变程度之和， 表示该像素的颜色值， 表示颜色簇的像素数目， 表示环境图像像素点的坐标， 表示第 类颜色的质心， 表示颜色簇数，为二分量，用于判断 是否属于第 类颜色， 表示第 个颜色簇的像素数目；
根据图像样本色彩提取结果，获取与图像颜色簇关联的所有颜色名称，计算图像样本第 个颜色类别的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；颜色簇色块的边界由像素颜色的平均值计算得到，并经过适度平滑，避免因简化边界造成的测算误差；
构造环境主导色特征参数，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
对环境主导色特征进行min‑max归一化处理，具体为：其中， 表示归一化前的特征值， 表示归一化后所述特征值的结果。
7.根据权利要求6所述的一种基于图像脑敏感数据的建成环境主导色测度系统，其特征在于：所述脑敏感度提取模块具体包括：脑电信号预处理单元，用于对原始脑电数据进行过滤和伪像矫正，并将幅值在10μV‑
100μV区间范围以外的数据作为坏导剔除，根据图像样本进行重分类与叠加平均；
脑电频带提取单元，用于提取八个导联的 频带、 频带、 频带的平均相对功率谱；
敏感度指数计算单元，用于从脑电特征中计算脑敏感度指数；
主导色敏感度获取单元，用于得到建成环境主导色的敏感度数值，作为环境主导色测度模型的训练数据。
8.根据权利要求6所述的一种基于图像脑敏感数据的建成环境主导色测度系统，其特征在于：所述主导色特征提取模块具体包括：样本图像处理单元，用于将所述图像样本进行数据维度转换；
颜色簇提取单元，用于将所述图像样本进行色彩识别与分割，得到图像样本的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；
主导色特征选择单元，用于构造环境主导色特征，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
特征参数计算单元，用于分别计算各主导色特征的参变量；
归一化单元，用于主导色特征能够编码为建成环境主导色测度模型的输入特征，使得环境主导色特征参数落在 区间内。
9.根据权利要求6所述的一种基于图像脑敏感数据的建成环境主导色测度系统，其特征在于：所述环境主导色测度模型训练模块具体包括：环境主导色测度模型构建单元，利用XGBoost决策树算法建立对于环境主导色敏感度和主导色特征的测度模型；
环境主导色测度模型训练单元，用于训练以环境主导色特征为影响指标的非线性回归模型；
环境主导色敏感度预测单元，用于将待预测的建成环境图像数据输入训练好的环境主导色测度模型，得到预测的建成环境主导色敏感度。
基于图像脑敏感数据的建成环境主导色测度方法和系统
技术领域
[0001] 本发明涉及城市品质测度技术领域，具体为基于图像脑敏感数据的建成环境主导色测度方法和系统。
背景技术
[0002] 建成环境主导色是环境品质的重要内容，通常能够有效地揭示建成环境的秩序和功能，被认为是区域定位和空间组织的关键线索，对环境提质增效有着重要作用。由于环境主导色自身具备的易感性和可识别性等特质，环境主导色可以有效运用于城市品质测度任务上，例如，在城市更新保护方面中，通过优化建成环境的主导色系，测量并改善杂乱的城市色彩现状，以此达到较为统一协调的老城主导色风貌；在环境规划设计方面，通过环境主导色品质测度建立环境标准色，及时调整色彩设计问题，提升建成环境的规划管理水平。
[0003] 近年来，为提高环境主导色对于环境品质的评估效果，越来越多的学者将图像数据与主导色预测联系起来。相关知识产权成果有：例如专利申请号为CN202110987218.0的一种基于色彩分布学习的街景图像评分方法，通过图像语义分割、图像实体色彩值计算、实体混合色评估、标签数据训练基于机器学习和色彩分布的街景图像评分方法；例如专利申请号为CN202110893036.7的一种城市街区建筑色彩和谐度评估方法与装置，通过建筑照片样本获取、照片取色分析、色彩区域划分、美景度评价的城市街区建筑色彩和谐度评估方法；例如专利申请号为CN201910833403.7的一种基于人工智算法的城市景观评价指标计算方法通过建立影响因子权重、采集城市评价图片集、景观色彩丰富度评分、因子目标评价函数计算的色彩景观评价指标计算方法。尽管基于图像数据的环境主导色测度方法研究已取得一定进展，但对融合多维度图像主导色特征与环境品质的非线性模型预测效果问题仍需进一步解决。与此同时，传统环境图像评价过程繁琐，整体流程周期过长，人工识别环境主导色品质成本较高，导致无法同步地对主导色信息进行反馈、影响环境主导色品质预测的精度和效率。
[0004] 因此，现有技术方法在进行环境主导色测度分析时，具有基础数据处理分析主观性、随机性，测度模型运行效率、精度及全局性等不足，难以支撑复杂的建成环境主导色测度研究、深入指导建成环境景观主导色品质解析，一种结合图像脑敏感数据进行的建成环境主导色测度方法亟待优化发展并实践应用，有助于实现建成环境主导色与环境品质测度的精确、多维度、整体分析，进而推动城市品质的提质增效。
发明内容
[0005] (一)解决的技术问题
[0006] 针对现有技术的不足，本发明提供了基于图像脑敏感数据的建成环境主导色测度方法和系统，解决了对融合多维度的图像色彩特征与环境品质的非线性模型预测的效果和传统环境图像评价过程繁琐，整体流程周期过长，导致无法同步地对主导色信息进行反馈、影响环境主导色品质预测的精度和效率的问题。
[0007] (二)技术方案
[0008] 为实现以上目的，本发明通过以下技术方案予以实现：
[0009] 一方面，提供了一种基于图像脑敏感数据的建成环境主导色测度方法，所述方法包括：
[0010] 获取建成环境图像样本对应脑电数据；
[0011] 基于脑电数据计算环境主导色的敏感度；
[0012] 根据建成环境图像样本提取主导色特征参数；
[0013] 搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练；
[0014] 将待分析的环境图像输入到训练好的模型中，得到预测的主导色敏感度结果。
[0015] 优选的，所述获取建成环境图像样本对应脑电数据包括对J位被试者在相同的实验室环境下进行I张建成环境图像样本脑电数据采集，得到I*J组脑电数据，每一组数据的(d)数据量均为n ，其中，d是每组数据的主导色特征维度，n是单次采集的脑电数据样本数量。
[0016] 优选的，所述基于脑电数据计算环境主导色的敏感度机体的具体步骤包括：
[0017] 选取枕叶区域O1、OZ、O2、POZ、PO3、PO4、PO7、PO8八个导联对于建成环境图像样本刺激的前后3秒脑电信号；
[0018] 通过原始脑电数据获取样本视觉刺激数据的差异波；
[0019] 对每一导脑电信号做短时傅里叶变换，分别提取预处理后脑电数据的α频带8～13Hz、β频带14～41Hz、θ频带4～8Hz功率谱密度；
[0020] 根据α频带、β频带、θ频带的平均相对功率谱计算脑电敏感度指数，从而得到图像样本建成环境主导色敏感度，计算过程如下：
[0021]
[0022] 其中，EFT表示脑电敏感度指数，1≤k≤8，表示所述的八个导联；Pθ(k)、Pα(k)、Pβ(k)分别是导联α、β、θ频带的平均相对功率谱；
[0023] 针对被试者个体差异的影响，对脑敏感度指数进行无量纲处理，具体为：
[0024]
[0025] 其中，Zj(i)表示第j个被试者对第i张图像样本无量纲处理后的脑敏感度，n表示图像样本的数量；
[0026] 建成环境主导色敏感度数值如下：
[0027]
[0028] 优选的，所述根据建成环境图像样本提取主导色特征参数的步骤具体包括：
[0029] 将样本图像{i1,i2,...,im}进行数据维度转换，令图像缩放后的尺寸设为1024×600像素；
[0030] 对图像进行色彩识别与分割，输出颜色簇划D＝{d1,d2,...,dk}，具体为：
[0031]
[0032]
[0033] 其中，S表示各个颜色簇畸变程度之和，Q(n)表示该像素的颜色值，N表示颜色簇的像素数目，n表示环境图像像素点的坐标，dk表示第k类颜色的质心，K表示颜色簇数，rnk为二分量，用于判断Q(n)是否属于第k类颜色，Tk表示第k个颜色簇的像素数目；
[0034] 根据图像样本色彩提取结果，获取与图像颜色簇关联的所有颜色名称，计算图像样本第{k1,k2,...,kj}个颜色类别的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；颜色簇色块的边界由像素颜色的平均值计算得到，并经过适度平滑，避免因简化边界造成的测算误差；
[0035] 构造环境主导色特征特征参数，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
[0036] 对环境主导色特征进行min‑max归一化处理，具体为：
[0037]
[0038] 其中，Hstd表示归一化前的特征值，Hint表示归一化后所述特征值的结果。
[0039] 优选的，所述搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练的具体步骤包括：
[0040] 将建成环境图像及其脑敏感数据转化为若干建成环境序列样本，使用XGBoost决策树算法搭建建成环境主导色测度模型，将75％建成环境样本数据进行训练，其余作为测试集；
[0041] 利用串联拼接方法融合所述8个维度的环境主导色特征，得到总体环境主导色特征Hall；
[0042] 将敏感度数据和主导色特征参数输入建成环境主导色测度模型，具体为：
[0043] Z＝{(Hi,yi)|i＝1,2,...,n}
[0044] 其中，Hi表示第i张图像样本的总体环境主导色特征，yi表示该图像样本的主导色敏感度数值，n表示图像样本的数量；
[0045] 对输入特征参数进行Kaiser‑Meyer‑Olkin检验和Bartlett球面判断。
[0046] 优选的，所述将待分析的环境图像输入到训练好的模型中，得到预测的主导色敏感度结果的步骤具体包括：
[0047] 将色相比例(HS)、饱和度比例(BS)、明度比例(VS)、最大颜色簇面积(MCA)、颜色簇面积复杂度(NPC)、颜色簇多样性(CDS)、颜色簇分割度(DPC)和相似颜色簇蔓延度(IPS)作为影响指标得到关于建成环境主导色预测的非线性回归模型为：
[0048]
[0049] 其中， 表示预测的主导色敏感度数据，
[0050] 为达到最终训练的权重更加平滑，以避免过拟合现象，采用的损失函数为：
[0051]
[0052] 其中， 表示所述模型回归树所有预测参数与真实参数之差的集合，表示测量预测参数与目标参数之差， 表示正则项优化函数，以避免过拟合，
T表示所述回归树的叶子结点数，ω表示每个叶子结点的得分，η和ρ表示需要调参的系数；
[0053] 计算所述模型的主导色特征重要性分数，具体为：
[0054]
[0055] 其中， 表示所述序列样本第i个主导色特征值的平均数， 分别表示所有阳性样本和阴性样本特征值的平均数，r表示第i个环境主导色特征对应的实例；F(i)越大，表示所述特征对主导色敏感度的影响较大，由此可以筛选重点景观色彩特征，综合建立环境主导色量化体系，以此提环境规划设计品质；
[0056] 计算所述模型的主导色特征权重，并根据特征权重进行环境主导色品质评估，具体处理过程如下：
[0057]
[0058] 其中， 表示所述建成环境序列样本第t个环境主导色特征的权重值， 表示所述模型回归树叶片的所有样本的梯度统计的总和； 表示所述模型回归树叶片所有样本的二阶统计的总和；
[0059] 环境主导色品质分数计算公式为：
[0060] Hquality＝n·w1H1+5·w2H2+4·w3H3+w4H4+w5H5++w6H6+w7H7+w8H8
[0061] 其中，n表示颜色簇的色相总数，w表示主导色特征的权重值，H表示主导色特征的参变量，将Hquality归一化得到最终的环境主导色品质分数。
[0062] 又一方面，提供了一种基于图像脑敏感数据的建成环境主导色测度系统，所述系统包括：
[0063] 数据采集处理模块，用于获取若干建成环境图像对应的脑电数据，转化为若干建成环境图像序列样本；
[0064] 脑敏感度提取模块，用于从所述脑电数据中提取脑敏感度指数，得到建成环境主导色的敏感度数值；
[0065] 主导色特征提取模块，用于从所述图像样本中进行图像色彩识别与分割，得到图像颜色簇和主导色特征参数；
[0066] 环境主导色测度模型训练模块，用于搭建建成环境主导色测度模型，将敏感度数据和主导色特征参数输入，利用XGBoost决策树算法进行训练；
[0067] 特征重要性识别模块，用于识别重要主导色特征，根据环境主导色特征选择表，建立综合环境主导色测度体系；
[0068] 品质量化评估模块，用于建成环境测度方法应用，根据主导色特征权重进行环境主导色品质评估。
[0069] 优选的，所述脑敏感度提取模块具体包括：
[0070] 脑电信号预处理单元，用于对原始脑电数据进行过滤和伪像矫正，并将幅值在10μV‑100μV区间范围以外的数据作为坏导剔除，根据图像样本进行重分类与叠加平均；
[0071] 脑电频带提取单元，用于提取八个导联的α频带、β频带、θ频带的平均相对功率谱；
[0072] 敏感度指数计算单元，用于从脑电特征中计算脑敏感度指数；
[0073] 主导色敏感度获取单元，用于得到建成环境主导色的敏感度数值，作为环境主导色测度模型的训练数据。
[0074] 所述主导色特征提取模块体包括：
[0075] 样本图像处理单元，用于将所述图像序列样本进行数据维度转换；
[0076] 颜色簇提取单元，用于将所述图像序列样本进行色彩识别与分割，得到图像样本的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；
[0077] 主导色特征选择单元，用于构造环境主导色特征，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
[0078] 特征参数计算单元，用于分别计算各主导色特征的参变量；
[0079] 归一化单元，用于主导色特征能够编码为建成环境主导色测度模型的输入特征，使得环境主导色特征参数落在[0,1]区间内。
[0080] 所述环境主导色测度模型训练模块具体包括：
[0081] 环境主导色测度模型构建单元，利用XGBoost决策树算法建立对于环境主导色敏感度和主导色特征的测度模型；
[0082] 特征融合单元，用于加速训练过程；
[0083] 环境主导色测度模型训练单元，用于训练以环境主导色特征为影响指标的非线性回归模型；
[0084] 环境主导色敏感度预测单元，用于将待预测的建成环境图像数据输入训练好的环境主导色测度模型，得到预测的建成环境主导色敏感度。
[0085] (三)有益效果
[0086] 本发明基于图像脑敏感数据的建成环境主导色测度方法和系统，从而有效的解决了对融合图像色彩特征与环境品质的非线性模型预测的效果和传统环境图像评价过程繁琐，整体流程周期过长，导致无法同步地对主导色信息进行反馈、影响环境主导色品质预测的精度和效率的问题。
附图说明
[0087] 图1为本发明基于图像脑敏感数据的建成环境主导色测度方法流程图；
[0088] 图2为本发明实施例的大脑视觉区域和枕叶区域采集导联电极示意图；
[0089] 图3为本发明实施例的建成环境图像样本脑敏感度分布图；
[0090] 图4为本发明实施例的XGBoost决策树方法路线图；
[0091] 图5为本发明实施例的建成环境主导色测度模型的预测结果示意图；
[0092] 图6为本发明实施例测度模型的预测环境敏感度与实际环境敏感度的拟合曲线图。
具体实施方式
[0093] 下面将结合本发明的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0094] 实施例一
[0095] 如图1所示，本实施例提供一种基于图像脑敏感数据的建成环境主导色测度方法，该方法具体包括：
[0096] 获取建成环境图像样本对应脑电数据。对J位被试者在相同的实验室环境下进行I(d)张建成环境图像样本脑电数据采集，得到I*J组脑电数据，每一组数据的数据量均为n ，其中，d是每组数据的主导色特征维度，n是单次采集的脑电数据样本数量。
[0097] 本实施例中，所有被试者按照年龄及性别比例进行选择，其生理与心理状态均为健康、生活环境类似，并签署知情同意书。利用E‑Prime实验操作系统搭建视觉刺激呈现和脑电数据收集系统，每张环境图像显示时长3s并播放三轮，实时采集被试者的原始脑电信号。
[0098] 基于脑电数据计算环境主导色的敏感度，该步骤具体包括：
[0099] 选取枕叶区域O1、OZ、O2、POZ、PO3、PO4、PO7、PO8八个导联(如图2所示)对于建成环境图像样本刺激的前后3秒脑电信号，该区域导联的脑电信号能较好反映环境色彩的视觉信息；
[0100] 为了提高脑电特征提取速度，减少冗余计算，通过原始脑电数据进行电极定位、滤波、独立成分分析、伪迹去除、基线校准和重分类处理，求得样本视觉刺激数据的差异波，差异波幅值大小代表了该图像样本对人的脑敏感度影响力；
[0101] 对每一导脑电信号做短时傅里叶变换，采用汉宁窗进行窗处理，分别提取预处理后脑电数据的α频带(8～13Hz)、β频带(14～41Hz)、θ频带(4～8Hz)功率谱密度；
[0102] 根据α频带、β频带、θ频带的平均相对功率谱计算脑电敏感度指数(FTG)，从而得到图像样本建成环境主导色敏感度(ATD)，计算过程如下：
[0103]
[0104] 式中，EFT表示脑电敏感度指数，1≤k≤8，表示所述的八个导联；Pθ(k)、Pα(k)、Pβ(k)分别是导联α、β、θ频带的平均相对功率谱；
[0105] 针对到被试者个体差异的影响，对脑敏感度指数进行无量纲处理,具体为：
[0106]
[0107] 式中，Zj(i)表示第j位被试者对第i张图像样本无量纲处理后的脑敏感度，n表示图像样本的数量；
[0108] 建成环境主导色敏感度数值如下：
[0109]
[0110] 主导色敏感度用以测量建成环境主导色品质，当敏感度数值EAT≥60，认定主导色敏感度的影响力较强；当40≤EAT＜60，认定主导色敏感度的影响力中等；0＜EAT＜40认定主导色敏感度的影响力较弱。
[0111] 实施例中，原始脑电数据的预处理和频带提取采用eego mylab的asa分析软件包，该软件具有较高的脑电过滤和伪像矫正速度，通过导入所需导联电极位，采用平均电极参考，将幅值在10μV‑100μV区间范围以外的数据作为坏导剔除，去除眼电、肌电的伪迹干扰，根据图像样本进行重分类与叠加平均，进而对其幅度和相位进行分析，抓取其中的α频带、β频带、θ频带平均相对功率谱，最终得到建成环境图像样本的主导色敏感度数值(如图3所示)。
[0112] 根据建成环境图像样本提取主导色特征参数，该步骤具体包括：
[0113] 将样本图像{i1,i2,...,im}进行数据维度转换，令图像缩放后的尺寸设为1024×600像素；
[0114] 对图像进行色彩识别与分割，输出颜色簇划D＝{d1,d2,...,dk}，具体为：
[0115]
[0116]
[0117] 式中，S表示各个颜色簇畸变程度之和，Q(n)表示该像素的颜色值，N表示颜色簇的像素数目，n表示环境图像像素点的坐标，dk表示第k类颜色的质心，K表示颜色簇数，rnk为二分量，用于判断Q(n)是否属于第k类颜色，Tk表示第k个颜色簇的像素数目；
[0118] 根据图像样本色彩提取结果，获取与图像颜色簇关联的所有颜色名称，计算图像样本第{k1,k2,...,kj}个颜色类别的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；颜色簇色块的边界由像素颜色的平均值计算得到，并经过适度平滑，避免因简化边界造成的测算误差；
[0119] 构造环境主导色特征特征参数(表1)，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度。
[0120] 表1：环境主导色特征选择
[0121]
[0122]
[0123]
[0124] 对环境主导色特征进行min‑max归一化处理，具体为：
[0125]
[0126] 其中，Hstd表示归一化前的特征值，Hint表示归一化后所述特征值的结果。
[0127] 本实施例中，令K取[4,6]，得到接近人们视觉空间的颜色簇(如图3所示)即图像样本的主导色，利用误差平方和(SSE)作为评价指标，SSE值越小表示数据点越接近于颜色簇的质心，即样本色彩的提取效果越好。颜色簇饱和度和明度的计算采用OpenCV自带的开源直方图估算器cv2.calcHist；颜色簇像素色块的面积和周长采用Canny边缘检测器，从而计算出每张图像的环境主导色特征参变量。通过归一化处理使得环境主导色特征参数均落在[0,1]区间内，能够编码为建成环境主导色测度模型的输入特征。
[0128] 搭建建成环境主导色测度模型，将敏感度数据和主导色特征作为输入进行训练，该步骤具体包括：
[0129] 将建成环境图像及其脑敏感数据转化为若干建成环境序列样本，使用XGBoost决策树算法搭建建成环境主导色测度模型(如图5所示)，将75％建成环境样本数据进行训练，其余作为测试集；
[0130] 利用串联拼接(Concat)方法融合8个维度的环境主导色特征，得到总体环境主导色特征Hall；
[0131] 将敏感度数据和主导色特征参数输入建成环境主导色测度模型，具体为：
[0132] Z＝{(Hi,yi)|i＝1,2,...,n}
[0133] 其中，Hi表示第i张图像样本的总体环境主导色特征，yi表示该图像样本的主导色敏感度数值，n表示图像样本的数量；
[0134] 对输入特征参数进行Kaiser‑Meyer‑Olkin检验和Bartlett球面判断，若数据结果KMO值＞0.5，同时Bartlett球度检验概率P值＜0.05，则可进行主导色测度模型参数设置。
[0135] 本实施例训练建成环境主导色测度模型时，利用随机搜索算法对决策树算法进行参数寻优，网络参数设置值如表2所示，进而根据模型评估指标进行超参数优化，并利用K‑2
fold交叉验证、决定系数(R)、平均绝对误差(MAE)和均方根误差(RMSE)进一步进行模型评
2
估(表3)，其中R越大表明模型的效果越好，MAE、RMSE的值越小表明模型预测越准确。
[0136] 为控制迭代速率，防止过拟合，利用参数learning_rate以控制迭代速率，同时采用LightGBM算法在保证精确度的前提下加速训练过程。
[0137] 表2：XGBoost决策树算法参数设置
[0138]
[0139] 表3：K‑fold交叉验证下建成环境主导色测度模型性能评估结果
[0140]
[0141]
[0142] 将待分析的环境图像输入到训练好的模型中,得到预测的主导色敏感度结果，该步骤具体包括：
[0143] 将相比例(HS)、饱和度比例(BS)、明度比例(VS)、最大颜色簇面积(MCA)、颜色簇面积复杂度(NPC)、颜色簇多样性(CDS)、颜色簇分割度(DPC)和相似颜色簇蔓延度(IPS)作为影响指标得到关于建成环境主导色预测的非线性回归模型为：
[0144]
[0145] 其中， 表示预测的主导色敏感度数据，
[0146] 为达到最终训练的权重更加平滑，以避免过拟合现象，采用的损失函数为：
[0147]
[0148] 式中， 表示所述模型回归树所有预测参数与真实参数之差的集合，表示测量预测参数与目标参数之差， 表示正则项优化函数，以避免过拟合，
T表示所述回归树的叶子结点数，ω表示每个叶子结点的得分，η和ρ表示需要调参的系数；
[0149] 本实施例中，当融合环境主导色特征时，可以将任一种特征的预测评价结果和多个样本图像主导色敏感度评价结果中输入预设的损失函数，得到环境主导色融合特征的预测结果和评价结果的误差，该误差即为预设的损失函数的输出，从而判断训练模型是否收敛。
[0150] 计算模型的主导色特征重要性分数，具体为：
[0151]
[0152] 式中， 表示所述样本第i个主导色特征值的平均数， 分别表示所有阳性样本和阴性样本特征值的平均数，r表示第i个环境主导色特征对应的实例；F(i)越大，表示特征对主导色敏感度的影响较大，由此可以筛选重点景观色彩特征，综合建立环境主导色量化体系，以此提环境规划设计品质。
[0153] 计算模型的主导色特征权重，并根据特征权重进行环境主导色品质评估，具体处理过程如下：
[0154]
[0155] 式中， 表示样本所述第t个环境主导色特征的权重值， 表示所述模型回归树叶片的所有样本的梯度统计的总和； 表示所述模型回归树叶片所有样本的二阶统计的总和；
[0156] 环境主导色品质分数计算过程如下所示：
[0157] Hquality＝n·w1H1+5·w2H2+4·w3H3+w4H4+w5H5++w6H6+w7H7+w8H8
[0158] 其中，n表示颜色簇的色相总数，w表示主导色特征的权重值，H表示主导色特征的参变量，将Hquality归一化得到最终的环境主导色品质分数。
[0159] 本实施例中，根据模型计算色彩特征重要性分数，其中最大颜色簇面积、颜色簇分割度、色相比例、颜色簇多样性、相似颜色簇蔓延度、饱和度比例、明度比例、颜色簇形状复杂度的分数依次为8486.848、4135.527、3665.604、1270.764、764.674、474.965、440.531、205.862(表4)。因此，结合图5和图6中的分析结果得出，在景观色彩规划设计时，应重点考虑景观色彩的最大颜色簇面积、颜色簇分割度、色相比例以及颜色簇多样性。建成环境样本的部分预测结果如图6所示，其中9.28％的建成环境图像样本的主导色敏感度(ATD)≥70，
78％的图像样本的主导色敏感度(ATD)≥40，均可被认为是获得一定关注的景观。在这样的情况下则基于主导色敏感度的高低(主导色敏感度越大说明配色效果越容易引起人们的关注和兴趣)，由此选出其中相对敏感度低的景观进行更新设计。
[0160] 表4：环境主导色特征权重及重要性分布
[0161]
[0162] 实施例二
[0163] 本实施例，提供了一种基于图像脑敏感数据的建成环境主导色测度系统，该系统包括：
[0164] 数据采集处理模块，用于获取若干建成环境图像对应的脑电数据，转化为若干建成环境图像序列样本；
[0165] 脑敏感度提取模块，用于从脑电数据中提取脑敏感度指数，得到建成环境主导色的敏感度数值；
[0166] 主导色特征提取模块，用于从图像样本中进行图像色彩识别与分割，得到图像颜色簇和主导色特征参数；
[0167] 环境主导色测度模型训练模块，用于搭建建成环境主导色测度模型，将敏感度数据和主导色特征参数输入，利用XGBoost决策树算法进行训练；
[0168] 特征重要性识别模块，用于识别重要主导色特征，根据环境主导色特征选择表，建立综合环境主导色测度体系；
[0169] 品质量化评估模块，用于建成环境测度方法应用，根据主导色特征权重进行环境主导色品质评估。
[0170] 脑敏感度提取模块具体包括：
[0171] 脑电信号预处理单元，用于对原始脑电数据进行过滤和伪像矫正，并将幅值在10μV‑100μV区间范围以外的数据作为坏导剔除，根据图像样本进行重分类与叠加平均；
[0172] 脑电频带提取单元，用于提取八个导联的α频带、β频带、θ频带的平均相对功率谱；
[0173] 敏感度指数计算单元，用于从脑电特征中计算脑敏感度指数；
[0174] 主导色敏感度获取单元，用于得到建成环境主导色的敏感度数值，作为环境主导色测度模型的训练数据。
[0175] 主导色特征提取模块体包括：
[0176] 样本图像处理单元，用于将图像序列样本进行数据维度转换；
[0177] 颜色簇提取单元，用于将图像序列样本进行色彩识别与分割，得到图像样本的饱和度、明度、亮度、通道以及颜色簇色块的面积和周长；
[0178] 主导色特征选择单元，用于构造环境主导色特征，包括色相比例、饱和度比例、明度比例、最大颜色簇面积、颜色簇面积复杂度、颜色簇多样性、颜色簇分割度和相似颜色簇蔓延度；
[0179] 特征参数计算单元，用于分别计算各主导色特征的参变量；
[0180] 归一化单元，用于主导色特征能够编码为建成环境主导色测度模型的输入特征，使得环境主导色特征参数落在[0,1]区间内。
[0181] 环境主导色测度模型训练模块具体包括：
[0182] 环境主导色测度模型构建单元，利用XGBoost决策树算法建立对于环境主导色敏感度和主导色特征的测度模型；
[0183] 特征融合单元，用于加速训练过程；
[0184] 环境主导色测度模型训练单元，用于训练以环境主导色特征为影响指标的非线性回归模型，训练时利用参数learning_rate以控制迭代速率，防止过拟合，同时采用LightGBM算法在保证精确度的前提下加速训练过程；
[0185] 环境主导色敏感度预测单元，用于将待预测的建成环境图像数据输入训练好的环境主导色测度模型，得到预测的建成环境主导色敏感度。
[0186] 需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
