碑湾没有碑的博客

原创 C++ Primer Plus (第6版)习题之循环和关系表达式

习题5.1#include<iostream>using namespace std;int main(){ int num1, num2; int sum = 0; cout << "Please enter two integers:"; cin >> num1 >> num2; for (int i = num1; i <= num2; i++) sum += i; cout << "The sum of al

原创 Python编程学习之基础数据类型

本文介绍Python语言的整型、浮点数、复数、布尔型和字符串这五种基本的数据类型。

原创 C++基础之循环语句

在程序设计过程中，有时候需要多次去执行相同的过程，例如打印“我爱我的祖国”这句话100次。如果我们像下面这样编写程序：cout<<"我爱我的祖国\n"；cout<<"我爱我的祖国\n"；cout<<"我爱我的祖国\n"；.....cout<<"我爱我的祖国\n"；不断的编写100条相同的语句，这会使得程序变得冗长。在C++编程语言中，通过引入一些循环语句来解决该问题，这些语句主要包含三种：for循环、while循环、do whie循环。本博客简单

原创 C++基础之数据类型（三）

在日常生活中，常常需要将一些信息存储在一个单元内，比如运动员的姓名、性别、工资、身高、体重等信息，这样一组数据是无法通过数组来实现存储的。数组在存储数据时，有一个必不可少的条件是这些数据的类型必须相同，不能将double类型与int类型的数据存储在同一个数组中。为了能够解决该问题，C++引入了一种新的数据类型，即结构。结构是一种比数组更灵活的数据格式，因为同一个结构可以存储多种类型的数据，即可以将double类型、int类型、char类型等存储在一起。结构是用户定义的类型，结构声明定义了这种类型的数据属性

原创 C++基础之数据类型（二）

我们常常发现简单的数据类型已经不能满足编程过程中所需要的数据要求。因此，在C++语言中提供了一些复合类型来表示这些数据，本博客主要是针对于数组进行说明。1 数组数组是一种数据格式，能够存储多个同类型的值。数组的创建可以通过声明语句来实现，数组声明时应当注意下列三点：存储在每个元素中的值的类型；数组名；数组中的元素数。声明数组的通用格式为：类型名 数组名[数组大小arraySize]## exampleint cards[10]float loans[50]char name[2

原创 C++ Primer Plus (第6版)习题之处理数据

习题3.1#include<iostream>using namespace std;int main(){ const int FeetTransInches = 12; int height_inches, height_feet, remain_inches; cout << "Please enter height(inches):__\b\b "; cin >> height_inches; height_feet = height_inch

原创 C++ Primer Plus (第6版)习题之开始学习C++

习题2.2#include<iostream>using namespace std;int main(){ const int Long_of_Yard = 220; int LongDistance; cout << "Please enter a distance (long): "; cin >> LongDistance; cout << LongDistance << " long is equal to " &lt

原创 C++基础之数据类型(一)

面向对象编程的本质是设计并扩展自己的数据类型，设计自己的数据类型就是让类型与数据匹配。在C++编写语言中，其数据类型大致包含了两种类型，其详细情况见下图：本博客主要对C++的基本数据类型进行相应的说明。1 整型整型是没有小数部分的数字，比如3, 20，-45等。在C++中整型大致包含了short int、int、long int等，不同的整型使用的内存量不同，内存量越大，可以表示整数值的范围越大。另外，有的类型可表示正值和负值，被称为有符号类型，还有一种只能表示正值和0，并不能表示负值，该类型为无符

原创 python中的条件分支与循环结构

1. 条件分支条件分支也可以成为选择结构，即根据某个条件决定执行不同部分的语句或语句块。在条件分支结构中，可以划分为单分支结构、二分支结构和多分支结构。1.1 单分支语句单分支语句可用if单分支语句来实现，其一般格式为：if 表达式： 语句块语句在执行过程中，若表达式的值为True，则执行语句中的语句块，否则，直接执行if语句的后续语句。if单分支语句流程图见下图1.例1：判断从键盘输入整数的奇偶性并输出结果。n = int(input("请输入一个整数："))flag = str

原创 图像处理一些大牛团队的网站

1、徐宗本研究方向： 智能信息处理、机器学习、数据建模基础理论研究。个人主页：http://gr.xjtu.edu.cn/web/zbxucn2、陈柯研究方向：图像处理的数学方法个人主页：https://www.liverpool.ac.uk/~cmchenke/3、Laurent D. Cohen研究方向：偏微分方程和变分法在图像处理和计算机视觉中的应用个人主页：...

原创 有用的一些网站

1、统计之都https://cosx.org/2、漫谈机器学习算法http://database.51cto.com/art/201407/444810.htm3、nuclear normhttps://math.stackexchange.com/questions/701062/derivative-of-the-nuclear-norm-with-respect-to-i...

