extractCenters.m 1.1 KB

123456789101112131415161718192021222324252627282930313233
  1. function fcenters=extractCenters(data,v0,sigma2,na)
  2. %extract response by segment centered around point v0, taking weights of surrounding pixels as gaussian with sigma2
  3. %nd is number of time samples
  4. %v0 is nr by 3 matrix of centers
  5. %na is the neighborhood size (in pixels, should work great for odd na)
  6. nr=size(v0,1);
  7. fprintf("Sigma %f\n",sigma2);
  8. [~,~,~,nd]=size(data);
  9. fcenters=zeros(nr,nd);
  10. na2=int8(0.5*(na+1));
  11. for fi=1:nr
  12. fs=0;
  13. fv=v0(fi,:);%zero beased coordinates
  14. idx0=int8(fv);
  15. idx=idx0;
  16. for i=1:na
  17. idx(1)=idx0(1)+i-na2;
  18. for j=1:na
  19. idx(2)=idx0(2)+j-na2;
  20. for k=1:na
  21. idx(3)=idx0(3)+k-na2;
  22. df=double(idx)-fv;
  23. fd=df*df';
  24. fw=exp(-0.5*fd/sigma2);
  25. fcenters(fi,:)=fcenters(fi,:)+fw*reshape(data(idx(1)+1,idx(2)+1,idx(3)+1,:),1,[]);
  26. fs=fs+fw;
  27. end
  28. end
  29. end
  30. fcenters(fi,:)=fcenters(fi,:)/fs;
  31. end
  32. end