I'm assuming your problem is with the radial dimension H. What I would do is calculate the median for each pixel such that its neighborhood is centered around 0.5 (the pixel). For example for a neighborhood size of [5 5]
colfilt(img_hsv(:,:,c), [5 5], 'sliding', @(x) mod(median(mod(x - x(13) + 0.5,1)) + x(13) - 0.5,1));
Its not perfect but gives decent results
