Project

General

Profile

Bug #2912 » testMatchImages.py

Anchi Cheng, 09/17/2014 09:38 AM

 
#!/usr/bin/env python
import time
import leginon.libCVwrapper
import pyami.mrc

# read the two images as numpy array
arraynew = pyami.mrc.read('./arraynew.mrc')
arrayold = pyami.mrc.read('./arrayold.mrc')

t0 = time.time()
# run libcv MatchImages
# feature minsize and maxsize are defined as in RCTAcquisition
# when the value is integer, it represents number of pixels in the feature
# when the value is a fractional float, it represents the percentage of total
# pixels in the image
minsize = 75
maxsize = 0.8
result = leginon.libCVwrapper.MatchImages(arrayold, arraynew, minsize, maxsize)
# timing
t1 = time.time()
print 'process time = %.1f seconds' % (t1 - t0)
(3-3/4)