Estoy tratando de entender este código:

d=edge(d,'canny',.6);
figure,
imshow(d,[])

ds = bwareaopen(d,40);
figure,
imshow(ds,[])

iout = d1;
BW=ds;

iout(:,:,1) = iout;
iout(:,:,2) = iout(:,:,1);
iout(:,:,3) = iout(:,:,1);
iout(:,:,2) = min(iout(:,:,2) + BW, 1.0);
iout(:,:,3) = min(iout(:,:,3) + BW, 1.0);

Entiendo qued es la imagen y se aplica el detector canny y se descuidan 40 píxeles. La imagen es escala de grises y el contorno se agrega a la imagen.

¿Puedes explicar las siguientes líneas? ¿Qué principio / algoritmo se usa aquí? Tengo problemas, especialmente con la parte de detección de contornos del código.