Abstract: A method of imaging a sample, the method comprising: (a) irradiating a sample with a pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 25 GHz to 100 THz; (b) determining a first parameter related to the amplitude of the radiation, which is either reflected from and/or transmitted by the sample, in the time domain; (c) calculating the value of the first parameter at a first time value relative to the value of the first parameter at a second time value which coincides with a physical feature of the dataset of the first parameter with respect to time; and (d) generating an image by plotting the value calculated in step (c) for different points of the sample.