Simulating intensity from surface plot?

So I'm trying to essentially show what the intensity distribution from a given 3D object would be, and have an intensity plot produced, but this is proving tricky.

given the following code:

then = linspace(0,3*pi/4,100);
phi = linspace(0,2*pi,100);
thes = linspace(0,-3*pi/4,100);
r = 21;
h = linspace(0,29 - 21/sqrt(2),55);
rc = linspace(0,21/sqrt(2),55);

[Tn,Pn] = meshgrid(then,phi);

ln = r.*sin(Tn).*cos(Pn);
yn = r.*sin(Tn).*sin(Pn);
bn = r.*cos(Tn) + 29;

[Ts,Ps] = meshgrid(thes,phi);

ls = r.*sin(Ts).*cos(Ps);
ys = r.*sin(Ts).*sin(Ps);
bs = -r.*cos(Ts) - 29;

[RC,PH] = meshgrid(rc,phi);
[H,PH] = meshgrid(h,phi);

lc = RC.*cos(PH);
yc = RC.*sin(PH);
bcn = H;
bcs = -H;

surf(ln,yn,bn,'FaceColor','m')

hold on

surf(ls,ys,bs,'FaceColor','m')

surf(lc,yc,bcn,'FaceColor','m')

surf(lc,yc,bcs,'FaceColor','m')

axis image

set(gca, 'XDir', 'reverse')
set(gca, 'FontSize', 16)

xlabel('Longitude (l) (degrees)')
ylabel('Depth (degrees)')
zlabel('Latitude (b) (degrees)')
title('Idealised, spherically symmetric Fermi bubbles')

camlight(180,0)

view(0,0)

which produces the following plot:

Is there any way in which I can make a plot showing what the intensity from this object would be, looking from this angle? Assuming the bubbles are uniform (in all ways, i.e. luminosity and specifically, surface brightness). I would be expecting the plot of the intensity to look Gaussian for example.

Fiddling around with functions like surfc hasn't yielded any results as of yet. Looking forward to hearing back from you. Thanks!