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:
