% dessin de la vigne de dessus

h=plot(1,1);
hold on
delete(h)
nbli=10;
nbco=20;
h=3;    % hauteur du piquet
oh=-1;   % hauteur de l'oeil
mur=5;
s=2;    % zoom

[x,y]=meshgrid(1:nbco,1:nbli);
y=y+mur+1;
plot(x,y,'ob')

% place de l'observateur
xo=4;
yo=3;
obs=plot(xo,yo,'or');

% projection sur le mur
for li=1:nbli
    for co=1:nbco
        xi=x(li,co);
        yi=y(li,co);
        dx=xo+(xi-xo).*(mur-yo)./(yi-yo);
        dx=s.*(dx-xo)+xo;
        dz1=mur+(oh.*(mur-yo))./(yi-yo);
        dz2=((h-oh).*(mur-yo))./(yi-yo);
        piq(li,co)=plot([dx dx],[dz1 dz1+dz2],'k');
    end
end

axis([0 nbco 0 nbli+mur+1])
axis equal
axis off
hold off

set(gcf,'color','w')
figure(gcf)
set(gcf,'doublebuffer','on')
set(gcf,'resize','off')
set(gcf,'windowButtonMotionFcn','clicvigne')



