function [JJ1,JJ2,JJ3,j1median,j2median,j3median]=swpbsdm(img1,img2,block)
% load image 
X = double(imread(img1,'bmp'));
XC = double(imread(img2,'bmp'));

[jj1r,jj2r,jj3r,JJ1R,JJ2R,JJ3R]=bsdm(X(:,:,1),XC(:,:,1),block);
[jj1g,jj2g,jj3g,JJ1G,JJ2G,JJ3G]=bsdm(X(:,:,2),XC(:,:,2),block);
[jj1b,jj2b,jj3b,JJ1B,JJ2B,JJ3B]=bsdm(X(:,:,3),XC(:,:,3),block);

clear X;
clear XC;


j1median=(jj1r(2)+jj1g(2)+jj1b(2))./3;
j2median=(jj2r(2)+jj2g(2)+jj2b(2))./3;
j3median=(jj3r(2)+jj3g(2)+jj3b(2))./3;
JJ1=(JJ1R+JJ1G+JJ1B)./3;
JJ2=(JJ2R+JJ2G+JJ2B)./3;
JJ3=(JJ3R+JJ3G+JJ3B)./3;