harris.tcl

# $Id$ #

#
#/** harris point related macro. \file harris.tcl  \ingroup gfm */
#

package require kernel
package require image
package require gfm

namespace eval gfm {
  proc computeHarrisPoints {
    {fileName}
    {reduction 2}
    {scale 1}
    {desiredNb 1000}
    {threshold 3.0}
    {doDisplay true}
    {doReturnVecPoints false}
    {doReturnSaveResult false}
    } {
    
    set image [image::JfrImage_loadImage $fileName]
    if {$reduction==1} {
        set image_r $image
    } else {
      # resize the image
      set image_r [image::new_JfrImage [expr round([$image getWidth]/$reduction)]  [expr round([$image getHeight]/$reduction)] [$image getDepth] [$image getnChannels]]
      $image resize $image_r
    }
    
    # Start the chrono
    kernel::tic;      
    set harrisdetector [gfm::new_HarrisDetector $scale 1.0 2.0 1 $threshold "true" $desiredNb]
    set harrispoints [gfm::HarrisDetector_detectIn $harrisdetector $image_r]
    
    putsDbg "Time to detect the harris points: [kernel::toc] ms"
    
    if $doDisplay {
      package require display
      display::show $image_r
      gfm::displayHarrisPoints $image_r $harrispoints
      display::dumpEps $image_r "image.eps"
    }
  }
  proc computeHarrisLaplacePoints {
    {fileName}
    {reduction 2}
    {threshold 3.0}
    {doDisplay true}
    {doReturnVecPoints false}
    {doReturnSaveResult false}
  } {
    
    set image [image::JfrImage_loadImage $fileName]
    if {$reduction==1} {
        set image_r $image
    } else {
      # resize the image
      set image_r [image::new_JfrImage [expr round([$image getWidth]/$reduction)]  [expr round([$image getHeight]/$reduction)] [$image getDepth] [$image getnChannels]]
      $image resize $image_r
    }
    
# Start the chrono
    kernel::tic;      
    set harrislaplacedetector [gfm::new_HarrisLaplaceDetector 2.0 1 $threshold ]
    set harrislaplacepoints [gfm::HarrisLaplaceDetector_detectIn $harrislaplacedetector $image_r]
    
    putsDbg "Time to detect the harris-laplace points: [kernel::toc] ms"
    
    if $doDisplay {
      package require display
      display::show $image_r
      gfm::displayHarrisPoints $image_r $harrislaplacepoints
    }
  }

}

package provide gfm 0.4