un programme pour appliquer la fonction de transformation suivante à une image en niveaux de gris

Question

Je veux appliquer à la suite fonction de transformation à une image en niveaux de gris, je sais comment l'appliquer à la fonction suivante,

ma question est comment puis-je appliquer un programme à la fonction de transformation suivante,

Code à ce jour,

clear; 
pollen = imread('Fig3.10(b).jpg'); 
u = double(pollen); 
[nx ny] = size(u) 
nshades = 256; 
r1 = 80; s1 = 10; % Transformation by piecewise linear function. 
r2 = 140; s2 = 245; 
for i = 1:nx 
for j = 1:ny 
if (u(i,j)< r1) 
uspread(i,j) = ((s1-0)/(r1-0))*u(i,j) 
end 
if ((u(i,j)>=r1) & (u(i,j)<= r2)) 
uspread(i,j) = ((s2 - s1)/(r2 - r1))*(u(i,j) - r1)+ s1; 
end 
if (u(i,j)>r2) 
uspread(i,j) = ((255 - s2)/(255 - r2))*(u(i,j) - r2) + s2; 
end 
end 
end 
hist= zeros(nshades,1); 
for i=1:nx 
for j=1:ny 
for k=0:nshades-1 
if uspread(i,j)==k 
hist(k+1)=hist(k+1)+1; 
end 
end 
end 
end 
plot(hist); 
pollenspreadmat = uint8(uspread); 
imwrite(pollenspreadmat, 'pollenspread.jpg'); 

Merci à l'avance

Answer 1

La figure dit que pour toutes les intensités qui sont entre A et B, ils doivent être mis à C. Tout ce que vous avez à faire est de modifier vos deux boucles for afin que, pour toute valeur comprise entre A et B, définissez l'emplacement de sortie sur C. Je suppose également que la gamme est inclusive. Vous pouvez simplement supprimer les première et dernière if conditions et utilisez celle du milieu:

for i = 1:nx 
    for j = 1:ny 
     if ((u(i,j)>=r1) && (u(i,j)<= r2)) 
      uspread(i,j) = C; 
     end 
    end 
end 

C est une constante que vous vous fixer. Habituellement, pour la segmentation, ce résultat est très élevé pour distinguer le premier plan de l'arrière-plan. Vous avez une image uint8 ici, donc C = 255; fonctionnerait.

Cependant, je vous recommande de réaliser une solution plus vectorisée. Évitez for boucles et utiliser logical indexation à la place:

uspread = u; 
uspread(u >= r1 & u <= r2) = C;