CN107064890A - 一种脉冲体制雷达海冰探测能力评估方法 - Google Patents
一种脉冲体制雷达海冰探测能力评估方法 Download PDFInfo
- Publication number
- CN107064890A CN107064890A CN201710232310.XA CN201710232310A CN107064890A CN 107064890 A CN107064890 A CN 107064890A CN 201710232310 A CN201710232310 A CN 201710232310A CN 107064890 A CN107064890 A CN 107064890A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- msub
- mfrac
- mtr
- 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
- 238000000034 method Methods 0.000 title claims abstract description 12
- 239000013535 sea water Substances 0.000 claims abstract description 25
- 230000010287 polarization Effects 0.000 claims abstract description 22
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 230000005540 biological transmission Effects 0.000 claims abstract description 8
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 230000035699 permeability Effects 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000005291 magnetic effect Effects 0.000 claims description 5
- 239000005457 ice water Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 230000005294 ferromagnetic effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003305 oil spill Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种脉冲体制雷达海冰探测能力评估方法,对于脉冲体制雷达,设其脉冲信号宽度为τ,电磁波传播速度为c0,则受脉宽限制的可探测海冰厚度为其中c0=3×108米/秒,θ2为海冰中电磁波透射角。海冰识别度用Dice表示,物理意义为雷达散射图像中海冰区域与海水区域的区分程度,水平极化和垂直极化时海冰识别度分别为 本发明基于海冰介电常数模型,考虑气‑冰界面与海水‑海冰界面之间电磁波多次散射,要评估因子为海冰可探测厚度和海冰识别度。
Description
技术领域
本发明属于海洋技术领域,具体涉及一种海冰可探测厚度及海冰识别度模拟方法。
背景技术
自上世纪80年代以来,由于全球变暖,海冰融化加快,北极海冰对全球气候的影响有放大作用,利用空基和天基雷达对极区海冰进行探测具有重要的科学意义和应用价值。
合成孔径雷达(简记为SAR)具有全天候、可穿透云层等优点,自1978年首次发射升空以来,被应用于对地观测特别是对海洋观测,包括海面风场、海浪参数提取、海流锋面检测、海表面溢油、海上目标识别、海冰运动及分类等应用领域。海洋表面各种洋面特征引起的海面粗糙度变化,对SAR海面成像图像雷达散射截面产生影响,使得SAR的海洋特征检测应用成为可能,因此对于不同介质的海面电磁散射也是一直以来的研究热点。其中海冰的SAR遥感监测越来越受到重视,研究热点包括海冰分类和海冰厚度探测。海冰厚度反演研究一直是气候模式研究所急需解决的技术瓶颈,最近的研究主要是利用海冰薄层的盐度和卤水体积随冰厚的变化经验关系进行海冰厚度探索研究,但该方法仅适用于小于薄冰厚度探测(小于0.4米)。目前对于海冰的探测能力尚缺乏定量理论研究,我们提出了一种脉冲体制雷达海冰探测能力评估模型,发明了一种脉冲体制雷达海冰探测能力评估方法。
发明内容
本发明的目的在于克服现有技术的不足,提供一种脉冲体制雷达海冰探测能力评估方法,是基于分层多次电磁散射机制的一种归一化雷达散射截面(NRCS)、可探测海冰厚度以及海冰识别度的模拟方法。
本发明采用的技术方案为:一种脉冲体制雷达海冰探测能力评估方法,包括以下步骤:
步骤1:介质1、2、3分别为空气、海冰和海水,它们的介电常数、磁导率和电导率分别为(ε1,μ1,σ1)、(ε2,μ2,σ2)、(ε3,μ3,σ3);一般来说,除了铁磁质外的介质,磁导率都近似等于真空中的磁导率,即有μ1=μ2=μ3=μ0。空气中ε1=ε0,σ1=0,其中真空中磁导率μ0=4π×10-7亨/米(H/m)、介电常数ε0=8.85419×10-12(F/m),海冰电导率为σ2=ωε″ε0,其中,电磁波角频率为ω=2πf,f为电磁波频率。海水电导率为:
式中,SW为海水的盐度,t为海水的温度。海水介电常数Debye模型函数为:
式中,
与温度和盐度相关的弛豫时间τ为:
ε*=a′+a″SW,式中,
步骤2:电磁波部分能量透射入海冰中传播,入射角为θ1,透射角为:
透入海水中的透射角为:
海冰中的电磁波衰减系数为:
步骤3:本发明方法中的受脉宽限制的可探测海冰厚度为:
其中电磁波传播速度c0=3×108米/秒。
步骤4:本发明中,海冰识别度用Dice表示,物理意义为雷达散射图像中海冰区域与海水区域的区分程度,水平极化(HH)和垂直极化(VV)时海冰识别度分别为:
其中,SHH,water,SVV,water分别为水平极化和垂直极化下海水的雷达散射系数。
分别水平极化和垂直极化下海冰的雷达散射系数。
其中,水平极化时冰-气界面的反射系数:
垂直极化时冰-气界面的反射系数:
其中,
有益效果:本发明基于海冰介电常数模型,考虑气-冰界面与海水-海冰界面之间电磁波多次散射,提出一种脉冲体制雷达海冰探测能力评估模型,主要评估因子为海冰可探测厚度和海冰识别度。
附图说明
图1为平面电磁波入射厚度为h的海冰时多次反射和透射示意图。
图2为雷达海冰识别度随雷达频率的变化。
图3为雷达海冰识别度随雷达入射角的变化。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明。
一种脉冲体制雷达海冰探测能力评估方法,包括以下步骤:
步骤1:如图1所示,介质1、2、3分别为空气、海冰和海水,它们的介电常数、磁导率和电导率分别为(ε1,μ1,σ1)、(ε2,μ2,σ2)、(ε3,μ3,σ3);一般来说,除了铁磁质外的介质,磁导率都近似等于真空中的磁导率,即有μ1=μ2=μ3=μ0。空气中ε1=ε0,σ1=0,其中真空中磁导率μ0=4π×10-7亨/米(H/m)、介电常数ε0=8.85419×10-12(F/m),海冰电导率为σ2=ωε″ε0,其中,电磁波角频率为ω=2πf,f为电磁波频率。海水电导率为:
式中,SW为海水的盐度,t为海水的温度。海水介电常数Debye模型函数为:
式中,
与温度和盐度相关的弛豫时间τ为:
ε*=a′+a″SW,式中,
步骤2:电磁波部分能量透射入海冰中传播,入射角为θ1,透射角为:
透入海水中的透射角为:
海冰中的电磁波衰减系数为:
步骤3:本发明方法中的受脉宽限制的可探测海冰厚度为:
其中电磁波传播速度c0=3×108米/秒。
步骤4:如图2和图3所示,本发明中,海冰识别度用Dice表示,物理意义为雷达散射图像中海冰区域与海水区域的区分程度,水平极化(HH)和垂直极化(VV)时海冰识别度分别为:
其中,SHH,water,SVV,water分别为水平极化和垂直极化下海水的雷达散射系数。
分别水平极化和垂直极化下海冰的雷达散射系数。
其中,水平极化时冰-气界面的反射系数:
垂直极化时冰-气界面的反射系数:
其中,
以上结合附图对本发明的实施方式做出详细说明,但本发明不局限于所描述的实施方式。对本领域的普通技术人员而言,在本发明的原理和技术思想的范围内,对这些实施方式进行多种变化、修改、替换和变形仍落入本发明的保护范围内。
Claims (1)
1.一种脉冲体制雷达海冰探测能力评估方法,其特征在于:包括以下步骤:
步骤1:介质1、2、3分别为空气、海冰和海水,它们的介电常数、磁导率和电导率分别为(ε1,μ1,σ1)、(ε2,μ2,σ2)、(ε3,μ3,σ3);磁导率都近似等于真空中的磁导率,即有μ1=μ2=μ3=μ0;空气中ε1=ε0,σ1=0,其中真空中磁导率μ0=4π×10-7H/m、介电常数ε0=8.85419×10- 12F/m,海冰电导率为σ2=ωε″ε0,其中,电磁波角频率为ω=2πf,f为电磁波频率;海水电导率为:
<mrow>
<mfenced open='' close=''>
<mtable>
<mtr>
<mtd>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mn>0.086374</mn>
<mo>+</mo>
<mn>0.030606</mn>
<mi>t</mi>
<mo>-</mo>
<msup>
<mrow>
<mn>0.0004121</mn>
<mi>t</mi>
</mrow>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>W</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0.077454</mn>
<mo>+</mo>
<mn>0.001687</mn>
<mi>t</mi>
<mo>+</mo>
<msup>
<mrow>
<mn>0.000019371</mn>
<mi>t</mi>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,SW为海水的盐度,t为海水的温度;海水介电常数Debye模型函数为:
<mrow>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mi>&infin;</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>&epsiv;</mi>
<mo>*</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mi>&infin;</mi>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>&epsiv;</mi>
<mo>*</mo>
</msup>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mi>&infin;</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;f&epsiv;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>&infin;</mi>
</msub>
<mo>=</mo>
<mn>6.4857</mn>
<mo>-</mo>
<mn>0.04203</mn>
<mi>t</mi>
<mo>-</mo>
<mn>0.006588</mn>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>0.0006492</mn>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
<mo>-</mo>
<mn>1.2328</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<mn>5.0433</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>8</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>5</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
与温度和盐度相关的弛豫时间τ为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mn>17.03</mn>
<mo>-</mo>
<mn>0.66651</mn>
<mi>t</mi>
<mo>+</mo>
<mn>5.1482</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1.2145</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>5.0325</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<mn>5.8272</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>7</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>5</mn>
</msup>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>W</mi>
</msub>
<mo>(</mo>
<mo>-</mo>
<mn>6.772</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>2.357</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msup>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>5.075</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>6.3983</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<mn>2.463</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>6</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>4</mn>
</msup>
<mo>-</mo>
<mn>3.0676</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>8</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>5</mn>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
ε*=a′+a″SW,式中,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>a</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mn>81.82</mn>
<mo>-</mo>
<mn>6.0503</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>t</mi>
<mo>-</mo>
<mn>3.1661</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>3.1097</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
<mo>-</mo>
<mn>1.1791</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<mn>1.4838</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>6</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>5</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>a</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mo>=</mo>
<mn>0.12544</mn>
<mo>+</mo>
<mn>9.4037</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<mi>t</mi>
<mo>-</mo>
<mn>9.5551</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>9.0888</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
<mo>-</mo>
<mn>3.6011</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>6</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<mn>4.71</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>8</mn>
</mrow>
</msup>
<msup>
<mi>t</mi>
<mn>5</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤2:电磁波部分能量透射入海冰中传播,入射角为θ1,透射角为:
<mrow>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>a</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
</mfrac>
</msqrt>
<msub>
<mi>sin&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
透入海水中的透射角为:
<mrow>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mi>a</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
</mfrac>
</msqrt>
<msub>
<mi>sin&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
海冰中的电磁波衰减系数为:
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>&omega;</mi>
<msqrt>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msubsup>
<mi>&sigma;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>&epsiv;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3:所述评估方法中的受脉宽限制的可探测海冰厚度为:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中电磁波传播速度c0=3×108米/秒;
步骤4:所述评估方法中,海冰识别度用Dice表示,物理意义为雷达散射图像中海冰区域与海水区域的区分程度,水平极化和垂直极化时海冰识别度分别为:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
<mo>,</mo>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
<mo>,</mo>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
<mo>,</mo>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
<mo>,</mo>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,SHH,water,SVV,water分别为水平极化和垂直极化下海水的雷达散射系数;
分别水平极化和垂直极化下海冰的雷达散射系数;
其中,水平极化时冰-气界面的反射系数:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>b</mi>
<mi>e</mi>
<mi>k</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>a</mi>
<mi>e</mi>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
垂直极化时冰-气界面的反射系数:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>g</mi>
<mi>l</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>c</mi>
<mi>g</mi>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710232310.XA CN107064890B (zh) | 2017-04-11 | 2017-04-11 | 一种脉冲体制雷达海冰探测能力评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710232310.XA CN107064890B (zh) | 2017-04-11 | 2017-04-11 | 一种脉冲体制雷达海冰探测能力评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064890A true CN107064890A (zh) | 2017-08-18 |
CN107064890B CN107064890B (zh) | 2019-10-25 |
Family
ID=59603197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710232310.XA Active CN107064890B (zh) | 2017-04-11 | 2017-04-11 | 一种脉冲体制雷达海冰探测能力评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064890B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679476A (zh) * | 2017-09-26 | 2018-02-09 | 南京大学 | 一种海冰类型遥感分类方法 |
RU2723437C1 (ru) * | 2019-09-03 | 2020-06-11 | Открытое акционерное общество "Авангард" | Способ обнаружения и высокоточного определения параметров морских ледовых полей и радиолокационная система для его реализации |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3665466A (en) * | 1970-03-20 | 1972-05-23 | Exxon Production Research Co | Determination of ice thickness |
WO2003048803A1 (fr) * | 2001-11-07 | 2003-06-12 | Communications Research Laboratory,Independent Administrative Institution | Procede d'observation de glace marine |
JP2005291782A (ja) * | 2004-03-31 | 2005-10-20 | National Institute Of Information & Communication Technology | Sarによる氷厚推定方法 |
CN101105395A (zh) * | 2007-08-01 | 2008-01-16 | 大连海事大学 | 雷达海冰厚度测量仪 |
CN203720349U (zh) * | 2014-02-20 | 2014-07-16 | 中船重工鹏力(南京)大气海洋信息系统有限公司 | 一种实现海洋雷达多种扫描方式的扫描系统 |
-
2017
- 2017-04-11 CN CN201710232310.XA patent/CN107064890B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3665466A (en) * | 1970-03-20 | 1972-05-23 | Exxon Production Research Co | Determination of ice thickness |
WO2003048803A1 (fr) * | 2001-11-07 | 2003-06-12 | Communications Research Laboratory,Independent Administrative Institution | Procede d'observation de glace marine |
JP2005291782A (ja) * | 2004-03-31 | 2005-10-20 | National Institute Of Information & Communication Technology | Sarによる氷厚推定方法 |
CN101105395A (zh) * | 2007-08-01 | 2008-01-16 | 大连海事大学 | 雷达海冰厚度测量仪 |
CN203720349U (zh) * | 2014-02-20 | 2014-07-16 | 中船重工鹏力(南京)大气海洋信息系统有限公司 | 一种实现海洋雷达多种扫描方式的扫描系统 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679476A (zh) * | 2017-09-26 | 2018-02-09 | 南京大学 | 一种海冰类型遥感分类方法 |
CN107679476B (zh) * | 2017-09-26 | 2020-10-09 | 南京大学 | 一种海冰类型遥感分类方法 |
RU2723437C1 (ru) * | 2019-09-03 | 2020-06-11 | Открытое акционерное общество "Авангард" | Способ обнаружения и высокоточного определения параметров морских ледовых полей и радиолокационная система для его реализации |
Also Published As
Publication number | Publication date |
---|---|
CN107064890B (zh) | 2019-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Moum et al. | Structure and generation of turbulence at interfaces strained by internal solitary waves propagating shoreward over the continental shelf | |
Callaghan et al. | Two regimes of laboratory whitecap foam decay: Bubble-plume controlled and surfactant stabilized | |
CN103439694B (zh) | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 | |
CN102798856B (zh) | 一种基于小波包的24GHz LFMCW雷达测距方法 | |
Haas | Dynamics versus thermodynamics: The sea ice thickness distribution | |
CN104991288A (zh) | 一种基于声呐电磁协同探测技术的海洋沉潜油检测系统及方法 | |
CN107064890B (zh) | 一种脉冲体制雷达海冰探测能力评估方法 | |
CN103605118A (zh) | 一种利用极地探冰雷达提取极地冰层位方法 | |
CN103323817B (zh) | 一种机载合成孔径雷达海面风矢量反演方法 | |
CN103630896A (zh) | 非线性调频尺度成像的方法 | |
CN102495293A (zh) | 基于系统辨识理论的机场道面介质层电磁特性的反演方法 | |
CN113686964A (zh) | 一种基于泄漏模态声波导特性的海冰厚度观测方法 | |
CN103105603B (zh) | 一种x波段测波雷达海流反演预处理方法 | |
CN106597546B (zh) | 一种探测河底表层泥沙物性的方法 | |
CN104280737B (zh) | 加权宽带时反算子分解声成像方法 | |
CN107300694A (zh) | 一种基于电磁波透射系数的未知墙体参数估计方法 | |
CN104482853B (zh) | 一种现场实时测量海上油膜厚度装置及方法 | |
CN114878696A (zh) | 一种弧形复合材料层合板分层损伤的识别方法 | |
Liu et al. | A feasibility study to identify ice types by measuring attenuation of ultrasonic waves for aircraft icing detection | |
CN103744119B (zh) | 一种建筑体裂缝的三维特征探测方法及系统 | |
WO2017041243A1 (en) | The calculation method of wave reflective index on interface | |
CN105334507B (zh) | 基于极化多特征的对海面漂浮雷达目标的检测方法 | |
CN102866171B (zh) | 雪层覆盖下海冰的后向电磁散射系数探测模块 | |
Anderson | Monitoring the marginal ice zone with HF radar | |
Duarte et al. | Ultrasonic time of flight estimation for wind speed measurement based on time-frequency domain using STFT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |