???
[pst-anamorphosis.git] / doc / pst-anamorphosis-doc-part1.tex
1 \listfiles
2 \documentclass[a4paper,fleqn]{article}
3 \usepackage[utf8]{inputenc}
4 \usepackage[T1]{fontenc}
5 \usepackage{amsmath,amssymb}
6 \usepackage[frenchb]{babel}
7 %\usepackage{mathptmx}
8 \usepackage[charter]{mathdesign}
9 \usepackage[margin=2.5cm]{geometry}
10 \usepackage{multicol}
11 \usepackage{listings}
12 \usepackage[svgnames]{xcolor}
13 \usepackage{showexpl}
14 \usepackage[nomessages]{fp}
15 \usepackage{xspace}
16 \usepackage{pst-plot,pst-solides3d,pst-anamorphosis-add,pst-3d}
17 \usepackage{pst-grad}
18 \usepackage[absolute,notitlepage]{pst-abspos}
19 \usepackage{url}
20
21 \psset{path=C:/Dokumente und Einstellungen/Besitzer/Desktop/bergen/bergen/}
22 %\def\epsRoot{C:/Dokumente und Einstellungen/Besitzer/Desktop/bergen/bergen/}
23
24 \renewcommand{\ttdefault}{lmtt}
25
26 \definecolor{syracuseGRIS}{HTML}{C1C1C1}
27 \definecolor{syracuseVERT}{HTML}{029235}
28
29 \definecolor{sepia}{rgb}{1,0.8,0.8}
30 \definecolor{grisclair}{rgb}{0.8,0.8,0.8}
31 \definecolor{BleuCiel}{cmyk}{0.2,0,0,0}
32 \definecolor{OrangePale}{cmyk}{0,0.2,0.4,0}
33
34 \lstset{%
35 language=[LaTeX]TeX,%
36 float=hbp,%
37 basicstyle=\ttfamily\small, %
38 texcsstyle=*\color{blue},%
39 identifierstyle=\color{black}, %
40 keywordstyle=\color{syracuseVERT}, %
41 otherkeywords={$, \{, \}, \[, \]},
42 stringstyle=\color{syracuseVERT}, %
43 commentstyle=\color{syracuseVERT}, %
44 backgroundcolor=\color{syracuseGRIS!30},%
45 columns=flexible, %
46 tabsize=4, %
47 frame=single, %
48 %frame=shadowbox, %
49 %rulesepcolor=\color{syracuseGRIS!30},%
50 extendedchars=true, %
51 showspaces=false, %
52 showstringspaces=false, %
53 numbers=left,
54 numbersep=0.8em,
55 numberstyle=\tiny, %
56 breaklines=true, %
57 breakautoindent=true,%
58 captionpos=b,%
59 xleftmargin=1em,%
60 sensitive=true,%
61 morekeywords=[6]{pspicture,center},%
62 keywordstyle=[6]\color{FireBrick},%
63 %morekeywords=[7]{(,)},%
64 %keywordstyle=[7]\color{syracuseVERT}
65 }
66
67 \renewcommand{\lstlistingname}{Source}
68 \renewcommand{\labelitemi}{$\bullet$}
69
70
71 \newcommand\cs[1]{\texttt{\char`\\#1}}
72 \newcommand\file[1]{\texttt{#1}}
73
74 \parindent0pt
75 \parskip5pt
76
77 \def\syracuseTitle{Les anamorphoses : pr\'{e}sentation th\'{e}orique}
78 %\def\syracuseGraphic{eiffel2}
79
80
81 \begin{document}
82
83 %% === BEGIN == Page de garde =================================================
84
85 \thispagestyle{empty}
86
87 \pstPutAbs(0,-29.7){%
88 \begin{pspicture}(0,0)(21,29.7)
89 \pspolygon[fillstyle=solid,linecolor=syracuseVERT,fillcolor=syracuseVERT](0,0)(10.5,14.85)(21,0)
90 \pspolygon[fillstyle=solid,linecolor=syracuseGRIS,fillcolor=syracuseGRIS](0,0)(21,29.7)(0,29.7)
91 \pspolygon[fillstyle=solid,linecolor=syracuseGRIS!50,fillcolor=syracuseGRIS!50](21,0)(10.5,14.85)(21,29.7)
92 \end{pspicture}
93 }
94 \pstPutAbs(2.5,-3.75){%
95 \includegraphics[scale=1]{pst-anamorphosis}
96 }
97 \pstPutAbs(2.5,-5.25){%
98 \LARGE \textbf{\syracuseTitle}
99 }
100 \pstPutAbs(2.5,-13.5){%
101 \begin{pspicture}(0,0)(12,12)
102 \rput(4,4){\includegraphics[height=8cm]{Eiffel}}
103 %\rput(4,4){\includegraphics[height=8cm]{\syracuseGraphic}}
104 %\psframe(0,0)(8,8)
105 \end{pspicture}
106 }
107 \pstPutAbs(12.5,-15){%
108 \parbox{0.4\textwidth}{\Large\raggedleft
109 {\LARGE\textbf{Contributeurs}}\\[0.2cm]
110 J\"{u}rgen \textsc{Gilg}\\
111 Manuel \textsc{Luque}\\
112 Jean-Michel \textsc{Sarlat}
113 }}
114 \vfill
115 \begin{center}
116 \textcolor{white}{\textbf{\today}}\\[0.3cm]
117 \textcolor{white}{\url{http://melusine.eu.org/syracuse/G/pst-anamorphosis/}}\\
118 \includegraphics[scale=0.4]{logo_syracuse}
119 \end{center}
120
121 %% == END == Page de garde ====================================================
122
123 \newpage
124
125 \section{L'anamorphose cylindrique}
126
127 On place \`{a} l'int\'{e}rieur du cylindre l'image telle qu'elle doit \^{e}tre vue par un observateur regardant dans le miroir cylindrique (on peut la placer \`{a} l'ext\'{e}rieur, mais il faut que les rayons lumineux rencontrent toujours le cylindre, il faut donc veiller aux dimensions). L'objet anamorphique est <<~l'objet d\'{e}form\'{e}~>> dont le miroir reconstituera les proportions r\'{e}elles.
128
129 Objet et image ob\'{e}issent aux lois de la r\'{e}flexion de l'optique g\'{e}om\'{e}trique :
130 \begin{itemize}
131 \item rayon incident et rayon r\'{e}fl\'{e}chi appartiennent \`{a} un m\^{e}me plan ;
132 \item rayon incident et rayon r\'{e}fl\'{e}chi sont sym\'{e}triques par rapport \`{a} la normale au miroir au point d'incidence.
133 \end{itemize}
134 \input{fig3d-anacyl.tex}
135 L'image non d\'{e}form\'{e}e (celle qui est vue dans le miroir) est plac\'{e}e, dans cet exemple, au centre du miroir. Un rayon incident partant de l'objet anamorphique se r\'{e}fl\'{e}chit sur le miroir et apr\`{e}s r\'{e}flexion parvient \`{a} l'{\oe}il de notre observateur. L'observateur a l'illusion que le rayon provient du point image. Il faut donc reconstruire math\'{e}matiquement la marche d'un tel rayon lumineux en partant de l'image dans le miroir.
136
137 L'observateur est suffisamment \'{e}loign\'{e} du miroir pour pouvoir \^{e}tre consid\'{e}r\'{e} comme ponctuel.
138
139 Soit $P$ un point de l'image(not\'{e} $A'$ dans le sch\'{e}ma ci-apr\`{e}s), $V$ l'{\oe}il de l'observateur. Tra\c{c}ons un droite $(PV)$ et d\'{e}terminons le point d'intersection $I$ avec le cylindre : c'est le point d'incidence.
140 \[
141 V(x_V,y_V,z_V)\quad\text{ et }\quad P(x_P,y_P,0)
142 \]
143 L'\'{e}quation param\'{e}trique de la droite $(PV)$ s'\'{e}crit $\overrightarrow{IV}=\rho\overrightarrow{PV}$:
144 \begin{equation}\label{eq:paracyl}
145 \left\lbrace
146 \begin{array}{lcl}
147 x_V-x_I&=&\rho(x_V-x_P)\\
148 y_V-y_I&=&\rho(y_V-y_P)\\
149 z_V-z_I&=&\rho(z_V-0)
150 \end{array}
151 \right.
152 \Longrightarrow
153 \left\lbrace
154 \begin{array}{lcl}
155 x_I&=&x_V(1-\rho)+\rho x_P\\
156 y_I&=&y_V(1-\rho)+\rho y_P\\
157 z_I&=&z_V(1-\rho)
158 \end{array}
159 \right.
160 \end{equation}
161 \begin{center}
162 \begin{pspicture}(-4,-0.5)(6.5,8)
163 \pnode(6,7){V}
164 \uput[0](V){$V$}
165 \pnode(3,6.5){S}
166 \pnode(-3,6.5){S'}
167 \pnode(3,2.8){I}
168 \pnode(1,0){P}
169 \pnode(5,0){P'}
170 \uput[45](P'){$P'$}
171 \pnode(3,0){G}
172 \pnode(-3,0){G'}
173 \pspolygon[fillstyle=solid,fillcolor=lightgray!30](G')(G)(S)(S')
174 \psaxes[ticks=none,labels=none]{->}(0,0)(-4,0)(6.5,7.5)
175 \uput[0](6.5,0){$x$}
176 \uput[90](0,7.5){$z$}
177 \uput[135](P){$P$}
178 \uput[-45](G){$G$}
179 \uput[135](I){$I$}
180 \psline(V)(P)
181 \psline(S)(G)
182 \rput(I){%
183 \psline[linestyle=dashed](3;0)%
184 \uput[0](3;0){$N$}%
185 \psarc[doubleline=true](0,0){1}{-90}{-54.46}%
186 \psarc[doubleline=true](0,0){1}{54.46}{90}
187 \uput[72](1.2;72){$\varepsilon$}%
188 \uput[-72](1.2;-72){$\varepsilon$}%
189 \rput{-90}(0,0){\psframe(0.3,0.3)}
190 }
191 \psline[linecolor=red](V)(I)(P')
192 \pcline[nodesepB=1.5,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(P')(I)
193 \pcline[nodesepB=2.5,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(I)(V)
194 \qdisk(I){2pt}
195 \qdisk(P){2pt}
196 \qdisk(V){2pt}
197 \qdisk(P'){2pt}
198 \qdisk(G){2pt}
199 \psline[linewidth=0.5\pslinewidth]{|<->|}(0,-0.2)(1,-0.2)
200 \uput[-90](0.5,-0.2){$r_P$}
201 \psline[linewidth=0.5\pslinewidth]{<->}(0,2.8)(3,2.8)
202 \uput[-90](1.5,2.8){$r_I=R$}
203 \end{pspicture}
204 \end{center}
205 Le point $I$ appartenant au cylindre, ses coordonn\'{e}es v\'{e}rifient la relation :
206 \begin{equation}\label{eq:cylindre}
207 x_I^2+y_I^2=R^2
208 \end{equation}
209 (\ref{eq:paracyl}) en (\ref{eq:cylindre}) et apr\`{e}s d\'{e}veloppement, on obtient l'\'{e}quation du second degr\'{e} en $\rho$:
210 \begin{gather*}
211 \left(x_V(1-\rho)+\rho x_P\right)^2+\left(y_V(1-\rho)+\rho y_P\right)^2=R^2\\
212 x_V^2(1-2\rho+\rho^2)+2(1-\rho)\rho x_Vx_P+\rho^2 x_P^2+y_V^2(1-2\rho+\rho^2)+2(1-\rho)\rho y_Vy_P+\rho^2 y_P^2=R^2\\
213 (x_V^2+y_V^2-2x_Vx_P-2y_Vy_P)\rho^2+2(-x_V^2+x_Vx_P-y_V^2+y_Vy_P)\rho+x_V^2+y_V^2=R^2
214 \end{gather*}
215 Comparaison avec
216 \[
217 a\rho^2+2b'\rho+c=0
218 \]
219 donne :
220 \[
221 \left\lbrace
222 \begin{array}{lcl}
223 a&=&(x_V-x_P)^2+(y_V-y_P)^2\\
224 2b'&=&x_Vx_P+y_Vy_P-x_V^2-y_V^2\\
225 c&=&x_V^2+y_V^2-R^2
226 \end{array}
227 \right.
228 \]
229 La r\'{e}solution de cette \'{e}quation nous donne les solutions classiques:
230 \[
231 \left\lbrace
232 \begin{array}{lcl}
233 \rho'&=&\dfrac{-b'+\sqrt{\Delta'}}{a}\\[0.5cm]
234 \rho''&=&\dfrac{-b'-\sqrt{\Delta'}}{a}
235 \end{array}
236 \right.
237 \qquad \Delta'=b'^2-ac
238 \]
239 On retiendra la plus petite valeur positive des deux, que par la suite j'appelle $\rho$.
240
241 $(IV)$ repr\'{e}sente le rayon r\'{e}fl\'{e}chi par le miroir. Le rayon incident est d\'{e}fini par la droite sym\'{e}trique de $(IV)$ par rapport \`{a} la normale au miroir en $I$. Je cherche le sym\'{e}trique de $V$, nomm\'{e} $V'$ par rapport \`{a} cette normale $(IN)$. Ce point $V'$ remplit deux conditions :
242 \begin{enumerate}
243 \item $\overrightarrow{IV}+\overrightarrow{IV'}=k\overrightarrow{IN}$
244 \item $\overrightarrow{VV'}.\overrightarrow{IN}=0$
245 \end{enumerate}
246 La normale $(IN)$ a pour vecteur directeur $\overrightarrow{IN}(x_I,y_I,0)$
247
248 La premi\`{e}re condition se traduit par :
249 \[
250 \left\lbrace
251 \begin{array}{lcl}
252 x_V-x_I+x_{V'}-x_I&=&kx_I\\
253 y_V-x_I+y_{V'}-y_I&=&ky_I\\
254 z_V-z_I+z_{V'}-z_I&=&0
255 \end{array}
256 \right.
257 \Longrightarrow
258 \left\lbrace
259 \begin{array}{lcl}
260 x_{V'}&=&kx_I+2x_I-x_V\\
261 y_{V'}&=&ky_I+2y_I-y_V\\
262 z_{V'}&=&2z_I-z_V
263 \end{array}
264 \right.
265 \]
266 La deuxi\`{e}me par :
267 \[
268 (x_{V'}-x_V)x_I+(y_{V'}-y_V)y_I=0
269 \]
270 En rempla\c{c}ant $x_{V'}$ et $y_{V'}$ tir\'{e}s de la premi\`{e}re condition dans la deuxi\`{e}me :
271 \begin{gather*}
272 k(x_I^2+y_I^2)+2x_I^2-2x_Vx_I+2y_I^2-2y_Vy_I=0\\
273 kR^2+2R^2=2(x_Vx_I+y_Vy_I)\\
274 k+2=\dfrac{2}{R^2}(x_Vx_I+y_Vy_I)
275 \end{gather*}
276 Les coordonn\'{e}es de $V'$ s'en d\'{e}duisent :
277 \[
278 \left\lbrace
279 \begin{array}{lcl}
280 x_{V'}&=&(k+2)x_I-x_V\\
281 y_{V'}&=&(k+2)y_I-y_V\\
282 z_{V'}&=&z_V(1-2\rho)
283 \end{array}
284 \right.
285 \]
286 Il reste \`{a} trouver l'intersection de $(IV')$ avec le plan horizontal $z=0$.
287
288 \'Equation param\'{e}trique de $(IV')$, $M$ \'{e}tant un point courant : $\overrightarrow{MV'}=\alpha\overrightarrow{IV'}$
289 \[
290 \left\lbrace
291 \begin{array}{lcl}
292 x_{V'}-x&=&\alpha(x_{V'}-x_I)\\
293 y_{V'}-y&=&\alpha(y_{V'}-y_I)\\
294 z_{V'}-z&=&\alpha(z_{V'}-z_I)
295 \end{array}
296 \right.
297 \]
298 $z=0\Longrightarrow \alpha=\dfrac{z_{V'}}{z_{V'}-z_I}$ soit
299 \[
300 \alpha=\dfrac{1-2\rho}{-\rho}
301 \]
302 En rempla\c{c}ant $\alpha$ par son expression, nous obtenons les coordonn\'{e}es du point $P'$ de l'objet anamorphique.
303 \[
304 \left\lbrace
305 \begin{array}{lcl}
306 x_{P'}&=&x_{V'}-\alpha(x_{V'}-x_I)\\
307 y_{P'}&=&y_{V'}-\alpha(y_{V'}-y_I)
308 \end{array}
309 \right.
310 \]
311 Cette s\'{e}rie de calculs doit \^{e}tre appliqu\'{e}e \`{a} tous les points de l'image <<~normale~>> afin d'obtenir l'objet anamorphique (d\'{e}form\'{e}) dont le miroir <<~redressera~>> la forme.
312
313 On notera que \textit{la cote de l'observateur $z_V$ n'intervient pas}. On le comprend ais\'{e}ment en faisant un dessin du plan vertical passant par l'\oe{}il et l'axe du cylindre miroir. La position de la projection horizontale \'{e}tant fix\'{e}e $(x_V,y_V)$, quelle que soit la valeur de $z_V$, $A$ et $A'$ \'{e}tant sym\'{e}triques par rapport \`{a} la g\'{e}n\'{e}ratrice du miroir appartenant au plan vertical choisi, si $A'$ est donn\'{e} alors $A$ est fix\'{e} quelque soit~$z_V$.
314
315 \newpage
316
317 \section{L'anamorphose conique}
318
319 Le principe est identique \`{a} celui de l'anamorphose cylindrique : imaginons un rayon lumineux provenant de l'objet <<~anamorphique~>>, se r\'{e}fl\'{e}chissant sur le miroir conique et parvenant \`{a} l'{\oe}il de l'observateur plac\'{e} au-dessus et dans l'axe du c\^{o}ne \`{a} une position suffisamment haute pour que l'observateur puisse \^{e}tre consid\'{e}r\'{e} comme ponctuel. Ainsi l'observateur aura l'illusion d'observer l'image reconstitu\'{e}e par le miroir conique. Image et objet sont dans le plan horizontal.
320 \input{fig3d-anacon.tex}
321 Il s'agit de d\'{e}terminer en premier, l'intersection de $(VP)$ avec le c\^{o}ne, qu'on appelle $I$ (point d'incidence). Les coordonn\'{e}es de $V$, $S$ et $P$ sont not\'{e}es : $V(0,0,z_V)$, $S(0,0,z_S)$ et $P(x_P,y_P,0)$.
322
323 L'\'{e}quation param\'{e}trique de la droite $(PV)$ s'\'{e}crit $\overrightarrow{IV}=\lambda\overrightarrow{PV}$:
324 \[
325 \left\lbrace
326 \begin{array}{lcl}
327 0-x_I&=&\lambda(0-x_P)\\
328 0-y_I&=&\lambda(0-y_P)\\
329 z_V-z_I&=&\lambda(z_V-0)
330 \end{array}
331 \right.
332 \Longrightarrow
333 \left\lbrace
334 \begin{array}{lcl}
335 x_I&=&\lambda x_P\\
336 y_I&=&\lambda y_P\\
337 z_I&=&(1-\lambda)z_V
338 \end{array}
339 \right.
340 \]
341 On pose :
342 \[
343 r_I^2=x_I^2+y_I^2\quad\textrm{et}\quad r_P^2=x_P^2+y_P^2\quad\textrm{et}\quad |\overrightarrow{OG}|=R
344 \]
345 Le point $I$ appartenant au c\^{o}ne, ses coordonn\'{e}es v\'{e}rifient la relation (th\'{e}or\`{e}me de Thal\`{e}s):
346 \begin{align*}
347 \frac{R}{z_S}&=\frac{r_I}{z_S-z_I}\\
348 \frac{R}{z_S}&=\frac{\lambda r_P}{z_S-(1-\lambda)z_P}\\
349 \lambda&=\frac{R(z_S-z_V)}{r_Pz_S-R z_V}
350 \end{align*}
351 \begin{center}
352 \begin{pspicture}(-4,-1)(6.5,7.75)
353 \pnode(0,10){V}
354 \uput[0](V){$V$}
355 \pnode(0,5){S}
356 \uput[0](S){$S$}
357 \pnode(1.5,2.5){I}
358 \pnode(2,0){P}
359 \pnode(4.545,0){P'}
360 \uput[45](P'){$P'$}
361 \psdots[dotstyle=|](P')
362 \pnode(3,0){G}
363 \pnode(-3,0){G'}
364 \pspolygon[fillstyle=solid,fillcolor=lightgray!30](G')(G)(S)
365 \psaxes[ticks=none,labels=none]{->}(0,0)(-4,0)(5.5,10.25)
366 \uput[0](5.5,0){$x$}
367 \uput[90](0,10.25){$z$}
368 \uput[45](P){$P$}
369 \uput[-45](G){$G$}
370 \uput[70](I){$I$}
371 \psline(V)(P)
372 \psline(S)(G)
373 \rput(I){%
374 \psline[linestyle=dashed](4;30.96)
375 \uput[15.5](4;30.1){$N$}
376 \rput{-59.036}(0,0){\psframe(0.5,0.5)}
377 \psarc[doubleline=true](0,0){1}{101.31}{120.964}
378 %\psarc[doubleline=true](0,0){1}{-79.09}{-59.036}
379 \psarc[doubleline=true](0,0){1.2}{-59.036}{-39.38}
380 \pnode(1.2;-50){I1}
381 \pnode(1.2;112){I2}
382 \uput[-50](I1){$\varepsilon$}
383 \uput[112](I2){$\varepsilon$}
384 }
385 \rput(P'){\pswedge[fillstyle=solid,fillcolor=lightgray](0,0){1}{140.175}{180}
386 \uput[160](1;160){$\varepsilon'$}}
387 \rput(S){\psarc[linecolor=gray,linewidth=2\pslinewidth](0,0){1}{-90}{-59.036}
388 \uput[-75](1;-75){$\theta$}}
389 \rput(V){\psarc(0,0){2}{-90}{-78.69}
390 \uput[-80](2;-85){$\beta$}}
391 \psline[linecolor=red](V)(I)(P')
392 \pcline[nodesepB=1.5,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(P')(I)
393 \pcline[nodesepB=4,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(I)(V)
394 \qdisk(I){2pt}
395 \qdisk(P){2pt}
396 \qdisk(S){2pt}
397 \qdisk(V){2pt}
398 \qdisk(P'){2pt}
399 \qdisk(G){2pt}
400 \psline[linewidth=0.5\pslinewidth]{|<->|}(0,-0.2)(2,-0.2)
401 \uput[-90](1,-0.2){$r_P$}
402 \psline[linewidth=0.5\pslinewidth]{|<->|}(0,-0.8)(3,-0.8)
403 \uput[-90](1.5,-0.8){$R$}
404 \psline[linewidth=0.5\pslinewidth]{<->}(1.5,2.5)(0,2.5)
405 \uput[-90](0.75,2.5){$r_I$}
406 \psline[linewidth=0.5\pslinewidth]{<->|}(-0.2,0)(-0.2,5)
407 \uput[180](-0.2,2.5){$z_S$}
408 \psline[linewidth=0.5\pslinewidth]{<->|}(-1,0)(-1,10)
409 \uput[180](-1,5){$z_V$}
410 \psline[linewidth=0.5\pslinewidth,linestyle=dashed](1.5,2.5)(1.5,0)
411 \end{pspicture}
412 \end{center}
413 Pour construire le rayon incident, $(P'I)$, $(IV)$ est le rayon r\'{e}fl\'{e}chi, d\'{e}terminons~$\varepsilon'$.
414
415 Un raisonnement g\'{e}om\'{e}trique \'{e}l\'{e}mentaire nous montre que :
416 \[
417 \varepsilon'=90^\circ-2\theta+\beta
418 \]
419 avec
420 \[
421 \beta=\arctan\frac{r_P}{z_V}
422 \]
423 et
424 \[
425 \theta=\arctan\frac{R}{z_S}\quad\textrm{demi-angle au sommet du c\^{o}ne }
426 \]
427 Dans le plan horizontal, les coordonn\'{e}es de $P'$ sont :
428 \[
429 \left\lbrace
430 \begin{array}{lcl}
431 x_{P'}&=&x_I+\dfrac{z_I}{\tan \varepsilon'}\\[0.5cm]
432 y_{P'}&=&y_I+\dfrac{z_I}{\tan \varepsilon'}\\
433 \end{array}
434 \right.
435 \]
436
437 \newpage
438
439 \section{L'anamorphose sph\'{e}rique}
440
441 On place \`{a} l'int\'{e}rieur de la demi-sph\`{e}re l'image telle qu'elle doit \^{e}tre vue par un observateur regardant dans le miroir sph\'{e}rique (on peut la placer \`{a} l'ext\'{e}rieur, mais il faut que les rayons lumineux rencontrent toujours la sph\`{e}re, il faut donc veiller aux dimensions). L'objet anamorphique est <<~l'objet d\'{e}form\'{e}~>> dont le miroir reconstituera les proportions r\'{e}elles.\par Objet et image ob\'{e}issent aux lois de la r\'{e}flexion de l'optique g\'{e}om\'{e}trique :
442 \begin{itemize}
443 \item rayon incident et rayon r\'{e}fl\'{e}chi appartiennent \`{a} un m\^{e}me plan ;
444 \item rayon incident et rayon r\'{e}fl\'{e}chi sont sym\'{e}triques par rapport \`{a} la normale au miroir au point d'incidence.
445 \end{itemize}
446 \input{fig3d-anasphere.tex}
447 L'image non d\'{e}form\'{e}e (celle qui est vue dans le miroir) est plac\'{e}e, dans cet exemple, au centre du miroir. Un rayon incident partant de l'objet anamorphique se r\'{e}fl\'{e}chit sur le miroir et apr\`{e}s r\'{e}flexion parvient \`{a} l'{\oe}il de notre observateur. L'observateur a l'illusion que le rayon provient du point image. Il faut donc reconstruire math\'{e}matiquement la marche d'un tel rayon lumineux en partant de l'image dans le miroir.
448
449 L'observateur est suffisamment \'{e}loign\'{e} du miroir pour pouvoir \^{e}tre consid\'{e}r\'{e} comme ponctuel.
450
451 Soit $P$ un point de l'image, $V$ l'{\oe}il de l'observateur. Tra\c{c}ons un droite $PV$ et d\'{e}terminons le point d'intersection $I$ avec la sph\`{e}re : c'est le point d'incidence.
452 \[
453 V(x_V,y_V,z_V)\quad\text{ et }\quad P(x_P,y_P,0)
454 \]
455 \begin{center}
456 \shorthandoff{!}
457 \begin{pspicture}(-4,-0.5)(8,5.5)
458 \pswedge[fillstyle=solid,fillcolor=lightgray!30](0,0){2}{0}{180}
459 \psaxes[labels=none,ticks=none]{->}(0,0)(-3,0)(7.5,5)
460 \uput[0](7.5,0){$x$}
461 \uput[90](0,5){$z$}
462 \pnode(1.5,0){P}
463 \pnode(0,0){O}\uput[dl](O){$O$}
464 \pnode(!/alpha 15 def /xI 2 alpha cos mul def /yI 2 alpha sin mul def /beta yI xI 1.5 sub atan def xI yI){I}
465 \psline(P)(I)
466 \rput(I){\pnode(!beta cos 5 mul beta sin 5 mul){V}}
467 \psline[linecolor=red](I)(V)
468 \rput(I){\psarc[doubleline=true](0,0){1}{!alpha}{!beta}
469 \psarc[doubleline=true](0,0){1.2}{!2 alpha mul beta sub}{!alpha}}
470 \uput{1.2}[!alpha beta add 2 div](I){$\varepsilon$}
471 \uput{1.4}[!3 alpha mul beta sub 2 div](I){$\varepsilon$}
472 \pnode(!yI neg xI 2 alpha mul beta sub tan mul add 2 alpha mul beta sub tan div 0){M}
473 \psline[linecolor=red](M)(I)
474 \pcline[nodesepB=0.5,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(M)(I)
475 \pcline[nodesepB=2,linecolor=red,arrowsize=0.175,arrowinset=0.075]{->>}(I)(V)
476 \psline[linewidth=0.5\pslinewidth]{<->}(O)(I)
477 \uput[75](1;15){$R$}
478 \psline[linewidth=0.5\pslinewidth]{<->}(O)(V)
479 \uput[75](3,2.7){$r_V$}
480 \psline[linewidth=0.5\pslinewidth]{|<->|}(0,-0.3)(1.5,-0.3)
481 \uput[-90](0.75,-0.3){$r_P$}
482 \rput(I){%
483 \psline[linestyle=dashed](5;15)
484 \uput[15](5;15){$N$}
485 \psline[linewidth=0.5\pslinewidth,linecolor=gray](1.5;105)
486 \uput[105](1.5;105){$T$}
487 \psline[linewidth=0.5\pslinewidth,linecolor=gray](1.25;-75)
488 \rput{-75}(0,0){\psframe(0.3,0.3)}
489 }
490 \uput[-45](P){$P$}
491 \uput[-45](M){$P'$}
492 \uput[70](I){$I$}
493 \uput[u](V){$V$}
494 \psdot(P)
495 \psdot(I)
496 \psdot(V)
497 \psdot(M)
498 \end{pspicture}
499 \end{center}
500 L'\'{e}quation param\'{e}trique de la droite $(PV)$ s'\'{e}crit $\overrightarrow{IV}=\lambda\overrightarrow{PV}$:
501 \begin{equation}\label{eq:para}
502 \left\lbrace
503 \begin{array}{lcl}
504 x_V-x_I&=&\lambda(x_V-x_P)\\
505 y_V-y_I&=&\lambda(y_V-y_P)\\
506 z_V-z_I&=&\lambda(z_V-0)
507 \end{array}
508 \right.
509 \Longrightarrow
510 \left\lbrace
511 \begin{array}{lcl}
512 x_I&=&x_V(1-\lambda)+\lambda x_P\\
513 y_I&=&y_V(1-\lambda)+\lambda y_P\\
514 z_I&=&z_V(1-\lambda)
515 \end{array}
516 \right.
517 \end{equation}
518 On pose :
519 \[
520 r_V^2=x_V^2+y_V^2+z_V^2\quad\textrm{et}\quad r_P^2=x_P^2+y_P^2
521 \]
522 Le point $I$ appartenant \`{a} la sph\`{e}re, ses coordonn\'{e}es v\'{e}rifient la relation :
523 \begin{equation}\label{eq:sphere}
524 x_I^2+y_I^2+z_I^2=R^2
525 \end{equation}
526 (\ref{eq:para}) en (\ref{eq:sphere})
527 \begin{gather*}
528 \left(x_V+\lambda(x_P-x_V)\right)^2+\left(y_V+\lambda(y_P-y_V)\right)^2+\left((1-\lambda)z_V\right)^2=R^2\\
529 x_V^2+2\lambda x_V(x_P-x_V)+\lambda^2(x_P-x_V)^2+y_V^2+2\lambda y_V(y_P-y_V) +\lambda^2(y_P-y_V)^2+(1-2\lambda+\lambda^2)z_V^2=R^2
530 \end{gather*}
531 Apr\`{e}s d\'{e}veloppement, on obtient l'\'{e}quation du second degr\'{e} en
532 $\lambda$:
533 \begin{gather*}
534 \lambda^2\left((x_P-x_V)^2+(y_P-y_V)^2+z_V^2\right)+2\lambda\left(x_V(x_P-x_V)+y_V(y_P-y_V)-z_V^2\right)+x_V^2+y_V^2+z_V^2-
535 R^2=0\\
536 \lambda^2\left(x_V^2+y_V^2+z_V^2+x_P^2+y_P^2-2x_Px_V-2y_Py_V\right)+2\lambda\left(-x_V^2-y_V^2-z_V^2+x_Px_V+y_Py_V\right)+x_V^2+y_V^2+z_V^2-
537 R^2=0
538 \end{gather*}
539 Comparaison avec
540 \[
541 a\lambda^2+2b'\lambda+c=0
542 \]
543 donne pour le coefficient $a$ de $\lambda^2$ :
544 \[
545 a=x_V^2+y_V^2+z_V^2+x_P^2+y_P^2-2x_Px_V-2y_Py_V
546 \]
547 Pour le coefficient $2b'$ de $\lambda$ :
548 \[
549 2b'=-x_V^2-y_V^2-z_V^2+x_Px_V+y_Py_V
550 \]
551 Pour le coefficient $c$ :
552 \[
553 c=x_V^2+y_V^2+z_V^2-R^2
554 \]
555 Alors
556 \[
557 a\lambda^2+2b'\lambda+c=0
558 \]
559 avec :
560 \[
561 \left\lbrace
562 \begin{array}{lcl}
563 a&=&x_V^2+y_V^2+z_V^2+x_P^2+y_P^2-2x_Px_V-2y_Py_V\\
564 2b'&=&-x_V^2-y_V^2-z_V^2+x_Px_V+y_Py_V\\
565 c&=&x_V^2+y_V^2+z_V^2-R^2
566 \end{array}
567 \right.
568 \]
569 \[\left\lbrace
570 \begin{array}{lcl}
571 a&=&r_V^2+r_P^2-2x_Px_V-2y_Py_V\\
572 2b'&=&-r_V^2+x_Px_V+y_Py_V\\
573 c&=&r_V^2-R^2
574 \end{array}
575 \right.
576 \]
577 La r\'{e}solution de cette \'{e}quation nous donne les solutions classiques:
578 \[
579 \left\lbrace
580 \begin{array}{lcl}
581 \lambda'&=&\dfrac{-b'+\sqrt{\Delta'}}{a}\\[0.5cm]
582 \lambda''&=&\dfrac{-b'-\sqrt{\Delta'}}{a}
583 \end{array}
584 \right.
585 \qquad \Delta'=b'^2-ac
586 \]
587 On retiendra la valeur positive.% et \texttt{Coeff1} dans le programme.
588
589 $(IV)$ repr\'{e}sente le rayon r\'{e}fl\'{e}chi par le miroir. Le rayon incident est d\'{e}fini par la droite sym\'{e}trique de $(IV)$ par rapport \`{a} la normale au miroir en $I$. Je cherche le sym\'{e}trique de $V$, nomm\'{e} $V'$ par rapport \`{a} cette normale $(IN)$. Ce point $V'$ remplit deux conditions :
590 \begin{enumerate}
591 \item $\overrightarrow{IV}+\overrightarrow{IV'}=k\overrightarrow{IN}$
592 \item $\overrightarrow{VV'}.\overrightarrow{IN}=0$
593 \end{enumerate}
594 La normale $(IN)$ a pour vecteur directeur $\overrightarrow{IN}(x_I,y_I,z_I)$
595
596 La premi\`{e}re condition se traduit par :
597 \[
598 \left\lbrace
599 \begin{array}{lcl}
600 x_V-x_I+x_{V'}-x_I&=&kx_I\\
601 y_V-x_I+y_{V'}-y_I&=&ky_I\\
602 z_V-z_I+z_{V'}-z_I&=&kz_I
603 \end{array}
604 \right.
605 \Longrightarrow
606 \left\lbrace
607 \begin{array}{lcl}
608 x_{V'}&=&kx_I+2x_I-x_V\\
609 y_{V'}&=&ky_I+2y_I-y_V\\
610 z_{V'}&=&kz_I+2z_I-z_V
611 \end{array}
612 \right.
613 \]
614 La deuxi\`{e}me par :
615 \[
616 (x_{V'}-x_V)x_I+(y_{V'}-y_V)y_I+(z_{V'}-z_V)z_I=0
617 \]
618 En rempla\c{c}ant $x_{V'}$, $y_{V'}$ et $z_{V'}$ tir\'{e}s de la premi\`{e}re condition dans la deuxi\`{e}me :
619 \begin{gather*}
620 k(x_I^2+y_I^2+y_I^2)+2x_I^2-2x_Vx_I+2y_I^2-2y_Vy_I+2z_I^2-2z_Vz_I=0\\
621 kR^2+2R^2=2(x_Vx_I+y_Vy_I+z_Vz_I)\\
622 k+2=\dfrac{2}{R^2}(x_Vx_I+y_Vy_I+z_Vz_I)
623 \end{gather*}
624 Les coordonn\'{e}es de $V'$ s'en d\'{e}duisent :
625 \[
626 \left\lbrace
627 \begin{array}{lcl}
628 x_{V'}&=&(k+2)x_I-x_V\\
629 y_{V'}&=&(k+2)y_I-y_V\\
630 z_{V'}&=&(k+2)z_I-z_V
631 \end{array}
632 \right.
633 \]
634 Il reste \`{a} trouver l'intersection de $(IV')$ avec le plan horizontal $z=0$.
635
636 \'Equation param\'{e}trique de $(IV')$, $M$ \'{e}tant un point courant : $\overrightarrow{MI}=\alpha\overrightarrow{V'I}$
637 \[
638 \left\lbrace
639 \begin{array}{lcl}
640 x_I-x&=&\alpha(x_I-x_{V'})\\
641 y_I-y&=&\alpha(y_I-y_{V'})\\
642 z_I-z&=&\alpha(z_I-z_{V'})
643 \end{array}
644 \right.
645 \]
646 $z=0\Longrightarrow \alpha=\dfrac{z_{I}}{z_I-z_{V'}}$.
647
648 En rempla\c{c}ant $\alpha$ par son expression, nous obtenons les coordonn\'{e}es du point $P'$ de l'objet anamorphique.
649 \[
650 \left\lbrace
651 \begin{array}{lcl}
652 x_{P'}&=&x_I-\alpha(x_I-x_{V'})\\
653 y_{P'}&=&y_I-\alpha(y_I-y_{V'})\\
654 z_{P'}&=&0
655 \end{array}
656 \right.
657 \]
658 Cette s\'{e}rie de calculs doit \^{e}tre appliqu\'{e}e \`{a} tous les points de l'image <<~normale~>> afin d'obtenir l'objet anamorphique (d\'{e}form\'{e}) dont le miroir <<~redressera~>> la forme.
659
660 \textbf{Remarque} : l'image doit se former du c\^{o}t\'{e} de l'observateur \`{a} l'int\'{e}rieur du miroir, plus pr\`{e}s du bord du miroir que du centre. Si on d\'{e}place le point $P$ vers $O$, il arrive un moment o\`{u} le rayon r\'{e}fl\'{e}chi part au-dessus de l'horizontale et ne rencontre plus le plan horizontal. L'anamorphose n'est plus possible et pour les calculs c'est le CRASH !
661
662 \newpage
663
664 \section{La perspective}
665
666 Dans le livre de Jurgis Baltru\v{s}a\"{\i}tis\footnote{ \textit{Anamorphoses : les perspectives d\'{e}prav\'{e}es} en livre de poche chez Flammarion.}, on trouve le principe de la <<~\textit{costruzione legittima}~>> avec un sch\'{e}ma de L\'{e}onard de Vinci (1492) et des sch\'{e}mas anamorphiques de Niceron (1658). Je cite page 58 :
667 \begin{quote}\itshape
668 <<~Rappelons en quelques mots quels ont \'{e}t\'{e} le proc\'{e}d\'{e}s utilis\'{e}s par les artistes pour l'ordonnancement de leurs tableaux en perspective normale. La premi\`{e}re ligne trac\'{e}e est celle de l'horizon \`{a} la hauteur de l'\oe{}il. Deux points y sont ensuite fix\'{e}s : au milieu le point principal vers o\`{u} convergent toutes les lignes droites parall\`{e}les qui s'\'{e}loignent en profondeur ; sur la m\^{e}me horizontale et \`{a} la m\^{e}me distance du point principal que l'\oe{}il, en face de la composition -- le point de distance, vers lequel convergent les diagonales.~>>
669 \end{quote}
670 \begin{center}
671 \begin{pspicture}(-5,-3)(5,14.25)
672 \psset{ua=2,F=10,D=4,type=perspective}
673 \psframe*[linecolor=BleuCiel](-5,10)(5,13)
674 \psframe*[linecolor=OrangePale](-5,2)(5,10)
675 \psline[arrowinset=0,arrowsize=2mm]{->}(0,14)
676 \psline[arrowinset=0,arrowsize=2mm]{->}(7,0)
677 \uput[0](7,0){$x$}
678 \uput[90](0,14){$y$}
679 \psline(-5,2)(5,2)
680 \pnode(4,10){F'}
681 \pnode(0,10){F}
682 \uput[ul](F){${F}$}
683 \uput[ur](F'){${F'}$}
684 \rput(0,12){\multido{\n=22.5+45.0}{8}{%
685 \psline[linecolor=yellow](1;\n)}%
686 \pscircle[linestyle=none,fillstyle=gradient,gradmidpoint=0,gradend=yellow,gradbegin=gray]{0.5}}
687 \multido{\i=-2+1}{5}{%
688 \pnode(! \i\space -2){A}
689 \pnode(! \i\space 2){B}
690 \psline(A)(B)
691 \pslineA(A)(B)
692 }%
693 \multido{\i=-2+1}{5}{%
694 \pnode(!-2 \i){A}
695 \pnode(!2 \i){B}
696 \pslineA[linecolor=blue](A)(B)
697 \psline[linecolor=blue](A)(B)
698 }%
699 \pnode(0,0){O}
700 \pnodeA(0,0){O'}
701 \pnode(-2,2){A}
702 \pnodeA(-2,2){A'}
703 \pnode(2,2){B}
704 \pnodeA(2,2){B'}
705 \pnode(-2,-2){C}
706 \pnodeA(-2,-2){C'}
707 \pnode(2,-2){D}
708 \pnodeA(2,-2){D'}
709 \pnode(-1,-1){M1}
710 \pnodeA(-1,-1){M1'}
711 \pnode(-1,1){M2}
712 \pnodeA(-1,1){M2'}
713 \pnode(-1,0){N1}
714 \pnodeA(-1,0){N1'}
715 \pnode(-1,2){P}
716 \pnodeA(-1,2){P'}
717 \pnode(-2,0){N2}
718 \pnodeA(-2,0){N2'}
719 \pnode(-2,1){N3}
720 \pnodeA(-2,1){N3'}
721 \pnodeA(-2,-1){S'}
722 \pnode(0,2){Q}
723 \pnode(1,2){R}
724 \psline(A)(F)(B)
725 \psline[linecolor=red](A)(F')
726 \psline[linecolor=red](P)(F)
727 \psline[linecolor=lightgray](P)(F')
728 \psline[linecolor=lightgray](Q)(F')
729 \psline[linecolor=lightgray](R)(F')
730 \psline[linecolor=lightgray](S')(F')
731 \psline[linecolor=lightgray](N2')(F')
732 \psline[linecolor=lightgray](N3')(F')
733 \uput[dl](A){${A}$}
734 \uput[dr](B){${B}$}
735 \uput[ul](A'){${A'}$}
736 \uput[ur](B'){${B'}$}
737 \uput[dl](C){${C}$}
738 \uput[dr](D){${D}$}
739 \uput[ul](P){${P}$}
740 \uput[u](P'){${P'}$}
741 \uput[ul](C'){${C'}$}
742 \uput[ur](D'){${D'}$}
743 \uput[dl](O){${O}$}
744 \uput[dr](O'){${O'}$}
745 \psline{<->|}(-3,2)(-3,10)
746 \uput[l](-3,6){$f$}
747 \pcline{|<->|}(F)(F')
748 \uput[u](2,10){$e$}
749 \psdots(A)(B)(C)(D)(A')(B')(C')(D')(P)(P')(F)(F')(O)(O')
750 \psdots[linecolor=blue](N1)(N1')
751 \uput[dl](N1){$\blue {N_1}$}
752 \uput[dl](N1'){$\blue {N_1'}$}
753 \psdots[linecolor=gray](M2)(M2')
754 \uput[dr](M2){\tiny $\gray {(X,Y)}$}
755 \uput[dr](M2'){\tiny $\gray {(x',y')}$}
756 \psdots[linecolor=red](M1')
757 \uput[dr](M1'){\tiny $\red {(\alpha_1,\beta_1)}$}
758 \end{pspicture}
759 \end{center}
760
761 \newpage
762
763 Exemples :
764 \begin{itemize}
765 \item $A\longrightarrow A'$
766 \item $B\longrightarrow B'$
767 \item $C\longrightarrow C'$
768 \item $D\longrightarrow D'$
769 \item $O\longrightarrow O'$
770 \item $M_1\longrightarrow M_1'$
771 \item $M_2\longrightarrow M'_2$
772 \end{itemize}
773 D\'{e}terminons les coordonn\'{e}es $(\alpha_1,\beta_1)$ de l'intersection de $(PF)$ avec $(AF')$.
774
775 Posons que les coordonn\'{e}es des points essentiels sont :
776 \begin{itemize}
777 \item $F(0,f)$
778 \item $F'(e,f)$
779 \item $A(-a,a)$
780 \item $B(a,a)$
781 \item $C(a,-a)$
782 \item $D(-a,-a)$
783 \item $P(X,a)$
784 \end{itemize}
785 \'Equation de $(AF')$ :
786 \[
787 \frac{y-f}{x-e}=\frac{a-f}{-a-e}\Longrightarrow x(a-f)+y(a+e)-a(f+e)=0
788 \]
789 \'Equation de $(PF)$ :
790 \[
791 \frac{x-0}{y-f}=\frac{X-0}{a-f}\Longrightarrow x(a-f)-yX+fX=0
792 \]
793 Intersection $(PF)\bigcap (AF')$
794 \[
795 \alpha_1=\frac{Xe}{X+a+e}\qquad\beta_1=\frac{a(f+e)+fX}{X+a+e}
796 \]
797 Si on prend maintenant, un point d'ordonn\'{e}e $Y\neq X$ par exemple $N_1$ dont l'image $N'_1$ se situe toujours sur $(PF)$, mais \`{a} l'intersection de $PF$ avec la parall\`{e}le \`{a} $x'Ox$ men\'{e}e par le point-image du point de coordonn\'{e}e $(Y,Y)$ (ici $O'$ qui est l'image de $O(0,0)$).
798
799 Il s'agit de d\'{e}terminer l'intersection de $(PF)$ avec la droite d'\'{e}quation :
800 \[
801 y=\beta_2=\frac{a(f+e)+fY}{Y+a+e}
802 \]
803 Apr\`{e}s calculs et simplifications, on trouve pour l'abscisse :
804 \[
805 \alpha_2=\frac{Xe}{Y+a+e}
806 \]
807 En r\'{e}sum\'{e} si dans le rep\`{e}re $Oxy$, on appelle $({\red X},{\red Y})$ les coordonn\'{e}es d'un point-objet et $({\blue x'},{\blue y'})$ les coordonn\'{e}es du point image dans la transformation \textit{anamorphose oblique} ou \textit{perspective}, les formules qui permettent de passer de l'objet \`{a} l'image s'\'{e}crivent :
808 \[
809 \left\lbrace
810 \begin{array}{lcl}
811 {\blue x'}&=&\displaystyle\frac{{\red X}e}{{\red Y}+a+e}\\[0.5cm]
812 {\blue y'}&=&\displaystyle\frac{a(f+e)+f{\red Y}}{{\red Y}+a+e}
813 \end{array}
814 \right.
815 \]
816 \end{document}

Licence Creative Commons Les fichiers de Syracuse sont mis à disposition (sauf mention contraire) selon les termes de la
Licence Creative Commons Attribution - Pas d’Utilisation Commerciale - Partage dans les Mêmes Conditions 4.0 International.