Timing Performance Evaluation of a Dual-Axis Rotational PET System According to NEMA NU 4-2008 Standards: A Simulation Study

This repository contains the dataset used to identify easyPET.3D’s optimal performance based on NEMA standards. In addition, explores the impact of different CTWs on PET image quality by comparing simulated electronics capable of a 300 ps CTW with a 40 ns CTW 

The dataset is provided in the easyPET format, which includes:

Header: Contains metadata about the acquisition, structured as a binary file.

 Version_binary,

header,

dates,

otherinfo,

acquisitionInfo,

stringdata,

systemConfigurations_info,

energyfactor_info,

peakMatrix_info

The header is  processed into a JSON format. e.g.

systemConfigurations_info = np.fromfile(binary_file, dtype='|S1', count=size_systemConfigurations[0]* 2).astype('|U1')systemConfigurations_info = systemConfigurations_info.tolist()systemConfigurations_info = ''.join(systemConfigurations_info)systemConfigurations_info = json.loads(systemConfigurations_info)

7 or 8 arrays: stores the primary data. Motor angles, Detector IDs, Photon Energies and time of arrival 

number_of_arrays = 7 # or 8 for simulationsreading_data = np.fromfile(binary_file, dtype=[('stepbot', np.int16), ('steptop', np.int16), ('idA', np.int16), ('idB', np.int16), ('EA', np.int16), ('EB', np.int16), ('time', np.float64), ('time_b', np.float64)])

step_bot = round(header[1], 4) / header[2]step_top = round(header[3], 4) / header[4]topRange = header[5]listMode = np.zeros((len(reading_data["stepbot"]), number_of_arrays))listMode[:, 0] = reading_data["EA"] # energias AlistMode[:, 1] = reading_data["EB"] # energias BlistMode[:, 2] = reading_data["idA"] # id AlistMode[:, 3] = reading_data["idB"] # id bif open_files_to_join: # This dont change the data in the way to merge the acquisitions files in just one listMode[:, 4] = (reading_data["stepbot"]) # *2 # bot angle remove factor of 2 because it is already saved that way listMode[:, 5] = (reading_data["steptop"] ) # top_angleelse: listMode[:, 4] = (reading_data["stepbot"]) * step_bot #*2 # bot angle remove factor of 2 because it is already saved that way listMode[:, 5] = (reading_data["steptop"] * step_top - topRange / 2) # top_anglelistMode[:, 6] = reading_data["time"] # timeif Version_binary == "Version 3": listMode[:, 7] = reading_data["time_b"] # time

15 December ---> Resolution test6 December ----> Sensitivity