raster.m

超声图像实时消噪实验程序

function r = raster( si, sj, di, dj, wi, wj )
% make a raster image, black lines on white background
r = ones( si, sj );
[x,y] = ndgrid( 1:si, 1:sj );
r( ismember( mod(x,di), 0:(wi-1) ) ) = 0;
r( ismember( mod(y,dj), 0:(wj-1) ) ) = 0;