Circular mean smoothing using BandMath

you can modify the maximum neighborhood size (see Band Math X with large neighborhood ) then your method with dotpr should work