集团网站建设特点 互联网课堂,做类似返利网的网站,wordpress支持字体,上饶网站设计公司《Sentinel-1 Level 1数据处理的详细算法定义》文档定义和描述了Sentinel-1实现的Level 1处理算法和方程#xff0c;以便生成Level 1产品。这些算法适用于Sentinel-1的Stripmap、Interferometric Wide-swath (IW)、Extra-wide-swath (EW)和Wave模式。 今天介绍的内容如下…《Sentinel-1 Level 1数据处理的详细算法定义》文档定义和描述了Sentinel-1实现的Level 1处理算法和方程以便生成Level 1产品。这些算法适用于Sentinel-1的Stripmap、Interferometric Wide-swath (IW)、Extra-wide-swath (EW)和Wave模式。 今天介绍的内容如下 Sentinel-1 Level 1数据处理的详细算法定义一 1 引言1.1目的1.2 范围1.3 关键术语1.4 文档结构 2.文档3.概述3.1 TOPSAR特定处理3.1.1 模式概述3.1.2 科学概述 3.2 双极化处理 4. 预处理算法4.1 原始数据分析 1 引言
1.1目的
本文档描述了Sentinel-1图像处理设施IPF用于生成Sentinel-1 Level 1L1产品的处理算法。这些算法适用于Sentinel-1获取模式Stripmap、干涉宽幅IW、超宽幅EW和Wave模式。在IW和EW模式中数据是使用地形观测渐进扫描SARTOPSAR技术获取的。在本文档中IW和EW将统称为TOPSAR模式。
1.2 范围
本文档满足了Sentinel-1 IPF可交付成果PAL1-1如Sentinel-1产品定义和仪器处理设施开发工作说明书SOW[A-1]和MDA向ESA提出的Sentinel-1 IPF合同变更通知第2号[A-5]中所述。本文档定义了满足Sentinel-1 IPF系统需求文档SRD[A-2]中要求的算法。
本文档展示了将在Sentinel-1 IPF软件中实现的Level 1处理算法和方程。文档[R-1]描述了Sentinel-1 IPF软件设计用于实现Level 1处理算法和方程。
本文档不尝试提供描述算法方程式的详细推导这些在其他文档中有讨论。本文档中讨论的许多Sentinel-1 IPF算法也在[R-5]中有更详细的描述。对于TOPSAR特定算法可以在[R-2]和[R-4]中找到更多的细节、论证和模式设计及算法描述。
由于IPF软件还将支持处理ENVISAT ASAR Stripmap和Wave模式数据本文档也介绍了ASAR特定的处理算法。注意除了预处理外IPF将对Sentinel-1和ENVISAT ASAR使用相同的算法只有相对较小的、局部的调整将反映在软件设计文档[R-1]中。
内部校准算法基于[A-3]的当前版本由[R-13]中包含的信息修改。
1.3 关键术语
• 仰角线仰角或慢时间逐样本变化的样本序列。
• 范围线范围或快时间逐样本变化的样本序列。雷达脉冲的回声是输入到IPF的原始范围线。
• 原始或信号数据SAR回声数据。为了与现有惯例兼容有时使用原始数据如在原始数据处理或原始数据校正中和有时使用信号数据如“信号数据文件”预处理阶段的输出。原始数据通常采用BAQ编码。然而在某些情况下“原始数据”被理解为BAQ解码后的数据如在原始数据校正中该操作是在BAQ解码数据上执行的。BAQ状态的信号原始数据将在每个上下文中明确。
1.4 文档结构
本文档的结构如下
• 第1节介绍目的、范围、关键术语和文档结构。
• 第2节列出适用和参考文档。
• 第3节提供Level 1产品类型和Sentinel-1 IPF支持的处理模块的概述。
• 第4节描述用于Level 0数据的预处理算法。
• 第5节介绍DCE处理算法。
• 第6节介绍SLC图像处理算法。
• 第7节介绍L1后处理算法。
• 第8节介绍支持处理切片数据的支持算法。
• 第9节描述多个处理模块通用的一些较低级别算法或算法。
• 附录A提供了文档中使用的主要符号列表。
• 附录B总结了2D数据的域这些数据是处理各个阶段的对象。
• 附录C提供了Sentinel-1处理器归一化方案的摘要。注意每个归一化步骤也在特定算法的上下文中呈现并在第1至8节中的一个中进行了描述。为了便于识别归一化步骤用缩写NORM标记。
• 附录D提供了文档中引用的辅助和内部参数的摘要。
2.文档
第二章是适用文档和参考文档暂不做翻译。
3.概述
Sentinel-1 SAR可以在以下四种标准获取模式之一下操作见图3-1
• Stripmap模式SM
• 干涉宽幅模式IW
• 超宽幅模式EW
• Wave模式WV Sentinel-1获取模式的简要描述可以在[A-4]中找到。尽管这四种模式都具有不同的属性但它们通常被归为两类单波束单条带的Stripmap模式SM和WV以及多条带TOPSAR模式IW和EW。请注意WV模式是一种采样的Stripmap模式在这种模式下数据的“缩影”会交替地在两个不同的入射角下获取即采用“蛙跳”模式。然而这种特殊的获取方式对所使用的处理算法没有影响每个缩影在单独处理时都被当作Stripmap数据对待。
Sentinel-1仪器处理设施IPF可以从这四种获取模式生成以下Level 1产品
• 斜距单视复数SLC
• 地面范围多视图检测GRD
GRD产品根据其分辨率进一步分类为
• 全分辨率FRSM模式
• 高分辨率HRSM、IW和EW模式
• 中分辨率MRSM、IW、EW和WV模式
每种产品类型的定义在[A-4]中给出。Sentinel-1 Level 1产品家族如图3-2所示。 为了生成这些产品Sentinel-1 IPF由多个软件组件组成如图3-3所示。完整的Sentinel-1 IPF软件设计在[R-1]中描述。 IPF支持处理双极化数据Stripmap或TOPSAR但由于两种极化数据的处理是相同的可能有一些微小的例外通常不会明确提及双极化情况。然而特定于双极化处理的算法方面在第3.2节中介绍。
在Sentinel-1 IPF中实现的算法在本文档中描述如下
• 预处理算法在第4节中描述。
• 多普勒中心频率估计算法在第5节中描述。
• SLC处理算法在第6节中描述。
• Level 1后处理算法在第7节中描述。
• 切片支持算法在第8节中描述。
• 常见和支持算法在第9节中描述。
请注意波模式Wave Mode的多普勒中心频率估计算法被封装在Level 1波模式SLC处理组件内而不是像其他模式那样作为一个单独的进程。
这些组中的算法可能依赖于常见或支持算法这些在第9节中描述。
除了Sentinel-1处理外IPF还支持处理ENVISAT ASAR图像模式IM和波模式WV数据。本文档中描述的Sentinel-1算法也适用于ASAR处理其中处理ASAR IM等同于Sentinel-1 SM处理ASAR WV等同于Sentinel-1 WV。Sentinel-1和ASAR处理算法之间的关键差异与输入原始数据的解码和校准有关这是由于Sentinel-1和ASAR仪器源数据包ISP数据格式的差异所致。
3.1 TOPSAR特定处理
Sentinel-1的IW和EW模式使用TOPSAR成像技术获取宽阔条带的数据分别由3个和5个子条带组成。
3.1.1 模式概述
TOPSAR成像是ScanSAR成像的一种形式其数据是通过对多个相邻子条带循环切换天线波束以突发方式获取的。然而在TOPSAR中对于每个突发波束在方位向上电子地从后向前进行转向。通过这种方式TOPSAR的突发长度比ScanSAR的突发要长得多。图3-4以图解方式表示了在TOPSAR模式下对同样多的子条带进行三个突发的获取强调了天线的旋转。 因此每个目标被照亮的时间较短但是整个天线足迹都会照亮目标因此在TOPSAR中分辨率的损失是通过缩小足迹来实现的而不是像ScanSAR那样通过切割足迹。这对图像质量有几个有益的影响例如在突发内实现均匀的噪声等效Sigma零NESZ和模糊度水平。
TOPSAR特定的概念和算法在文档中描述如下
• 天线转向率以及由于转向产生的直流DC速率在第9.14节中描述。
• 数据所需的方位预处理和后处理在第6.2.4节和第6.4节中描述。
• 在执行以下算法之前对数据进行去斜率处理
o 如第5.2.1节中描述的精细直流DC估计。
o 如第7.2节中描述的地面范围检测GRD方位处理。
下一节提供了TOPSAR原理的数学概述。
3.1.2 科学概述
要更好地理解TOPSAR模式的原理一种方便直观的方法是引入天线转向率和“收缩因子”的概念。 首先考虑在传统Scan-SAR中使用固定天线的突发获取方位天线图案对获取的SAR回波的权重作为方位时间 的函数可以表示为 其中
( L ) 天线长度( \lambda ) 波长( V_r ) 有效传感器速度见9.10节( R_0 ) 零多普勒传感器-目标距离( \theta(\eta) ) 方位角
方位角由下式定义 考虑到现在获取TOPSAR突发的情况方位天线图案的权重函数 Gaap(η) 变为 其中
• kθ 天线转向速度
• Tb 突发长度
请注意为了本节的目的方位天线图案可以被建模为一个sinc函数忽略了由于其转向产生的栅瓣。
结合方程3-2和3-3我们得到以下表达式 比较方程(3-4)和方程(3-1)可以注意到目标被转向天线照亮的地面足迹等同于一个固定天线的足迹但是被“收缩因子”缩小了 因此从方位分辨率的角度来看TOPSAR模式可以被视为一种条带图其等效天线长度是实际长度的k倍因此分辨率是实际分辨率的k倍。
图3-5展示了针对TOPSAR-IW1k4得到的GaapTOPS(η)与相同突发长度的Scan-SAR获取的Gaap SCAN(η)显然k0之间的比较。可以注意到天线图案上的“收缩”效应在TOPSAR情况下方位图案的旁瓣可以看作是由于旋转产生的效果而在Scan-SAR情况下只能看到主瓣。 3.2 双极化处理
双极化数据处理的步骤本质上与单极化数据处理相同通常每种极化都是以相同的方式处理的。然而在处理双极化数据时有一些额外的考虑因素 每种极化都会生成单独的图像并在处理每种极化时应用极化特定的校正和校准因子例如校准脉冲、天线图案等。 必须使用相同的估计多普勒中心频率来处理两种极化以确保两个结果图像能够准确共注册。因此将使用仅从两种极化图像之一估计出的DC。特别是将使用共极化通道因为交叉极化通道可能更嘈杂因此可靠性较低。 交叉极化通道的内部校准复制品重建和PG计算取决于共极化通道的值。
注意通过使用正确的参考复制品进行距离压缩和正确的漂移补偿隐含地应用了通道间的校正如第4.2.1.4节和第4.2.1.5节所述。
4. 预处理算法
Sentinel-1 IPF的预处理模块有四个主要的算法组成部分
如第4.1节所述的原始数据分析。如第4.2节所述的内部校准。如第4.3节所述的下行链路头验证。如第4.4节所述的地形高度函数。
请注意这里描述的算法是专门针对Sentinel-1预处理模块的而不是ASAR预处理模块。ASAR预处理模块执行类似的操作但是专门针对处理ASAR数据实现的。例如ASAR原始数据也必须被解码但它使用的是一种不同的编码算法。
ASAR校准算法在[R-6]中描述原始数据解码算法在[R-8]中描述。
4.1 原始数据分析
为了对原始信号数据的I和Q通道进行校正需要进行原始数据分析。传统的原始数据校正例如在ENVISAT-ASAR和RADARSAT-2的情况下应用包括另见第9.2节
I/Q偏置去除I/Q增益不平衡校正I/Q非正交性校正
然而对于Sentinel-1仪器的接收模块在数字域内执行解调因此不再需要I/Q增益不平衡和I/Q非正交性校正。
对于ASAR数据的原始数据校正所必需的原始数据分析在[R-6]中定义。由于IPF也支持ASAR数据的处理为了完整性在本节中复制了ASAR原始数据分析方案。
尽管对于Sentinel-1来说I/Q增益不平衡和I/Q非正交性校正并不需要但它们将作为可选功能提供通过使用配置输入参数来实现。不过不管校正标志如何本节描述的原始数据分析将对ASAR和Sentinel-1数据执行并且结果将对两种数据都进行报告。
请注意原始数据分析是在BAQ解码后的Level 0数据的一个子集上执行的。为此只有待分析的子集在这个阶段会被BAQ解码。有关BAQ解码的描述以及对整个Level 0数据集进行BAQ解码的阶段请参阅第9.1节。
