一种多光谱影像标定方法技术

本发明专利技术提供一种多光谱影像标定方法，包括步骤：S1，测量三种校准目标的第一反射率；S2，获得多光谱影像；S3，利用所述多光谱影像获得对应三种校准目标和对应待标定地物目标的像元的亮度值；S4，根据所述多光谱影像将所述第一反射率重采样以得到第二反射率；S5，根据所述第二反射率和对应三种校准目标的像元的亮度值构建标定方程；S6，利用所述标定方程和对应待标定地物目标的像元的亮度值获得待标定目标的标定反射率。本发明专利技术的多光谱影像标定方法，通过采用价格便宜，操作方便的黑布、灰布和白布作为校准目标，利用分段线性方程进行标定，标定过程简单易行，且标定结果的精度较高。

A method of multi spectral image calibration

The invention provides a multi spectral image calibration method, comprising the steps of: S1, the first three reflectivity measurement calibration target; S2, for multi spectral image; S3, using the multi spectral image to obtain brightness corresponding to the three kinds of calibration target and corresponding calibration target pixel value; according to the S4. The first reflectance spectral image resampling to get second S5, reflectance; calibration equation was constructed according to the luminance of the second reflectivity and three corresponding calibration target pixel value; S6 calibration equation and the corresponding to the target pixel brightness calibration object is obtained to be calibrated using the calibration target reflectivity. Calibration method for multi spectral images of the invention, through the use of cheap and convenient operation of the cloth, grey cloth and cloth as the calibration target, using piecewise linear calibration equation, the calibration process is simple, high precision and calibration results.

【技术实现步骤摘要】

本专利技术涉及无人机遥感技术，尤其涉及多光谱影像标定方法。
技术介绍
无人机遥感技术能够快速获取国土资源、自然环境等空间的遥感信息，现已获得越来越广泛的应用。近年来，随着材料科学的飞速发展和野外环境研究的需要，多光谱数据，尤其是对于反射率的标定成为研究的重要一步。无人机遥感的飞行高度低，受大气传播干扰极小，可以忽略大气影响。通过利用无人机遥感技术获得多光谱数据，进而对地物的反射率进行标定，现已提出了多种方法。现有技术中常用的利用线性方程进行标定的方法，，其缺点是校准精度不好，另外，校准目标一般选取漫反射标准灰度参照板，其价格昂贵，并且由于其展开范围大，导致操作不方便。
技术实现思路
本专利技术要解决的技术问题是，提供一种校准精度好，操作方便的多光谱影像标定方法。为了解决上述问题，本专利技术提供了一种多光谱影像标定方法，包括以下步骤：S1，测量三种校准目标的第一反射率；S2，获得多光谱影像；S3，利用所述多光谱影像获得对应三种校准目标和对应待标定地物目标的像元的亮度值；S4，根据所述多光谱影像将所述第一反射率重采样以得到第二反射率；S5，根据所述第二反射率和对应三种校准目标的像元的亮度值构建标定方程；S6，利用所述标定方程和对应待标定地物目标的像元的亮度值获得待标定目标的标定反射率。优选地，三种校准目标分别是黑布、白布和灰布。优选地，S5具体为，根据黑布、白布和灰布的第二反射率以及分别对应黑布、白布和黑布的像元的亮度值计算黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移，并利用黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移构建分段式线性方程。优选地，所述分段式线性方程具体为其中，R*为标定反射率，a1、b1分别是利用黑布和灰布计算出的增益和偏移，a2、b2分别是利用灰布和白布计算出的增益和偏移，X为对应待标定地物目标的像元的亮度值，X灰为灰布的亮度值。优选地，S1中利用手持地物光谱仪测量第一反射率，S2中利用多光谱相机获得多光谱影像。优选地，所述第一反射率和所述第二反射率的采样波长位于400nm至1000nm的范围内。优选地，所述第一反射率的采样波长为400nm至1000nm，采样间隔为1nm，所述第二反射率的采样波段包括6个波段，其中心波长分别是490nm、550nm、680nm、720nm、800nm和900nm，其波长宽度为10nm；S4具体为，将1nm的采样间隔的第一反射率重采样到对应于6个波段的采样间隔的第二反射率。与现有技术相比，本专利技术的有益效果在于：通过采用价格便宜，操作方便的黑布、灰布和白布作为校准目标，利用分段线性方程进行标定，标定过程简单易行，且标定结果的精度较高。附图说明图1是本专利技术实施例的多光谱影像标定方法的流程图；图2是手持地物光谱仪测量黑布、灰布和白布三种校准目标的反射率的曲线图。具体实施方式下面结合附图和具体实施例对本专利技术作进一步详细描述，但不作为对本专利技术的限定。图1是本专利技术实施例的多光谱影像标定方法的流程图。如图1所示，本专利技术的多光谱影像标定方法，包括以下步骤：S1，测量三种校准目标的第一反射率；S2，获得多光谱影像；S3，利用所述多光谱影像获得对应三种校准目标和对应待标定地物目标的像元的亮度值；S4，根据所述多光谱影像将所述第一反射率重采样以得到第二反射率；S5，根据所述第二反射率和对应三种校准目标的像元的亮度值构建标定方程；S6，利用所述标定方程和对应待标定地物目标的像元的亮度值获得待标定目标的标定反射率。其中，步骤S3与步骤S4的顺序可以互换。通过上述方法，利用有限的校准目标便可以得到多光谱影像的大量像元所对应的待标定地物目标的标定反射率。其中，在S1步骤中，可以利用手持地物光谱仪测量第一反射率，步骤S2中可以利用多光谱相机获得多光谱影像。本实施例中，选取的三种校准目标分别是黑布、白布和灰布，黑布、白布和灰布可以是无纺布材料，其价格便宜，携带方便。本实施例中，步骤S5可以具体为，根据黑布、白布和灰布的第二反射率以及分别对应黑布、白布和黑布的像元的亮度值计算黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移，并利用黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移构建分段式线性方程。具体的，分段式线性方程可以为其中，R*为标定反射率，a1、b1分别是利用黑布和灰布计算出的增益和偏移，a2、b2分别是利用灰布和白布计算出的增益和偏移，X为对应待标定地物目标的像元的亮度值，X灰为灰布的亮度值。根据上述分段式线性方程，在步骤S6中，将对应待标定地物目标的像元的亮度值与作为校准目标的灰布的亮度值进行比较，并将对应待标定地物目标的像元的亮度值带入对应的方程中以计算该像元对应待标定地物目标的标定反射率。本实施例中，为了获得较好的标定结果，第一反射率和第二反射率的采样波长优选位于400nm至1000nm的范围内。例如，第一反射率的采样波长为400nm至1000nm，采样间隔为1nm。第二反射率的采样波段包括6个波段，其中心波长分别是490nm、550nm、680nm、720nm、800nm和900nm，其波长宽度为10nm。当然，第二反射率的采样波段、中心波长和波长宽度可以根据具体的情况而变化，一般的，中心波长选择颜色为蓝、绿、红、近红外所对应的波长范围内。步骤S4具体为，将1nm的采样间隔的第一反射率应用积分法重采样到对应于6个波段的采样间隔的第二反射率。下面以具体数据进行说明。利用手持地物光谱仪测量黑布、灰布和白布三种校准目标的反射率数据，具体如图2。其中，手持地物光谱仪记录的是三种校准目标在波长为325nm-1075nm的反射率，其包括了第一反射率的优选的波长范围400nm-1000nm，采样间隔为1nm。多光谱相机的6个波段的中心波长分别是，490nm、550nm、680nm、720nm、800nm和900nm，波长宽度为10nm。将手持地物光谱仪的1nm采样间隔的光谱数据重采样到多光谱相机对应于波段的采样间隔中，即对第一反射率进行重采样到第二反射率，得到的重采样结果如表1所示。表1三种校准目标的第二反射率中心波长490nm550nm680nm720nm800nm900nm黑布(R黑)0.0410380.0420810.0433990.0443610.0458270.046427灰布(R灰)0.2317160.2114110.2148320.2066690.3058340.320635白布(R白)0.9510610.9269060.9037140.9093310.8992590.900016在多光谱影像上得到对应黑布、灰布、白布的像元的亮度值，如表2所示。表2多光谱影像中的对应黑布、灰布、白布的像元的亮度值中心波长490nm550nm680nm720nm800nm900nm黑布(DN黑)29.5714317.2857118.142866.2857142.52.5灰布(DN灰)121.8571117.142997.2857178.1428681.4285765.85714白布(DN白)225.4667244239.6667226.8667197154.5333根据表1和表2的数据，通过公式可计算黑布与灰布之间的增益a1和偏移b1以及灰布与白布之间本文档来自技高网...

