makehalfdomeimage.m

来自「beamlet变化的工具箱」· M 代码 · 共 49 行

function image = MakeHalfDomeImage(n,r, theta,dec)
% MakeHalfDome -- Make synthetic image of half dome
%	Usage
%		image = MakeHalfDomeImage(n,r, theta,dec)
%	Inputs
%		n		extent of delivered image, e.g. 256
%		r		parameter of halfspace mutilation
%		theta	parameter of halfspace mutilation
%		dec		parameter of anti-jaggie decimation,
%				optional, default 2		
%	Outputs
%		image	Anti-aliased version of mutilated gaussian
%
%	Description
%		The Gaussian exp(-r^2) is mutilated by the
%	halfspace indicator with parameters (r,theta)
%	

	if nargin < 4,
		dec = 2;
	end
	fac = 2^(dec-1);
	x   = (-fac*n):(fac*n-1);
	y   = x;
	[xx,yy] = meshgrid(x,y);
	BigImage = exp( - (xx./(fac*n/2)).^2 - (yy./(fac*n/2)).^2) ...
		.* ( (xx .* cos(theta) + yy .* sin(theta)) > (fac .*r) );
	image = HCoarsen2d(BigImage,dec);

% Copyright (c) 1999 David L. Donoho

%	Revision History:
%		1.0	09/05/99	 initial version DLD	
%		1.1	09/13/99	 changed line 27 to "> (fac .* r));"	

%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%