假设我有一个元素为 0 和 1 的数组s。据我了解，Hoshen-Kopelman 算法s将是（在 3 维中d=3）：

label=labelmatrix(bwconncomp(s,2*d));

如果我想做同样的事情，但现在s会有周期性的边界条件（特别是 3 维，如果重要的话），我该怎么做？

编辑：在搜索了更多之后，我偶然发现了这个：

https://uk.mathworks.com/matlabcentral/fileexchange/66079-cc2periodic?s_tid=mwa_osa_a

编辑：如果这是题外话，那么我想尽可能将其移至堆栈溢出。