import utility.utility_bokeh as ut_bokeh import utility.utility_detection as ut_detection from webService.backend.coreSketchDyBinService import sketchDyBinService from csv import DictWriter windowSize = 20 kmax = 9 initialBin = 3 durations = [500, 60, 100, 200] # height = 3 heights = [1] if __name__ == '__main__': for duration in durations: for height in heights: listFile = ut_light.getListLight(height=height, duration=duration, isFullPath=True) listAns = ut_light.getListLight(height=height, duration=duration, isFullPath=True, folderType="answer") output3D_htmlFile = "{}3D_html\\".format( ut_light.getFullPath(height=height, duration=duration)) output3D_JPGFile = "{}3D_JPG\\".format( ut_light.getFullPath(height=height, duration=duration)) outputdyFile = "{}dy_html\\".format( ut_light.getFullPath(height=height, duration=duration)) for mainFile, ansFile in zip(listFile, listAns): file_name = os.path.basename(mainFile) mainRaw = [float(item) for item in ut.txt_to_list(mainFile)]
import pywt import numpy as np import matplotlib.pyplot as plt import utility.utility_light_curve as ut_lc heights = [1,3] durations = [60,100,200,500] i = 38 if __name__ == '__main__': for height in heights: for duration in durations: listFiles = ut_lc.getListLight(height=height, duration=duration) for lc_file in listFiles[i:i+1]: data = ut_lc.getDataFromFile(fileName=lc_file, height=height, duration=duration) fig, axs = plt.subplots(2, 1) axs[0].plot(data["timestamp"], data["instances"]) coef, freqs=pywt.cwt(data["instances"],np.arange(1,129),'gaus1') # plt.subplot(2, 1, 2) axs[1].matshow(coef) # doctest: +SKIP axs[1].autoscale() plt.show() # doctest: +SKIP plt.clf()
print("a") indexFromMaxBin = [*range(start_point, end_point, 1)] return indexFromMaxBin if __name__ == '__main__': # data = pd.Series([1, 8, 9, 10, 9]) # print(grubbs.max_test_outliers(data, alpha=0.05)) # result = grubbs.max_test_outliers(data, alpha=0.05) for window_size in window_sizes: for I in IS: for L in LS: print("start L {}, I {} ,window_size {}".format( L, I, window_size)) files_list = ut_lc.getListLight(height=L, duration=I) rows = [] for index, file_name in enumerate(files_list): print("file : {}".format(file_name)) rawInstances, transList = ut_service.getDataLC_test( file_name, L=L, I=I) timestamp = [*range(0, len(rawInstances), 1)] FdataMDF_dy = ut_gen.genListDyBin(instances=rawInstances, timestamp=timestamp, windowSize=window_size, initialBin=3) # FdataMDF_Fix = ut_gen.genFixBin(binSize=binSize, instances=rawInstances,isExtract= False) gb_core = gb_test() for instance in FdataMDF_dy["window"]: gb_core.add_element(instance.get_representation()) for alpha in alphas: