diff --git a/dat/AldoJF250k.npz b/dat/AldoJF250k.npz new file mode 100644 index 0000000..17cc6e7 --- /dev/null +++ b/dat/AldoJF250k.npz @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:e4fc74cb0bbc8edffc76d3b80687c8b11bfc0c07815ebc04767db2aba23c0132 +size 27788624 diff --git a/dat/AldoJF250k_000000.dat b/dat/AldoJF250k_000000.dat new file mode 100644 index 0000000..3ac7857 --- /dev/null +++ b/dat/AldoJF250k_000000.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:ee5133916a3612ef3350314527aa4653356b845f14d07c6e2dd59a63222a0910 +size 4194432 diff --git a/dat/AldoJF250k_000001.dat b/dat/AldoJF250k_000001.dat new file mode 100644 index 0000000..0e02ec5 --- /dev/null +++ b/dat/AldoJF250k_000001.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:d1d9f70d9124c58b231212149121ca5a4535221d23c852226464d3a8cbd438b9 +size 4194432 diff --git a/dat/AldoJF250k_000002.dat b/dat/AldoJF250k_000002.dat new file mode 100644 index 0000000..164661f --- /dev/null +++ b/dat/AldoJF250k_000002.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:e890e23e7dd7ef2500f960b4bc4ab79cd33656bb4fdd5aa9835f3b5284d5a676 +size 4194432 diff --git a/dat/AldoJF250k_000003.dat b/dat/AldoJF250k_000003.dat new file mode 100644 index 0000000..fef2404 --- /dev/null +++ b/dat/AldoJF250k_000003.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:32a5b51f154a3c519b2472858d42b060ea68bf40910370f828a05eb3325fa669 +size 4194432 diff --git a/dat/AldoJF250k_000004.dat b/dat/AldoJF250k_000004.dat new file mode 100644 index 0000000..e2edbdb --- /dev/null +++ b/dat/AldoJF250k_000004.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:9b0ee1e850b7859bb8b656c62167cee2a8d1d4f9232cbe1973ecaae4626c0cf7 +size 4194432 diff --git a/dat/AldoJF250k_000005.dat b/dat/AldoJF250k_000005.dat new file mode 100644 index 0000000..d0cd9e2 --- /dev/null +++ b/dat/AldoJF250k_000005.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:2e87bde338e4bda8e283b71e88dd29995a16c90053ec5925bcd10bb31583a51d +size 4194432 diff --git a/dat/AldoJF250k_000006.dat b/dat/AldoJF250k_000006.dat new file mode 100644 index 0000000..8fe69f3 --- /dev/null +++ b/dat/AldoJF250k_000006.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:120ca3d34039efa5a9be87438cbc78757821516ee372f142e43c6ed49c35708d +size 2621520 diff --git a/dat/AldoJF500k.npz b/dat/AldoJF500k.npz index 7912439..03e475b 100644 --- a/dat/AldoJF500k.npz +++ b/dat/AldoJF500k.npz @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:a47abc96fd62fb4f9c23504389c81e7ad2e40f452767f89257e799a1a7f828bc +oid sha256:a4169ba17ed74e479cefa0f41b4b8f0672e8a88e1f3add845a92fe91761afba0 size 25166720 diff --git a/dat/AldoJF65k.npz b/dat/AldoJF65k.npz new file mode 100644 index 0000000..aa59243 --- /dev/null +++ b/dat/AldoJF65k.npz @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:86f656fd2850bfad41e11b62306a60cb47cafd74573f4a31849c23c4444c0e3d +size 14813456 diff --git a/dat/AldoJF65k_000000.dat b/dat/AldoJF65k_000000.dat new file mode 100644 index 0000000..03be452 --- /dev/null +++ b/dat/AldoJF65k_000000.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:643f5ab54dc5eda3dd6072f9c9817719895010e845e9211c884ed603d6db7b33 +size 2359584 diff --git a/dat/AldoJF65k_000001.dat b/dat/AldoJF65k_000001.dat new file mode 100644 index 0000000..62493c3 --- /dev/null +++ b/dat/AldoJF65k_000001.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:3f667ff4e708254b8738876e23449d1bb48e302ef8e40efc9ad06ea1d78fc603 +size 2359584 diff --git a/dat/AldoJF65k_000002.dat b/dat/AldoJF65k_000002.dat new file mode 100644 index 0000000..26144f4 --- /dev/null +++ b/dat/AldoJF65k_000002.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:666a8cdb4fd7cf65aead31c4e4080f3e1ee973b5c4fb2fe1f2c90cf54ff1c6d1 +size 2359584 diff --git a/dat/AldoJF65k_000003.dat b/dat/AldoJF65k_000003.dat new file mode 100644 index 0000000..4960612 --- /dev/null +++ b/dat/AldoJF65k_000003.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:359c7a7d6bc770ebedbccf2c8b05aa6bafaa8fd24ee05e006d928e103ee0953f +size 2359584 diff --git a/dat/AldoJF65k_000004.dat b/dat/AldoJF65k_000004.dat new file mode 100644 index 0000000..9a757eb --- /dev/null +++ b/dat/AldoJF65k_000004.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:36edcf89c46207831231dfe22f4c839502da4313f4385d3e3cc5255661866098 +size 2359584 diff --git a/dat/AldoJF65k_000005.dat b/dat/AldoJF65k_000005.dat new file mode 100644 index 0000000..7e1effd --- /dev/null +++ b/dat/AldoJF65k_000005.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:de5e961dbeef2e9de998782d97c2e348cdc7befd0a9c82fc862252f63eb62618 +size 2359584 diff --git a/dat/AldoJF65k_000006.dat b/dat/AldoJF65k_000006.dat new file mode 100644 index 0000000..892a097 --- /dev/null +++ b/dat/AldoJF65k_000006.dat @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:df942b9c7511be6ad3a8b21ccc8ac4f8dcf3402b43537eac72ef8d21145887a8 +size 655440 diff --git a/dat/write_dat_files.py b/dat/write_dat_files.py index 0e7f70a..b1950ca 100644 --- a/dat/write_dat_files.py +++ b/dat/write_dat_files.py @@ -5,36 +5,130 @@ import numpy as np import tifffile -img = tifffile.imread("AldoJF500k.tiff") +def write_dat_file_module(): + # img = tifffile.imread("AldoJF500k.tiff") -# Filename is name_000000.dat -# Header is frame number and bunch id 8+8 bytes -# Data is 2 bytes per pixel + # # Filename is name_000000.dat + # # Header is frame number and bunch id 8+8 bytes + # # Data is 2 bytes per pixel + # header = np.zeros(2, dtype=np.uint64) + # file_index = 0 + # frame_number = 1 + # frames = [] + # headers = [] + # for file_index in range(4): + # with open(f"AldoJF500k_{file_index:06d}.dat", "wb") as f: + # for i in range(1,8,1): + # frame = np.random.normal(img, np.sqrt(img)) + # frame = np.clip(frame, 0, 65535) #avoid negative values + # frame = frame.astype(np.uint16) + # header[0] = frame_number + # header[1] = frame_number**2 + # frame_number += 1 + # header.tofile(f) + # frame.tofile(f) + # headers.append(header.copy()) + # frames.append(frame.copy( )) -header = np.zeros(2, dtype=np.uint64) -file_index = 0 -frame_number = 1 -frames = [] -headers = [] -for file_index in range(4): - with open(f"AldoJF500k_{file_index:06d}.dat", "wb") as f: - for i in range(1,8,1): - frame = np.random.normal(img, np.sqrt(img)) - frame = np.clip(frame, 0, 65535) #avoid negative values - frame = frame.astype(np.uint16) - header[0] = frame_number - header[1] = frame_number**2 - frame_number += 1 - header.tofile(f) - frame.tofile(f) - headers.append(header.copy()) - frames.append(frame.copy( )) + # if file_index == 3 and i == 3: + # break - if file_index == 3 and i == 3: - break + # headers = np.array(headers) + # frames = np.array(frames) + # print(header.shape, frames.shape) + with np.load('AldoJF500k.npz') as f: + headers = f["headers"] + frames = f["frames"] -headers = np.array(headers) -frames = np.array(frames) -print(header.shape, frames.shape) + new_header = np.zeros(headers.shape[0], dtype = [('framenum', np.uint64), ('bunchid', np.uint64)]) + for i in range(headers.shape[0]): + new_header[i][0] = headers[i][0] + new_header[i][1] = headers[i][1] -np.savez('AldoJF500k.npz', headers=headers, frames=frames) \ No newline at end of file + np.savez('AldoJF500k.npz', headers=new_header, frames=frames) + +def write_dat_file_half_module(): + # img = tifffile.imread("AldoJF500k.tiff") + # img = img[0:256,:] + + # header = np.zeros(2, dtype=np.uint64) + # file_index = 0 + # frame_number = 1 + # frames = [] + # headers = [] + # for file_index in range(7): + # with open(f"AldoJF250k_{file_index:06d}.dat", "wb") as f: + # for i in range(1,9,1): + # frame = np.random.normal(img, np.sqrt(img)) + # frame = np.clip(frame, 0, 65535) #avoid negative values + # frame = frame.astype(np.uint16) + # header[0] = frame_number + # header[1] = frame_number**2 + # frame_number += 1 + # header.tofile(f) + # frame.tofile(f) + # headers.append(header.copy()) + # frames.append(frame.copy( )) + + # if file_index == 6 and i == 5: + # break + + # headers = np.array(headers) + # frames = np.array(frames) + # print(header.shape, frames.shape) + with np.load('AldoJF250k.npz') as f: + headers = f["headers"] + frames = f["frames"] + + new_header = np.zeros(headers.shape[0], dtype = [('framenum', np.uint64), ('bunchid', np.uint64)]) + for i in range(headers.shape[0]): + new_header[i][0] = headers[i][0] + new_header[i][1] = headers[i][1] + np.savez('AldoJF250k.npz', headers=new_header, frames=frames) + + +def write_dat_file_single_chip(): + # img = tifffile.imread("AldoJF500k.tiff") + # x = 580 + # y = 220 + # img = img[y:y+256, x:x+256] + # header = np.zeros(2, dtype=np.uint64) + # file_index = 0 + # frame_number = 1 + # frames = [] + # headers = [] + # for file_index in range(7): + # with open(f"AldoJF65k_{file_index:06d}.dat", "wb") as f: + # for i in range(1,19,1): + # frame = np.random.normal(img, np.sqrt(img)) + # frame = np.clip(frame, 0, 65535) #avoid negative values + # frame = frame.astype(np.uint16) + # header[0] = frame_number + # header[1] = frame_number**2 + # frame_number += 1 + # header.tofile(f) + # frame.tofile(f) + # headers.append(header.copy()) + # frames.append(frame.copy( )) + + # if file_index == 6 and i == 5: + # break + + # headers = np.array(headers) + # frames = np.array(frames) + # print(header.shape, frames.shape) + with np.load('AldoJF65k.npz') as f: + headers = f["headers"] + frames = f["frames"] + + new_header = np.zeros(headers.shape[0], dtype = [('framenum', np.uint64), ('bunchid', np.uint64)]) + for i in range(headers.shape[0]): + new_header[i][0] = headers[i][0] + new_header[i][1] = headers[i][1] + np.savez('AldoJF65k.npz', headers=new_header, frames=frames) + + +if __name__ == "__main__": + # write_dat_file_module() + write_dat_file_half_module() + write_dat_file_single_chip() \ No newline at end of file