我想在cocoa应用程序中显示图像RGB颜色直方图。请建议使用目标c或任何可用于实现此目的的第三方库来实现此目的。
答案 0 :(得分:2)
这是一个问题,因为 RGB 颜色是 3D 空间所以他们的直方图会导致 4D 情节,这是我们真的没有理解。
因此,解决方案是将 4D 绘图以某种方式转换为 3D 绘图。这可以通过用具有某种含义的东西对颜色进行分类来完成。我不会推测和描述我正在使用的东西。我使用 HSV 色彩空间并忽略 V 值。这样我就失去了很多色调信息,但它仍然足以为我的目的描述颜色。这就是它的样子:
您还可以使用不同V
的更多绘图来覆盖更多颜色。有关详细信息,请参阅:
无论如何,你可以使用任何渐变排序或任何形状完整的你的情节。
如果您想要纯 RGB ,那么您可以调整它并使用 RGB 立方体表面或将其映射到球体上并忽略(0,0,0)
的长度(使用单位)矢量)这样的事情:
因此,如果R,G,B
位于<0,1>
,则将其转换为<-1,+1>
,然后计算球面坐标(忽略半径),并获得2个变量,而不是3个可用作情节( 2D 地球基地或 3D 球体......)。
这里的C ++代码如何执行此操作(由HSV直方图制作):
picture pic0,pic1,pic2,zed;
const int na=360,nb=180,nb2=nb>>1; // size of histogram table
int his[na][nb];
DWORD w;
int a,b,r,g,x,y,z,l,i,n;
double aa,bb,da,db,dx,dy,dz,rr;
color c;
pic2=pic0; // copy input image pic0 to pic2
for (a=0;a<na;a++) // clear histogram
for (b=0;b<nb;b++)
his[a][b]=0;
for (y=0;y<pic2.ys;y++) // compute it
for (x=0;x<pic2.xs;x++)
{
c=pic2.p[y][x];
r=c.db[picture::_r]-128;
g=c.db[picture::_g]-128;
b=c.db[picture::_b]-128;
l=sqrt(r*r+g*g+b*b); // convert RGB -> spherical a,b angles
if (!l) { a=0; b=0; }
else{
a=double(double(na)*acos(double(b)/double(l))/(2.0*M_PI));
if (!r) b=0; else b=double(double(nb)*atan(double(g)/double(r))/(M_PI)); b+=nb2;
while (a<0) a+=na; while (a>=na) a-=na;
if (b<0) b=0; if (b>=nb) b=nb-1;
}
his[a][b]++; // update color usage count ...
}
for (n=0,a=0;a<na;a++) // max probability
for (b=0;b<nb;b++)
if (n<his[a][b]) n=his[a][b];
// draw the colored RGB sphere and histogram
zed =pic1; zed .clear(9999); // zed buffer for 3D
pic1.clear(0); // image of histogram
da=2.0*M_PI/double(na);
db=M_PI/double(nb);
for (aa=0.0,a=0;a<na;a++,aa+=da)
for (bb=-M_PI,b=0;b<nb;b++,bb+=db)
{
// normal
dx=cos(bb)*cos(aa);
dy=cos(bb)*sin(aa);
dz=sin(bb);
// color of surface (darker)
rr=75.0;
c.db[picture::_r]=double(rr*dx)+128;
c.db[picture::_g]=double(rr*dy)+128;
c.db[picture::_b]=double(rr*dz)+128;
c.db[picture::_a]=0;
// histogram center
x=pic1.xs>>1;
y=pic1.ys>>1;
// surface position
rr=64.0;
z=rr;
x+=double(rr*dx);
y+=double(rr*dy);
z+=double(rr*dz);
if (zed.p[y][x].dd>=z){ pic1.p[y][x]=c; zed.p[y][x].dd=z; }
// ignore lines if zero color count
if (!his[a][b]) continue;
// color of lines (bright)
rr=125.0;
c.db[picture::_r]=double(rr*dx)+128;
c.db[picture::_g]=double(rr*dy)+128;
c.db[picture::_b]=double(rr*dz)+128;
c.db[picture::_a]=0;
// line length
l=(xs*his[a][b])/(n*3);
for (double xx=x,yy=y,zz=z;l>=0;l--)
{
if (zed.p[y][x].dd>=z){ pic1.p[y][x]=c; zed.p[y][x].dd=z; }
xx+=dx; yy+=dy; zz+=dz; x=xx; y=yy; z=zz;
if (x<0) break; if (x>=xs) break;
if (y<0) break; if (y>=ys) break;
}
}
pic0
,输出图像为pic1
(直方图)pic2
是pic0
(旧代码的残余)zed
是用于3D显示的Zed缓冲区,避免了Z排序...... 我将自己的图片类用于图片,因此有些成员是:
xs,ys
图片大小(以像素为单位)
p[y][x].dd
是(x,y)位置的像素,为32位整数类型
clear(color)
- 清除整个图像
resize(xs,ys)
- 将图片大小调整为新分辨率
由于球体是3D物体,您应该向其添加旋转,以便所有表面都能及时显示(或使用鼠标或其他任何方式旋转)......