повышение контрастности изображения с помощью окрестности

Привет, я хочу повысить контрастность изображения, используя значения соседних пикселей.

Пусть изображение рассматривается как u0. Затем я хочу улучшить изображение, используя формулу

Здесь M1 — минимумы, а M2 — максимумы u0 среди соседних пикселей, Mg — максимальное значение уровня серого исходного изображения. Район, взятый для моей операции, — 9х9. uN — новое сгенерированное изображение (изображение с контрастным усилением).

Я пробовал следующий код, но не уверен, прав ли я или нет.

%Generate a contrast enhanced image
tic
clear all; close all;
I = imread('4.jpg');
I = imresize(I,[128 128]);
if size(I,3)== 3
            P = rgb2gray(uint8(I));
            P = double(P);
elseif size(I,3) == 2
            P = 0.5.*(double(I(:,:,1))+double(I(:,:,2)));
else
            P = double(I);
end
ssize=9;
mg=max(P(:));

f1 = @(x) min(x(:));
m1 = nlfilter(P,[9 9],f1);
f2 = @(x) max(x(:));
m2 = nlfilter(P,[9 9],f2);

P_op=((P-m1)./(m2-m1)).*mg;
subplot(2,1,1),imagesc(P,[0 255]);colormap(gray);axis off;
subplot(2,1,2),imagesc(P_op,[0 255]);colormap(gray);axis off;
toc

Некоторые из результатов, которые я получил, показаны ниже:

Может ли кто-нибудь проверить и сказать мне, правильный ли мой код или нет? Я и сам не уверен. Также, пожалуйста, скажите мне, есть ли лучший способ сделать это. Заранее спасибо, ребята.

ОТРЕДАКТИРОВАННЫЙ вопрос Я перечитал работу и ** должен применить функцию скользящего окна только к нескольким указанным пикселям. **

Пиксели, которые я должен применить, выясняются в этом методе. Обнаружен начальный контур изображения (показан на изображениях красным цветом). Затем вокруг контура на заданном расстоянии проводится полоса. Функция скользящего окна должна применяться только к тем пикселям, которые находятся в пределах узкой полосы для исходных изображений

Даю изображения, и начальные контуры, и изображения полос.

Пиксели, отмеченные белым, являются моими указанными пикселями, к которым должна быть применена функция скольжения. Можно ли применять nfilter по таким критериям? Пожалуйста помоги. Я уточню дополнительно, если мой вопрос неверен.