tag 标签: 系统经管大学堂:名校名师名课

相关帖子

版块 作者 回复/查看 最后发表
用定理证明器构造条件计划 外文文献专区 kedemingshi 2022-3-8 0 411 kedemingshi 2022-3-8 18:49:20
[分享]世界统计与分析全才SAS系统应用指南(上、下册) attachment 数据分析与数据科学 gxg365 2006-6-17 610 53160 hyq2003 2020-10-30 06:10:17
基于MATLAB的系统分析与设计-小波分析 attachment MATLAB等数学软件专版 lee2010wl 2013-1-23 9 3855 butterfly1980 2019-7-30 17:41:53
[下载]Aptech GAUSS 9.0.0.1108 (数学和统计系统) attachment Gauss专版 xuehe 2008-6-24 42 24959 tianwk 2019-3-18 18:49:43
[免费 下载]DPS数据处理系统教程(2006最新版) attachment 计量经济学与统计软件 yangyang9189 2007-10-31 156 51497 岳岳岳岳岳 2016-10-10 10:38:27
[公告]人大经济论坛新版系统测试 站务与外事 angelboy 2009-5-3 115 39176 leon_9930754 2016-4-22 06:34:17
丰田管理系统 模式全书 attachment 运营管理(物流与供应链管理) 未醉的狐狸 2013-2-6 1 4615 shlywtz 2015-9-5 20:46:55
[分享]复杂适应系统中智能体的行为规则及程序实现 数据交流中心 tybjsky 2005-11-27 9 4481 yifeng1025217 2015-4-16 16:05:55
information systems management in practice信息系统管理实践课件ppt attachment 管理信息系统 simmba 2013-1-25 6 4848 tigerbaby 2015-3-12 09:57:43
海氏岗位价值评价系统 attach_img 人力资源管理 rxy111666 2013-2-18 5 7486 antiloo 2014-4-2 12:48:38
马克威分析系统单机版的学员版(4.0) attachment 数据分析与数据挖掘 事为 2009-3-7 17 6710 frienddxm 2013-4-2 11:05:59
论坛问卷系统正式启用www.wenjuan.cc域名 attach_img 调查问卷专版 admin 2013-1-18 53 10675 /狂奔的蜗牛 2013-2-27 22:44:36
怎么删除帖子 - [阅读权限 30] 学术道德监督 kobecass 2013-2-19 6 536 QYZ2012 2013-2-19 11:16:17
怎么删除帖子呢? 教师之家与经管教育 kobecass 2013-2-19 3 2980 jinxi118 2013-2-19 10:09:26
[轉貼] 现代SCM系统功能特點 运营管理(物流与供应链管理) Toyotomi 2013-2-15 0 4997 Toyotomi 2013-2-15 00:14:43
bys 命令老出错? Stata专版 peyzf 2013-1-27 4 8360 peyzf 2013-1-27 10:12:28
SAS 9.1.3 在Vista系统下的安装笔记 SAS专版 gzndxf 2008-10-10 17 8751 kyolee2006 2012-11-19 12:06:05
【国内资源】[4-12]中国统计数据应用支持系统 国内外文献账号区 yinhezhiwang 2009-4-12 22 4872 日尧光 2012-11-6 12:02:22
【国内文献】[2-23]中国统计数据应用支持系统 国内外文献账号区 yzhongb 2009-2-23 4 2034 544964487 2011-11-15 18:28:58

相关日志

