Generally, least-squares rational function estimation-based approach can be used to estimate the state-space model {A, B, C ,D} using data that can be in the frequency domain or the time domain.
If you are satisfied with the estimation result, then you can proceed with the LQR design.