U_o = imread('wook_L30_G30.bmp');
U_o = double(U_o(:,:,1))/255;  % U_o (ox, oy) ~ not image size
index_z2=1;

for z2 = z2_under_image + z2_image : z2_series_step : z2_upper_image + z2_image
    
    
    final_image(:,:,index_z2)=conv2(abs(final_h(:,:,index_z2)),U_o,'same');
    
    
    index_z2
    index_z2=index_z2+1;
    
    
end

%final image should be normalized too

for as = 1:601
    final_image(:,:,as)=final_image(:,:,as)/sum(sum(final_image(:,:,as)));
end


slice=zeros(grid_size,601);
slice(:)=final_image(grid_size/2,:,:);
imagesc(abs(slice));