The 20*log (1+32p/1p) part is not correct. This is somewhat mixing with the dB scale. However the numerical result for 20*log .. is not far off the corret factor without the 20*log part. So by lucky coincidence the mistake does not change much.
The calculation assumes that the ouput is sampled not much faster than some 50 kHz, so it does not include much higher frequency noise. Chances are this is OK for the application.
There would be a small contribution from low frequency 1/f noise of the OP. So kind of add the 0.1-10 Hz RMS noise to the final result. This would be negligible in this case with high TIA gain and high BW, but could be an issue in other cases.