发明名称 一种基于Cloude分解的极化SAR图像分类方法
摘要 本发明公开了一种基于Cloude特征分解的极化SAR图像无监督分类方法,包括:读入一幅待分类的极化SAR图像,对图像中的每个像素点进行Cloude分解,得到散射熵H和散射角α;计算散射熵H和散射角α的统计直方图,并对直方图进行分割,得到分割阈值;将得到的分割阈值作为散射熵H和散射角α特征构成的二维平面的划分点,对极化SAR图像进行初始划分为m类;将得到的初始分类的类中心和类别数输入Wishart分类器,得到所述极化SAR图像的分类结果。本发明对H、α参数进行直方图分割获得划分的阈值,改进H/α分类面的初始类中心,使分类更合理,提高了分类器的性能。
申请公布号 CN104123563B 申请公布日期 2017.04.12
申请号 CN201410341457.9 申请日期 2014.07.17
申请人 电子科技大学 发明人 张扬;蒋霞;郝恩义
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 周永宏
主权项 一种基于Cloude分解的极化SAR合成孔径雷达图像分类方法,其特征是,包括如下步骤:步骤S1:读入一幅待分类的极化SAR图像,对图像中的每个像素点进行Cloude分解,得到散射熵H和散射角α;步骤S2:计算散射熵H和散射角α的统计直方图,并对直方图进行分割,得到分割阈值,步骤S2具体包含以下分步骤:步骤S21:计算所述散射熵H和散射角α的统计直方图;步骤S22:对得到统计直方图进行峰值搜索,得到初始峰值,进一步去掉小峰值和相邻峰值,得到明显峰值,步骤S22还包括以下分步骤:步骤S221:找初始峰值,p<sub>0</sub>={i,s(i)|s(i)≥s(i‑1)&amp;s(i)≥s(i+1),2≤i≤256},p<sub>0</sub>是初始峰值的集合,当i=256时,则s(i+1)=s(1),s(i)表示i点出现的频数,i表示量化的刻度(0~1);步骤S222:去除小峰值,对于集合p<sub>0</sub>,当峰值小于第一阈值时,将这个峰值从集合p<sub>0</sub>中去掉;步骤S223:去除相邻的峰值,当两个峰值之间的距离小于第二阈值时,去除较小的峰值;步骤S23:将得到的明显峰值输入到FCM模糊C均值聚类分割器中,得到分割阈值,步骤S23还包括以下分步骤:步骤S231:初始化直方图的聚类中心<img file="FDA0001170518540000011.GIF" wi="41" he="56" />为步骤S22产生的峰值,设定门限ε,且初始化t=1,m=2,步骤S232:计算模糊隶属度,<img file="FDA0001170518540000012.GIF" wi="563" he="187" />其中,c表示聚类中心的总个数;<img file="FDA0001170518540000013.GIF" wi="38" he="59" />表示第i个聚类中心;<img file="FDA0001170518540000014.GIF" wi="75" he="80" />表示第t次迭代计算时的第i个聚类中心与第j个数据点间的欧几里德距离;<img file="FDA0001170518540000015.GIF" wi="71" he="77" />表示第t次迭代计算时的第i个聚类中心与第k个数据点间的欧几里德距离;<img file="FDA0001170518540000016.GIF" wi="65" he="71" />表示第t次迭代计算的模糊隶属度矩阵中第i 行j列的模糊隶属度,如果<img file="FDA0001170518540000021.GIF" wi="174" he="72" />则<img file="FDA0001170518540000022.GIF" wi="441" he="74" />其中,<img file="FDA0001170518540000023.GIF" wi="291" he="99" />f<sub>j</sub>是映射到直方图的横坐标刻度,1≤j≤256,1≤i≤S,S表示明显峰值个数;步骤S233:计算直方图的聚类中心<img file="FDA0001170518540000024.GIF" wi="130" he="70" /><img file="FDA0001170518540000025.GIF" wi="916" he="214" />其中,h(f<sub>j</sub>)是f<sub>j</sub>的频数,<img file="1.GIF" wi="14" he="26" />为灰度级,1≤j≤256,1≤i≤S,S表示明显峰值个数;步骤S234:如果<img file="FDA0001170518540000026.GIF" wi="326" he="92" />迭代停止,否则t=t+1,重复步骤S232到步骤S234,得到收敛的聚类中心即峰值和隶属度矩阵,根据最大隶属度原则找到分割阈值;步骤S3:将步骤S2得到的阈值作为散射熵H和散射角α特征构成的二维平面的划分点,对极化SAR图像进行初始划分为m类,得到初始分类的类中心,所述类别数m由步骤S2得到的阈值决定;步骤S4:将步骤S3得到的类中心和类别数输入到Wishart分类器中,得到所述极化SAR图像的分类结果。
地址 611731 四川省成都市高新区(西区)西源大道2006号