发明名称 一种地表粗糙度参与阴影模型的遥感图像反立体校正方法
摘要 本发明提供一种地表粗糙度参与阴影模型的遥感图像反立体校正方法,该方法是对光学遥感红、绿、蓝(RGB)三色图像进行RGB色彩空间向HIS(饱和度、彩色亮度和色相)色彩空间变换,获得彩色亮度I(Intensity)、色相H(Hue)和饱和度S(Saturation)三分量,并以线性标准化地表粗糙度<img file="2014100392879100004DEST_PATH_IMAGE002.GIF" wi="22" he="20" />作为太阳西北方向地形阴影图<img file="2014100392879100004DEST_PATH_IMAGE004.GIF" wi="40" he="20" />(ShadedReliefModel)的权重,以<img file="2014100392879100004DEST_PATH_IMAGE006.GIF" wi="41" he="20" />作为彩色亮度<img file="2014100392879100004DEST_PATH_IMAGE008.GIF" wi="14" he="18" />权重,求和构成新的彩色亮度分量<img file="2014100392879100004DEST_PATH_IMAGE010.GIF" wi="17" he="18" />,最后基于新的彩色亮度数据<img file="767361DEST_PATH_IMAGE010.GIF" wi="17" he="18" />,由HIS色彩空间反变换至RGB色彩空间,实现对遥感图像的反立体现象的校正。同传统的基于地形阴影模型反立体校正方法相比较,该方法能根据地形复杂程度的不同,进行反立体校正,因而具有反立体校正后图像光谱信息损失少的优点。
申请公布号 CN103810706A 申请公布日期 2014.05.21
申请号 CN201410039287.9 申请日期 2014.01.27
申请人 鲁东大学 发明人 王涛;何福红;蒋卫国;顾丽娟;张振华
分类号 G06T7/00(2006.01)I;G06T15/80(2011.01)I 主分类号 G06T7/00(2006.01)I
代理机构 烟台双联专利事务所(普通合伙) 37225 代理人 矫智兰
主权项 1. 一种地表粗糙度参与阴影模型的遥感图像反立体校正方法,其特征在于,包括以下步骤:(1)对光学遥感图像红、绿、蓝三波段(RGB三波段)数据和DEM(数字高程模型)数据进行配准,统一坐标系统为高斯-克吕格平面直角坐标系,并作重采样和裁剪工作,使两者具有相同的空间分辨率和图幅范围;(2)基于光学遥感图像红、绿、蓝三波段数据,进行红、绿、蓝色彩空间向彩色亮度I(Intensity)、色相H(Hue)和饱和度S(Saturation)色彩空间变换,获得彩色亮度<img file="2014100392879100001DEST_PATH_IMAGE002.GIF" wi="14" he="18" />、色相<img file="DEST_PATH_IMAGE004.GIF" wi="20" he="18" />和饱和度<img file="DEST_PATH_IMAGE006.GIF" wi="16" he="20" />三分量数据;(3)根据遥感图像成像时刻(年月日时分秒)及图像中心点经纬度坐标,计算成像时刻太阳方位角<img file="DEST_PATH_IMAGE008.GIF" wi="14" he="22" />和太阳高度角<img file="DEST_PATH_IMAGE010.GIF" wi="17" he="16" />;(4)基于DEM数据,计算地表粗糙度,获得地表粗糙度数据<img file="DEST_PATH_IMAGE012.GIF" wi="18" he="16" />;(5)对地表粗糙度数据<img file="664959DEST_PATH_IMAGE012.GIF" wi="18" he="16" />做线性标准化处理,获得线性标准化地表粗糙度数据<img file="DEST_PATH_IMAGE014.GIF" wi="22" he="20" />;(6)基于研究区DEM数据,以遥感图像成像时刻太阳方位角<img file="582100DEST_PATH_IMAGE008.GIF" wi="14" he="22" />加<img file="DEST_PATH_IMAGE016.GIF" wi="16" he="16" />(180°)和遥感图像成像时刻太阳高度角<img file="361837DEST_PATH_IMAGE010.GIF" wi="17" he="16" />为参数计算研究区正立体地形阴影数据,即太阳西北方向地形阴影图<img file="DEST_PATH_IMAGE018.GIF" wi="40" he="20" />(Shaded Relief Model);(7)太阳西北方向地形阴影图<img file="570095DEST_PATH_IMAGE018.GIF" wi="40" he="20" />和彩色亮度<img file="384468DEST_PATH_IMAGE002.GIF" wi="14" he="18" />分量数据分别以<img file="269247DEST_PATH_IMAGE014.GIF" wi="22" he="20" />和1-<img file="536280DEST_PATH_IMAGE014.GIF" wi="22" he="20" />作为权重,计算获得新的彩色亮度<img file="DEST_PATH_IMAGE020.GIF" wi="17" he="18" />分量数据(式1):<img file="2014100392879100001DEST_PATH_IMAGE022.GIF" wi="176" he="22" />(1);(8)基于HIS-&gt;RGB色彩空间变换算法,将(2)获得的色相<img file="545300DEST_PATH_IMAGE004.GIF" wi="20" he="18" />和饱和度<img file="214179DEST_PATH_IMAGE006.GIF" wi="16" he="20" />以及(7)获得的新的彩色亮度<img file="207543DEST_PATH_IMAGE020.GIF" wi="17" he="18" />三分量数据反变换至RGB色彩空间,此时,遥感图像已经正立体化,实现对原始遥感RGB图像的反立体现象的校正。
地址 264025 山东省烟台市芝罘区经旗中路186号