Методы сложения поправочных итераций при калибровке

теоретические и практические аспекты колориметрии, системы управления цветом
Ответить

Какой вариант сложения итераций при калибровке вы используете?

Можно выбрать до 1 варианта ответа

 
 
Результаты голосования

Аватара пользователя
mihas
Администратор
Сообщения: 1368
Зарегистрирован: 18 авг 2004, 16:58
Откуда: Москва
Контактная информация:

Методы сложения поправочных итераций при калибровке

Сообщение mihas »

Итерации или повторное применение математических операций - это последовательное шаг за шагом приближение к идеалу в калибровке. Итерация отличается от повторения тем, что функция итерации не вызывает заново сама себя. Откат к линейным пластинам всякий раз при начале калибровки - это не итерации, это повторения одного и того же. Итерации используют тот же математический инструментарий, но всякий раз с новыми данными, возникшими на предыдущем этапе работы цикличной функции.
Когда мы откалибровали процесс (посчитаем на примере офсетного процесса) но по результату увидели, что идеал не достигнут, хотя мы к нему и приблизились, мы можем продолжить калибровку, но не повторить ее с нуля, а перейти к новому шагу, оперируя данными с предыдущей итерации. Таких итераций может быть бессчетное количество, офсет постоянно немного "дышит", поэтому подновлять калибровочную, не откатываясь далеко назад в самое начало, - обычное рутинное дело. Такие калибровки я выполняю для нескольких печатных машин приблизительно раз в месяц на протяжении ~10 лет, иногда они носят чисто косметический характер и почти не меняют размер точки на пластине, иногда поправки к точке ощутимые. Так например смена времен года с жары на похолодание и наоборот влечет как правило за собой заметные поправки. Смена краски (марки или партии) может также повлечь за собой необходимость существенно подновить калибровку новой итерацией. Когда через несколько лет в документе Excel перестает хватать столбцов на все итерации - я переношу последнюю результирующую колонку в новый документ (в виде значений, а не формул) и продолжаю прибавлять итерации далее уже в новом документе.

Очевидно, что суммировать отклонения в нелинейных величинах приращения к точке - TVI или DotGain - нельзя. А вот оперировать суммами линейных величин размеров точки в файле или на пластине можно смело. Для перехода от нелинейного приращения точки в печати к поправке линейной величины точки на пластине есть семейство удобных калькуляторов: для приведения кривой TVI под ISO 12647-2:2004 и ГОСТ Р 54766-2011, для кривой TVI под ISO 12647-2:2013, для компенсации сразу по всем секциям к любому из указанных стандартов (интегрирован также со спектральным калькулятором и может на выходе помимо таблиц записывать device link profile в формате icc).

