由于我在用SPH方法仿真时用的是FORTRAN语言,并且没有找到直接输出vtk文件的代码,因此偷懒通过MATLAB转换一下数据. 用到的Matlab子程序可通过一下链接找到. Matlab处理数据导出Paraview可读的vtk文件(一) 假设我每个粒子输出的格式按照下面输出的 x(i), y(i), z(i), vx(i), vy(i), vz(i), p(i), rho(i), phase(i) 含义为每个粒子的坐标,速度,压强,密度,相. vtkwrite(‘output.vtk’,'un

Paraview是一个开源的可视化软件. 用到matlab子程序从这里下载 或者到博客末尾复制粘贴 子程序名为 vtkwrite 示例1: load mri D = squeeze(D); vtkwrite('mri.vtk', 'structured_points', 'mri', D) 示例2:云图 load wind [cu,cv,cw] = curl(x, y, z, u, v, w); div = divergence(x, y, z, u, v, w); vtkwrite('wind

matlab算法转为c语言后,影响c语言效率的关键在于multiword的产生,基于此会有multiword加减法和乘除法,极大消耗资源,减少甚至消除multiword很重要,需注意的是:算法中尽量减少除法使用,如果用到乘法,减少乘法项的个数,尽量都弄成两个数的乘法,对于矩阵算法尤其要注意乘除法的使用.矩阵算法可以在最后要得到结果时再用除法或乘法.多用加减法. 多步算法变为定点算法后很可能误差指数倍增长.

MyYuanLaiPic = imread('e:/image/matlab/Cluo.jpg');%读取RGB格式的图像 MyFirstGrayPic = rgb2gray(MyYuanLaiPic);%用已有的函数进行RGB到灰度图像的转换 [rows , cols , colors] = size(MyYuanLaiPic);%得到原来图像的矩阵的参数 MidGrayPic = zeros(rows , cols);%用得到的参数创建一个全零的矩阵,这个矩阵用来存储用下面的方法产生的灰度图

今天是第一天写博客,哈哈哈!把完成的数字图像作业放上来和大家一起分享一下! 如果有什么问题,希望大家和我多多交流 1518234852@qq.com width=512; height=512; img = fopen('RAW.raw');% 打开原图像%同文件夹下 m = fread(img,[width height],'uchar');% uchar为无符号字符型 m = uint8(m);% unit8表示无符号整数,范围0-255 mt= m';% 将点阵转置(这是由于BMP格式图像

直接保存为txt文件: 可以用fprintf函数,来代替save函数 比如现在我有一个变量a=[0.1223   345.4544] 如果我想保存它的话,可以用下面的程序: fid = fopen('a.txt','wt'); fprintf(fid,'%6.4f   %8.4f\n',a); fclose(fid); 然后打开a.txt,你会发现其中是0.1223   345.4544如果用save a.txt -ASCII a 保存的话,a.txt中的内容是   1.2234500e-001

前言:本程序是我去年实现论文算法时所做.主要功能为标记切割肝脏区域.时间有点久,很多细节已经模糊加上代码做了很多注释,因此在博客中不再详述. NOTE: 程序分几大段功能模块,仔细阅读,对解决医学图像还是有一定的借鉴意义 想借鉴本文的一定要仔细阅读代码和注释,中间有人机交互部分,空跑会抛异常 .dcm数据,我放到了我的百度云盘,有兴趣的可以下载,实测一下代码.dcm数据连接 clc,clear img_1=dicomread('10011.dcm');%读取dcm文件 (所谓的灰度值) meta

容易忘记的小知识点: Matlab程序首行程序: clear:close all:clc:程序运行开始最好清理下之前残留的各种数据,不然matlab可能会出错. clear(清理变量空间),close all(关掉所有显示的图像figure),clc(清理命令窗口中的显示) Matlab程序换行(一个空格加上3个.即 ...然后直接下一行就可以继续写代码) Matlab下程序计时 如下:tic和toc可以在程序中使用多个,计算各个程序字块的运行时间.还可以嵌套使用... tic;%计时开始 %%

将srcimg文件下的bmp文件转为jpg图像,存放在dstimg文件夹下 str = 'srcimg'; dst = 'dstimg'; file=dir([str,'\*.bmp']); :length(file) img = imread( [str,'\',file(i).name]); imwrite(img, [dst,'\',file(i).name(1:end-4),'.jpg']); end

PLY格式是比较流行的保存点云Point Cloud的格式,可以用MeshLab等软件打开,而VTK是医学图像处理中比较常用的格式,可以使用VTK库和ITK库进行更加复杂的运算处理.我们可以使用ParaView软件对VTK格式文件进行预览和简单处理,ParaView也可以打开PLY格式,但是就没有texture了,而且我们如果直接用ParaView导出VTK格式也没有texture的,这不是我们想要的结果.MeshLab虽然可以打开有texture的PLY文件,但是却不支持导出VTK格式,那么我

常用功能的实现 获取当前脚本所在目录 current_script_dir = fileparts(mfilename('fullpath')); % 结尾不带'/' 常用函数的说明 bsxfun matlab函数 bsxfun浅谈 bsxfun是一个matlab自版本R2007a来就提供的一个函数,作用是"applies an element-by-element binary operation to arrays a and b, with singleton expansion enab

转载:http://blog.sina.com.cn/s/blog_6c41e2f30101559d.html ****************假设某图像数据A(uint8格式)********************** A = 235  200  89  20 >> double(A)                 %返回与原矩阵数值相同但类型为double的矩阵: ans = 235   200 89   20 >> im2double(A)           %返回矩阵

阅读目录 前言 Matlab的开发环境与简单介绍 Matlab的常见命令 Matlab的灵魂-矩阵操作 Matlab的.m或.fig的编程与技巧 从大二开始接触到matlab,讲真,这是一个我觉得很良心的一个平台,之所以这样认为是因为matlab有个help功能,这里可以查询到所有的函数以及关键字的使用方法,但是可惜是英文版的,就比较难找到自己想要的函数功能,长时间里自己都是在百度里将自己所需要的函数功能说出,然后百度,一般都能找到,我觉得很有用的一个matlab论坛平台是MATLAB中文论坛,

0 系统+软件版本 系统:CentOS 6.7 x64, 内核 2.6.32-573.el6.x86_64软件:Matlab R2015b(包括威锋网和东北大学ipv6下载的资源,都测试过) 1 脚本运行“陷阱“ 1.1 未知的“陷阱” 首先,这个程序在Matlab R2013a中可以完美运行,这个“陷阱“在是新安装的R2015b上才出现的. 说它是“陷阱“,是因为脚本文件涉及到大量的数据处理,比如一个几百次的循环,它可能在执行某一句的时候就失去响应了,可能是一个循环,也可能是单句,仿佛掉进了一

1 设计内容 本设计是基于FPGA的音频信号FIR低通滤波,根据要求,采用Matlab对WAV音频文件进行读取和添加噪声信号.FFT分析.FIR滤波处理,并分析滤波的效果.通过Matlab的分析验证滤波效果后,将叠加噪声信号的音频信号输出到txt文件里.然后使用Matlab语言编写滤波器模块和测试模块,通过Modelsim软件读取txt文件的数据,将数据送入滤波模块,最后将滤波的结果输出到txt文件里,最后用Matlab将处理的结果从txt文件读出.显示.FFT分析用Verilog设计的FIR滤

上上篇的博客已经提到如何配置环境,即如何在asp.net中调用matlab生成的dll文件.这篇博客打算做个笔记,那就是matlab和C#数据类型如何转换.随着需求的增加,我会不断增加新的类型转换. (1)C#二维数组转为Matlab矩阵 //这段代码主要就是将C#的二维数组转为Matlab的矩阵 double[,] frontDrawDataArray = new double[250, 290]; double[,] sideDrawDataArray = new double[250, 2

最近的项目需要matlab和C的混合编程,经过一番努力终于完成了项目要解决的问题.现在就将Mex的一些经验总结一下,当然只是刚刚开始,以后随着学习的深入继续添加.首先讲讲写Mex的一些常规规定,然后我们会重点关注混合编程中最难解决数据的问题--结构到底如何转换,并且后面会重点说一下自己的程序. 一.Mex的结构 先看一个简单的程序(该程序保存在matlab主目录下名字是mexDemon.cpp,或者在主目录下新建一个.cpp文件): #include "mex.h" //加入头文件,该

一.常用对象操作:除了一般windows窗口的常用功能键外.1.!dir 可以查看当前工作目录的文件. !dir& 可以在dos状态下查看.2.who 可以查看当前工作空间变量名, whos 可以查看变量名细节.3.功能键:功能键 快捷键 说明方向上键 Ctrl+P 返回前一行输入方向下键 Ctrl+N 返回下一行输入方向左键 Ctrl+B 光标向后移一个字符方向右键 Ctrl+F 光标向前移一个字符Ctrl+方向右键 Ctrl+R 光标向右移一个字符Ctrl+方向左键 Ctrl+L 光标向左移

安装 1. 在http://www.csie.ntu.edu.tw/~cjlin/ 中下载libsvm 2. 按照http://zjhello123.blog.163.com/blog/static/25495143201251792625811/ 中的安装方法安装libsvm 3. 设置路径. 如果不设置的话会报错:Error using svmtrain (line 233) Y must be a vector or a character array. 我开始设置了也报错,后来重新mak

将matlab中数据输出保存为txt或dat格式 总结网上各大论坛,主要有三种方法. 第一种方法:save(最简单基本的) 具体的命令是:用save *.txt -ascii x x为变量 *.txt为文件名,该文件存储于当前工作目录下,再打开就可以 打开后,数据有可能是以指数形式保存的. 例子: a =[17 24 1 8 15;23 5 7 14 16 ;4 6 13 20 22 ;10 12 19 21 3 ;11 18 25 2 9 ]: save afile.txt -ascii a