% void CLASS cam_xyz_coeff (float rgb_cam[3][4], double cam_xyz[4][3])
% {
% double cam_rgb[4][3], inverse[4][3], num;
% int i, j, k;
%
% for (i=0; i < colors; i++) /* Multiply out XYZ colorspace /
% for (j=0; j < 3; j++)
% for (cam_rgb[i][j] = k=0; k < 3; k++)
% cam_rgb[i][j] += cam_xyz[i][k]
xyz_rgb[k][j];
%
% for (i=0; i < colors; i++) { /* Normalize cam_rgb so that /
% for (num=j=0; j < 3; j++) /
cam_rgb * (1,1,1) is (1,1,1,1) */
% num += cam_rgb[i][j];
% for (j=0; j < 3; j++)
% cam_rgb[i][j] /= num;
% pre_mul[i] = 1 / num;
% }
% pseudoinverse (cam_rgb, inverse, colors);
% for (i=0; i < 3; i++)
% for (j=0; j < colors; j++)
% rgb_cam[i][j] = inverse[j][i];
% }

% void CLASS pseudoinverse (double (in)[3], double (out)[3], int size)
% {
% double work[3][6], num;
% int i, j, k;
%
% for (i=0; i < 3; i++) {
% for (j=0; j < 6; j++)
% work[i][j] = j == i+3;
% for (j=0; j < 3; j++)
% for (k=0; k < size; k++)
% work[i][j] += in[k][i] * in[k][j];
% }
% for (i=0; i < 3; i++) {
% num = work[i][i];
% for (j=0; j < 6; j++)
% work[i][j] /= num;
% for (k=0; k < 3; k++) {
% if (k==i) continue;
% num = work[k][i];
% for (j=0; j < 6; j++)
% work[k][j] -= work[i][j] * num;
% }
% }
% for (i=0; i < size; i++)
% for (j=0; j < 3; j++)
% for (out[i][j]=k=0; k < 3; k++)
% out[i][j] += work[j][k+3] * in[i][k];
% }

% const double xyz_rgb[3][3] = { /* XYZ from RGB */
% { 0.412453, 0.357580, 0.180423 },
% { 0.212671, 0.715160, 0.072169 },
% { 0.019334, 0.119193, 0.950227 } };
clc
clear all
cmsrc0 = [1.895316, -0.688256,-0.20706;
-0.180953,1.625399,-0.444446;
-0.0000147,-0.524012,1.524159;
0,0,0]

cam_xyz = [0.617188, -0.140625, -0.0390625;
-0.460938, 1.19531, 0.226563;
-0.0859375, 0.1953113, 0.545875;
0,0,0
];

cam_rgb = zeros(4,3);
inverse = zeros(4,3);
xyz_rgb = [0.412453, 0.357580, 0.180423,0;
0.212671, 0.715160, 0.072169,0;
0.019334, 0.119193, 0.950227,0];
cam_rgb = cam_xyz * xyz_rgb;
for i =1:3
cam_rgb(i,:) = cam_rgb(i,:) / sum(cam_rgb(i,:));
end
cam_rgb = cam_rgb(:,1:3);
colors = 3;

x11 = cam_rgb(1,1);
x12 = cam_rgb(1,2);
x13 = cam_rgb(1,3);
x21 = cam_rgb(2,1);
x22 = cam_rgb(2,2);
x23 = cam_rgb(2,3);
x31 = cam_rgb(3,1);
x32 = cam_rgb(3,2);
x33 = cam_rgb(3,3);
work = zeros(3,6);
work(1,1) = x11 x11 + x21 x21 + x31 * x31;
work(1,2) = x11 x12 + x21 x22 + x31 * x32;
work(1,3) = x11 x13 + x21 x23 + x31 * x33;
work(2,1) = x11 x12 + x21 x22 + x31 * x32;
work(2,2) = x12 x12 + x22 x22 + x32 * x32;
work(2,3) = x12 x13 + x22 x23 + x32 * x33;
work(3,1) = x11 x13 + x21 x23 + x31 * x33;
work(3,2) = x12 x13 + x22 x23 + x32 * x33;
work(3,3) = x13 x13 + x23 x23 + x33 * x33;
work(:,4:6) = [1,0,0;0,1,0;0,0,1];

% i =1
num1 = work(1,1);
work(1,1) = work(1,1) / num1;
work(1,2) = work(1,2) / num1;
work(1,3) = work(1,3) / num1;
work(1,4) = work(1,4) / num1;
work(1,5) = work(1,5) / num1;
work(1,6) = work(1,6) / num1;

% k = 2
num2 = work(2,1);
work(2,1) = work(2,1) - num2 * work(1,1);
work(2,2) = work(2,2) - num2 * work(1,2);
work(2,3) = work(2,3) - num2 * work(1,3);
work(2,4) = work(2,4) - num2 * work(1,4);
work(2,5) = work(2,5) - num2 * work(1,5);
work(2,6) = work(2,6) - num2 * work(1,6);