【技术保护点】
一种多光谱影像标定方法，其特征在于，包括以下步骤：S1，测量三种校准目标的第一反射率；S2，获得多光谱影像；S3，利用所述多光谱影像获得对应三种校准目标和对应待标定地物目标的像元的亮度值；S4，根据所述多光谱影像将所述第一反射率重采样以得到第二反射率；S5，根据所述第二反射率和对应三种校准目标的像元的亮度值构建标定方程；S6，利用所述标定方程和对应待标定地物目标的像元的亮度值获得待标定目标的标定反射率。

【技术特征摘要】
1.一种多光谱影像标定方法，其特征在于，包括以下步骤：S1，测量三种校准目标的第一反射率；S2，获得多光谱影像；S3，利用所述多光谱影像获得对应三种校准目标和对应待标定地物目标的像元的亮度值；S4，根据所述多光谱影像将所述第一反射率重采样以得到第二反射率；S5，根据所述第二反射率和对应三种校准目标的像元的亮度值构建标定方程；S6，利用所述标定方程和对应待标定地物目标的像元的亮度值获得待标定目标的标定反射率。2.根据权利要求1所述的多光谱影像标定方法，其特征在于，三种校准目标分别是黑布、白布和灰布。3.根据权利要求2所述的多光谱影像标定方法，其特征在于，S5具体为，根据黑布、白布和灰布的第二反射率以及分别对应黑布、白布和黑布的像元的亮度值计算黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移，并利用黑布与灰布之间的增益和偏移以及灰布与白布之间的增益和偏移构建分段式线性方程。4.根据权利要求3所述的多光谱影像标定方法，其特征在于，...

【专利技术属性】
技术研发人员：邓磊，张爱武，郝向磊，毛智慧，
申请(专利权)人：首都师范大学，
类型：发明
国别省市：北京;11