《基于图像技术的人工湿地植物生长信息提取方法.pdf》由会员分享,可在线阅读,更多相关《基于图像技术的人工湿地植物生长信息提取方法.pdf(10页珍藏版)》请在专利查询网上搜索。
1、10申请公布号CN102022984A43申请公布日20110420CN102022984ACN102022984A21申请号200910070434822申请日20090915G01B11/24200601G01B11/28200601G06T7/60200601G06F9/4420060171申请人南开大学地址300071天津市南开区卫津路94号72发明人黄岁樑唐婷芳子54发明名称基于图像技术的人工湿地植物生长信息提取方法57摘要一种基于计算机图像技术的人工湿地植物生长及形态变化的信息提取方法。将定期采集的人工湿地植物图像文件数码照片导入程序,显示在计算机屏幕上,程序将屏幕坐标转换成实际坐。
2、标;用户通过人机互动的软件操作获取图像中植物轮廓通过图像边缘识别实现及植株大小通过植株投影面积表征;将同组图像信息相对时间轴表示,可绘制出人工湿地植物生长曲线。本发明的优点在于结合数码成像原理,采用计算机图像技术,在VISUALC环境下设计开发软件,对人工湿地植物生长进行监测,对人工湿地植物不造成任何破坏。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书6页附图2页CN102022998A1/1页21一种基于计算机图像技术的,对采集的人工湿地植物图像进行处理以提取植物生长过程及形态变化信息的方法,其特征在于采用数码相机,定期对人工湿地植物进行图像采集,用户通过。
3、人机互动的软件操作,获取图像中植物轮廓及植株大小植物投影面积表征,以及将同组图像信息相对时间轴表示绘制出人工湿地植物生长曲线。2根据权利要求1所述的人工湿地植物生长信息提取方法,其特征在于,所述对人工湿地植物进行有背景并参照物式图像采集采用了这样的方法选取尺寸远大于植株大小的正红纯色矩形板作背景,矩形板左上角位置水平钉上150CM20CM或根据需要纯白色参照矩形条,右下角位置垂直钉上100CM20CM或根据需要纯白色参照矩形条,图像采集时将背景板尽量紧贴近植物,并保持背景板底边平行于地面,竖边与地面垂直。3根据权利要求1所述的人工湿地植物生长信息提取方法,其特征在于,所述植物轮廓提取采用了这样。
4、的方法自主识别与人机互动相结合。建立的CDIB位图操作类库中包含各种算法的图像边缘检测函数,用户可选取最合适的函数命令计算机进行图像边缘自主识别;所得图像轮廓图包含多条轮廓线及一些不相关信息,用户在线提取模式下,通过鼠标点击所需提取的轮廓线中任意一点,软件从这一点出发进行迭代搜索将植物轮廓线一次提取出来。4根据权利要求1所述的人工湿地植物生长信息提取方法,其特征在于,所述植株大小定量采用了这样的方法通过软件计算将图像的RGB色彩属性模式转化为HSV色彩属性模式,测试若干像片后选取HSV颜色模式下H分量色相在120附近即区间75,165中的点为绿色,二重循环扫描图像获得图像中绿颜色像素点集,即为。
5、植株投影面积像素集;结合坐标转换因子将像素块面积换算为实际投影面积。5、根据权利要求4所述的人工湿地植物生长信息提取方法,其特征在于,所述坐标转换采用了这样的方法坐标转换模式下,鼠标分别点击水平参照矩形条/垂直参照矩形条的两端点,在对话框中输入实际长度15CM/实际长度10CM,得到电脑屏幕上图像显示的水平标度/垂直标度,即得到像素点间距离的换算因子,像素块面积换算因子即为水平标度与垂直标度的乘积。6根据从2到5中任一权利要求所述的人工湿地植物生长信息提取方法,其特征在于基于计算机图像技术,所述图形图像文件为BMP位图文件。权利要求书CN102022984ACN102022998A1/6页3基。
6、于图像技术的人工湿地植物生长信息提取方法技术领域0001本发明涉及的基于计算机图像技术的人工湿地植物生长及形态变化信息提取方法及程序,实现对人工湿地植物的“无损伤式”监测。背景技术0002人工湿地是20世纪50年代产生,7080年代发展起来的一种抗冲击负荷能力强、系统配置可塑性强、投资及运行费用低、生态效益显著的新型污水处理技术。湿地植物是人工湿地床体重要组成部分之一;对于人工湿地植物生长的研究,传统方法仅停留在污水处理前后将整株的植物挖出称干重,或者将其截成根、茎和叶三段测污染物富集量。这种方法不能监测植物生长的过程。0003本发明的目的在于提供这样一种人工湿地植物生长过程监测的方法通过数码。
7、相机成像,定期对人工湿地植物进行图像采集;将图像文件导入程序,显示在计算机屏幕上,用户通过人机互动的软件操作获取图像中植物轮廓通过图像边缘识别实现及植株大小通过植株投影面积表征;将同组图像信息相对时间轴表示,可绘制出人工湿地植物生长曲线。发明内容0004本发明提出的目的通过以下方案来实现1图像采集。采用数码相机,针对所研究的人工湿地植物进行有背景并参照物式图像采集,以避开背景杂色对信息提取的影响,并可尽量减免照相过程中由于拍摄远近、焦距调节等因素带来的成像缩小比例不一致问题。2图像导入及轮廓提取。用操作DIB的WIN32API函数在计算机屏幕上显示带有数据信息的图像位图文件;计算机通过图像边缘。
8、检测函数自主识别得到图像轮廓图,用户通过鼠标点击植物轮廓线中的任意一点可将整个植物轮廓从错综复杂的轮廓线中提取出来,并将信息列表转化成TXT文件输出。3植株大小定量。利用照片中参照物对图像进行坐标转换,得到像素点间长度换算因子及像素块面积换算因子;植株大小通过植株投影面积表征,植株投影面积通过图像中绿色像素数统计并结合像素块面积换算因子换算而得。0005所述对人工湿地植物进行有背景并参照物式图像采集采用了这样的方法选取尺寸远大于植株大小的正红纯色矩形板作背景,矩形板左上角位置水平钉上150CM20CM纯白色参照矩形条,右下角位置垂直钉上100CM20CM纯白色参照矩形条,图像采集时将背景板尽量。
9、紧贴近植物,并保持背景板底边平行于地面,竖边与地面垂直。0006所述植物轮廓提取采用了这样的方法自主识别与人机互动相结合。建立的CDIB位图操作类库中包含各种算法的图像边缘检测函数,用户可选取最合适的函数命令计算机进行图像边缘自主识别;所得图像轮廓图包含多条轮廓线及一些不相关信息,用户可以在线提取模式下,通过鼠标点击所需提取的轮廓线中任意一点,软件从这一点出说明书CN102022984ACN102022998A2/6页4发进行迭代搜索将植物轮廓线一次提取出来。0007所述坐标转换采用了这样的方法坐标转换模式下,鼠标分别点击水平参照矩形条/垂直参照矩形条的两端点,在对话框中输入实际长度15CM/。
10、实际长度10CM,得到电脑屏幕上图像显示的水平标度/垂直标度,即得到像素点间距离的换算因子,像素块面积换算因子即为水平标度与垂直标度的乘积。0008所述植株大小定量采用了这样的方法通过软件计算将图像的RGB色彩属性模式转化为HSV色彩属性模式,测试若干像片后选取HSV颜色模式下H分量色相在120附近即区间75,165中的点为绿色,二重循环扫描图像获得图像中绿颜色像素点集,即为植株投影面积像素集;结合坐标转换因子将像素块面积换算为实际投影面积。0009与现有的方法相比,本发明的优点在于0010本发明针对在科研及生产实践中,人工湿地植物生长的研究需要将植物连根拔出或收割地上部分从而破坏人工湿地结构。
11、这一实际情况,结合数码相机成像,采用计算机图像技术,在VISUALC环境下设计开发软件,使研究者从过去“要么破坏人工湿地,要么只研究起点和终点的植物状况”中解放出来,解决该问题。附图说明0011下面结合附图和并非特定的实施例来对本发明作更详尽的说明。0012图1是本发明所述方法的流程图。0013图2是HSV色彩属性图中本发明所取绿色区间所处的位置示意图。0014图3是对五组用本发明所述方法采集的人工湿地植物图像进行信息提取所得的植物生长曲线图。具体实施方式0015图像采集与图像处理是本发明的两个要点,下面将以按上述有背景参照物式方法采集的人工湿地植物图像转化为BMP位图文件为例来对本发明做出更。
12、详细的说明。0016在图像采集部分中,按上述发明内容中所述做好背景及参照物,采用普通数码相机,选取一定的拍摄距离,正对植株拍摄即可。由于软件的普适性,对图像的尺寸大小没有严格要求,将JEPG格式的图像文件正确导入计算机中转化为24位位图格式便可读入软件进行处理。0017WINDOWS操作系统下有两大类位图格式GDIGRAPHICSDEVICEINTERFACE,图形设备接口位图和DIBDEVICEINDEPENDENTBITMAPS,设备无关位图。由于DIB在程序设计方面提供许多超越GDI位图的优势DIB本身携带有自己的颜色信息,使颜色调配管理变得更加简单;任何运行WINDOWS操作系统的计算。
13、机均可处理存储于BMP文件或可执行文件内部资源中的DIB;尽管TIFF、GIF或JPEG等格式的图像文件在转换后也可以被使用,但却只有采用DIB格式的文件被WIN32API所直接支持;虽然在WIN32中允许在不同进程之间通过传递GDI位图句柄到剪贴板进行位图的传递,但在后台所进行的操作却是借助DIB来完成的将设备相关的GDI位图转换成一个DIB,并将此DIB复制说明书CN102022984ACN102022998A3/6页5到一段共内存,因此,本软件图像操作部分的程序实现代码均以DIB作为处理对象。0018在图像处理部分中,将注意力集中在图像处理的算法实现上。封装WINDOWS系统提供的一些重。
14、要的DIB访问函数,建立一个实现DIB主要操作功能的类CDIB,在进行图像处理编程时,只需通过CDIB类提供的公有型PUBLIC接口函数调用,即可实现与之相关的功能。由于图像边缘识别功能块有多种算法,将这些函数独立出来封装成CBIANYUANJIANCEDIB类。0019CDIB类中主要的功能函数有五个LOADFILE、SAVEFILE、SAVELINE、HSV、GETAREA,分别用于读取位图、保存位图、提取植物边缘信息、RGB模式转换为HSV模式以及统计植物像素面积。CBIANYUANJIANCEDIB类中包含的边缘识别函数有GRADIENT、LAPLACIAN、GUASSLAPLACIA。
15、N、KRISCH、SOBEL、ROBERT及PREWITT,分别对应边缘检测的梯度算子、拉普拉斯算子、高斯拉普拉斯算子、克瑞斯算子、卷积算子、罗伯特算子和普瑞维特算子。0020到目前为止,MICROSOFT一直没有将有关DIB的API函数封装进MFC,本软件封装到CDIB类中的位图读写函数VOIDLOADFILECONSTCHARDIBFILENAME和VOIDSAVEFILECONSTCSTRINGFILENAME以及封装至CBIANYUANJIANCEDIB类中的图像边缘检测函数,主要参阅参考文献VC图像处理程序设计杨淑莹编著,清华大学出版社,2005年1月第2版,ISBN97878108。
16、24507提供的相关程序代码。0021SAVELINE函数负责从已进行边缘识别的图像中,将用户指定的植物轮廓提取出来,并将所提取的轮廓信息以列表形式写入TXT文件中。函数原型在DIBH中声明如下0022INTSAVELINECPOINTPOINT;0023软件在轮廓提取模式BOOL型变量BCURVE为TRUE下,用户通过鼠标点击所需提取的轮廓,该鼠标单击事件驱动CPINFOEVIEW类中的ONLBUTTONDOWNUINTNFLAGS,CPOINTPOINT函数运行,调用SAVELINE函数时ONLBUTTONDOWN函数中实参POINT传递给SAVELINE函数的形参POINT,进行定点指定。
17、的轮廓一次提取。函数返回值为INT型,为0表示保存曲线成功;大于0时表示取点失败,用户重新取点;小于0表示取点失败,用户不重新取点。取点成功与否主要取决于用户在鼠标点击时,是否正确点在了轮廓线上,若偏离轮廓线在函数允许范围内3像素偏移量,程序会自动纠正找到最近的轮廓点并开始一次提取,若偏离超出允许范围,则程序提示需重新取点。定点轮廓一次提取功能实现函数程序及注释如下说明书CN102022984ACN102022998A4/6页600240025说明书CN102022984ACN102022998A5/6页70026上述代码中斜体标出轮廓上所有点一次提取使用的是递归搜索法,通过调用LINEPOI。
18、NT函数实现,程序及注释如下00270028HSV函数负责将24位位图中的颜色信息由RBG模式转化为HSV模式,返回FLOAT型的像素点H分量。函数原型在DIBH中声明如下0029FLOATHSVCONSTINTR,CONSTINTG,CONSTINTB;0030其中R、G、B分别为RGB色彩模式下像素点的红、绿、蓝三原色分量,先将说明书CN102022984ACN102022998A6/6页8其从0255的数值转化到0010区间内,如果三个因数中最大的值为MAX,最小的值为MIN,对应像素点的H、S、V值的计算方程组如下。00310032计算所得的色相H值范围为0360如果H大于360,用3。
19、60除H值的余数作为色相值,饱和度S和明度V值范围为0010,再将其转化为0100,即为HSV模式下的分量值。0033若MAXMIN,H0,即没有色彩,是纯灰色。0034若MAX0,S0,也是没有色彩。0035若V0,即没有亮度,是纯黑色。0036GETAREA函数负责统计植物像素面积。图像中绿颜色像素点的获取主要基于HSV颜色模式下,H分量在120附近的点取区间75,165为各种绿色这一原理,故只需用二重循环对整个位图矩阵进行扫描,调用HSV获取H分量进行筛选,记录下所有符合要求的点即可。用户可通过菜单栏中的坐标转换及对话框提示完成实际坐标与屏幕坐标之间的换算,进而得到实际投影面积。0037。
20、在点击参考矩形条时,用户点击偏移量若超出3像素的误差范围,程序会以提醒用户重新点击参考点以确保计算数据的精度。0038本发明所述的基于图像技术的人工湿地植物生长信息提取方法,不仅仅限于说明书和实施方式中所列运用,它完全可以适用于各种适合本发明之领域,对于熟悉本领域的人员而言可容易地实现另外的优点和进行修改,因此在不背离权利要求及等同范围所限定的一般概念的精神和范围的情况下,本发明并不限于特定的细节、代表性的设备和这里示出与描述的图示示例。说明书CN102022984ACN102022998A1/2页9图1图2说明书附图CN102022984ACN102022998A2/2页10图3说明书附图CN102022984A。