Окрестность в пределах евклидова расстояния R1
может быть получена следующим образом:
[x,y] = meshgrid(-R1:R1,-R1:R1);
r = x.^2 + y.^2;
k1 = r < R1.^2;
Свертка с этим ядром приводит для каждого выходного пикселя к сумме пикселей в его окрестности:
sum1 = conv2(img, k1);
(где img
- ваше изображение - убедитесь, что оно имеет тип с плавающей запятой, если оно целочисленного типа, возможно, что результат суммы не может быть представлен правильно).
Вы можете повторить тот же код выше для R2
, ведущий к sum2
, а затем вычислить ваш коэффициент:
result = sum1 ./ sum2;