wassname/scikit-image
1
0
Fork 0
mirror of https://github.com/wassname/scikit-image.git synced 2026-07-21 12:50:27 +08:00
Code Issues Packages Projects Releases Wiki Activity

Remove duplicate subtraction

Browse Source
This commit is contained in:
Johannes Schönberger
2012-08-30 17:13:30 +02:00
parent 4cd1f8798b
commit 146d5a3f5b
1 changed files with 10 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files
skimage/_shared/interpolation.pyx
+10 -11
View File
@@ -123,28 +123,27 @@ cdef inline double bicubic_interpolation(double* image, int rows, int cols,
"""
cdef int r0 = <int>r
cdef int c0 = <int>c
cdef int r0 = <int>r - 1
cdef int c0 = <int>c - 1
if r < 0:
r0 -= 1
if c < 0:
c0 -= 1
# scale position to range [0, 1]
cdef double xr = (r - r0 + 1) / 3
cdef double xc = (c - c0 + 1) / 3
cdef double xr = (r - r0) / 3
cdef double xc = (c - c0) / 3
cdef double fc[4], fr[4]
cdef int pr, pc
for pr in range(r0 - 1, r0 + 3):
# row-wise cubic interpolation
for pr in range(r0, r0 + 4):
for pc in range(c0, c0 + 4):
fc[pc - c0] = get_pixel(image, rows, cols, pr, pc, mode, cval)
fr[pr - r0] = cubic_interpolation(xc, fc)
# do row-wise cubic interpolation
for pc in range(c0 - 1, c0 + 3):
fc[pc + 1 - c0] = get_pixel(image, rows, cols, pr, pc, mode, cval)
fr[pr + 1 - r0] = cubic_interpolation(xc, fc)
# do cubic interpolation for interpolated values of each row
# cubic interpolation for interpolated values of each row
return cubic_interpolation(xr, fr)