Mercurial > hg > octave-image

changeset 355:7f70d7848175

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix indexing (by Dongik Shin)
author hauberg
date Sat, 13 Dec 2008 12:55:23 +0000
parents 594055014a19
children 83b86e861886
files inst/immaximas.m
diffstat 1 files changed, 4 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/inst/immaximas.m
+++ b/inst/immaximas.m
@@ -80,12 +80,13 @@
   [sub{1:nd}] = ind2sub(s, ind);
   if (!isempty(ind))
     w = 1; # Width that we look out on each side of the feature point to fit a local parabola
-
+    ws = w*cumprod([1; s(:)]);
+    
     ## We fit a parabola to the points in each dimension
     for d = 1:nd
       ## Indices of points above, below, left and right of feature point
-      indminus1 = max(ind-w,1);
-      indplus1  = min(ind+w,numel(immx));
+      indminus1 = max(ind-ws(d), 1);
+      indplus1  = min(ind+ws(d), numel(immx));
 
       ## Solve quadratic
       c = im(ind);