发明名称 一种用于宽带光谱成像系统的最优化训练样本选择方法
摘要 本发明公开了一种用于宽带光谱成像系统的最优化训练样本选择方法,包括搭建实际宽带光谱成像系统,标定光谱成像系统的各项参数;准备光谱成像常用样本集并测量其光谱反射率信息;利用实际光谱成像系统拍摄样本集,依据系统标定信息对样本集数字响应值校正;利用伪逆光谱重建方法重建样本集,以光谱均方根误差最小化为原则,从样本集中逐个挑选训练样本组成训练样本集;每挑选一个训练样本加入到训练样本集时,将样本集中所有样本遍历一次,直至训练样本集对总样本集的光谱重建误差开始收敛于某最小值Jmin时,锁定最优训练样本集,实现宽带光谱成像最优训练样本选择。该方法具有显著降低光谱成像复杂性、提高光谱重建精度等技术优势。
申请公布号 CN104359556A 申请公布日期 2015.02.18
申请号 CN201410653613.5 申请日期 2014.11.14
申请人 武汉大学 发明人 万晓霞;梁金星;刘强;朱时良;李焕;李俊锋
分类号 G01J3/28(2006.01)I;G06F19/00(2011.01)I 主分类号 G01J3/28(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 薛玲
主权项 一种用于宽带光谱成像系统的最优化训练样本选择方法,其特征在于,包括以下步骤:步骤1:搭建实际宽带光谱成像系统;步骤2:对光谱成像系统的暗电流水平、光源稳定性、光照不均匀性、各通道线性化响应程度进行标定;步骤3:准备光谱成像常用的样本集Θ,包括标准色靶或自制色靶,利用专业测量设备测量获取样本集Θ的光谱反射率信息r;步骤4:用标定后的真实光谱成像系统对样本集Θ进行拍摄,并根据步骤2对光谱成像系统的标定结果对拍摄的样本集进行校正,获得颜色样本集校正后的数字响应信号值d;步骤5:利用伪逆方法作为训练样本选择时的光谱重建方法,从选择第一个训练样本开始,针对样本集Θ中的每一个样本r<sub>i</sub>,利用伪逆方法计算得到样本对应的转换矩阵T;然后对样本集Θ进行光谱重建;以光谱重建误差(root mean square,RMS)作为评价指标,挑选出光谱重建误差最小的一个样本作为第一个训练样本s<sub>1</sub>;此时,训练样本集中的第一个样本Ω<sub>1</sub>被选定;以此类推,在挑选剩余的第2个到第k个训练样本;所述的转换矩阵TWEI:T=r<sub>i</sub>·PINV(d<sub>i</sub>)   (式一);式一中,d<sub>i</sub>表示与r<sub>i</sub>对应的第i个样本的系统响应信号值;所述的对样本集Θ进行光谱重建的公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>r</mi><mo>^</mo></mover><mo>=</mo><mi>T</mi><mo>&CenterDot;</mo><mi>d</mi></mrow>]]></math><img file="FDA0000608375000000011.GIF" wi="189" he="85" /></maths>   (式二);<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mi>RMS</mi><mo>=</mo><mi>E</mi><mo>{</mo><mo>|</mo><mo>|</mo><mover><mi>r</mi><mo>^</mo></mover><mo>-</mo><mi>r</mi><mo>|</mo><mo>|</mo><mo>}</mo></mrow>]]></math><img file="FDA0000608375000000012.GIF" wi="483" he="166" /></maths>   (式三);<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>s</mi><mn>1</mn></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><msub><mi>r</mi><mi>i</mi></msub><mo>&Element;</mo><mi>&Theta;</mi></mrow></munder><mi>J</mi></mrow>]]></math><img file="FDA0000608375000000013.GIF" wi="315" he="109" /></maths>   (式四);Ω<sub>1</sub>={r<sub>1</sub>}   (式五);所述的挑选剩余的第2个到第k个训练样本的计算公式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>T</mi><mo>=</mo><mo>{</mo><msub><mi>&Omega;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&cup;</mo><msub><mi>r</mi><mi>i</mi></msub><mo>}</mo><mo>&CenterDot;</mo><mi>PINV</mi><mrow><mo>(</mo><msub><mi>d</mi><mrow><mo>{</mo><msub><mi>&Omega;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&cup;</mo><msub><mi>r</mi><mi>i</mi></msub><mo>}</mo></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mover><mi>r</mi><mo>^</mo></mover><mo>=</mo><mi>T</mi><mo>&CenterDot;</mo><mi>d</mi></mtd></mtr><mtr><mtd><mi>J</mi><mo>=</mo><mi>E</mi><mo>{</mo><mo>|</mo><mo>|</mo><mover><mi>r</mi><mo>^</mo></mover><mo>-</mo><mi>r</mi><mo>|</mo><mo>|</mo><mo>}</mo></mtd></mtr><mtr><mtd><msub><mi>s</mi><mi>k</mi></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><msub><mi>r</mi><mi>i</mi></msub><mo>&Element;</mo><mi>&Theta;</mi></mrow></munder><mi>J</mi></mtd></mtr><mtr><mtd><msub><mi>&Omega;</mi><mi>k</mi></msub><mo>=</mo><msub><mi>&Omega;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&cup;</mo><mo>{</mo><msub><mi>s</mi><mi>k</mi></msub><mo>}</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000608375000000021.GIF" wi="717" he="589" /></maths>   (式六);步骤6:针对步骤5中的训练样本选择方法,直到满足式七所表示的条件时,完成训练样本选择,此时确定的训练样本集Ω<sub>k</sub>对样本集Θ光谱重建误差开始达到最小,开始收敛;当继续增加满足式六的训练样本s<sub>k+1</sub>,s<sub>k+2</sub>,…,s<sub>m</sub>时,训练样本集Ω<sub>k+1</sub>,Ω<sub>k+2</sub>,…,Ω<sub>m</sub>对样本集Θ的光谱重建误差仅有极其微弱减小趋势,基本收敛于最小值J<sub>min</sub>;其中:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munder><mrow><mi>arg</mi><mi>min</mi></mrow><msub><mi>s</mi><mi>k</mi></msub></munder><msub><mi>J</mi><mi>k</mi></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><msub><mi>s</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></munder><msub><mi>J</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><msub><mi>s</mi><mrow><mi>k</mi><mo>+</mo><mn>2</mn></mrow></msub></munder><msub><mi>J</mi><mrow><mi>k</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>=</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>=</mo><msub><mi>J</mi><mi>min</mi></msub></mrow>]]></math><img file="FDA0000608375000000022.GIF" wi="1186" he="118" /></maths>   (式七);每挑选一个新的样本加入到训练样本集Ω时,都要将样本集Θ中所有的样本遍历一次,其中包括已经存在于训练样本集Ω中的样本;步骤7:锁定满足步骤6中式七所示的规则和条件的最优训练样本集Ω<sub>k</sub>,完成宽带光谱成像系统最优训练样本选择。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学