suncn-CSDN博客

转载直线拟合之最小二乘法 Ax+By+C=0

最小二乘法拟合直线，对于大多数人都不是很陌生，直线方程 y=kx+b，令sum((y-kx-b)^2)，分别对k和b求导便可求出最佳的参数，但是如果遇到诸如“X=5”这样的直线方程，又该如何呢？这时候就不再适用了，参考如下描述http://mathforum.org/library/drmath/view/51809.html，这样拟合出来的才是直线的一般方程。Let the line s

2014-02-25 14:49:37 6575

转载旋转矩阵及旋转向量相互转化 Rodrigues矩阵

处理三维旋转问题时，通常采用旋转矩阵的方式来描述。一个向量乘以旋转矩阵等价于向量以某种方式进行旋转。除了采用旋转矩阵描述外，还可以用旋转向量来描述旋转，旋转向量的长度（模）表示绕轴逆时针旋转的角度（弧度）。旋转向量与旋转矩阵可以通过罗德里格斯（Rodrigues）变换进行转换。算法过程如下：式中，norm为求向量的模。反变换也可以很容易的通过如下公式实现：

2014-02-11 13:33:57 18565

原创 matlab 复制一个行向量为多个相同的行向量，用于求向量和矩阵距离

x=[1 2 3 ];out=repmat(x,[3,1]);out=[1 2 3;1 2 3; 1 2 3];

2013-11-22 14:58:44 15338 1

转载 matlab文件关联

