quiver plot on an image
显示 更早的评论
Hello
I am trying to quiverplot on an image, but quiver plot shows up very small on the upper left corner (attached here)
I have to use 'imshow' instead of 'imagesc' becasue I want the image in black and white and colored quiver vectors. When I assigne colormap for 'imagesc' it messes the quiver vectors.
Is there any way to superimpose quiver plot on imshow function?
thank you
采纳的回答
"Is there any way to superimpose quiver plot on imshow function?"
Yes.
% read an image:
img = imread('peppers.png');
% imshow the image:
imshow(img)
% make a quiver plot:
hold on
[m,n,~] = size(img);
[x,y] = meshgrid(1:n,1:m);
quiver(x,y,randn(m,n),randn(m,n))
% zoom in to see some details:
xlim([101 115])
ylim([101 110])
It looks like the x and y you use in quiver are incorrect is why your quiver plot is so small.
4 个评论
Reza Zh
2024-3-1
Here's an example where the image and x and y are the same sizes as yours:
% an image the size of your image:
img = imresize(imread('peppers.png'),[2304 2401]);
size(img)
figure
im = imshow(img);
By default the XData of a plotted image is [1 size(img,2)] and the YData is [1 size(img,1)], so the pixel coordinates are 1:size(img,2) in the x-dimension and 1:size(img,1) in the y-dimension.
im.XData
im.YData
% some arbitrary x and y of the same size as yours:
x = linspace(2,78,31); % x goes from 2 to 78
y = linspace(3,14,29); % y goes from 3 to 14
[x,y] = meshgrid(x,y);
size(x)
You either have to (1) scale your x to span from 1 to size(img,2) and scale your y to span from 1 to size(img,1), or (2) you have to set your image's XData and YData to match the limits of your quiver's x and y.
Option 1: rescale the x and y to the location of the image:
minx = min(x(:));
maxx = max(x(:));
miny = min(y(:));
maxy = max(y(:));
[m,n,~] = size(img);
X = (x-minx)./(maxx-minx).*(n-1)+1;
Y = (y-miny)./(maxy-miny).*(m-1)+1;
Check the min and max of X and Y (these match the image now):
min(X(:))
max(X(:))
min(Y(:))
max(Y(:))
% make a quiver plot:
hold on
quiver(X,Y,randn(size(x)),randn(size(x)))
Option 2: adjust the image XData and YData:
% new figure:
figure
% show the image in the new figure:
im = imshow(img);
% now adjust the XData and YData (makes the image tiny in this case):
im.XData = [min(x(:)) max(x(:))];
im.YData = [min(y(:)) max(y(:))];
% make a quiver plot:
hold on
quiver(x,y,randn(size(x)),randn(size(x)))
% the axes limits are still for the image's original
% XData and YData, so adjust those too:
set(gca(),'XLimMode','auto','YLimMode','auto') % setting the limit modes to auto
% xlim(im.XData+[-0.5 0.5]) % alternatively, explicitly set the limits to
% ylim(im.YData+[-0.5 0.5]) % match the new XData and YData
Reza Zh
2024-3-1
Voss
2024-3-1
You're welcome!
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Image Processing Toolbox 的更多信息
Translated by 
