In this article results of complete modelling of electromagnetic field distribution in a near-field scanning probe microscope (NSOM) are presented. It is shown, that using finite difference in time domain method the NSOM signal can be computed for real tip-sample geometry. The results of such a simulation can be used to predict presence and distribution of topography related artefacts in NSOM images.