技术标签: localization review WCE wireless capsule endoscopy medical robot
Efficient magnetic localization and orientation technique for capsule endoscopy [1]
Paper Link
Authors: Chao Hu, etc.
2005,International Journal of Information Acquisition
本文是将磁极子模型用于WCE定位领域的一篇重要、基础、经典的文章。将magnetic dipole model用于WCE定位,而magnetic dipole model将单一磁铁看作一个磁极子进行建模,给定磁体的位置和磁矩朝向,模型能计算出磁体外一点的磁场。magnetic dipole model如下式所示:
b = μ 0 M 4 π ∣ ∣ r ∣ ∣ 3 ( 3 ( m ^ ⋅ r ^ ) r ^ − m ^ ) \textbf{b}=\frac{\mu_{0}M}{4\pi||\textbf{r}||^{3}}\left(3(\widehat{\textbf{m}}\cdot\widehat{\textbf{r}})\widehat{\textbf{r}}-\widehat{\textbf{m}}\right) b=4π∣∣r∣∣3μ0M(3(m ⋅r )r −m )
这里的 μ 0 \mu_{0} μ0是真空磁导率(常数), μ 0 = 4 π × 1 0 − 7 \mu_{0}=4\pi\times10^{-7} μ0=4π×10−7; M M M是一个与磁体体积材料有关的常量, M = V ⋅ M 0 M=V \cdot M_{0} M=V⋅M0,V是磁体的体积, M 0 M_{0} M0是磁体的材料,矫顽力xxx kA/m; r \textbf{r} r表示磁体中心点到磁体外一点的向量, r = ( p s x , p s y , p s z ) T − ( p c x , p c y , p c z ) T \textbf{r}=(p_{sx},p_{sy},p_{sz})^{T}-(p_{cx},p_{cy},p_{cz})^{T} r=(psx,psy,psz)T−(pcx,pcy,pcz)T; m ^ \widehat{\textbf{m}} m
表示磁体的单位磁矩方向; r ^ \widehat{\textbf{r}} r
表示 r \textbf{r} r的单位向量(方向)。
若共有 N N N个sensor,那么优化如下非线性目标函数:
m i n ∑ s = 1 N ∣ ∣ b s a m p l e ( s ) − b m o d e l ( s ) ∣ ∣ , s = 1 , 2 , ⋯ , N . min\sum_{s=1}^{N}||\textbf{b}_{sample}^{(s)}-\textbf{b}_{model}^{(s)}||, s=1,2,\cdots,N. mins=1∑N∣∣bsample(s)−bmodel(s)∣∣,s=1,2,⋯,N.
最后的优化结果就是磁体的位置 ( p c x , p c y , p c z ) T (p_{cx},p_{cy},p_{cz})^{T} (pcx,pcy,pcz)T和磁矩朝向 m ^ \widehat{\textbf{m}} m 。
在这篇文章中,也对多个非线性算法对此问题的求解性能进行了分析,最终发现LM算法的性能最优,也成为了纯磁定位通常选用的方法。
作者采用了三轴折线图来表示一段时间的位置和朝向误差,以及在3D图和俯视2D图中,直观显示位置误差,最终位置和角度误差为 5.6 m m 5.6mm 5.6mm, 4.26 % 4.26\% 4.26%。
Efficient linear algorithm for magnetic localization and orientation in capsule endoscopy [2]
Paper Link
Authors: Chao Hu, etc.
2005,Engineering in Medicine and Biology 27th Annual Conference
在前一篇paper中,使用了非线性优化算法,计算量大,耗时长。作者在这篇文章中,提出了一种线性解析算法,一步一步根据采样磁场反推出磁体的位置和朝向。具体方法为,列出磁极子模型:
b = μ 0 M 4 π ∣ ∣ r ∣ ∣ 3 ( 3 ( m ^ ⋅ r ^ ) r ^ − m ^ ) \textbf{b}=\frac{\mu_{0}M}{4\pi||\textbf{r}||^{3}}\left(3(\widehat{\textbf{m}}\cdot\widehat{\textbf{r}})\widehat{\textbf{r}}-\widehat{\textbf{m}}\right) b=4π∣∣r∣∣3μ0M(3(
文章浏览阅读6.3k次。日志默认info级别debug日志不会打印,但是会执行日志填充的数据例如:logger.debug("日志输出",2*10); 1. 2*10会先执行出结果,然后继续往下走2. 在ch.qos.logback.classic.Logger#filterAndLog_1方法中判断是否符合级别要求是否需要输出3.如图:..._debug中的计算是否在info级别也会跑
文章浏览阅读1.4k次。Similarly to the previous example, let us apply our calibration engine onto the data that comes with the originalcalibration toolbox of Heikkil� from the University of Oulu. Once again. do not bothe_non-planar calibration
文章浏览阅读1w次,点赞10次,收藏63次。物联网常用的网络协议:MQTT、AMQP、HTTP、CoAP、LwM2M物联网设备间沟通的语言,就是网络协议。设备间想相互交流,通信双方必须使用同一种“语言”。比如说你和中国人问好说’你好‘、日本人问好要说‘こんにちは’、和英国人问好要说‘hello’.说起网络协议,你可能马上就想到了 HTTP 协议。是的,在日常的 Web 开发中,我们总是需要跟它打交道,因为 HTTP 协议是互联网的主流网络协议。类似地,应用在互联网中的网络协议,还有收发电子邮件的 POP3 、SMTP 和 IMAP 协议,以及_lmm2m和mqtt
文章浏览阅读7.4k次,点赞4次,收藏20次。这篇博文简要记录一下使用MKL函数库计算一般矩阵的特征值与特征向量对形如对称矩阵或是埃尔米特等特殊矩阵有其对应的子程序,在这里先不涉及。有需求的可以自行查阅MKL官方文档下面给出本次示例代码:代码使用f95接口。f77借口参数太多,笔者太懒<不过懒惰是创新的原动力^_^>program testGeev use lapack95 implicit..._fortran求矩阵特征值
文章浏览阅读147次。学习内容来自:Numpy Tutorial文章目录Array SlicingArray IndexingMathematical ManipulationBroadcastingImage Processing基本的用法课程里面说的挺详细了。 特别记录一些需要关注的点。Array Slicing使用固定数字进行array寻址会导致数组降维。y = np.random.random((3,..._np.imresize
文章浏览阅读355次。题目阅览 观察数字:12321,123321 都有一个共同的特征,无论从左到右读还是从右向左读,都是相同的。这样的数字叫做:回文数字。 本题要求你找到一些5位或6位的十进制数字。满足如下要求: 该数字的各个数位之和等于输入的整数。 输入格式 一个正整数 n (10<n<100), 表示要求满足的数位和。 输出格式若干行,每行包含一个满足要求的5位或6位整数。 数字按从小到大的顺序排列。 如果没有满足条件的,输出:-1样例输入144样例输出199899_c++蓝桥杯 回文数
文章浏览阅读6.2k次,点赞3次,收藏13次。需要的pom文件 <dependency> <groupId>com.google.zxing</groupId> <artifactId>core</artifactId> <version>3.1.0</version>_java扫二维码进入自己制作的网页
文章浏览阅读650次。在遥感图像处理中,我们经常需要将多波段遥感影像拆分成多个单波段图像,以便进行各种分析和后续处理。本篇博客将介绍一个用Python编写的程序,该程序可以读取多波段遥感影像,将其拆分为单波段图像,并保存为单独的文件。本程序使用GDAL库来处理遥感影像数据,以及NumPy库来进行数组操作。结果如下图所示,选中的影像为输入的多波段影像,其他影像分别为拆分后的多波段影像。_一个多波段影像分解成多个单波段影像
文章浏览阅读5.1k次,点赞2次,收藏4次。0前言一直用的好好的移动硬盘突然不显示了,前段时间因为比较忙,一直没顾得上管它,趁这个假期,好好捅咕了一番,总算是弄好了,特此将解决的过程记录如下:1.问题描述 1.我的移动硬盘在其他人的电脑上能够正常显示和使用 2.其他移动硬盘在我电脑上能够正常的显示和使用 3.在我的电脑上,该移动硬盘,既不显示盘符,磁盘管理 又不显示该磁盘2.问题分析1.我的移动硬盘能够在其他人电脑上_电脑无法显示移动硬盘
文章浏览阅读1k次。Kernel initialization. Part 10.在原文的基础上添加了5.10.13部分的源码解读。End of the linux kernel initialization processThis is tenth part of the chapter about linux kernel initialization process and in the previous part we saw the initialization of the RCU and stopped o_linux 标志着kernel启动完成
文章浏览阅读5.3k次,点赞5次,收藏23次。Scala语言概述:Scala语言是一门以Java虚拟机为运行环境,支持面向对象和函数式编程的静态语言,java语言是面向对象的,所以代码写起来就会相对比较模块儿,而函数式编程语言相对比较简洁_scala安装及环境配置
文章浏览阅读2.4k次。“他来听我的演唱会,门票换了手铐一对”。最近歌神张学友变阿SIR,演唱会上频频抓到罪犯,将人脸识别技术又一次推到了大众的视线中。要说人脸识别技术的爆发,当属去年9月份苹果iPhone x的发布,不再需要指纹,只需要扫描面部就可以轻松解锁手机。任何技术一旦进入智能手机这个消费市场,尤其是被苹果这个标志性的品牌采用,就意味着它将成为一种趋势,一个智能设备的标配。在智能手机快速崛起的这几年,其密码锁..._人脸识别发展历史