织田信长-CSDN博客

原创【OpenCV】基于kmeans的细胞检测方法

问题是这样的，有一幅经过二值化处理之后的图像，我们希望统计其中细胞的个数，和不同粘连情况的细胞个数，比如，下图中有1个细胞组成连通区域的，也有2个细胞组成连通区域的，也有更多个细胞组成连通区域的，我们希望分别统计不同的情况。我想出的一种可行的方法是这样的：通过图像形态学的处理erode，将一些邻接的细胞分离开来，并减少单个像素的噪声干扰计算其中的连通域计算每一个连通域的面积

2013-12-15 21:04:36 1445

原创【OpenCV】基于Python的图像高斯平滑和椒盐噪声处理

最近要做一个Project，是使用TV来对添加了Gaussian和Salt&Pepper噪声的图像进行恢复，前期的任务是生成噪声污染的图像。噪声图像的生成采用对图像进行高斯平滑，之后在随机的对图像添加椒盐噪声。里面用到的几个opencv的方法：random.random_integers:产生范围内的随机整数cv2.GaussianBlur:对图像进行高斯滤波源代码：

2013-11-05 22:23:49 1707 1

原创【数字图像处理】使用kmeans算法对TrueColor图片进行优化

实验的主要内容是将truecolor的图片通过一个优化算法得到其256色的最优表示。本实验采用kmean做算法对像素的色彩进行聚类的计算，分类得到一个色彩数为256的CodeBook，和一个包含有CodeBook索引的矩阵，原图片的颜色值都从CodeBook中获取，colorMap矩阵保存其中像素颜色的索引值，这样通过索引恢复出一个近似的图片文件，达到比传统的GIF图像更好的效果。

2013-10-31 21:23:27 1164

原创【StatLearn】统计学习中knn算法的实验（2）

接着统计学习中knn算法实验（1）的内容Problem：Explore the data before classification using summary statistics or visualizationPre-process the data (such as denoising, normalization, feature selection, …) Try ot

2013-10-31 21:05:26 1098

原创【StatLearn】统计学习中knn算法的实验（1）

Problem：Develop a k-NN classifier with Euclidean distance and simple votingPerform 5-fold cross validation, find out which k performs the best (in terms of accuracy)Use PCA to reduce the dimensi

2013-10-31 20:38:04 1538

原创【Python】计算GPA

最近开始我的研究生生涯了，周围都在谈论GPA这种高深的东西，像我这种工科学校背景的孩子一向是个喜欢在键盘上耕耘的人，对考高分没有多大兴趣（其实是没有那个本事）。不过我也想知道我自己本科的GPA是多少来着，网上的各种工具都不是很理想的，决定自己写吼。from numpy import *from scipy import *import csvi def USTC_GPA(r):

2013-09-08 13:07:39 1430

原创【Python】Coding the Matrix:Week 6 Secret Sharing Lab

这个Problem 2 有点困惑的，我通过程序获得了一组3对线性独立的向量，但是任然不过，过会再继续，但是代码的原理应该是可以的。choose_rand_vector()返回一个符合要求3个tuple的listchoose_rand_vector(randNum=4,vecLength=6):choose_pair()验证所有的vector组合是否满足6个向量相互线性独立的要求def

2013-08-27 13:38:14 927

原创【Python】Coding the Matrix:Week 5 Perspective Lab

这个Lab的内容光是说明就有7页之巨，我反复看了很久才看懂一点点，Lab主要完成的是从不同坐标系表示之间变换的方法。原始的图片，从Camera basis的表示转换成WhiteBoard basis的表示里面的Problem 3是难点，Problem 4我没有完成，因为还缺少之前的代码，暂时不写。注意Problem 3中的vector h不能通过print(h)来获得，

2013-08-26 18:06:20 1211 2

原创【Python】Coding the Matrix:Week 5: Dimension Homework 5

这一周的作业，刚压线写完。Problem3 没有写，不想证明了。从Problem 9 开始一直到最后难度都挺大的，我是在论坛上看过了别人的讨论才写出来的，挣扎了很久。Problem 9在给定的基上分解向量，里面调用了hw4的一些函数，通过solve函数获得矩阵方程的解Problem 10判断矩阵是不是可逆的，注意判断矩阵是不是square的Problem 11和Problem 12 都

2013-08-26 12:07:48 989

