如何根据四个旋转和缩放的点绘制椭圆?

时间:2015-09-04 14:53:53

标签: c++ math geometry ellipse drawellipse

我从一个完美的圆圈开始,然后我需要旋转并自由缩放它而不会松散椭圆特征:

Circle distortion

我想知道是否有办法描述一个像仅基于余弦和正弦的扭曲椭圆。

我有这段代码:

    float centerX = 100;
    float centerY = 100;
    float radiusX = 100;
    float radiusY = 100;
    float rotation = PI;
    float res = 30;
    for (int i = 0; i < res; i++) {

        //find the point in space for the circle resolution
            float angle = (float)i/res * TWO_PI;

            float x = radiusX * cos(angle);
            float y = -radiusY * sin(angle);

        //rotate the point
            float s = sin(rotation);
            float c = cos(rotation);

            // translate point to origin:
            p->x -= centerX;
            p->y -= centerY;

            float xnew = p->x * c - p->y * s;
            float ynew = p->x * s + p->y * c;

            // translate point back
            x = xnew + centerX;
            y = ynew + centerY;

        //apply X Y
    }

此代码只能表示忽略旋转和比例之间关系的相同椭圆:

enter image description here

1 个答案:

答案 0 :(得分:0)

我使用的方法与我多年前在这里提出的方法相同,但是那个时候扭曲一个方块是有用的: What's the correct way to draw a distorted plane in OpenGL?

以下是代码:

//setup the four coordinates
float topX = 20;
float topY = 10;

float bottomX = 30;
float bottomY = 20;

float rightX = 30;
float rightY = 20;

float leftX = 30;
float leftY = 20;

//calculate the horizontal radius
float distHorX = rightX-leftX;
float distHorY = rightY-leftY;
float radiusX = sqrt(distHorX * distHorX + distHorY * distHorY)/2.f;

//calculate the vertical radius
float distVerX = topX-bottomX;
float distVerY = topY-bottomY;
radiusY = sqrt(distVerX * distVerX + distVerY * distVerY)/2.f;

float res = 30;
for (int i = 0; i < res; i++) {

        float angle = (float)i/res * TWO_PI;

        float x = radiusX * cos(angle);
        float y = -radiusY * sin(angle);

        //corvert the circle inside a square to a square inside a circle
        //it is a magical number I have found to convert that proportion
        x /= 0.705069124;
        y /= 0.705069124;

        //transform the points based on that four coordinates
        float pctx = (x + radiusX) / (radiusX*2);
        float pcty = (y + radiusY) / (radiusY*2);

        float linePt0x = (1-pcty)* topX + pcty * leftX;
        float linePt0y = (1-pcty)* topY + pcty * leftY;
        float linePt1x = (1-pcty)* rightX + pcty * bottomX;
        float linePt1y = (1-pcty)* rightY + pcty * bottomY;
        float ptx      = (1-pctx) * linePt0x + pctx * linePt1x;
        float pty      = (1-pctx) * linePt0y + pctx * linePt1y;

        //apply X Y
        x = ptx;
        y = pty;
    }