椭圆弧旋转后计算角度问题

lzz0517

嗯，领导说不行，必须要这样搞，所以才这么麻烦，如果没有块的话，什么问题也没有了！呵呵，只是没办法啊！

highflybird

实际上你的这个椭圆已经很接近圆了，不如干脆做圆弧处理不更简单吗？没必要做椭圆处理啊？

lzz0517

这个只是为了更通用，有些椭圆弧接近圆弧，但有些不那么接近圆弧，所以就要这样计算，就不会出现画出来的跟CAD上显示的不一致！

highflybird

关于椭圆的OCS是这样计算的：
设椭圆中心坐标为cen，长轴矢量是maj,法线矢量是nrm --这里是(0 0 -1)
那么你这个题目中由长轴矢量 的单位矢量 XDir是：(0.134264 0.990946 0.0)
那么短轴的单位矢量是： nrm跟 XDir的叉积：  (0.990946 -0.134264 0.0)
因而这个椭圆自身的OCS变换矩阵mat是：
|0.134264 0.990946  0.0|
|0.990946 -0.134264 0.0|
|0        0         -1.|
注意了，在计算的时候矩阵要倒置一下。即：
|XDir.x  YDir.x  ZDir.x|
|XDir.y  yDir.y  ZDir.y|
|XDir.z  yDir.z  ZDir.z|
因而对于椭圆的起始角度可以算出：
由初始点： (340.122 * cos(2.35439) 340.122*0.9999999999977266*sin(2.35439) 0) = (-240.069 240.936 0)
进行上面的矩阵变换得到新的点。
(mat x (-240.069 240.936 0)) = (206.522 -270.244 0.0)
最后加上椭圆的中心，相当于平移: (3217176.6902 532617.0907 0)
这样就得到坐标的端点。
(3217383.2119 532346.8465 0)
这个坐标是块内坐标，
然后把这个坐标通过图块的变换矩阵变换到WCS(世界坐标系)这样就得到解了。
对于椭圆的长轴矢量，也可以由椭圆的中心加上长轴矢量，就是长轴端点，然后把椭圆的中心和长轴端点的块内坐标
再次乘图块的变换矩阵，分别得到世界坐标系统的椭圆中心和长轴端点坐标，把WCS下的长轴端点坐标减去椭圆中心坐标
这样就得到椭圆的WCS下的长轴矢量。
不知道这样说你是否理解了？

highflybird

另外要说的一点，一个图块炸开后和没炸开后，里面的东西有可能不一样的。譬如你这个图就是这样。
没炸开前的椭圆的参数和炸开后的椭圆的参数不同。

lzz0517

如果要求出块炸开的椭圆弧的起始角度和终止角度，没有长轴矢量坐标的话，好像没别的办法可以求出角度了，因为我用炸开块后的椭圆弧的参数进行求角度，得到的角度就是CAD上面显示的角度，用的办法如下：http://blog.csdn.net/ymain/article/details/6225722；其中的参数为长轴矢量坐标，起点矢量坐标，这两个参数代入getRotateAngle函数求出起始角度，长轴矢量坐标，端点矢量坐标，这两个参数代入getRotateAngle函数求出端点角度。

问题是现在不知道怎么求出这个长轴矢量坐标！

如果用余弦定理求夹角，也还是要知道起点或端点与圆心的连线跟圆心左边的长轴端点坐标的夹角还是右边的夹角，所以要求角度总绕不开长轴矢量坐标！

highflybird

lzz0517 发表于 2012-6-6 10:52
如果要求出块炸开的椭圆弧的起始角度和终止角度，没有长轴矢量坐标的话，好像没别的办法可以求出角度了，因 ...

我现在已经在70楼做了修改。请回看。

highflybird

另外，如果不是用CAD编程里面的函数，那么你用C#如何创建一个椭圆的？
能否把函数原型写出来？或者这段代码写出来？

lzz0517

我用来画倾斜椭圆弧函数如下：
void drawellipsearc(double x,double y,double sa,double ea,double rotateangle,double R,double r)
{
    if(sa>ea)
        ea+=2*PI;
   List<PointF> listpoint=new List<PointF>();
   for(int i=0;i<65;i++)
{
    double angle=sa+i*(ea-sa)/64;
   double x=R*cos(angle);
   double y=r*sin(angle);
   double nx=x*cos(rotateangle)-y*sin(rotateangle)+x;
   double ny=y*cos(rotateangle)+x*sin(rotateangle)+y;
   PointF p=new PointF((float)nx,(float)ny);
   listpoint.add(p);
}
listpoint.add(listpoint[0]);
Pen pen = new Pen(Color.Black, 3);
Graphics.DrawLines(pen, listpoint.ToArray());}

lzz0517

highflybird 发表于 2012-6-6 10:59
我现在已经在70楼做了修改。请回看。

你说那个是不是用来求长轴矢量的呢，我今天看了一下没看明白怎么用！