[发明专利]一种SAR幅度图像均匀相对误差量化方法

摘要

本发明公开了一种SAR幅度图像均匀相对误差量化方法，该方法包括以下步骤：将SAR幅度图像中非零数据先变换到对数域；将变换后的对数域数据进行均匀量化，并将量化后的数据编码存储；对SAR幅度图像中零值数据进行单独量化编码，并存储；根据量化编码存储结果对原SAR幅度图像数据进行量化解码。本发明给出的均匀相对误差量化方法为采用了零值独立编码和对数变换后的均匀量化方法。该方法的相对量化误差符合均匀分布。通过使用实际SAR数据的实施例试验显示，均匀相对误差量化方法较均匀量化方法可使用更多的量化等级、具有更小的量化误差，且量化后图像具备更好的可视效果。

主权项

1.一种SAR幅度图像均匀相对误差量化方法，其特征在于，所述方法包括以下步骤：步骤A、将SAR幅度图像中非零数据先变换到对数域；步骤B、将变换后的对数域数据进行均匀量化，并将量化后的数据编码存储；步骤C、对SAR幅度图像中零值数据进行单独量化编码，并存储；步骤D、根据量化编码存储结果对原SAR幅度图像数据进行量化解码；所述步骤B中SAR幅度图像中不包含零值数据的量化编码；其SAR幅度图像中不包含零值数据的量化编码过程包括：令A表示SAR幅度图像中像素点的散射幅度值，将所有像素点的散射幅度值都变为dB值，即U＝10×log₁₀(A)   (1)其中log₁₀表示以10为底的对数，U为散射幅度值A变换后的对数数据；对U采用16bit的均匀量化，包括：对U进行归一化，即将其值域变为[0,1]，利用整幅图像全部像素点U值的最小值min(U)和最大值max(U)进行如下变换X＝(U‑min(U))/[max(U)‑min(U)]   (2)然后将X对应为[0,1,2,…,65535]量化等级，具体可使用下面的表达式实现Y＝round(65535×X)   (3)其中round表示四舍五入取整运算；步骤C中原SAR幅度图像中包含零值数据的量化编码，其SAR幅度图像中含零值数据的量化编码过程包括：令A表示SAR幅度图像中的非零值的散射幅度值，将非零散射幅度值变换为dB值，即U＝10×log₁₀(A)   (4)其中log₁₀表示以10为底的对数；对U采用16bit不包括零值的均匀量化，包括：对U进行归一化，即将其值域变为[0,1]，利用整幅图像全部非零像素点U值的最小值min(U)和最大值max(U)进行如下变换X＝(U‑min(U))/[max(U)‑min(U)]   (5)然后将X对应为[1,2,…，65535]量化等级，具体可使用下面的表达式实现Y＝round(65534×X)+1   (6)其中round(·)表示四舍五入运算取整运算；对图像中的z个零值像素点的Y值编码为16bit数的零值。