MATLAB SPMD和nlfilter - 更改窗口大小

Question

我使用John Burkardt的CONTRAST2_SPMD的修改版本在大型灰度图像上运行3x3方差计算。这样做的好处是我可以使用PCT在本地机器上使用8个内核，但是我想改变窗口大小（目前为3x3），但是我想尝试几种窗口大小。我如何修改下面的代码以允许可变窗口大小？ （我将使用奇数，方形窗口，例如3x3,5x5,7x7,9x9等）

function y = parwinvar (x) 

%************************************************************************** 
% 
% 
%% PARWINVAR uses MATLAB's SPMD command for parallel windowed variance. 
% 
% Discussion: 
% Calculates windowed standard deviation (squared to get variance). 
% Based on CONTRAST2_SPMD by John Burkardt. 
% 
% 
% Parameters: 
% 
% Input, image X, the initial black and white image. 
% Output, image Y, the contrast-enhanced black and white image. 
% 
% 
% 
% Open the matlabpool 
% 
    matlabpool open local 8 
% 
% Check image is grayscale, if not covert it 
% 
    if ndims(x)>2 
     x=rgb2gray(x); 
    end 
% 
% Since the image is black and white, it is a 2D array. 
% Hence, it will be distributed by columns. 
% 
    xd = distributed (x); 
% 
% Have each worker enhance the contrast in its portion of the picture. 
% 
% You will see lots of meaningless error messages, because NLFILTER 
% wants to put out a "waitbar" telling you to wait. But the workers 
% don't have an associated display. 
% 
    spmd 

    xl = getLocalPart (xd); 
% 
% In order to use LABSENDRECEIVE, we need to reference the previous 
% and next labindex. 
% 
    if (labindex ~= 1) 
     previous = labindex - 1; 
    else 
     previous = numlabs; 
    end 

    if (labindex ~= numlabs) 
     next = labindex + 1; 
    else 
    next = 1; 
end 
% 
% Each worker sends its first column to the previous worker, 
% and receives corresponding data from the next worker. 
% 
column = labSendReceive (previous, next, xl(:,1)); 

if (labindex < numlabs) 
    xl = [ xl, column ]; 
end 
% 
% Each worker sends its last column to the next worker, 
% and receives corresponding data from the previous worker. 
% 
column = labSendReceive (next, previous, xl(:,end)); 

if (1 < labindex) 
    xl = [ column, xl ]; 
end 
% 
% Now do the enhancement. 
% We can only do a 3x3 neighborhood, because we're only sending 
% one column to the left and right. 
% 
xl = nlfilter (xl, [3,3], @std2); 
% 
% Now strip off the extra columns. 
% 
    if (1 < labindex) 
     xl = xl(:,2:end); 
    end 

    if (labindex < numlabs) 
     xl = xl(:,1:end-1); 
    end 

    xl = uint8 (xl); 

    end 
% 
% We are working with a black and white image, so we can simply 
% concatenate the submatrices to get the whole object. 
% 
    y = [ xl{:} ]; 
    y = (y .* y); % square to get variance 
% Close matlabpool 
matlabpool close 
    return 
end 

在此先感谢您的帮助。

Answer 1

您需要修改labSendReceive行以发送更多数据。尤其是，这样的事情：

% 
% Each worker sends its first column to the previous worker, 
% and receives corresponding data from the next worker. 
% 
column = labSendReceive (previous, next, xl(:,1)); 

需求，成为：

% 
% Each worker sends its first N columns to the previous worker, 
% and receives corresponding data from the next worker. 
% 
columns = labSendReceive (previous, next, xl(:,1:N)); 

其中N是(blocksize-1)/2。在发送右侧“幻象单元格”时以及删除“额外列”时，会做出相应的更改。