Each reflect each boundary value outwards This mimics what the skimage.filters.gaussian_filter(... mode="nearest") does.
(x, axis, depth)
| 247 | |
| 248 | |
| 249 | def nearest(x, axis, depth): |
| 250 | """Each reflect each boundary value outwards |
| 251 | |
| 252 | This mimics what the skimage.filters.gaussian_filter(... mode="nearest") |
| 253 | does. |
| 254 | """ |
| 255 | left = ( |
| 256 | (slice(None, None, None),) * axis |
| 257 | + (slice(0, 1),) |
| 258 | + (slice(None, None, None),) * (x.ndim - axis - 1) |
| 259 | ) |
| 260 | right = ( |
| 261 | (slice(None, None, None),) * axis |
| 262 | + (slice(-1, -2, -1),) |
| 263 | + (slice(None, None, None),) * (x.ndim - axis - 1) |
| 264 | ) |
| 265 | |
| 266 | l = repeat(x[left], depth, axis=axis) |
| 267 | r = repeat(x[right], depth, axis=axis) |
| 268 | |
| 269 | l, r = _remove_overlap_boundaries(l, r, axis, depth) |
| 270 | |
| 271 | return concatenate([l, x, r], axis=axis) |
| 272 | |
| 273 | |
| 274 | def constant(x, axis, depth, value): |