58 lines
1.3 KiB
C++
58 lines
1.3 KiB
C++
#include<iostream>
|
|
#include<cstdio>
|
|
#include<cmath>
|
|
#define H 1210
|
|
#define W 1240
|
|
|
|
using namespace std;
|
|
|
|
int main()
|
|
{
|
|
unsigned char* grey;
|
|
// 存灰度值的数组
|
|
grey = new unsigned char[H * W];
|
|
FILE* filetest;
|
|
//读取文件
|
|
fopen_s(&filetest, "E:\\Digital image process\\lab-1\\第一次上机实验数据\\test.raw", "rb");
|
|
fread(grey, sizeof(unsigned char), H * W, filetest);
|
|
|
|
int i, j;
|
|
int* Pi_number = new int[256];
|
|
for (i = 0; i < 256; i++) {
|
|
Pi_number[i] = 0;
|
|
}
|
|
for (i = 0; i < H; i++) {
|
|
for (j = 0; j < W; j++) {
|
|
Pi_number[grey[i * W + j]]++;
|
|
}
|
|
}
|
|
float* Pi = new float[256];
|
|
for (i = 0; i < 256; i++) {
|
|
Pi[i] = (float)Pi_number[i] / (H * W);
|
|
}
|
|
|
|
//计算均值和标准差
|
|
float mean = 0;
|
|
for (i = 0; i < 256; i++) {
|
|
mean += (float)i * Pi[i];
|
|
}
|
|
|
|
float variance = 0;
|
|
for (i = 0; i < 256; i++) {
|
|
variance += Pi[i] * (i - mean) * (i - mean);
|
|
}
|
|
float stdvar = sqrt(variance);
|
|
|
|
FILE* fileout;
|
|
|
|
fopen_s(&fileout, "E:\\Digital image process\\lab-1\\第一次上机实验数据\\result.txt", "w");
|
|
|
|
fprintf(fileout, "完成人：测绘2401 24110901003 毕爽爽 \n");
|
|
fprintf(fileout, "均值：%f，标准差：%f\n", mean, stdvar);
|
|
for (i = 0; i < 256; i++) {
|
|
fprintf(fileout, "灰度值为%d的像素个数：%d,频率：%f\n", i, Pi_number[i], Pi[i]);
|
|
}
|
|
|
|
fclose(filetest);
|
|
fclose(fileout);
|
|
} |