transforming matrix after change of variables

1 次查看(过去 30 天)
Hi,
I have a matrix M(r, t), which I would like to transform into a matrix N(x, y, t), where x = r*cos(theta), y = r*sin(theta). I would like to use this to generate a heat map using image() where it plots x against y with the color corresponding to the entry in the matrix. Is there a way to do this?
Thanks for the help!
  2 个评论
the cyclist
the cyclist 2015-11-24
What is stored in M? The value at a given (r,t) location? Is it a column vector? Do you have the corresponding values of r and t stored in corresponding vectors?
María Jesús
María Jesús 2015-11-24
Yes, M has a value at each (r,t), it is a 100x100 array. r and t are each a 1x100 array

请先登录,再进行评论。

回答(2 个)

the cyclist
the cyclist 2015-11-24
Assuming you are using t and theta interchangeably here, then I think it should just be
x = r.*cos(t);
y = r.*sin(t);
figure
imagesc(x,y,M)
  1 个评论
María Jesús
María Jesús 2015-11-24
They're meant to be different. Theta runs from 0 to 2pi and t from 0 to 0.2, although the vectors do have the same size.

请先登录,再进行评论。


Walter Roberson
Walter Roberson 2015-11-25
image() and imagesc() and imshow() only deal with rectangular matrices aligned with the X/Y axes. When you pass x or y coordinates into image() or imagesc(), only the first and last values are paid attention to.
If you want an image with nonlinear coordinates you have a few choices:
  1. pcolor() or surf() with the coordinate grids and the image data. pcolor() is exactly equivalent to surf() viewed from above. Note that the coordinates are taken to be of vertices and color is interpolated based upon the color of the vertices
  2. use patch() in texturemap mode to draw the image stretched on to the surface
  3. use a gridded interpolant and sample the function at a rectangular grid that is then an image that can be used with image() or imagesc()
  4. probably not for this case but in some cases, use one of the image transforms
  5. you might be able to use iradon()

类别

Help CenterFile Exchange 中查找有关 Data Distribution Plots 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by