Проиллюстрируем графически. На первом рисунке представлена измеренная и посчитанная градационная характеристика офсетной печати синей кривой, целевая кривая А стандарта красной кривой, и поправка к размеру точки на пластине зеленой кривой:
Измеренное в печати приращение к точке TVI и калибровочная кривая (зеленым) для приближения к идеалу стандарта.
Измеренное в печати приращение к точке TVI и калибровочная кривая (зеленым) для приближения к идеалу стандарта.
• 42.89 КБ • 3650 просмотров
На втором рисунке показано, как изменилась печать после применения к пластинам поправочной (калибровочной) зеленой кривой верхнего рисунка. Печать заметно приблизилась к идеалу - красной целевой кривой стандарта.
Последовательное приближение к идеалу. Измеренное в печати приращение к точке TVI и калибровочная второй итерации.
Последовательное приближение к идеалу. Измеренное в печати приращение к точке TVI и калибровочная второй итерации.
• 39.42 КБ • 3650 просмотров
Заметим, что печать, показанная синей кривой на втором графике, выполнена не с линейных пластин, а с пластин со внесенной поправкой к размеру точки на пластине в разделе calibration рипа. Но отсчет приращения синим цветом кривой по вертикальной оси выполнен по всегда неизменной шкале файла, это важно понимать, и мы чуть позже снова к шкалам градуировки обратимся. Нам осталось лишь прибавить новую поправку (уже почти линейную зеленую кривую второго графика) к предыдущей зеленой кривой первого графика, или сложить столбцы D+H из таблицы ниже по тексту. И простые операции сложения линейных величин - поправок к размеру точки на пластине - нам в помощь. В данном сложении поправок используется тот же принцип, что во всех рипах, когда линеаризационная кривая складывается с калибровочной кривой для лостижения заданной точки на пластине.
В таблице представлены наглядно такие вычисления, такую таблицу можно дальше продолжать вправо хоть до бесконечности, внося все новые и новые итерации и не откатываясь назад к первым колонкам AB.
•••• A •••••••• B •••••••• C •••••••• D •••••••• E •••••••• F •••••••• G •••••••• H •••••••• I •••••••• J ••••
TVTVI 1TVDot to A (1)TVTVI 2TVDot to A (2)TVResult Dot+
10001000100010001000
981.29980.42981.35980.3980.72
953951953.34950.3951.3
905.53901.6906.1900.57902.17
808.94802.698010.01801.25803.94
7010.94703.127012.33701.55704.67
6011.56603.286013.78601.09604.37
5011.04503.535013.89500.68504.21
4013.17400.14012.93400.29400.39
3012.9230-1.653011.1530-0.1630-1.81
2010.0120-1.64207.32200.3920-1.25
106.5810-1.89104.910-0.5410-2.43
54.875-1.7552.145-0.035-1.78
22.792-0.8521.142-0.22-1.05
0000000000
Прокомментируем содержимое колонок (имена колонок отмечены знаком ••••):
AB - процент заполнения площади растровой точкой, градуировка шкалы или Tone Value (TV сокращенно) и приращение полутонов в печати первой итерации Tone Value Increase или TVI, колориметрически нелинейно рассчитанное. Данные представлены на первом рисунке синей кривой.
CD - рассчитанная калибровочная кривая или поправка к размеру точки на пластине. Расчет выполнен в упомянутом калькуляторе калибровки с применением интерполяции сплайнами Эрмита (Hermite) и с выбранной в поле "Целевая" целевой кривой A(14) стандарта ISO и ГОСТ. Данные столбцов представлены на первом рисунке зеленой кривой. Красной кривой представлена на рисунке целевая кривая А(14) из стандарта (в таблице ее нет, но точные табличные значения могут быть получены тут).
EF - неизменные TV контрольной шкалы в файле в колонке E и приращение TVI в результате печати с применением калибровочной кривой CD.
GH - вторая итерация поправки к размеру точки на пластине, представленная зеленой кривой на втором рисунке.
IJ - сумма двух итераций, в колонке J просто линейно суммируются колонки линейных значений поправке к точке первой и второй итерации D и H, итого J = D + H.

