Simulating intensity from surface plot?

3 次查看（过去 30 天）
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!

类别

Help CenterFile Exchange 中查找有关 Surface and Mesh Plots 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by