High-linear energy transfer (LET) ion irradiation of cell nuclei induces complex and severe DNA lesions, and foci of repair proteins are formed densely along the ion trajectory. To efficiently discriminate the densely distributed/overlapping foci along the ion trajectory, a focus recognition algorithm called FociPicker3D based on a local fraction thresholding technique was developed. We analyzed high-resolution 3D immunofluorescence microscopic focus images and obtained the kinetics and spatial development of γ-H2AX, 53BP1 and phospho-NBS1 foci in BJ1-hTERT cells irradiated with 55 MeV carbon ions and compared the results with the dynamics of double-strand break (DSB) distributions simulated using the PARTRAC model. Clusters consisting of several foci were observed along the ion trajectory after irradiation. The spatial dynamics of the protein foci supports that the foci clusters are not formed by neighboring foci but instead originate from the DSB cluster damage induced by high-LET radiations.