加速 Octave / Matlab 绘图

Question

Gnoivce 和 Hartmut 为这段代码提供了很多帮助，但运行需要一段时间。 bar 命令中的 CData 属性在我使用的 Octave 4.0-4.2.1 版本中似乎没有实现。解决此问题的方法是单独绘制所有单个条形图并为每个单独的条形图设置单独的颜色。人们帮助了我，让我走到了这一步，但是情节显示需要 5 分钟 有人知道加快速度的方法吗？

以下代码运行s： marbles.jpg 下面使用的图像文件：

clear all,clf reset,tic,clc 
rgbImage = imread('/tmp/marbles.jpg');
hsvImage = rgb2hsv(rgbImage);         % Convert the image to HSV space
hPlane = 360.*hsvImage(:, :, 1);      % Get the hue plane scaled from 0 to 360
binEdges = 0:360;                 %# Edges of histogram bins
N = histc(hPlane(:),binEdges);    %# Bin the pixel hues from above

C = colormap(hsv(360));  %# create an HSV color map with 360 points
stepsize = 1; % stepsize 1 runs for a while...

for n=binEdges(2:stepsize:end)   %# Plot the histogram, one bar each
    if (n==1), hold on, end
    h=bar(n,N(n));
    set(h,'FaceColor',C(n,:));   %# set the bar color individually
end

axis([0 360 0 max(N)]);           %# Change the axes limits
set(gca,'Color','k');             %# Change the axes background color
set(gcf,'Pos',[50 400 560 200]);  %# Change the figure size
xlabel('HSV hue (in degrees)');   %# Add an x label
ylabel('Bin counts');             %# Add a y label
fprintf('\nfinally Done-elapsed time -%4.4fsec- or -%4.4fmins- or -%4.4fhours-\n',toc,toc/60,toc/3600);

5 分钟后创建的绘图：

查看原题

Answer 1

我猜循环是您代码中花费这么长时间的瓶颈？您可以删除循环并通过调用 bar, then call set 修改 hggroup 对象及其子对象 patch 来创建绘图：

h = bar(binEdges(1:end-1), N(1:end-1), 'histc');   % hggroup object
set(h, 'FaceColor', 'flat', 'EdgeColor', 'none');
hPatch = get(h, 'Children');                       % patch object
set(hPatch, 'CData', 1:360, 'CDataMapping', 'direct');

在 Octave 4.0.3 中使用此修复程序重复您的代码会立即为我呈现：

Answer 2

正如我在评论中建议的那样，我会使用图像（在我的系统上为您的图像花费 0.12 秒）。

编辑：更多评论，修复小错误，允许创建步长 > 1

的容器

img_fn = "17S9PUK.jpg";
if (! exist (img_fn, "file"))
  disp ("downloading image from imgur.com...");
  fflush (stdout);
  urlwrite ("http://i.imgur.com/17S9PUK.jpg", "17S9PUK.jpg");
endif

rgbImage = imread (img_fn);

## for debugging so the matrixes fit on screen
if (0)
  pkg load image
  rgbImage = imresize (rgbImage, [6 8]);
endif

hsvImage = rgb2hsv(rgbImage);
hPlane = 360 .* hsvImage(:, :, 1); 

## create bins, I've choosen 2 step to "smooth" the result
binEdges = 1:2:360;
N = histc (hPlane(:), binEdges)';

cm = permute (hsv (numel (binEdges)), [3 1 2]);

## Create an image with x = hue
img = repmat (cm, max(N), 1);

## Create binary mask which is used to black "img" dependent on N
sp = sparse (N(N > 0), (1:360)(N > 0), true, max(N), numel (binEdges));
mask = full (cumsum (flipud (sp)));

## extend mask in depth to suppress RGB
mask = repmat (mask, [1 1 3]);

## use inverted mask to "black out" pixels < N
img(logical (1 - flipud (mask))) = 0;

## show image
image (binEdges, 1:max(N), img)
set (gca, "ydir", "normal");
xlabel('HSV hue (in degrees)');
ylabel('Bin counts');

## print it for Whosebug
print ("out.png")

与上面相同，但 bin 宽度为 1（运行时间为 0.167423 秒。）

binEdges = 1:360;

加速 Octave / Matlab 绘图

Speeding up Octave / Matlab plot

matlab

plot

octave