详情

原文传递一种用于矿物信息提取的高光谱影像处理方法

专利名称：	一种用于矿物信息提取的高光谱影像处理方法
摘要：	本发明涉及一种用于地开石、黄钾铁钒、利蛇纹石、菱锰矿、明矾石、石膏、水铝矿、阳起石、叶蜡石信息提取的高光谱影像处理方法，特别是涉及一种降低其他地物或噪声的影响，减少了处理数据量，用于地开石、黄钾铁钒、利蛇纹石、菱锰矿、明矾石、石膏、水铝矿、阳起石、叶蜡石信息提取的高光谱影像处理方法。本发明依次包括如下步骤：步骤1、获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；步骤2、对地面反射率的影像数据进行重采样，提取不同波段的影像；步骤3、比较不同波段影像的像元值，并标记比较结果；步骤4、根据步骤3确定地开石、黄钾铁钒、利蛇纹石、菱锰矿、明矾石、石膏、水铝矿、阳起石、叶蜡石的丰度值。
专利类型：	发明专利
申请人：	核工业北京地质研究院
发明人：	杨燕杰;邱骏挺;刘德长;赵英俊
专利状态：	有效
申请日期：	1900-01-20T20:00:00+0805
发布日期：	1900-01-20T15:00:00+0805
申请号：	CN201911327172.9
公开号：	CN111157459A
代理机构：	核工业专利中心
代理人：	闫兆梅
分类号：	G01N21/25;G06T5/00;G;G01;G06;G01N;G06T;G01N21;G06T5;G01N21/25;G06T5/00
申请人地址：	100029 北京市朝阳区小关东里十号院
主权项：	1.一种用于矿物信息提取的高光谱影像处理方法，其特征在于：依次包括如下步骤：步骤(1)、获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；步骤(2)、对地面反射率的影像数据进行重采样，提取不同波段的影像；步骤(3)、比较不同波段影像的像元值，并标记比较结果；步骤(4)、根据步骤(3)确定地开石、黄钾铁钒、利蛇纹石、菱锰矿、明矾石、石膏、水铝矿、阳起石、叶蜡石的丰度值。 2.根据权利要求1所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(2)中，提取影像的波段包括980nm，1040nm，1145nm，1175nm，1325nm，1340nm，1355nm，1415nm，1430nm，1445nm，1460nm，1490nm，1520nm，1535nm，1550nm，1595nm，1625nm，1700nm，1715nm，1745nm，1760nm，1775nm，1790nm，1805nm，1820nm，1940nm，1940nm，1955nm，1985nm，2015nm，2045nm，2060nm，2075nm，2090nm，2105nm，2120nm，2135nm，2150nm，2165nm，2180nm，2195nm，2210nm，2225nm，2230nm，2240nm，2255nm，2270nm，2285nm，2300nm，2315nm，2330nm，2345nm，2360nm，2375nm，2390nm，2405nm，2420nm。 3.根据权利要求1所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(3)包括如下步骤，当1145nm波段影像的每个像元值都大于1175nm波段影像相对应的像元值，判断满足条件一；获得1145nm波段影像减去1175nm波段影像的结果影像H1；当1325nm波段影像的每个像元值都大于1445nm波段影像相对应的像元值，判断满足条件二；获得1325nm波段影像减去1445nm波段影像的结果影像H2；当1340nm波段影像的每个像元值都大于1355nm波段影像相对应的像元值，判断满足条件三；获得1340nm波段影像减去1455nm波段影像的结果影像H3；当1355nm波段影像的每个像元值都大于980nm波段影像相对应的像元值，判断满足条件四；获得1355nm波段影像减去980nm波段影像的结果影像H4；当1355nm波段影像的每个像元值都大于1415nm波段影像相对应的像元值，判断满足条件五；获得1355nm波段影像减去1415nm波段影像的结果影像H5；当1355nm波段影像的每个像元值都大于1445nm波段影像相对应的像元值，判断满足条件六；获得1355nm波段影像减去1445nm波段影像的结果影像H5；当1415nm波段影像的每个像元值都大于1430nm波段影像相对应的像元值，判断满足条件七；获得1415nm波段影像减去1430nm波段影像的结果影像H7；当1415nm波段影像的每个像元值都大于1460nm波段影像相对应的像元值，判断满足条件八；获得1415nm波段影像减去1460nm波段影像的结果影像H8；当1445nm波段影像的每个像元值都大于1415nm波段影像相对应的像元值，判断满足条件九；获得1445nm波段影像减去1415nm波段影像的结果影像H9；当1460nm波段影像的每个像元值都大于1445nm波段影像相对应的像元值，判断满足条件十；获得1460nm波段影像减去1445nm波段影像的结果影像H10；当1490nm波段影像的每个像元值都大于1445nm波段影像相对应的像元值，判断满足条件十一；获得1490nm波段影像减去1445nm波段影像的结果影像H11；当1490nm波段影像的每个像元值都大于1520nm波段影像相对应的像元值，判断满足条件十二；获得1490nm波段影像减去1520nm波段影像的结果影像H12；当1520nm波段影像的每个像元值都大于1490nm波段影像相对应的像元值，判断满足条件十三；获得1520nm波段影像减去1490nm波段影像的结果影像H13；当1535nm波段影像的每个像元值都大于1445nm波段影像相对应的像元值，判断满足条件十四；获得1535nm波段影像减去1445nm波段影像的结果影像H14；当1595nm波段影像的每个像元值都大于1550nm波段影像相对应的像元值，判断满足条件十五；获得1595nm波段影像减去1550nm波段影像的结果影像H15；当1625nm波段影像的每个像元值都大于1550nm波段影像相对应的像元值，判断满足条件十六；获得1625nm波段影像减去1550nm波段影像的结果影像H16；当1700nm波段影像的每个像元值都大于1040nm波段影像相对应的像元值，判断满足条件十七；获得1700nm波段影像减去1040nm波段影像的结果影像H17；当1700nm波段影像的每个像元值都大于1745nm波段影像相对应的像元值，判断满足条件十八；获得1700nm波段影像减去1745nm波段影像的结果影像H18；当1715nm波段影像的每个像元值都大于1760nm波段影像相对应的像元值，判断满足条件十九；获得1715nm波段影像减去1760nm波段影像的结果影像H19；当1775nm波段影像的每个像元值都大于1940nm波段影像相对应的像元值，判断满足条件二十；获得1775nm波段影像减去1940nm波段影像的结果影像H20；当1790nm波段影像的每个像元值都大于2015nm波段影像相对应的像元值，判断满足条件二十一；获得1790nm波段影像减去2015nm波段影像的结果影像H21；当1805nm波段影像的每个像元值都大于1775nm波段影像相对应的像元值，判断满足条件二十二；获得1805nm波段影像减去1775nm波段影像的结果影像H22；当1820nm波段影像的每个像元值都大于1745nm波段影像相对应的像元值，判断满足条件二十三；获得1820nm波段影像减去1745nm波段影像的结果影像H23；当1820nm波段影像的每个像元值都大于1940nm波段影像相对应的像元值，判断满足条件二十四；获得1820nm波段影像减去1940nm波段影像的结果影像H24；当1820nm波段影像的每个像元值都大于1955nm波段影像相对应的像元值，判断满足条件二十五；获得1820nm波段影像减去1955nm波段影像的结果影像H25；当2045nm波段影像的每个像元值都大于1985nm波段影像相对应的像元值，判断满足条件二十六；获得2045nm波段影像减去1985nm波段影像的结果影像H26；当2045nm波段影像的每个像元值都大于2015nm波段影像相对应的像元值，判断满足条件二十七；获得2045nm波段影像减去2015nm波段影像的结果影像H27；当2045nm波段影像的每个像元值都大于2075nm波段影像相对应的像元值，判断满足条件二十八；获得2045nm波段影像减去2075nm波段影像的结果影像H28；当2060nm波段影像的每个像元值都大于1940nm波段影像相对应的像元值，判断满足条件二十九；获得2060nm波段影像减去1940nm波段影像的结果影像H29；当2060nm波段影像的每个像元值都大于2210nm波段影像相对应的像元值，判断满足条件三十；获得2060nm波段影像减去2210nm波段影像的结果影像H30；当2075nm波段影像的每个像元值都大于2120nm波段影像相对应的像元值，判断满足条件三十一；获得2075nm波段影像减去2120nm 波段影像的结果影像H31；当2105nm波段影像的每个像元值都大于2090nm波段影像相对应的像元值，判断满足条件三十二；获得2105nm波段影像减去2090nm波段影像的结果影像H32；当2105nm波段影像的每个像元值都大于2165nm波段影像相对应的像元值，判断满足条件三十三；获得2105nm波段影像减去2165nm波段影像的结果影像H33；当2120nm波段影像的每个像元值都大于2165nm波段影像相对应的像元值，判断满足条件三十四；获得2120nm波段影像减去2165nm波段影像的结果影像H34；当2135nm波段影像的每个像元值都大于2165nm波段影像相对应的像元值，判断满足条件三十五；获得2135nm波段影像减去2165nm波段影像的结果影像H35；当2135nm波段影像的每个像元值都大于2180nm波段影像相对应的像元值，判断满足条件三十六；获得2135nm波段影像减去2180nm波段影像的结果影像H36；当2150nm波段影像的每个像元值都大于2255nm波段影像相对应的像元值，判断满足条件三十七；获得2150nm波段影像减去2255nm波段影像的结果影像H37；当2195nm波段影像的每个像元值都大于2165nm波段影像相对应的像元值，判断满足条件三十八；获得2195nm波段影像减去2165nm波段影像的结果影像H38；当2195nm波段影像的每个像元值都大于2180nm波段影像相对应的像元值，判断满足条件三十九；获得2195nm波段影像减去2180nm波段影像的结果影像H39；当2195nm波段影像的每个像元值都大于2210nm波段影像相对应的像元值，判断满足条件四十；获得2195nm波段影像减去2210nm波段影像的结果影像H40；当2195nm波段影像的每个像元值都大于2225nm波段影像相对应的像元值，判断满足条件四十一；获得2195nm波段影像减去2225nm波段影像的结果影像H41；当2225nm波段影像的每个像元值都大于2120nm波段影像相对应的像元值，判断满足条件四十二；获得2225nm波段影像减去2120nm波段影像的结果影像H42；当2225nm波段影像的每个像元值都大于2255nm波段影像相对应的像元值，判断满足条件四十三；获得2225nm波段影像减去2255nm波段影像的结果影像H43；当2225nm波段影像的每个像元值都大于2300nm波段影像相对应的像元值，判断满足条件四十四；获得2225nm波段影像减去2300nm波段影像的结果影像H44；当2225nm波段影像的每个像元值都大于2360nm波段影像相对应的像元值，判断满足条件四十五；获得2225nm波段影像减去2360nm波段影像的结果影像H45；当2240nm波段影像的每个像元值都大于2210nm波段影像相对应的像元值，判断满足条件四十六；获得2240nm波段影像减去2210nm波段影像的结果影像H46；当2240nm波段影像的每个像元值都大于2255nm波段影像相对应的像元值，判断满足条件四十七；获得2240nm波段影像减去2255nm波段影像的结果影像H47；当2240nm波段影像的每个像元值都大于2315nm波段影像相对应的像元值，判断满足条件四十八；获得2240nm波段影像减去2315nm波段影像的结果影像H48；当2255nm波段影像的每个像元值都大于2210nm波段影像相对应的像元值，判断满足条件四十九；获得2255nm波段影像减去2210nm波段影像的结果影像H49；当2270nm波段影像的每个像元值都大于2300nm波段影像相对应的像元值，判断满足条件五十；获得2270nm波段影像减去2300nm波段影像的结果影像H50；当2285nm波段影像的每个像元值都大于2300nm波段影像相对应的像元值，判断满足条件五十一；获得2285nm波段影像减去2300nm波段影像的结果影像H51；当2285nm波段影像的每个像元值都大于2315nm波段影像相对应的像元值，判断满足条件五十二；获得2285nm波段影像减去2315nm波段影像的结果影像H52；当2330nm波段影像的每个像元值都大于2405nm波段影像相对应的像元值，判断满足条件五十三；获得2330nm波段影像减去2405nm 波段影像的结果影像H53；当2330nm波段影像的每个像元值都大于2255nm波段影像相对应的像元值，判断满足条件五十四；获得2330nm波段影像减去2255nm波段影像的结果影像H54；当2330nm波段影像的每个像元值都大于2315nm波段影像相对应的像元值，判断满足条件五十五；获得2330nm波段影像减去2315nm波段影像的结果影像H55；当2330nm波段影像的每个像元值都大于2390nm波段影像相对应的像元值，判断满足条件五十六；获得2330nm波段影像减去2390nm波段影像的结果影像H56；当2330nm波段影像的每个像元值都大于2405nm波段影像相对应的像元值，判断满足条件五十七；获得2330nm波段影像减去2405nm波段影像的结果影像H57；当2345nm波段影像的每个像元值都大于2315nm波段影像相对应的像元值，判断满足条件五十八；获得2345nm波段影像减去2315nm波段影像的结果影像H58；当2345nm波段影像的每个像元值都大于2330nm波段影像相对应的像元值，判断满足条件五十九；获得2345nm波段影像减去2330nm波段影像的结果影像H59；当2360nm波段影像的每个像元值都大于2375nm波段影像相对应的像元值，判断满足条件六十；获得2360nm波段影像减去2375nm波段影像的结果影像H60；当2360nm波段影像的每个像元值都大于2390nm波段影像相对应的像元值，判断满足条件六十一；获得2360nm波段影像减去2390nm波段影像的结果影像H61；当2360nm波段影像的每个像元值都大于2405nm波段影像相对应的像元值，判断满足条件六十二；获得2360nm波段影像减去2405nm波段影像的结果影像H62；当2375nm波段影像的每个像元值都大于2390nm波段影像相对应的像元值，判断满足条件六十三；获得2375nm波段影像减去2390nm波段影像的结果影像H63；当2405nm波段影像的每个像元值都大于2360nm波段影像相对应的像元值，判断满足条件六十四；获得2405nm波段影像减去2360nm波段影像的结果影像H64；当2405nm波段影像的每个像元值都大于2390nm波段影像相对应的像元值，判断满足条件六十五；获得2405nm波段影像减去2390nm波段影像的结果影像H65；当2420nm波段影像的每个像元值都大于2390nm波段影像相对应的像元值，判断满足条件六十六；获得2420nm波段影像减去2390nm波段影像的结果影像H66。 4.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件五、条件九、条件三十六、条件三十九、条件四十、条件四十六、条件四十七、条件五十一、条件五十六、条件六十五的像元，获得上述像元范围内的H5、H9、H36、H39、H40、H46、H47、H51、H56、H65相加的和H69，H69值越高，代表像元内的地开石丰度值越高。 5.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件四、条件八、条件十五、条件二十一、条件三十一、条件四十二、条件四十三、条件五十四、条件五十七的像元，获得上述像元范围内的H4、H8、H15、H21、H31、H42、H43、H54、H57相加的和H70，H70值越高，代表像元内的黄钾铁钒丰度值越高。 6.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件三、条件二十五、条件二十六、条件四十一、条件四十四、条件五十九、条件六十的像元，获得上述像元范围内的H3、H25、H26、H41、H44、H59、H60相加的和H71。H71值越高，代表像元内的利蛇纹石丰度值越高。 7.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件十四、条件二十四、条件二十七、条件三十五、条件四十五、条件六十四的像元，获得上述像元范围内的H14、H24、H27、H35、H45、H64相加的和H72，H72值越高，代表像元内的菱锰矿丰度值越高。 8.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件一、条件二、条件十、条件十八、条件二十三、条件二十四、条件二十九、条件三十、条件五十三的像元，获得上述像元范围内的H1、H2、H10、H18、H23、H24、H29、H30、H53相加的和H73，H73值越高，代表像元内的石膏丰度值越高。 9.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件六、条件十一、条件十二、条件十六、条件二十、条件三十七的像元，获得上述像元范围内的H6、H11、H12、H16、H20、H37相加的和H74，H74值越高，代表像元内的阳起石丰度值越高。 10.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件二十八、条件三十二、条件三十三、条件三十八、条件四十八、条件五十五、条件六十三的像元，获得上述像元范围内的H28、H32、H33、H38、H48、H55、H63相加的和H75，H75值越高，代表像元内的叶蜡石丰度值越高。 11.根据权利要求3所述的一种用于矿物信息提取的高光谱影像处理方法，其特征在于：所述步骤(4)中，选择同时满足条件七、条件十三、条件十九、条件二十二、条件二十四、条件三十四、条件三十九、条件四十九、条件五十二、条件六十二的像元，获得上述像元范围内的H7、H13、H19、H22、H24、H34、H39、H49、H52、H62相加的和H76。H76值越高，代表像元内的明矾石丰度值越高。
所属类别：	发明专利