% k = 3
num3 = work(3,1);
work(3,1) = work(3,1) - num3 * work(1,1);
work(3,2) = work(3,2) - num3 * work(1,2);
work(3,3) = work(3,3) - num3 * work(1,3);
work(3,4) = work(3,4) - num3 * work(1,4);
work(3,5) = work(3,5) - num3 * work(1,5);
work(3,6) = work(3,6) - num3 * work(1,6);

% i = 2
num1 = work(2,2);
work(2,1) = work(2,1) / num1;
work(2,2) = work(2,2) / num1;
work(2,3) = work(2,3) / num1;
work(2,4) = work(2,4) / num1;
work(2,5) = work(2,5) / num1;
work(2,6) = work(2,6) / num1;

num2 = work(1,2);
work(1,1) = work(1,1) - num2 * work(2,1);
work(1,2) = work(1,2) - num2 * work(2,2);
work(1,3) = work(1,3) - num2 * work(2,3);
work(1,4) = work(1,4) - num2 * work(2,4);
work(1,5) = work(1,5) - num2 * work(2,5);
work(1,6) = work(1,6) - num2 * work(2,6);

num3 = work(3,2);
work(3,1) = work(3,1) - num3 * work(2,1);
work(3,2) = work(3,2) - num3 * work(2,2);
work(3,3) = work(3,3) - num3 * work(2,3);
work(3,4) = work(3,4) - num3 * work(2,4);
work(3,5) = work(3,5) - num3 * work(2,5);
work(3,6) = work(3,6) - num3 * work(2,6);

% i = 3
num1 = work(3,3);
work(3,1) = work(3,1) / num1;
work(3,2) = work(3,2) / num1;
work(3,3) = work(3,3) / num1;
work(3,4) = work(3,4) / num1;
work(3,5) = work(3,5) / num1;
work(3,6) = work(3,6) / num1;

num2 = work(1,3);
work(1,1) = work(1,1) - num2 * work(3,1);
work(1,2) = work(1,2) - num2 * work(3,2);
work(1,3) = work(1,3) - num2 * work(3,3);
work(1,4) = work(1,4) - num2 * work(3,4);
work(1,5) = work(1,5) - num2 * work(3,5);
work(1,6) = work(1,6) - num2 * work(3,6);

num3 = work(2,3);
work(2,1) = work(2,1) - num3 * work(3,1);
work(2,2) = work(2,2) - num3 * work(3,2);
work(2,3) = work(2,3) - num3 * work(3,3);
work(2,4) = work(2,4) - num3 * work(3,4);
work(2,5) = work(2,5) - num3 * work(3,5);
work(2,6) = work(2,6) - num3 * work(3,6);

out = zeros(3,3);
out(1,1) = work(1,4) * cam_rgb(1,1) + work(1,5) * cam_rgb(1,2) + work(1,6) * cam_rgb(1,3);
out(1,2) = work(2,4) * cam_rgb(1,1) + work(2,5) * cam_rgb(1,2) + work(2,6) * cam_rgb(1,3);
out(1,3) = work(3,4) * cam_rgb(1,1) + work(3,5) * cam_rgb(1,2) + work(3,6) * cam_rgb(1,3);
out(2,1) = work(1,4) * cam_rgb(2,1) + work(1,5) * cam_rgb(2,2) + work(1,6) * cam_rgb(2,3);
out(2,2) = work(2,4) * cam_rgb(2,1) + work(2,5) * cam_rgb(2,2) + work(2,6) * cam_rgb(2,3);
out(2,3) = work(3,4) * cam_rgb(2,1) + work(3,5) * cam_rgb(2,2) + work(3,6) * cam_rgb(2,3);
out(3,1) = work(1,4) * cam_rgb(3,1) + work(1,5) * cam_rgb(3,2) + work(1,6) * cam_rgb(3,3);
out(3,2) = work(2,4) * cam_rgb(3,1) + work(2,5) * cam_rgb(3,2) + work(2,6) * cam_rgb(3,3);
out(3,3) = work(3,4) * cam_rgb(3,1) + work(3,5) * cam_rgb(3,2) + work(3,6) * cam_rgb(3,3);
out'

% work = zeros(3,6);
% for i =1:3
% for j = 1:6
% if j == i+3
% work(i,j) = 1;
% else
% work(i,j) = 0;
% end
% end
%
% for j = 1:3
% for k = 1:colors
% work(i,j) = work(i,j) + cam_rgb(k,i) * cam_rgb(k,j);
% end
% end
% end
%
% for i =1:3
% num = work(i,i);
% for j = 1:6
% work(i,j) = work(i,j) / num;
% end
% for k = 1:3
% if(k == i)
% continue;
% end
% num = work(k,i);
% for j = 1:6
% work(k,j) = work(k,j) - work(i,j) * num;
% end
% end
% end

