专利名称：星载斜视sar原始数据实用化频域压缩比特分配方法
技术领域：
本发明涉及星载斜视SAR原始数据压缩技术领域，是一种星载斜视SAR原始数据实用化频域压缩比特分配方法，以解决星载斜视SAR多普勒中心非零以及二维频谱的能量分布沿方位向存在倾斜的特性导致传统频域压缩算法失效的问题。
背景技术：
星载合成孔径雷达(SAR)频域类压缩算法的主要思想是将SAR原始数据变换到二维频域，利用频域的去相关性以及二维频谱的能量分布特性，在满足传输数据率的前提下提高信噪比，代表性的压缩算法包括FFT-BAQ (Fast Fourier Transform-Block Adaptive Quantization)、DCT-BAQ (Discrete Cosine Transform) > WHT-BAQ (ffalsh-Hadamard Transform)。沿距离向，SAR原始数据的频谱能量主要集中在发射调频信号带宽之内；沿方位向，受天线方向图影响，多普勒谱能量在多普勒中心频率附近最高。因此，传统频域压缩算法的比特分配方式通常有两种第一，采用固定的比特分配方式，将距离向发射信号带宽之外的频谱置零；由系统参数计算多普勒中心频率，并将靠近多普勒中心频率的高能量区分配较多比特数，远离多普勒中心频率的低能量区分配较少比特数。第二，估算多普勒中心频率，沿方位向将多普勒中心频率附近的高能量区分成一块，并分配较多比特数；远离多普勒中心频率的低能量区分配较少比特数。但对于斜视SAR，多普勒中心不为零，而是与斜视角大小以及天线指向误差有关， 当对SAR方位向回波作FFT时，方位向频谱将因FFT的周期性而反折，多普勒中心位置不固定。因此在进行频域压缩时，若采用固定的比特分配方式，会造成压缩性能的严重恶化，必须对多普勒中心频率进行估计，但目前的多普勒中心估计算法运算量较大，星上无法实现。 此外，斜视SAR的多普勒中心与距离频率有关，二维频谱的能量沿方位向存在倾斜，若仅对方位向分块进行比特分配，则同一分块内的能量分布不均衡，会造成信噪比的下降。因此，如何选择一种实用化的自适应频域比特分配方法是解决问题的关键。
发明内容
本发明的目的是公开一种星载斜视SAR原始数据实用化频域压缩比特分配方法， 解决斜视SAR多普勒中心非零以及二维频谱的能量分布沿方位向存在倾斜的问题。为达到上述目的，本发明的技术解决方案是—种星载斜视SAR原始数据实用化频域压缩比特分配方法，其包括A)根据传输数据率、分块数以及距离向过采样系数确定比特分配方案；B)比特分配时，分别计算各分块数据幅度均值，并由大到小排序；C)无需大运算量的多普勒中心估计，由各块数据幅度均值大小区分高低能量区， 并进行比特分配。所述的原始数据实用化频域压缩比特分配方法，其所述A)步，包括Al)根据距离向过采样系数确定距离向频谱的弃置区，并将其置零；
A2)将有效频谱区域沿方位向和距离向分成MXN块，M和N均为大于0的整数，分块方式与BAQ压缩分块方式相同；A3)根据传输数据率、分块数以及距离向过采样系数计算出不同压缩比的分配比例。所述的原始数据实用化频域压缩比特分配方法，其所述似)步由于传统频域压缩区分高低能量区时采用的分块方式在斜视SAR的情况下失效， 针对斜视SAR 二维频谱具有多普勒中心非零以及能量分布沿方位向存在倾斜的特性，分别将有效频谱区域沿方位向和距离向分成MXN块，M和N均为大于0的整数，分块方式与BAQ 压缩分块方式相同，从而更好的区分高低能量区，提高算法性能。所述的原始数据实用化频域压缩比特分配方法，其所述Α； )步按照系统传输数据率的要求，每个采样点实部或虚部的平均编码比特数为R bits/sample ；距离向过采样率为α。s"；沿方位向距离向的分块数为MXN，M、N均为大于0 的整数；二维频谱沿距离向的能量主要分布在发射信号带宽之内，占距离向频谱总长度的比例为nr = l/aoSjrX100% (1)假设进行pbitBAQ压缩的数据块个数为A，进行qbitBAQ压缩的数据块个数为B， 且ρ > q，A、B均为大于等于0的整数，则比特分配方案的计算公式如下A+B = MXN(2)maxji ρ-Α + ^Β\ηΓ<λ(3)
^ Il A + B J /r J由公式⑵和(3)计算出A和B的值，即为比特分配方案。所述的原始数据实用化频域压缩比特分配方法，其所述C)步各分块数据幅度均值由大到小排序，前A个幅度均值较大的数据块为高能量区， 分配P bit;后B个幅度均值较小的数据块为低能量区，分配q bit。本发明的有益效果是，在满足传输速率要求的前提下，针对星载斜视SAR原始数据二维频谱的特性，采用了自适应比特分配的方式，既可以得到优于传统频域压缩算法的信噪比，又可以满足星上实现对运算量的要求。
图1斜视SAR原始数据频域压缩实用化比特分配流程图。图2斜视SAR原始数据二维频谱幅度仿真图。图3是本发明的斜视SAR原始数据二维频谱比特分配示意图。
具体实施例方式本发明的一种星载斜视SAR原始数据实用化频域压缩比特分配方法，由传输数据率、分块数以及距离向过采样系数确定比特分配方案；比特分配时，分别对斜视SAR 二维频谱进行距离向方位向分块，计算各块幅度均值大小并排序；按幅度均值大小以及比特分配方案对各块数据分配不同的比特数，无需进行多普勒中心估计。具体操作包括根据已知的传输数据率、分块数以及距离向过采样系数计算出不同压缩比特的分配比例；根据距离向过采样率，将斜视SAR回波的二维频谱在距离向发射带宽以外的部分置零；沿距离向方位向对频谱进行分块，计算各块幅度均值并按大小排序； 按幅度均值大小以及不同压缩比特的分配比例对各块数据进行BAQ压缩。其中幅度均值较大的为高能量区，分配较多比特；幅度均值较小的为低能量区，分配较少比特。下面结合附图和实施例对本发明进一步说明。图1为斜视SAR原始数据频域压缩实用化比特分配流程图。假设按照系统传输数据率的要求，每个采样点实部或虚部的平均编码比特数为R bits/sample ；距离向过采样率为0。3,,;沿方位向距离向的分块数为1父仏将斜视SAR原始数据作二维FFT变换到二维频域，则二维频谱沿距离向的能量主要分布在发射信号带宽之内，占距离向频谱总长度的比例为nr = l/aoSjrX100% (1)将距离向频谱的弃置区(发射信号带宽以外的部分)置零，并将发射信号带宽之内的频谱沿方位向和距离向平均分成MXN块。假设压缩比特分配方式为分配pbit的数据块个数为A，分配qbit的数据块个数为B，且ρ > q，则A+B = MXN(2)
权利要求
1.一种星载斜视SAR原始数据实用化频域压缩比特分配方法，其特征在于包括A)根据传输数据率、分块数以及距离向过采样系数确定比特分配方案；B)比特分配时，分别计算各分块数据幅度均值，并由大到小排序；C)无需大运算量的多普勒中心估计，由各块数据幅度均值大小区分高低能量区，并进行比特分配。
2.根据权利要求1所述的原始数据实用化频域压缩比特分配方法，其特征在于，所述 A)步，包括Al)根据距离向过采样系数确定距离向频谱的弃置区，并将其置零；A2)将有效频谱区域沿方位向和距离向分成MXN块，M和N均为大于0的整数，分块方式与BAQ压缩分块方式相同；A3)根据传输数据率、分块数以及距离向过采样系数计算出不同压缩比的分配比例。
3.根据权利要求2所述的原始数据实用化频域压缩比特分配方法，其特征在于，所述 A2)步由于传统频域压缩区分高低能量区时采用的分块方式在斜视SAR的情况下失效，针对斜视SAR 二维频谱具有多普勒中心非零以及能量分布沿方位向存在倾斜的特性，分别将有效频谱区域沿方位向和距离向分成MXN块，M和N均为大于0的整数，分块方式与BAQ压缩分块方式相同，从而更好的区分高低能量区，提高算法性能。
4.根据权利要求2所述的原始数据实用化频域压缩比特分配方法，其特征在于，所述 A3)步按照系统传输数据率的要求，每个采样点实部或虚部的平均编码比特数为Rbits/ sample ；距离向过采样率为α。_ ；沿方位向距离向的分块数为MXN，Μ、N均为大于0的整数；二维频谱沿距离向的能量主要分布在发射信号带宽之内，占距离向频谱总长度的比例为nr = l/aoSjrX100% (1)假设进行PbitBAQ压缩的数据块个数为A，进行qbitBAQ压缩的数据块个数为B，且ρ > q，A、B均为大于等于0的整数，则比特分配方案的计算公式如下
5.根据权利要求1所述的原始数据实用化频域压缩比特分配方法，其特征在于，所述 C)步各分块数据幅度均值由大到小排序，前A个幅度均值较大的数据块为高能量区，分配 Pbit ；后B个幅度均值较小的数据块为低能量区，分配qbit。
全文摘要
本发明公开了一种星载斜视SAR原始数据实用化频域压缩比特分配方法，涉及数据压缩技术，其根据传输数据率、分块数及距离向过采样系数确定比特分配方案；比特分配时，分别计算各分块数据幅度均值，并由大到小排序；无需大运算量的多普勒中心估计，由各块数据幅度均值大小区分高低能量区，并进行比特分配。本发明方法在满足传输数据率要求的前提下，比传统频域压缩比特分配方法具有更好的适应性和量化信噪比，解决了星载斜视合成孔径雷达(SAR)多普勒中心非零以及二维频谱的能量分布沿方位向存在倾斜的特性导致传统频域压缩算法失效的问题。
文档编号G01S13/90GK102298138SQ20101021842
公开日2011年12月28日 申请日期2010年6月25日 优先权日2010年6月25日
发明者华斌, 李信, 祁海明 申请人:中国科学院电子学研究所