te += 1
#estimateNoise(im,mode='basic'):
tc += 1
try:
mean, std = estimateNoise(im, mode='basic')
nm = makeNoiseMap(im.shape, mean, std)
print nm.shape
except:
print 'Test failed (%i):' % tc, sys.exc_info()[0]
te += 1
#constructModel(bvals,coeffs,size):
tc += 1
try:
shapeDict = fileio.readLageurreCoeffs('../data/testHermite.pkl')
rx = np.array(range(0, shapeDict['size'][0]),
dtype=float) - shapeDict['xc'][0]
ry = np.array(range(0, shapeDict['size'][1]),
dtype=float) - shapeDict['xc'][1]
xx, yy = shapelet.xy2Grid(rx, ry)
bvals = decomp.genBasisMatrix(shapeDict['beta'], shapeDict['norder'],
shapeDict['phi'], xx, yy)
mdl = constructModel(bvals, shapeDict['coeffs'], shapeDict['size'])
print mdl.shape
except:
print 'Test failed (%i):' % tc, sys.exc_info()[0]
te += 1
#polarCoeffImg(coeffs,nmax):
tc += 1
print 'Error: Unknown mode'
return np.nan
if __name__ == "__main__":
print '============================================'
print 'Testing measure module:'
print '============================================'
import fileio
import sys
tc=0
te=0
#load precomputed shapelet coeffs (polar and cartesian)
hermDict=fileio.readHermiteCoeffs('../data/testHermite.pkl')
laDict=fileio.readLageurreCoeffs('../data/testLageurre.pkl')
#polarIDs(nmax):
tc+=1
try:
print polarIDs(5)
print cartIDs(5)
except:
print 'Test failed (%i):'%tc, sys.exc_info()[0]
te+=1
#compute flux (cartesian)
#compute flux (polar)
tc+=1
try:
print flux(hermDict['coeffs'],hermDict['beta'],hermDict['norder'],mode=hermDict['mode'])
return np.nan
if __name__ == "__main__":
print '============================================'
print 'Testing measure module:'
print '============================================'
import fileio
import sys
tc = 0
te = 0
#load precomputed shapelet coeffs (polar and cartesian)
hermDict = fileio.readHermiteCoeffs('../data/testHermite.pkl')
laDict = fileio.readLageurreCoeffs('../data/testLageurre.pkl')
#polarIDs(nmax):
tc += 1
try:
print polarIDs(5)
print cartIDs(5)
except:
print 'Test failed (%i):' % tc, sys.exc_info()[0]
te += 1
#compute flux (cartesian)
#compute flux (polar)
tc += 1
try:
print flux(hermDict['coeffs'],
