Despite three decades of effort, predicting accurately extraordinary transmission through subwavelength hole arrays has proven challenging. This lack of quantitative design and modelling able to take into account the inherent complexity of high frequency instrumentation has prevented the development of practical high-performance components based on this phenomenon. This manuscript resorts to the Method of Moments to provide, not only such missing quantitative prediction, but also a theoretical framework to understand and shed more light on the far-field and near-field physics of the extraordinary terahertz (THz) transmission through subwavelength hole arrays under different illumination and detection conditions. An excellent agreement between numerical and experimental results with various illumination and detection setups is obtained, demonstrating the suitability of this computationally efficient modelling tool to predict the response of extraordinary transmission structures in practical situations.