% out = zeros(3,3);
% for i = 1:colors
% for j =1:3
% out(i,j) =0;
% for k = 1:3
% out(i,j) = out(i,j) + work(j,k+3) * cam_rgb(i,k);
% end
% end
% end
%
% out = out'

ceshi的更多相关文章

  1. http://www.spasvo.com/ceshi/open/kyxncsgj/Jmeter/

    http://www.spasvo.com/ceshi/open/kyxncsgj/Jmeter/

  2. 解决Swap file &quot;.ceshi.c.swp&quot; already exists!问题

    关于swp文件:使用vi,常常能够看到swp这个文件,那这个文件是怎么产生的呢.当你打开一个文件,vi就会生成这么一个.(filename)swp文件以备不測,假设你正常退出,那么这个.(filena ...

  3. Fatal error: Call to undefined function Think\C() in /var/www/html/ceshi.hzheee.com/think/ThinkPHP/Library/Think/Think.class.php on line 334 这个问题解决

    当APP_DEBUG为true时,包含图中这个文件,文件中又引导包含这些库文件,可以看出安装thinkphp3.2.3时ThinkPHP/Common/下是functions.php,把它改成func ...

  4. laotech老师唠科mac 深入浅出MAC OS X ceshi ruguokeyi

    laotech老师唠科mac 深入浅出MAC OS X http://study.163.com/plan/planLearn.htm?id=1637004#/learn/resVideo?lesso ...

  5. Jquery的点击事件,三句代码完成全选事件

    先来看一下Js和Jquery的点击事件 举两个简单的例子 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN&q ...

  6. Jquery(1)

    鼠标点击事件: <input type="button" value="测试" onclick="test()" /> < ...

  7. JQuery实现表格的增加行和删除行

    利用JQuery实现datatables插件的增加和删除行操作 在学习过程中遇到了这个利用JQuery对表格行的增加和删除,特记录下来以供初学者参考. 下面是主要的代码: <meta http- ...

  8. Jquary入门(添加 修改 表单元素)+ JSON+弹框

    字符串拼接 计算机语言 都是 对  数据的处理(获取/修改数据)  添加元素  除了  固定的方法添加 以外 都是   字符串拼接(拼接成固定格式即可执行). 1.表单添加元素  append() 方 ...

  9. Jquery事件

    在JS里加事件  <input type="button" value="测试" onclick="test()"/>跟一个函数 ...

随机推荐

  1. 中国计算机学会CCF推荐国际学术会议

    中国计算机学会推荐国际学术会议 (计算机系统与高性能计算) 一.A类 序号 会议简称 会议全称 出版社 网址 1 ASPLOS Architectural Support for Programmin ...

  2. 一个ListView怎么展示两种样式

    private class MyBaseMsgAdapter extends BaseAdapter { //获取数据适配器中条目类型的总数,修改成两种(纯文本,输入+文字) @Override pu ...

  3. C# 深入浅出 异步(八)

    C#异步调用学习链接:从C#5.0说起:再次总结C#异步调用方法发展史

  4. 函数式 js 接口实现原理,以及 lodash/fp 模块

    函数式 js 接口 之前在 youtube 上看到一个技术视频,讲“underscore.js的接口为什么不好用”,以及什么样的接口更好用.演讲者是 lodash.js 的作者,他提出了一种“全面函数 ...

  5. linux -目录结构

    摘自:http://www.comptechdoc.org/os/linux/usersguide/linux_ugfilestruct.html 这个目录结构介绍是我目前看到介绍最全的,有时间在翻译 ...

  6. JS判断是不是手机浏览器浏览网站的网页,并自动跳转

    现在智能手机上网越来越普遍了,为了获得用户体验增加网站流量,你有必要为你的网站增加一个访问端设备的判断功能,若发现是手机用户访问,则直接跳转到手机站,通过百度的APP site,很容易就可实现这功能. ...

  7. Android6.0权限组申请

    void checkPermission() { final List<String> permissionsList = new ArrayList<>(); if (Bui ...

  8. C++指针参数引用

    粘个代码占位置,以后有时间把指针函数,函数指针都补上 #include <iostream> using namespace std; void freePtr1(int* p1){ /* ...

  9. Owin Self Host

    http://owin.org/ Owin 定义了webserver和webapplication之间的标准接口,目标就是为了解耦webapplication对webserver的依赖, 就是说以后可 ...

  10. Python~字典

    if not isinstance(x, (int, float)): raise TypeError('bad operand type')   range() raw_input(‘birth’) ...