CN104361602A - 一种基于modis图像的水体颜色检测方法及装置 - Google Patents
一种基于modis图像的水体颜色检测方法及装置 Download PDFInfo
- Publication number
- CN104361602A CN104361602A CN201410693596.8A CN201410693596A CN104361602A CN 104361602 A CN104361602 A CN 104361602A CN 201410693596 A CN201410693596 A CN 201410693596A CN 104361602 A CN104361602 A CN 104361602A
- Authority
- CN
- China
- Prior art keywords
- coordinates system
- coordinate
- angle
- water body
- chromaticity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/001—Industrial image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明实施例提供一种基于MODIS图像的水体颜色检测方法及装置,其中方法包括:获取反映水体颜色的具有RGB波段组合的MODIS图像;将MODIS图像的每个象元的RGB值转换成三刺激值;基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;计算每个象元的色度坐标在预定坐标系中的新坐标;计算各象元的新坐标相对于预定坐标系的X轴的夹角;根据预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别;根据所确定的颜色主波长与水色级别,确定水体颜色。本发明实施例实现了采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的适用。
Description
技术领域
本发明涉及遥感图像颜色参量提取技术领域,更具体地说,涉及一种基于MODIS(Moderate-resolution Imaging Spectroradiometer,中分辨率成像光谱仪)图像的水体颜色检测方法及装置。
背景技术
水体颜色是太阳光线与水体中的光活性物质在可见光范围内作用的直接结果,与水体水质有着密切的关系。水体颜色是水质的常规监测指标之一,因此对水体颜色的进行检测,对于水质监测而言具有重要意义。
本发明的发明人研究发现:目前利用遥感手段对水体颜色进行宏观大范围长时间的遥感检测成为了一种需求;而水体中三种主要的光活性物质:浮游植物、悬浮物和有色可溶性有机物共同决定了水体的颜色;在水体中,三种物质表现出各种可能的组合,三种物质浓度独立变化且变化范围很大,水体光学特性复杂多变,使水体水色遥感成为一个多变量非线性的复杂问题;因此现有技术还没有一种能够适用于采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的方案。
可见,提供一种水体颜色检测方法,以适用采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测,成为本领域技术人员亟需解决的技术问题。
发明内容
有鉴于此,本发明实施例提供一种基于MODIS图像的水体颜色检测方法及装置,以实现适用采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的目的。
为实现上述目的,本发明实施例提供如下技术方案:
一种基于MODIS图像的水体颜色检测方法,包括:
获取反映水体颜色的具有RGB波段组合的MODIS图像;
将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,所述预定坐标系的Y轴与色度图坐标系的X轴平行,正向方向一致;
计算各象元的新坐标相对于预定坐标系的X轴的夹角;
根据预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别;
根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
其中,所述预定的夹角与颜色主波长的对应关系的确定过程包括:
选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2;
确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长点在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标;
通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标;
计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角;
确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系。
其中,所述预定的夹角与水色级别的对应关系的确定过程包括:
确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标;
通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标;
计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角;
确定各水色级别对应的夹角范围,该夹角范围即该水色级别与上下两个水色级别夹角角分线确定的夹角范围,落入这个角度范围的夹角值对应该水色级别,建立各水色级别与对应的夹角的对应关系。
其中,所述获取反映水体颜色的具有RGB波段组合的MODIS图像包括:
获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
其中,所述将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值包括:
基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
其中,所述预定坐标系的构建过程包括:
在CIE标准色度图中,基于色度坐标为(1/3,1/3)的等能白点构建二维左手坐标系,所构建的二维左手坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致,所构建的二维左手坐标系作为所述预定坐标系。
本发明实施例还提供一种基于MODIS图像的水体颜色检测装置,包括:
图像获取模块,用于获取反映水体颜色的具有RGB波段组合的MODIS图像;
三刺激值转换模块,用于将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
色度坐标计算模块,用于基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
新坐标计算模块,用于通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,所述预定坐标系的Y轴与色度图坐标系的X轴平行,正向方向一致;
夹角计算模块,用于计算各象元的新坐标相对于预定坐标系的X轴的夹角;
颜色波长级别确定模块,用于根据预定的夹角与颜色主波长的对应关系,确定所述夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定所述夹角的水色级别;
水体颜色确定模块,用于根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
其中,所述装置还包括:
第一对应关系确定模块,用于选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2,确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标,通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标,计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角,确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系;
第二对应关系确定模块,用于确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标,通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标,计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角,确定各水色级别对应的夹角,建立各水色级别与对应的夹角的对应关系。
其中,所述图像获取模块包括:
图像数据获取单元,用于获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
其中,所述三刺激值转换模块包括:
转换执行单元,用于基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
基于上述技术方案,本发明实施例提供的基于MODIS图像的水体颜色检测方法,通过获取反映水体颜色的具有RGB波段组合的MODIS图像,将MODIS图像的每个象元的RGB值转换成色度学中的三刺激值,计算各象元的三刺激值在色度图坐标系中的色度坐标,并且通过建立X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致的预定坐标系,从而通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,计算各象元的新坐标相对于预定坐标系的X轴的夹角;进而通过预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别,最终根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。本发明实施例提供的基于MODIS图像的水体颜色检测方法,利用遥感数据对水体水质进行宏观监测,发挥了卫星遥感环境监测的优势;同时基于MODIS反射率遥感图像,建立了水体颜色参量提取的方法,使遥感图像上的水体颜色能够定量化表示,实现了采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的适用。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的基于MODIS图像的水体颜色检测方法的流程图;
图2为本发明实施例提供的确定夹角与颜色主波长的对应关系的方法流程图;
图3为本发明实施例提供的确定夹角与水色级别的对应关系的方法流程图;
图4为本发明实施例提供的基于MODIS图像的水体颜色检测方法的另一流程图;
图5为本发明实施例提供的基于MODIS图像的水体颜色检测装置的结构框图;
图6为本发明实施例提供的基于MODIS图像的水体颜色检测装置的另一结构框图;
图7为本发明实施例提供的图像获取模块的结构框图;
图8为本发明实施例提供的三刺激值转换模块的结构框图;
图9为本发明实施例提供的基于MODIS图像的水体颜色检测装置的再一结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的基于MODIS图像的水体颜色检测方法的流程图,该方法可应用于具有数据处理能力的电子设备,参照图1,该方法可以包括:
步骤S100、获取反映水体颜色的具有RGB波段组合的MODIS图像;
可选的,MODIS图像可以认为是MODIS拍摄的反映水体颜色的图像数据;
可选的,RGB波段组合主要是指红光(R)范围波段,绿光(G)范围波段和蓝光(B)范围波段的组合;本发明实施例可通过MODIS获取地表反射率水体区域图像,并且该地表反射率水体区域图像中可包含RGB波段组合。
步骤S110、将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
可选的,本发明实施例可基于CIE(International Commission onIllumination,国际照明委员会)标准色度学将MODIS图像中每个象元的RGB值(即红光范围波段值,绿光范围波段值和蓝光范围波段值)转换成色度学中的三刺激值X,Y,Z;
步骤S120、基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
可选的,色度图坐标系可以是CIE标准色度学对应的色度坐标系;通过三刺激值在色度图坐标系中的坐标转换算法,可计算各象元的三刺激值所对应在色度图坐标系中的色度坐标。
步骤S130、通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标;
其中,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致;
可选的,本发明实施例可通过两个坐标系之间的坐标转换算法,计算每个象元在色度图坐标系中的色度坐标,在预定坐标系中的新坐标。
步骤S140、计算各象元的新坐标相对于预定坐标系的X轴的夹角;
若象元在新坐标系中的坐标为(x,y),则本发明实施例可计算(x,y)与预定坐标系的X轴的夹角α值。
步骤S150、根据预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别;
可选的,本发明实施例可计算色度图的颜色主波长在预定坐标系中对应的坐标,确定出色度图的颜色主波长在预定坐标系中对应的坐标,与预定坐标系的X轴的夹角,从而得到夹角与颜色主波长的对应关系,进而确定出步骤S140所计算的各夹角的颜色主波长;由于各颜色主波长均可对应相应的颜色,因此通过确定出的各夹角的颜色主波长,本发明实施例可确定出MODIS图像所反映的水体颜色分布情况。
可选的,本发明实施例可通过水体Forel-Ule比色表中各水色级别在预定坐标系中对应的坐标,确定出各水色级别在预定坐标系中对应的坐标,与预定坐标系的X轴的夹角,从而得到夹角与水色级别的对应关系,进而确定出步骤S140所计算的各夹角的水色级别;通过确定出的各夹角的水色级别,本发明实施例可确定出MODIS图像所反映的水体的水色级别分布情况。
步骤S160、根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
可选的,本发明实施例可仅以所确定的颜色主波长,确定MODIS图像所反映的水体颜色;也可仅以所确定的水色级别确定MODIS图像所反映的水体颜色;显然,本发明实施例也可采用颜色主波长与水色级别相结合的方式确定MODIS图像所反映的水体颜色。
本发明实施例提供的基于MODIS图像的水体颜色检测方法,通过获取反映水体颜色的具有RGB波段组合的MODIS图像,将MODIS图像的每个象元的RGB值转换成色度学中的三刺激值,计算各象元的三刺激值在色度图坐标系中的色度坐标,并且通过建立X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致的预定坐标系,从而通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,计算各象元的新坐标相对于预定坐标系的X轴的夹角;进而通过预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别,最终根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。本发明实施例提供的基于MODIS图像的水体颜色检测方法,利用遥感数据对水体水质进行宏观监测,发挥了卫星遥感环境监测的优势;同时基于MODIS反射率遥感图像,建立了水体颜色参量提取的方法,使遥感图像上的水体颜色能够定量化表示,实现了采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的适用。
可选的,图2示出了确定夹角与颜色主波长的对应关系的方法流程图,值得注意的是,在前文所述的夹角与颜色主波长的对应关系的表达含义的基础上,图2所示方法仅为实现夹角与颜色主波长的对应关系的表达含义的一种可选方式;参照图2,该方法可以包括:
步骤S200、选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2;
可选的,设定纳米波长范围可以是380nm(纳米)-700nm,其中包含321个整数纳米波长,则可对应321个颜色主波长。
步骤S210、确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长点在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标;
可选的,本发明实施例可计算380nm-700nm中的321个整数纳米波长,在色度图坐标系中的色度坐标,其中,一个整数纳米波长对应一个颜色主波长,得到321个颜色主波长所对应的321个色度坐标。
色度图坐标系可以为图1所示的与CIE标准色度学对应的色度坐标系。
步骤S220、通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标;
可选的,预定坐标系可以为图1所示的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致的坐标系。
步骤S230、计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角;
步骤S240、确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系。
可选的,本发明实施例可建立321个颜色主波长与对应的夹角的对应关系,根据该对应关系可构建颜色主波长角度查找表,可选的,按照此方法构建的颜色主波长查找表,从380nm到700nm对应的角度是逐波长递增的。
本发明实施例在得到各整数纳米波长对应的夹角后,也就得到了整数纳米波长表示的颜色主波长对应的夹角,可建立出各颜色主波长与夹角的对应关系,从而作为图1所示方法中表示的预定的夹角与颜色主波长的对应关系,确定出MODIS图像的各象元的新坐标相对于预定坐标系的X轴的夹角,所对应的颜色主波长。可选的,本发明实施例可基于主波长角度查找表,查找MODIS图像每个象元对应的角度值最近的一个角度,这个角度对应的颜色主波长值即为该象元的颜色主波长值,输出颜色主波长图像。
可选的,图3示出了确定夹角与水色级别的对应关系的方法流程图,值得注意的是,在前文所述的夹角与水色级别的对应关系的表达含义的基础上,图3所示方法仅为实现夹角与水色级别的对应关系的表达含义的一种可选方式;参照图3,该方法可以包括:
步骤S300、确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标;
可选的,水体Forel-Ule比色表为Forel-Ule目视比色法的所采用的比色表(Forel-Ule目视比色法为现有的水质常规监测算法),水体Forel-Ule比色表表达了21个水色级别。本发明实施例可计算水体Forel-Ule比色表表达的21个水色级别在色度图坐标系中的色度坐标。
步骤S310、通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标;
步骤S320、计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角;
步骤S330、确定各水色级别对应的夹角,建立各水色级别与对应的夹角的对应关系。
可选的,本发明实施例可建立水体Forel-Ule比色表表达的21个水色级别与对应的夹角的对应关系,根据该对应关系可构建水色分级角度查找表,可选的,按照此方法构建的水色分级角度查找表,从FU第1级别到第22级别对应的角度是逐渐递增的。
可选的,本发明实施例基于水色分级角度查找表,可查找MODIS图像的每个象元的角度最近的一个角度,这个角度对应的FU分级值(水色分级值)即为该象元的FU分级(水色分级);并且只有当角度值在[-150°,60°]区间内时,才进行查找分级,否则不能进行查找FU分级。
本发明实施例在得到各水色级别与对应的夹角的对应关系后,可将该对应关系作为图1所示方法中表示的预定的夹角与水色级别的对应关系,确定出MODIS图像的各象元的新坐标相对于预定坐标系的X轴的夹角,所对应的水色级别。
本发明实施例提取的水体颜色FU分级值,基于常规水质监测使用的Forel-Ule比色表中的21种水体颜色,能够直观且定量地反映水质状况,实现了基于卫星遥感图像的水体颜色分级,并且为与历史常规监测数据衔接对比提供了可能。
可选的,反映水体颜色的具有RGB波段组合的MODIS图像,可以是包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;可选的,第1波段的波长范围可以为620nm-670nm,第4波段的波长范围可以为545nm-565nm,第3波段的波长范围可以为459nm至479nm;可以看出,第1波段的波长范围位于红光范围,第4波段的波长范围位于绿光范围,第3波段的波长范围位于蓝光范围,该MODIS第1波段,第4波段和第3波段的组合即是MODIS图像RGB真彩色合成的波段组合。
对应的,在将MODIS图像的每个象元的RGB值转换成色度学中的三刺激值时,本发明实施例可基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
可选的,本发明实施例构建的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致的预定坐标系,可以是在CIE标准色度图中,基于色度坐标为(1/3,1/3)的等能白点所构建的二维左手坐标系。
可选的,图4示出了本发明实施例提供的基于MODIS图像的水体颜色检测方法的另一流程图,图4所示方法可作为本发明实施例提供的优选方法,参照图4,该方法可以包括:
步骤S400、获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;
其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
步骤S410、基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
步骤S420、基于各象元的三刺激值,计算每个象元在CIE标准色度学对应的色度图坐标系中的色度坐标;
步骤S430、在CIE标准色度图中,基于色度坐标为(1/3,1/3)的等能白点所构建的二维左手坐标系中,通过坐标转换,计算每个象元的色度坐标在所述二维左手坐标系中的新坐标;
其中,所构建的二维左手坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致。
步骤S440、计算各象元的新坐标相对于所述二维左手坐标系的X轴的夹角;
步骤S450、根据预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别;
其中,本发明实施例可计算色度图的颜色主波长在所述二维左手坐标系中对应的坐标,确定出色度图的颜色主波长在所述二维左手坐标系中对应的坐标,与所述二维左手坐标系的X轴的夹角,从而得到夹角与颜色主波长的对应关系;本发明实施例可通过水体Forel-Ule比色表中各水色级别在所述二维左手坐标系中对应的坐标,确定出各水色级别在所述二维左手坐标系中对应的坐标,与所述二维左手坐标系的X轴的夹角,从而得到夹角与水色级别的对应关系。
步骤S460、根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
可选的,本发明实施例所检测的水体可以是II类水体,可以是内陆水体;显然,本发明实施例所检测的水体类型也可不受限制。
本发明实施例提供的基于MODIS图像的水体颜色检测方法,利用遥感数据对水体水质进行宏观监测,发挥了卫星遥感环境监测的优势;同时基于MODIS反射率遥感图像,建立了水体颜色参量提取的方法,使遥感图像上的水体颜色能够定量化表示,实现了采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的适用。
下面对本发明实施例提供的基于MODIS图像的水体颜色检测装置进行介绍,下文描述的基于MODIS图像的水体颜色检测装置与上文描述的基于MODIS图像的水体颜色检测方法可相互对应参照。
图5为本发明实施例提供的基于MODIS图像的水体颜色检测装置的结构框图,参照图5,该装置可以包括:
图像获取模块100,用于获取反映水体颜色的具有RGB波段组合的MODIS图像;
三刺激值转换模块200,用于将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
色度坐标计算模块300,用于基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
新坐标计算模块400,用于通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,所述预定坐标系的Y轴与色度图坐标系的X轴平行,正向方向一致;
夹角计算模块500,用于计算各象元的新坐标相对于预定坐标系的X轴的夹角;
颜色波长级别确定模块600,用于根据预定的夹角与颜色主波长的对应关系,确定所述夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定所述夹角的水色级别;
水体颜色确定模块700,用于根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
可选的,图6示出了本发明实施例提供的基于MODIS图像的水体颜色检测装置的另一结构框图,结合图5和图6所示,该装置还可以包括:
第一对应关系确定模块800,用于选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2,确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长点在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标,通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标,计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角,确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系;
第二对应关系确定模块900,用于确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标,通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标,计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角,确定各水色级别对应的夹角,建立各水色级别与对应的夹角的对应关系。
可选的,图7示出了本发明实施例提供的图像获取模块100的一种可选结构,参照图7,图像获取模块100可以包括:
图像数据获取单元110,用于获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
可选的,图8示出了本发明实施例提供的三刺激值转换模块200的一种可选结构,参照图8,三刺激值转换模块200可以包括:
转换执行单元210,用于基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
可选的,图9示出了本发明实施例提供的基于MODIS图像的水体颜色检测装置的再一结构框图,结合图5和图9所示,该装置还可以包括:
坐标系建立模块1000,用于在CIE标准色度图中,基于色度坐标为(1/3,1/3)的等能白点构建二维左手坐标系,所构建的二维左手坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致,所构建的二维左手坐标系作为所述预定坐标系。
本发明实施例提供的基于MODIS图像的水体颜色检测装置,利用遥感数据对水体水质进行宏观监测,发挥了卫星遥感环境监测的优势;同时基于MODIS反射率遥感图像,建立了水体颜色参量提取的方法,使遥感图像上的水体颜色能够定量化表示,实现了采用遥感手段对水体颜色进行宏观大范围长时间的遥感检测的适用。
同时,本发明实施例提取的水体颜色FU分级值,基于常规水质监测使用的Forel-Ule比色表中的21种水体颜色,能够直观且定量地反映水质状况,实现了基于卫星遥感图像的水体颜色分级,并且为与历史常规监测数据衔接对比提供了可能。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种基于MODIS图像的水体颜色检测方法,其特征在于,包括:
获取反映水体颜色的具有RGB波段组合的MODIS图像;
将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,所述预定坐标系的Y轴与色度图坐标系的X轴平行,正向方向一致;
计算各象元的新坐标相对于预定坐标系的X轴的夹角;
根据预定的夹角与颜色主波长的对应关系,确定各夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定各夹角的水色级别;
根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
2.根据权利要求1所述的基于MODIS图像的水体颜色检测方法,其特征在于,所述预定的夹角与颜色主波长的对应关系的确定过程包括:
选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2;确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长点在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标;
通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标;
计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角;确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系。
3.根据权利要求1所述的基于MODIS图像的水体颜色检测方法,其特征在于,所述预定的夹角与水色级别的对应关系的确定过程包括:
确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标;
通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标;
计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角;
确定各水色级别对应的夹角,建立各水色级别与对应的夹角的对应关系。
4.根据权利要求1-3任一项所述的基于MODIS图像的水体颜色检测方法,其特征在于,所述获取反映水体颜色的具有RGB波段组合的MODIS图像包括:
获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
5.根据权利要求4所述的基于MODIS图像的水体颜色检测方法,其特征在于,所述将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值包括:
基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
6.根据权利要求1所述的基于MODIS图像的水体颜色检测方法,其特征在于,所述预定坐标系的构建过程包括:
在CIE标准色度图中,基于色度坐标为(1/3,1/3)的等能白点构建二维左手坐标系,所构建的二维左手坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,Y轴与色度图坐标系的X轴平行,正向方向一致,所构建的二维左手坐标系作为所述预定坐标系。
7.一种基于MODIS图像的水体颜色检测装置,其特征在于,包括:
图像获取模块,用于获取反映水体颜色的具有RGB波段组合的MODIS图像;
三刺激值转换模块,用于将所述MODIS图像的每个象元的RGB值转换成色度学中的三刺激值;
色度坐标计算模块,用于基于各象元的三刺激值,计算每个象元在色度图坐标系中的色度坐标;
新坐标计算模块,用于通过坐标转换,计算每个象元的色度坐标在预定坐标系中的新坐标,所述预定坐标系的X轴与色度图坐标系的Y轴平行,正向方向一致,所述预定坐标系的Y轴与色度图坐标系的X轴平行,正向方向一致;
夹角计算模块,用于计算各象元的新坐标相对于预定坐标系的X轴的夹角;
颜色波长级别确定模块,用于根据预定的夹角与颜色主波长的对应关系,确定所述夹角的颜色主波长,及根据预定的夹角与水色级别的对应关系,确定所述夹角的水色级别;
水体颜色确定模块,用于根据所确定的颜色主波长与水色级别,确定MODIS图像所反映的水体颜色。
8.根据权利要求7所述的基于MODIS图像的水体颜色检测装置,其特征在于,还包括:
第一对应关系确定模块,用于选取色度图外围线的设定纳米波长范围,所述设定纳米波长范围包含n个整数纳米波长,一个整数纳米波长表示色度图外围线的一个颜色主波长,n大于2,确定所述设定纳米波长范围中,色度图外围线上各整数纳米波长点在色度图坐标系中的色度坐标,得到n个色度坐标,一个整数纳米波长对应一个色度坐标,通过坐标转换,计算各整数纳米波长对应的色度坐标在预定坐标系中的新坐标,计算各整数纳米波长对应的新坐标,相对于预定坐标系的X轴的夹角,确定各整数纳米波长对应的夹角,建立整数纳米波长表示的颜色主波长,与对应的夹角的对应关系;
第二对应关系确定模块,用于确定水体Forel-Ule比色表中各水色级别在色度图坐标系中的色度坐标,通过坐标转换,计算各水色级别对应的色度坐标在预定坐标系中的新坐标,计算各水色级别对应的新坐标,相对于预定坐标系的X轴的夹角,确定各水色级别对应的夹角,建立各水色级别与对应的夹角的对应关系。
9.根据权利要求7或8所述的基于MODIS图像的水体颜色检测装置,其特征在于,所述图像获取模块包括:
图像数据获取单元,用于获取包含MODIS第1波段,第4波段和第3波段共3个波段的地表反射率水体区域图像数据;其中,第1波段的波长范围为620纳米至670纳米,位于红光范围,第4波段的波长范围为545纳米至565纳米,位于绿光范围,第3波段的波长范围为459纳米至479纳米,位于蓝光范围。
10.根据权利要求9所述基于MODIS图像的水体颜色检测装置,其特征在于,所述三刺激值转换模块包括:
转换执行单元,用于基于CIE标准色度学,将MODIS图像每个象元的第1波段值,第4波段值,第3波段值转换成色度学中的三刺激值X,Y,Z。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693596.8A CN104361602B (zh) | 2014-11-26 | 2014-11-26 | 一种基于modis图像的水体颜色检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693596.8A CN104361602B (zh) | 2014-11-26 | 2014-11-26 | 一种基于modis图像的水体颜色检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104361602A true CN104361602A (zh) | 2015-02-18 |
CN104361602B CN104361602B (zh) | 2017-12-08 |
Family
ID=52528860
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410693596.8A Active CN104361602B (zh) | 2014-11-26 | 2014-11-26 | 一种基于modis图像的水体颜色检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361602B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105488488A (zh) * | 2015-12-10 | 2016-04-13 | 中国科学院遥感与数字地球研究所 | 城市黑臭水体遥感识别方法及装置 |
CN107167431A (zh) * | 2017-05-26 | 2017-09-15 | 中国科学院遥感与数字地球研究所 | 一种基于光谱指数模型的黑臭水体识别方法及系统 |
CN110068821A (zh) * | 2019-05-22 | 2019-07-30 | 广东华宫水利水电建设工程有限公司 | 一种基于遥感的小型水利工程水面面积及蓄水量推求方法 |
CN113109331A (zh) * | 2021-04-12 | 2021-07-13 | 生态环境部卫星环境应用中心 | 水体色度快速检测方法和装置 |
US20230154081A1 (en) * | 2019-07-04 | 2023-05-18 | Zhejiang University | Method for reconstructing geostationary ocean color satellite data based on data interpolating empirical orthogonal functions |
CN117152634A (zh) * | 2023-09-12 | 2023-12-01 | 广东省水利水电科学研究院 | 基于色度指数的多源卫星影像浮水植物识别方法和系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798384A (zh) * | 2012-07-03 | 2012-11-28 | 天津大学 | 一种基于压缩采样的海洋遥感图像水色水温监测方法 |
CN102914505A (zh) * | 2012-09-27 | 2013-02-06 | 国家海洋局第二海洋研究所 | 利用modis水色遥感数据的海洋赤潮藻种判别方法 |
US20130064420A1 (en) * | 2011-09-13 | 2013-03-14 | United States Government, Department of the Navy, Secretary of the Navy, Naval Research Laboratory | Automated system and method for optical cloud shadow detection over water |
CN103389494A (zh) * | 2013-07-24 | 2013-11-13 | 中国科学院南海海洋研究所 | 一种ii类水体水色遥感数据大气校正新方法 |
CN104122213A (zh) * | 2014-06-30 | 2014-10-29 | 苏州科技学院 | 一种测量水质色度的方法 |
CN104132942A (zh) * | 2014-08-11 | 2014-11-05 | 江苏恒创软件有限公司 | 一种基于无人机的水体富营养化检测与处理方法 |
-
2014
- 2014-11-26 CN CN201410693596.8A patent/CN104361602B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130064420A1 (en) * | 2011-09-13 | 2013-03-14 | United States Government, Department of the Navy, Secretary of the Navy, Naval Research Laboratory | Automated system and method for optical cloud shadow detection over water |
CN102798384A (zh) * | 2012-07-03 | 2012-11-28 | 天津大学 | 一种基于压缩采样的海洋遥感图像水色水温监测方法 |
CN102914505A (zh) * | 2012-09-27 | 2013-02-06 | 国家海洋局第二海洋研究所 | 利用modis水色遥感数据的海洋赤潮藻种判别方法 |
CN103389494A (zh) * | 2013-07-24 | 2013-11-13 | 中国科学院南海海洋研究所 | 一种ii类水体水色遥感数据大气校正新方法 |
CN104122213A (zh) * | 2014-06-30 | 2014-10-29 | 苏州科技学院 | 一种测量水质色度的方法 |
CN104132942A (zh) * | 2014-08-11 | 2014-11-05 | 江苏恒创软件有限公司 | 一种基于无人机的水体富营养化检测与处理方法 |
Non-Patent Citations (2)
Title |
---|
BING ZHANG等: "Monitoring water quality of urban water supply sources using optical remote sensing", 《2009 URBAN REMOTE SENSING JOINT EVENT》 * |
GERHARD MEISTER等: "Corrections to the MODIS Aqua Calibration Derived From MODIS Aqua Ocean Color Products", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105488488A (zh) * | 2015-12-10 | 2016-04-13 | 中国科学院遥感与数字地球研究所 | 城市黑臭水体遥感识别方法及装置 |
CN105488488B (zh) * | 2015-12-10 | 2018-12-21 | 中国科学院遥感与数字地球研究所 | 城市黑臭水体遥感识别方法及装置 |
CN107167431A (zh) * | 2017-05-26 | 2017-09-15 | 中国科学院遥感与数字地球研究所 | 一种基于光谱指数模型的黑臭水体识别方法及系统 |
CN107167431B (zh) * | 2017-05-26 | 2019-07-05 | 中国科学院遥感与数字地球研究所 | 一种基于光谱指数模型的黑臭水体识别方法及系统 |
CN110068821A (zh) * | 2019-05-22 | 2019-07-30 | 广东华宫水利水电建设工程有限公司 | 一种基于遥感的小型水利工程水面面积及蓄水量推求方法 |
US20230154081A1 (en) * | 2019-07-04 | 2023-05-18 | Zhejiang University | Method for reconstructing geostationary ocean color satellite data based on data interpolating empirical orthogonal functions |
US11790580B2 (en) * | 2019-07-04 | 2023-10-17 | Zhejiang University | Method for reconstructing geostationary ocean color satellite data based on data interpolating empirical orthogonal functions |
CN113109331A (zh) * | 2021-04-12 | 2021-07-13 | 生态环境部卫星环境应用中心 | 水体色度快速检测方法和装置 |
CN117152634A (zh) * | 2023-09-12 | 2023-12-01 | 广东省水利水电科学研究院 | 基于色度指数的多源卫星影像浮水植物识别方法和系统 |
CN117152634B (zh) * | 2023-09-12 | 2024-06-04 | 广东省水利水电科学研究院 | 基于色度指数的多源卫星影像浮水植物识别方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104361602B (zh) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104361602A (zh) | 一种基于modis图像的水体颜色检测方法及装置 | |
CN101874401B (zh) | 用于测量物体活力的单芯片图像传感器 | |
KR100552695B1 (ko) | 칼라영상의 색조절방법 및 장치 | |
CN105488488B (zh) | 城市黑臭水体遥感识别方法及装置 | |
US8577170B2 (en) | Shadow detection in a single image | |
CN105338262A (zh) | 一种热成像图像处理方法及装置 | |
CN101770646B (zh) | 基于Bayer RGB图像的边缘检测方法 | |
US9036047B2 (en) | Apparatus and techniques for image processing | |
CN103514585B (zh) | 一种图像处理方法 | |
US20150312541A1 (en) | Image pickup device | |
US20050219380A1 (en) | Digital camera for determining chromaticity coordinates and related color temperature of an object and method thereof | |
CN113176211B (zh) | 一种基于hsv颜色空间模型变换的日间云检测方法 | |
JP2000171304A (ja) | カラ―照明色温度の検出方法及び該装置 | |
AU2003244341B2 (en) | Optical sensor, optical sensing method, program, and recording medium | |
KR20080031702A (ko) | 원격탐사영상 구름식별 방법 | |
CN115984711A (zh) | 一种基于卫星遥感的非蓝藻水华监测方法和系统 | |
JP2008292195A (ja) | 物体識別装置及びプログラム | |
CN106815602B (zh) | 一种基于多层次特征描述的跑道fod图像探测方法及装置 | |
US7071965B2 (en) | Artistic color gamut display | |
JP2010272934A (ja) | カラーバランス調整装置 | |
CN104954627A (zh) | 一种信息处理方法及电子设备 | |
CN110390698B (zh) | 一种灰度传感器模组及其数据处理方法和机器人 | |
Crichton et al. | Skin chromaticity gamuts for illumination recovery | |
CN101977328B (zh) | 一种基于补色模式的颜色插值方法 | |
CN201440274U (zh) | 一种目标检测设备及其使用的图像采集装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |