void afisare(double lungime, int nivel)
{
CPunct p(0,0);
Im3(lungime*3, 0, p);
Im3(lungime, nivel, p);
}
void Im3(double lungime, int nivel, CPunct p)
{
double x,y;
p.getxy(x,y);
CPunct p1(x-lungime/2,y+lungime/2);
CVector v(1,0);
v.deseneaza(p1,lungime);
p1=CPunct(x+lungime/2,y+lungime/2);
v=CVector(0,-1);
v.deseneaza(p1,lungime);
p1=CPunct(x+lungime/2,y-lungime/2);
v=CVector(-1,0);
v.deseneaza(p1,lungime);
p1=CPunct(x-lungime/2,y-lungime/2);
v=CVector(0,1);
v.deseneaza(p1,lungime);
if(nivel!=0)
{
Im3(lungime/3,nivel-1,CPunct(x-lungime,y));
Im3(lungime/3,nivel-1,CPunct(x-lungime,y-lungime));
Im3(lungime/3,nivel-1,CPunct(x-lungime,y+lungime));
Im3(lungime/3,nivel-1,CPunct(x+lungime,y));
Im3(lungime/3,nivel-1,CPunct(x+lungime,y-lungime));
Im3(lungime/3,nivel-1,CPunct(x+lungime,y+lungime));
Im3(lungime/3,nivel-1,CPunct(x,y-lungime));
Im3(lungime/3,nivel-1,CPunct(x,y+lungime));
}
}
void bench_im_rle(Pt2di sz,INT NbLabel,Pt2di SzBox,INT NbLissage)
{
typedef EliseRle::tIm tIm;
bool ModeV8 = NRrandom3() > 0.5;
Im2D<tIm,INT> Im1(sz.x,sz.y);
Im2D<tIm,INT> Im2(sz.x,sz.y);
Im2D<tIm,INT> Im3(sz.x,sz.y);
Im2D<tIm,INT> ImOri(sz.x,sz.y);
Im2D<tIm,INT> ImBox1(sz.x,sz.y,0);
Im2D<tIm,INT> ImBox2(sz.x,sz.y,0);
Fonc_Num f = Iconv(frandr()*20*NbLabel)%NbLabel;
for (INT k=0 ; k<NbLissage ; k++)
f = label_maj(f,NbLabel,Pt2di(5,5));
ELISE_COPY(Im1.all_pts() , f-2, Im1.out()|Im2.out()|Im3.out()|ImOri.out() );
ELISE_COPY(Im1.border(1) , 0 , Im1.out()|Im2.out()|Im3.out()|ImOri.out() );
tIm ** data1 = Im1.data();
tIm ** data2 = Im2.data();
Neighbourhood V8 = Neighbourhood::v8();
Neighbourhood V4 = Neighbourhood::v4();
U_INT1 ColBig = NbLabel * 2 +1;
EliseRle::tContainer Rles;
INT cptSmall1 =0;
INT cptSmall2 =0;
INT cptBig1 =0;
INT cptBig2 =0;
ELISE_COPY(select(Im1.all_pts(),Im1.in()<0),0,Im1.out());
BenchConc aBc(NbLabel,false);
BenchConc * aBcPtr = &aBc;
if (NRrandom3() >0.5)
aBcPtr = new BenchConc(NbLabel,true);
ELISE_COPY
(
Im3.all_pts(),
BoxedConc(Im3.in(0),SzBox,ModeV8,aBcPtr,(INT)(9+SzBox.y*2.9*NRrandom3())),
Im3.out()
);
for (INT y=0; y<sz.y ; y++)
{
for (INT x=0; x<sz.x ; x++)
{
tIm v1 = data1[y][x];
Neighbourhood V = ModeV8 ? V8 : V4;
Liste_Pts_INT2 l2(2);
if ((v1 >0) && (v1<=NbLabel))
{
Pt2di p0,p1;
ELISE_COPY
(
conc
(
Pt2di(x,y),
Im1.neigh_test_and_set (V,v1,ColBig,ColBig+1)
),
Virgule(FX,FY),
(l2|p0.VMin()|p1.VMax())
);
p1 += Pt2di(1,1);
if ((p1.x-p0.x<=SzBox.x) && (p1.y-p0.y<=SzBox.y))
{
cptSmall1++;
ELISE_COPY(l2.all_pts(),v1+NbLabel,Im1.out());
}
else
{
cptBig1++;
}
ELISE_COPY(rectangle(p0,p1),(ImBox1.in()+v1)%256,ImBox1.out());
}
tIm v2 = data2[y][x];
if ((v2 >0) && (v2<=NbLabel))
{
Box2di Box = EliseRle::ConcIfInBox
(
Pt2di(x,y),
Rles,
data2,
v2,ColBig,
ModeV8,
SzBox
);
if (! Rles.empty())
{
cptSmall2++;
EliseRle::SetIm(Rles,data2, v2+NbLabel);
}
else
{
cptBig2++;
}
ELISE_COPY
(
rectangle(Box._p0,Box._p1),
(ImBox2.in()+v2)%256,
ImBox2.out()
);
}
}
}
ELISE_COPY(select(Im1.all_pts(),ImOri.in()<0),ImOri.in(),Im1.out());
BENCH_ASSERT(cptSmall2 == cptSmall1);
BENCH_ASSERT(cptBig2 == cptBig1);
INT Dif1;
ELISE_COPY
(
Im1.all_pts(),
Abs(Im1.in()-Im2.in()),
VMax(Dif1)
);
BENCH_ASSERT(Dif1==0);
ELISE_COPY
(
Im1.all_pts(),
Abs(ImBox1.in()-ImBox2.in()),
VMax(Dif1)
);
BENCH_ASSERT(Dif1==0);
ELISE_COPY
(
Im1.all_pts(),
Abs(Im1.in()-Im3.in()),
VMax(Dif1)
);
BENCH_ASSERT(Dif1==0);
}
