Как насчет этого?
Сделайте некоторые данные:
>> im = [rand(200,200),rand(200,200)*2];
Отобразите его и выберите ROI, используя imfreehand
:
>> imagesc(im)
>> h = imfreehand;
Сделать маску из ROI:
>> mask = createMask(h);
Просто чтобы показать, что произошло, выделите область за пределами ROI NaN и отобразите:
>> im(~mask) = NaN;
>> figure; imagesc(im)
Рассчитать среднее значение пикселей в области интереса (которое больше 0,5, поскольку я выбрал область с одной стороны от середины):
>> mean(im(mask))
ans =
0.83782
Надеюсь, это поможет!