Seed to Voxel Analysis Multiple Seeds FDR Correction

Hello! I'm currently working on a Seed to Voxel Analysis with 5 different seeds.
I am aware that one way to correct for multiple comparisons is to use Bonferroni correction and multiply the cluster-size p-FWE or p-FDR with the number of seeds (5 in my case) and see if it still survives the threshold value (0.05). However, I'm looking for a more liberal FDR approach and this is what I've got so far:

For all 5 seeds, I have extracted all the cluster-size p-uncorrected values that survive only the voxel threshold p<0.001 in Advanced Family-wise Error Control settings (that would be many more than the clusters that show up in Random Field Theory parametric statistics). After transforming all the cluster-size p-uncorrected into a single vector in R, I have used the "Benjaminie-Hochberg" method to transform them into p-FDR values. To me, this seems like the correct way to do it, as I get the same cluster-size p-FDR values on RFT when I try this method with p-values from one seed.

I'd very much appreciate it if you could please let me know if this is the correct way to do it; if not, I'd appreciate it if you could tell me other ways to compute the seed-corrected p-FDR values.

Thanks!

Soo Park