- 博客(196)
- 收藏
- 关注
原创 噪声
https://blog.csdn.net/matrix_space/article/details/105745560/https://blog.csdn.net/lpsl1882/article/details/51339742
2020-12-03 14:03:32
373
原创 clang - msvc
想体验下openmp后续的一键GPU优化新特性,可惜VS只支持openmp2.0。https://devblogs.microsoft.com/cppblog/clang-llvm-support-in-visual-studio/然后装了个clang体验,结果neon_2_sse.h这个文件各种编译通不过。error报各种未定义标识符:__fp16 __builtin_ia32_shufps __builtin_ia32_pshufd等等网上也找不到相关解决方案,看来只能等后续V...
2020-11-10 17:00:22
581
原创 interpolation method
Cubic Key:cubic B-spline:cubic o-Moms:Thévenaz P, Blu T, Unser M. Image interpolation and resampling[J]. Handbook of medical imaging, processing and analysis, 2000, 1(1): 393-420R. Keys (1981). "Cubic convolution interpolation for digit...
2020-10-15 16:48:11
871
原创 Bicubic interpolation
网上的公式基本指的是Cubic convolution interpolation。看了下这篇文章,主要就是提出了一个 双三次卷积插值算法,这个算法精度介于线性插值与三次样条插值之间,并且是三阶近似。——————————————————————————————————————————————————————插值kernel公式:1. 插值使用特定区间分段三次多项式,并且插值的权重kernel是对称的 (为什么?),那么形式如下,8个未知数待定2.插值kernel满足整点...
2020-10-13 17:21:58
908
2
原创 bayer AHD插值
最近做RAW,遇到一些issue涉及到raw的伪彩,为了有希望解决这个问题,去看了下AHD paper,希望能学到一些技巧有助于帮助伪彩的解决。这里做个记录。AHD插值step1. G的插值: paper里用最优化过程导出了一个filter。filter 长度为5,对于R0,G1,R2,G3,R4,运用该filter可得到G2的值。这个定点化后的filter的理解(性质): 若进行变形可得到 a. G2 = weight(Gi-Ri)+R2,也就是说R2加上一...
2020-09-28 09:17:58
1656
原创 omp
int layer_for = 2;int each_task_size_for = 100000 / 20;int num_threads = omp_get_max_threads();// int num_threads_limit = omp_get_thread_limit();//编不过?printf("多线程:%d\t%d\n", num_threads);omp_set_num_threads(num_threads);#pragma omp parallel num_th..
2020-08-05 09:56:16
265
原创 超高速指数模糊算法
网上有一篇非常快的模糊算法实现:SSE图像算法优化系列五:超高速指数模糊算法的实现和优化(10000*10000在100ms左右实现)。名字叫 Exponential blur(指数模糊)。在搜材料时搜到只搜到了代码:blur.cpp两个链接的核心代码就是 *zR += (alpha * ((R << zprec) - *zR)) >> aprec;事实上可以看出代码实现的就是Deriche filter里提到的公式:a first order I...
2020-07-28 09:45:12
530
原创 问题:矩阵分解角度
一个链接:里面展示了一堆问题以及矩阵分解角度来描述看待该问题:The Advanced Matrix Factorization Jungle
2020-07-23 15:10:39
173
原创 harris角点,亚像素 以及 hdrplus的亚像素对齐
最近在弄对齐,对齐里有个亚像素对齐,然后opencv里又有个角点亚像素refine,然后就到了角点...,这一系列关联下去,就去学了个角点,然后这些基本推导都非常类似... 看上去只要是扯到图像结构啥的,都会涉及到梯度啊,局部极小,然后taylor,hessian,特征值一系列都会有所涉及,所以原理推导基本都很类似。 1.harris角点认为局部移动导致变化强烈的点就是角点,然后公式里talor展开合并得到一个二次型,矩阵就是类hessian。二次型矩阵是实对称矩阵A,可以正交分解成...
2020-07-17 09:05:34
645
原创 网上一个网站(線代啟示錄)
最近在做对齐相关的工作,查资料时偶然查到了一个网站,感觉帮助很大,做个记录線代啟示錄1、一元多次多项式求根直接转化为 对伴随矩阵求特征值2、里面有许多专题:傅里叶变换,奇异值分解,奇异值与特征值联系,KKT条件,马尔科夫等等,都值得一读pageRank:PageRank傅里叶变换:快速傅立葉轉換 傅里叶变换应该还有取模形式(精度),但是任意模数的实现资料有待查找KKT条件:KKT条件母函数:母函数一个有趣的问题,矩阵乘链的计算顺序,看上去就是动...
2020-07-13 15:51:59
1084
原创 知乎上的一个数值计算专栏
数值计算实验室查资料过程中偶然发现了知乎上看到的一篇文章,并由此看到了相应的一个专栏,看了好几篇,太赞了!正是我一直迫切需要的,数值知识梳理总结理解都有,还有推导过程以及点评等,太牛叉了以后很多知识都可以直接在此看,若需要更详细再进一步查,可省很多时间。...
2020-07-03 15:25:41
355
原创 tonemap SR
一种快速简便优秀的全局曲线调整与局部信息想结合的非线性彩色增强算法(多图深度分析和探索)SSE图像算法优化系列二十:一种快速简单而又有效的低照度图像恢复算法。低对比度图像增强及图像边缘检测方法SSE图像算法优化系列二十八:深度优化局部拉普拉斯金字塔滤波器。论文笔记:Handheld Multi-Frame Super-ResolutionHMSR - 根据图像梯度方向生成滤波kernel权重Kernel Regression for Image Processing and R.
2020-07-03 15:02:52
318
原创 性能优化相关以及一些困惑
最近稍微看了些CPU性能相关,大致浅显地了解了下“缓存”“cache line”“false sharing”“缓存分块算法”对性能的影响。在理论了解后,打算实际应用时产生了困惑:1.缓存大小已知后,实际编码时怎样简单快速地感知当前代码在CPU上的缓存表现?如果使用缓存分块算法,怎样确定计算具体分块参数?这些疑惑网上没找到解答,网上多数都停留于理论介绍上,实际如何应用的材料还没看到,快速操作推荐方法没怎么看到。 看来只能日后自己在代码中尝试做下实验进行实际感受了。2.既然性能还与缓存相关,那么给定.
2020-06-30 16:03:35
151
原创 poisson快速算法
泊松融合是个经典问题,网上已经有很多材料介绍了,自己以前也有做过笔记:poisson方程组矩阵前段时间,项目中用了泊松融合,为了能在产品使用,对性能有一定要求,通用的解法基本是不满足要求的。后来网上看到了一个关于poisson的优化材料,code实现了下,性能优化非常明显,可见数学优化还是非常重要。材料链接:Solving the Discrete Poisson Equation using Jacobi, SOR, Conjugate Gradients, and the FFT各种算法复杂
2020-06-29 15:05:38
1775
原创 neon笔记
一个sse2neon:https://github.com/DLTcollab/sse2neonneon指令的使用介绍与范例:https://blog.csdn.net/EmSoftEn/article/details/51718763Neon Intrinsics各函数介绍
2020-06-24 14:20:57
232
原创 Hessian笔记
Hessian笔记最近做了一些DIS光流的工作。过程中文献里提到了一个Hessian矩阵,但其表达式又与二阶导不同,所以这里是个疑惑点,做个记录。如链接:http://image.sciencenet.cn/olddata/kexue.com.cn/upload/blog/file/2010/9/2010929122517964628.pdf同时还有: “根据费舍尔信息矩阵的定义,梯度的外积矩阵是 Hessian 矩阵的一个渐近无偏估计:”而在图像中,梯度的外积矩阵会遇到的比较多,需要注意两者是有区
2020-06-19 10:06:27
264
原创 Recursive filter笔记记录
Recursive gaussian filter是gaussian filter的一种近似实现,可以参看这篇survey: A Survey of Gaussian Convolution AlgorithmsPascal Getreuer其中之一是Vliet–Young–Verbeek :1.Recursive implementati...
2019-11-07 14:56:17
1404
原创 ARM和x86区别
https://blog.csdn.net/qq_34385566/article/details/80182441ARM X86 Neon SSE指令集各不同不兼容
2019-07-01 15:32:02
983
原创 dp四边形优化
四边形优化资料: https://blog.csdn.net/noiau/article/details/72514812 题目链接 : postoffice已过: https://cn.vjudge.net/problem/POJ-1160 #include<iostream>#define maxlen 2001+90#define maxva...
2018-12-18 23:01:49
262
转载 yuv rgb理解笔记
(转)视频与图像RGB/YUV格式详解但是对于视频捕获和编解码等应用来讲,这样的表示方式数据量太大了。需要想办法在不太影响感觉的情况下,对原始数据的表示方法进行更改,减少数据量。 无论中间处理过程怎样,最终都是为了展示给人观看,这样的更改,也是从人眼睛的特性出发,和发明RGB三原色表示方法的出发点是一样的。 于是我们使用Y,Cb,Cr模型来表示颜色。Ia...
2018-11-13 14:23:50
333
原创 傅里叶变换
很多东西学了后就忘光了,淡淡的忧伤。本科数学分析学了后,教材保留带到研究生阶段。研究生毕业竟然扔了,现在到用到的时候,感觉亏死了。还好还是留有一些印象,在网上找了些资料又回顾了一些。 但还是觉得好亏,怎么能扔了呢!!!傅里叶级数展开:点击打开链接傅里叶级数转化为复数形式:点击打开链接算法导论里是讲傅里叶变换的应用在多项式相乘上。里面嵌套了DFT的如何快速求...
2017-11-13 17:19:49
437
原创 gamma矫正 记录
http://blog.csdn.net/candycat1992/article/details/46228771读完后,也就是说现在的rgb图像0-255只是一种采样计数。原本可以1代表0.1亮度,2代表0.2。但是由于人眼对于亮度暗的地方更加敏感,所以可以在暗处采样点更多,两处少一点采样,这样就能有效提高图像的显示效果。这样一来,采样点与真实亮度之间的关系就不是线性的了,而...
2017-07-26 14:12:27
382
转载 dos文件(夹)复制命令:copy和xcopy
http://www.cnblogs.com/KevinJasmine/p/4159234.html
2017-03-29 14:30:49
779
原创 pro文件
QT += coreQT -= guiTARGET = opencv300TestCONFIG += consoleCONFIG -= app_bundleTEMPLATE = appDESTDIR = binSOURCES += main.cpp \# ImageFactory.cpp#OPENCV_PATH= F:/8rfuz/
2017-03-27 13:17:40
652
转载 【转】C Runtime Library 的来历, API, STL, MFC, ATL关系
http://www.cppblog.com/dawnbreak/archive/2009/08/11/92906.html
2017-03-24 16:19:18
344
转载 PASCAL VOC数据集分析
http://blog.csdn.net/zhangjunbob/article/details/52769381
2017-03-22 10:04:46
487
原创 faster-rcnn相关
http://blog.csdn.net/sinat_30071459/article/details/53202977
2017-03-21 10:06:47
285
原创 机器人SLAM算法漫谈
http://mp.weixin.qq.com/s?__biz=MzIzNDM2OTMzOQ==&mid=2247485029&idx=1&sn=bfa57a1113734a9ab639477e8f70f91a&chksm=e8f63e31df81b727f5a40627e277e299b0419c71b33e805f1a80363d358eefca0531c92362b5&mpshare=1&s
2017-03-14 17:01:08
534
转载 记 整数划分(integer partition)
好几次了,怕以后找不到了,还是记录一下网址吧。http://www.cnblogs.com/xubenben/p/3664959.htmlhttp://duoerbasilu.iteye.com/blog/1606001
2017-02-20 16:43:37
680
原创 Sort函数
关于sort函数的几种排序方法:http://www.cppblog.com/mzty/archive/2005/12/15/1770.html
2017-01-09 01:25:28
278
原创 51nod 1608 管道安装
51nod 链接 点击打开链接codeforces 115C 链接 点击打开链接本题求放置水管的方式的总数。管道就4个,而且管道的出口都是一边竖直一边水平。分析后可以发现,水平方向与竖直方向可以分离开来进行考虑。对于特定的一种放置管道的方式,可以按如下方式进行转化为0,1矩阵:只考虑竖直: 令向上开口为0,向下为1,则化为0,1的一个矩阵。 上下方向0,1
2016-12-30 14:29:35
390
原创 BSP代码 日志记录
主要记录一些觉得可能存在的问题。最近以为代码有个潜在的bug,只是因为特例的特殊,所以出现概率较小。 所以还是记录一下比较好,记忆太差。存在的bug是这样的:当建成的BSP树对某个polygon A进行空间划分切割时,如果出现共面情形,并且该 划分超平面上同时有正负2种polygon时,并且A同时与这2种polygon有重叠时。这时候该polygonA哪些
2016-12-17 16:15:43
673
原创 BSP理解笔记
最近看BSP的早期原始论文,里面好几部分都花了好多心思才转过弯来,也没个图示差点就没反应过来。一堆超平面 可以把3d空间划分成一堆的凸体,这些凸体的 包围面,其实就是sub 子超平面。 而上面就是在 已知超平面的法向量的条件下,怎么求得这些sub子超平面的法向量。 具体规则就是 根据子超平面 所在的凸体空间“内外”属性以及这些凸体与超平面的位置关系 来共同
2016-12-15 14:56:59
431
原创 51nod 1588 幸运树
51nod 1588 幸运树codeforces lucky tree这一题第一眼看去就认为是树形dp,然后推出公式,时间复杂度O(n)*check lucky的代价,最多9。并且可以对公式化简合并,得到简洁的形式。然后就没继续想下去了,直接套模板了。然后才发现解题报告里有另一种方法,考虑了一下,该方法才更接近问题的本质,将问题转化为遍历联通块了,即可得到最终答
2016-12-14 11:25:30
663
原创 51nod 异或凑数
此题是一道有关 异或的秩的题目,已经有解题报告了,不过感觉解题报告里没看到 严格的证明,所以记录下自己的思路。这道题要算给定区间内能否存在数异或成给定的数K。首先数目达到500000,所以要考虑线性算法,当然nlogn也是可以的,不过一般都是希望O(n)。一开始没什么思路,不过后来发现给定区间里的选择的数的数目不定,并且并不要求连续,还觉得比较困难。想到异或有性质就是 秩
2016-12-08 11:50:07
457
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人