原创【Algorithm】Random Contraction Algorithm

写

2013-08-15 17:56:13 387

原创【Python】Coding the Matrix:Week 3 The Matrix

写Homework之前需要写Matrix的一些函数，也就是完成Matrix Class Homework。这篇文章也就先写到这里先。这个类主要实现矩阵的加法，转置，矩阵和向量的乘法，矩阵和矩阵的乘法。from vec import *def getitem(M, k): "Returns the value of entry k in M. The value of k sh

2013-08-10 12:59:52 790

原创【Python】Coding the Matrix: Week 3: The Matrix

Matrix Class Homework

2013-08-10 12:53:14 116

原创【Python】Coding the Matrix:Week 2 The Vector Space

这周的作业还比较少，Python的语言确实比较精简，所以如果代码写的优雅的话，代码量可是很少的哦！# version code 761# Please fill out this stencil and submit using the provided submission script.from vec import Vecfrom vec import *import iter

2013-08-08 11:43:12 966

原创【OpenMesh】Training：Getting Started with OpenMesh

内容出自236329 - Digital Geometry Processing, Spring2013 - Announcements上的OpenMesh Assignment感觉比较好，适合刚开始学习的同学。练习要求：这次的练习主要是写两个函数，calc_valences()和color_coding()，函数主要完成了计算一个Vertex的valence以及根据这个valence给图

2013-08-06 21:57:03 1789

原创【OpenGL】细分三角形

#include typedef GLfloat point2[2];point2 v[]={{-1.0,-0.58},{1.0,-0.58},{0.0,1.15}};int n;void mydisplay(){ glClear(GL_COLOR_BUFFER_BIT); glColor3f(0.5,0.5,0.4); glBegin(GL_LINE); glVertex2

2013-07-30 20:44:10 690

原创【OpenGL】细分三角形

#include typedef GLfloat point2[2];point2 v[]={{-1.0,-0.58},{1.0,-0.58},{0.0,1.15}};int n;void mydisplay(){ glClear(GL_COLOR_BUFFER_BIT); glColor3f(0.5,0.5,0.4); glBegin(GL_LINE); glVertex2

2013-07-30 20:41:37 716

翻译【OpenMesh】存储自定义属性

OpenMesh的专有OM格式允许存储和恢复自定义属性,以及标准的属性。我们必须使用指定的自定义属性，像下面的代码显示的一样。OpenMesh::VPropHandleT vprop_float;mesh.add_property(vprop_float, "vprop_float");我们给顶点注册一个浮点属性，名为"vprop_float"。属性的名字需要符合一些规定，如果我们想要长

2013-07-22 22:16:40 1233

翻译【OpenMesh】使用IO::Options

这个例子展示了：如何控制Mesh::IO::read_mesh()行为,如何控制Mesh::IO::write_mesh()行为。OpenMesh::IO::Options可以用作读写网格。这个类通过开关相应的bit位实现控制读写器的行为。这个类提供接口实现开、关和修改相应的位。我们区分一下的控制位mode bit-控制二进制读写Options::BinaryOptions

2013-07-20 22:56:24 1577

翻译【OpenMesh】删除几何元素

这个小例子演示如何移除网格中的面和顶点。我们使用第一个例子中的正方体。如果我们希望我们的网格能够移除顶点，面或者边，我们必须为网格类扩展默认特性。顶点，面和Halfedge需要OpenMesh::Attributes::Status属性，这个属性用于当一个元素被删除的时候保持其已删除标记。定义网格特性的结构如下：struct MyTraits : public OpenMesh::

2013-07-20 19:07:17 1426

翻译【OpenMesh】使用特性扩展网格

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/tutorial_07.html例子显示：如何使用特性扩展实体的性能在之前的指南中我们使用属性和改变Point，Normal，TexCoord，和Color数据类型。但是我们可以使用特性（Trait）做更多。我们可以改变网格实体比如顶点，面，边和Halfedge的性

2013-07-20 18:45:28 1116

翻译【OpenMesh】使用网格的属性和特征

例子主要展示如何改变位置，法向量，颜色和纹理的数据类型。在之前的指南中我们学习使用标准属性，通过调用适合的请求方法。不像自定义属性，用户通过传递数据类型到句柄来指定数据类型（比如，MyMesh::FPropHandleT），标准属性的数据类型定义为网格特征。我们可以和特征一起定制和扩展网格数据结构。我们通过两方面做到这一点。改变位置（Position），法向量（Normal），颜色（C

