1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
| G = zeros(nz,nx,ny);
T1 = ones(WL,WX,WY);
T2 = ones(WL,WX,WY);
T3 = ones(WL,WX,WY);
wl = floor(WL/2);
wx = floor(WX/2);
wy = floor(WY/2);
for i = 1:nz-WL+1
for j = 1:nx-WX+1
for k = 1:ny-WY+1
num = sum(sum(sum(((X(i:i+WL-1,j:j+WX-1,k:k+WY-1)-(T1*mean(mean(mean(X(i:i+WL-1,j:j+WX-1,k:k+WY-1)))))) .* ...
(Y(i:i+WL-1,j:j+WX-1,k:k+WY-1)-(T2*mean(mean(mean(Y(i:i+WL-1,j:j+WX-1,k:k+WY-1))))))))));
den = sum(sum(sum((X(i:i+WL-1,j:j+WX-1,k:k+WY-1)-(T3*mean(mean(mean(X(i:i+WL-1,j:j+WX-1,k:k+WY-1)))))))).^2);
G(i+wl,j+wx,k+wy)=num./den;
end
end
end |
Partager