![Intersectionpoints 3D => path3 与贝塞尔三角形的交点尚未实现](https://linux22.com/image/449832/Intersectionpoints%203D%20%3D%3E%20path3%20%E4%B8%8E%E8%B4%9D%E5%A1%9E%E5%B0%94%E4%B8%89%E8%A7%92%E5%BD%A2%E7%9A%84%E4%BA%A4%E7%82%B9%E5%B0%9A%E6%9C%AA%E5%AE%9E%E7%8E%B0.png)
我在 W10 PC 上使用 Asymptote 2.80 版。我尝试重新使用旧的(2019 年)代码来绘制球体与圆锥体的交点。2019 年生成的绘图已合并。?
我在 2019 年就这样做过,但现在编译时出现以下错误消息:
if(s.triangular) abort(triangular); ^ C:\Program Files\Asymptote2_80/three_surface.asy: 1296.25: path3 与贝塞尔三角形的交点尚未实现。
代码如下——绘制了所有的构造步骤。
\begin{figure}[h]
\centering
\begin{asy}
size(10cm);
import graph3;
currentprojection=orthographic((5,1,2),up=Z);
settings.render=4; // Pose pb si =0 : Le plan devant la sphere est absent
settings.prc=false;
currentlight=Viewport;
real k=1;
path3 Directrice = rotate(30,Y)*unitcircle3;
draw(Directrice,orange);
triple S=(0,-1.5,2);
draw(Label ("S",green),S);
path3 Directrice = rotate(60,-Y)*unitcircle3;
path3 Directrice = shift(1.5,0,0)*Directrice;
triple P1 = relpoint(Directrice,0);
triple P2 = relpoint(Directrice,0.25);
triple P3 = relpoint(Directrice,0.5);
triple P4 = relpoint(Directrice,0.75);
// t.x => point m de D -- t.y => facteur k
triple f(pair t) {
triple Point =relpoint(Directrice,t.x);
real X= S.x+t.y*(Point.x-S.x);
real Y= S.y+t.y*(Point.y-S.y);
real Z= S.z+t.y*(Point.z-S.z);
//dot((X,Y,Z),orange);
return (X,Y,Z);
}
surface s1=surface(f,(0,0),(1,1),36,36,Spline);
draw(s1,orange+opacity(.7));
surface sph=shift(0,0,0)*scale3(1)*unitsphere;
draw(sph,blue);
limits((0,0,0),(k+2,k+2,k+2));
xaxis3("$x$",Arrow3);
yaxis3("$y$",Arrow3);
zaxis3("$z$",Arrow3);
draw(Directrice,red);
dot(Directrice^^S,orange);
// size(p) = C'est le nombre de points qui définissent le path p
for(int i=0; i<size(Directrice);++i) {
draw(point(Directrice,i)--S,dashed);
}
/////// INTERSECTION -- On recherche les points d'intersection de 100 segments [Sm] avec la sphère.
triple[][] inter;
path3 ch1,ch2;
int n=100;
for(int i=0; i<n; ++i){
path3 p=S--relpoint(Directrice,i/n);
inter.push(intersectionpoints(p,sph,0.001));
draw(p,yellow);
}
for(int i=0; i<inter.length; ++i){
if (inter[i].length!=0){
ch1=ch1..inter[i][0];
dot(inter[i][0],2bp+green);
}
if(inter[i].length==2){
ch2=ch2..inter[i][1];
dot(inter[i][1],2bp+red);}
}
draw(ch1..cycle,1bp+green);
draw(ch2..cycle,1bp+red);
// Il reste à fermer la boucle : du 1er point de construction de ch1 au 1er point de construction ce ch2, et idem pour le dernier point
draw(point(ch1,0)..point(ch2,0),5bp+orange );
draw(point(ch1,(size(ch1)-1) )..point(ch2,(size(ch2)-1) ) ,5bp+orange );
\end{asy}
\caption{Projection gnomonique 3D \label{Fig_Proj_Gnomo_3D} }
\end{figure}
我可以要求重新调整哪些交点(什么与什么之间)
这个想法是先在与北极相切的平面上画出地面平行线的球面投影,然后再在与赤道相切的平面上画出,最后在北纬 50 度的平面上画出切线:欢迎任何不使用“交点”的解决方案