sift.tcl

package require image
package require display
package require gfm

namespace eval sift {
    proc points2 { fileName {reduction 2}} {

        
   set image [image::Image_loadImage $fileName]
    if {$reduction==1} {
        set image_r $image
    } else {
      # resize the image
      set image_r [image::Image_resize $image [expr 1.0 / $reduction ] ]      
    }
    
    # Start the chrono
    kernel::tic;      
    set detector [sift::new_SiftDetector ]
    set points [sift::SiftDetector_detectIn $detector $image_r]
      display::show $image_r
      gfm::displayHarrisPoints $image_r $points

    }

    proc points {image {siftScale 0.5} {displayScale 0.5}} {
    # This macro extract SIFT points from an image

    set pointsVec [sift::createSiftPointsVector]
    sift::demoWrapper $image $siftScale $pointsVec
    putsDbg "Nb points: [sift::sizeSiftPointsVector $pointsVec]"

    set imFullSize [image::JfrImage_loadImage $image]
    set imResized [$imFullSize resize $displayScale]
    image::delete_JfrImage $imFullSize
    display::show $imResized
    image::delete_JfrImage  $imResized

    sift::displayPoints $imResized $pointsVec $siftScale $displayScale
    }

    proc match {image1 image2 {siftScale 0.5} {displayScale 0.5} } {
    # This macro extracts SIFT points from image1, add the points to a database.
    # Then it extracts SIFT points from image2 and looks for them in the database
    set pointsVec1 [sift::createSiftPointsVector]
    set pointsVec2 [sift::createSiftPointsVector]
    
    
    sift::demoWrapper $image1 $siftScale $pointsVec1
    sift::demoWrapper $image2 $siftScale $pointsVec2

    set imFullSize [image::JfrImage_loadImage $image1]
    set imResized1 [$imFullSize resize $displayScale]
    image::delete_JfrImage $imFullSize
    display::show $imResized1 "image 1"
    sift::displayPoints $imResized1 $pointsVec1 $siftScale $displayScale

    set imFullSize [image::JfrImage_loadImage $image2]
    set imResized2 [$imFullSize resize $displayScale]
    image::delete_JfrImage $imFullSize
    display::show $imResized2 "image 2"
    sift::displayPoints $imResized2 $pointsVec2 $siftScale $displayScale
    
    set db [sift::new_DataBase]
    $db processPoints $pointsVec1 1
    $db processPoints $pointsVec2 0

    displayIds $imResized1 $pointsVec1 $siftScale $displayScale
    displayIds $imResized2 $pointsVec2 $siftScale $displayScale
    
    }

    proc displayPoints {image pointsVec siftScale displayScale {color "red"}} {
    set scale [expr $displayScale/$siftScale]
    for {set i 0} {$i < [sizeSiftPointsVector $pointsVec]} {incr i} {
      set pt [sift::getSiftPointsVector $pointsVec $i]
      #       putsDbg [sift::print $pt]
      display::drawCross $image [expr [$pt cget -u]*$scale] [expr [$pt cget -v]*$scale] 3 1 $color {jafar sift_point}
    }
    }

    proc displayIds {image pointsVec siftScale displayScale {color "green"}} {
    set scale [expr $displayScale/$siftScale]
    for {set i 0} {$i < [sizeSiftPointsVector $pointsVec]} {incr i} {
      set pt [sift::getSiftPointsVector $pointsVec $i]
      if {[$pt cget -id] != $::sift::DBKey_NoId} {
        display::drawText $image [expr [$pt cget -u]*$scale-8] [expr [$pt cget -v]*$scale] "[$pt cget -id]" $color "c" {jafar sift_id}
        display::drawSquare $image [expr [$pt cget -u]*$scale] [expr [$pt cget -v]*$scale] 5 1 $color {jafar sift_id}
      }
    }
    }

}

package provide sift 0.2