# on whether or not averaging was done.
#
# Hazen 10/13
#
import os
import sys
import storm_analysis.sa_library.parameters as params
import storm_analysis.sa_utilities.std_analysis as std_analysis
# Setup
if (len(sys.argv) == 4):
input_file = sys.argv[1]
output_file = sys.argv[2]
parameters = params.Parameters(sys.argv[3])
else:
print(
"usage: track_average_correct.py <input_list.bin> <output_list.bin> <params>"
)
exit()
# Tracking
print("Tracking")
std_analysis.tracking(input_file, parameters)
# Averaging
print("Averaging")
did_averaging = False
if (parameters.radius > 0.0):
did_averaging = True
src_directory = os.path.dirname(__file__)
if(len(sys.argv)!=5):
print("usage: cs_analysis <dax_file> <params_file> <hres_file> <bin_file>")
exit()
movie_data = datareader.inferReader(sys.argv[1])
#
# FIXME:
#
# This should also start at the same frame as hres in the event of a restart.
#
i3_file = writeinsight3.I3Writer(sys.argv[4])
params = parameters.Parameters(sys.argv[2])
#
# Load the a matrix and setup the homotopy image analysis class.
#
a_mat_file = params.a_matrix
print("Using A matrix file:", a_mat_file)
a_mat = setup_A_matrix.loadAMatrix(a_mat_file)
image = movie_data.loadAFrame(0)
htia = homotopy_imagea_c.HomotopyIA(a_mat,
params.epsilon,
image.shape)
#