本发明属于道路影像信息提取技术领域,尤其涉及一种高分辨率遥感影像道路信息提取系统。
背景技术:
城市道路交通系统主要由路网系统、交通流和交通监测、控制系统构成。路网系统是交通流运行的载体,交通检测和控制系统指导车流有序通过路网中的冲突点。路网系统和交通监测和控制系统好比沟渠和沟渠之间的控制阀门,通过阀门的控制,控制道路上的交通流有序的通过路网节点,最终使个体车辆到达其行驶的目的地。遥感是指非接触的,远距离的探测技术。一般指运用传感器/遥感器对物体的电磁波的辐射、反射特性的探测。遥感是通过遥感器这类对电磁波敏感的仪器,在远离目标和非接触目标物体条件下探测目标地物。然而,现有城市道路影像提取技术,获取的影像清晰度不高,画面昏暗不清、对比度下降、颜色饱和度下降;同时采集的影像信息处理效不高。
综上,现有技术存在的问题是:现有城市道路影像提取技术,获取的影像清晰度不高,画面昏暗不清、对比度下降、颜色饱和度下降;同时采集的影像信息处理效不高。
技术实现要素:
针对现有技术存在的问题,本发明提供了一种高分辨率遥感影像道路信息提取系统。
本发明是这样实现的,一种高分辨率遥感影像道路信息提取系统,所述高分辨率遥感影像道路信息提取系统包括:
影像采集模块,与数据处理模块连接,用于通过卫星获取遥感图像,并发送给数据处理模块;
所述影响采集模块sar图像变化检测方法包括:
(1)对同一地域不同时间的两幅sar原始图像分别进行去噪处理,得到去噪后图像i1和去噪后图像i2;
(2)分别利用对数比算子和均值比算子,对得到的去噪后图像i1和去噪后图像i2进行运算,得到对数比值差异图xl和均值比值差异图xm,计算公式如下:
其中,μ1,μ2分别表示i1和i2的局部均值;
(3)对得到的对数比值差异图xl和均值比值差异图xm进行特征描述,得到对应的特征描述矩阵xl和特征描述矩阵xm;
(4)对得到的特征描述矩阵xl和特征描述矩阵xm,利用各列向量的欧氏距离度量其相似度,分别构造对应的相似度矩阵,亲和矩阵w1和w2;
(5)对得到的亲和矩阵w1和w2,利用基于亲和矩阵融合谱聚类方法对其进行聚类;
(6)使用k-means聚类方法将输出的特征向量f进行聚类,将得到的像素聚类结果,按照对数比值差异图xl或均值比值差异图xm中对应像素的排列顺序,恢复为原始图像的尺寸,得到最终的变化检测结果图,输出结果;
数据处理模块,与影像采集模块、无线通信模块、显示模块连接,用于对影像采集模块采集的数据进行分析处理;
所述无线通信模块基于用户的位置信息以及当前时间段内统计到的视频请求信息,对用户进行聚类处理包括用户的位置信息用当前的位置坐标来描述:
li=(xi,yi);
其中xi,yi分别表示用户i的横纵坐标值,对于用户i,构建一个内容请求频率向量:
ni=(ni,1,ni,2,...,ni,c);
其中ni,c表示用户i请求内容c的次数,每个用户对应一个内容请求向量,该向量反映了用户的内容请求偏好;
基于用户的位置信息和内容请求偏好信息对用户进行聚类,具有相似内容请求偏好且位置相近的用户分到一个多播组,使用余弦相似度准则来计算两个用户间的相似度,用如下公式计算:
其中β是一个0-1之间的权重系数;
使用k-means聚类方法,对小区内所有的用户d进行聚类,ui={li,ni}表示用户i的聚类信息,聚类的目的是将原始用户分成c类d={d1,…,dc},数学模型上是对下式求最小值:
其中γk为用户群的中心;
所述基于用户的位置信息以及当前时间段内统计到的视频请求信息,对用户进行聚类处理具体步骤如下:
步骤一,从d中随机取c个用户,作为c个用户群的中心;
步骤二,根据相似度的计算公式,计算剩下的用户到c个用户群中心的相似度,将用户划分到相似度最高的用户群;
步骤三,根据聚类结果,更新c个用户群的中心γk={lk,nk},用如下公式:
其中mi是一个0-1之间的权重系数,重复步骤二和步骤三,直到聚类中心不再发生变化;
所述根据用户聚类结果,根据每个用户群的位置信息,计算出每个用户群中心位置的水平方位角和垂直仰角具体包括:
采用有源天线波束赋形模型,基站对每个用户群有一个特定波束,即对每个用户群设置一个特定的电子下倾角和垂直半功率带宽的波束,基站坐标为原点o(0,0,hbs),用户群k的质心为γk,位置坐标为(xk,yk,zk),垂直仰角和水平方位角为
基于聚类后的用户群位置信息,用户群质心的水平方向角和垂直仰角通过下面的公式求出:
显然,垂直仰角和水平方位角的取值范围为θ1∈(0,π),
所述基站天线波束实现对用户群的精确对准具体包括:
步骤一,将调整波束的电子下倾角、电子水平角和半功率带宽,使波束的辐射方向对准用户群的中心位置,使半功率带宽范围覆盖用户群中的所有用户,基站到用户的下倾角和水平角将调整为:
其中,
步骤二,确定波束宽度,用户群的覆盖区域为一个圆点在用户群中心的圆形,则该圆形区域的半径为该用户群中离中心位置最远的用户与中心的距离,即:
其中(xk,yk)为用户群k的中心γk的坐标,则第k个波束的垂直半功率带宽为:
所述采用有源天线阵列的天线模型,并确定基站到用户的信道增益模型具体包括:
步骤一,根据每个用户的位置信息以及所属的用户群波束,计算出每个用户位置的实际水平方位角和垂直仰角,计算出用户i相对于基站的水平方位角
步骤二,有源天线阵列的天线模型:
3d天线增益模型采用3gpp标准中提出的有源天线阵列辐射模型,天线增益模型表示如下:
其中,
其中,θetilt表示天线波束的下倾角,
步骤三,基站到用户的信道增益模型,采用多播信道增益模型,在一个多播组中的用户以相同的速率接受数据,基站的传输速率超过了该群中的某个用户的最大承受速率,则这个用户将不能正常解码该数据,基站以用户群中最小的速率传输数据,因此用户群k中基站到用户的等效信道增益等于该用户群中用户的最差信道增益,即:
其中
其中,f和pl分别表示快衰落和路径损耗,
所述提出用户群分簇算法,根据用户群的位置信息,对用户群进行分簇处理具体包括:
基于图论的知识对用户群进行分簇,定义波束间的干扰图g=(v,e),其中v表示波束的集合,作为干扰图的顶点,e表示波束间的干扰系数,作为干扰图的边,定义指示函数e(vk,vm)(k≠m)指示波束k和波束m间的干扰:
其中ok和om分别表示用户群k和用户群m的半径,rth表示两个波束间干扰忽悠不计的门限距离,另外,定义e(vk,vk)=0,表示波束自身不存在干扰,根据指示函数,构建一个二值干扰矩阵:
定义波束的干扰度:
当dg(vk)=0时,称vk为零度节点;
分簇的具体步骤如下:
步骤一,用顶点集合v构建干扰矩阵ag,初始化迭代因子h=1,孤立节点集合
步骤二,找到所有的零度节点vk,更新s=s∪vk;剩余节点集合记为φ1=v-s;
步骤三,分簇:a)
步骤四,用节点集合bh重新构建ag≠0,更新节点集合φh+1=bh,更新迭代因子h=h+1,执行步骤(3);如果ag=0或者|bh|=1,如果|bh|=1,则φh+1=bh;
步骤五,将孤立节点集合s分配到最少节点的一簇中;
经过用户群的分簇处理后,用户群d={d1,…,dk,…,dc}经过分簇算法被划分为φ={φ1,…,φh,…},φh表示第h个用户群簇,每个簇中的总的用户传输速率为:
系统总的吞吐量为所有用户群簇的传输速率之和:
其中
条件(2)表示一个载波只能分配给一个用户群簇,同一簇中的用户群共享一个载波资源,不同簇中的用户群不可以复用;
所述基于最大化吞吐量的载波分配算法具体步骤如下:
步骤一,根据公式:
计算每个簇中的用户在载波n上的总传输速率;
步骤二,为了最大化系统的吞吐量,找到获得最大速率的载波和用户群簇,首先分配该载波给该用户群簇,根据公式:
将载波n分配给用户群簇φh获得最大的传输速率,载波n分配给簇φh的频谱利用率最高,所以将载波n分配给用户群簇φh;
步骤三,将载波n从载波集合f中移除,同时,将用户群簇φh从集合φ中移除;
步骤四,重复执行步骤二和步骤三,直到载波集合或者用户群簇集合成为空集;
无线通信模块,与数据处理模块连接,与无线基站无线连接,用于将数据处理模块的处理的数据信息通过无线形式发送给无线基站,无线基站将数据信息发送给远程服务器进行存储;
所述无线基站宽带无线全双工mimo通信系统回波自干扰自适应抑制方法将迭代滤波法的迭代思想运用到回波自干扰抑制方法中,通过自适应迭代算法实现对自干扰信号的准确估计及系统误码率ber性能的提升;具体步骤如下:
步骤一:近端通信节点的接收信号
tr(n)=hfe(n)tfe(n)+hne(n)tne(n)+w(n);
其中,
步骤二:在接收端通过利用归一化最小均方误差nlms算法对混有自干扰、信道噪声的接收信号进行自干扰抑制,定义算法的代价函数为:
其中,min表示取最小值,n表示第n时刻,e(ene(n))2表示近端误差信号
步骤三:设置采用归一化最小均方误差nlms算法进行自干扰抑制的相关初始值:
令初始迭代次数k=1,并设置最大迭代次数k及根据近端输入信号的自相关矩阵设置收敛步长因子μne,自适应滤波器的初始化权值矢量αne(0)以及滤波器的长度m,开始迭代过程,分别设置k=25、m=11、
步骤四:根据公式
其中j=n1,…,nt,nt表示发送天线总数目,m为自适应滤波器的长度,αne(n)在表示n时刻的权值矢量,
j<nt,则令j=j+1,估计下一接收天线上的估计信号
j=nt,则前进至步骤五;
步骤五:根据下式,更新n时刻的权值矢量
如下式更新下一时刻的权值矢量:
其中,j=1,…,nt,nt表示发送天线总数目,ε表示的是权值矢量αne(n)在迭代过程中的调整因数,仿真中设置的大小为0.001,
步骤六:根据最佳权值矢量αne(n)以及公式:
由下式得近端估计信号
其中,j=1,…,nt,nt表示发送天线总数目,αne(n)表示n时刻权值矢量,
步骤七:从总体接收信号tr(n)中滤除估计出的回波自干扰信号,以获得来自远端节点的有用传输信号,将该信号送入后续的mimo译码检测单元,以获得对远端发送信号的准确估计,具体包括:
第一步,从接收信号tr(n)中减去回波自干扰估计信号
第二步,将信号tes(n)送入后续的mimo译码检测单元,以获得对远端发送信号tfe(n)的估计;
显示模块,与数据处理模块连接,用于显示高清的遥感影像。
进一步,所述数据处理模块,包括特征提取模块、特征分类模块、图像高清模块、图像转换模块;
特征提取模块模块,用于对遥感图像特征信息进行提炼;
特征分类模块,用于对遥感图像特征信息进行分类;
图像高清模块,用于对遥感图像进行高清处理;
图像转换模块,用于对遥感图像进行尺寸大小处理。
进一步,所述图像高清模块图像处理方法如下:
首先,获取原始图像以及原始图像加入环境光底色后的降质图像;
然后,获取原始图像和降质图像的色度差异;将原始图像和降质图像从第一颜色空间转换到第二颜色空间;以及获取原始图像和降质图像在第二颜色空间中的色度差异;
最后,利用色度差异对原始图像进行色度补偿以得到补偿后图像。
本发明通过图像高清模块对原始图像进行色度补偿例如对颜色饱和度进行补偿或者对原始图像色域范围进行扩展,因此能够提高图像色彩鲜艳度,显示设备的显示效果大大提高;通过特征提取模块模块和特征分类模块可以对图像进行提炼、分类,大大提高采集影像的处理效率。
附图说明
图1是本发明实施例提供的高分辨率遥感影像道路信息提取系统结构示意图;
图中:1、影像采集模块;2、数据处理模块;3、无线通信模块;4、显示模块;5、无线基站;6、远程服务器。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下。
下面结合附图对本发明的结构作详细的描述。
如图1所示,本发明实施例提供的高分辨率遥感影像道路信息提取系统包括:影像采集模块1、数据处理模块2、无线通信模块3、显示模块4、无线基站5、远程服务器6。
影像采集模块1,与数据处理模块2连接,用于通过卫星获取遥感图像,并发送给数据处理模块2;
数据处理模块2,与影像采集模块1、无线通信模块3、显示模块4连接,用于对影像采集模块1采集的数据进行分析处理;
无线通信模块3,与数据处理模块2连接,与无线基站5无线连接,用于将数据处理模块2的处理的数据信息通过无线形式发送给无线基站5,无线基站5将数据信息发送给远程服务器6进行存储;
显示模块4,与数据处理模块2连接,用于显示高清的遥感影像。
高分辨率遥感影像道路信息提取系统,其特征在于,所述高分辨率遥感影像道路信息提取系统包括:
影像采集模块,与数据处理模块连接,用于通过卫星获取遥感图像,并发送给数据处理模块;
所述影响采集模块sar图像变化检测方法包括:
(1)对同一地域不同时间的两幅sar原始图像分别进行去噪处理,得到去噪后图像i1和去噪后图像i2;
(2)分别利用对数比算子和均值比算子,对得到的去噪后图像i1和去噪后图像i2进行运算,得到对数比值差异图xl和均值比值差异图xm,计算公式如下:
其中,μ1,μ2分别表示i1和i2的局部均值;
(3)对得到的对数比值差异图xl和均值比值差异图xm进行特征描述,得到对应的特征描述矩阵xl和特征描述矩阵xm;
(4)对得到的特征描述矩阵xl和特征描述矩阵xm,利用各列向量的欧氏距离度量其相似度,分别构造对应的相似度矩阵,亲和矩阵w1和w2;
(5)对得到的亲和矩阵w1和w2,利用基于亲和矩阵融合谱聚类方法对其进行聚类;
(6)使用k-means聚类方法将输出的特征向量f进行聚类,将得到的像素聚类结果,按照对数比值差异图xl或均值比值差异图xm中对应像素的排列顺序,恢复为原始图像的尺寸,得到最终的变化检测结果图,输出结果;
所述无线通信模块基于用户的位置信息以及当前时间段内统计到的视频请求信息,对用户进行聚类处理包括用户的位置信息用当前的位置坐标来描述:
li=(xi,yi);
其中xi,yi分别表示用户i的横纵坐标值,对于用户i,构建一个内容请求频率向量:
ni=(ni,1,ni,2,...,ni,c);
其中ni,c表示用户i请求内容c的次数,每个用户对应一个内容请求向量,该向量反映了用户的内容请求偏好;
基于用户的位置信息和内容请求偏好信息对用户进行聚类,具有相似内容请求偏好且位置相近的用户分到一个多播组,使用余弦相似度准则来计算两个用户间的相似度,用如下公式计算:
其中β是一个0-1之间的权重系数;
使用k-means聚类方法,对小区内所有的用户d进行聚类,ui={li,ni}表示用户i的聚类信息,聚类的目的是将原始用户分成c类d={d1,…,dc},数学模型上是对下式求最小值:
其中γk为用户群的中心;
所述基于用户的位置信息以及当前时间段内统计到的视频请求信息,对用户进行聚类处理具体步骤如下:
步骤一,从d中随机取c个用户,作为c个用户群的中心;
步骤二,根据相似度的计算公式,计算剩下的用户到c个用户群中心的相似度,将用户划分到相似度最高的用户群;
步骤三,根据聚类结果,更新c个用户群的中心γk={lk,nk},用如下公式:
其中mi是一个0-1之间的权重系数,重复步骤二和步骤三,直到聚类中心不再发生变化;
所述根据用户聚类结果,根据每个用户群的位置信息,计算出每个用户群中心位置的水平方位角和垂直仰角具体包括:
采用有源天线波束赋形模型,基站对每个用户群有一个特定波束,即对每个用户群设置一个特定的电子下倾角和垂直半功率带宽的波束,基站坐标为原点o(0,0,hbs),用户群k的质心为γk,位置坐标为(xk,yk,zk),垂直仰角和水平方位角为
基于聚类后的用户群位置信息,用户群质心的水平方向角和垂直仰角通过下面的公式求出:
显然,垂直仰角和水平方位角的取值范围为θ1∈(0,π),
所述基站天线波束实现对用户群的精确对准具体包括:
步骤一,将调整波束的电子下倾角、电子水平角和半功率带宽,使波束的辐射方向对准用户群的中心位置,使半功率带宽范围覆盖用户群中的所有用户,基站到用户的下倾角和水平角将调整为:
其中,
步骤二,确定波束宽度,用户群的覆盖区域为一个圆点在用户群中心的圆形,则该圆形区域的半径为该用户群中离中心位置最远的用户与中心的距离,即:
其中(xk,yk)为用户群k的中心γk的坐标,则第k个波束的垂直半功率带宽为:
所述采用有源天线阵列的天线模型,并确定基站到用户的信道增益模型具体包括:
步骤一,根据每个用户的位置信息以及所属的用户群波束,计算出每个用户位置的实际水平方位角和垂直仰角,计算出用户i相对于基站的水平方位角
步骤二,有源天线阵列的天线模型:
3d天线增益模型采用3gpp标准中提出的有源天线阵列辐射模型,天线增益模型表示如下:
其中,
其中,θetilt表示天线波束的下倾角,
步骤三,基站到用户的信道增益模型,采用多播信道增益模型,在一个多播组中的用户以相同的速率接受数据,基站的传输速率超过了该群中的某个用户的最大承受速率,则这个用户将不能正常解码该数据,基站以用户群中最小的速率传输数据,因此用户群k中基站到用户的等效信道增益等于该用户群中用户的最差信道增益,即:
其中
其中,f和pl分别表示快衰落和路径损耗,
所述提出用户群分簇算法,根据用户群的位置信息,对用户群进行分簇处理具体包括:
基于图论的知识对用户群进行分簇,定义波束间的干扰图g=(v,e),其中v表示波束的集合,作为干扰图的顶点,e表示波束间的干扰系数,作为干扰图的边,定义指示函数e(vk,vm)(k≠m)指示波束k和波束m间的干扰:
其中ok和om分别表示用户群k和用户群m的半径,rth表示两个波束间干扰忽悠不计的门限距离,另外,定义e(vk,vk)=0,表示波束自身不存在干扰,根据指示函数,构建一个二值干扰矩阵:
定义波束的干扰度:
当dg(vk)=0时,称vk为零度节点;
分簇的具体步骤如下:
步骤一,用顶点集合v构建干扰矩阵ag,初始化迭代因子h=1,孤立节点集合
步骤二,找到所有的零度节点vk,更新s=s∪vk;剩余节点集合记为φ1=v-s;
步骤三,分簇:a)
步骤四,用节点集合bh重新构建ag≠0,更新节点集合φh+1=bh,更新迭代因子h=h+1,执行步骤(3);如果ag=0或者|bh|=1,如果|bh|=1,则φh+1=bh;
步骤五,将孤立节点集合s分配到最少节点的一簇中;
经过用户群的分簇处理后,用户群d={d1,…,dk,…,dc}经过分簇算法被划分为φ={φ1,…,φh,…},φh表示第h个用户群簇,每个簇中的总的用户传输速率为:
系统总的吞吐量为所有用户群簇的传输速率之和:
其中
条件(2)表示一个载波只能分配给一个用户群簇,同一簇中的用户群共享一个载波资源,不同簇中的用户群不可以复用;
所述基于最大化吞吐量的载波分配算法具体步骤如下:
步骤一,根据公式:
计算每个簇中的用户在载波n上的总传输速率;
步骤二,为了最大化系统的吞吐量,找到获得最大速率的载波和用户群簇,首先分配该载波给该用户群簇,根据公式:
将载波n分配给用户群簇φh获得最大的传输速率,载波n分配给簇φh的频谱利用率最高,所以将载波n分配给用户群簇φh;
步骤三,将载波n从载波集合f中移除,同时,将用户群簇φh从集合φ中移除;
步骤四,重复执行步骤二和步骤三,直到载波集合或者用户群簇集合成为空集。
所述无线基站宽带无线全双工mimo通信系统回波自干扰自适应抑制方法将迭代滤波法的迭代思想运用到回波自干扰抑制方法中,通过自适应迭代算法实现对自干扰信号的准确估计及系统误码率ber性能的提升;具体步骤如下:
步骤一:近端通信节点的接收信号
tr(n)=hfe(n)tfe(n)+hne(n)tne(n)+w(n);
其中,
步骤二:在接收端通过利用归一化最小均方误差nlms算法对混有自干扰、信道噪声的接收信号进行自干扰抑制,定义算法的代价函数为:
其中,min表示取最小值,n表示第n时刻,e(ene(n))2表示近端误差信号
步骤三:设置采用归一化最小均方误差nlms算法进行自干扰抑制的相关初始值:
令初始迭代次数k=1,并设置最大迭代次数k及根据近端输入信号的自相关矩阵设置收敛步长因子μne,自适应滤波器的初始化权值矢量αne(0)以及滤波器的长度m,开始迭代过程,分别设置k=25、m=11、
步骤四:根据公式
其中j=n1,…,nt,nt表示发送天线总数目,m为自适应滤波器的长度,αne(n)在表示n时刻的权值矢量,
j<nt,则令j=j+1,估计下一接收天线上的估计信号
j=nt,则前进至步骤五;
步骤五:根据下式,更新n时刻的权值矢量
如下式更新下一时刻的权值矢量:
其中,j=1,…,nt,nt表示发送天线总数目,ε表示的是权值矢量αne(n)在迭代过程中的调整因数,仿真中设置的大小为0.001,
步骤六:根据最佳权值矢量αne(n)以及公式:
由下式得近端估计信号
其中,j=1,…,nt,nt表示发送天线总数目,αne(n)表示n时刻权值矢量,
步骤七:从总体接收信号tr(n)中滤除估计出的回波自干扰信号,以获得来自远端节点的有用传输信号,将该信号送入后续的mimo译码检测单元,以获得对远端发送信号的准确估计,具体包括:
第一步,从接收信号tr(n)中减去回波自干扰估计信号
第二步,将信号tes(n)送入后续的mimo译码检测单元,以获得对远端发送信号tfe(n)的估计。
本发明提供的数据处理模块2,包括特征提取模块、特征分类模块、图像高清模块、图像转换模块;
特征提取模块模块,用于对遥感图像特征信息进行提炼;
特征分类模块,用于对遥感图像特征信息进行分类;
图像高清模块,用于对遥感图像进行高清处理;
图像转换模块,用于对遥感图像进行尺寸大小处理。
本发明提供的图像高清模块图像处理方法如下:
首先,获取原始图像以及原始图像加入环境光底色后的降质图像;
然后,获取原始图像和降质图像的色度差异;将原始图像和降质图像从第一颜色空间转换到第二颜色空间;以及获取原始图像和降质图像在第二颜色空间中的色度差异;
最后,利用色度差异对原始图像进行色度补偿以得到补偿后图像。
本发明通过影像采集模块1取遥感图像,并发送给数据处理模块2,数据处理模块2对影像采集模块1采集的数据进行分析处理;然后由无线通信模块3将数据处理模块2的处理的数据信息通过无线形式发送给无线基站5,无线基站5将数据信息发送给远程服务器6进行存储,通过显示模块4可以显示高清的遥感影像。
以上所述仅是对本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰,均属于本发明技术方案的范围内。