938 |
return psnr_y; |
return psnr_y; |
939 |
} |
} |
940 |
|
|
941 |
|
|
942 |
|
float sse_to_PSNR(long sse, int pixels) |
943 |
|
{ |
944 |
|
if (sse==0) |
945 |
|
return 99.99F; |
946 |
|
|
947 |
|
return 48.131F - 10*(float)log10((float)sse/(float)(pixels)); // log10(255*255)=4.8131 |
948 |
|
|
949 |
|
} |
950 |
|
|
951 |
|
long plane_sse(uint8_t * orig, |
952 |
|
uint8_t * recon, |
953 |
|
uint16_t stride, |
954 |
|
uint16_t width, |
955 |
|
uint16_t height) |
956 |
|
{ |
957 |
|
int diff, x, y; |
958 |
|
long sse=0; |
959 |
|
|
960 |
|
for (y = 0; y < height; y++) { |
961 |
|
for (x = 0; x < width; x++) { |
962 |
|
diff = *(orig + x) - *(recon + x); |
963 |
|
sse += diff * diff; |
964 |
|
} |
965 |
|
orig += stride; |
966 |
|
recon += stride; |
967 |
|
} |
968 |
|
return sse; |
969 |
|
} |
970 |
|
|
971 |
/* |
/* |
972 |
|
|
973 |
#include <stdio.h> |
#include <stdio.h> |
1064 |
} |
} |
1065 |
|
|
1066 |
|
|
|
#define ABS(X) (((X)>0)?(X):-(X)) |
|
1067 |
float |
float |
1068 |
image_mad(const IMAGE * img1, |
image_mad(const IMAGE * img1, |
1069 |
const IMAGE * img2, |
const IMAGE * img2, |
1229 |
} |
} |
1230 |
|
|
1231 |
} |
} |
1232 |
|
|