利用阴影指数和方位搜索法检测Landsat 8 OLI影像中云影

利用阴影指数和方位搜索法检测Landsat 8 OLI影像中云影 Detection of cloud shadow in Landsat 8 OLI image by shadow index and azimuth search method 王凌 WANG Ling 1 first-author 赵庚星 ZHAO Genxing 1 姜远茂 JIANG Yuanmao 2 朱西存 ZHU Xicun 1 常春艳 CHANG Chunyan 1 王明媚 WANG Mingmei 3 山东农业大学资源与环境学院土肥资源高效利用国家工程实验室 山东农业大学资源与环境学院土肥资源高效利用国家工程实验室 College of Resources and Environment,Shandong Agricultural University/National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources College of Resources and Environment,Shandong Agricultural University/National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources 山东农业大学园艺科学与工程学院 山东农业大学园艺科学与工程学院 College of Horticulture Science and Engineering,Shandong Agricultural University College of Horticulture Science and Engineering,Shandong Agricultural University 山东省乳山市国土资源局 山东省乳山市国土资源局 Land and Resources Bureau of Rushan City Land and Resources Bureau of Rushan City Landsat 8 OLI影像已成为重要的数据源,但受云及云影的影响较大,降低了数据的可用性,因此,快速识别云及云影,为后续的数据恢复有着积极的应用价值。通过构建云指数（CI）、归一化暗像元指数（NDPI）和比值阴影指数（RSI）,采用阈值法和方位角搜索法,以两景Landsat 8 OLI影像（一景试验影像,另一景验证影像）为例进行云及云影检测。每类随机选取200个样点进行精度分析,结果表明:CI可快速区分OLI影像中的云区与非云区,厚云样本点正确识别率达到99%;NDPI与归一化植被指数（NDVI）构建的比值阴影指数RSI放大了水体、云影与其他阴影间的差异,更便于区分;方位搜索合理设置搜索方位角和搜索距离,简化了云影与云的相对关系模型,可准确区分水体与云影,两者的正确识别率都超过93%,弥补了阈值法的局限性。本方法可行快捷,为OLI影像的后续应用提供了基础,可有效提高其利用精度。 Landsat 8 OLI images have become important data sources;however,they are usually covered by clouds and cloud shadows,which reduce data availability.Therefore,a rapid method for detecting clouds and cloud shadows in a single image is necessary for followup data recovery applications.Threshold setting has been a commonly used method;however,it is difficult to use because the same threshold value usually indicates different objectives and varies across different data sources,as well as across different period images.Thus,the relative position relationship between a cloud shadow and a cloud is the key to detect a cloud shadow.First,a cloud index（CI） with a cirrus band and coastal/aerosol band was established to distinguish the thin and thick cloud pixels by setting thresholds.Second,a normalized dark pixel index（NDPI） with coastal/aerosol band and short-wave infrared 2 band was established.Furthermore,a ratio of shadow index（RSI） was established based on the NDPI and normalized differential vegetation index.RSI was used to distinguish the potential cloud shadow pixels in thin-cloud and non-cloud areas also by using thresholds.Third,an azimuth search method based on the solar azimuth angle and an appropriate searching distance were adopted to detect the cloud shadow from the potential cloud shadow pixels in the images.Two Landsat 8 OLI images were selected for the case study;the image captured on April 22,2014 was used for the test,whereas the image captured on July 12,2015 was used for validation.For each type（e.g.,thick cloud,water,cloud shadow,or other kinds of shadow）,200 random sampling points were used to assess the detection accuracy.Results showed that CI could quickly distinguish cloud from on-cloud pixels by statistics according to the cloud coverage ratio in the header files,with thresholds of CI ≥ 0.0011 and CI ≥ 0.0048 in the test and validation OLI images for the thick cloud type.Both user accuracy rates of detection for the thick cloud samples were over 99%.RSI could enhance the difference among water,cloud shadow,and other kinds of shadow,thereby facilitating the differentiation among the different types.The pixels with 0.45 ≤ RSI < 0.76 in the test image and 0.36 ≤ RSI < 0.76 in the validation image belong to potential cloud shadows.Using the solar azimuth angle at the time of imaging as the searching azimuth angle and a reasonable searching distance（500-2200 m in the test image and270-800 m in the validation image）,the azimuth searching method simplified the models of the relative position relationship of a cloud shadow to a cloud and accurately distinguished cloud shadow from water and other shadows.The detection precision levels for both test and validation reached more than 93%,which compensated for the limitation of the threshold method.The proposed shadow detection method combines threshold with simplified relative relationship model and leverages band difference.The method is feasible and rapid when applied to a single image,further improving the utilization accuracy of OLI images. 比值阴影指数 云指数 方位搜索 云影检测 OLI影像 ratio of shadow index cloud index azimuth search cloud shadow detection OLI image TP751 国家青年科学基金（编号:41301482）;国家自然科学基金（编号:41271369）;山东省自主创新及成果转化专项（编号:2014ZZCX07402）;山东农业大学博士后基金（编号:76390）;2013年度山东省高等学校青年骨干教师国内访问学者项目~~ 2021-06-16 6 20