This paper proposes a parallelization method for voxel-based 3D shape reconstruction. The voxel-based technique offers accurate 3D shape reconstruction by integrating information obtained from multiple images; however, it requires huge computation time for the reconstruction. Our main approach is parallel processing by voxel space division, which divides the voxel space by planes along with the lines passing through a virtual view point and assigns the process for each divided voxel space to a processing element (PE) individually. This makes it possible to reduce the amount of inter-PE communication. In addition to this, we introduce a method to decrease the amount of computation by controlling iteration count of iterative process for each voxel. Experimental results show that our method more efficiently reduces the computation time than simple horizontal domain division method.