Error in tile extraction?
In gsg_get_tm_coords, lines 70/71, indices to extract tiles from the tomogram are computed based on the tile-start-indices (bs
; coming from a non-overlapping tiling of the tomogram or the masked region thereof) like this:
% Extraction start
es = bs - repmat(floor(o.padsize/2),[n_tiles,3]);
ee = es + repmat(tilesize-1,[n_tiles,1]);
o.padsize
and tilesize
(which is o.tilesize
) are computed beforehand from the size of the template in gsg_determine_tile_size
like this:
% Decide padding size on scoring algorithm
switch s.scoring_fcn
case 'flcf'
o.padsize = o.tmpl_size;
case 'scf'
o.padsize = o.tmpl_size.*2;
otherwise
error('ACHTUNG!!! Invalid scoring function!!!');
end
% Set tilesize
o.tilesize = o.patchsize + o.tmpl_size;
If I am not mistaking, this means that the extraction is only symmetric around the tile for the 'flcf' scoring function. Is this intended or a bug in gsg_determine_tile_size
? Not being symmetric is not a problem per-se but can lead to odd effects of tiles on the boundary basically being filled with mean-values and thus computed for nothing and wasting resources + some data on the other side of the tomogram might be skipped from the correlation.