主权项 |
一种用于宽带光谱成像系统的最优化训练样本选择方法,其特征在于,包括以下步骤:步骤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>·</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>∈</mo><mi>Θ</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>Ω</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>∪</mo><msub><mi>r</mi><mi>i</mi></msub><mo>}</mo><mo>·</mo><mi>PINV</mi><mrow><mo>(</mo><msub><mi>d</mi><mrow><mo>{</mo><msub><mi>Ω</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>∪</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>·</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>∈</mo><mi>Θ</mi></mrow></munder><mi>J</mi></mtd></mtr><mtr><mtd><msub><mi>Ω</mi><mi>k</mi></msub><mo>=</mo><msub><mi>Ω</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>∪</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>,完成宽带光谱成像系统最优训练样本选择。 |