原创 python中的列表数组

在C/C++编程语言中，同种类型的数据会通过数组的形式存放在一起，并借助于下标方式对其进行索引。python语言中也存在相类似的概念来存放一组数据，该操作被称为列表。1 创建列表列表的创建与普通变量的创建一样，通过中括号将一堆数据括起来，并采用逗号对数据进行隔开。不同于C/C++语言中的数组，python的列表可以将不同数据类型存放在一起。例如：number = [1,2,3,4,5]mix = [1,'love', 3.56,[1,2,4]]empty = [] # 空列表的创建2 元

原创 python的数据类型与操作符

在数据运算过程中，常常会使用不同的数据类型（整数、小数、字符等）与运算操作，本文简单说明一下关于python语言中数据类型与操作符。1. 数据类型python中的数据类型大致包含了整型（int）、浮点数（float）、布尔类型（bool）。1.1 整型与浮点数python中的整型与长整型进行了无缝结合，它的长度只限于计算机的虚拟内存总数。整型与浮点数最大的区别大于是否存在小数点。1.2 布尔类型布尔类型包含了True和False，其分别来表示“真”和“假”。另外，布尔类型也可以当作整数来处理，即

原创 python中的random模块的使用

在程序设计过程中，我们可能需要引入一些随机数据，本文简单的说明一下python语言中关于random模块的使用。1、random.random()表示随机产生0到1之间的浮点数import randomprint(random.random())#-----#0.304439830825960162、random.randint(a,b)表示随机产生[a,b]之间的一个整数import randomprint(random.randint(2,5))#-----#

原创 扩散加权磁共振成像(DWI)

扩散加权磁共振成像(diffusion-weighted imaging, DWI)是20世纪90年代初中期发展起来的MRI新技术。DWI能够无创地测量生物组织内水分子的扩散速度。水分子扩散是一个三维过程，在理想环境下，水分子沿各个方向的扩散速度使相同的，称为各向同性扩散。然而，在人体组织内，水分子受到各种因素的影响（如脑白质纤维，由于髓鞘阻碍了水分子垂直于纤维方向的扩散，水分子表现为沿纤维束方向...

原创 矩阵的Kronecker积、Khatri-Rao积、Hadamard积

1、矩阵的Kronecker积Kronecker积也称为克罗内克积，是任意大小矩阵的运算，使用符号其表示为：若A为大小m*n的矩阵，B为大小p*q的矩阵，则A与B的克罗内克积是一个大小为mn*pq的矩阵，其表述为： 其具体形式为： ...

转载 课程笔记4--图像K空间理解

K空间的数据分布实际上是图像空间中数据的二维傅立叶变换结果。K空间中的数据点和图像空间中的数据点并不是一一对应的。一个K空间中的数据点对应了图像空间中所有数据点的一部分信息。事实上，K空间中的数据正是图像空间中的数据作二维傅立叶变换的结果（图1），也就是说，我们的“大脑图像”可以被看作是由一系列频率、相位、方向各异的二维正弦波叠加而成的，而K空间的数据正表示了图像的正弦波组成。因此，为了理解如...

原创 如何使用MATLAB去给图像添加噪声

图像去噪是数字图像处理中最为关键的环节，去噪效果的好坏也会影响到图像的后续处理，比如图像分割、特征提取、目标识别等。数字图像在产生和传输过程中都可能会受到噪声的破坏，一般有高斯噪声、脉冲噪声、泊松噪声等。在进行数字图像去噪实验时，我们经常使用MATLAB工具去实现，本文主要讲解的是如何使用MATLAB对上述噪声进行模拟仿真。1、高斯噪声高斯噪声，从字面含义上可以看出，其概率密度函数服从高斯...

转载 拉格朗日乘数法

阅读目录1. 拉格朗日乘数法的基本思想 2. 数学实例 3. 拉格朗日乘数法的基本形态 4. 拉格朗日乘数法与KKT条件　　拉格朗日乘数法（Lagrange Multiplier Method）之前听数学老师授课的时候就是一知半解，现在越发感觉拉格朗日乘数法应用的广泛性，所以特意抽时间学习了麻省理工学院的在线数学课程。新学到的知识一定要立刻记录下来，希望对各位博友有些许帮助。回到顶...

转载 机器学习中的数学(5)-强大的矩阵奇异值分解(SVD)及其应用

前言： 上一次写了关于PCA与LDA的文章，PCA的实现一般有两种，一种是用特征值分解去实现的，一种是用奇异值分解去实现的。在上篇文章中便是基于特征值分解的一种解释。特征值和奇异值在大部分人的印象中，往往是停留在纯粹的数学计算中。而且线性代数或者矩阵论里面，也很少讲任何跟特征值与奇异值有关的应用背景。奇异值分解是一个有着很明显的物理意义的一种方法，它可以将一个比较复杂的矩阵用更小更简单的...

