Отрегулируйте гистограмму N-D изображения в соответствии с гистограммой опорного изображения
B = imhistmatchn(A,ref)A и возвращает выходное изображение B гистограмма которого приблизительно совпадает с гистограммой опорного изображения ref. Оба A и ref должны быть изображениями в оттенках серого, но они не должны иметь одинаковый тип данных, размер или количество измерений.
B = imhistmatchn(A,ref,nbins)nbins равномерно разнесенные ячейки в пределах соответствующего диапазона для данного типа данных изображения. Возвращенное изображение B имеет не более nbins дискретные уровни.
Если тип данных изображения:
single или double, диапазон гистограмм - [0, 1].
uint8, диапазон гистограмм - [0, 255].
uint16, диапазон гистограмм - [0, 65535].
int16, диапазон гистограмм - [-32768, 32767].
Загрузите N-D изображение в оттенках серого в рабочую область. Также загрузите изображение в градациях серого, чтобы получить эталонную гистограмму.
load mri D load mristack
Отображение исходного объема в виде фрагментов.
figure montage(D,'DisplayRange',[]) title('Original 3-D Image')
Изменение формы ссылки в виде стека фрагментов оттенков серого для отображения.
ref = reshape(mristack,[256,256,1,21]);
Отображение опорного объема в виде фрагментов. Для правильного отображения на экране эталонный объем уменьшается в 0,5 раза с помощью imresize.
ref_downsized = imresize(ref,0.5); figure montage(ref_downsized,'DisplayRange',[]) title('Reference 3-D Image')
Сопоставить гистограмму D к гистограмме полного размера ref.
Dmatched = imhistmatchn(D,ref);
Просмотрите выходные данные. Обратите внимание, что уровни яркости выходных данных более точно соответствуют опорному изображению, чем исходному изображению.
figure montage(Dmatched,'DisplayRange',[]) title('Histogram Matched MRI')
histeq | imadjust | imhist | imhistmatch