another question has come up about hist3()

Question

jgillis16 2015-6-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/225498-another-question-has-come-up-about-hist3

评论： Star Strider 2015-7-5

OK. So, I have access to the statistics toolbox now. And, I was wondering how I would go about plotting the hist3 function for the same code from:

http://www.mathworks.com/matlabcentral/answers/224847-how-to-create-bins-that-are-then-used-to-create-surf-function

Because [N,C]=hist3(...) doesn't actually display a plot, and having something along the lines of the Seamount:Data Point Density Histogram and Intensity plot (link is below) would be extremely helpful.

http://www.mathworks.com/help/stats/hist3.html?refresh=true

Thanks again!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2015-6-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/225498-another-question-has-come-up-about-hist3#answer_184161

在 MATLAB Online 中打开

This works, but I’m not terribly happy with the result:

[N,C] = hist3([C12 S12], {0:359; -89:89});
figure(1)
bar3(N)
grid on
xtk = linspace(-89,89, 11);
set(gca, 'XTick',xtk+90,'XTickLabel',xtk)
xlabel('DEC (-89 +89)')
ylabel('RA (1 360)')
axis tight

11 个评论
显示 9更早的评论隐藏 9更早的评论

Star Strider 2015-6-28

在 MATLAB Online 中打开

I prefer to use mesh rather than bar3, although for a histogram, bar3 is more appropriate.

With this code:

fid = fopen( 'jgillis16 GWGCCatalogrm.txt');
trashLine = fgets(fid); %Skips the first line
data = textscan(fid, '%f%s%f%f%s%f%f%f%f%f%f%f%f%f%f%f%f%f', 'Delimiter', '|', 'TreatAsEmpty','~');
fclose(fid);
M1=data;
GalList.pgc = data{1};
GalList.name = data{2};
GalList.ra = data{3};
GalList.dec = data{4};
GalList.major = data{7};
GalList.abs_mag = data{14};
GalList.dist = data{15};
GalList.mass =(10.^( (-20.8 - GalList.abs_mag)./ 2.5 )) ; % Converts Absolute Magnitude into (blue luminosity) mass in units of Milky Way (luminosity) mass (how many potential forming stars in region using the blue luminosity which is rough estimate of mass).
GalList.ra = GalList.ra.*15;
GalList.dec = GalList.dec ;
GalaxyList = GalList;
C1 = GalList.ra;
S1 = GalList.dec;
C12 = round(C1);
S12 = round(S1);
[N,C] = hist3([C12 S12], {0:359; -89:89});
figure(1)
mesh(N)
grid on
xtk = linspace(-89,89, 11);
set(gca, 'XTick',xtk+90,'XTickLabel',xtk)
xlabel('DEC (-89 +89)')
ylabel('RA (1 360)')
axis tight

I get this plot:

Star Strider 2015-7-1

在 MATLAB Online 中打开

Your data are not easy to analyse. This time, I believe it’s due to inconsistencies between your two data files and my ignorance of astronomy.

I cannot get the DEC to display correctly on these plots. I invite you to look over this code to see if you can see the problem.

The code works. The problem is likely that I don’t understand your data and how to work with it. Please be specific as to what data from what files you want to analyse.

fid = fopen( 'jgillis16 NewGalaxy.txt');
trashLine = fgets(fid); %Skips the first line
data = textscan(fid, '%f%s%f%f%s%f%f%f%f%f%f%f%f%f%f%f%f%f', 'Delimiter', '|', 'TreatAsEmpty','~');
fclose(fid);
M2=data;
C12_NG = round(M2{3});
S12_NG = round(M2{4});
fid = fopen( 'jgillis16 GWGCCatalogrm.txt');
trashLine = fgets(fid); %Skips the first line
data = textscan(fid, '%f%s%f%f%s%f%f%f%f%f%f%f%f%f%f%f%f%f', 'Delimiter', '|', 'TreatAsEmpty','~');
fclose(fid);
M1=data;
C12_GW = round(M1{3});
S12_GW = round(M1{4});
[NNG,CNG] = hist3([C12_NG S12_NG], {0:359; -89:89});    % NewGalaxy
[NGW,CGW] = hist3([C12_GW S12_GW], {0:359; -89:89});    % GWGCCatalogrm
Add = NNG + NGW;                                        % Summed Histograms
Rat = NNG ./ NGW;                                       % Ratio: NG/GW
Pct = NNG ./ Add;                                       % Percent: NG/Total
Q1 = CNG{1};                                            % Check Values
Q2 = CNG{2};                                            % Check Values
[XCtrs, YCtrs] = meshgrid(CNG{2}, CNG{1});
figure(1)
surf(XCtrs, YCtrs, Pct)
grid on
xtk = linspace(-89,89, 9);
set(gca, 'XTick',xtk,'XTickLabel',xtk)
ytk = linspace(0, 359, 10);
set(gca, 'YTick',ytk, 'YTickLabel',ytk)
qxt = get(gca, 'XTick');
qyt = get(gca, 'YTick');
xlabel('DEC (-89 +89)')
ylabel('RA (1 360)')
zlabel('Fraction:  NewGalaxy/Total')
axis tight
shading interp
Q3 = XCtrs(1:36:360, 1:18:179);                         % Check Values
Q4 = YCtrs(1:36:360, 1:18:179);                         % Check Values
NNG = accumarray([C12_NG+1, S12_NG+90], 1);
NGW = accumarray([C12_GW+1, S12_GW+90], 1);
Add = NNG + NGW;                                        % Summed Histograms
Rat = NNG ./ NGW;                                       % Ratio: NG/GW
Pct = NNG ./ Add;                                       % Percent: NG/Total
figure(2)
surf(Pct)
grid on
xtk = linspace(-89,89, 9);
set(gca, 'XTick',xtk,'XTickLabel',xtk)
ytk = linspace(0, 359, 10);
set(gca, 'YTick',ytk, 'YTickLabel',ytk)
qxt = get(gca, 'XTick');
qyt = get(gca, 'YTick');
xlabel('DEC (-89 +89)')
ylabel('RA (1 360)')
zlabel('Fraction:  NewGalaxy/Total')
axis tight
shading interp

jgillis16 2015-7-5

OK. Could you perhaps attach the plot you produce with the updated code provided?

I am moving the morphology question to another string: http://it.mathworks.com/matlabcentral/answers/228451-applying-a-fourth-variable-to-a-stem3-plot

Star Strider 2015-7-5

I would, but since there are discrepancies between the corresponding values in the two files, it doesn’t make any sense. That’s the reason I didn’t post it here. If you look though your files to be sure the declinations and right-ascensions all span the same ranges (-89 +89 and 0 359 as appropriate), I will run my code and plot it.

I posted an Answer to your other Question. It tangentially relates to this one.

请先登录，再进行评论。