diff options
Diffstat (limited to 'SD-VBS/benchmarks/tracking/src')
11 files changed, 0 insertions, 600 deletions
diff --git a/SD-VBS/benchmarks/tracking/src/matlab/Filtering.m b/SD-VBS/benchmarks/tracking/src/matlab/Filtering.m deleted file mode 100755 index d42ebb7..0000000 --- a/SD-VBS/benchmarks/tracking/src/matlab/Filtering.m +++ /dev/null | |||
| @@ -1,30 +0,0 @@ | |||
| 1 | function imageOut = Filtering(imageIn, rows, cols, kernel, kernelSize) | ||
| 2 | |||
| 3 | imageOut = zeros(rows, cols);%initalize output image to all zeros | ||
| 4 | imageIn = double(imageIn);%convert to double to allow image arithmetic | ||
| 5 | |||
| 6 | intialCol = ((kernelSize+1)/2); | ||
| 7 | endCol = round(cols - ((kernelSize+1)/2)); | ||
| 8 | halfKernel = (kernelSize-1)/2; | ||
| 9 | |||
| 10 | initialRow = ((kernelSize+1)/2); | ||
| 11 | endRow = (rows - ((kernelSize+1)/2)); | ||
| 12 | |||
| 13 | %% Start 1-D filtering row-wise first. | ||
| 14 | |||
| 15 | for i=initialRow:endRow | ||
| 16 | for j=initialCol:endCol | ||
| 17 | imageOut(i,j) = sum(imageIn(i,j-halfKernel:j+halfKernel).*kernel)/sum(kernel);%actual filtering step | ||
| 18 | end | ||
| 19 | end | ||
| 20 | |||
| 21 | %% Start 1-D filtering col-wise first. | ||
| 22 | |||
| 23 | % kernelT = kernel'; | ||
| 24 | % for i=initialRow:endRow | ||
| 25 | % for j=initialCol:endCol | ||
| 26 | % imageOut(i,j) = sum(imageOut(i-halfKernel:i+halfKernel,j).*kernelT)/sum(kernel);%kernel to be transposed for performing multiplcation | ||
| 27 | % end | ||
| 28 | % end | ||
| 29 | % | ||
| 30 | % %imshow(uint8(imageOut)); | ||
diff --git a/SD-VBS/benchmarks/tracking/src/matlab/calcAreaSum.m b/SD-VBS/benchmarks/tracking/src/matlab/calcAreaSum.m deleted file mode 100755 index f8dd3da..0000000 --- a/SD-VBS/benchmarks/tracking/src/matlab/calcAreaSum.m +++ /dev/null | |||
| @@ -1,48 +0,0 @@ | |||
| 1 | %calcAreaSum: | ||
| 2 | % sizeX = cols | ||
| 3 | % sizeY = rows | ||
| 4 | function ret = calcAreaSum(src, sizeX, sizeY, winSize, dataDir) | ||
| 5 | |||
| 6 | src = double(src); | ||
| 7 | nave = winSize; | ||
| 8 | nave_half = floor((nave+1)/2); | ||
| 9 | a1=zeros(1,sizeX+nave); | ||
| 10 | |||
| 11 | for i=1:sizeY | ||
| 12 | %pull out one row | ||
| 13 | for j=1:sizeX | ||
| 14 | a1(j+nave_half)=src(i,j); | ||
| 15 | end | ||
| 16 | |||
| 17 | a1sum=0; | ||
| 18 | %sum up values within a window | ||
| 19 | for k=1:nave | ||
| 20 | a1sum = a1sum+a1(k); | ||
| 21 | end | ||
| 22 | |||
| 23 | for j=1:sizeX | ||
| 24 | ret(i,j) = a1sum; | ||
| 25 | a1sum = a1sum + a1(j+nave) - a1(j); | ||
| 26 | end | ||
| 27 | end | ||
| 28 | |||
| 29 | a1=zeros(1,sizeY+nave); | ||
| 30 | for i=1:sizeX | ||
| 31 | |||
| 32 | %pull out one col | ||
| 33 | for j=1:sizeY | ||
| 34 | a1(j+nave_half)=ret(j,i); | ||
| 35 | end | ||
| 36 | |||
| 37 | a1sum=0; | ||
| 38 | %sum up values within a window | ||
| 39 | for k=1:nave | ||
| 40 | a1sum = a1sum+a1(k); | ||
| 41 | end | ||
| 42 | |||
| 43 | for j=1:sizeY | ||
| 44 | ret(j,i) = a1sum; | ||
| 45 | a1sum = a1sum + a1(j+nave) - a1(j); | ||
| 46 | end | ||
| 47 | end | ||
| 48 | |||
diff --git a/SD-VBS/benchmarks/tracking/src/matlab/calcGoodFeature.m b/SD-VBS/benchmarks/tracking/src/matlab/calcGoodFeature.m deleted file mode 100755 index b83faa6..0000000 --- a/SD-VBS/benchmarks/tracking/src/matlab/calcGoodFeature.m +++ /dev/null | |||
| @@ -1,28 +0,0 @@ | |||
| 1 | %calcGoodFeature: | ||
| 2 | |||
| 3 | function [lambda, tr, det,c_xx, c_xy, c_yy] = calcGoodFeature(dX, dY, sizeX, sizeY, winSize, dataDir) | ||
| 4 | |||
| 5 | for i=1:sizeY | ||
| 6 | for j=1:sizeX | ||
| 7 | xx(i,j)=dX(i,j)*dX(i,j); | ||
| 8 | xy(i,j)=dX(i,j)*dY(i,j); | ||
| 9 | yy(i,j)=dY(i,j)*dY(i,j); | ||
| 10 | end | ||
| 11 | end | ||
| 12 | |||
| 13 | c_xx=calcAreaSum(xx, sizeX, sizeY, winSize,dataDir); | ||
| 14 | c_xy=calcAreaSum(xy, sizeX, sizeY, winSize,dataDir); | ||
| 15 | c_yy=calcAreaSum(yy, sizeX, sizeY, winSize,dataDir); | ||
| 16 | |||
| 17 | for i=1:sizeY | ||
| 18 | for j=1:sizeX | ||
| 19 | tr(i,j)= c_xx(i,j)+c_yy(i,j); | ||
| 20 | det(i,j)= c_xx(i,j)*c_yy(i,j)-c_xy(i,j)*c_xy(i,j); | ||
| 21 | % if( tr(i,j) == 0 ) | ||
| 22 | % lambda(i,j) = 0; | ||
| 23 | % else | ||
| 24 | lambda(i,j)=det(i,j)/(tr(i,j) + 0.00001); | ||
| 25 | % end | ||
| 26 | end | ||
| 27 | end | ||
| 28 | |||
diff --git a/SD-VBS/benchmarks/tracking/src/matlab/calcPyrLKTrack.m b/SD-VBS/benchmarks/tracking/src/matlab/calcPyrLKTrack.m deleted file mode 100755 index c965c2b..0000000 --- a/SD-VBS/benchmarks/tracking/src/matlab/calcPyrLKTrack.m +++ /dev/null | |||
| @@ -1,108 +0,0 @@ | |||
| 1 | function [newFPnt, valid] = calcPyrLKTrack(iP, iDxP, iDyP, jP, fPnt, nFeatures, winSize, accuracy_th, max_iter) | ||
| 2 | |||
| 3 | cellDims = size(iP); | ||
| 4 | GOOD_FEATURE_LAMBDA_TH = accuracy_th; | ||
| 5 | |||
| 6 | for i=1:(cellDims(1)) | ||
| 7 | curImgDims = size(iP{i}); | ||
| 8 | imgDims(i,1)= curImgDims(1); | ||
| 9 | imgDims(i,2)= curImgDims(2); | ||
| 10 | |||
| 11 | end | ||
| 12 | |||
| 13 | pLevel = cellDims(1); | ||
| 14 | |||
| 15 | rate=[1, 0.5, 0.25, 0.125, 0.0625, 0.03125]; | ||
| 16 | winSizeSq=4*winSize*winSize; | ||
| 17 | iPatch=cell(1, winSizeSq); | ||
| 18 | jPatch=cell(1, winSizeSq); | ||
| 19 | iDxPatch=cell(1,winSizeSq); | ||
| 20 | iDyPatch=cell(1,winSizeSq); | ||
| 21 | |||
| 22 | valid(1:nFeatures) = 1; | ||
| 23 | newFPnt = zeros(2,nFeatures); | ||
| 24 | |||
| 25 | for i=1:nFeatures | ||
| 26 | |||
| 27 | dX=0; | ||
| 28 | dY=0; | ||
| 29 | |||
| 30 | %% x is rows here and y is cols | ||
| 31 | |||
| 32 | x=fPnt(1,i)*rate(pLevel+1); %half size of real level | ||
| 33 | y=fPnt(2,i)*rate(pLevel+1); | ||
| 34 | |||
| 35 | for level=pLevel:-1:1 | ||
| 36 | |||
| 37 | x = x+x; | ||
| 38 | y = y+y; | ||
| 39 | dX = dX+dX; | ||
| 40 | dY = dY+dY; | ||
| 41 | imgSize(1)=imgDims(level,1); %y,x | ||
| 42 | imgSize(2)=imgDims(level,2); %y,x | ||
| 43 | |||
| 44 | c_xx = 0; | ||
| 45 | c_xy = 0; | ||
| 46 | c_yy = 0; | ||
| 47 | |||
| 48 | %when feature goes out to the boundary. | ||
| 49 | |||
| 50 | if ((x-winSize)<1 || (y-winSize)<1 || (y+winSize+1)>imgSize(1) || (x+winSize+1)>imgSize(2)) | ||
| 51 | %winSize+1due to interpolation | ||
| 52 | %error or skip the level?? | ||
| 53 | valid(i) = 0; | ||
| 54 | break; | ||
| 55 | end | ||
| 56 | |||
| 57 | |||
| 58 | iPatch = getInterpolatePatch(iP{level}, imgSize(1), imgSize(2), x, y, winSize); | ||
| 59 | iDxPatch = getInterpolatePatch(iDxP{level}, imgSize(1), imgSize(2), x, y, winSize); | ||
| 60 | iDyPatch = getInterpolatePatch(iDyP{level}, imgSize(1), imgSize(2), x, y, winSize); | ||
| 61 | |||
