CN103954223B - 基于双摄像头毫米级生物3d打印机出丝宽度测量方法 - Google Patents
基于双摄像头毫米级生物3d打印机出丝宽度测量方法 Download PDFInfo
- Publication number
- CN103954223B CN103954223B CN201410185499.8A CN201410185499A CN103954223B CN 103954223 B CN103954223 B CN 103954223B CN 201410185499 A CN201410185499 A CN 201410185499A CN 103954223 B CN103954223 B CN 103954223B
- Authority
- CN
- China
- Prior art keywords
- image
- printer
- represent
- edge
- axis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于双摄像头毫米级生物3D打印机出丝宽度测量方法。在完成机械安装之后,用户启动生物3D打印机,测量系统利用双摄像头实时采集图像,再进行深度图像处理,经过滤波、锐化及边缘提取、二值化和计算后能够得到较为准确的打印出丝宽度。两路摄像头处理性能优势明显,不仅克服了单摄像头图像采集存在图像视觉死角的缺点,而且在一定程度上补偿了光源造成的视觉误差,实时计算获得出丝的宽度,一改原有生物3D打印设备缺少反馈的缺点。且该发明专利硬件可移植性强,能够作为辅助工具应用于绝大部分3D打印机。产品具有广阔的前景与市场。
Description
技术领域
本发明属于生物3D打印快速成形领域,具体涉及一种基于双摄像头毫米级生物3D打印机出丝宽度测量方法。
背景技术
3D打印技术逐渐进入各个下游应用领域,生物3D打印是最前沿的研究领域。作为一项前沿制造技术,“3D打印”已经逐步应用于航天军工,模具制造,动漫制作、文化创意等多个领域,随着技术的发展,它的应用领域还在不断扩展。生物3D打印技术是使用3D打印的方法成型生物材料,特别是细胞材料,用来制造人工的组织、器官等一系列生物医疗领域的产品。它是3D打印研究中最前沿的领域。现阶段,生物3D打印的应用主要包括:细胞打印、组织工程支架和植入物打印。
组织工程支架和植入物的生物3D打印逐渐成熟,骨组织的打印市场空间巨大。生物3D打印借助CT、ECT技术获取人体模型器官模型,其技术牵涉到3D模型的建立,有限元的分析等。最后指导3D打印设备喷射生物相容性材料,形成所需要的结构。3D打印在构建植入物的微观结构方面相对传统工艺有很明显的优势。在美国,仅骨移植修复材料的市场空间就达200亿美金。此外,一些血管支架等领域的应用也在逐渐发展。
生物3D打印现在桌面级别的通常是采用热熔堆积技术,即将塑料先融化然后挤压,经过一个零点几毫米的喷口挤出,按照电脑设计好的轨迹运动,最后成型。专业的打印机,能够保证吐丝匀称性,但造型时间长。高速生物3D打印机,在提高打印速度的同时损失了打印质量,尤其对于拐角的打印,由于机械结构的原因,出丝宽度大于正常值。
保证打印质量是增材制造的核心,而打印质量和硬件软件都有着密切的关系。现市场上同类仪器产品,均采用了开环成型控制系统,这就意味着,但即使精度再高的生物3D打印机,打印时都是无法实时控制出丝的宽度,更谈不上通过测量实际丝的宽度来修正这一系统误差。
目前市场上同类产品针对这一无法避免的系统误差,一般采取如下手段进行修正克服:
1.大数据修正,即将获各组数据与实际测量宽度进行比较,推导生成经验修正公式加以修正。修正公式基于大量的测量数据之上,工程量巨大,且该修正公式普适性较差。
2.模型简化法,即将复杂的行程转换为简单的行程,降低打印难度。一般需要打印质量精细点的通常设置层高为0.3mm,如果需粗糙的可以设置成0.4mm,如此以来,其在提高打印出丝宽度的同时损失了物品的逼真度。
发明内容
本发明针对现有技术的不足,提出了一种基于双摄像头毫米级生物3D打印机出丝宽度测量方法,利用两个彩色高清摄像头OV7670,进行对出丝过程进行实时拍摄,同时获得实时出丝宽度数据。
本发明方法具体是:
步骤一:以打印机喷头工作原点建立xyz三维坐标轴,即在初始状态下,喷头中心与z轴重合,x轴与工作台平面短边共线,y轴与工作台平面长边共线;两个摄像头将安装在成型室工作台导轨上,第一摄像头在zoy平面,其拍摄方向对准喷头,与xoy平面平行的平面保持θ1的工作角度,25°≤θ1≤35°;第二摄像头在zox平面,其拍摄方向对准喷头,与xoy平面平行的平面保持θ2的工作角度,25°≤θ2≤35°;喷头在x轴的运动分量在第一摄像头中能观测到;喷头在y轴的运动分量在第二摄像头中能观测到;成型工作台面在初识状态下距离成型室底面500mm,其后在成型过程中,随着每一层成型结束,成型工作台面将下降固定高度,但其在成型过程中始终和xoy平面保持平行;
步骤二:打开系统电源按钮,初始化各个模块;
步骤三:调节两个彩色高清摄像头的焦距,调节至画面最佳状态,将实时拍摄图像显示在3.5寸TFT液晶屏上;
步骤四:启动生物3D打印机,喷头在伺服电机的带动下开始运动成型,此时喷头开始出丝;由于喷头在工作台面的轨迹运动,在液晶屏中将会产生相对的移动;
步骤五:将图像用傅立叶变换、沃尔什变换进行预处理;
步骤六:将步骤五处理后的图像再进行深度图像处理,处理步骤包括滤波、锐化及边缘提取、二值化和计算;
滤波:拍摄的图像边缘具有一定的模糊性,在边缘提取之前,对图像进行中值滤波的处理,较好地去除图像中的噪声并避免了图像模糊的加剧;
锐化及边缘提取:经过滤波处理后的图像,采用基于梯度的边缘检测法,其一阶微分算子采用Roberts交叉微分算子,二阶微分算子采用Laplacian算子和Wallis算子;
二值化:在对图进行完边缘提取后,在基于整幅图像的灰度值统计平均的基础上选择一个阈值Th来对提取的边缘图片进行二值化处理,将图像边缘赋值为黑色并进行存储;
计算:将得到的照片数据传入控制器进行处理;生物3D打印机在造型时,会打印正方形的支撑支架;在分析过程中将牵涉到四个量:第一个量,待测丝直径图像上的像素距离:其中xa表示图像上待测丝边缘上点A的x轴坐标,ya表示图像上待测丝边缘上点A的y轴坐标;同理,xb表示图像上待测丝边缘上点B的x轴坐标,yb表示图像上待测丝边缘上点B的y轴坐标,该量通过图像分析获得;第二个量,正方形支架边长图像像素距离,其中xc表示A、B点所在的正方形支架的一个顶点C的x轴坐标,yc表示正方形支架顶点C的y轴坐标;同理,xd表示点C相邻的正方形支架顶点D的x轴坐标,yd表示顶点D的y轴坐标;第三个量,正方形支架边长实际距离RCD,该量是受上位机软件控制,故其准确值是已知的;第四个量,即我们所测量出丝的宽度RAB;根据比例关系:
本发明的有益效果在于:两路摄像头处理性能优势明显,不仅克服了单摄像头图像采集有死角的缺点,而且在一定程度上补偿了光源造成的视觉误差,能够实时计算获得出丝的宽度,一改原有生物3D打印设备缺少反馈的缺点。且该发明专利硬件可移植性强,能够作为辅助工具应用于绝大部分3D打印机。产品具有广阔的前景与市场。
附图说明
图1是本发明使用流程图;
图2是摄像头、喷头物理位置关系与安装图;
图3是经过图像滤波后的图像;
图4是经过多次滤波后用于计算的图像。
具体实施方式
本专利着重偏重于基于双摄像头毫米级生物3D打印机出丝宽度测量方法的硬件实现部分和简单的软件算法,详细的软件提取算法将在另一专利中做详细说明阐述。下面结合附图对本发明作进一步描述。如图1所示,本发明方法具体步骤是:
步骤一:以打印机喷头工作原点建立xyz三维坐标轴,即在初始状态下,喷头3中心与z轴重合,x轴与成型室门垂直,y轴与成型室门平行。两个摄像头将安装在成型室工作台导轨5上,第一摄像头1在zoy平面,其拍摄方向对准喷头,且与经过喷头、与xoy平面平行的平面保持θ1的工作角度,25°≤θ1≤35°;第二摄像头2在zox平面,其拍摄方向对准喷头,且与经过喷头、与xoy平面平行的平面保持θ2的工作角度,25°≤θ2≤35°;喷头在x轴的运动分量在第一摄像头1中能观测到;喷头在y轴的运动分量在第二摄像头2中能观测到;成型工作台4面在初识状态下距离成型室底面500mm,其后在成型过程中,随着每一层成型结束,成型工作台面将下降固定高度,但其在成型过程中始终和xoy平面保持平行。该步骤为机械安装步骤,在今后的设备使用中可以不需要重复该步骤。具体安装示意图见图2。
步骤二:打开系统电源按钮,初始化各个模块,将实时拍摄图像显示在3.5寸TFT液晶屏上。
步骤三:调节两个彩色高清摄像头的焦距,调节至画面最佳状态。
步骤四:启动生物3D打印机,喷头在伺服电机的带动下开始运动成型,此时喷头开始出丝。由于喷头在工作台面的轨迹运动,在液晶屏中将会产生相对的移动。
步骤五:将图像用傅立叶变换、沃尔什变换进行预处理。
步骤六:将步骤五处理后的图像再进行深度图像处理,处理步骤包括滤波、锐化及边缘提取、二值化和计算。
滤波:拍摄的图像边缘具有一定的模糊性,在边缘提取之前,对图像进行中值滤波的处理,较好地去除图像中的噪声并避免了图像模糊的加剧,如图3所示。
锐化及边缘提取:经过滤波处理后的图像,采用基于梯度的边缘检测法,其一阶微分算子采用Roberts交叉微分算子,二阶微分算子采用Laplacian算子和Wallis算子。
二值化:在对图进行完边缘提取后,在基于整幅图像的灰度值统计平均的基础上选择一个阈值Th来对提取的边缘图片进行二值化处理,将图像边缘赋值为黑色并进行存储;
锐化及边缘提取:经过滤波处理后的图像,存在一定程度的模糊,需要对图像进行锐化处理,可以采用基于梯度的边缘检测法,其一阶微分算子可采用Roberts交叉微分算子、Sobel算子和Priwitt微分算子,二阶微分算子可采用Laplacian算子、Wallis算子。锐化处理主要是加强图像的边缘,便于将图像从背景中准确地提取出来。
计算:将得到的照片数据传入控制器进行处理。生物3D打印机在造型时,会打印正方形的支撑支架。在分析过程中将牵涉到四个量:第一个量,待测丝宽度图像像素距离:其中xa表示图像上待测丝边缘上点A的x轴坐标,ya表示图像上待测丝边缘上点A的y轴坐标;同理,xb表示图像上待测丝边缘上点B的x轴坐标,yb表示图像上待测丝边缘上点B的y轴坐标,该量通过图像分析获得。第二个量,正方形支架边长图像像素距离,其中xc表示图像上正方形支架顶点C的x轴坐标,yc表示图像上正方形支架顶点C的y轴坐标;同理,xd表示图像上与点C相邻的正方形支架顶点D的x轴坐标,yd表示图像上与点C相邻的正方形支架顶点D的y轴坐标,该量通过图像分析获得。
第三个量,正方形支架边长实际距离RCD,该量是受上位机软件控制,故其准确值是已知的。第四个量,即我们所测量出丝的宽度RAB。结合处理后的图像图4,根据计算公式,可得:
待测丝宽度图像像素距离:
正方形支架边长图像像素距离:
支架间隙实际距离RCD可通过上位机修改,一般取2.5mm。
根据比例公式:
计算得RAB=0.88mm。经过实际测试,该方法的测量精度误差在10%。
Claims (1)
1.基于双摄像头毫米级生物3D打印机出丝宽度测量方法,其特征在于,该方法具体包括以下步骤:
步骤一:以打印机喷头工作原点建立xyz三维坐标轴,即在初始状态下,喷头中心轴线与z轴重合,x轴与成型室工作台短边共线,y轴与成型室工作台长边共线;两个彩色高清摄像头将安装在成型室工作台导轨上,第一彩色高清摄像头在zoy平面,其拍摄方向对准喷头,与xoy平面平行的平面保持θ1的工作角度,25°≤θ1≤35°;第二彩色高清摄像头在zox平面,其拍摄方向对准喷头,与xoy平面平行的平面保持θ2的工作角度,25°≤θ2≤35°;喷头在x轴的运动分量在第一彩色高清摄像头中能观测到;喷头在y轴的运动分量在第二彩色高清摄像头中能观测到;成型室工作台在初始状态下距离成型室底面500mm,其后在成型过程中,随着每一层成型结束,成型室工作台将下降固定高度,但其在成型过程中始终和xoy平面保持平行;
步骤二:打开系统电源按钮,初始化各个模块;
步骤三:调节两个彩色高清摄像头的焦距,调节至画面最佳状态,将实时拍摄图像显示在3.5寸TFT液晶屏上;
步骤四:启动生物3D打印机,喷头在伺服电机的带动下开始运动成型,此时喷头开始出丝;由于喷头在成型室工作台的轨迹运动,在液晶屏中将会产生相对的移动;
步骤五:将图像用傅立叶变换、沃尔什变换进行预处理;
步骤六:将步骤五处理后的图像再进行深度图像处理,处理步骤包括滤波、锐化及边缘提取、二值化和计算;
滤波:拍摄的图像边缘具有一定的模糊性,在边缘提取之前,对图像进行中值滤波的处理,较好地去除图像中的噪声并避免了图像模糊的加剧;
锐化及边缘提取:经过滤波处理后的图像,采用基于梯度的边缘检测法,其一阶微分算子采用Roberts交叉微分算子,二阶微分算子采用Laplacian算子和Wallis算子;
二值化:在对图进行完边缘提取后,在基于整幅图像的灰度值统计平均的基础上选择一个阈值Th来对提取的边缘图片进行二值化处理,将图像边缘赋值为黑色并进行存储;
计算:将得到的照片数据传入控制器进行处理;生物3D打印机在造型时,会打印正方形的支撑支架;在分析过程中将牵涉到四个量:第一个量,待测丝直径图像上的像素距离:其中xa表示图像上待测丝边缘上点A的x轴坐标,ya表示图像上待测丝边缘上点A的y轴坐标;同理,xb表示图像上待测丝边缘上点B的x轴坐标,yb表示图像上待测丝边缘上点B的y轴坐标,该量通过图像分析获得;第二个量,正方形支架边长图像像素距离, 其中xc表示A、B点所在的正方形支架的一个顶点C的x轴坐标,yc表示正方形支架顶点C的y轴坐标;同理,xd表示点C相邻的正方形支架顶点D的x轴坐标,yd表示顶点D的y轴坐标;第三个量,正方形支架边长实际距离RCD,该量是受上位机软件控制,故其准确值是已知的;第四个量,即我们所测量出丝的宽度RAB;根据比例关系:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410185499.8A CN103954223B (zh) | 2014-05-04 | 2014-05-04 | 基于双摄像头毫米级生物3d打印机出丝宽度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410185499.8A CN103954223B (zh) | 2014-05-04 | 2014-05-04 | 基于双摄像头毫米级生物3d打印机出丝宽度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103954223A CN103954223A (zh) | 2014-07-30 |
CN103954223B true CN103954223B (zh) | 2017-01-11 |
Family
ID=51331533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410185499.8A Active CN103954223B (zh) | 2014-05-04 | 2014-05-04 | 基于双摄像头毫米级生物3d打印机出丝宽度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103954223B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105338307B (zh) * | 2015-08-24 | 2019-02-15 | 西安电子科技大学 | 一种3d打印远程监控装置及其方法 |
US10766070B2 (en) * | 2015-08-28 | 2020-09-08 | Materialise N.V. | Self supporting in additive manufacturing |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303226A (zh) * | 2008-06-25 | 2008-11-12 | 南昌航空大学 | 一种基于最大连通域的线路板线宽测量方法 |
CN103213280A (zh) * | 2013-03-29 | 2013-07-24 | 杭州电子科技大学 | 基于线宽测量的智能快速成形流量控制方法 |
CN103613272A (zh) * | 2013-11-12 | 2014-03-05 | 杭州电子科技大学 | 一种玻璃光纤拉丝设备优化生产方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6276402A (ja) * | 1985-09-30 | 1987-04-08 | Kobe Steel Ltd | 板幅測定方法 |
JP2007240503A (ja) * | 2005-04-08 | 2007-09-20 | Hitachi Kokusai Electric Inc | 顕微鏡撮像装置及び寸法測定装置 |
-
2014
- 2014-05-04 CN CN201410185499.8A patent/CN103954223B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303226A (zh) * | 2008-06-25 | 2008-11-12 | 南昌航空大学 | 一种基于最大连通域的线路板线宽测量方法 |
CN103213280A (zh) * | 2013-03-29 | 2013-07-24 | 杭州电子科技大学 | 基于线宽测量的智能快速成形流量控制方法 |
CN103613272A (zh) * | 2013-11-12 | 2014-03-05 | 杭州电子科技大学 | 一种玻璃光纤拉丝设备优化生产方法 |
Non-Patent Citations (3)
Title |
---|
《组织工程与人体器官快速成形制造》;岳秀艳等;《华北煤炭医学院学报》;20010731;第3卷(第4期);第405-406页 * |
《面向RP的CT图像处理技术综述》;姜睿智等;《CT理论与应用研究》;20080331;第17卷(第1期);第14-20页 * |
《高分子生物材料快速成形构建组织工程支架的研究》;索海瑞等;《中国组织工程研究与临床康复》;20080205;第12卷(第6期);第1101-1105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103954223A (zh) | 2014-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107672180B (zh) | 一种基于逆向工程技术的3d打印精度检测方法 | |
CN104881874B (zh) | 基于二元四次多项式畸变误差补偿的双远心镜头标定方法 | |
CN104677305B (zh) | 一种基于十字结构光的物体表面三维重建方法和系统 | |
CN102136142B (zh) | 基于自适应三角形网格的非刚性医学图像配准方法 | |
CN101571950B (zh) | 基于各向同性扩散和稀疏表示的图像修复方法 | |
CN104732900B (zh) | 像素缺陷检测方法和装置 | |
CN103134900B (zh) | 一种高超声速飞行器的热防护结构的烧蚀测试方法 | |
CN103471910A (zh) | 一种基于随机点跟踪的金属材料断裂伸长率智能测试方法 | |
CN103247058B (zh) | 一种基于误差分布式多层网格的快速光流场计算方法 | |
CN101984916A (zh) | 基于数字图像处理技术的血管管径测量方法 | |
CN104089586A (zh) | 发动机曲轴轴颈形状误差的图像检测装置和图像检测方法 | |
CN108444413A (zh) | 陶瓷墙地砖平整度检测装置及方法 | |
CN103954223B (zh) | 基于双摄像头毫米级生物3d打印机出丝宽度测量方法 | |
DE102012215215A1 (de) | Verfahren und Vorrichtung zum Erkennen von Abweichungen einer Oberfläche eines Objekts | |
CN106530255A (zh) | 一种汽泡群数据测量方法 | |
CN105227802A (zh) | 一种面向最佳质量图像快速扫描的自适应控制方法及系统 | |
CN115760837A (zh) | 一种基于深度神经网络的晶体质量评估方法及系统 | |
CN109145913A (zh) | 一种非均匀背景离散相灰度图像二值化处理方法 | |
EP4385712A1 (en) | 3d printing control method and apparatus, and electronic device | |
CN107525768A (zh) | 一种dna倍体分析设备的质量控制方法 | |
CN107219619A (zh) | 三维折射率层析显微成像系统及其方法 | |
CN109272575A (zh) | 一种提高数字切片扫描仪建模速度的方法 | |
CN115760893A (zh) | 一种基于核相关滤波算法的单雾滴粒径和速度测量方法 | |
CN106018198B (zh) | 一种气泡粒径的反演计算方法 | |
CN114659463A (zh) | 一种植物表型采集装置及其采集方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |