# 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) #