Bezjeove krive
Bezjeove krive su parametarske krive. Njihovu najbolju primenu pronalazimo u računarskoj grafici.
Zadaju se određenim brojem kontrolnih tačaka čijim se pomeranjem menja oblik krive. Dodavanjem kontrolnih tačaka se povećava stepen krive, ali i složenost izračunavanja, zbog čega se najčešće koriste krive malog stepena. Takozvanim „lepljenjem“ krivih malog stepena, dobijaju se Bezjeovi splajnovi (složene krive).
Pol Kastelžo i Pjer Bezje su ih prvi koristili pri dizajnu automobila krajem pedesetih godina dvadesetog veka[1].
Definicija Bezjeove krive
[uredi | uredi izvor]Rekurzivna definicija
[uredi | uredi izvor]Ako Bezjeovu krivu određenu pomoću tački označimo kao . Onda:
- i
Eksplicitna definicija
[uredi | uredi izvor]Neka je Bezjeova kriva, stepena , zadata pomoću tačaka ravni. Možemo je predstaviti sledećim polinomom:
Tačke pomoću kojih je definisana kriva nazivaju se kontrolne tačke, a polinomi Bernstajnovi polinomi, tj. bazne funkcije.
Tačke obrazuju kontrolni poligon.
Polinomijalni oblik
[uredi | uredi izvor]Primenom binomne teoreme na definiciju krive dobijamo:
gde je:
Izvod
[uredi | uredi izvor]Izvod krive stepena je:
Specijalni slučajevi
[uredi | uredi izvor]Najveću praktičnu primenu imaju Bezjeove krive drugog i trećeg stepena.
Linearna Bezjeova kriva
[uredi | uredi izvor]Bezjeova kriva stepena 1, tj. kriva određena sa samo dve tačke i , jeste zapravo duž , zadata sledećom parametrizacijom:
Kvadratna Bezjeova kriva
[uredi | uredi izvor]
Kvadratna Bezjeova kriva zadata je sledećom parametrizacijom:
Bezjeova kriva drugog stepena je deo parabole. Dokaz leži u osobini invarijantnosti afinih preslikavanja.
Izvod:
Kubna Bezjeova kriva
[uredi | uredi izvor]Kriva je zadata pomoću tačaka . Počinje u , ide prema i stiže u dolazeći iz smera . Uglavnom ne prolazi kroz i ; ove tačke služe da obezbede dodatne informacije (razdaljina između i određuje „koliko daleko“ i „koliko brzo“ se kriva kreće prema pre nego što skrene ka ).
Kubna Bezjeovakriva zadata je sledećom paramtrizacijom:
Prvi izvod:
Drugi izvod:
Osobine Bezjeovih krivih
[uredi | uredi izvor]- Najveći stepen Bezjeove krive definisane sa kontrolnih tačaka je .
- Početna tačka krive je , tj.
- Vektor je kolinearan tangentnom vektoru krive u tački dok je vektor kolinearan tangentnom vektoru krive .
- Osobina nenegativnosti: za možemo videti da su svi Bernstajnovi polinomi , tj. nenegativni su.
- Osobina konveksnog omotača: Bezjeova kriva koja je određena pomoću kontrolnih tačka pripada konveksnom omotaču tih tačka. Centar mase tačaka sa pridruženim pozitivnim masama je tačka . Na slici je prikazana kriva 6. stepena, a konveksni omotač njenih 7 kontrolnih tačaka je označen roze bojom. Izuzev dve krajnje tačke koje su na rubu, čitava kriva leži unutar omotača.
Konveksni omotač Bezjeove krive | Osobina manje varijacije |
- Osobina manje varijacije: nijedna prava ne preseca Bezijerovu krivu više puta nego što seče kontrolnu poligonsku liniju.
- Osobina afine invarijantnosti: slika Bezjeove krive pri afinom preslikavanju je određena slikama njenih kontrolnih tačaka pri tom istom preslikavanju.
Bezjeova kriva na proizvoljnom segmentu
[uredi | uredi izvor]Osnovni domen Bezjeove krive je segment , ali ga možemo proširiti na proizvoljni segment. U tom slučaju je potrebno izvršiti smenu promenljive na sledeći način:
Sada je promenljiva definisana na intervalu i potrebno ju je zameniti u bazne funkcije kako bi kriva bila definisana na segmentu :
Pomeranje kontrolnih tačaka
[uredi | uredi izvor]Pretpostavimo da je kontrolna tačka pomerena za vektor .
Nova kriva je određena tačkama i njena jednačina je:
Nova kriva se predstavlja kao zbir početne krive i dodatnog člana Tačka na novoj krivoj je dobijena translacijom tačke na početnoj krivoj u pravcu na rastojanju . Osim krajnjih tačaka, sve tačke početne krive se pomeraju na nove pozicije.
Pomeranje bilo koje kontrolne tačke Bezjeove krive menja krivu globalno.
Početna kriva na slici je kriva 8. stepena sa 9 kontrolnih tačaka i označena je crnom bojom. Ako kontrolnu tačku pomerimo na novu poziciju, crna kriva se menja u crvenu.
Matrična reprezentacija[2]
[uredi | uredi izvor]Kvadratne Bezjeove krive
[uredi | uredi izvor]
Kubne Bezjeove krive
[uredi | uredi izvor]
Gde je
Odakle se analogijom dobija uopštenje:
Podela Bezjeove krive na dva dela korišćenjem matrične reprezentacije
[uredi | uredi izvor]Kvadratna Bezjeova kriva
[uredi | uredi izvor]Razmotrimo sad podelu krive na dva dela pravljenjem reza u tački .
Prva kriva:
,
gde je data sa:
Temena kontrolnog poligona dobijamo iz:
Za temena su data sa:
Druga kriva je definisana sa:
,
gde je:
.
Kontrolni poligon je određen sledećim tačkama:
.
tj. za :
Kubna Bezjeova kriva
[uredi | uredi izvor]Razmotrimo sad podelu krive na dva dela pravljenjem reza u tački .
Prva kriva:
,
gde je data sa:
Temena kontrolnog poligona dobijamo iz:
Za temena su data sa:
Druga kriva je definisana sa:
,
gde je:
.
Kontrolni poligon je određen sledećim tačkama:
tj. za :
De Kastelžoov algoritam
[uredi | uredi izvor]Ako želimo da odredimo tačku na krivoj, možemo jednostavno zameniti parametar u jednačinu krive. Ova metoda je numerički nestabilna i zato se koristi jednostavnija metoda: De Kastelžoov algoritam, čija je osnovna ideja u iterativnoj podeli kontrolnih duži u odnosu .
Podela Bezjeove krive
[uredi | uredi izvor]Smisao podele krive je da je „isečemo“ u nekoj tački, tako da ona ostane podeljena na dva dela od kojih je svaki i dalje Bezjeova kriva. Početni skup kontrolnih tačaka se odbacuje, a rezultujuće Bezjeove krive se zadaju pomoću svojih novih kontrolnih tačkama.
Algoritam podele bazira se na De Kastelžoovom algoritmu.
Podela krive ima veliki broj primena: dizajn krivih, renderovanje Bezjeove krive, određivanje preseka dve krive... Ako imamo neku krivu koja nam ne odgovara, možemo je podeliti na nezadovoljavajući i zadovoljavajući deo i onda se skoncentrisati na promenu nezadovoljavajućeg dela.
Kriva se može deliti proizvoljan broj puta, a ako želimo da se podeoni preseci krive glatko spajaju, onda granične tačke i njihove dve susedne tačke moraju biti kolinearne.
Isto to važi ako želimo spojiti dve krive u jednu. Kako je crtanje krivih većeg stepena zahtevno, mogu se iskoristiti krive manjeg stepena i onda ih glatko povezivati.
Povećanje stepena krive
[uredi | uredi izvor]Računarstvo
[uredi | uredi izvor]Kodiranje[3]
[uredi | uredi izvor]Neki programski jezici omogućavaju crtanje Bezjeovih krivih ako su poznate kontrolne tačke kojima su određene.
Na primer, u html-u napravimo canvas tag:
<canvas id="myCanvas" width="500" height="300"></canvas>
U okviru JavaScript-a definišemo canvas preko promenljivih:
var c=document.getElementById("myCanvas");
var ctx=c.getContext("2d");
Sada možemo koristiti sledeći kod koji iscrtava krivu stepena 3, gde i označavaju i koordinate -te tačke.
ctx.beginPath();
ctx.moveTo(<math>P_0x,P_0y</math>);
ctx.bezierCurveTo(<math>P_1x,P_1y,P_2x,P_2y,P_3x,P_3y</math>);
ctx.stroke();
Ako želimo krivu drugog stepena, potrebno je pozvati odgovarajuću funkciju, tj.:
ctx.beginPath();
ctx.moveTo(<math>P_0x,P_0y</math>);
ctx.quadraticCurveTo(<math>P_1x,P_1y,P_2x,P_2y</math>);
ctx.stroke();
Naredni kod je jednostavan i praktičan primer koji prikazuje kako se crta kubna Bezjeova kriva u programskom jeziku C. Ovaj algoritam je linearan jer je u praksi brži i lakši za implementaciju, za razliku od rekurzivnog De Kastelžoovog algoritma koji se često pominje u računarskoj grafici.
/******************************************************
Kod za generisanje kubne Bezjeove krive
*******************************************************/
typedef struct
{
float x;
float y;
}
Point2D;
/******************************************************
cp je niz koji sadrži 4 elementa gde:
cp[0] je početna tačka
cp[1] je prva kontrolna tačka
cp[2] je druga kontrolna tačka
cp[3] je krajnja tačka
t je parametar, 0 <= t <= 1
*******************************************************/
Point2D PointOnCubicBezier( Point2D* cp, float t )
{
float ax, bx, cx;
float ay, by, cy;
float tSquared, tCubed;
Point2D result;
/* računa koeficijente polinoma*/
cx = 3.0 * (cp[1].x - cp[0].x);
bx = 3.0 * (cp[2].x - cp[1].x) - cx;
ax = cp[3].x - cp[0].x - cx - bx;
cy = 3.0 * (cp[1].y - cp[0].y);
by = 3.0 * (cp[2].y - cp[1].y) - cy;
ay = cp[3].y - cp[0].y - cy - by;
/* računa tačku krive u zavisnosti od parametra t */
tSquared = t * t;
tCubed = tSquared * t;
result.x = (ax * tCubed) + (bx * tSquared) + (cx * t) + cp[0].x;
result.y = (ay * tCubed) + (by * tSquared) + (cy * t) + cp[0].y;
return result;
}
/*****************************************************************************
Argument funkcije ComputeBezier je niz struktura Point2D sa tačkama krive
generisanim uz pomoć kontrolnih tačaka cp. Pozivalac mora da obezbedi dovoljno
memorije za rezultat čija je veličina <sizeof(Point2D) * brojTačaka>
******************************************************************************/
void ComputeBezier( Point2D* cp, int numberOfPoints, Point2D* curve )
{
float dt;
int i;
dt = 1.0 / ( numberOfPoints - 1 );
for( i = 0; i < numberOfPoints; i++ )
curve[i] = PointOnCubicBezier( cp, i*dt );
}
Grafika
[uredi | uredi izvor]Bezjeove krive se naširoko koriste za modelovanje glatkih krivih. Afine transformacije kao što su translacija i rotacija mogu biti primenjene na kontrolne tačke krive.
Najčešće se koriste kvadratna i kubna kriva, jer krive većeg stepena troše mnogo više resursa računara daljom obradom. Kada su potrebni složeniji oblici, krive nižeg reda se spajaju i formiraju složenu Bezjeovu krivu, tzv. Bezjeov splajn. U jezicima vektorske grafike se složena Bezijerova kriva zove još i putanja. Kako bi se garantovala glatkoća krive, kontrolna tačka u kojoj se dve krive spajaju mora da bude na liniji između dve kontrolne tačke sa obe strane.
Standard za vektorsku grafiku naziva se SVG, a postoje programi za vektorsku grafiku: Adobe Illustrator, CorelDraw, Inkscape...
Najčešće korišćeni metod za rasterizovanje Bezjeove krive je rekurzivna podela, u kojoj se proveravaju kontrolne tačke, kako bi se otkrilo da li je kriva dovoljno glatka. U slučaju da nije, kriva se deli parametarski na dva segmenta i ista rekurzivna procedura se sprovodi nad oba dela.
Animacije
[uredi | uredi izvor]U Adobe Flash i Synfig animacijama, Bezijerove krive se koriste za konstruisanje. Recimo da korisnik želi da se njegov objekat kreće određenom putanjom. On će tu putanju zadati kao Bezjeovu krivu, a aplikacija će na osnovu unosa napraviti potrebne okvire po kojima će se objekat kretati.
Konstrukcija krivih
[uredi | uredi izvor]Linearne krive
[uredi | uredi izvor]Parametar kod linearne Bezjeove krive može biti protumačen kao „koliko je daleko od krećući se ka . Na primer, kada je , je na pola puta od ka . Promenom parametra od 0 do 1, ispisuje se prava linija.
Kvadratne krive
[uredi | uredi izvor]Dok parametar varira od 0 do 1, moguće je uočiti sledeće:
- Tačku koja varira od do i time ispisuje linearnu Bezjeovu krivu
- Tačku koja varira od do i time ispisuje linearnu Bezjeovu krivu
- Tačku koja je linearna interpolacija tačaka i i time ispisuje kvadratnu Bezjeovu krivu
Krive višeg stepena
[uredi | uredi izvor]- Kubna Bezjeova kriva ima tri tačke koje ispisuju linearne krive i tačke koje opisuju kvadratne Bezjeove krive.
- Kriva stepena četiri ima četiri tačke za opis linearnih krivih, tri tačke koje opisuju kvadratne krive i dve tačke koje opisuju kubne Bezjeove krive.
- Analogijom konstruišemo i krivu petog stepena, kao i krive većih stepena
Reference
[uredi | uredi izvor]- ^ Farin, Gerald E.; Hoschek, Josef; Myung-Soo Kim (2002). Handbook of Computer Aided Geometric Design. Elsevier. стр. 4—6. ISBN 978-0-444-51104-1.
- ^ http://www.idav.ucdavis.edu/education/CAGDNotes/Matrix-Cubic-Bezier-Curve/Matrix-Cubic-Bezier-Curve.html Архивирано на сајту Wayback Machine (4. март 2016) Matrična reprezentacija, 29.01.2016.
- ^ http://codegolf.stackexchange.com/questions/21178/animated-drawing-of-a-b%C3%A9zier-curve Programerski kodovi, 29.01.2016.
Literatura
[uredi | uredi izvor]- Farin, Gerald E.; Hoschek, Josef; Myung-Soo Kim (2002). Handbook of Computer Aided Geometric Design. Elsevier. стр. 4—6. ISBN 978-0-444-51104-1.
- Šukilović, T. Vukmirović, S. (2015). Geometrija za informatičare. Beograd, Matematički fakultet u Beogradu. стр. 139—149. ISBN 978-86-7589-106-2.
Video snimci
[uredi | uredi izvor]- Animacija Bezijerovih krivih
- Animacija Bezijerovih krivih
- Crtanje hobotnice pomoću Bezijerovih krivih
- GeoGebra
- Bezijerove krive u PowerPoint-u
- OpenGL
- Bezijerove krive u računarskoj grafici
- JavaScript tutorial za Bezijerove krive
- Splajnovi
- Mali čas Bezijerovih krivih
- Racionalne Bezijerove krive
- Inkscape tutorial
- Inkscape tutorial: nacrtati macu pomoću Bezijerovih krivih