В принципе можно воспользоваться альтернативным сложением итераций, оно подключает к расчетам виртуальную нелинейную пластину по оси градуировки, представим эту математику в новой таблице:
•••• A •••••••• B •••••••• C •••••••• D •••••••• E •••••••• F •••••••• G •••••••• H •••••••• I •••••••• J •••••••• K •••••••• L ••••
TVTVI 1TVDot to A (1)TVTVI 2TVDot to A (2)Summ C+DSumm E+FTVResult Dot+
10001000100010001001001000
981.29980.42981.35980.398.4299.35980.66
953951953.34950.39698.34951.25
905.53901.6906.1900.5791.696.1902.1
808.94802.698010.01801.2582.6990.01803.84
7010.94703.127012.33701.5573.1282.33704.63
6011.56603.286013.78601.0963.2873.78604.35
5011.04503.535013.89500.6853.5363.89504.29
4013.17400.14012.93400.2940.152.93400.47
3012.9230-1.653011.1530-0.1628.3541.1530-1.82
2010.0120-1.64207.32200.3918.3627.3220-1.25
106.5810-1.89104.910-0.548.1114.910-2.43
54.875-1.7552.145-0.033.257.145-1.77
22.792-0.8521.142-0.21.153.142-0.97
000000000000
ABCDEFGH - повтор колонок первой таблицы.
I - суммируются колонки C и D - получаем значения градуировки на пластине после первой калибровки, а не линейный TV файла шкалы.
J - суммируются колонки E и F - градуировки шкалы и приращения этих значений во второй печати, получаем абсолютное значение TVI (TV+TVI) когда величина приращения суммируется со значением точки в файле шкалы. Специально отмечаем ту особенность, что в столбце J данные приращения представлены по оси градуировки TV файла, тогда как столбец градуировки I содержит данные TV пластины. То есть приращения, посчитанные относительно файла шкалы, помещаются не на шкалу градуировки файла, а на шкалу градуировки пластин. Тем не менее, линеаризационная кривая в альтернативном сложении не учитывается, поэтому TV пластин колонки I в расчетах все равно носят виртуальный абстрактный характер, реальная сумма кривых линеаризации и калибровки в рипе не совпадает с физически измеренной геометрией площади точки на пластине.
KL - вычисленная в калибровочном калькуляторе поправка к размеру точки из колонок IJ на входе калькулятора. К удобным заданным целочисленным TV оси в колонке K от произвольных TV пластины колонки I нас возвращает интерполяция для любой оси в калибровочном калькуляторе (среднее окно ввода формата данных).
Интерфейс программы калибровки и форматирование данных на выходе по среднему окну
Интерфейс программы калибровки и форматирование данных на выходе по среднему окну
• 113.25 КБ • 3650 просмотров
В обеих таблицах столбцы D,H,L посчитаны в калибровочном калькуляторе cielab.xyz/dgcor/ на JavaScript, и я не могу привести математику этих столбцов в Excel к сожалению: в расчетах используются вложенные друг в друга циклы и алгоритмы не простой кусочно-линейной, а сложной сплайновой интерполяции, и то и другое просто как функцию ячейки не пропишешь в этом прекрасном табличном редакторе, понадобится как минимум VB и большущий кусок кода на нем.

Если мы сравним в двух таблицах результирующую последнюю колонку калибровки по двум итерациям - мы не увидим большой разницы, максимальные различия в 0.08% размера точки. Это многократно меньше, чем минимальная дискретность 8-битного вывода с шагом 0.39, то есть на оттиске мы не увидим разницу всего в 0.08% точки никогда. В рип обычно вносят поправки с точностью до десятых, а не до сотых, но и тут разница в 0.1% в 4 раза меньше дискретности 8-битного вывода, поэтому разницу в одну десятую процента точки мы увидим вряд ли.
У первого и второго вариантов сложения итераций есть свои поклонники, каждая из них содержит свои плюсы и минусы. Первый вариант сложения менее громоздкий и прост в администрировании, второй вариант больше подвержен осцилляциям в критично-сложных случаях при использовании методов интерполяции, критичных к осцилляциям, таких как кубическая интерполяция (интерполяция кубическими сплайнами), интерполяция многочленом Лагранжа и подобными. Сплайны Эрмита и Акимы менее подвержены этому недостатку - осцилляциям - незапланированным выбросам на кривой. И во всех случаях каких-то запредельно нереалистичных данных мы можем "словить" ту или иную ошибку и в первом и во втором способе сложения итерационных поправок.

Полезный совет: во избежание ошибок как при расчете результирующей калибровочной, так и во избежание ошибок рипа (многочлены Лагранжа там запросто могут нам осцилляции выдать на пластине при сложении калибровочной кривой с линеаризационной): шкала от 0 до 100 просто с шагом 10 подойдет для калибровки, но можно и уточнить также как в некоторых рипах значения точек 2, 5, 95, 98. Попытки зачастить на шкале TV в светах вроде 1-2-3-4 ни к чему хорошему не приведут. Такие данные можно промерить и использовать в расчетах, но такие данные не нужно "в лоб" прописывать в рип. Вопроса выбора оптимальных шкал для печати с разным числом патчей и их расположением, вопрос усреднения TVI таких шкал с использованием интерполяции разбирался по ссылке, точность калибровки заметно повышается.