2013-07-20 14:37:23 2022 2

翻译【OpenMesh】使用标准属性

这个例子展示：如何添加标准属性如何获取和设置标准属性的值正如之前所说，我们可以通过属性绑定数据到网格实体中。OpenMesh提供一连串的标准属性。不像自定义属性，标准属性有一些特定的功能和不同的接口，这就是这篇指南所要讨论的。下面的表格列举了所有可能的标准属性和能够使用这些属性适合的实体。 VertexFaceEdgeHalfedge

2013-07-19 23:49:57 2127

翻译【OpenMesh】使用STL算法

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/tutorial.html因为OpenMesh迭代器几乎与STL的迭代器一致，意味着可以对网格使用STL的算法。下列的示例代码展示了如何使用STL算法中的for_each，因为这样写可读性更好而比手写的循环效率更高。我们定义一个类提供平滑算法，因此定义一个可重用的组件。

2013-07-19 22:15:14 1295

翻译【OpenMesh】使用自定义属性

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/tutorial.html这个例子展示：如何添加和移除自定义属性如何取得和设置自定义属性的值再上一个例子中我们计算了每一个顶点的重心并将它们保存在Array中。如果我们将数据保存在网格中并由OpenMesh管理数据，这将会跟方便并且不易出错。如果我们连接这样的特性将

2013-07-19 21:58:07 4216

翻译【OpenMesh】使用迭代器和循环机

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/tutorial.html这个例子展现：如何使用迭代器如何使用循环机这个例子是第一个版本的平滑器。我们将会介绍迭代器（iterators）和循环机（circulators）。这两个概念提供线性枚举的功能，比如网格中所有的顶点，围绕一个顶点循环，枚举所有的1-ri

2013-07-19 18:20:31 1258

翻译【OpenMesh】创建一个正方体

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/tutorial.html这个例程演示了：如何声明MyMesh如何添加顶点和面到Mesh中如何使用IO函数写入网格数据每一个程序第一步都是定义MyMesh类型。OpenMesh支持常见的多边形网格（面由可变顶点数的多边形组成）以及特殊的三角网格（所有面都是三角形的）

2013-07-19 17:56:08 1482

翻译【OpenMesh】如何使用OpenMesh创建项目

在这篇指南中我们将会解释如何使用OpenMesh创建项目并且使用CMake构建系统。我们假定你已经下载了OpenMesh的源代码并安装了CMake构建工具。按照下面的几步添加你自己的应用到构造树中：到OpenMeshRoot/src/OpenMesh/Apps中创建一个新的目录，名为"MyPwnProject"现在创建一个新的文件名为"CMakeLists.txt"包含下列代码：

2013-07-19 17:21:41 1470 1

翻译【OpenMesh】指南（代码示例）

OpenMesh库是处理多边形网格强有力的库。因为OpenMesh的内在创建结构允许用户创建用户定制的网格类型，这种网格类型针对特定需求的应用。用户可以应用自己的数据结构来表示顶点，边，和面，或者便捷地使用预先定义的数据结构。另外，OpenMesh提供动态特性允许用户运行时连接或者中断一个数据。这篇文档介绍OpenMesh的基本概念。在这个指南中，我们通过代码示例介绍OpenMesh库

2013-07-19 17:02:19 1106

翻译【OpenMesh】概念类层次结构

因为这里没有OpenMesh类和大量采用C++模板的类库，我们展示了OpenMesh::TriMesh_ArrayKernelT继承图，作为所有网格类型的代表。请注意！大多数继承关系已经实现了模板参数化！因此有一些继承关系没有在参考继承图中文档化。构建内核BaseKernel定义了基础的属性操作比如 add/remove/access。接下啦，AttribKernelT添加

2013-07-16 16:31:56 969

翻译【OpenMesh】一些基本操作:翻转和折叠

这一节中你将学到一些OpenMesh中早已提供的基础操作。内容包括三角形网格边的翻转以及通过连接邻接的顶点边缘折叠。三角形网格的翻转（Flipping edges）考虑到两个邻接面的三角形网格中存在两种不同内部边缘设置。调用函数OpenMesh::TriConnectivity::flip(EdgeHandle _eh)将翻转特定的边到其相反的方向，如下图所述。下面的代码段

