123456789101112 |
- function map=gauss_smoothing_2D(map, sigma)
- % filt=fspecial('gaussian',[ceil(3*sigma(1)) ceil(3*sigma(2))]);%,sigma(1));
- l1=ceil(3*sigma(1));
- l2=ceil(3*sigma(2));
- [X,Y]=meshgrid(-l1:l1,-l2:l2);
- filt=exp(-.5*((X/sigma(1)).^2+(Y/sigma(2)).^2));
- filt=filt/sum(filt(:));
- [m1,m2]=size(map);
- % map=repmat(map,3,3); % angular data can be tricky
- map=nanconv(map,filt,'same','edge','2d' );
- % map=map(m1+(1:m1),m2+(1:m2));
|