wassname/scikit-image
1
0
Fork 0
mirror of https://github.com/wassname/scikit-image.git synced 2026-07-15 11:25:53 +08:00
Code Issues Packages Projects Releases Wiki Activity

fix bilinear interpolation when position outside of image

Browse Source
This commit is contained in:
Johannes Schönberger
2012-08-20 22:48:34 +02:00
parent a6bdda8291
commit bfa5578979
1 changed files with 10 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files
skimage/feature/_texture.pyx
+10 -5
View File
@@ -69,7 +69,7 @@ def _glcm_loop(np.ndarray[dtype=np.uint8_t, ndim=2,
cdef _bilinear_interpolation(np.ndarray[double, ndim=2] image,
np.ndarray[double, ndim=2] coords,
np.ndarray[double, ndim=1] output,
double r0=0, double c0=0):
double r0=0, double c0=0, double cval=0):
cdef double r, c, dr, dc
cdef int i, minr, minc, maxr, maxc
@@ -82,10 +82,15 @@ cdef _bilinear_interpolation(np.ndarray[double, ndim=2] image,
maxc = <int>ceil(c)
dr = r - minr
dc = c - minc
top = (1 - dc) * image[minr, minc] + dc * image[minr, maxc]
bottom = (1 - dc) * image[maxr, minc] + dc * image[maxr, maxc]
output[i] = (1 - dr) * top + dr * bottom
if (
minr < 0 or maxr >= image.shape[0]
or minc < 0 or maxc >= image.shape[1]
):
output[i] = cval
else:
top = (1 - dc) * image[minr, minc] + dc * image[minr, maxc]
bottom = (1 - dc) * image[maxr, minc] + dc * image[maxr, maxc]
output[i] = (1 - dr) * top + dr * bottom
@cython.boundscheck(False)