We present a method, based on a single scattering model, to calculate the attenuation coefficient of each pixel in optical coherence tomography (OCT) depth profiles. Numerical simulations were used to determine the model's response to different depths and attenuation coefficients. Experiments were performed on uniform and layered phantoms with varying attenuation coefficients. They were measured by a 1300 nm OCT system and their attenuation coefficients were evaluated by our proposed method and by fitting the OCT slope as the gold standard. Both methods showed largely consistent results for the uniform phantoms. On the layered phantom, only our proposed method accurately estimated the attenuation coefficients. For all phantoms, the proposed method largely reduced the variability of the estimated attenuation coefficients. The method was illustrated on an in-vivo retinal OCT scan, effectively removing common imaging artifacts such as shadowing. By providing localized, per-pixel attenuation coefficients, this method enables tissue characterization based on attenuation coefficient estimates from OCT data.