发明名称 基于二维压缩感知的遥感影像高效采集和增量更新方法
摘要 基于二维压缩感知的遥感影像高效采集和增量更新方法,属于遥感和压缩感知技术领域,提供了一种将传统影像采集方法和压缩感知采集方法融为一体的遥感影像高效采集和增量更新方法。本发明所述的方法:首先,以传统方法采集时相t1的遥感影像;然后,以压缩感知方法采集时相t1和时相t2的测量数据;再以时相t1和时相t2测量数据的差值,基于二维压缩感知模型,利用变化区域的时空连续性先验信息,高效重构二维稀疏信号;最后,以重构结果和时相t1的遥感影像实现时相t2的遥感影像增量更新。更新时相t2的遥感影像仅需2-3倍于变化区域的测量数据,这是传统遥感影像采集方法无法做到的。本发明为压缩感知理论在遥感领域的实际应用奠定了基础。
申请公布号 CN103017738A 申请公布日期 2013.04.03
申请号 CN201210547673.X 申请日期 2012.12.18
申请人 程涛 发明人 程涛;朱国宾
分类号 G01C11/00(2006.01)I 主分类号 G01C11/00(2006.01)I
代理机构 代理人
主权项 1.基于二维压缩感知的遥感影像高效采集和增量更新方法,其特征是:所述方法的过程为:步骤一:以传统方法采集时相t1的遥感影像数据<img file="955096DEST_PATH_IMAGE001.GIF" wi="25" he="25" />(L表示<img file="372302DEST_PATH_IMAGE001.GIF" wi="25" he="25" />中列的总数);以压缩感知方法通过测量矩阵<img file="340258DEST_PATH_IMAGE002.GIF" wi="18" he="18" />(<img file="436390DEST_PATH_IMAGE003.GIF" wi="69" he="21" />,<img file="616836DEST_PATH_IMAGE004.GIF" wi="52" he="20" />)采集时相t1和t2的测量数据<img file="652663DEST_PATH_IMAGE005.GIF" wi="24" he="25" />和<img file="791520DEST_PATH_IMAGE006.GIF" wi="24" he="25" />(其中<img file="640527DEST_PATH_IMAGE007.GIF" wi="73" he="25" />,<img file="155822DEST_PATH_IMAGE008.GIF" wi="74" he="25" />);计算<img file="547620DEST_PATH_IMAGE009.GIF" wi="28" he="18" />(其中<img file="857379DEST_PATH_IMAGE010.GIF" wi="94" he="25" />);采用分段正交匹配追踪(Stagewise orthogonal matching pursuit,StOMP)算法(或其他重构算法,如MP、OMP、ROMP、StOMP、SP和CoSaMP)通过<img file="662524DEST_PATH_IMAGE011.GIF" wi="210" he="26" />(其中<img file="247089DEST_PATH_IMAGE012.GIF" wi="14" he="21" />表示<img file="86869DEST_PATH_IMAGE013.GIF" wi="28" he="18" />的列号,<img file="974054DEST_PATH_IMAGE014.GIF" wi="60" he="22" />)重构二维稀疏信号<img file="532074DEST_PATH_IMAGE013.GIF" wi="28" he="18" />(其中<img file="654751DEST_PATH_IMAGE015.GIF" wi="178" he="25" />),重构结果记为<img file="614617DEST_PATH_IMAGE016.GIF" wi="34" he="25" />;步骤二:设定临界值<img file="203861DEST_PATH_IMAGE017.GIF" wi="61" he="22" />,设定迭代误差<img file="983598DEST_PATH_IMAGE018.GIF" wi="66" he="22" />;判断<img file="283867DEST_PATH_IMAGE016.GIF" wi="34" he="25" />各列的稀疏度<img file="98239DEST_PATH_IMAGE019.GIF" wi="58" he="20" />,如果是,属于确定区域;否则,属于不确定区域;步骤三:计算成块分布的不确定区域数量<img file="451860DEST_PATH_IMAGE020.GIF" wi="26" he="25" />及各块不确定区域列的数量<img file="718894DEST_PATH_IMAGE021.GIF" wi="33" he="25" />;步骤四:判断<img file="324318DEST_PATH_IMAGE022.GIF" wi="52" he="25" />,如果是,设定迭代次数<img file="993197DEST_PATH_IMAGE024.GIF" wi="9" he="18" />的初始值为1;否则,执行步骤十一;步骤五:提取不确定区域<img file="517719DEST_PATH_IMAGE024.GIF" wi="9" he="18" />和不确定区域<img file="272049DEST_PATH_IMAGE024.GIF" wi="9" he="18" />的列数<img file="681165DEST_PATH_IMAGE021.GIF" wi="33" he="25" />;设定迭代次数<img file="204550DEST_PATH_IMAGE012.GIF" wi="14" he="21" />的初始值为1;步骤六:提取不确定区域<img file="899973DEST_PATH_IMAGE024.GIF" wi="9" he="18" />的左邻列<img file="141599DEST_PATH_IMAGE025.GIF" wi="49" he="26" />;根据左邻列<img file="682302DEST_PATH_IMAGE025.GIF" wi="49" he="26" />的支撑集<img file="466718DEST_PATH_IMAGE026.GIF" wi="48" he="26" />扩充<img file="67464DEST_PATH_IMAGE027.GIF" wi="32" he="20" />,得到不确定区域列的支撑集<img file="61964DEST_PATH_IMAGE028.GIF" wi="29" he="26" />;判断<img file="406358DEST_PATH_IMAGE029.GIF" wi="84" he="28" />,如果是,执行步骤七;否则,执行步骤八;步骤七:以支撑集<img file="543816DEST_PATH_IMAGE028.GIF" wi="29" he="26" />计算列<i>j</i>的稀疏解<img file="315463DEST_PATH_IMAGE030.GIF" wi="125" he="28" />;判断<img file="797260DEST_PATH_IMAGE031.GIF" wi="158" he="33" />和<img file="679765DEST_PATH_IMAGE032.GIF" wi="206" he="30" />,如果是,<img file="438774DEST_PATH_IMAGE033.GIF" wi="54" he="25" />,执行步骤八;否则使<img file="381322DEST_PATH_IMAGE034.GIF" wi="58" he="22" />,判断<img file="350415DEST_PATH_IMAGE035.GIF" wi="54" he="25" />,如果是,执行步骤十,否则返回执行步骤六;步骤八:提取不确定区域<img file="36611DEST_PATH_IMAGE024.GIF" wi="9" he="18" />的右邻列<img file="650126DEST_PATH_IMAGE036.GIF" wi="49" he="26" />;根据右邻列<img file="763576DEST_PATH_IMAGE036.GIF" wi="49" he="26" />的支撑集<img file="219965DEST_PATH_IMAGE037.GIF" wi="48" he="26" />扩充<img file="444273DEST_PATH_IMAGE027.GIF" wi="32" he="20" />,得到不确定区域列的新的支撑集<img file="177874DEST_PATH_IMAGE028.GIF" wi="29" he="26" />;判断<img file="196645DEST_PATH_IMAGE029.GIF" wi="84" he="28" />,如果是,执行步骤九;否则,执行步骤十;步骤九:以支撑集<img file="609172DEST_PATH_IMAGE028.GIF" wi="29" he="26" />计算列的稀疏解<img file="168329DEST_PATH_IMAGE030.GIF" wi="125" he="28" />;判断<img file="349912DEST_PATH_IMAGE031.GIF" wi="158" he="33" />和<img file="710224DEST_PATH_IMAGE032.GIF" wi="206" he="30" />,如果是,执行步骤十;否则使<img file="344468DEST_PATH_IMAGE038.GIF" wi="57" he="22" />,判断<img file="441737DEST_PATH_IMAGE039.GIF" wi="34" he="22" />,如果是,执行步骤十,否则返回执行步骤八;步骤十:使<img file="743405DEST_PATH_IMAGE040.GIF" wi="50" he="20" />,判断<img file="776083DEST_PATH_IMAGE041.GIF" wi="46" he="25" />,如果是,退出迭代1,否则返回执行步骤五;步骤十一:输出:重构信号<img file="897623DEST_PATH_IMAGE042.GIF" wi="36" he="25" />,并标出重构失败的列;<img file="798583DEST_PATH_IMAGE043.GIF" wi="101" he="25" />。
地址 545006 广西壮族自治区柳州市城中区东环大道268号广西工学院教工信箱30号