武汉理工大学数字图像处理英文课件 Chapter2（1）.pptVIP

Ellipses Polygon Function Polygon Demo * function img = ellipse(sz, cen, rad, theta, thres) x = 1:sz(2); y = 1:sz(1); x = ones(sz(1),1)*x; y = y*ones(1,sz(2)); theta = theta*pi/180; p = cos(theta)*(x-cen(1)) - sin(theta)*(y-cen(2)); q = sin(theta)*(x-cen(1)) + cos(theta)*(y-cen(2)); r = sqrt((p/rad(1)).^2+(q/rad(2)).^2); img = rthres; im = ellipse([200 200],[100 100],[80 10],20,1); imshow(im); imwrite(im,ellipse.tif); function z = polygon(x,y,p) if nargin3 theta = (2*pi/3)*(0:2) -pi/2; p = [cos(theta) sin(theta)]; end n = size(p,1); z = ones(size(x)); for i=1:n, u = p(mod(i,n)+1,:) - p(i,:); v = p(mod(i+1,n)+1,:) - p(i,:); if ((u(1)*v(2)-u(2)*v(1))0) u = -u; end s = u(1)*(y-p(i,2))-u(2)*(x-p(i,1)); z = z.*(s0); end pi Pi+1 Pi+2 (x, y) Test whether (x, y) is on the same side of the line defined by the points pi pi+1 as is the point pi+2 Repeat for each side of the polygon. Points that pass all tests are set to 1. theta = (pi/3)*(0:5); p = [cos(theta) sin(theta)]; z = polygon(x,y,p); imshow(z); imwrite(z0.5,hexagon.tif); [x y] = meshgrid(linspace(-1.5,1.5,251)); z = polygon(x,y); imshow(z); imwrite(z0.5,triangle.tif); *