An important part of the issue is what the detector can resolve.
Wiki puts the angular resolution of the human eye at 1 arcminute. The dark adapted pupil is about 1 cm, the theoretical resolution would be
https://en.wikipedia.org/wiki/Dawes'_limit about 10 arceconds. Taking the more conservative estimate, we covert 10 arc seconds to radians and get 5*10^-5 radians.
THis means that optimistically, the human eye might be able to distinguish objects of ##r \, \theta## = where r = (4 light years) and ##\theta## = (5*10^-5). The actual answer is probaby closer to an arcminute that the theoretical maximum , though.
This optimisitc estimate gives ~ 2*10^12 meters of resolution. The radius of the sun is 7*10^8 meters. This is much lower than what the eye is capable of resolving
This justifies mathematically what I assumed earlier without a detailed proof (I rather thought it was obvious) - the human eye cannot resolve the sun at this distance, and that we are more than justified in using the plane wave approximation. The paper explicitly mentions they use the plane wave approximation.
To use the plane wave approximation simply need to normalize the radiated energy in terms of meter^2, rather than the "solid angle". We have to do this because our detector cannot resolve said angle - we can't use geometric optics, we need to take the diffraction limit into account.
Thus we multiply the power emitted / unit area of the black body, by the surface area of the sun (to get the total power), multiply it by (1 meter^2) [if we are using MKS units] area and divide by the area of a sphere of radius of 4 light years to get the power / meter^2
This gives us the intensity in power / meter^2.
(as per
http://en.wikipedia.org/wiki/Intensity_(physics))
AFTER we have done the above, we can , if we like, convert this to a classical EM wave of amplitude E in a vacuum (see the above wiki)
I = \frac{c \epsilon_o}{2} \left| E^2 \right|
The approach that Einstein used to transform this, and probably the least mathematically demanding one is to look at how E transforms. Because E&M radition is transverse, we need to look only at how the transverse component of E transforms.
The answer can be looked up in a classic E&M textbook or on Wiki (but probably not randomly guessed at):
## {{E}_{\bot}}'= \gamma \left( \mathbf {E}_{\bot} + \mathbf{ v} \times \mathbf {B} \right)##
We can see right away that E DOES depend on our velocity. It's mistake to assume that E doesn't transform - E changes with our velocity in a manner described by the above equation.
The magnetic field B = (1/c)E (I believe, I took the trouble of working this out, but I don't have a reference) of the plane wave will contribute to our electric field in the new frame. if we work out the details one should get the standard result that E transforms as ##\gamma \left( \beta + 1\right)## which is the relativistic doppler factor.
This right away gives us the result that I (as expressed in power/meter^2) which is proportional to E, transforms as the square of the above factor.
If one wants to calculate when the sun is swamped out by black body radiation , one can probably estimate it by looking at the amount of energy radiated by CMB by the solid angle resolution of the human eye (this should be proportioal to the angular resolution squared, I'm not sure if there is any proportionality constant or not). One then compares this with the amount of energy emitted by the sun itself, which acts as a point source, that we calculated above.