分享 CT系统参数标定问题研究
斯诺帝亚2 2012-8-21 10:22
摘 要 本文通过对体模系统和投影点坐标的分析,提出了较为充分合理的假设,并对数据进行了拟合处理。先后建立了三个模型(理想状态下的模型,只考虑载物台倾斜的模型,考虑载物台倾斜与探测器的 俯仰角的模型) 来求解 CT 系统参数标定的问题。 在对问题一的分析中, 只考虑在理想情况下:即载物台水平和探测器无偏转,并不考虑机械误差的情况下各个参数之间的关系 ,以 X 光放射源为原点建立坐标系, 沿水平方向射向探测器的水平直线为 X 轴,竖直方向为 Y 轴, 建立模型一, 绘出侧视图如一所示 。运用几何知识作图,通过建立相似图形的等比例方程来研究确定几个系统参数之间的关系( X 射线源的位置参数、载物中心线的位置参数和探测器的位置参数)。同时求出载物台半径等参数,最后讨论了引入像素、载物台偏角、观测器俯仰角和数据不确定度对误差的影响。 在分析问题二时,考虑到所给椭圆点集可能有偏差,于是首先运用拟合的方法对数据行处理,得到更为确信的数据并由此求出修正后的“椭圆”方程(观察均方根误差,得知此方程能符合条件,可以使用)。由于考虑多个误差因素将使模型计算变得极为复杂,在此不同时考虑。假设 探测器无偏转, 只有载物台倾斜,,建立模型二。易知倾斜的载物台与安装没有偏差的水平观测器可以等效为水平的载物台与只纵向倾斜的观测器。然后根据拟合后的数据求出两类椭圆长轴斜率。从而求出载物台倾斜角为 0.1146︒ ,并依据此进一步修正数据。 在模型二修正数据的基础上,考虑探测器的俯角 因素(易知 不影响模型二的求解)建立模型三。由修正后的图知两“椭圆”长轴长等。由上方长轴长大于下方长轴长知探测器向前倒。经过一定数据近似后,再根据等比例关系可求得探测器俯角 为 。 近一步推导可得坐标及光源到旋转台的旋转轴的距离为 ,光源在成像平面上的水平直射点,光源到探测器水平直射点的间距为 ,。 最后在问题三,对模型三的机械误差进行了分析计算。机械误差主要来源为像素、载物台偏角,探测器的俯仰角,圆心距离的误差。并进行了灵敏度分析。 关键词:坐标修正 载物台偏角探测器俯仰角 机械误差 一、问题的重述 扫描物件围绕某一固定转轴旋转,每隔一定角度采集一张图像,然后根据采集的图像采用 3D 图像重建算法即可将原始 3D 物件重建出来。在实际中,由 X 射线源发出 X 射线,经过扫描物体衰减后照射在探测器上,探测器根据接收到的光子数的计数实现光电转换,从而形成灰度图像。一般平板探测器大小为 3000 × 2000 像素,每个像素为 0.127mm 。 在 CT 系统安装过程中往往存在机械误差,而这些误差对于物件重建的准确性往往是致关重要的,实际中就需要对安装好的 CT 系统进行参数标定,系统参数主要包括 X 射线源的位置、载物中心线的位置和探测器的位置参数等。对于 CT 系统参数标定精度要求是不够的(通常 CT 重建过程中系统参数的机械误差不允许超过一个像素)。通过扫描已知参数的体模,分析投影数据,估计系统参数的方法。针对工业 CT 要求解决如下问题: ( 1 )建立合适的坐标系,正确描述 CT 系统的各种参数和机械误差,并建模分析这些参数的关系和可能的机械误差。 ( 2 )通常采用轴承钢球作为 CT 参数标定体模。因为它具有各向投影一致,边缘清晰等优点,在里采用两个钢球实现 CT 的参数标定实验 . 两个钢球置于有机玻璃管中( X 射线容易透过有机玻璃,容易后期球心获取的图像处理),钢球直径为 8mm ± 0.0008mm ,两球心距离约为 100mm ± 1mm 。将有机玻璃管固定在旋转载物台上,载物台携带有机玻璃管以均匀速度旋转,与此同时 X 射线源发出的射线,探测器采集数据(旋转一周,等间距采集 180 张)。两个钢球球心投影的轨迹是两个椭圆。试根据实验数据估计该 CT 系统的参数值,给出 CT 系统的标定。 ( 3 )在参数标定过程中可能存在多种可能的机械误差,对参数标定可能的误差进行分析。 二、问题的分析 问题一:本题误差主要来源于工业 CT 系统安装过程。除了装置本身的由于在装系统时,可能会导致探测器扭转、前后俯仰,载物台的倾斜等,都会导致机械误差的存在。我们在计算时会考虑这些因素 , 并建立相应的模型。所以,在标定时应该给出其机械误差大小。可以根据体模的投影数据,利用几何知识来求解 主要包括 X 射线源的位置、载物中心线的位置和探测器的位置并求出他们的关系。对于载物台来说,由于载物台要进行旋转,所以其倾斜方式只有一种。 问题二:根据分析,载物台的倾斜可以等效为探测器投影绕一个点在竖直方向旋转(图三所示)。由附录一中图分析可知,,两个椭圆的长轴不是水平的,这是由于载物台的倾斜造成的;上椭圆与下椭圆的长轴不一样长,这是由于探测器的前后俯仰造成的;“椭圆”图形左右不对称,是由于探测器的扭转造成的。我们分别建立模型将数据一步一步进行处理,得到机械误差来源,再进行修正,从而给出准确的 CT 参数标定。 问题三:本题中考虑的机械误差来源有两圆心距离误差,载物台倾斜,探测器前俯后仰,利用问题二中所得数据,可求得总机械误差。 三、模型的假设与符号说明 1. 模型的假设 ( 1 )考虑到此处光子的波动不明显,故忽略光的衍射现象,以光的直线传播进行处理; ( 2 )忽略 X 光辐射源的体积和形状,将其视为质点; (3) 一些题目中未说明但实际易得的参数,如载物台的半径等,假设可以测量得到,此处不考虑具体方法及相关误差; ( 4 ) 数据能客观反映实际情况; ( 5 )载物台与有机玻璃管始终垂直; ( 6 )两小球圆心视为质点; ( 7 )在数据采集过程中,放射源点,载物台,探测器参数不再变化。 2. 符号说明 : X 光射线源到载物台中心垂线的距离; : X 光射线源到探测器的距离; :探测器俯仰偏转后下方椭圆最下沿点到 Y 轴距离; : 探测器俯仰偏转后上方椭圆最上沿点到 Y 轴距离; :有机玻璃管和小球的轨迹圆半径; :两个小球重心之间的距离; :椭圆间最近两点的距离; :椭圆间最远两点的距离; :上投影椭圆短轴长度; :上投影椭圆短轴长度; :投影的竖直方向最大长度; :载物台的倾斜角; :探测器偏转的俯仰角; :大椭圆的长轴长度; :小椭圆的长轴长度; :两椭圆中心之间的距离; :当偏差为一个像素时候,探测器的俯仰角度; :当偏差为一个像素时候 的变化量; :当偏差为一个像素时候 的变化量; :当偏差为一个像素时候 的变化量 四、模型的建立和求解 问题一 : 1 、模型一的建立 首先考虑在理想情况下,即载物台水平和探测器无偏转和不考虑机械误差的情况下各个参数之间的关系。建立坐标系时选定 X 光辐射源为坐标原点,沿水平方向射向探测器的水平直线为 X 轴,竖直方向为 Y 轴,绘出侧视图如一所示。 图一理想情况下坐标系的建立 图二理想模型下的俯视图 由平面几何中相似图形成比例可以得出下列方程: (1) ; (2) ; (3) 若考虑引入机械误差,则可能存在以下几种机械误差: ( 1 )载物台与地面有一定倾斜角,此角会使投影偏转。 ( 2 )探测器地面间存在一定俯仰角,此偏角的存在会使原本应是椭圆的投影产生一定偏差; ( 3 ) x 放射点与载物台轴线成面与探测器成角不成 90 度,会使原本应是椭圆的投影产生一定偏差; ( 4 )像素不够大引起的误差,一般平板探测器大小为 3000 2000 ,每个像素为 0.127 ,若两个投影点之间的水平或竖直距离小于 0.127 ,在探测器上就会区分不出,例如以下 2 组数据: 上方球心坐标 ( X , Y ):( 1415 , 1639 ),( 1375 , 1639 ); ( 5 )题目中给出的数据(小球半径、球心距离等)不精确产生的误差; 问题二 : 1、 模型二的建立 暂 只考虑载物台的倾斜角 的情况下,建立此模型。 图三探测器的偏转 倾斜的载物台与安装没有偏差的观测器可以等效为水平的载物台与只纵向倾斜的观测器,如图二 由于数据可能有较大的误差,我们用 matlab 对数据进行拟合。 图四:修正前的图形 图五:拟合修正后的图形 再由图 中所示方法求出长轴短轴大小及对应各个切点。 图六:椭圆的切线与椭圆轴的关系 再根据对应切点确定长轴斜率 tanx=(lxj-rxj)/(lx-rx); % 下方椭圆长轴斜率 tans=(lsj-rsj)/(ls-rs); % 上方椭圆长轴斜率 tan=(tans+tanx)/2; 得 根据相应切点求出圆心轨迹所成圆之圆心,水平射于探测器的点横纵坐标 xs0=(ls+rs)/2; ys0=(lsj+rsj)/2; % 上方圆心轨迹圆圆心水平射于探测器的点横纵坐标 xx0=(lx+rx)/2; yx0=(lxj+rxj)/2; % 下方圆心轨迹圆圆心水平射于探测器的点横纵坐标 bili=(gs-ds)/((gs-ds)+(gx-dx)); 根据附件一中所给数据,计算出 x 放射点水平射于探测器的横纵坐标 y0=ds-(gx-ds)*bili;%x 放射点水平射于探测器的纵坐标 x0=dsi-(gxi-dsi)*bili; %x 放射点水平射于探测 图七 修正的原理图 再根据对应切点确定长轴斜率 h=abs ( sqrt(tan^2+1)-(tan*(x(i,j)-ls)-y(i,j)+lsj)/sqrt(tan^2+1) ) ; % 上方椭圆上任意一点至上方椭圆长轴的距离 y(i,j)=y0+(ys0-y0)*sqrt(tan^2+1)-(tan*(x(i,j)-ls)-y(i,j)+lsj)/sqrt(tan^2+1); % 修正后的纵坐标 x(i,j)=x0+fuhao*hh; % 修正后的纵横坐标 修正后图形如图五。详细数据见附录一。 2、 模型三的建立 我们可以建立在模型二的基础上考虑探测器的俯仰角度 ,如图八所示,由于载物台的偏转角度 与探测器俯仰角度 不影响,所以可以把他们结合起来计算。考虑他们偏转角的误差。 由图八和图可知,探测器向前偏转了 ,导致上椭圆和下椭圆大小不一样,根据几何知识,可以求得方程 (4) 和方程 (5) : 于是,我们就可以求出在探测器前后俯仰时候的偏转角 与体模投影点坐标之间的关系。找出偏转角 与其他参数之间的关系,进行修正。得到准确的载物中心线的位置和探测器的位置参数,较精确对 CT 进行标定。 图八模型三主视图 图九模型三的左视图 而 与 的值非常相近,即 非常小。相似三角形对应边要小于对应 , 但由与 非常小,相差不大,二者可约等于。 则由图三和图四可以得到下列方程如下: ; (4) ; (5) 非常小。则可以将 , , 的长度看成一样。即式子可改为: ; (6) ; (7) 于是,我们就可以求出在探测器前后俯仰角 。 由模型三的修正,我们可以得到准确的数据,就可以求得 : 带入模型一的方程,由模型一, X 射线源的位置在原点,即坐标为 ,由射线源到载物中心线的距离为 , 即位置坐标为 ,由射线源到探测器的距离为 , 即位置坐标为 , 我们还求出在载物台上物体绕载物台中心线旋转的半径,由于小球在载物台边缘,则旋转半径约等于载物台半径 问题三 : 根据前面的描述,误差的来源主要有三种:像素、偏转角和球心距。 对于 CT 系统参数标定精度要求是其机械误差不超过一个像素,而由式( 1 )( 2 )( 5 )知球心距只与 有直接联系,与像素无直接联系,故此处主要研究偏转角与投影像素之间的关系。有式 (5) 可进一步推出它们的关系如下: ; ( 8 ) ; ( 9 ) ; ( 10 ) 计算( 8 )得到 ; ; 计算( 10 )得到 ; ; 比较两次结果,取最小值 ,即要保证投影在探测板上的点位误差不超过一个相位,探测板的俯仰角浮动范围最大不能超过 。 五、模型的推广与改进方向 考虑三维直角坐标系,能更精确的求出机械误差,具有更好的实际应用。考虑探测器的扭转角对体模投影图像的影响。考虑多种误差同时影响 。 六、模型的优缺点 优点: 1、 建立多个模型,模型的建立由易到难,比较有层次,能更清晰的分析问题,找到它们之间的关系; 2、 对体模的投影点坐标进行数据的修正,计算的结果更加准确; 3、 考虑到的误差比较合理,有较大的实际意义; 缺点: 1、 对探测器的扭转角考虑不全面,仅在理论上进行分析。 2、 由于像素造成的机械误差比较大,像素精度不高,此方法解决这个问题效果欠佳; 七、参考文献 姜启源 . 数学模型 . 北京:高等教育出版社, 1999. 韩中庚 . 数学建模方法及其应用 . 北京:高等教育出版社, 2009. 八、附录 附录一:投影点坐标的修正 zhen=zeros(2000,3000); for i=1:180 for j=1:2 zhen(y(i,j),x(i,j))=1; end end % 矩阵之生成 rs=1000; lx=1000; rx=1000; for i=1:3000 for j=1:2000 if zhen(j,i)==1 if ils if j1400 ls=i; lsj=j; end end if ilx if j500 lx=i; lxj=j; end end if irs if j1400 rs=i; rsj=j; end end if irx if j500 rx=i; rxj=j; end end end end end gs=1800; ds=1800; for i=1:3000 for j=1:2000 if zhen(j,i)==1 if j1000 if jds ds=j; dsi=i; end end end if zhen(j,i)==1 if jgs gs=j; gsi=i; end end end end gx=200; dx=200; for i=1:3000 for j=1:2000 if zhen(j,i)==1 if jdx dx=j; dxi=i; end end if zhen(j,i)==1 if j1000 if jgx gx=j; gxi=i; end end end end end tanx=(lxj-rxj)/(lx-rx); % 下方椭圆长轴斜率 tans=(lsj-rsj)/(ls-rs); % 上方椭圆长轴斜率 tan=(tans+tanx)/2; xs0=(ls+rs)/2; ys0=(lsj+rsj)/2; % 上方圆心轨迹圆圆心水平射于探测器的点横纵坐标 xx0=(lx+rx)/2; yx0=(lxj+rxj)/2; % 下方圆心轨迹圆圆心水平射于探测器的点横纵坐标 bili=(gs-ds)/((gs-ds)+(gx-dx)); y0=ds-(gx-ds)*bili; %x 放射点水平射于探测器的纵坐标 x0=dsi-(gxi-dsi)*bili; %x 放射点水平射于探测器的横坐标 for j=1:2 for i=1:180 if (y(i,j)1500) h=(tan*(x(i,j)-xs0)-y(i,j)+ys0)/sqrt(tan^2+1); hh=sqrt((x(i,j)-xs0)^2+(y(i,j)-ys0)^2-h^2); shu=(x(i,j)-dsi)*(ds-gx)/(dsi-gxi)+ds-y(i,j); fuhao=abs(shu)/shu; y(i,j)=y0+(ys0-y0)*sqrt(tan^2+1)-(tan*(x(i,j)-ls)-y(i,j)+lsj)/sqrt(tan^2+1); x(i,j)=x0+fuhao*hh; end end end % 载物台偏角修正方程 for j=1:2 for i=1:180 if (y(i,j)500) h=(tan*(x(i,j)-xx0)-y(i,j)+yx0)/sqrt(tan^2+1); hh=sqrt((x(i,j)-xx0)^2+(y(i,j)-yx0)^2-h^2); shu=(x(i,j)-dsi)*(ds-gx)/(dsi-gxi)+dx-y(i,j); fuhao=abs(shu)/shu; y(i,j)=y0-(y0-yx0)*sqrt(tan^2+1)-(tan*(x(i,j)-lx)-y(i,j)+lxj)/sqrt(tan^2+1); x(i,j)=x0+fuhao*hh; end end end % 载物台偏角修正方程 ·½³ for j=1:180 for i=1:2 x(j,i)=round(x(j,i)); y(j,i)=round(y(j,i)); end end zhen=zeros(2000,3000); for i=1:180 for j=1:2 zhen(y(i,j),x(i,j))=1; end end % 修正后矩阵之生成
9 次阅读|0 个评论
分享 A银行业务核心系统使用总结(待更新)
gznshql 2012-8-8 00:03
使用的人性化不足之处 1.输入方面提示不够人性化,比如如何选中,某些地方不能选中复制 2.输入出错时,不能定位到精确位置给出相应提示 3.当进入某个菜单选择项时,一定要选中某项,否则无法退出 4.根据某项内容回显的数据,当修改该项内容时,回显的内容不会联动改变(此处为Bug) 5.在功能界面上缺少一步退出功能界面的选项,必须连续多按几次ESC键,容易出错 6.可以通过增加身份证识别和读取系统,提高业务办理速度。
个人分类: 银行IT|19 次阅读|0 个评论
分享 系统一和系统二的故事
热度 1 zhangjinami 2012-8-3 14:43
系统一是人的潜意识,反应速度比较快,准确度不高,可以培养 系统二是人的思考,反应速度慢,准确度高,可以锻炼
19 次阅读|1 个评论
分享 交通银行首推3D网上银行系统
web3d 2012-5-30 08:22
随着互联网的发展和安全性能的加强,越来越多的公司、企业、工厂开始了网上电子商务。而传统的电子商务已经走向暮年,不再具备吸引力。而在这个处处要求创新的世界里,跟随着发展必然没有空间,这样散失了发展的意义,而如今时代潮流把目光移向了3D技术,也是在这样的时代大潮下,互联网电子商务也谋求了变化,开始向3D商务转变。 在时代的发展下,网上银行是电子商务的依托,虽然他们具备最先天的优势。但是他们还是需要进步,不然等待他们的也只是淘汰。这不如今交通银行首次在网上开展了网上银行3D商务,把整个银行都搬到了网络上,这样只要你登录交通银行的电子3D商务银行,就可以让你享受到和真实生活中银行服务大厅一样的效果。这不得不说是一次大胆的创新,也说明了我国对虚拟现实技术的重视,3D商务普及必将是未来电子商务发展的必然趋势。 3D网上银行系统,网上银行3D商务。 交通银行推出的这款3D网上银行系统,是结合了虚拟现实技术进行现实生活中的交通银行服务系统的模拟进行3D模型制作。然后,在通过3D技术把所有的模型嫁接、植入构成一个完整的系统。上传上网络上。由此,你就可以登录该3D网上银行系统进行和日常生活中一样的,犹如身临其境般地办理查询、转账、理财和贵金属产品购买等多项网上银行业务。 3D网上银行系统,所包含的系统和功能。 目前为止,就交通银行推出的这款网上银行3D商务来看,它主要包含的系统业务区域,基本和现实生活中的交通银行大厅一样。分别有:大厅、新产品体验区、转账交易区、财富管理中心、贵金属交易区五大功能区。而在这五大功能区中,最为热点和技术科技性的是贵金属交易区。在这个贵金属交易区,你可以实现实物黄金贵金属的360度全方位的立体观测。最后,确定是否购买。 除此之外,交通银行的此款网上银行3D商务更具备强大的交互式功能,你可以和周围的其中用户沟通,此外,在你办理业务的时候,你还可以看到在你的区域有4个虚拟客户服务经理,任你挑选,为你完成业务办理。 来源于纳金网
0 个评论

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-25 07:10