У меня еще пара вопрсов 
1. При использовании разбиения на треугольники на некоторых границах возникают пробелы.
На приведенной картинке видно эти артефакты в областях искажения. С начало я во всех циклах использовал только целые числа. При таком подходе возникали горизонтальные искажения. Не пробелы, а именно искажения, т.е. рисовались не те цвета. Когда же я заменил все циклы на вещественные, то получил разрывы текстур на некоторых границах. В чем может быть загвоздка?
Вообще пишут, что "При реализации алгоритма удобно рассчитывать границы по какому-либо алгоритму растеризации отрезка". Если это может исправить данную ситуацию, то можно пример, пожалуйста
2. Можно по подробнее о сглаживании Просто при больших искажениях получается достаточно видимый алиасинг (идет большая пикселизация).
Вот код, который я накатал (заранее извиняюсь ):
| Код |
/*
* OffsetMesh - класс вершин, имеет методы GetNvY() - количество вершин по Y, GetNvX() - количество вершин по X,
* GetVretex(i, j, vex) - возвращает вершину по ее индексам в массиве.
*/
void DistortPictureWithTriangle(IN OffsetMesh* Mesh, IN Gdiplus::Bitmap* Texture, OUT Gdiplus::Bitmap* Picture)
{
#define GRID_SIZE 9 // постоянная сетки
int i1, j1, i2, j2;
PointF vx1, vx2, vx3; // вершины треугольника, которые содержат относительные смещения (относительно прямоугольной сетки)
PointF p12, p22, p32; //
PointF p11, p21, p31;
Graphics grph(Picture);
grph.Clear(Color(100, 100, 100));
// for each vertex in mesh (для каждой вершины из массива)
for (int j = 0; j < Mesh->GetNvY() - 1; j++)
{
for (int i = 0; i < Mesh->GetNvX() - 1; i++)
{
/************************************************************************/
/* Triangle 1
/*
/* p2 *---* p3
/* | /
/* | /
/* |/
/* * p1
/*
/************************************************************************/
// извлекаем вершины треугольника из массива
Mesh->GetVretex(i, j + 1, &vx1);
Mesh->GetVretex(i, j, &vx2);
Mesh->GetVretex(i + 1, j, &vx3);
// вычисляем реальное расположение вершин на картинке (т.к. величина ячейки сетки константа, то это сделать легко)
i1 = i * GRID_SIZE;
i2 = (i + 1) * GRID_SIZE;
j1 = j * GRID_SIZE;
j2 = (j + 1) * GRID_SIZE;
// undeformed points
p11.X = (REAL)i1;
p11.Y = (REAL)j2;
p21.X = (REAL)i1;
p21.Y = (REAL)j1;
p31.X = (REAL)i2;
p31.Y = (REAL)j1;
// deformed points
p12.X = (REAL)(i1 + vx1.x);
p12.Y = (REAL)(j2 + vx1.y);
p22.X = (REAL)(i1 + vx2.x);
p22.Y = (REAL)(j1 + vx2.y);
p32.X = (REAL)(i2 + vx3.x);
p32.Y = (REAL)(j1 + vx3.y);
DrawTriangle(p11, p21, p31, p12, p22, p32);
/************************************************************************/
/* Triangle 2
/*
/* * p2
/* /|
/* / |
/* / |
/* p1 *---* p3
/*
/************************************************************************/
// извлекаем вершины треугольника из массива
m_Mesh->GetVretex(i, j + 1, &vx1);
m_Mesh->GetVretex(i + 1, j, &vx2);
m_Mesh->GetVretex(i + 1, j + 1, &vx3);
i1 = i * GRID_SIZE;
i2 = (i + 1) * GRID_SIZE;
j1 = j * GRID_SIZE;
j2 = (j + 1) * GRID_SIZE;
// undeformed points
p11.X = (REAL)i1;
p11.Y = (REAL)j2;
p21.X = (REAL)i2;
p21.Y = (REAL)j1;
p31.X = (REAL)i2;
p31.Y = (REAL)j2;
// deformed points
p12.X = (REAL)(i1 + vx1.x);
p12.Y = (REAL)(j2 + vx1.y);
p22.X = (REAL)(i2 + vx2.x);
p22.Y = (REAL)(j1 + vx2.y);
p32.X = (REAL)(i2 + vx3.x);
p32.Y = (REAL)(j2 + vx3.y);
DrawTriangle(p11, p21, p31, p12, p22, p32);
}
}
}
void DrawTriangle(Gdiplus::Bitmap* Texture, Gdiplus::Bitmap* Picture, PointF p11, PointF p21, PointF p31, PointF p12, PointF p22, PointF p32)
{
REAL b11, b12;
REAL b21, b22;
REAL a11, a12;
REAL a21, a22;
REAL x1, y1, x2, y2, x2min, x2max;
REAL Det, UnDet;
PointF *p2min, *p2mid, *p2max;
int width, height;
Gdiplus::Color pix;
// save picture size
width = Texture->GetWidth();
height = Texture->GetHeight();
// coefficients (x2,y2 -> u,v)
b11 = p22.X - p12.X;
b12 = p32.X - p12.X;
b21 = p22.Y - p12.Y;
b22 = p32.Y - p12.Y;
// coefficients (u,v -> x1,y1)
a11 = p21.X - p11.X;
a12 = p31.X - p11.X;
a21 = p21.Y - p11.Y;
a22 = p31.Y - p11.Y;
// Det
Det = b11 * b22 - b12 * b21;
UnDet = 1 /Det;
// calculate scan bounds
p2min = p12.Y <= p22.Y ? (p12.Y <= p32.Y ? &p12 : &p32) : (p22.Y <= p32.Y ? &p22 : &p32);
p2max = p12.Y > p22.Y ? (p12.Y > p32.Y ? &p12 : &p32) : (p22.Y > p32.Y ? &p22 : &p32);
if ((p2max != &p12) && (p2min != &p12)) p2mid = &p12;
if ((p2max != &p22) && (p2min != &p22)) p2mid = &p22;
if ((p2max != &p32) && (p2min != &p32)) p2mid = &p32;
if (p2min->Y != p2mid->Y)
for (y2 = p2min->Y; y2 <= p2mid->Y; y2++)
{
x2min = (p2min->X + (p2mid->X - p2min->X) / (p2mid->Y - p2min->Y) * (y2 - p2min->Y));
x2max = (p2min->X + (p2max->X - p2min->X) / (p2max->Y - p2min->Y) * (y2 - p2min->Y));
if (x2min > x2max)
{
for (x2 = x2min; x2 >= x2max; x2--)
{
x1 = UnDet * ((a11 * b22 - a12 * b21) * (x2 - p12.X) + (a12 * b11 - a11 * b12) * (y2 - p12.Y)) + p11.X;
y1 = UnDet * ((a21 * b22 - a22 * b21) * (x2 - p12.X) + (a22 * b11 - a21 * b12) * (y2 - p12.Y)) + p11.Y;
// bounds checking
if (x1 < 0 || x1 > width || y1 < 0 || y1 > height ||
x2 < 0 || x2 > width || y2 < 0 || y2 > height) continue;
Texture->GetPixel((int)(x1), (int)(y1), &pix);
Picture->SetPixel((int)(x2), (int)(y2), pix);
}
}
else
{
for (x2 = x2min; x2 <= x2max; x2++)
{
x1 = UnDet * ((a11 * b22 - a12 * b21) * (x2 - p12.X) + (a12 * b11 - a11 * b12) * (y2 - p12.Y)) + p11.X;
y1 = UnDet * ((a21 * b22 - a22 * b21) * (x2 - p12.X) + (a22 * b11 - a21 * b12) * (y2 - p12.Y)) + p11.Y;
// bounds checking
if (x1 < 0 || x1 > width || y1 < 0 || y1 > height ||
x2 < 0 || x2 > width || y2 < 0 || y2 > height) continue;
Texture->GetPixel((int)(x1), (int)(y1), &pix);
Picture>SetPixel((int)(x2), (int)(y2), pix);
}
}
}
if (p2max->Y != p2mid->Y)
for (y2 = p2mid->Y; y2 <= p2max->Y; y2++)
{
x2min = (p2mid->X + (p2max->X - p2mid->X) / (p2max->Y - p2mid->Y) * (y2 - p2mid->Y));
x2max = (p2min->X + (p2max->X - p2min->X) / (p2max->Y - p2min->Y) * (y2 - p2min->Y));
if (x2min > x2max)
{
for (x2 = x2min; x2 >= x2max; x2--)
{
x1 = UnDet * ((a11 * b22 - a12 * b21) * (x2 - p12.X) + (a12 * b11 - a11 * b12) * (y2 - p12.Y)) + p11.X;
y1 = UnDet * ((a21 * b22 - a22 * b21) * (x2 - p12.X) + (a22 * b11 - a21 * b12) * (y2 - p12.Y)) + p11.Y;
// bounds checking
if (x1 < 0 || x1 > width || y1 < 0 || y1 > height ||
x2 < 0 || x2 > width || y2 < 0 || y2 > height) continue;
Texture->GetPixel((int)(x1), (int)(y1), &pix);
Picture->SetPixel((int)(x2), (int)(y2), pix);
}
}
else
{
for (x2 = x2min; x2 <= x2max; x2++)
{
x1 = UnDet * ((a11 * b22 - a12 * b21) * (x2 - p12.X) + (a12 * b11 - a11 * b12) * (y2 - p12.Y)) + p11.X;
y1 = UnDet * ((a21 * b22 - a22 * b21) * (x2 - p12.X) + (a22 * b11 - a21 * b12) * (y2 - p12.Y)) + p11.Y;
// bounds checking
if (x1 < 0 || x1 > width || y1 < 0 || y1 > height ||
x2 < 0 || x2 > width || y2 < 0 || y2 > height) continue;
Texture->GetPixel((int)(x1), (int)(y1), &pix);
Picture->SetPixel((int)(x2), (int)(y2), pix);
}
}
}
}
|
|
-> 
,