专利名称：水下无阴影声纳遥感正射数字图像的计算机生成方法
技术领域：
本发明涉及一种声纳遥感数字图像处理方法，特别是一种声纳遥感数字图像上阴影消除和阴影中像元遥感值恢复的原理方法。
背景技术：
水下地形对声纳遥感信息的影响与改正(包括水下山体阴影、水体背景阴影的消除)一直都是声纳遥感数字图像处理与应用中的难题。现有声纳遥感数字图像处理与应用技术与正射影像处理技术，均没有从根本上解决这些重要的问题。经过查新搜索国内外均没有相关文献和专利报道，本项发明从理论和实践上较好地解决了这一问题。
发明内容
本发明的目的在于针对现有技术存在的缺陷，提供一种水下无阴影声纳遥感正射数字图像的计算机生成方法。
为达到上述发明目的，本发明的构思是本项发明是在GIS支持下建立和利用水中声纳直射声波、水体散射声波照度、水下地形与声纳遥感数字图像相互之间的定量关系，分解原始声纳遥感数字图像为声纳直射声波和水体散射声波遥感数字图像。在此基础上利用数字地形图计算像元声纳直射声波变换系数，水体散射声波变换系数和像元地面辐射参数(散射/直射比，散射/辐射比)进行声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换(如将水下起伏地形上的声纳遥感数字图像转换为水平地面上相应声纳遥感数字图像)。
同时进行水下地形阴影和水体背景阴影的判断与消除，最后完成原始声纳遥感数字图像的水下无阴影正射影像的计算机生成。
声纳遥感数字图像水下无阴影正射影像的计算机生成必须使生成的该图像的每个像元遥感值均严格定量满足声纳遥感图像的点成像原理。
根据上述的发明构思，本发明采用下述技术方案一种水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于首先，建立声纳遥感信息与水下地形、声纳直射声波和水体散射声波之间定量关系的数理模型；继之解析声纳遥感数字图像生成声纳直射声波遥感图像和水体散射声波遥感图像、计算像元声纳直射声波变换系数，水体散射声波变换系数和像元水下地面辐射参数，即散射/直射比，散射/辐射比；建立声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换模型，分别进行声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换，消除水下起伏地形上的辐射差异以及这种差异对声纳遥感数字图像的影响；然后判断和消除遥感数字图像的水下山体与水体背景阴影；在此基础上最后完成声纳遥感数字图像水下无阴影正射影像的计算机生成。
其具体操作步骤如下(1)、声纳遥感数字图像与水下数字地形图配准。
(2)、声纳遥感数字图像的水体修正。
(3)、声纳遥感数字图像上各点声纳高度角、方位角的计算。
(4)、声纳遥感数字图像上水下山体与水体背景阴影的判断。
(5)、观察测量或者计算声纳遥感数字图像上各点水下水平地面上的散射/辐射比。
(6)、计算声纳遥感数字图像上各点水下地面声纳直射声波水下地形影响系数。
(7)、计算声纳遥感数字图像上各点水下地面水体散射声波水下地形影响系数。
(8)、解析声纳遥感图像，生成声纳直射声波遥感图像和水体散射声波遥感图像。
(9)、声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水下水平地形变换。
(10)、声纳遥感数字图像上水下山体与水体背景阴影的消除。
(11)、声纳遥感数字图像水下无阴影正射影像的计算机生成。
上述的骤(1)中的声纳遥感数字图像与数字地形图配准是GIS常规操作。
上述的步骤(2)中的声纳遥感数字图像的水体修正的步骤如下利用声纳遥感图片相同水下地面上相邻两非同类像元DN1和DN2的地面反射率γ1和γ2，按下式计算可以得到离散分布的像元水体程辐射遥感值。
DN(a)＝γ1*(DN1-DN2)/(γ1-γ2)(1)或者对声纳遥感数字图像上水下山体和水体背景阴影中的水体像元遥感值进行分析也可以得到离散分布的像元水体程辐射遥感值；对上述结果，再采用插值算法得到每个像元水体程辐射遥感值；对声纳遥感数字图像上各像元点遥感值DNij进行减去水体程辐射遥感值DNAij的运算，即DNij-DNAij，便完成了声纳遥感数字图像的水体修正。
其中，DNij代表声纳遥感数字图像上各像元点遥感值，DNAij代表水体程辐射遥感值，DNij-DNAij代表声纳遥感数字图像的水体修正值。
上述的步骤(3)中的声纳遥感数字图像上各点声纳高度角、方位角的计算按下式进行计算θij＝arcsin(sinφ*sinδ+cosφ*cosδ*cost(i，j)) (2)A(i，j)＝arcsin(sinθij*sinφ-sinδ)/cosθij*cosφ) (3)δ＝arcsin(sinθ*sinφ-cosθ*cosφ*cosA) (4)t＝arcsin(cosθij*sinA/cosδ)+Δλ (5)θ，A、δ声纳遥感图片注记中声纳下点声纳高度角、方位角和声纳赤角；λ、φ分别为声纳下点地理经纬度。Δλ为像元点对声纳下点的经度增量。
上述的步骤(4)中的声纳遥感数字图像上水下山体与水体背景阴影的判断方法如下利用与声纳遥感数字图像配准的数字地形图和声纳遥感图片像元的声纳的位置参数，即高度角、方位角，进行水下山体与水体背景阴影的判断，其判断的准则是在声纳照射方向该像元点的最大水下地形、水体背景、高度角等于或者大于该像元点的声纳高度角即DH(i，j)≥θij，则该像元点为阴影；反之则不是阴影。
上述的步骤(5)中的观察测量或者计算声纳遥感数字图像上各点水下水平地面上的直、散射辐射比方法如下若声纳遥感图片水下水平地面上有相邻2个同类像元X和Y，像元Y位于阴影中，DN(y)和DN(x)、DNa分别是其遥感值和水体程辐射遥感值，则该处水下水平地面上的散、直射辐射比L＝((DN(y)-DNa)/(DN(x)-DN(y)) (6)上述的步骤(6)中的计算声纳遥感数字图像上各点水下地面直射声波水下地形影响系数按下式计算Fij＝1-tgαij·ctgθij·cosωij，ωij＝ALij-Aij(7)其中，αij代表像元水下地表坡角，ALij代表水下地表坡向，θij代表声纳遥感图片注记中声纳下点声纳高度角，Aij代表声纳方位角，ωij代表声纳方位角与地面坡向之夹角。
上述的步骤(7)中的计算声纳遥感数字图像上各点水下地面水体散射声波水下地形影响系数，即声纳遥感图片像元水下地面水体散射声波照射立体角与水下水平地面上的水体散射声波照射立体角2π之比按下式计算Gij=1-2nπΣk=1nβk,]]>式中n＝2π/Δt (8)其中，βk为第K个方位上的最大平均遮蔽角即立体角。Δt为方位角步长，为360°/n。n为方位划分数，即n为按计算精度要求将360°划分的范围个数。
上述的步骤(8)中的解析声纳遥感图像；生成声纳直射声波遥感图像和水体散射声波遥感图像，方法如下声纳直射声波遥感图像表达式DNSij=(DNij-DNAij)·Fij(Fij+Gij·Lij)---(9)]]>水体散射声波遥感图像表达式DNDij=(DNij-DNAij)·Gij·LijFij+Gij·Lij---(10)]]>i＝1，2…m；j＝1，2…n其中，DNSij代表像元遥感数值中直射声波分量，DNij代表像元遥感数值，DNAij代表水体程辐射遥感值，Lij代表水下地面反射亮度，Fij代表数字地形图上各点水下地面直射声波水下地形影响系数，Gij代表数字地形图上各点水下地面水体散射声波水下地形影响系数，DNDij代表像元遥感数值中散射声波分量，LSij代表水下地面声波谱反射亮度，LDij代表水下地面辐射照度。
上述的步骤(9)中的声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水下地形变换，采用以下公式计算DN′Sij=DNSij/Fij=(DNij-DNAij)/(Fij+Gij·Lij)---(11)]]>DN′Dij=DNDij/Gij=(DNij-DNAij)·Lij/(Fij+Gij·Lij)---(12)]]>i＝1，2…m；j＝1，2…n其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNSij′代表经水体修正后的声纳像元遥感数据中直射声波分量，DNDij′代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
上述的步骤(10)中的声纳遥感数字图像上水下山体与水体背景阴影的消除，采用声纳遥感数字图像上的阴影像元进行以下数字运算DN′ij=(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)---(13)]]>其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNAij代表水体程辐射遥感值。
上述的步骤(11)中的声纳遥感数字图像水下无阴影正射影像的计算机生成，使用以下数字运算完成DNij′=DNSij′+DNDij′---(14)]]>DN′ij=(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)---(15)]]>其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNSij′代表经水体修正后的声纳像元遥感数据中直射声波分量，DNDij′代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
上述的像元水体程辐射DAij，通过数字图像中相关信息中提取或水下地面观测数据计算得出；上述的像元水下地面声纳高度θij、方位角ALij由声纳遥感图片注记给出的相关信息计算得出；上述的像元水下地面坡度αij、水下坡向Aij以及像元经、纬度由声纳遥感图片对应DTM给出。
本发明与现有技术相比，具有如下显而易见的突出实质性特点和显著优点由于声纳遥感成像是在水下声波条件之下，水下声波成份包括声纳直射声波(方向性)分量和水体散射声波(基本上是各向同性)，其相对强度以及在水下起伏地表上的再分配方式，对声纳遥感数据的影响有很大差异。将水体修正后声纳遥感图像解析，分离成声纳直射声波遥感图像和水体散射声波遥感图像，分别进行声纳直射声波和水体散射声波的变换，形成对应的声纳直射声波和水体散射声波的水下地形变换遥感数字图像，最后将二者合成水下声波条件下的水下地形变换声纳遥感数字图像。其理论、方法具有原创性、开拓性。同时，所生成的声纳直射声波和水体散射声波遥感数字图像以及它们的水下地形变换数字图像都是水下无法直接获取的，在理论上和实践上均有重要价值的声纳遥感信息。
由于声纳遥感数字图像的水下地形变换消除了水下地形起伏带来的水下地表声波辐射差异以及水体背景与水下山体阴影的影响，在声纳遥感信息中突出了水下地物波谱特征，使其更具可比性。从根本上改善了声纳遥感数字图像的质量具有重要的科研和应用前景。
本发明可广泛应用于声纳遥感技术的定量研究与应用，如资源环境的声纳遥感调查研究与成像、军事中的伪装与反伪装、水下阴影地物的识别与声纳遥感数字图像的模式识别、水下地物波谱特性的声纳遥感反演、声纳遥感数字图像的虚拟与仿真等方面。
图1为本发明一个实施例的水下无阴影声纳遥感正射数字图像计算机生成流程图。
图2为青岛渤海湾2004.10.10经过水体修正后的声纳遥感数字图像(1024×1024)。
图3为与图2声纳遥感图片匹配的对应数字地形图(立体)。
图4为图2声纳遥感图片解析后生成的声纳直射声波遥感数字图像。
图5为图2声纳解析后生成的水体散射声波遥感数字图像。
图6为水下地形直射声波水下地形变换系数的线性拉伸图。
图7为水体散射声波水下地形变换系数的线性拉伸图。
图8为声纳直射声波水下地形修正系数所成的图像。
图9为水体散射声波水下地形修正系数所成的图像。
图10为计算机生成的水下无阴影声纳遥感数字正射影像。
具体实施例方式
本发明的一个优选实施例结合附图详细说明如下本例为青岛渤海湾2004.10.10的声纳遥感图片以计算机生成无阴影声纳遥感正射数字图像。
参见图1，本水下无阴影声纳正射数字图像的计算机生成方法是首先，建立声纳遥感信息与水下地形、声纳直射声波和水体散射声波之间定量关系的数理模型；继之解析声纳遥感数字图像生成声纳直射声波遥感图像和水体散射声波遥感图像、计算像元声纳直射声波变换系数，水体散射声波变换系数和像元地面辐射参数，即散射/直射比，散射/辐射比；建立声纳直射声波和水体散射声波遥感数字图像的水平地形变换模型，分别进行声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换，消除水下起伏地形上的辐射差异以及这种差异对声纳遥感数字图像的影响；然后判断和消除遥感数字图像的水下山体与水体背景阴影；在此基础上最后完成声纳遥感数字图像水下无阴影正射影像的计算机生成。
其具体操作步骤如下(1)、声纳遥感数字图像与水下数字地形图配准。
(2)、声纳遥感数字图像的水体修正。
(3)、声纳遥感数字图像上各点声纳高度角、方位角的计算。
(4)、声纳遥感数字图像上水下山体与水体背景阴影的判断。
(5)、观察测量或者计算声纳遥感数字图像上各点水平地面上的散射/辐射比。
(6)、计算声纳遥感数字图像上各点水下地面声纳直射声波水下地形影响系数。
(7)、计算声纳遥感数字图像上各点水下地面水体散射声波水下地形影响系数。
(8)、解析声纳遥感图像，生成声纳直射声波遥感图像和水体散射声波遥感图像。
(9)、声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水平地形变换。
(10)、声纳遥感数字图像上水下山体与水体背景阴影的消除。
(11)、声纳遥感数字图像无阴影正射影像的计算机生成。
声纳遥感数字图像的水下地形变换在有多点同步实时辐射观测资料(散射/直射比)、声纳遥感数字图像高精度水体修正、水下典型地物反射率和相应高精度的数字地形图等条件下有最佳的水下地形变换效果。如果上述条件不能满足，仅仅利用一般精度的数字地形图进行声纳遥感数字图像水下地形变换也可以取得令人满意的较好效果。
上述方案中各步骤说明如下1、声纳遥感数字图像与水下数字地形图配准声纳遥感数字图像与水下数字地形图(TEM)配准系GIS常规操作。
2、声纳遥感数字图像的水体修正获得像元水体程辐射遥感值的方法有两种利用声纳遥感图片相同水下地面上相邻两非同类像元DN1和DN2的水下地面反射率γ1和γ2，进行以下计算DN(a)＝γ1*(DN1-DN2)/(γ1-γ2)(1)可以得到离散分布的像元水体程辐射遥感值。
对声纳遥感数字图像上水下山体和水体背景阴影中的水体像元遥感值进行分析也可以得到离散分布的像元水体程辐射遥感值。
对上述结果，再采用插值算法可得到每个像元水体程辐射遥感值。对声纳遥感数字图像上各像元点遥感值进行减去水体程辐射遥感值的运算，即DNij-DNAij，便完成了声纳遥感数字图像的水体修正。
其中，DNij代表声纳遥感数字图像上各像元点遥感值，DNAij代表水体程辐射遥感值，DNij-DNAij代表声纳遥感数字图像的水体修正值。
3、声纳遥感数字图像上各点声纳高度角、方位角的计算按下式进行计算θij＝arcsin(sinφ*sinδ+cosφ*cosδ*cost(i，j)) (2)A(i，j)＝arcsin(sinθij*sinφ-sinδ)/cosθij*cosφ) (3)δ＝arcsin(sinθ*sinφ-cosθ*cosφ*cosA) (4)t＝arcsin(cosθij*sinA/cosδ)+Δλ (5)θ，A、δ声纳遥感图片注记中声纳下点声纳高度角、方位角和声纳赤角；λ、φ分别为声纳下点地理经纬度。Δλ为像元点对声纳下点的经度增量。
4、声纳遥感数字图像上水下山体与水体背景阴影的判断利用与声纳遥感数字图像配准的数字地形图和声纳遥感图片像元的声纳的位置参数，即高度角、方位角，进行水下山体与水体背景阴影的判断，其判断的准则是在声纳照射方向该像元点的最大水下地形、水体背景、高度角等于或者大于该像元点的声纳高度角即DH(i，j)≥θij，则该像元点为阴影；反之则不是阴影。
5、观测或计算声纳遥感数字图像上各点水下水平地面上的直、散射辐射比若声纳遥感图片水平地面上有相邻2个同类像元X和Y，像元Y位于阴影中，DN(y)和DN(x)、DNa分别是其遥感值和水体程辐射辐射遥感值，则该处水下水平地面上的散、直射辐射比L＝((DN(y)-DNa)/(DN(x)-DN(y))(6)6、计算数字地形图上各点地面直射声波水下地形影响系数Fij＝1-tgαij·ctgθij·cosωij，ωij＝ALij-Aij(7)其中，αij代表像元水下地表坡角，ALij代表水下地表坡向，θij代表声纳遥感图片注记中声纳下点声纳高度角，Aij代表声纳方位角，ωij代表声纳方位角与地面坡向之夹角。
7、计算数字地形图上各点水下地面水体散射声波水下地形影响系数，即声纳遥感图片像元水下地面水体散射声波照射立体角与水平地面上的水体散射声波照射立体角2π之比Gij=1-2nπΣk=1nβk,]]>式中n＝2π/Δt (8)其中，βk为第K个方位上的最大平均遮蔽角即立体角。Δt为方位角步长，为360°/n。n为方位划分数，即n为按计算精度要求将360°划分的范围个数。
8、解析声纳遥感图像，生成声纳直射声波遥感图像和水体散射声波遥感图像，方法如下声纳直射声波遥感图像表达式DNSij=(DNij-DNAij)·Fij(Fij+Gij·Lij)---(9)]]>水体散射声波遥感图像表达式DNDij=(DNij-DNAij)·Gij·LijFij+Gij·Lij---(10)]]>i＝1，2…m；j＝1，2…n其中，DNSij代表像元遥感数值中直射声波分量，DNij代表像元遥感数值，DNAij代表水体程辐射遥感值，Lij代表水下地面反射亮度，Fij代表数字地形图上各点地面直射声波水下地形影响系数，Gij代表数字地形图上各点地面水体散射声波水下地形影响系数，DNDij代表像元遥感数值中散射声波分量，LSij代表水下地面声波谱反射亮度，LDij代表水下地面辐射照度。
9、声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水下地形变换，采用以下公式计算DN′Sij=DNSij/Fij=(DNij-DNAij)/(Fij+Gij·Lij)---(11)]]>DN′Dij=DNDij/Gij=(DNij-DNAij)·Lij/(Fij+Gij·Lij)---(12)]]>i＝1，2…m；j＝1，2…n其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNSij′代表经水体修正后的声纳像元遥感数据中直射声波分量，DNDij′代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
10、声纳遥感数字图像上水下山体与水体背景阴影的消除，采用声纳遥感数字图像上的阴影像元进行以下数学运算DN′ij=(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)---(13)]]>其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNAij代表水体程辐射遥感值。
11、声纳遥感数字图像无阴影正射影像的计算机生成，使用以下数字运算完成DNij′=DNSij′+DNDij′---(14)]]>DN′ij=(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)---(15)]]>其中，DNij代表像元遥感数值，DNij′代表经水体修正后的声纳像元遥感数据，DNSij′代表经水体修正后的声纳像元遥感数据中直射声波分量，DNDij′代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
上述的像元水体程辐射DNAij，通过数字图像中相关信息中提取或水下地面观测数据计算得出；上述的像元水下地面声纳高度角θij、方位角ALij由声纳遥感图片注记给出的相关信息计算得出；上述的像元水下地面坡度αij、坡向Aij以及像元经、纬度由声纳遥感图片对应DTM给出。声纳遥感数字图像的水下地形变换在有多点实时辐射观测资料(散射/辐射比)、声纳遥感数字图像高精度水体修正、水下典型地物反射率和相应高精度的数字地形图等条件下有最佳的水下地形变换效果。如果上述条件不能满足，仅仅利用一般精度的数字地形图进行声纳遥感数字图像水下地形变换也可以取得令人满意的较好效果。
由于声纳遥感成像是在水下声波条件之下，水下声波成份包括声纳直射声波(方向性)和水体散声波(基本上是各向同性)分量，其相对强度以及在水下起伏地表上的再分配方式以及对声纳遥感数据的影响都有很大差异。将水体修正后声纳遥感图像解析，分离成声纳直射声波遥感图像和水体散射声波遥感图像，分别进行声纳直射声波和水体散射声波的变换以形成对应的声纳直射声波和水体散射声波的水下地形变换遥感数字图像，最后将二者合成水下声波条件下的地形变换声纳遥感数字图像。其理论、方法具有原创性、和开拓性。同时，所生成的声纳直射声波和水体散射声波遥感数字图像以及它们的水下地形变换数字图像都是水下无法直接获取的，在理论上和实践上均有重要价值的遥感信息。
声纳遥感数字图像的水下地形变换消除了水下地形起伏带来的地表声波辐射差异以及水体背景与水下山体阴影的影响，在声纳遥感信息中突出了水下地物波谱特征，使其更具可比性。从根本上改善了声纳遥感数字图像的质量具有重要的科研和应用前景，与现有技术相比具有的显著的优点与积极的效果。
图2为青岛渤海湾2004.10.10经过水体修正后的声纳遥感数字图像(1024×1024)；图3为与该声纳遥感图片匹配的对应数字地形图(立体)；图4、图5分别为该声纳遥感图片解析后生成的声纳直射声波和水体散射声波遥感数字图像声纳直射声波遥感数字图像(图4)为水下地形影响的声纳直射声波下的遥感数字图像上有明显的阴影分布、山脊两边直射声波辐射有明显差异；水体散射声波遥感数字图像(图5)为水下地形影响的低照度水体散射声波下的遥感数字图像。
图6、图7分别水下地形直射声波地形变换系数和水体散射声波水下地形变换系数的可视化，它们有着十分清晰的物理和地学意义图6为水下地形直射声波地形变换系数的线性拉伸。表达成像瞬间声纳直射声波辐射在水下起伏地面上的真实归一化(水平地面上声纳直射声波辐射照度为1)分布。
图7为水体散射声波水下地形变换系数的线性拉伸。表达成像瞬间水体散射声波辐射在水下起伏地面上的真实归一化(水平地面上水体散射声波辐射照度为1)分布。
直射声波水下地形变换系数和水体散射声波水下地形变换系数及其可视化在数字沙盘、地图晕染和数字可视化等方面有广泛应用。
图8、图9分别为声纳直射声波和水体散射声波遥感数字图像所对应的水平地形变换的图像声纳直射声波(非偏振声波)遥感数字图像的水下地形变换图像(图8)上阴影已经消除，原来隐藏在阴影中的水体清晰可见，山脊两边直射声波辐射差异已不明显，水下起伏山地已变成平原。与原来的声纳遥感数字图像相比声纳直射声波遥感数字图像的水下地形变换图像显然有利于计算机的自动识别与分类，其图像质量有重大的改善和提高。
水体散射声波遥感数字图像的水下地形变换图像(图9)为了消除水下地形影响在各向均匀的水体散射声波(偏振声波)条件下的成像的遥感数字图像。
图10为计算机生成的水下无阴影声纳遥感数字正射影像(图8和图9的合成)。
声纳遥感数字图像的水下地形变换在有多点同步实时辐射观测资料(散射/直射比)、声纳遥感数字图像高精度水体修正、水下典型地物反射率和相应高精度的数字地形图等条件下有最佳的水下地形变换效果。如果上述条件不能满足，仅仅利用一般精度的数字地形图进行声纳遥感数字图像地形变换也可以取得令人满意的较好效果。
权利要求
1.一种水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于首先，建立声纳遥感信息与水下地形、声纳直射声波和水体散射声波之间定量关系的数理模型；继之解析声纳遥感数字图像生成声纳直射声波遥感图像和水体散射声波遥感图像、计算像元声纳直射声波变换系数，水体散射声波变换系数和像元水下地面辐射参数，即散射/直射比，散射/辐射比；建立声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换模型，分别进行声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换，消除水下起伏地形上的辐射差异以及这种差异对声纳遥感数字图像的影响；然后判断和消除遥感数字图像的水下山体与水体背景阴影；在此基础上最后完成声纳遥感数字图像水下无阴影正射影像的计算机生成；其具体操作步骤如下(1)、声纳遥感数字图像与水下数字地形图配准；(2)、声纳遥感数字图像的水体修正；(3)、声纳遥感数字图像上各点声纳高度角、方位角的计算；(4)、声纳遥感数字图像上水下山体与水体背景阴影的判断；(5)、观察测量或者计算声纳遥感数字图像上各点水下水平地面上的散射/辐射比；(6)、计算声纳遥感数字图像上各点水下地面声纳直射声波水下地形影响系数；(7)、计算声纳遥感数字图像上各点水下地面水体散射声波水下地形影响系数；(8)、解析声纳遥感图像，生成声纳直射声波遥感图像和水体散射声波遥感图像；(9)、声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水下水平地形变换；(10)、声纳遥感数字图像上水下山体与水体背景阴影的消除；(11)、声纳遥感数字图像水下无阴影正射影像的计算机生成；
2.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(2)中的声纳遥感数字图像的水体修正的步骤如下(1)利用声纳遥感图片相同水下地面上相邻两非同类像元DN1和DN2的地面反射率γ1和γ2，按下式计算可以得到离散分布的像元水体程辐射遥感值；DN(a)＝γ1*(DN1-DN2)/(γ1-γ2)，或者对声纳遥感数字图像上水下山体和水体背景阴影中的水体像元遥感值进行分析得到离散分布的像元水体程辐射遥感值；(2)对上述结果，再采用插值算法得到每个像元水体程辐射遥感值；(3)对声纳遥感数字图像上各像元点遥感值DNij进行减去水体程辐射遥感值DNAij的运算，即DNij-DNAij，便完成了声纳遥感数字图像的水体修正，其中，DNij代表声纳遥感数字图像上各像元点遥感值，DNAij代表水体程辐射遥感值，DNij-DNAij代表声纳遥感数字图像的水体修正值。
3.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(3)中的声纳遥感数字图像上各点声纳高度角、方位角的计算按下式进行计算θij＝arcsin(sinφ*sinδ+cosφ*cosδ*cost(i，j))，A(i，j)＝arcsin(sinθij*sinφ-sinδ)/cosθij*cosφ)，δ＝arcsin(sinθ*sinφ-cosθ*cosφ*cosA)，t＝arcsin(cosθij*sinA/cosδ)+Δλ，其中θ，A、δ分别为声纳遥感图片注记中声纳下点声纳高度角、方位角和声纳赤角；λ、φ分别为声纳下点地理经纬度，Δλ为像元点对声纳下点的经度增量。
4.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(4)中的声纳遥感数字图像上水下山体与水体背景阴影的判断方法如下利用与声纳遥感数字图像配准的数字地形图和声纳遥感图片像元的声纳的位置参数，即高度角、方位角，进行水下山体与水体背景阴影的判断，其判断的准则是在声纳照射方向该像元点的最大水下地形、水体背景、高度角等于或者大于该像元点的声纳高度角即DH(i，j)≥θij，则该像元点为阴影；反之则不是阴影。
5.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(5)中的观测或计算声纳遥感数字图像上各点水下水平地面上的直、散射辐射比方法如下若声纳遥感图片水下水平地面上有相邻2个同类像元X和Y，像元Y位于阴影中，DN(y)和DN(x)、DNa分别是其遥感值和水体程辐射遥感值，则该处水下水平地面上的散、直射辐射比L＝(DN(y)-DNa)/(DN(x)-DN(y))。
6.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(6)中的计算数字地形图上各点水下地面直射声波水下地形影响系数按下式计算Fij＝1-tgαij·ctgθij·cosωij，ωij＝ALij-Aij，其中，αij代表像元水下地表坡角，ALij代表水下地表坡向，θij代表声纳遥感图片注记中声纳下点声纳高度角，Aij代表声纳方位角，ωij代表声纳方位角与水下地面坡向之夹角。
7.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(7)中的计算数字地形图上各点水下地面水体散射声波水下地形影响系数，即声纳遥感图片像元地面水体散射声波照射立体角与水平地面上的水体散射声波照射立体角2π之比按下式计算Gij=1-2nπΣk=1nβk,]]>式中n＝2π/Δt，其中，βk为第K个方位上的最大平均遮蔽角即立体角，Δt为方位角步长，为360°/n，n为方位划分数，即n为按计算精度要求将360°划分的范围个数。
8.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(8)中的解析声纳遥感图像，生成声纳直射声波遥感图像和水体散射声波遥感图像的方法如下声纳直射声波遥感图像表达式DNSij=(DNij-DNAij)·Fij(Fij+Gij·Lij),]]>水体散射声波遥感图像表达式DNDij=(DNij-DNAij)·Gij·LijFij+Gij·Lij,]]>i＝1，2…m；j＝1，2…n，其中，DNSij代表像元遥感数值中直射声波分量，DNij代表像元遥感数值，DNAij代表水体程辐射遥感值，Lij代表水下地面反射亮度，Fij代表数字地形图上各点水下地面直射声波水下地形影响系数，Gij代表数字地形图上各点水下地面水体散射声波水下地形影响系数，DNDij代表像元遥感数值中散射声波分量，LSij代表水下地面声波谱反射亮度，LDij代表水下地面辐射照度。
9.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(9)中的声纳直射声波遥感数字图像和水体散射声波遥感数字图像的水下地形变换，采用以下公式计算DN′Sij＝DNSij/Fij＝(DNij-DNAij)/(Fij+Gij·Lij)，DN′Dij＝DNDij/Gij＝(DNij-DNAij)·Lij/(Fij+Gij·Lij)，i＝1，2…m；j＝1，2…n，其中，DNij代表像元遥感数值，DN′ij代表经水体修正后的声纳像元遥感数据，DN′Sij代表经水体修正后的声纳像元遥感数据中直射声波分量，DN′Dij代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
10.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(10)中的声纳遥感数字图像上水下山体与水体背景阴影的消除，采用声纳遥感数字图像上的阴影像元进行以下数字运算DN′ij＝(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)，其中，DNij代表像元遥感数值，DN′ij代表经水体修正后的声纳像元遥感数据，DNAij代表水体程辐射遥感值。
11.根据权利要求1所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于所述的步骤(11)中的声纳遥感数字图像水下无阴影正射影像的计算机生成，使用以下数字运算完成DN′ij＝DN′Sij+DN′Dij，DN′ij＝(DNij-DNAij)·(1+Lij)/(Fij+Gij·Lij)，其中，DNij代表像元遥感数值，DN′ij代表经水体修正后的声纳像元遥感数据，DN′Sij代表经水体修正后的声纳像元遥感数据中直射声波分量，DN′Dij代表经水体修正后的声纳像元遥感数据中散射声波分量，DNDij代表像元遥感数值中散射声波分量，DNSij代表像元遥感数值中直射声波分量，DNAij代表水体程辐射遥感值。
12.根据权利要求2或3或6所述的水下无阴影声纳遥感正射数字图像的计算机生成方法，其特征在于像元水体程辐射DNAij，通过数字图像中相关信息中提取或水下地面观测数据计算得出；所述的像元地面声纳高度θij、方位角ALij由声纳遥感图片注记给出的相关信息计算得出；所述的像元水下地面坡度αij、坡向Aij以及像元经、纬度由声纳遥感图片对应DTM给出。
全文摘要
本发明涉及一种水下无阴影声纳遥感正射数字图像的计算机生成方法。本方法是首先，建立声纳遥感信息与水下地形、声纳直射声波和水体散射声波之间定量关系的数理模型；继之解析声纳遥感数字图像生成声纳直射声波遥感图像和水体散射声波遥感图像、计算像元声纳直射声波变换系数，水体散射声波变换系数和像元地面辐射参数，即散射/直射比，散射/辐射比；建立声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换模型，分别进行声纳直射声波和水体散射声波遥感数字图像的水下水平地形变换，消除水下起伏地形上的辐射差异以及这种差异对声纳遥感数字图像的影响；然后判断和消除声纳遥感数字图像的水下山体与水体背景阴影；在此基础上最后完成声纳遥感数字图像水下无阴影正射影像的计算机生成。本发明的理论、方法具有原创性、开拓性，可广泛应用于声纳遥感技术的定量研究和应用。
文档编号G01S15/00GK101034476SQ20071003870
公开日2007年9月12日 申请日期2007年3月29日 优先权日2007年3月29日
发明者师彪, 李先华, 张凌峰, 唐振军, 魏为民 申请人:上海大学