[发明专利]一种用于石膏信息提取的高光谱影像处理方法

摘要

本发明属于高光谱影像处理方法，特别是涉及一种用于石膏信息提取的高光谱影像处理方法。它包括步骤一，预处理；步骤二采样，对波段在1145nm，1175nm，1325nm，1445nm，1460nm，1700nm，1745nm，1820nm，1940nm，2060nm,2210nm，2330nm，2405nm的图像采样；步骤三判断，和步骤四计算。本发明本方法的效果在于只用了13个波段，相对于高光谱影像SASI全波段101个波段，需要处理的数据量减少了87%，并且由于是计算机自动一步提取，减少了主成分变换、端元波谱的选择等操作步骤，运算速度可以提高了7倍以上。由于去除了大部分对信息提取关系不大的波段，减少其他物质或噪声对其光谱的干扰，提高了信息提取的精度。对高光谱影像数据中石膏信息的快速提取具有较好的作用和意义。

主权项

一种用于石膏信息提取的高光谱影像处理方法，其特征在于：包括下述步骤步骤一，预处理获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；步骤二：采样对步骤一得到的影像数据进行采样，对波段在1145nm，1175nm，1325nm，1445nm，1460nm，1700nm，1745nm，1820nm，1940nm，2060nm,2210nm，2330nm，2405nm的图像采样，并依次记录为b1～b13，每次采样得到的都是一幅灰度图，图中每一个像元的值都是其灰度值，即b1为波段1145nm的采样数据，b2为波段1175nm的采样数据，以此类推；步骤三：判断进行下面一系列判断，并记录结果a1＝(b1大于b2)；a2＝(b3大于b4)；a3＝(b4小于b5)；a4＝(b6大于b7)；a5＝(b7小于b8)；a6＝(b8大于b9)；a7＝(b9小于b10)；a8＝(b10大于b11)；a9＝(b12大于b13)；上述判断是针对每次判断图像的相应像元进行判断的，当判断结果为“是”时，记录判断结果为1，否则记录结果为0，本步骤结束后得到a1～a9，共9个矩阵；步骤四：计算用下述公式进行计算b14b14＝b1+b3+b5+b6+2*b8+2*b10+b12‑b2‑2*b4‑2*b7‑2*b9‑b11‑b13所述的*表示相乘；用下面公式计算a0a0＝a1*a2*a3*a4*a5*a6*a7*a8*a9*b14a0为石膏信息的丰度图；上述所有计算均为相应像元计算，即使用不同图像的相应像元计算。