OldClaudio” post=98957Visto che l’appetito vie mangiando, saresti anche capace di disegnare una lossodromia fra due punti dati?
Forse siamo un po’ OT e per questo chiedo scusa a @Squall. Questo regalo è tutto per @OldClaudio.
Allo stesso modo dell’ortodromia ho definito un pic [tt]lossodromia[/tt] che genera la curva lossodromica tra due punti assegnati sulla sfera. Il codice usa solo i comandi di TikZ, ma per disegnare la sfera mi sono appoggiato a pgfplots, che non è necessario. Il pic utilizza l’istruzione [tt]plot[/tt] di TikZ con i valori di default (25 campioni), ma è possibile passare al comando le opzioni [tt]smooth[/tt], [tt]hobby[/tt] oppure [tt]samples=
Il cuore è il comando:
`\pic [partenza={latitudine A}{longitudine A}, arrivo={latitudine B}{longitudine B}] {lossodromia};`
nel quale, come di consueto, la latitudine 0 corrisponde al piano [xy] e la longitudine 0 al semipiano [xz] di ascisse positive. Non è possibile scegliere come punto di partenza o di arrivo uno dei due poli, ma è possibile prendere un punto di latitudine molto vicina (ad esempio 89.98°). Non è nemmeno possibile prendere due punti sullo stesso meridiano, ma in tal caso la lossodromia coincide con l’ortodromia e si può usare quel codice. Volendo si può pensare di affinare il codice per includere questi casi particolari.
Anche qui è possibile utilizzare la chiave [tt]raggio[/tt] che è un alias per la chiave [tt]scale[/tt] e accetta solo numeri adimensionali che vengono interpretati come fattore di scala della lunghezza del versore dell’asse delle ascisse.
`\documentclass
{standalone}
\usepackage{pgfplots}
\pgfplotsset{compat=1.11}
% componenti dei versori degli assi coordinati
% rispetto alla base della tela
\pgfmathsetlengthmacro{\xx}{-2cm*sqrt(3)/2}
\pgfmathsetlengthmacro{\xy}{-1cm}
\pgfmathsetlengthmacro{\yx}{2cm*sqrt(3)/2}
\pgfmathsetlengthmacro{\yy}{-1cm}
\pgfmathsetlengthmacro{\zx}{0cm}
\pgfmathsetlengthmacro{\zy}{2cm}
\tikzset{
latA/.initial=0,
lonA/.initial=0,
latB/.initial=0,
lonB/.initial=0,
partenza/.style 2 args={latA=#1, lonA=#2},
arrivo/.style 2 args={latB=#1, lonB=#2},
posizione/.style 2 args={latA=#1, lonA=#2},
raggio/.initial=1,
raggio/.code={\tikzset{/tikz/scale=#1}},
declare function={
s(\x) = ln( (1+tan(\x/2)) / (1-tan(\x/2))) );
lat(\x) = 2*atan( tanh( (\m*(\x-\lonA) + s(\latA))/2) ) );
},
lossodromia/.pic={
% riferimento http://homepages-fb.thm.de/boergens/marken/marke062/marke062_loxo.pdf
\pgfmathsetmacro{\latA}{\pgfkeysvalueof{/tikz/latA}}
\pgfmathsetmacro{\lonA}{\pgfkeysvalueof{/tikz/lonA}}
\pgfmathsetmacro{\latB}{\pgfkeysvalueof{/tikz/latB}}
\pgfmathsetmacro{\lonB}{\pgfkeysvalueof{/tikz/lonB}}
\pgfmathsetmacro{\m}{(s(\latB)-s(\latA)) / (\lonB-\lonA)}
\draw [variable=\t, domain=\lonA:\lonB, pic actions]
plot (
{abs(\t)!=90 ? cos(\t)*cos(lat(\t)) : 0},
{abs(\t)!=90 ? sin(\t)*cos(lat(\t)) : 0},
{abs(\t)!=90 ? sin(lat(\t)) : sin(\t)});
},
coordinata/.pic={
\pgftransformshift{
\pgfpointspherical{90-\pgfkeysvalueof{/tikz/lonA}}{\pgfkeysvalueof{/tikz/latA}}{1}
}
\pgfnode{coordinate}{center}{}{#1}{}
}
}
\begin{document}
\begin{tikzpicture}
\begin{axis} [x={(\xx,\xy)}, y={(\yx,\yy)}, z={(\zx,\zy)}, hide axis]
% disegna i meridiani e i paralleli
\addplot3 [surf, z buffer=sort, fill opacity=0.6, white, faceted color=blue!40,
samples=19, samples y=37, variable=\u, variable y=\v,
domain=0:180, y domain=0:360]
({cos(u)*sin(v)}, {sin(u)*sin(v)}, {cos(v)});
% disegna il contorno della sfera
\draw [blue!40] (0,0,0) circle ({2cm/cos(atan(1/sqrt(2)))});
% disegna la lossodromia
\pic [partenza={-10}{20}, arrivo={89.98}{360}, smooth, samples=50, red] {lossodromia};
% definisce le coordinate dei punti di partenza e arrivo
\pic [posizione={-10}{20}, x={(\xx,\xy)}, y={(\yx,\yy)}, z={(\zx,\zy)}] {coordinata=A};
\pic [posizione={89.98}{360}, x={(\xx,\xy)}, y={(\yx,\yy)}, z={(\zx,\zy)}] {coordinata=B};
% traccia i pallini in corrispondenza degli estremi della curva
\fill (A) circle (1pt)
(B) circle (1pt);
% disegna gli assi coordinati
\draw [-latex] (1,0,0) — (1.5,0,0);
\draw [-latex] (0,1,0) — (0,1.5,0);
\draw [-latex] (0,0,1) — (0,0,1.5);
\end{axis}
\end{tikzpicture}
\end{document}`
[attachment]C:\fakepath\lossodromia.png[/attachment]
Il codice, come anche quello dell’ortodromia, al momento presenta un bug: quando usato all’interno dell’ambiente [tt]axis[/tt] di pgfplots per disegnare la sfera passando come opzioni i versori degli assi coordinati, i pallini (e nel caso dell’ortodromia la curva), potrebbero essere tracciati fuori asse per alcune combinazioni di assi coordinati. Il problema nasce dal fatto che il sistema di riferimento viene specificato due volte: una per l’ambiente [tt]axis[/tt] e l’altra per i pic. Se qualcuno mi aiuta a risolverlo, gli sono grato.
Ciao