cwd=pwd; cd([matlabroot '\toolbox\matlab\winfun\private']);% fileassoc('add','.m') ;fileassoc('add', {'.m','.mat','.mdl','.fig','.p','.mlprj','.mexw64'})cd(cwd); disp('Changed Windows file a

2013-11-05 09:59:57 1378 1

在利用cvFindContours查找轮廓时，时常要统计轮廓的个数，cvFindContours(src,storage,&count,sizeof(CvContour),CV_RETR_EXTERNAL,CV_CHAIN_APPROX_SIMPLE,cvPoint(0,0));刚开始，我是利用count->total统计个数的，但是结果却不对，其实这个结果返回的个数就是轮廓的个数。及int nu

2013-08-19 09:35:02 2813 3

原创计算机视觉共享代码备份

转载至http://www.cvchina.info/codes/1，通用库opencv http://www.opencv.org.cn/index.php/%E9%A6%96%E9%A1%B52，局部不便特征vlfeat http://www.vlfeat.org/index.htmlferns http://cvlab.epfl.ch/software/ferns/ind

2013-07-10 14:32:11 1150

原创 libsvm 之svmpredict 输出概率用法

［p1,p2,p3］=svmpredict(label_test,testdata,svmstruct,'-b 1');这个想必大家都用过，p1输出预估的类别，p2准确率，p3不同类的分类概率，但是这里面却有个小陷阱，之前我的理解，p3的一行代表不同类的概率，且是按照顺序排列的，既列的索引对应类别的索引，这次小师弟问我，才注意到并不是按照顺序的，而是在svmtrain训练的模型中既model.La

2013-06-06 15:37:54 18294 4

原创 linux 环境下java调用matlab小结

由于工作需要，需在linux环境下用java调用matlab，此中遇到诸多问题，现记录下来，以备不日之需1，安装虚拟机，我用的是vmware workstation 9.10版本，linux版本用到的是centOS6.0 64位系统，安装开发者界面，这样可以免除安装eclipse的繁琐，但是也有问题，这个后面介绍。2，安装matlab，此版本我用的是matlab2011b版本，这一切相信大

2013-03-02 14:27:57 2402 1

原创 labview访问数据库（labview sql tookit+sql2000）

链接方式无非两种常见的方式，odbc和udl文件，写下此文章仅记录我工作中遇到的问题，链接方式网上已经比较多了，不再熬述，在链接数据时会提示对象名无效，此时应将数据库的所有者跟表名中间加“.”链接起来就可以了。

2013-01-22 10:44:06 2250 1

原创吐槽一下opencv&xcode

近些天在改写代码，将原来matlab的代码改写成c的，这就用到opencv了，有个mac，就行在自己电脑上搭环境，需要安装xcode4.5，opencv，刚开始用cmake安装的opencv2.4.2版本的，用macport安装的，在这个网站上http://opencv.willowgarage.com/wiki/Mac_OS_X_OpenCV_Port用第一中方法安装，可后来一直出错（没有卸载，

2012-09-28 15:38:47 2034

原创 p 分位法实现图像二值化 matlab

最近在做字符识别，由于对于粘连情况识别无效，所以一直在考虑一种能够解决字符粘连的方法，而最为关键的是二值化过程，这种方法不能很好的分割字符，但是也是一种思路，写下来以备以后不时之需。代码很少，clc;clear;close allga=imread('1.bmp');if ndims(ga)==3 ga=rgb2gray(ga);endprc=prctile(ga,

2012-07-02 16:58:39 4809

原创图像主方向

参考文献：利用图像边缘信息估算图像纹理主方向matlab实现：function [mdex] = MainDirec(a)% function: the main direction of the image% description: the direction edge is [0,180],and the direction are divided%

2011-12-20 16:36:22 3335 1

原创医学图像数据库

DICOM sample image sets：http://pubimage.hcuge.ch:8080/ 上面不乏有ct、ads等等，(*^__^*) 嘻嘻……看有没有需要的吧@，

2010-12-13 11:52:00 1287

原创 matlab 生成EXE文件命令

已经有m文件和fig文件时，生成exe可执行文件方法： 在matlab的command窗口中输入 mcc -B sgl gui.m (gui是自己写的m文件) 2.将上步生成的exe和CTF文件以及MCRinstaller.exe这三个文件放到同一文件夹下一起拷到待运行的机器 (MCRinstaller.exe这个文件在<matlab path >/MATLAB/R2006a/toolbox/compiler/deploy/win32下可以找到此时仍需matlab所必需的动

2010-12-03 12:04:00 1630

原创解决“两个CU文件冲突”expected an expression

最近在学习cuda编程，刚刚起步，未免有很多困难，也想把自己的每一步的历程记录下来，将每一个所犯的错误记下来，以便将来不会犯同一个错误。这个程式是做的“GPU-高性能运算之CUDA”一本书的第二个例子，matrixassign，编译出现了如下错误 1>e:/CUDA/Practice/MatrixAssign/MatrixAssign/matrix_assign_kernel.cu(17): error: expected an expression 1>e:/CUDA/Pr

2010-12-02 10:49:00 5994 2

转载使用MATLAB进行USB摄像头的编程

1、查询USB2.0Camera 的具体参数(imaqhwinfo) 2、创建视频输入对象(videoinput) 3、图像预览和显示(preview、stoppreview、closepreview和image) 4、获取视频图像(getsnapshot) 5、图像获取设备的获取和设置(get和set) 6、关闭视频对象(delete) 在正式讲解之前，我想说明下几个个在图像获取工具箱中的术语： 图像获取设备：

2010-12-01 11:55:00 10445

原创 CUDA3.2+VS2008+CUDA_VS_WIZARD配置详细说明及第一个入门案例

摸索了好多天，终于完成了自己cuda环境的配置，可以编写简单的程序，但是还没有办法编译sdk中的例子，应该还是环境配置有问题吧。还待进一步努力，现将我的配置说明分享给大家，也是我从网路上寻得，希望大家少走弯路吧！在次感谢“http://1662170622.blogbus.com/logs/81511439.html”，因为我是按照他的配置来写的，一次成功。多谢！我的系统是win7 32位系统的。VS是2008版本的，CUDA_SDK和CUDAD_TOOLKIT都是3.2版本的。第一步：下载安装所

2010-12-01 10:13:00 1629 3

转载 debug release 的区别

关于Debug和Release之本质区别的讨论 一、Debug 和 Release 编译方式的本质区别 Debug 通常称为调试版本，它包含调试信息，并且不作任何优化，便于程序员调试程序。Release 称为发布版本，它往往是进行了各种优化，使得程序在代码大小和运行速度上都是最优的，以便用户很好地使用。 Debug 和 Release 的真正秘密，在于一组编译选项。下面列出了分别针对二者的选项（当然除此之外还有其他一些，如/Fd /Fo，但区别并不

2010-11-29 14:13:00 648

转载 Gray-level Co-occurrence Matrix 灰度共生矩阵

灰度共生矩阵就是一种通过研究灰度的空间相关特性来描述纹理的常用方法。 灰度共生矩阵是对图像上保持某距离的两象素分别具有某灰度的状况进行统计得到的。 取图像(N×N)中任意一点（x，y）及偏离它的另一点（x+a，y+b），设该点对的灰度值为（g1，g2）。令点（x，y）在整个画面上移动，则会得到各种（g1，g2）。值，设灰度值的级数为，则（g1，g2）。的组合共有 k2种。对于整个画面，统计出每一种（g1，g2）值出现的次数，然后排列成一个方阵，在用（g1，g2）出现的总次数将

2010-11-25 14:35:00 2712

转载图像的灰度共生矩阵GLCM(matlab 函数帮助)

Gray-level co-occurrence matrix from an image 图像的灰度共生矩阵 灰度共生矩阵是像素距离和角度的矩阵函数，它通过计算图像中一定距离和一定方向的两点灰度之间的相关性，来反映图像在方向、间隔、变化幅度及快慢上的综合信息。 使用方法： glcm = graycomatrix(I) glcms = graycomatrix(I,param1,val1,param2,val2,...) [glcms,

2010-11-25 14:34:00 18708 1

转载灰度共生矩阵（GLCM）

共生矩阵用两个位置的象素的联合概率密度来定义，它不仅反映亮度的分布特性，也反映具有同样亮度或接近亮度的象素之间的位置分布特性，是有关图象亮度变化的二阶统计特征。它是定义一组纹理特征的基础。 一幅图象的灰度共生矩阵能反映出图象灰度关于方向、相邻间隔、变化幅度的综合信息，它是分析图象的局部模式和它们排列规则的基础。 设f(x,y)为一幅二维数字图象，其大小为M×N，灰度级别为Ng,则满足一定空间关系的灰度共生矩阵为 P(i,j)=#｛(x1,y1),(x2,y2)∈M×

2010-11-25 14:31:00 12615 3

linxue博客

原创关于直线镜像矩阵生成