Thread Subject: how to plot 2D histogram of medical image

I have two a images T1 and T2 from magnetic resonans session in DICOM standard and I have to plot 2-dimensional histogram off this image.I tried to do it gor example like that:
[A B] = meshgrid(1:N,1:N);
figure;
plot3(A,B,image1,'. k','MarkerSize',3);
hold on;
plot3(A,B,image2,'. k','MarkerSize',3);

but it's 3D histogram

I tried also something like that:

figure;
load seamount;
wym=[N,M];
hold on;
hist3(wym);
n = hist3(wym);
n1=n;
n1(size(n,1)+1,size(n,2)+1)=0;
xb=linspace(min(wym(:,1)),max(wym(:,1)),size(n,1)+1);
yb=linspace(min(wym(:,2)),max(wym(:,2)),size(n,1)+1);
h=pcolor(xb,yb,n1);
set(h,'zdata',obraz1(size(n1))*max(max(n)))
colormap(hot)
grid on
view(3);

but it also 3D histogram and I don't know how to do it in 2D...

On Oct 27, 12:42 pm, "Magda Kie?basa" <magda13...@wp.pl> wrote:
> I have two a images T1 and T2 from magnetic resonans session in DICOM standard and I have to plot 2-dimensional histogram off this image.I tried to do it gor example like that:
> [A B] = meshgrid(1:N,1:N);
> figure;
> plot3(A,B,image1,'. k','MarkerSize',3);
> hold on;
> plot3(A,B,image2,'. k','MarkerSize',3);
>
> but it's 3D histogram
>
> I tried also something like that:
>
> figure;
> load seamount;
> wym=[N,M];
> hold on;
> hist3(wym);
> n = hist3(wym);
> n1=n;
> n1(size(n,1)+1,size(n,2)+1)=0;
> xb=linspace(min(wym(:,1)),max(wym(:,1)),size(n,1)+1);
> yb=linspace(min(wym(:,2)),max(wym(:,2)),size(n,1)+1);
> h=pcolor(xb,yb,n1);
> set(h,'zdata',obraz1(size(n1))*max(max(n)))
> colormap(hot)
> grid on
> view(3);
>
> but it also 3D histogram and I don't know how to do it in 2D...

Do a search on histxy in this Usenet group.
You will find a couple of algorithms there that date back to 2002.