The disclosure relates to a system and method for image reconstruction. The method may include the steps of: obtaining raw data corresponding to radiation rays within a volume, determining a radiation ray passing a plurality of voxels, grouping the voxels into a plurality of subsets such that at least some subset of voxels are sequentially loaded into a memory, and performing a calculation relating to the sequentially loaded voxels. The radiation ray may be determined based on the raw data. The calculation may be performed by a plurality of processing threads in a parallel hardware architecture. A processing thread may correspond to a subset of voxels.