I have a partial answer to your question. I found this reasoning in
https://www.amazon.com/dp/0716703440/?tag=pfamazon01-20 so don't hesitate to read the page 72 if you find my explanation too messy, or simply wrong.
The geodesic deviation equation is
[tex]\frac{D^{2}\xi_{\alpha}}{d\tau^{2}} = -R_{\alpha\beta\gamma\delta}u^{\beta}\xi^{\gamma}u^{\delta}[/tex] in any coordinate system.
The idea is to use the method of least squares to determine the components of the Riemann tensor.
We define the 'error' between the observed value and the predicted value for the Kth particle by
[tex]\ddot{\xi}_{\alpha}^{\phantom{\alpha}K} + R_{\alpha\beta\gamma\delta}u^{\beta,K}\xi^{\gamma,K}u^{\delta,K} = \delta a_{\alpha}^{\phantom{\alpha}K}[/tex]
The expression you now have to minimize is the sum of the squares of the [itex]\delta a_{\alpha}^{\phantom{\alpha}K}[/itex]
[tex]S = \sum_{K}\eta^{\alpha\beta}\delta a_{\alpha}^{\phantom{\alpha}K}\delta a_{\beta}^{\phantom{\beta}K}[/tex]
In this expression, everything is regarded as known, except the components [itex]R_{\alpha\beta\gamma\delta}[/itex]. Minimize with respect to these components. It gives a system of 20 equations (the 20 independent components out of the 256).
It should be possible to determine the minimum number of particles by checking for which value of K the determinant of the 20 equations system 'stops to vanish'. However, this is not very elegant and might prove difficult to achieve (I haven't tried).