Результатом и первого и второго варианта сложения калибровочных итераций будет зеленая калибровочная кривая на рисунке, ее то мы и внесем в рип в его раздел calibration в соответствии с двумя последними результирующими столбцами предложенных таблиц:
Две итерации калибровки, сложенные вместе в зеленой кривой.
Две итерации калибровки, сложенные вместе в зеленой кривой.
• 44.78 КБ • 3650 просмотров
В зависимости от формы кривой TVI мы можем в расчетах использовать как интерполяцию функции TVI, так и ее различные аппроксимации (10 вариантов сглаживания или приближения в калькуляторе), управляя соответственно результатом калибровки. Тут также может пригодиться альтернативный алгоритм представления данных по шкале TV пластины на второй итерации калибровки:
Интерполяция и аппроксимация функции TVI печати
Интерполяция и аппроксимация функции TVI печати
• 81.26 КБ • 3650 просмотров
В первом варианте сложения итераций мы всегда остаемся в рамках градационной шкалы файла. Во втором варианте по градационной шкале у нас появляются нелинейные точки на пластине по оси X, тогда как соотносятся они все же с приращением относительно линейных полутонов файла по оси Y. Поскольку в большинстве случаев практического применения того или иного сложения итерациями линейных величин поправки к размеру точки на пластине большой разницы по результату между алгоритмами сложения нет - вы можете выбрать любой из предложенных вариантов, который кажется вам идейно ближе. Переходить в расчетах на абстрактную градационную шкалу пластины примерно с тем же результатом (обычная максимальная погрешность в 0.08% при использовании сплайнов Эрмита для интеполяции), или суммировать линейные величины (относительные поправки к точке) - решать вам. Абстрактность и виртуальность градационной шкалы пластины альтернативного метода состоит в том, что не учитывает реальное сложение кривой калибровочной с кривой линеаризационной, значения полутонов пластины не соответствуют реально измеренным, так как сторонники альтернативного метода принципиально отвергают возможность суммирования линейных величин поправки к размеру точки, то есть теоретически отвергают ровно то, что все равно практически делает любой рип с их данными, суммируя кривые из разделов линеаризейшен и калибрейшен.

В заключение приведем разницу между простым расчетом сложения поправок и его чуть усложненной альтернативой в таблице и на графике. Надо также отметить что третья, четвертая, пятая и так далее итерации внесения поправок к размеру точки будут различаться между представленными способами исчисления итераций по-прежнему на ту же небольшую величину - не более нескольких сотых размера точки.
•••• A •••••••• B ••••
TVdifference
1000
980.06
950.05
900.07
800.1
700.04
600.02
50-0.08
40-0.08
300.01
200
100
5-0.01
2-0.08
00
Разница в поправке к размеру точки между вариантами сложения на графике
Разница в поправке к размеру точки между вариантами сложения на графике
• 19.02 КБ • 3650 просмотров
Не забывайте, что отдельные вычисления по красочным секциям мы проводим в одном калькуляторе, но есть и другой, где все расчеты производятся быстрее в 4 раза по 4 секциям разом. И он интегрирован со спектральным калькулятором, что удобно: спектральный предоставляет данные TVI, посчитанные из спектральных замеров, в нужном формате, напрямую передавая их между страницами с калькуляторами.

Короткие инструкции по калибровке на Youtube, Lesson 1 and 2.
Вложения
Аватара пользователя
mihas
Администратор
Сообщения: 1368
Зарегистрирован: 18 авг 2004, 16:58
Откуда: Москва
Контактная информация:

Re: Методы сложения поправочных итераций при калибровке

Сообщение mihas »

Немного конкретики про ввод калибровочных в различные рипы приводится в этой новой теме.
В этой теме подобраны все необходимы ссылки по процессу калибровки.
Ответить

Вернуться в «Колориметрия - наука о цвете: теория и практика»