В результате использования данного алгоритма формируется поток битов, который подается в линию связи. При этом распределение имеющейся квоты битов происходит по стандартам алгоритма SPIHT. Очевидно, что в случае зашумленных данных, алгоритм SPIHT на финальной стадии использует нулевую зону, что равносильно грубой пороговой обработке вейвлет-коэффициентов. Таким образом, вероятны появления артефактов звона в восстановленном РСАизображении, снижающие эффект от предварительной перестановки ПОД, из-за ошибок квантования, появляющихся на стадии кодирования. Следовательно, необходимо рассмотреть вопрос о снижении ошибок квантования при SPIHTкодировании опорного РСА-изображения. 3.5. Методика расчета квантователя. Распределение квоты бит, выполняемое в рамках SPIHT-кодирования, происходит с учетом принятия гипотезы о равномерном квантовании с высоким разрешением. Равномерное квантование не всегда обеспечивает наименьшую ошибку, особенно, при низких скоростях кодирования для вейвлет-коэффициентов с малой амплитудой. Таким образом, необходимо учитывать особенности распределения вейвлет-коэффициентов в каждой субполосе, то есть учитывать их плотность вероятности при расчете интервала квантования. Одним из вариантов является применение многоуровневого квантователя Ллойда-Макса [112, 113]. Пусть с помощью какого-либо из предложенных выше алгоритмов получено распределение квоты бит: ./ R, /=1 (3.8) где оптимальные/ 1,..., J Тогда число уровней квантователя дляу-й субполосы равно: Ll round j = I--.J (3.9) Пусть вычислены также оценки вейвлет-коэффициентов без учета квантования по уровню Д/1 дпп’{и +<р(^;\тш)? j h-, J ? то есть после пороговой обработки с 144 и,[./] |
-203 3.6. Методы и алгоритмы расчета квантователя Распределение квоты бит, выполненное в предыдущих пунктах данной главы для разных функций пороговой обработки вейвлет-коэффициентов, происходит с учетом принятия гипотезы о равномерном квантовании с высоким разрешением. Равномерное квантование не всегда обеспечивает наименьшую ошибку, особенно, при низких скоростях кодирования для вейвлеткоэффициентов с малой амплитудой, но которые оказались выше порога. Таким образом, необходимо учитывать особенности распределения вейвлеткоэффициентов в каждой субполосе, то есть учитывать их плотность вероятности при расчете интервала квантования. Одним из вариантов является применение многоуровневого квантователя Ллойда-Макса/155, 156/. Пусть с помощью какого-либо из предложенных выше алгоритмов получено распределение квоты бит: Rc = R(B) = £ а fi] , (3.6.1) 7=1 где Ъ* оптимальные значения бюджета бит для каждой субполосы, j = 1,...,./. Тогда число уровней квантователя дляj-й субполосы равно: Lj =round(2“Л ), j = (3.6.2) Пусть вычислены также оценки вейвлет-коэффициентов без учета квантования* по уровню Wy = Wy +9(Wy^,T ), j = , то есть после пороговой^ ^ f 9 обработки с найденной величиной порога. Тогда в пределах каждой субполосы имеется собственное распределение оценок значимых вейвлет-коэффициентов И ) для построения квантователя применяется плотность вероятности Pw (wY) 5 вычисляемая по зашумленным вейвлет-коэффициентам, использование оценок wY, позволяет рассчитать оптимальный квантователь Ллойда-Макса с учетом изменений динамического диапазона вейвлет-коэффициентов, происходящих |