Вредоносное ПО (malware) - это назойливые или опасные программы,...
Исходя из полученных выше общих формул, можно указать конкретные способы определения координат центров тяжести тел.
1. Симметрия. Если однородное тело имеет плоскость, ось или центр симметрии (рис.7), то его центр тяжести лежит соответственно в плоскости симметрии, оси симметрии или в центре симметрии.
Рис.7
2. Разбиение. Тело разбивается на конечное число частей (рис.8), для каждой из которых положение центра тяжести и площадь известны.
Рис.8
3.Метод отрицательных площадей. Частный случай способа разбиения (рис.9). Он применяется к телам, имеющим вырезы, если центры тяжести тела без выреза и вырезанной части известны. Тело в виде пластинки с вырезом представляют комбинацией сплошной пластинки (без выреза) с площадью S 1 и площади вырезанной части S 2 .
Рис.9
4.Метод группировки. Является хорошим дополнением двух последних методов. После разбиения фигуры на составные элементы часть их бывает удобно объединить вновь, чтобы затем упростить решение путем учета симметрии этой группы.
Центры тяжести некоторых однородных тел.
1) Центр тяжести дуги окружности. Рассмотрим дугу АВ радиуса R с центральным углом . В силу симметрии центр тяжести этой дуги лежит на оси Ox (рис. 10).
Рис.10
Найдем координату по формуле . Для этого выделим на дуге АВ элемент ММ’ длиною , положение которого определяется углом . Координата х элемента ММ’ будет . Подставляя эти значения х и dl и имея в виду, что интеграл должен быть распространен на всю длину дуги, получим:
где L - длина дуги АВ , равная .
Отсюда окончательно находим, что центр тяжести дуги окружности лежит на ее оси симметрии на расстоянии от центра О , равном
где угол измеряется в радианах.
2) Центр тяжести площади треугольника. Рассмотрим треугольник, лежащий в плоскости Oxy , координаты вершин которого известны: A i (x i ,y i ), (i = 1,2,3). Разбивая треугольник на узкие полоски, параллельные стороне А 1 А 2 , придем к выводу, что центр тяжести треугольника должен принадлежать медиане А 3 М 3 (рис.11) .
Рис.11
Разбивая треугольник на полоски, параллельные стороне А 2 А 3 , можно убедиться, что он должен лежать на медиане А 1 М 1 . Таким образом, центр тяжести треугольника лежит в точке пересечения его медиан , которая, как известно, отделяет от каждой медианы третью часть, считая от соответствующей стороны.
В частности, для медианы А 1 М 1 получим, учитывая, что координаты точки М 1 - это среднее арифметическое координат вершин А 2 и А 3:
x c = x 1 + (2/3)∙(x М 1 - x 1) = x 1 + (2/3)∙[(x 2 + x 3)/2-x 1 ] = (x 1 + x 2 +x 3)/3.
Таким образом, координаты центра тяжести треугольника представляют собой среднее арифметическое из координат его вершин:
x c =(1/3)Σx i ; y c =(1/3)Σy i .
3) Центр тяжести площади кругового сектора. Рассмотрим сектор круга радиуса R с центральным углом 2α, расположенный симметрично относительно оси Ox (рис.12) .
Очевидно, что y c = 0, а расстояние от центра круга, из которого вырезан этот сектор, до его центра тяжести можно определить по формуле:
Рис.12
Проще всего этот интеграл вычислить, разбивая область интегрирования на элементарные секторы с углом d φ. С точностью до бесконечно малых первого порядка такой сектор можно заменить треугольником с основанием, равным R ×d φ и высотой R . Площадь такого треугольника dF =(1/2)R 2 ∙d φ, а его центр тяжести находится на расстоянии 2/3R от вершины, поэтому в (5) положим x = (2/3)R ∙cosφ. Подставляя в (5) F = αR 2 , получим:
С помощью последней формулы вычислим, в частности, расстояние до центра тяжести полукруга .
Подставляя в (2) α = π/2, получим: x c = (4R )/(3π) ≅ 0,4R .
Пример 1. Определим центр тяжести однородного тела, изображённого на рис. 13.
Рис.13
Тело однородное, состоящее из двух частей, имеющих симметричную форму. Координаты центров тяжести их:
Объёмы их:
Поэтому координаты центра тяжести тела
Пример 2. Найдем центр тяжести пластины, согнутой под прямым углом. Размеры – на чертеже (рис.14).
Рис.14
Координаты центров тяжести:
Площади:
|
Рис.15
В этой задаче удобнее разделить тело на две части: большой квадрат и квадратное отверстие. Только площадь отверстия надо считать отрицательной. Тогда координаты центра тяжести листа с отверстием:
координата 
так как тело имеет ось симметрии (диагональ).
Пример 4. Проволочная скобка (рис.16) состоит из трёх участков одинаковой длины l .
Рис.16
Координаты центров тяжести участков:
Поэтому координаты центра тяжести всей скобки:
Пример 5. Определить положение центра тяжести фермы, все стержни которой имеют одинаковую погонную плотность (рис.17).
Напомним, что в физике плотность тела ρ и его удельный вес g связаны соотношением: γ= ρg , где g - ускорение свободного падения. Чтобы найти массу такого однородного тела, нужно плотность умножить на его объем.
Рис.17
Термин «линейная» или «погонная» плотность означает, что для определения массы стержня фермы нужно погонную плотность умножить на длину этого стержня.
Для решения задачи можно воспользоваться методом разбиения. Представив заданную ферму в виде суммы 6 отдельных стержней, получим:
где L i длина i -го стержня фермы, а x i , y i - координаты его центра тяжести.
Решение этой задачи можно упростить, если сгруппировать 5 последних стержней фермы. Нетрудно видеть, что они образуют фигуру, имеющую центр симметрии, расположенный посредине четвертого стержня, где и находится центр тяжести этой группы стержней.
Таким образом, заданную ферму можно представить комбинацией всего двух групп стержней.
Первая группа состоит из первого стержня, для нее L 1 = 4 м, x 1 = 0 м, y 1 = 2 м. Вторая группа стержней состоит из пяти стержней, для нее L 2 = 20 м, x 2 = 3 м, y 2 = 2 м.
Координаты центра тяжести фермы находим по формуле:
x c = (L 1 ∙x 1 + L 2 ∙x 2)/(L 1 + L 2) = (4∙0 + 20∙3)/24 = 5/2 м;
y c = (L 1 ∙y 1 + L 2 ∙y 2)/(L 1 + L 2) = (4∙2 + 20∙2)/24 = 2 м.
Отметим, что центр С лежит на прямой, соединяющей С 1 и С 2 и делит отрезок С 1 С 2 в отношении: С 1 С /СС 2 = (x c - x 1)/(x 2 - x c ) = L 2 / L 1 = 2,5/0,5.
Вопросы для самопроверки
Что называется центром параллельных сил?
Как определяются координаты центра параллельных сил?
Как определить центр параллельных сил, равнодействующая которых равна нулю?
Каким свойством обладает центр параллельных сил?
По каким формулам вычисляются координаты центра параллельных сил?
Что называется центром тяжести тела?
Почему силы притяжения Земле, действующие на точку тела, можно принять за систему параллельных сил?
Запишите формулу для определения положения центра тяжести неоднородных и однородных тел, формулу для определения положения центра тяжести плоских сечений?
Запишите формулу для определения положения центра тяжести простых геометрических фигур: прямоугольника, треугольника, трапеции и половины круга?
Что называют статическим моментом площади?
Приведите пример тела, центр тяжести которого расположен вне тела.
Как используются свойства симметрии при определении центров тяжести тел?
В чем состоит сущность способа отрицательных весов?
Где расположен центр тяжести дуги окружности?
Каким графическим построением можно найти центр тяжести треугольника?
Запишите формулу, определяющую центр тяжести кругового сектора.
Используя формулы, определяющие центры тяжести треугольника и кругового сектора, выведите аналогичную формулу для кругового сегмента.
По каким формулам вычисляются координаты центров тяжести однородных тел, плоских фигур и линий?
Что называется статическим моментом площади плоской фигуры относительно оси, как он вычисляется и какую размерность имеет?
Как определить положение центра тяжести площади, если известно положение центров тяжести отдельных ее частей?
Какими вспомогательными теоремами пользуются при определении положения центра тяжести?
Алгоритм разработан группой экспертов в области фотографии (Joint Photographic Expert Group) специально для сжатия 24-битных и полутоновых изображений в 1991 году. Этот алгоритм не очень хорошо сжимает двухуровневые изображении, но он прекрасно обрабатывает изображения с непрерывными тонами, в которых близкие пикселы обычно имеют схожие цвета. Обычно глаз не в состоянии заметить какой-либо разницы при сжатии этим методом в 10 или 20 раз.
Алгоритм основан на ДКП, применяемом к матрице непересекающихся блоков изображения, размером 8х8 пикселей. ДКП раскладывает эти блоки по амплитудам некоторых частот. В результате, получается матрица, в которой многие коэффициенты, как правило, близки к нулю, которые можно представить в грубой числовой форме, т.е. в квантованном виде без существенной потери в качестве восстановления.
Рассмотрим работу алгоритма подробнее. Предположим, что сжимается полноцветное 24-битное изображение. В этом случае получаем следующие этапы работы.
Шаг 1. Переводим изображение из пространства RGB в пространство YCbCr с помощью следующего выражения:
Отметим сразу, что обратное преобразование легко получается путем умножения обратной матрицы на вектор , который по существу является пространством YUV:
.
Шаг 2. Разбиваем исходное изображение на матрицы 8х8. Формируем из каждой три рабочие матрицы ДКП – по 8 бит отдельно для каждой компоненты. При больших степенях сжатия блок 8х8 раскладывается на компоненты YCbCr в формате 4:2:0, т.е. компоненты для Cb и Cr берутся через точку по строкам и столбцам.
Шаг 3. Применение ДКП к блокам изображения 8х8 пикселей. Формально прямое ДКП для блока 8х8 можно записать в виде
где 
. Так как ДКП является «сердцем»
алгоритма JPEG, то желательно
на практике вычислять его как можно быстрее. Простым подходом для ускорения
вычислений является заблаговременное вычисление функций косинуса и сведения
результатов вычисления в таблицу. Мало того, учитывая ортогональность функций
косинусов с разными частотами, вышеприведенную формулу можно записать в виде
.
Здесь является матрицей, размером 8х8 элементов, описывающая 8-ми мерное пространство, для представления столбцов блока в этом пространстве. Матрица является транспонированной матрицей и делает то же самое, но для строк блока . В результате получается разделимое преобразование, которое в матричном виде записывается как
Здесь - результат ДКП, для вычисления которого требуется операций умножения и почти столько же сложений, что существенно меньше прямых вычислений по формуле выше. Например, для преобразования изображения размером 512х512 пикселей потребуется арифметических операций. Учитывая 3 яркостных компоненты, получаем значение 12 582 912 арифметических операций. Количество умножений и сложений можно еще больше сократить, если воспользоваться алгоритмом быстрого преобразования Фурье. В результате для преобразования одного блока 8х8 нужно будет сделать 54 умножений, 468 сложений и битовых сдвигов.
В результате ДКП получаем матрицу , в которой коэффициенты в левом верхнем углу соответствуют низкочастотной составляющей изображения, а в правом нижнем – высокочастотной.
Шаг 4. Квантование. На этом шаге происходит отбрасывание части информации. Здесь каждое число из матрицы делится на специальное число из «таблицы квантования», а результат округляется до ближайшего целого:
.
Причем для каждой матрицы Y, Cb и Cr можно задавать свои таблицы квантования. Стандарт JPEG даже допускает использование собственных таблиц квантования, которые, однако, необходимо будет передавать декодеру вместе со сжатыми данными, что увеличит общий размер файла. Понятно, что пользователю сложно самостоятельно подобрать 64 коэффициента, поэтому стандарт JPEG использует два подхода для матриц квантования. Первый заключается в том, что в стандарт JPEG включены две рекомендуемые таблицы квантования: одна для яркости, вторая для цветности. Эти таблицы представлены ниже. Второй подход заключается в синтезе (вычислении на лету) таблицы квантовании, зависящей от одного параметра , который задается пользователем. Сама таблица строится по формуле
На этапе квантования осуществляется управление степенью сжатия, и происходят самые большие потери. Понятно, что задавая таблицы квантования с большими коэффициентами, мы получим больше нулей и, следовательно, большую степень сжатия.
С квантованием связаны и специфические эффекты алгоритма. При больших значениях шага квантования потери могут быть настолько велики, что изображение распадется на квадраты однотонные 8х8. В свою очередь потери в высоких частотах могут проявиться в так называемом «эффекте Гиббса», когда вокруг контуров с резким переходом цвета образуется волнообразный «нимб».
Шаг 5. Переводим матрицу 8х8 в 64-элементный вектор при помощи «зигзаг»-сканирования (рис. 2).
Рис. 2. «Зигзаг»-сканирование
В результате в начале вектора, как правило, будут записываться ненулевые коэффициенты, а в конце образовываться цепочки из нулей.
Шаг 6. Преобразовываем вектор с помощью модифицированного алгоритма RLE, на выходе которого получаем пары типа (пропустить, число), где «пропустить» является счетчиком пропускаемых нулей, а «число» - значение, которое необходимо поставить в следующую ячейку. Например, вектор 1118 3 0 0 0 -2 0 0 0 0 1 … будет свернут в пары (0, 1118) (0,3) (3,-2) (4,1) … .
Следует отметить, что первое число преобразованной компоненты , по существу, равно средней яркости блока 8х8 и носит название DC-коэффициента. Аналогично для всех блоков изображения. Это обстоятельство наводит на мысль, что коэффициенты DC можно эффективно сжать, если запоминать не их абсолютные значения, а относительные в виде разности между DC коэффициентом текущего блока и DC коэффициентом предыдущего блока, а первый коэффициент запомнить так, как он есть. При этом упорядочение коэффициентов DC можно сделать, например, так (рис. 3). Остальные коэффициенты, которые называются AC-коэффициентами сохраняются без изменений.
Шаг 7. Свертываем получившиеся пары с помощью неравномерных кодов Хаффмана с фиксированной таблицей. Причем для DC и AC коэффициентов используются разные коды, т.е. разные таблицы с кодами Хаффмана.
Рис. 3. Схема упорядочения DC коэффициентов
Рис. 4. Структурная схема алгоритма JPEG
Процесс восстановления изображения в этом алгоритме полностью симметричен. Метод позволяет сжимать изображения в 10-15 раз без заметных визуальных потерь.
При разработке данного стандарта руководствовались тем, что данный алгоритм должен был сжимать изображения довольно быстро – не более минуты на среднем изображении. Это в 1991 году! А его аппаратная реализация должна быть относительно простой и дешевой. При этом алгоритм должен был быть симметричным по времени работы. Выполнение последнего требования сделало возможным появление цифровых фотоаппаратов, снимающие 24 битные изображения. Если бы алгоритм был несимметричен, было бы неприятно долго ждать, пока аппарат «перезарядится» - сожмет изображение.
Хотя алгоритм JPEG и является стандартом ISO, формат его файлов не был зафиксирован. Пользуясь этим, производители создают свои, несовместимые между собой форматы, и, следовательно, могут изменить алгоритм. Так, внутренние таблицы алгоритма, рекомендованные ISO, заменяются ими на свои собственные. Встречаются также варианты JPEG для специфических приложений.
После вычисления всех коэффициентов DCT их необходимо проквантовать. На этом шаге происходит отбрасывание части информации (небольшие потери происходят и на предыдущем шаге из-за конечной точности вычислений на компьютере). Каждое число из матриц коэффициентов DCT делится на специальное число из «таблицы квантования», а результат округляется до ближайшего целого. Как уже отмечалось, необходимо иметь три такие таблицы для каждой цветовой компоненты. Стандарт JPEG допускает использование четырех таблиц, и пользователь может выбрать любую из этих таблиц для квантования компонентов цвета. Все 64 числа из таблицы квантования являются параметрами JPEG. В принципе, пользователь может поменять любой коэффициент для достижения большей степени сжатия. На практике весьма сложно экспериментировать с таким большим числом параметров, поэтому программное обеспечение JPEG использует два подхода:
1. Таблица квантования, принятая по умолчанию. Две такие таблицы, одна для компоненты светимости (и для градации серого цвета), а другая для хроматических компонент, являются результатом продолжительного исследования со множеством экспериментов, проделанных комитетом JPEG. Они являются частью стандарта JPEG и воспроизведены в табл. 3.50. Видно, как коэффициенты QC таблиц растут при движении из левого верхнего угла в правый нижний угол. В этом отражается сокращение коэффициентов DCT, соответствующих высоким пространственным частотам.
2.
Вычисляется простая таблица коэффициентов квантования, зависящая от параметра , который задается
пользователем. Простые выражения типа 
гарантируют убывание коэффициентов из
левого верхнего угла в правый нижний.
Светимость
Если квантование сделано правильно, то в блоке коэффициентов DCT останется всего несколько ненулевых коэффициентов, которые будут сконцентрированы в левом верхнем углу матрицы. Эти числа являются выходом алгоритма JPEG, но их следует еще сжать перед записью в выходной файл. В литературе по JPEG это сжатие называется «энтропийным кодированием», детали которого будут разбираться в § 3.7.5. Три технических приема используется при энтропийном кодировании для сжатия целочисленных матриц 8x8.
3. 64 числа выстраиваются одно за другим как при сканировании зигзагом (см. рис. 3.5а). В начале стоят ненулевые числа, за которыми обычно следует длинный хвост из одних нулей. В файл выводятся только ненулевые числа (после надлежащего кодирования) за которыми следует специальный код ЕОВ (end-of-block, конец блока). Нет необходимости записывать весь хвост нулей (можно также сказать, что ЕОВ кодирует длинную серию нулей).
Пример
: В табл. 3.51 приведен список
гипотетических коэффициентов DCT, из которых только 4 не равны
нулю. При зигзагообразном упорядочении этих чисел получается последовательность
коэффициентов: 
Табл. 3.51. Квантованные коэффициенты.
А как написать подпрограмму для считывания элементов матрицы по зигзагу? Простейший способ состоит в ручном прослеживании этого пути и в записи результата в массив структур zz, в котором каждая структура состоит из пары координат клеток, через которые проходит зигзагообразный путь (см. рис. 3.52).
Если компоненты структуры zz обозначить zz.r и zz.с, то путь по зигзагу можно совершить с помощью следующего цикла
4. Ненулевые коэффициенты преобразования сжимаются по методу Хаффмана (см. § 3.7.5).
5. Первое из этих чисел (коэффициент DC, см. стр. 145) обрабатывается отдельно от других чисел (коэффициентов АС).
Рис. 3.52. Координаты зигзагообразного пути.
- Tutorial
UPD. Был вынужден убрать моноширинное форматирование. В один прекрасный день хабрапарсер перестал воспринимать форматирование внутри тегов pre и code. Весь текст превратился в кашу. Администрация хабра не смогла мне помочь. Теперь неровно, но хотя бы читабельно.
Вам когда-нибудь хотелось узнать как устроен jpg-файл? Сейчас разберемся! Прогревайте ваш любимый компилятор и hex-редактор, будем декодировать это:
Специально взял рисунок поменьше. Это знакомый, но сильно пережатый favicon Гугла:
Сразу предупреждаю, что описание упрощено, и приведенная информация не полная, но зато потом будет легко понять спецификацию.
Даже не зная, как происходит кодирование, мы уже можем кое-что извлечь из файла.
- маркер начала. Он всегда находится в начале всех jpg-файлов.
Следом идут байты
. Это маркер, означающий начало секции с комментарием. Следующие 2 байта
- длина секции (включая эти 2 байта). Значит в следующих двух
- сам комментарий. Это коды символов ":" и ")", т.е. обычного смайлика. Вы можете увидеть его в первой строке правой части hex-редактора.
Немного теории
Очень кратко по шагам:
Давайте подумаем, в каком порядке могут быть закодированы эти данные. Допустим, сначала полностью, для всего изображения, закодирован канал Y, затем Cb, потом Cr. Все помнят загрузку картинок на диал-апе. Если бы они кодировались именно так, нам бы пришлось ждать загрузки всего изображения, прежде чем оно появится на экране. Так же будет неприятно, если потерятся конец файла. Вероятно, существуют и другие весомые причины. Поэтому закодированные данные располагаются поочередно, небольшими частями.
Напоминаю, что каждый блок Y ij , Cb ij , Cr ij - это матрица коэффициентов ДКП, закодированная кодами Хаффмана. В файле они располагаются в таком порядке: Y 00 Y 10 Y 01 Y 11 Cb 00 Cr 00 Y 20
Чтение файла
После того, как мы извлекли комментарий, будет легко понять, что:- Файл поделен на секторы, предваряемые маркерами.
- Маркеры имеют длину 2 байта, причем первый байт .
- Почти все секторы хранят свою длину в следующих 2 байта после маркера.
FF D8 FF FE 00 04 3A 29 FF DB 00 43 00 A0 6E 78
FF FF FF FF FF FF FF FF FF FF FF FF FF FF DB 00
43 01 AA B4 B4 F0 D2 F0 FF FF FF FF FF FF FF FF
FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF
FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF
FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF
FF FF FF C0 00 11 08 00 10 00 10 03 01 22 00 02
11 01 03 11 01 FF C4 00 15 00 01 01 00 00 00 00
00 00 00 00 00 00 00 00 00 00 03 02 FF C4 00 1A
10 01 00 02 03 01 00 00 00 00 00 00 00 00 00 00
00 01 00 12 02 11 31 21 FF C4 00 15 01 01 01 00
00 00 00 00 00 00 00 00 00 00 00 00 00 00 01 FF
C4 00 16 11 01 01 01 00 00 00 00 00 00 00 00 00
00 00 00 00 11 00 01 FF DA 00 0C 03 01 00 02 11
03 11 00 3F 00 AE E7 61 F2 1B D5 22 85 5D 04 3C
82 C8 48 B1 DC BF FF D9
Маркер : DQT - таблица квантования.
FF DB
00 43
00 A0 6E 78
8C 78 64 A0 8C 82 8C B4 AA A0 BE F0 FF FF F0 DC
DC F0 FF FF FF FF FF FF FF FF FF FF FF FF FF FF
FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF FF
FF FF FF FF FF FF FF FF FF FF FF FF FF
Заголовок секции всегда занимает 3 байта. В нашем случае это . Заголовок состоит из:
Длина: 0x43 = 67 байт
Длина значений в таблице: 0 (0 - 1 байт, 1 - 2 байта)
[_0]
Идентификатор таблицы: 0
Оставшимися 64-мя байтами нужно заполнить таблицу 8x8.
Приглядитесь, в каком порядке заполнены значения таблицы. Этот порядок называется zigzag order: 
Маркер : SOF0 - Baseline DCT
Этот маркер называется SOF0, и означает, что изображение закодировано базовым методом. Он очень распространен. Но в интернете не менее популярен знакомый вам progressive-метод, когда сначала загружается изображение с низким разрешением, а потом и нормальная картинка. Это позволяет понять что там изображено, не дожидаясь полной загрузки. Спецификация определяет еще несколько, как мне кажется, не очень распространенных методов.FF C0
00 11
08 00 10 00 10 03 01 22 00 02
11 01 03 11 01
Длина: 17 байт.
Precision: 8 бит. В базовом методе всегда 8. Как я понял, это разрядность значений каналов.
Высота рисунка: 0x10 = 16
Ширина рисунка: 0x10 = 16
Количество компонентов: 3. Чаще всего это Y, Cb, Cr.
1-й компонент:
Идентификатор: 1
Горизонтальное прореживание (H 1): 2
[_2]
Вертикальное прореживание (V 1): 2
Идентификатор таблицы квантования: 0
2-й компонент:
Идентификатор: 2
Горизонтальное прореживание (H 2): 1
[_1]
Вертикальное прореживание (V 2): 1
3-й компонент:
Идентификатор: 3
Горизонтальное прореживание (H 3): 1
[_1]
Вертикальное прореживание (V 3): 1
Идентификатор таблицы квантования: 1
Теперь посмотрите, как определить насколько прорежено изображение. Находим H max =2 и V max =2 . Канал i будет прорежен в H max /H i раз по горизонтали и V max /V i раз по вертикали.
Маркер : DHT (таблица Хаффмана)
Эта секция хранит коды и значения полученные кодированием Хаффмана .FF C4
00 15
00 01 01 00 00 00 00
00 00 00 00 00 00 00 00 00 00 03 02
длина: 21 байт.
класс: 0 (0 - таблица DC коэффициэнтов, 1 - таблица AC коэффициэнтов).
[_0]
идентификатор таблицы: 0
Длина кода Хаффмана: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Количество кодов:
Количество кодов означает количество кодов такой длины. Обратите внимание, что секция хранит только длины кодов, а не сами коды. Мы должны найти коды сами. Итак, у нас есть один код длины 1 и один - длины 2. Итого 2 кода, больше кодов в этой таблице нет.
С каждым кодом сопоставлено значение, в файле они перечислены следом. Значения однобайтовые, поэтому читаем 2 байта.
- значение 1-го кода.
- значение 2-го кода.
Построение дерева кодов Хаффмана
Мы должны построить бинарное дерево по таблице, которую мы получили в секции DHT. А уже по этому дереву мы узнаем каждый код. Значения добавляем в том порядке, в каком указаны в таблице. Алгоритм прост: в каком бы узле мы ни находились, всегда пытаемся добавить значение в левую ветвь. А если она занята, то в правую. А если и там нет места, то возвращаемся на уровень выше, и пробуем оттуда. Остановиться нужно на уровне равном длине кода. Левым ветвям соответствует значение 0 , правым - 1 .Замечание:
Не нужно каждый раз начинать с вершины. Добавили значение - вернитесь на уровень выше. Правая ветвь существует? Если да, идите опять вверх. Если нет - создайте правую ветвь и перейдите туда. Затем, с этого места, начинайте поиск для добавления следующего значения.
Деревья для всех таблиц этого примера:
UPD (спасибо anarsoul):
В узлах первого дерева (DC, id =0) должны быть значения 0x03 и 0x02
В кружках - значения кодов, под кружками - сами коды (поясню, что мы получили их, пройдя путь от вершины до каждого узла). Именно такими кодами (этой и других таблиц) закодировано само содержимое рисунка.
Маркер : SOS (Start of Scan)
Байт в маркере означает - «ДА! Наконец-то то мы перешли непосредственно к разбору секции закодированного изображения!». Однако секция символично называется SOS.
FF DA
00 0C
03 01 00 02 11
03 11 00 3F 00
Длина заголовочной части (а не всей секции): 12 байт.
Количество компонентов сканирования. У нас 3, по одному на Y, Cb, Cr.
1-й компонент:
Номер компонента изображения: 1 (Y)
Идентификатор таблицы Хаффмана для DC коэффициэнтов: 0
[_0]
Идентификатор таблицы Хаффмана для AC коэффициэнтов: 0
2-й компонент:
Номер компонента изображения: 2 (Cb)
[_1]
3-й компонент:
Номер компонента изображения: 3 (Cr)
Идентификатор таблицы Хаффмана для DC коэффициэнтов: 1
[_1]
Идентификатор таблицы Хаффмана для AC коэффициэнтов: 1
Данные компоненты циклически чередуются.
На этом заголовочная часть заканчивается, отсюда и до конца (маркера ) закодированные данные.
0
Нахождение DC-коэффициента.
1.
Читаем последовательность битов (если встретим 2 байта , то это не маркер, а просто байт )
. После каждого бита сдвигаемся по дереву Хаффмана (с соответствующим идентификатором) по ветви 0 или 1, в зависимости от прочитанного бита. Останавливаемся, если оказались в конечном узле.
10
1011101110011101100001111100100
2.
Берем значение узла. Если оно равно 0, то коэффициент равен 0, записываем в таблицу и переходим к чтению других коэффициентов. В нашем случае - 02. Это значение - длина коэффициента в битах. Т. е. читаем следующие 2 бита, это и будет коэффициент.
10
10
11101110011101100001111100100
3. Если первая цифра значения в двоичном представлении - 1, то оставляем как есть: DC_coef = значение. Иначе преобразуем: DC_coef = значение-2 длина значения +1 . Записываем коэффициент в таблицу в начало зигзага - левый верхний угол.
Нахождение AC-коэффициентов.
1.
Аналогичен п. 1, нахождения DC коэффициента. Продолжаем читать последовательность:
10
10
1110
1110011101100001111100100
2.
Берем значение узла. Если оно равно 0, это означает, что оставшиеся значения матрицы нужно заполнить нулями. Дальше закодирована уже следующая матрица. Первые несколько дочитавших до этого места и написавших об этом мне в личку, получат плюс в карму. В нашем случае значение узла: 0x31.
Первый полубайт: 0x3 - именно столько нулей мы должны добавить в матрицу. Это 3 нулевых коэффициэнта.
Второй полубайт: 0x1 - длина коэффициэнта в битах. Читаем следующий бит.
10
10
1110
1
110011101100001111100100
3. Аналогичен п. 3 нахождения DC-коэффициента.
Как вы уже поняли, читать AC-коэффициенты нужно пока не наткнемся на нулевое значение кода, либо пока не заполнится матрица.
В нашем случае мы получим:
10
10
1110
1
1100
11
101
10
0
0
0
1
11110
0
100
и матрицу:
Вы заметили, что значения заполнены в том же зигзагообразном порядке?
Причина использования такого порядка простая - так как чем больше значения v и u, тем меньшей значимостью обладает коэффициент S vu в дискретно-косинусном преобразовании. Поэтому, при высоких степенях сжатия малозначащие коэффициенты обнуляют, тем самым уменьшая размер файла.
[-4 1 1 1 0 0 0 0] [ 5 -1 1 0 0 0 0 0]
[ 0 0 1 0 0 0 0 0] [-1 -2 -1 0 0 0 0 0]
[ 0 -1 0 0 0 0 0 0] [ 0 -1 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [-1 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[-4 2 2 1 0 0 0 0]
[-1 0 -1 0 0 0 0 0]
[-1 -1 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
Ой, я забыл сказать, что закодированные DC-коэффициенты - это не сами DC-коэффициенты, а их разности между коэффициентами предыдущей таблицы (того же канала)! Нужно поправить матрицы:
DC для 2-ой: 2 + (-4) = -2
DC для 3-ой: -2 + 5 = 3
DC для 4-ой: 3 + (-4) = -1
[-2 1 1 1 0 0 0 0] [ 3 -1 1 0 0 0 0 0] [-1 2 2 1 0 0 0 0]
………
Теперь порядок. Это правило действует до конца файла.
… и по матрице для Cb и Cr:
[-1 0 0 0 0 0 0 0]
[ 1 1 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
Так как тут только по одной матрице, DC-коэфициенты можно не трогать.
Вычисления
Квантование
Вы помните, что матрица проходит этап квантования? Элементы матрицы нужно почленно перемножить с элементами матрицы квантования. Осталось выбрать нужную. Сначала мы просканировали первый компонент, его компонента изображения = 1. Компонент изображения с таким идентификатором использует матрицу квантования 0 (у нас она первая из двух). Итак, после перемножения:
[ 0 120 280 0 0 0 0 0]
[ 0 -130 -160 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
Аналогично получаем еще 3 матрицы Y-канала…
[-320 110 100 160 0 0 0 0] [ 480 -110 100 0 0 0 0 0]
[ 0 0 140 0 0 0 0 0] [-120 -240 -140 0 0 0 0 0]
[ 0 -130 0 0 0 0 0 0] [ 0 -130 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [-140 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[-160 220 200 160 0 0 0 0]
[-120 0 -140 0 0 0 0 0]
[-140 -130 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
… и по матрице для Cb и Cr.
[-170 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 180 210 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
[ 0 0 0 0 0 0 0 0] [ 0 0 0 0 0 0 0 0]
Обратное дискретно-косинусное преобразование
Формула не должна доставить сложностей*. S vu - наша полученная матрица коэффициентов. u - столбец, v - строка. s yx - непосредственно значения каналов.
*Вообще говоря, это не совсем правда. Когда я смог декодировать и отобразить на экране рисунок 16x16, я взял изображение размером 600x600 (кстати, это была обложка любимого альбома Mind.In.A.Box - Lost Alone). Получилось не сразу - всплыли различные баги. Вскоре я мог любоваться корректно загруженной картинкой. Только очень огорчала скорость загрузки. До сих пор помню, она занимала 7 секунд. Но это и неудивительно, если бездумно пользоваться приведенной формулой, то для вычисления одного канала одного пикселя потребуется нахождения 128 косинусов, 768 умножений, и сколько-то там сложений. Только вдумайтесь - почти тысяча непростых операций только на один канал одного пиксела! К счастью, тут есть простор для отимизации (после долгих экспериментов уменьшил время загрузки до предела точности таймера 15мс, и после этого сменил изображение на фотографию в 25 раз большей площадью. Возможно, напишу об этом отдельной статьей).
Напишу результат вычисления только первой матрицы канала Y (значения округлены):
[ 87 72 50 36 37 55 79 95]
[-10 5 31 56 71 73 68 62]
[-87 -50 6 56 79 72 48 29]
И 2-х оставшихся:
Cb Cr
[ 60 52 38 20 0 -18 -32 -40] [ 19 27 41 60 80 99 113 120]
[ 48 41 29 13 -3 -19 -31 -37] [ 0 6 18 34 51 66 78 85]
[ 25 20 12 2 -9 -19 -27 -32] [-27 -22 -14 -4 7 17 25 30]
[ -4 -6 -9 -13 -17 -20 -23 -25] [-43 -41 -38 -34 -30 -27 -24 -22]
[ -37 -35 -33 -29 -25 -21 -18 -17] [-35 -36 -39 -43 -47 -51 -53 -55]
[ -67 -63 -55 -44 -33 -22 -14 -10] [ -5 -9 -17 -28 -39 -50 -58 -62]
[ -90 -84 -71 -56 -39 -23 -11 -4] [ 32 26 14 -1 -18 -34 -46 -53]
[-102 -95 -81 -62 -42 -23 -9 -1] [ 58 50 36 18 -2 -20 -34 -42]
- О, пойду-ка поем!
- Да я вообще не въезжаю, о чем речь.
- Раз значение цветов YCbCr получены, осталось преобразовать в RGB, типа так: YCbCrToRGB(Y ij , Cb ij , Cr ij) , Y ij , Cb ij , Cr ij - наши полученные матрицы.
- 4 матрицы Y, и по одной Cb и Cr, так как мы прореживали каналы и 4 пикселям Y соответствует по одному Cb и Cr. Поэтому вычислять так: YCbCrToRGB(Y ij , Cb , Cr )
YCbCr в RGB
R = Y + 1.402 * CrG = Y - 0.34414 * Cb - 0.71414 * Cr
B = Y + 1.772 * Cb
Не забудьте прибавить по 128. Если значения выйдут за пределы интервала , то присвоить граничные значения. Формула простая, но тоже отжирает долю процессорного времени.
Вот полученные таблицы для каналов R, G, B для левого верхнего квадрата 8x8 нашего примера:
255 248 194 148 169 215 255 255
255 238 172 115 130 178 255 255
255 208 127 59 64 112 208 255
255 223 143 74 77 120 211 255
237 192 133 83 85 118 184 222
177 161 146 132 145 162 201 217
56 73 101 126 144 147 147 141
0 17 76 126 153 146 127 108
231 185 117 72 67 113 171 217
229 175 95 39 28 76 139 189
254 192 100 31 15 63 131 185
255 207 115 46 28 71 134 185
255 241 175 125 112 145 193 230
226 210 187 173 172 189 209 225
149 166 191 216 229 232 225 220
72 110 166 216 238 231 206 186
255 255 249 203 178 224 255 255
255 255 226 170 140 187 224 255
255 255 192 123 91 138 184 238
255 255 208 139 103 146 188 239
255 255 202 152 128 161 194 232
255 244 215 200 188 205 210 227
108 125 148 172 182 184 172 167
31 69 122 172 191 183 153 134
