I,d)*radii(d) ).

 This function has an additional parameter - flag. If flag==1, then all
 locations that are not strictly the biggest in their window are
 suppressed (ie, keeps max only if it is UNIQUE).  This can be useful for
 large flat regions -- nonMaxSupr(ones(30),3) does no suppression since
 all values are equal in each window, but nonMaxSuprList(ones(30),3) w
 flag==1 suppresses all locations. If flag==2 ties are broken randomly (in
 fact this allows for slightly faster execution).

 USA