2017年11月_billbliss

原创 RANSAC LMedS 详细分析

RANSAC的函数接口参照opencv来说主要需要3-4个参数（第四个不是必须的）1. 误差阈值ransacThreshold：区分inlier和outliner的依据2. 置信度confidence：设置之后代表RANSAC采样n次过程中会出现（至少一次）采样点数据集中的点都为内点的概率这个值设置的太大，会增加采样次数。太小，会使结果不太理想。3. 最大采样迭代次数maxIte

2017-11-21 14:38:12 8867 2

转载线性参数估计方法之比较 LS、WLS、Ransac LS、LMedS

四种参数估计算法的简单描述： LS既是最小方差，以此作为目标函数求解参数估计值的方法称为最小方差估计，SVD分解是目前解决这一问题的最有效手段。 WLS既是带权重的最小方差，其思想是将每个输入的样本赋予权值，初始时每个样本的权值相等，然后使用所有带权重的样本估计模型的参数，得到参数后，计算每个样本与模型的偏差，再根据偏差决定样本的新权重，偏差越大则权重越小，然后重复模型参数估计

2017-11-21 10:37:31 6393 1

转载矩阵分解 SVD分解

1.前言一般提到特征值分解(eigenvalue decomposition)或者奇异值分解(singular value decomposition)，大多数同学脑海里的第一反应就是一大堆矩阵以及数学计算方法。确实，学校学习阶段，不管是学线性代数或者矩阵分析，对于这部分内容，或者说绝大部分内容，老师一上来都是吧啦吧啦给你一堆定理推论或者公理，然后就是哗啦哗啦一堆公式出来，告诉你怎么计算。

2017-11-20 09:46:20 28654 3

转载张氏标定相机标定 - 1

这里先分析一下张氏标定的闭式解（解析解）1. 预备知识=> 1.1. 从像素坐标系(u,v) 到世界坐标系(Xw,Yw,Yw)公式如下：=> 1.2. 符号规定( Notation )为了和张正友教授的论文相统一，现在把公式符号统一一下。=> 1.3. 推导中的数学基础第一点 : 旋转向量 R 为正交

2017-11-17 14:43:22 1132 2

转载矩阵分解 LDL^T分解

分解实际问题中，当求解方程组的系数矩阵是对称矩阵时，则用下面介绍的分解法可以简化程序设计并减少计算量。从定理可知，当矩阵A的各阶顺序主子式不为零时，A有唯一的Doolittle分解A= LU。矩阵U的对角线元素Uii 不等于0，将矩阵U的每行依次提出，下面将U分解为定理：若对称矩阵A的各阶顺序主子式不为零时，则A可以唯

2017-11-17 14:14:17 26148 1

转载矩阵分解 Cholesky分解

Cholesky分解法又叫平方根法，是求解对称正定线性方程组最常用的方法之一。对于一般矩阵，为了消除LU分解的局限性和误差的过分积累，采用了选主元的方法，但对于对称正定矩阵而言，选主元是不必要的。定理：若对称正定，则存在一个对角元为正数的下三角矩阵，使得成立。假设现在要求解线性方程组，其中为对称正定矩阵，那么可通过下面步骤求解（1）求的

2017-11-17 11:36:08 35707 1

转载矩阵分解三角分解(LU分解)

三角分解(LU分解)在线性代数中， LU分解(LU Decomposition)是矩阵分解的一种，可以将一个矩阵分解为一个单位下三角矩阵和一个上三角矩阵的乘积（有时是它们和一个置换矩阵的乘积）。LU分解主要应用在数值分析中，用来解线性方程、求反矩阵或计算行列式。本质上，LU分解是高斯消元的一种表达方式。首先，对矩阵A通过初等行变换将其变为一个上三角矩阵。对于学习过线性代数的同学来说

2017-11-17 11:32:17 94271 6

转载整理：一些主流的编程竞赛网站对比

今天把各个主流的编程网站仔细对比了一下，下面从各个角度对比一下这些编程竞赛网站。分类和排名网站地址类别PRATRQF主办方备注TopCoder综合633,5864167

2017-11-15 14:35:16 5504

原创 win下使用shell使用文件夹中的文件生成文件列表

使用windows中shell命令批处理实现：使用文件夹中的文件生成文件列表举例：将opencv库文件中debug lib生成一个列表，release lib生成一个列表DIR *320d.lib /B >dbglist.txtDIR *320.lib /B >rlslist.txt

2017-11-13 11:10:37 582

转载 SSE设计简介2

（1）MMX指令集首先要提到MMX指令集，MMX指令集是在SSE之前的，后来的SSE指令集覆盖了MMX指令集的内容，现在的大多数CPU也都支持SSE指令集了，SSE指令集之后还有SSE2、SSE3、SSE4等，最新的Intel处理器支持AVX指令集。（2）SIMDsingle instruction multiple data，单指令流多数据流，也就是说一次运算指令可以执行多个数据流，

2017-11-09 19:58:06 327

原创 SSE加速实战之二阶范数计算

这里使用了纯C，SSE c++加速版本，SSE 汇编加速版本测试求二维向量二阶范数的计算：平方和求平方根结果使用纯C反而更快环境WIN10 + VS2013结果是arrayCalcCPP 0.040msarrayCalcSSE 0.212msarrayCalcSSEASM 0.225msvoid arrayCalcCPP(F32 *p1, F32 *p2,

2017-11-09 19:13:25 1224

原创 SSE设计简介

SSE程序设计详细介绍包含的头文件：所有的SSE指令函数和__m128数据类型都在xmmintrin.h文件中定义：#include 因为程序中用到的SSE处理器指令是由编译器决定，所以它并没有相关的.lib库文件。数据分组（Data Alignment）由SSE指令处理的每一个浮点数数组必须把其中需要处理的数每16个字节（128位二进制）分为一组。一个静态数组（sta

2017-11-09 19:06:24 656

转载 neon指令速查

以下是neon的指令速查包含部分的neon和SSE的对应关系#ifndef __ARM_NEON__#error You must enable NEON instructions (e.g. -mfloat-abi=softfp -mfpu=neon) to use arm_neon.h#endif/*(1)、正常指令：生成大小相同且类型通常与操作数向量相同的结果向量；

2017-11-09 11:54:00 8048

原创并行加速实战双边滤波器

之前分析了二维中值滤波器的并行加速由于二维中值滤波器是控制密集型的滤波器（排序操作），所以SSE加速不太明显这次选用了计算密集型的双边滤波器针对双边滤波器在5*5的滤波核下的运算速度做优化和分析以下会有主区域、全图、主循环、完整（初始化+主循环）的概念1. 由于双边滤波的滤波半径为2+1，所以不能忽略图像四周边界的区域了。所以，以下会对主区

2017-11-03 18:09:49 2338

billbliss的专栏