专利名称：具有相干斑特性的任意形状sar海面溢油图像的模拟方法
技术领域：
本发明属于微波海洋遥感技术领域，具体地讲就是涉及具有相干斑特性的任意形状SAR海面溢油图像的模拟方法。
背景技术：
海洋与人类的经济活动息息相关。然而海洋溢油污染具有影响海域范围广、持续时间长、对海洋生物和生态环境破坏大的特点。由于合成孔径雷达(SAR, syntheticaperture radar)具有全天时、全天候、高分辨率和大观测范围的优点，它已成为海洋溢油 监测的一种有利工具，然而高虚警率问题一直是SAR溢油监测系统的难点，尽管已提出了多种改进办法，但该问题仍未得到彻底的解决。相干斑噪声和海面溢油区域形状的复杂性是影响SAR溢油监测系统的高虚警率的重要因素。SAR海面溢油图像的模拟将有助于改进溢油识别方法的研究，从而降低SAR溢油监测系统的高虚警率。从过程上来看，SAR海面溢油图像的模拟会涉及四个方面海面的建模、溢油与海面的相互作用、溢油海面与电磁波的交互作用及SAR成像处理。海面模型有Longuet-Higgins法、分形法、多尺度模型等；Marangoni效应是溢油与海面的相互作用的理论基础；而溢油海面与电磁波的交互作用属于随机粗糙面电磁散射的研究范畴，分为解析法和数值法两大类；SAR成像处理算法有距离-多普勒算法、波数域算法及Chirp-Scaling 算法等。文献[I]采用多尺度海面模型和Marangon效应，建立了溢油海面模型，然后计算溢油海面的雷达截面，最后经过SAR成像处理，获得溢油海面的SAR仿真图像，该方法的不足之处在于计算过程复杂、效率低，另一方面，海面溢油区域的形状没有考虑。参考文献[I] G. Franceschetti, A. Iodice,D. Riccio. SAR Raw Signal Simulationof Oil Slicks in Ocean Environments. IEEE Trans. Geosci. RemoteSensing, 2002，40(9):1935-1949
发明内容
为了克服现有SAR海面溢油图像模拟方法的不足，本发明提供了一种具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，包括以下步骤I.设置SAR图像的参数，计算SAR图像所对应海域的面积A ；SAR图像的参数包括距离向像素点数M、方位向像素点数N、距离分辨率P r、方位分辨率P a ；计算SAR图像所对应海域面积A的方法是:A = MXNX PaX Pr2.根据海况信息，确定SAR图像中每个像素点的归一化灰度值I' (m，n)，且满足K分布；
首先，根据海况信息，确定K分布的尺度参数a和形状参数V ;然后产生具有尺度参数a和形状参数V的K分布的随机数并存储于矩阵/丨.:.、.(/ , ) = 其中，m为SAR图像在距离向上像素点的索引数,它取值为m = 1，2，…，M ；n为SAR图像在方位向上像素点的索引数，它取值为n = I, 2,…，N ;用matIab程序中自带的函数gamrnd()产生MXN的Gamma(伽玛)分布的随机数并存储于矩阵G1 (m，n)，其尺度参数为1，形状参数为I ;用matlab程序中自带的函数gamrnd ()产生MXN的Gamma分布的随机数并存储于矩阵Gv+1(m，n)，其尺度参数为1，形状参数为v+1 ;对矩阵P' a；v(m, n)进行归一化处理后，得到Pa，v(m, n);则SAR图像中每个像素点的归一化灰度值I' (m, n) = Pa，v(m, n)；3.根据相干斑噪声，对SAR图像中每个像素点的灰度值进行修正，得到含有相干斑噪声的SAR海面溢油图像的归一化灰度值I (m，n)；根据相干斑噪声，确定指数分布的参数；然后用matlab程序中自带的函数exprndO产生MXN指数分布的随机数并存储于矩阵B (m，n);含有相干斑噪声K分布的随机数为C' (m，n) =P' a，v(m, n)B(m, n),并对其进行归一化处理,得到C(m, n);则含有相干 斑噪声的SAR海面溢油图像的归一化灰度值I (m, n) = C(m, n)；4.设置SAR图像中溢油面积Atjil和溢油的初始位置p (x0, y。)；5.根据溢油面积，计算溢油区域内像素点的数目；根据溢油面积和SAR图像参数，计算溢油区域内像素点的数目Z = I
PaP16.根据溢油的初始位置和随机游走算法，在SAR图像中生成具有任意形状的溢油区域RoiI ；根据溢油区域的面积Atjil,确定随机游走窗口 W的尺寸wXw ；通常
W = JIoor ^xlO-6，其中，floor表示下取整;步长$ =，随机游走窗口 W的游
U PaPr JUJ
走方向dir = {0,45°，90°，135°，180°，225° }，初始化溢油像素点数j = 0 ;在溢油的初始位置p (x0, y0)处，设置随机游走窗口 W，并把随机游走窗口所在像素点的索引数记录在二维数组Q中，计算当前溢油像素点数j ；随机游走窗口 W按照均匀分布向dir游走一步，且游走步长为S，添加当前随机游走窗口 W像素点的索引数于二维数组Q中，若索引数已写入该数组，则不重复记录，并重新计算当前溢油像素点数j，若j < L，则重复进行该步骤；否则，对二维数组Q按照距离向的索引数进行排序，在SAR图像中生成具有任意形状的溢油区域Rtjil，其索引数已保存在二维数组Q中；7.根据溢油信息，对溢油区域Rtjil中的像素点的灰度值进行修改；根据溢油的信息，确定溢油的均值和方差，然后用matlab程序中自带的函数randn ()产生L个高斯分布的随机数并存储于数组D⑴，I = I, 2，...，L ;然后用matlab程序中自带的函数exp rnd()产生L个指数分布的随机数并存储于数组E (I)中，含有相干斑噪声高斯分布的随机数存储于数组F' (I) =E(I)D(I),并对其归一化处理后，得到数组F(I)，则含有相干斑噪声溢油区域Rtjil的像素点的灰度值已存储于数组F(I)中；按照二维数组Q中的索引数和数组F(l)，逐个修改SAR海面溢油图像中溢油区域Rtjil中像素点的归一化灰度值。
本发明是基于图像的统计分析方法来模拟SAR海面溢油图像。在该方法中，考虑了 SAR成像时产生的相干斑噪声的因素；同时，对溢油区域的几何形状没有增加任何限制因素。另一方面，海况和溢油的影响是通过统计分布来考虑的。因此，本发明具有相干斑噪声和溢油区域任意形状的特点，且适用范围较广。
图I为本发明的流程图。
具体实施例方式下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用 于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。图I为本发明的流程图。具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，包括以下步骤I.设置SAR图像的参数，计算SAR图像所对应海域的面积A ；SAR图像的参数包括距离向像素点数M、方位向像素点数N、距离分辨率P r、方位分辨率P a ；计算SAR图像所对应海域面积A的方法是:A = MXNX p aX P ^2.根据海况信息，确定SAR图像中每个像素点的归一化灰度值I' (m，n)，且满足K分布；首先，根据海况信息，确定K分布的尺度参数a和形状参数V ;然后产生具有尺度参数a和形状参数V的K分布的随机数并存储于矩阵/=...(/ ，《) = 2ayj(il., (/ ,//).其中为SAR图像在距离向上像素点的索引数,它取值为m = 1，2，…，M ;n为SAR图像在方位向上像素点的索引数，它取值为n= I, 2，…，N ;用!1^1:1313程序中自带的函数gamrndO产生MXNGamma分布的随机数并存储于矩阵G1On, n)，其尺度参数为1，形状参数为I ;用matlab程序中自带的函数gamrndO产生MXN的Gamma分布的随机数并存储于矩阵Gv+1 (m, n),其尺度参数为1，形状参数为v+1 ;对矩阵P' a,v(m, n)进行归一化处理后，得到Pa，v(m，n);则SAR图像中每个像素点的归一化灰度值r (m，n) =Pa，v(m，n)；3.根据相干斑噪声，对SAR图像中每个像素点的灰度值进行修正，得到含有相干斑噪声的SAR海面溢油图像的归一化灰度值I (m，n)；根据相干班噪声，确定指数分布的参数；然后用matlab程序中自带的函数exprndO产生MXN指数分布的随机数并存储于矩阵B (m，n);含有相干斑噪声K分布的随机数为C' (m，n) =P' a，v(m, n)B(m, n),并对其进行归一化处理,得到C(m, n);则含有相干斑噪声的SAR海面溢油图像的归一化灰度值I (m, n) = C(m, n)；4.设置SAR图像中溢油面积Atjil和溢油的初始位置p (x0, y0)；5.根据溢油面积，计算溢油区域内像素点的数目；根据溢油面积和SAR图像参数，计算溢油区域内像素点的数目i = i
PaPr6.根据溢油的初始位置和随机游走算法，在SAR图像中生成具有任意形状的溢油区域Ron ；根据溢油区域的面积Atjil,确定随机游走窗口 W的尺寸wXw ；通常w = floor J念xlO-6)，其中，floor表示下取整；步长$ = floor{^ ,随机游走窗口 W的游
走方向dir = {0，45° ,90° , 135°，180°，225° }，初始化溢油像素点数j = O。在溢油的初始位置p (x0, y0)处，设置随机游走窗口 W，并把随机游走窗口所在像素点的索引数记录在二维数组Q中，计算当前溢油像素点数j ;随机游走窗口 W按照均匀分布向dir游走一步,且游走步长为s,添加当前随机游走窗口 W像素点的索引数于二维数组Q 中，若索引数已写入该数组，则不重复记录，并重新计算当前溢油像素点数j，若j < L，则重复进行该步骤；否则，对二维数组Q按照距离向的索引数进行排序，在SAR图像中生成了具有任意形状的溢油区域Rtjil，其索引数已保存在二维数组Q中；7.根据溢油信息，对溢油区域Rtjil中的像素点的灰度值进行修改；根据溢油的信息，确定溢油的均值和方差，然后用matlab程序中自带的函数randn ()产生L个高斯分布的随机数并存储于数组D⑴，I = I, 2，...，L ;然后用matlab程序中自带的函数exp rnd()产生L个指数分布的随机数并存储于数组E (I)中，含有相干斑噪声高斯分布的随机数存储于数组F' (I) =E(I)D(I),并对其归一化处理后，得到数组F(I)，则含有相干斑噪声溢油区域Rtjil的像素点的灰度值已存储于数组F(I)中；按照二维数组Q中的索引数和数组F(l)，逐个修改SAR海面溢油图像中溢油区域Rtjil中像素点的归一化灰度值。
权利要求
1.一种具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于包括以下步骤 步骤I :设置SAR图像的参数，计算SAR图像所对应海域的面积A ； 步骤2 :根据海况信息，确定SAR图像中每个像素点的归一化灰度值I' (m，n)，且满足K分布； 步骤3 :根据相干斑噪声，对SAR图像中每个像素点的灰度值进行修正，得到含有相干斑噪声的SAR海面溢油图像的归一化灰度值I (m，n)； 步骤4 :设置SAR图像中溢油面积Atjil和溢油的初始位置p (X。，y0)； 步骤5 :根据溢油面积，计算溢油区域内像素点的数目； 步骤6 :根据溢油的初始位置和随机游走算法，在SAR图像中生成具有任意形状的溢油区域RoiI ； 步骤7 :根据溢油信息，对溢油区域Rtjil中的像素点的灰度值进行修改。
2.根据权利要求I所述的具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于步骤I按如下过程进行 SAR图像的参数包括距离向像素点数M、方位向像素点数N、距离分辨率Pr、方位分辨率Pa; 计算SAR图像所对应海域面积A的方法是A = MXNXpaXp r0
3.根据权利要求I所述的具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于步骤2按如下过程进行 首先，根据海况信息，确定K分布的尺度参数a和形状参数V ;然后产生具有尺度参数a和形状参数V的K分布的随机数并存储于矩阵(队")=(/；'/,/ )(；,. Amji)，其中，m为SAR图像在距离向上像素点的索引数，它取值为m = 1，2，...，M ；n为SAR图像在方位向上像素点的索引数，它取值为n = I, 2,…，N ;用mat Iab程序中自带的函数gamrnd ()产生MXN的Gamma分布的随机数并存储于矩阵G1On, n)，其尺度参数为1，形状参数为I ;用!^丨1&13程序中自带的函数gamrnd()产生MXN的Ga_a分布的随机数并存储于矩阵Gv+1 (m, n),其尺度参数为1，形状参数为v+1 ;对矩阵P' a，v(m，n)进行归一化处理后，得到Pa，v(m，n);则SAR图像中每个像素点的归一化灰度值I' (m, n) = Pa； v (m, n)。
4.根据权利要求3所述的具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于步骤3按如下过程进行 根据相干斑噪声，确定指数分布的参数；然后用matlab程序中自带的函数exp rnd()产生MXN指数分布的随机数并存储于矩阵B(m，n);含有相干斑噪声K分布的随机数为C' (m，n) =P' a，v(m, n)B(m, n),并对其进行归一化处理,得到C(m, n);则含有相干斑噪声的SAR海面溢油图像的归一化灰度值I (m, n) = C (m, n)。
5.根据权利要求2所述的具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于步骤5按如下过程进行 根据溢油面积和SAR图像参数，计算溢油区域内像素点的数目i
6.根据权利要求2所述的具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，其特征在于步骤6按如下过程进行根据溢油区域的面积Atjil,确定随机游走窗口 W的尺寸wXw ;取
7.根据权利要求6所述的具有相干斑特性的任意形状SAR海面溢油图像的模 拟方法，其特征在于步骤7按如下过程进行 根据溢油的信息，确定溢油的均值和方差，然后用matlab程序中自带的函数randn ()产生L个高斯分布的随机数并存储于数组D (I)，I = 1，2，…，L ;然后用matlab程序中自带的函数exp rnd()产生L个指数分布的随机数并存储于数组E (I)中，含有相干斑噪声高斯分布的随机数存储于数组F' (I) =E(I)D(I),并对其归一化处理后，得到数组F(l)，则含有相干斑噪声溢油区域Rtjil的像素点的灰度值已存储于数组F(I)中；按照二维数组Q中的索引数和数组F(I)，逐个修改SAR海面溢油图像中溢油区域Rtjil中像素点的归一化灰度值。
全文摘要
本发明公开了一种具有相干斑特性的任意形状SAR海面溢油图像的模拟方法，所采用的方法是设置SAR图像的参数，计算SAR图像所对应海域的面积A；根据海况信息，确定SAR图像中每个像素点的归一化灰度值，且满足K分布；根据相干斑噪声，对SAR图像中每个像素点的灰度值进行修正，得到含有相干斑噪声的SAR海面溢油图像的归一化灰度值；设置SAR图像中溢油面积和溢油的初始位置，并计算溢油区域内像素点的数目；根据溢油的初始位置和随机游走算法，在SAR图像中生成具有任意形状的溢油区域；最后，根据溢油信息，对溢油区域中的像素点的灰度值进行修改。本发明具有相干斑噪声和溢油区域任意形状的特点，且适用范围较广。
文档编号G01S13/90GK102798860SQ20121024496
公开日2012年11月28日 申请日期2012年7月13日 优先权日2012年7月13日
发明者杨永红, 奚彩萍, 凌霖 申请人:江苏科技大学