转载 机器学习中导数最优化方法(基础篇)

1. 前言熟悉机器学习的童鞋都知道，优化方法是其中一个非常重要的话题，最常见的情形就是利用目标函数的导数通过多次迭代来求解无约束最优化问题。实现简单，coding 方便，是训练模型的必备利器之一。这篇博客主要总结一下使用导数的最优化方法的几个基本方法，梳理梳理相关的数学知识，本人也是一边写一边学，如有问题，欢迎指正，共同学习，一起进步。2. 几个数学概念1) 梯度（一阶导数）...

转载 漫谈：机器学习中距离和相似性度量方法

在机器学习和数据挖掘中，我们经常需要知道个体间差异的大小，进而评价个体的相似性和类别。最常见的是数据分析中的相关分析，数据挖掘中的分类和聚类算法，如 K 最近邻（KNN）和 K 均值（K-Means）等等。根据数据特性的不同，可以采用不同的度量方法。一般而言，定义一个距离函数 d(x,y), 需要满足下面几个准则：1) d(x,x) = 0 // 到自己的...

原创 MRI图像噪声分析

MRI图像噪声的最为主要的来源是检测对象和电路元件，产生的噪声主要为热噪声，有时候还存在生理学噪声。磁共振设备采集到的原始数据是K空间数据，其包含实部和虚部这两部分的信号，且两者之间的相位差为90°。实部和虚部分别带有均值为0，方差相同且独立的加性高斯白噪声。对K空间数据进行傅里叶变换，得到图像域数据。傅里叶变换不改变噪声形式，因此我们得到的图像域数据是复数信号，其表达式为： ...

转载 BrainWeb: Simulated Brain Database使用说明

阅读目录(Content)BrainWeb: Simulated Brain Database使用说明 一、文件命名规则 二、文件下载格式：.rawb 三、离散标签命名规则 回到顶部(go to top)BrainWeb: Simulated Brain Database使用说明作者：凯鲁嘎吉-博客园http://www.cnblogs.com/kai...

转载 散粒噪声 shot noise 泊松噪声

在raw image中，主要的噪声为两种，高斯噪声和散粒噪声，其中，高斯噪声是与光强没有关系的噪声，无论像素值是多少，噪声的平均水平（一般是0）不变。另一种是散粒噪声，因为其符合泊松分布，又称为泊松噪声，下图可见，泊松噪声随着光强增大，平均噪声也增大。什么是散粒噪声？散粒噪声=泊松噪声=shot noise=poisson noiseShot noise存在的根本原因是因为光是由离...

原创 CT剂量和降低剂量的方法

X射线穿过介质会发生相互作用，此时能量将从X射线光子转移到介质上。因此，X射线辐射的剂量描述的是在与介质相互作用之中单位质量介质内能量转移多少的度量。主要包含两种：吸收剂量（absorbed dose）和等效剂量（equivalent dose）。吸收剂量指单位质量介质吸收能量的多少，单位是Gray（Gy）。吸收剂量描述人体吸收了多少能量。为了能够量化不同电离辐射在相同的吸收剂量下对生物组织造成...

原创 CT中的噪声和伪影

CT图像中存在影响图像质量的噪声和伪影。1、CT图像中的噪声主要有三种来源。（1）量子噪声。量子噪声服从泊松分布，是一种实验观测中的读出噪声。当观测中数量有限的X射线光子数量少到能够引发探测器上数据读出中出现可观测到的统计涨落时，这种读出的统计涨落被称作量子噪声。这种噪声也是低剂量CT中面临的主要噪声 。量子噪声主要是由X射线流以及探测器上探测到的X射线光子数决定的，主要受制于两个因素：扫描方...

原创 如何使用MATLAB实现标记感兴趣区域并将其放大

在图像处理中，常常会遇到这样一个问题：如何将图像中某一小块使用红色边框标记，然后再将其放大放置于左下角？之前在网上查找了很多关于此类方法的解决方法，但是都失败了。最后自己编写了一个小的程序来实现改功能。其代码如下：Img = imread('Lena.png');Local_Zoom = ShowEnlargedPatch(Img, 241, 426, 70, 70, 2, 2);imwri...

原创 如何采用Matlab来实现图片保存时白色边框的消除？

