def xmedian(ref,mwid): temp=np.isnan(ref) tmean=np.nanmean(ref) ref[temp]=tmean ref2=cv2.blur(ref,(mwid,mwid)) ref[temp]=ref2[temp] tempx=np.uint8(255*ref) return cv2.medianBlur(tempx,mwid)/255.0