lklmin.m

图像对齐（inage registration）工具箱

function [odp,odp2,odp3]=lklmin(matice,dr,dc)
%computes local minima in submatrices sized  drxdc
%odp=lokalminima(matice,dr,dc)

vel=size(matice);

pom=reshape(matice,dr,vel(1)*vel(2)/dr);
[pom2,pom3]=min(pom,[],1);
pom=reshape(pom2,vel(1)/dr,vel(2));	% minima z jednotlivych dilcich sloupku
odp2=reshape(pom3,vel(1)/dr,vel(2));	% souradnice

pom2=reshape(pom',dc,(vel(2)*vel(1)/dr)/dc);
pom4=reshape(odp2',dc,(vel(2)*vel(1)/dr)/dc);
[pom,pom3]=min(pom2,[],1);
idx=0:dc:(vel(2)*vel(1)/dr-1);
idy=pom3+idx;
odp2=pom4(idy);

odp=reshape(pom,vel(2)/dc,vel(1)/dr)';	% minima z jednotlivych dilcich sloupku
odp3=reshape(pom3,vel(2)/dc,vel(1)/dr)';	% souradnice
odp2=reshape(odp2,vel(2)/dc,vel(1)/dr)';	% souradnice