一种用于水铝矿信息提取的高光谱影像处理方法技术

本发明专利技术属于高光谱影像处理方法，特别是涉及一种用于水铝矿信息提取的高光谱影像处理方法。它包括：步骤一，预处理；步骤二：采样，对波段在1355nm，1445nm，1490nm，1520nm，1550nm，1625nm，1775nm，1940nm，2150nm，2255nm的图像采样；步骤三：判断，和步骤四：计算。本发明专利技术本方法的效果在于：只用了10个波段，相对于高光谱影像SASI全波段101个波段，需要处理的数据量减少了90％，并且由于是计算机自动一步提取，减少了主成分变换、端元波谱的选择等操作步骤，运算速度可以提高了10倍以上。由于去除了大部分对信息提取关系不大的波段，减少其他物质或噪声对其光谱的干扰，提高了信息提取的精度。对高光谱影像数据中水铝矿信息的快速提取具有较好的作用和意义。

【技术实现步骤摘要】

本专利技术属于高光谱影像处理方法，特别是涉及。
技术介绍
当前的高光谱遥感影像的水铝矿信息提取方法主要是光谱全波段匹配或是部分连续波段的光谱匹配，具体算法有光谱角、混合解调滤波等，由于地表的物质组成很少是由单一矿物组成的，这些方法在信息提取的过程易受其他地物波谱或噪声的影响，提取信息精度相对较低。其次现有的光谱提取方法人工操作步骤多，增加了人为判断误差。第三是高光谱数据波段多，数据量大，现有的方法处理时间长，降低了数据处理的速度和应用规模。因此，如何在水铝矿信息提取的过程中减少其他地物或噪声的影响、人工操作步骤和处理数据量，成为当前高光谱遥感影像处理的前沿之一。
技术实现思路
本专利技术要解决的技术问题是提供一种降低其他地物或噪声的影响，减少了处理的数据量的用于水铝矿信息提取的高光谱影像处理方法。通过对高光谱影像重采样，提取特定波段，进行一系列判断和计算，计算出影像范围内的不同区域水铝矿的丰度值。为解决上述技术问题，本专利技术提供，依次包括:步骤一，预处理获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；步骤二:采样对步骤一得到的影像数据进行采样，对波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm, 1625nm, 1775nm, 1940nm, 2150nm, 2255nm 的图像采样，并依次记录为 bl ?blO,即 bl为波段1355nm的采样数据，b2为波段1445nm的采样数据，以此类推；步骤三:判断按进行下面一系列判断，并记录结果al = (bl 大于 b2)；a2 = (b2 小于 b3)；a3 = (b3 大于 b4)；a4 = (b5 小于 b6)；a5 = (b7 大于 b8)；a6 = (b9 大于 blO)；上述判断是针对每次判断图像的相应像元进行判断的；本步骤结束后得到al?a6，共6个矩阵；步骤四:计算用下述公式进行计算bllbll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8_bl0所述的*表示相乘；用下面公式计算aOaO = al*a2*a3*a4*a5*a6*bll上述所有计算均为相应像元计算。如上所述的，其中，所述的步骤三中的判断是针对每次判断图像的相应像元进行判断的。如上所述的，其中，所述的对图像的相应像元进行判断是按照下述过程进行的:取需要判断图像的对应像元，判断像元灰度是否符合判断标准，如果符合，判定该像元位置的结果为“ 1 ”，否则判定该像元位置的结果为“0”，然后判断下一个对应像元，直到完成对整个图像的判断。本专利技术本方法的效果在于:只用了 10个波段，相对于高光谱影像SASI全波段101个波段，需要处理的数据量减少了 90 %，并且由于是计算机自动一步提取，减少了主成分变换、端元波谱的选择等操作步骤，运算速度可以提高了 10倍以上。由于去除了大部分对信息提取关系不大的波段，减少其他物质或噪声对其光谱的干扰，提高了信息提取的精度。对高光谱影像数据中水铝矿信息的快速提取具有较好的作用和意义。【具体实施方式】本专利技术依次包括下述步骤:步骤一，预处理获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；本步骤所述的进行大气校正是本领域公知的处理方法。步骤二:采样对步骤一得到的影像数据进行采样，对波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm, 1625nm, 1775nm, 1940nm, 2150nm, 2255nm 的图像采样，并依次记录为 bl ?blO,即 bl为波段1355nm的采样数据，b2为波段1445nm的采样数据，以此类推。每次采样得到的都是一幅灰度图，图中每一个像元的值都是其灰度值，即bl为一幅灰度图，bl图像的(1，1)点的值为灰度值，其余点依此类推，其余采样图也依次类推。步骤三:判断按进行下面一系列判断，并记录结果al = (bl 大于 b2)；a2 = (b2 小于 b3)；a3 = ((b3 大于 b4)；a4 = (b5 小于 b6)；a5 = (b7 大于 b8)；a6 = (b9 大于 blO)；上述判断是针对每次判断图像的相应像元进行判断的，以al = (bl大于b2)为例，取bl图像的某像元(例如(1，1)点)的灰度值，与b2图像的相应像元(当bl图像取(1，1)点，则b2图像也必须取(1，1)点)的灰度值，然后根据判断规则“bl大于b2”判断，当判断结果为“是”时，记录判断结果为1，否则记录结果为0。因此当al = (bl大于b2)判断完毕时，得到的al是与bl矩阵相同大小的矩阵，其中每个点的值是根据判断规则得到的判断结果(即每个点的值是0或1)。其它判断也按照类似的规则进行。本步骤结束后得到al?a6，共6个矩阵。步骤四:计算用下述公式进行计算bllbll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8_bl0 所述的 * 表不相乘。用下面公式计算aOaO = al*a2*a3*a4*a5*a6*bll上述所有计算均为相应像元计算，即使用不同图像的相应像元计算。以bll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8-bl0 公式为例，当计算点(x，y)时，取 bl、b3、b6、b7、b9、b2、b4、b5、b8、bl0的点(x, y)的灰度值参与计算,得到的结果是bll的点(x, y)的值。又例如 aO = al*a2*a3*a4*a5*a6*bll,当计算点(x, y)时，取 al、a2、a3、a4、a5、a6、bll 的点(x，y)的值参与计算，得到的结果是aO的点(X, y)的值。计算得到的aO就是水铝矿信息的丰度图，即图像中某区域的数值越大表示该区域水铝矿的丰度越高。【主权项】1.，其特征在于:包括下述步骤 步骤一，预处理 获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据； 步骤二:采样 对步骤一得到的影像数据进行采样，对波段在1355nm，1445nm, 1490nm, 1520nm, 1550nm，1625nm, 1775nm, 1940nm, 2150nm, 2255nm的图像采样，并依次记录为bl?blO,即bl为波段2045nm的采样数据，b2为波段2075nm的采样数据，以此类推； 步骤三:判断 按进行下面一系列判断，并记录结果 al = (bl 大于 b2)； a2 = (b2 小于 b3)； a3 = ((b3 大于 b4)； a4 = (b5 小于 b6)； a5 = (b7 大于 b8)； a6 = (b9 大于 blO)； 上述判断是针对每次判断图像的相应像元进行判断的； 本步骤结束后得到al?a6，共6个矩阵； 步骤四:计算 用下述公式进行计算bllbll = bl+2*b3+b6+b7+b9-2*b2-b4-b5-b8-bl0 所述的*表示相乘； 用下面公式计算aOaO = al*a2*a3*a4*a5*a6*bll 上述所有计算均为相应像元计算。2.如权利要求1所述的，其特征在于:所述的步骤三中的判断是针对每次判断本文档来自技高网...

【技术保护点】
一种用于水铝矿信息提取的高光谱影像处理方法，其特征在于：包括下述步骤 步骤一，预处理 获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据； 步骤二：采样 对步骤一得到的影像数据进行采样，对波段在1355nm,1445nm,1490nm,1520nm,1550nm,1625nm,1775nm,1940nm,2150nm，2255nm的图像采样，并依次记录为b1～b10，即b1为波段2045nm的采样数据，b2为波段2075nm的采样数据，以此类推； 步骤三：判断 按进行下面一系列判断，并记录结果 a1＝(b1大于b2)； a2＝(b2小于b3)； a3＝((b3大于b4)； a4＝(b5小于b6)； a5＝(b7大于b8)； a6＝(b9大于b10)； 上述判断是针对每次判断图像的相应像元进行判断的； 本步骤结束后得到a1～a6，共6个矩阵； 步骤四：计算 用下述公式进行计算b11 b11＝b1+2*b3+b6+b7+b9‑2*b2‑b4‑b5‑b8‑b10 所述的*表示相乘； 用下面公式计算a0 a0＝a1*a2*a3*a4*a5*a6*b11 上述所有计算均为相应像元计算。...

【技术特征摘要】

【专利技术属性】
技术研发人员：杨燕杰，赵英俊，武鼎，宗自华，张川，陆冬华，钱坤，车永飞，叶发旺，李瀚波，高仉生，吴文欢，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京;11