| 题目要求 |
| 程序代码 |
| 结果图片 |
题目要求：
生成一个较小方差得随机图像（用一个随机变量，使其大部分数值的差异不超过3，其中大部分数值接近0），将此图像载入一个绘制程序，在图像上画一些辐射状相交于一点的线条。
a、对图像使用双边滤波。
程序代码：
1 // OpenCVExerciseTesting.cpp : 定义控制台应用程序的入口点。
2 //
3 //D:\\Work\\Work_Programming\\Source\\Image\\lena.jpg
4
5
6 #include "stdafx.h"
7 #include <cv.h>
8 #include <highgui.h>
9 #include <iostream>
10 #include <math.h>
11 #include <stdlib.h>
12 using namespace cv;
13 using namespace std;
14 //函数声明-->--->-->--->-->--->-->--->//
15
16 double gaussrand(double E = 0.0, double V = 1.0);
17 void InitialzieImageToRandomValue(IplImage * img);
18
19 //<--<--<--<--<--<--<--<--<--函数声明//
20
21 int _tmain(int argc, _TCHAR* argv[])
22 {
23
24 const char * fileName1 = "D:\\Work\\Work_Programming\\Source\\Image\\OpenCVExerciseResult\\Source_5-8_2.jpg";
25 IplImage * src1 = cvLoadImage(fileName1, CV_LOAD_IMAGE_GRAYSCALE);
26 assert(src1);
27
28 //IplImage * img = cvCreateImage(cvSize(500, 500), IPL_DEPTH_8U, 1);
29 //assert(img);
30 //cvZero(img);
31 //InitialzieImageToRandomValue(img);
32
33 cvNamedWindow("原始图像", 0);
34 cvNamedWindow("题目_a", 0);
35 //cvNamedWindow("题目_b", 0);
36 //cvNamedWindow("题目_c", 0);
37 ////cvNamedWindow("题目_d", 0);
38
39 cvShowImage("原始图像", src1);
40 /*char * imageName = "D:\\Work\\Work_Programming\\Source\\Image\\OpenCVExerciseResult\\Result_5-8.jpg";
41 cvSaveImage(imageName, img);*/
42
43 //---------------------------a：开始--------------------------------//
44
45 IplImage * image_BIL = cvCloneImage(src1);
46 IplImage * temp_image = cvCloneImage(src1);
47 cvZero(temp_image);
48 cvZero(image_BIL);
49
50 cvSmooth(src1, temp_image, CV_BILATERAL, 3, 3, 11, 11);
51
52 int j = 0;
53 while (j < 100)
54 {
55 //注意:使用CV_BILATERAL平滑类型时cvSmooth参数中的最后两个11
56 cvSmooth(temp_image, image_BIL, CV_BILATERAL, 7, 7, 11, 11);
57 //cvSmooth(temp_image, image_BIL, CV_BILATERAL, 3, 3);
58
59 cvCopyImage(image_BIL, temp_image);
60 j++;
61 }
62 cvShowImage("题目_a", image_BIL);
63 char * imageName = "D:\\Work\\Work_Programming\\Source\\Image\\OpenCVExerciseResult\\Result_5-8_Step2.jpg";
64 cvSaveImage(imageName, image_BIL);
65 //---------------------------a：结束--------------------------------//
66
67 cvWaitKey(0);
68
69 cvReleaseImage(&src1);
70 cvReleaseImage(&image_BIL);
71 //cvReleaseImage(&cleandiff);
72 //cvReleaseImage(&dirtydiff);
73
74 cvDestroyWindow("原始图像");
75 cvDestroyWindow("题目_a");
76 //cvDestroyWindow("题目_b");
77 //cvDestroyWindow("题目_c");
78 ////cvDestroyWindow("题目_d");
79
80
81 return 0;
82 }
83
84 //生成的高斯分布随机数序列的期望为0.0，方差为1.0。若指定期望为E，方差为V，则只需增加:X = X * V + E;
85 //double gaussrand()
86 //{
87 // static double V1, V2, S;
88 // static int phase = 0;
89 // double X;
90 //
91 // if (phase == 0) {
92 // do {
93 // double U1 = (double)rand() / RAND_MAX;
94 // double U2 = (double)rand() / RAND_MAX;
95 //
96 // V1 = 2 * U1 - 1;
97 // V2 = 2 * U2 - 1;
98 // S = V1 * V1 + V2 * V2;
99 // } while (S >= 1 || S == 0);
100 //
101 // X = V1 * sqrt(-2 * log(S) / S);
102 // }
103 // else
104 // X = V2 * sqrt(-2 * log(S) / S);
105 //
106 // phase = 1 - phase;
107 //
108 // return X;
109 //}
110
111 //声称期望为E，方差为V的随机数
112 double gaussrand(double E , double V )
113 {
114 static double V1, V2, S;
115 static int phase = 0;
116 double X;
117
118 if (phase == 0) {
119 do {
120 double U1 = (double)rand() / RAND_MAX;
121 double U2 = (double)rand() / RAND_MAX;
122
123 V1 = 2 * U1 - 1;
124 V2 = 2 * U2 - 1;
125 S = V1 * V1 + V2 * V2;
126 } while (S >= 1 || S == 0);
127
128 X = V1 * sqrt(-2 * log(S) / S);
129 }
130 else
131 X = V2 * sqrt(-2 * log(S) / S);
132
133 phase = 1 - phase;
134
135 X = X * V + E;
136 return X;
137 }
138
139 void InitialzieImageToRandomValue(IplImage * img)
140 {
141 int imgWidth = img->width;
142 int imgHeight = img->height;
143
144 if (img->nChannels == 1)
145 {
146 for (int h = 0; h < imgHeight; ++h)
147 {
148 for (int w = 0; w < imgWidth; ++w)
149 {
150 int val = gaussrand(0.0, 3.0);
151 if (val >= 0)
152 {
153 cvSetReal2D(img, w, h, val);
154 }
155 else
156 {
157 w--;
158 }
159 }
160 }
161 }
162 }
结果图片：
时间： 2024-10-10 13:26:04