do_regions_file.m

phow is the code for extracting pyramid features from image

function supportRegion(par,image)

% elipse1 = sprintf('%f\t%f\t%f\n',0.08,0.00000001,0.08);
% elipse2 = sprintf('%f\t%f\t%f\n',0.008,0.00000001,0.008);
% elipse3 = sprintf('%f\t%f\t%f\n',0.002,0.00000001,0.002);
% elipse4 = sprintf('%f\t%f\t%f\n',0.0008,0.00000001,0.0008);

%cornerness scale angle point_type laplacian extremum_type mi11 mi12 mi21
%mi22 
elipse1 = sprintf('1000 10.000 0.0000 0 10 0 1.000000 0.000000 0.000000 1.000000\n');
elipse2 = sprintf('1000 15.000 0.0000 0 10 0 1.000000 0.000000 0.000000 1.000000\n');
elipse3 = sprintf('1000 20.000 0.0000 0 10 0 1.000000 0.000000 0.000000 1.000000\n');
elipse4 = sprintf('1000 25.000 0.0000 0 10 0 1.000000 0.000000 0.000000 1.000000\n');

s=sprintf('%s%s',par.directoriImatges,image);
I = imread(s);
    
fitxer = sprintf('%s%s.support',par.directoriSupport,image);
f = fopen(fitxer,'w');
fprintf(f,'1\n');
     
nPoints = (size(I,1)/par.M)*(size(I,2)/par.M)*4;

% for i=par.M:par.M:size(I,1)-par.M
%     for j=par.M:par.M:size(I,2)-par.M
%         nPoints = nPoints+4;
%     end
% end

%fprintf(f,'%d\n',nPoints);

fprintf(f,'100\n');
    
for i=par.M:par.M:size(I,1)-par.M
    for j=par.M:par.M:size(I,2)-par.M
        c = sprintf('%d\t%d\t%s',j,i,elipse1);
        fprintf(f,'%s',c);
        c = sprintf('%d\t%d\t%s',j,i,elipse2);
        fprintf(f,'%s',c);
        c = sprintf('%d\t%d\t%s',j,i,elipse3);
        fprintf(f,'%s',c);
        c = sprintf('%d\t%d\t%s',j,i,elipse4);
        fprintf(f,'%s',c);
    end
end
fclose(f);