Head scatter factors for high energy photon beams from linear accelerators can be modeled using a two-source model consisting of focal and extrafocal radiation. The focal radiation can be approximated as a point source, and the distribution of the extrafocal radiation is a two-dimensional (2D) radial symmetric function. Various methods, including analytical, Monte Carlo, and empirical trial functions, have been used to determine the radial symmetric function of extrafocal radiation distribution. This article describes a method for directly determining the extrafocal radiation distribution without assuming any empirical trial function. The extrafocal radiation distribution is determined with measured head scatter factors for rectangular fields defined by the lower jaw (X) fixed at 40 cm and the upper jaw (Y) varying from 3 to 40 cm. The derivatives of the measured head scatter factors, with respect to the Y jaw position projected in the plane of extrafocal radiation, are proportional to the one-dimensional (1D) projection (also called the line spread function) of the extrafocal radiation distribution. Two methods are used to solve the radial function of extrafocal radiation from the 1D projection. The first method uses a 2D filtered backprojection algorithm, originally developed for parallel beam computed tomography reconstruction, to directly derive the radial dependence of the extrafocal radiation distribution. The method has been applied to 6 and 18 MV photon beams from a Siemens linear accelerator and has been tested by comparing measured and calculated head scatter factors for square and rectangular fields. The second method uses a Fourier transform followed by a Fourier-Bessel transform to solve the problem. The distributions of extrafocal radiation derived from these two methods are virtually identical.