1 |
gforget |
1.1 |
|
2 |
|
|
%use the field produced in sample_proccessing/example_bin_average.m |
3 |
|
|
fld=obsmap; fld(find(isnan(fld)))=0; fld(find(mygrid.hFacCsurf==0))=NaN; |
4 |
|
|
|
5 |
|
|
%choose smoothing scale: here 3 X grid spacing |
6 |
|
|
distXC=3*mygrid.DXC; distYC=3*mygrid.DYC; |
7 |
|
|
|
8 |
|
|
%do the smoothing: |
9 |
|
|
obsmap_smooth=diffsmooth2D(fld,distXC,distYC); |
10 |
|
|
|
11 |
|
|
%display results: |
12 |
|
|
OBS_smooth=convert2array(obsmap_smooth); |
13 |
|
|
|
14 |
|
|
figure; set(gcf,'Units','Normalized','Position',[0.3 0.1 0.4 0.6]); |
15 |
|
|
imagescnan(OBS_smooth','nancolor',[1 1 1]*0.8); axis xy; caxis([-1 1]*0.1); colorbar; |
16 |
|
|
|