//--------------------------------------------------------------------------------------------------------
// function to compute dx, dy and dt for motion estimation
//--------------------------------------------------------------------------------------------------------
void OpticalFlow::getDxs(DImage &imdx, DImage &imdy, DImage &imdt, const DImage &im1, const DImage &im2)
{
//double gfilter[5]={0.01,0.09,0.8,0.09,0.01};
double gfilter[5]={0.02,0.11,0.74,0.11,0.02};
//double gfilter[5]={0,0,1,0,0};
if(1)
{
//DImage foo,Im;
//Im.Add(im1,im2);
//Im.Multiplywith(0.5);
////foo.imfilter_hv(Im,gfilter,2,gfilter,2);
//Im.dx(imdx,true);
//Im.dy(imdy,true);
//imdt.Subtract(im2,im1);
DImage Im1,Im2,Im;
im1.imfilter_hv(Im1,gfilter,2,gfilter,2);
im2.imfilter_hv(Im2,gfilter,2,gfilter,2);
Im.copyData(Im1);
Im.Multiplywith(0.4);
Im.Add(Im2,0.6);
//Im.Multiplywith(0.5);
//Im1.copyData(im1);
//Im2.copyData(im2);
Im.dx(imdx,true);
Im.dy(imdy,true);
imdt.Subtract(Im2,Im1);
}
else
{
// Im1 and Im2 are the smoothed version of im1 and im2
DImage Im1,Im2;
im1.imfilter_hv(Im1,gfilter,2,gfilter,2);
im2.imfilter_hv(Im2,gfilter,2,gfilter,2);
//Im1.copyData(im1);
//Im2.copyData(im2);
Im2.dx(imdx,true);
Im2.dy(imdy,true);
imdt.Subtract(Im2,Im1);
}
imdx.setDerivative();
imdy.setDerivative();
imdt.setDerivative();
}
//--------------------------------------------------------------------------------------------------------
// function to compute dx, dy and dt for motion estimation
//--------------------------------------------------------------------------------------------------------
void OpticalFlow::getDxs(DImage &imdx, DImage &imdy, DImage &imdt, const DImage &im1, const DImage &im2)
{
// Im1 and Im2 are the smoothed version of im1 and im2
DImage Im1,Im2;
double gfilter[5]= {0.05,0.2,0.5,0.2,0.05};
im1.imfilter_hv(Im1,gfilter,2,gfilter,2);
im2.imfilter_hv(Im2,gfilter,2,gfilter,2);
//Im1.copyData(im1);
//Im2.copyData(im2);
Im2.dx(imdx,true);
Im2.dy(imdy,true);
imdt.Subtract(Im2,Im1);
imdx.setDerivative();
imdy.setDerivative();
imdt.setDerivative();
}