Bezjeove krive su parametarske krive. Njihovu najbolju primenu pronalazimo u računarskoj grafici.

Bezijerova kriva trećeg stepena čiji poligon nije prost
Animacija Bezijerove krive rađena u programskom jeziku Ruby

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

уреди

Rekurzivna definicija

уреди

Ako Bezjeovu krivu određenu pomoću tački   označimo kao  . Onda:

  i
 

Eksplicitna definicija

уреди

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

уреди

Primenom binomne teoreme na definiciju krive dobijamo:

 

gde je:

 

Izvod krive stepena   je:

 

Specijalni slučajevi

уреди

Najveću praktičnu primenu imaju Bezjeove krive drugog i trećeg stepena.

Linearna Bezjeova kriva

уреди

Bezjeova kriva stepena 1, tj. kriva određena sa samo dve tačke   i  , jeste zapravo duž  , zadata sledećom parametrizacijom:

 


Kvadratna Bezjeova kriva

уреди
 
Bezijerova kriva drugog stepena


 
Bezjeova kriva drugog stepena je deo parabole

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

уреди
 
Bezjeova kriva trećeg stepena

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

уреди
  • 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

уреди

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

уреди
 
Pomeranje kontrolne tačke Bezjeove krive

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]

уреди

Kvadratne Bezjeove krive

уреди

 


Kubne Bezjeove krive

уреди

  


Gde je  


Odakle se analogijom dobija uopštenje:


  



Podela Bezjeove krive na dva dela korišćenjem matrične reprezentacije

уреди

Kvadratna Bezjeova kriva

уреди

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

уреди

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

уреди

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

уреди

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

уреди

Računarstvo

уреди

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 );
  }

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

уреди

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.


 
Animacija Bezijerove krive u programskom jeziku Java
 
Animacija Bezjeove krive u programskom jeziku Python

Konstrukcija krivih

уреди

Linearne krive

уреди

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.

 
Animacija linearne Bezjeove krive, t je iz [0,1]

Kvadratne krive

уреди

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


 
Konstrukcija kvadratne Bezjeove krive
 
Animacija kvadratne Bezjeove krive, t je iz [0,1]


Krive višeg stepena

уреди
  • Kubna Bezjeova kriva ima tri tačke   koje ispisuju linearne krive i tačke   koje opisuju kvadratne Bezjeove krive.


 
Konstrukcija kubne Bezjeove krive
 
Animacija kubne Bezjeove krive, t je iz [0,1]
  • 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.


 
Konstrukcija Bezjeove krive stepena 4
 
Animacija Bezjeove krive stepena 4, t je iz [0,1]
  • Analogijom konstruišemo i krivu petog stepena, kao i krive većih stepena
 
Animacija Bezijerove krive petog stepena

Reference

уреди
  1. ^ Farin, Gerald E.; Hoschek, Josef; Myung-Soo Kim (2002). Handbook of Computer Aided Geometric Design. Elsevier. стр. 4—6. ISBN 978-0-444-51104-1. 
  2. ^ 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.
  3. ^ http://codegolf.stackexchange.com/questions/21178/animated-drawing-of-a-b%C3%A9zier-curve Programerski kodovi, 29.01.2016.

Literatura

уреди

Video snimci

уреди

Vidi još

уреди

Spoljašnje veze

уреди