simplematching.tcl

# $Id$ #

package require image
package require gfm

namespace eval gfm {
  proc computeSimpleMatch {
    {fileName1}
    {fileName2}
    {reduction 2}
    {scale 1}
    {desiredNb 1000}
    {threshold 3.0}
    {doDisplay true}
  } {
    set ref [image::JfrImage_loadImage $fileName1]
    set match [image::JfrImage_loadImage $fileName2]

    if {$reduction==1} {
      set ref_r $ref
      set match_r $match
    } else {
    # resize the ref image
      set ref_r [image::new_JfrImage \
        [expr round([$ref getWidth]/$reduction)]  \
        [expr round([$ref getHeight]/$reduction)] \
        [$ref getDepth] [$ref getnChannels]]
      $ref resize $ref_r

    # resize the match image  
      set match_r [image::new_JfrImage \
        [expr round([$match getWidth]/$reduction)]  \
        [expr round([$match getHeight]/$reduction)] \
        [$match getDepth] [$match getnChannels]]
      $match resize $match_r
    }

    set harrisdetector [gfm::new_HarrisDetector $scale 1.0 2.0 1 $threshold "true" $desiredNb]
#    set znccfactory [gfm::new_ZnccDescriptorFactory 4]
    set simplematcher [gfm::new_SimpleMatcher 0.9 20 20 ]
    kernel::tic
    set engine [gfm::new_Engine $harrisdetector $simplematcher ]
    set result [gfm::Engine_computeMatch $engine $ref_r $match_r ]
    putsDbg "Time to match: [kernel::toc] ms"
    if {$doDisplay } {
      package require display
      display::show $ref_r
      displayHarrisPoints $ref_r [Result_points_ref_get $result ]
      display::show $match_r
      displayHarrisPoints $match_r [Result_points_match_get $result ]
      displayMatch $ref_r $match_r [Result_matches_get $result ]
}
  }

}

package provide gfm 0.4