Измерьте пространственную частотную характеристику с помощью графика Imatest eSFR
sharpnessTable = measureSharpness(chart)sharpnessTable включает частоту для каждого ROI, в котором ответ спадает до 50% начальных и пиковых значений.
sharpnessTable = measureSharpness(chart,Name,Value)
[ также возвращает средний швейцарский франк вертикальных и горизонтальных ROI, с помощью входных параметров любого из предыдущих синтаксисов.sharpnessTable,aggregateSharpnessTable] = measureSharpness(___)
Считайте изображение графика eSFR в рабочую область.
I = imread('eSFRTestImage.jpg');Создайте esfrChart объект, затем отобразите график с аннотациями ROI. 60 наклонных ROI ребра помечены зелеными числами.
chart = esfrChart(I); displayChart(chart,'displayColorROIs',false,... 'displayGrayROIs',false,'displayRegistrationPoints',false)
Измерьте резкость края в ROI 25-28 и возвратите измерения в sharpnessTable. Включайте измерения MTF70 и MTF30 путем определения 'percentResponse' аргумент пары "имя-значение".
sharpnessTable = measureSharpness(chart,'ROIIndex',25:28,'PercentResponse',[70 30])
sharpnessTable=4×9 table
ROI slopeAngle confidenceFlag SFR comment MTF70 MTF70P MTF30 MTF30P
___ __________ ______________ ____________ ____________ ____________________________________________ ____________________________________________ ________________________________________ ________________________________________
25 4.2391 true {88x5 table} {0x0 double} 0.060224 0.058884 0.051757 0.058482 0.060224 0.058884 0.051757 0.058482 0.10732 0.11267 0.11185 0.11178 0.10732 0.11267 0.11185 0.11178
26 5.0254 true {88x5 table} {0x0 double} 0.18561 0.18597 0.18667 0.18597 0.18561 0.18597 0.18667 0.18597 0.26028 0.26236 0.26224 0.26241 0.26028 0.26236 0.26224 0.26241
27 4.7131 true {88x5 table} {0x0 double} 0.070042 0.069599 0.066239 0.06951 0.070042 0.069599 0.066239 0.06951 0.2169 0.21863 0.22086 0.21904 0.2169 0.21863 0.22086 0.21904
28 4.8996 true {88x5 table} {0x0 double} 0.19485 0.20418 0.19669 0.20168 0.19485 0.20418 0.19658 0.20168 0.26447 0.27531 0.26004 0.27188 0.26447 0.27531 0.25998 0.27188
Выберите четвертую строку в таблице резкости, которая соответствует ROI 28. Отобразите график швейцарского франка ROI.
idx = 4; plotSFR(sharpnessTable(idx,:))
Распечатайте MTF70 и измерения MTF30 ROI. Сравните измерения с графиком.
Измерение MTF70 красных и синих цветовых каналов немного меньше, чем 0,2, в то время как измерение MTF70 зеленых каналов и каналов яркости немного больше, чем 0,2. Эти измерения соглашаются с визуальным осмотром графика швейцарского франка, на котором значение швейцарского франка 0,7 происходит на пространственных частотах приблизительно 0,2 пары линии на пиксель.
mtf70 = sharpnessTable.MTF70(idx,:)
mtf70 = 1×4
0.1949 0.2042 0.1967 0.2017
Измерение MTF30 синего цветового канала заметно меньше, чем измерение MTF30 других цветовых каналов. Это измерение соглашается с визуальным осмотром графика швейцарского франка, на котором кривая швейцарского франка синего канала понижается более быстро, чем другие каналы.
mtf30 = sharpnessTable.MTF30(idx,:)
mtf30 = 1×4
0.2645 0.2753 0.2600 0.2719
Наклонные ребра на правильно ориентированном графике под углом 5 градусов горизонтали или вертикальные. Измерения резкости не точны, когда ориентация ребра значительно отклоняется от 5 градусов.
Резкость выше к центру изображенной области и уменьшений к периферии. Горизонтальная резкость обычно выше, чем вертикальная резкость.
Алгоритм измерения швейцарского франка на основе работы [2]
[3] Питера Бернса. Во-первых, measureSharpness определяет положение ребра с субпиксельным разрешением для каждой строки развертки, или строку или столбец пиксельного перпендикуляра к ребру, в ROI. Например, каждая строка пикселей является строкой развертки для почти вертикального ребра. Затем measureSharpness выравнивает и составляет в среднем строки развертки, чтобы создать сверхдискретизированный профиль интенсивности ребра. Функция берет производную профиля интенсивности и применяет функцию работы с окнами. Возвращенное измерение швейцарского франка является абсолютным значением преобразования Фурье оконной производной.
[1] Imatest. "Esfr". http://www.imatest.com/mathworks/esfr/.
[2] Записи, Питер. "Slanted-Edge MTF для Анализа Цифрового фотоаппарата и Сканера". Общество Обработки изображений Науки и техники; Продолжения Обработки изображений, Качества изображения, Системной Конференции по Захвату изображения. Портленд, Орегон, март 2000, стр 135–138.
[3] Записи, Питер. "sfrmat3: оценка швейцарского франка для цифровых фотоаппаратов и сканеров". URL: http://losburns.com/imaging/software/SFRedge/sfrmat3_post/index.html.