Mercurial Hosting > traffic-intelligence
view python/utils.py @ 16:9d6831cfe675
added tools for speed
author | Nicolas Saunier <nico@confins.net> |
---|---|
date | Wed, 18 Nov 2009 18:36:52 -0500 |
parents | e7bbe8465591 |
children | 5a21d2cfee44 |
line wrap: on
line source
#! /usr/bin/env python ''' Generic utilities.''' #from numpy import * #from pylab import * __metaclass__ = type commentChar = '#'; def openCheck(filename, option = 'r', quit = False): '''Open file filename in read mode by default and checks it is open >>> f = openCheck('non_existant_file.txt') File non_existant_file.txt could not be opened. ''' try: return open(filename, option) except IOError: print 'File %s could not be opened.' % filename if quit: from sys import exit exit() return None def readline(f): '''Modified readline function to skip comments.''' s = f.readline() while (len(s) > 0) and s.startswith(commentChar): s = f.readline() return s.strip() def removeExtension(filename, delimiter = '.'): '''Returns the filename minus the extension (all characters after last .) >>> removeExtension('test-adfasdf.asdfa.txt') 'test-adfasdf.asdfa' >>> removeExtension('test-adfasdf') 'test-adfasdf' ''' i = filename.rfind(delimiter) if i>0: return filename[:i] else: return filename def listfiles(dirname, extension, remove = False): '''Returns the list of files with the extension in the directory dirname If remove is True, the filenames are stripped from the extension''' from os import listdir tmp = [f for f in listdir(dirname) if f.endswith(extension)] tmp.sort() if remove: return [removeExtension(f, extension) for f in tmp] else: return tmp def removeFile(filename): '''Deletes the file while avoiding raising an error if the file does not exist''' if (os.path.exists(filename)): os.remove(filename) def invertHomography(homography): 'Returns an inverted homography' invH = inv(homography) invH /= invH[2,2] return invH def project(homography, p): '''Returns the coordinates of the projection of the point p through homography''' from numpy.core._dotblas import dot from numpy.core.multiarray import array from numpy.lib.function_base import insert if (homography!=None) and (len(homography)>0): pAugmented = insert(array(p), [2],[1], axis=0) projected = dot(homography, pAugmented) projected[0] /= projected[2] projected[1] /= projected[2] else: projected = p return projected[:2] def printPoint(x,y): return '(%f,%f)'%(x,y) if __name__ == "__main__": import doctest import unittest #suite = doctest.DocFileSuite('tests/ubc_utils.txt') suite = doctest.DocTestSuite() unittest.TextTestRunner().run(suite) #doctest.testmod() #doctest.testfile("example.txt")