2013-07-16 15:50:02 2802 2

翻译【OpenMesh】文件读取和写入

这一节讲的是从文件中读取网格数据和写入网格数据到文件的方法。响应的函数定义在名字空间OpenMesh::MeshIO中。这一节分成三部分。第一步给出例子讲述如何使用IOManager，第二步讲述IOManager的背景信息，第三步展示如何添加模块到IOManager。第一步IOManager快速开始复制以下代码到工程中。注意：如果你静态链接OpenMesh，你必须添加OM_

2013-07-15 17:50:08 2470 2

翻译【OpenMesh】操作网格

介绍在这篇指南中你会学到如何使用OpenMesh类库操纵网格。在之前的一章中（参考Mesh Iterators and Circulators），你已经学了如何迭代顶点，边，halfedge和面，以及循环遍历某些结构比如1-ring。所以这篇指南中我们将会关注如何高效的使用halfedge数据结构和其中有用的特性比如说边界标志（boundary flag）。我们假定你已经熟悉OpenMesh中

2013-07-15 16:45:48 2952

翻译【OpenMesh】网格的迭代器和循环机

原文出处：http://openmesh.org/Documentation/OpenMesh-Doc-Latest/mesh_iterators.html迭代器Mesh提供线性迭代（枚举定点，halfedge，边和面）。这些功能提供更加便捷的网格导航。每一个迭代器XYZIter迭代也存在常数型的迭代ConstXYZIter。所有的迭代器在名字空间OpenMesh::Itera

2013-07-15 00:44:11 2495 1

翻译【OpenMesh】Halfedge数据结构

这一节描述了下列用于存储mesh实体，定点，边，面和连接信息的数据结构。现在有很多种数据结构可以用来表示多边形网格。在这一节末引用的论文有详细的比较。在这个项目中使用的数据结构叫做halfedge data structrue。face-based数据结构存储连接信息在面上引用其中的顶点和邻接元素，edge-based 数据结构将连接信息放到边上。每一条边引用它的两个顶点，属于包含这条边的

2013-07-14 00:40:16 2782 1

原创【Interactive 3D Graphics】Exercise: Head Light

题目的要求是在场景中加入点光源。 var point = new THREE.PointLight(0xFF0000,3.0); point.position.set(0,0,100); scene.add(point); var point = new THREE.PointLight(0xFF0000,3.0); point.position.set(0,0,10

2013-07-07 11:41:06 577

原创【Interactive 3D Graphics】Exercise: Set a Directional Light

题目是在场景中加入直射光，DirectionalLight我加入红光，设置好位置之后加入场景中。 var light=new THREE.DirectionalLight(0xFF0000,1.5); light.position.set(-200,200,-400); scene.add(light);

2013-07-07 11:15:12 1121

原创【Interactive 3D Graphics】Exercise: Diffuse Material

开始先说说MeshLambertMaterial和MeshBasicMaterial。MeshBasicMaterial:A material for drawing geometries in a simple shaded (flat or wireframe) way.MeshLambertMaterial:A material for non-shiny (Lambertian

2013-07-06 19:32:32 672

原创【Interactive 3D Graphics】Exercise: Vertex Attributes

题目主要是写补全一些代码来实现彩色的三角形。完成的效果如下：代码很简单，只是因为是JavaScript的缘故，真的很容易打错变量。 geometry.vertices.push(new THREE.Vector3(100,0,0)); geometry.vertices.push(new THREE.Vector3(0,100,0)); geometry.vertices.

2013-07-06 17:57:51 610

原创【OpenCV】双边滤波

双边滤波（Bilateral filter）是一种非线性的滤波方法，是结合图像的空间邻近度和像素值相似度的一种折衷处理，同时考虑空域信息和灰度相似性，达到保边去噪的目的。具有简单、非迭代、局部的特点。 Bilateral Filtering for Gray and Color Images可以看看这里的解析，里面讲的很详细。 C++的函数原型:voidbilateralFilte

2013-06-12 15:56:17 181

原创【OpenCV】中值滤波

中值滤波法是一种非线性平滑技术，它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值. 函数原型：void cvSmooth(const CvArr*src, CvArr* dst, int smoothtype=CV_GAUSSIAN, intparam1=3, int param2=0, double param3=0, doubleparam4=0)i

2013-06-12 15:24:29 184