莖叶图和箱线图在图像处理中的应用.docVIP

茎叶图和箱线图在图像处理中的应用 1 在MATLAB中获取图像数据 打开MATLAB6.1，在command window中输入命令，以获取图像中的灰度数据。本文以matlab中自带图片blood1.tif图片为示例图片。图片原图如图1所示。 图1 原图 在command window中输入如下命令： i=imread(blood1.tif) （注意此处不加分号） 随后在command window中会出现一个数字矩阵，如图2所示 图2 数字矩阵 将这个矩阵全部复制下来，留着数据备用。 2 在SAS中绘制图像的茎叶图和箱线图 打开sas8.2软件，在Editor中输入以下程序： title image process : Steam-leaf and Box plot; data imagedata; input gray; cards; 101 113 126 118 148 193 213 206 193 176 165 160 155 148 148 155 96 106 113 118 138 196 218 206 193 171 160 155 148 148 148 148 56 80 85 90 118 171 196 201 193 176 165 155 148 148 148 148 46 48 55 63 101 123 138 181 188 181 171 155 148 148 148 148 46 46 46 48 80 96 113 126 176 181 176 165 155 148 148 148 46 46 46 46 48 68 90 113 131 176 188 176 165 160 155 155 46 46 46 46 46 46 68 96 118 160 181 188 171 165 165 165 …………………………………………………………….. …………………………………………………………….. (将在matlab中获取的图像灰度数据复制于此处) ; run; proc univariate plot ;run; 此时在sas系统中output窗口中就会有数据出现，如图3所示： 图3 图3给出了图像数据的偏度（skewness）、峰度（kurtosis）、级差（range）等量化的数据。 在图4中给出了箱线图，由于数据量在观测区间上超过48，故给出了一个水平直方图。 图 4 水平直方图和箱线图 此时选择sas菜单的soluations－analysis-interactive data analysis 会出现如下对话框（如图5所示）： 图 5 对话框 在work工作区中选择imagedata，这是我们刚才在sas程序中当前工作区的数据，然后单击open，随即出现下一个数据窗口，如图6所示： 图6 数据表格 在表格中选中gray ，使其变黑，之后选择菜单analyze－Box plot/Mosaic plot 便得到箱线图，如图7所示： （a） （b） 图7 箱线图 其中，a、b是选中不同部分时数据的显示。 蚅肂蒄薅羃膁膃螀衿膀芆薃螅腿莈螈蚁膈薀薁肀膇芀蒄羆膆莂虿袂膆蒅蒂螈膅膄蚈蚄芄芇蒁羂芃荿蚆袈节蒁葿螄芁芁蚄螀芀莃薇聿芀蒅螃羅艿薈薅袁芈芇螁螇袄莀薄蚃羄蒂蝿羂羃膂薂羈羂莄螈袄羁蒆蚀螀羀蕿蒃肈罿芈虿羄羈莁蒁袀肈蒃蚇螆肇膃蒀蚂肆芅蚅肁肅蒇蒈羇肄蕿螃袃肃艿薆蝿膃蒀蚂肆芅蚅肁肅蒇蒈羇肄蕿螃袃肃艿薆蝿艿薆蝿膃蒀蚂肆芅蚅肁肅蒇蒈羇肄蕿螃袃肃艿薆蝿肅蒇蒈羇肄蕿螃袃肃艿薆蝿肃艿薆蝿