在做图像处理时，常常会遇到保存所绘制的图像或曲线后，所保存的图片周围有一圈白色的边，如何将这些白色边框消除掉？ 本文阐述了如何使用matlab来完成该任务。1、采用imshow()函数显示图像在图像处理过程中，经常使用imshow()函数显示图像，并采用imwrite()函数或者saves()函数进行保存，但是这常常面临着一个问题，所保存的图片都存在白色边框。我们可以采用下列函数来消除：[M...

原创 如何使用Matlab绘制图形并标记图例

当我们在写cscd、EI或者SCI论文时，常常需要放置一些图片，通常情况下需要对图片进行一些处理，对图像横纵坐标、曲线进行标记。本文阐述如何使用Matlab来实现这些功能。1、如何使用Matlab绘制相关曲线图在Matlab中，常采用plot函数来绘制二维图像。其函数格式如下：plot(x,y) 如果x,y都属于向量，则它们必须具有相同的长度，plot函数将以x为横轴，绘制出y。如果x,y都...

原创 医学图像 成像机理

1、X射线影像1.1 X射线的产生 X射线也被称为伦琴射线，是1895年德国科学家伦琴发现的具有很高能量、肉眼看不见但能够穿透不同物质且使荧光物质发光的射线。X射线是真空管内高速行进的电子流轰击钨靶时产生的。其发生过程是X射线管灯丝供电、加热，在阴极附近产生自由电子，当向X射线管两极提供高压电时，阴极与阳极间的电势差陡增，电子以高速由阴极向阳极行进，轰击阳极钨靶而发生能量转换，其...

转载 MATLAB实现dicom转换为png格式

clcclear%从dicom文件中读取数据%首先获得指定文件夹下的所有文件的文件名directory = uigetdir('','选择指定文件夹：');dirs=dir(directory);%dirs结构体类型,不仅包括文件名，还包含文件其他信息。dircell=struct2cell(dirs)'; %类型转化，转化为元组类型filenames=dircell(:,1) ;%...

原创 变换域去噪技术

变换域去噪的基本思想是先将图像从空间域变换至某个变换域中，随后再对变换域中的稀疏进行处理，最后通过反变换回到原始空间域来达到去除噪声的目的。典型的方法有小波变换去噪和多尺度分析去噪。1.小波变换去噪 小波变换去噪是将图像变换到小波域，利用图像和噪声在小波变换下的不同特性，对小波稀疏进行处理，通过最大限度地保留图像信号的系数而减少噪声系数来除去噪声。常见的小波去噪方...

原创 小波阈值去燥

转载 低秩恢复算法（图像去噪）

转载自：http://blog.sina.com.cn/s/blog_bfb629e50102xhsw.html                                                                                   低秩矩阵恢复算法 概述近几年，低秩矩阵恢复（LRMR）广泛用于图像处理用途图像恢复，比如去噪、去模糊等。一...

原创 稀疏模型求解优化算法

       稀疏模型求解，是给定字典或固定基D，求解对图像或信号x的最逼近原始真实图像的稀疏系数α的过程。当前，关于稀疏求解问题的优化算法大致归纳为：针对范数最小化问题的贪婪算法，是直接对范数进行求解；针对范数最小化问题的最优化算法（其包括范数正则化算法，是范数求解法和范数求解法的折中），是把范数松弛到范数进行求解；统计优化算法，是从数学期望中估计图像的最稀疏解。范数近似求解算法     ...

原创 稀疏性度量

       对于特定的稀疏表示方法，其优劣一般采用表示系数的稀疏程度来度量。图像的稀疏性一般采用范数来定义，即不为零的系数数量越少，越稀疏。然而，在图像处理的实际应用中，通常受到噪声干扰，采用范数来度量稀疏性的效果会变差。在数学分析理论中，范数常被用来度量稀疏性，其能较好地度量非零系数的个数与系数误差。而且，当0&lt;p&lt;1时，与存在一定的关系，当p→0时，             ...

原创 稀疏表示理论问题优化模型

       稀疏表示是针对不确定方程组的求解问题而提出的。对不确定方程组Ax=b,A∈，且n&lt;m，在A为满秩的条件下，方程数目小于未知数的个数，方程组有无数组解。给定不同的正则约束函数，不确定性方程组的解也会有所不同。稀疏表示模型是将解的系数性作为不确定性方程组一种约束，使方程组有唯一解。       稀疏表示的本质思想是：引入到信号表示理论中，从基或字典中寻求最小原子的线性组合来表示...

原创 图像噪声的分类与模型

        噪声是干扰和妨碍人类认知和理解信息的重要因素，而图像噪声则是图像中干扰和妨碍人类认识和理解图像信息的重要因素。由于噪声本身具有不可预测性，可以将它当做一种随机误差（这种误差只有通过概率统计的方法来识别）。因此，图像噪声可以视为一种多维随机过程，可以选择随机过程的概率分布函数和概率密度函数来作为对图像噪声进行描述的方法。1 图像噪声的分类        图像在采集和传输的过程当中...

转载 keras小数据集扩充

@原文地址：http://blog.keras.io/building-powerful-image-classification-models-using-very-little-data.html参考译文地址：http://keras-cn.readthedocs.io/en/latest/blog/image_classification_using_very_little...