The NMR spectrum of FeSe shows a dramatic broadening on cooling towards the bulk nematic phase at $T_s=90$ K, due to the formation of a quasi-static, short-range-ordered nematic domain structure. However, a quantitative understanding of the NMR broadening and its relationship to the nematic susceptibility is still lacking. Here, we show that the temperature and pressure dependence of the broadening is in quantitative agreement with the mean-field Edwards-Anderson parameter of an Ising-nematic model in the presence of random-field disorder introduced by non-magnetic impurities. Furthermore, these results reconcile the interpretation of NMR and Raman spectroscopy data in FeSe under pressure.