#include"cv.h"
#include"cxcore.h"
#include"highgui.h"
#include<iostream>
#include<math.h>
#include <stdio.h>
using namespace std;
int main(int argc, char** argv)
{
IplImage* src = cvLoadImage( "E:\\matlab\\test6.bmp" );
IplImage * I = cvCreateImage( cvGetSize(src), 8, 1);
IplImage * SVLM = cvCreateImage( cvGetSize(src), 8, 1);
IplImage * gama = cvCreateImage( cvGetSize(src), 8, 1);
IplImage * O = cvCreateImage( cvGetSize(src), 8, 1);
IplImage * E = cvCreateImage( cvGetSize(src), 8, 1);
IplImage * S = cvCreateImage( cvGetSize(src), 8, 1);
double P;
double a=0.5;
double f=0.7;
//float gama;
IplImage* r = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* g = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* b = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* R = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* R1 = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* G = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* G1 = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* B = cvCreateImage( cvGetSize(src), 8, 1 );
IplImage* B1 = cvCreateImage( cvGetSize(src), 8, 1 );
cvCvtPixToPlane( src, r, g, b, 0 );
IplImage* final = cvCreateImage( cvGetSize(src), 8, 3 );
//cvSetImageCOI(src,1);
//cvCopy(src,BlueChannel); //提取蓝色
//cvSetImageCOI(src,2);
//cvCopy(src,GreenChannel); //提取绿色
//cvSetImageCOI(src,3);
//cvCopy(src,RedChannel); //提取红色
cvCvtColor(src,I,CV_RGB2GRAY); //提取亮度
IplImage * L = cvCreateImage( cvGetSize(I), 8, 1);
cvSmooth( I , L , CV_GAUSSIAN , 3 ,3,0,0);
CvScalar MeanScalar;
CvScalar StandardDeviationScalar;
cvAvgSdv(I,&MeanScalar,&StandardDeviationScalar);
//printf("StandardDeviationScalar:%f\n",StandardDeviationScalar.val[0]);
if(StandardDeviationScalar.val[0]<=40)
P=2;
else if(40<StandardDeviationScalar.val[0]<80)
P=-0.025*StandardDeviationScalar.val[0]+3;
else if(StandardDeviationScalar.val[0]>80)
P=1;
//printf("P:%f\n",P);
IplImage* I2 = 0 ;
CvSize I2_cvsize; //目标图像尺寸
I2_cvsize.width = I->width / 2; //目标图像的宽为源图象宽的scale倍
I2_cvsize.height = I->height / 2; //目标图像的高为源图象高的scale倍
I2 = cvCreateImage( I2_cvsize, 8, 1); //构造目标图象
cvResize(I, I2, CV_INTER_LINEAR); //缩放源图像到目标图像,线性插值
IplImage * L_2 = cvCreateImage( cvGetSize(I2), 8, 1);
cvSmooth( I2 , L_2 , CV_GAUSSIAN , 3 ,3,0,0);
IplImage * L2 = cvCreateImage( cvGetSize(I), 8, 1);
cvResize(L_2, L2, CV_INTER_LINEAR);
IplImage* I3 = 0 ;
CvSize I3_cvsize; //目标图像尺寸
I3_cvsize.width = I->width / 4; //目标图像的宽为源图象宽的scale倍
I3_cvsize.height = I->height / 4; //目标图像的高为源图象高的scale倍
I3 = cvCreateImage( I3_cvsize, 8, 1); //构造目标图象
cvResize(I, I3, CV_INTER_LINEAR); //缩放源图像到目标图像,线性差值
IplImage * L_3 = cvCreateImage( cvGetSize(I3), 8, 1);
cvSmooth( I3 , L_3 , CV_GAUSSIAN , 3 ,3,0,0);
IplImage * L3 = cvCreateImage( cvGetSize(I), 8, 1);
cvResize(L_3, L3, CV_INTER_LINEAR);
IplImage* I4 = 0 ;
CvSize I4_cvsize; //目标图像尺寸
I4_cvsize.width = I->width / 8; //目标图像的宽为源图象宽的scale倍
I4_cvsize.height = I->height / 8; //目标图像的高为源图象高的scale倍
I4 = cvCreateImage( I4_cvsize, 8, 1); //构造目标图象
cvResize(I, I4, CV_INTER_LINEAR); //缩放源图像到目标图像,线性差值
IplImage * L_4 = cvCreateImage( cvGetSize(I4), 8, 1);
cvSmooth( I4 , L_4 , CV_GAUSSIAN , 3 ,3,0,0);
IplImage * L4 = cvCreateImage( cvGetSize(I), 8, 1);
cvResize(L_4, L4, CV_INTER_LINEAR);
int height = I->height;
int width = I->width;
int step = I->widthStep;
//double f1,f2;
for(int row=0;row<height;row++){
uchar* ptr=(uchar*)(SVLM->imageData+row*step);
uchar* ptr1=(uchar*)(L->imageData+row*step);
uchar* ptr2=(uchar*)(L2->imageData+row*step);
uchar* ptr3=(uchar*)(L3->imageData+row*step);
uchar* ptr4=(uchar*)(L4->imageData+row*step);
uchar* ptr5=(uchar*)(gama->imageData+row*step);
uchar* ptr6=(uchar*)(O->imageData+row*step);
uchar* ptr7=(uchar*)(I->imageData+row*step);
uchar* ptr8=(uchar*)(E->imageData+row*step);
uchar* ptr9=(uchar*)(S->imageData+row*step);
uchar* ptr10=(uchar*)(r->imageData+row*step);
uchar* ptr11=(uchar*)(R->imageData+row*step);
uchar* ptr12=(uchar*)(g->imageData+row*step);
uchar* ptr13=(uchar*)(G->imageData+row*step);
uchar* ptr14=(uchar*)(b->imageData+row*step);
uchar* ptr15=(uchar*)(B->imageData+row*step);
for(int col=0;col<width;col++){
ptr[col]=(ptr1[col]+ptr2[col]+ptr3[col]+ptr4[col])/4;
ptr5[col]=(double)pow(0.54,(128.0-(double)ptr[col])/128.0);
ptr6[col]=(double)255*pow((double)(ptr7[col]/255.0),0.5);//(double)ptr5[col]
//if(col==0&row==0) printf("P:%e\n",ptr6[col]);
ptr8[col]=(double)pow(((double)(ptr[col])*1.0)/((double)(ptr6[col])*1.0),P);
//if(col==0&row==0) printf("P:%e\n",ptr8[col]);
ptr9[col]=(double)255.0*pow(ptr6[col]/255.0,ptr8[col]);
ptr11[col]=(double)ptr9[col]*((double)(ptr10[col])*1.0)/((double)(ptr6[col])*1.0)*f;
//if(ptr11[col]==0) printf("%d %d ",row,col);
ptr13[col]=(double)ptr9[col]*((double)(ptr12[col])*1.0)/((double)(ptr6[col])*1.0)*f;
ptr15[col]=(double)ptr9[col]*((double)(ptr14[col])*1.0)/((double)(ptr6[col])*1.0)*f;
}
}
//cvDiv(r,O,R1,1.0);
//cvMul(R1,S,R,f);
//cvDiv(g,O,G1,1.0);
//cvMul(G1,S,G,f);
//cvDiv(b,O,B1,1.0);
//cvMul(B1,S,B,f);
cvCvtPlaneToPix(R,G,B,0,final);
//cvNamedWindow( "r", 1 );
//cvNamedWindow( "g", 1 );
//cvNamedWindow( "b", 1 );
//cvNamedWindow( "I", 1 );
//cvNamedWindow( "I2", 1 );
//cvNamedWindow( "I3", 1 );
//cvNamedWindow( "I4", 1 );
//cvNamedWindow( "L", 1 );
//cvNamedWindow( "L2", 1 );
//cvNamedWindow( "L3", 1 );
//cvNamedWindow( "L4", 1 );
//cvNamedWindow( "SVLM", 1 );
//cvNamedWindow( "S", 1 );
//cvNamedWindow( "O", 1 );
//cvNamedWindow( "E", 1 );
//cvNamedWindow( "R", 1 );
//cvNamedWindow( "G", 1 );
//cvNamedWindow( "B", 1 );
cvNamedWindow( "final", 1 );
//cvShowImage( "r", r );
//cvShowImage( "g", g );
//cvShowImage( "b", b);
//cvShowImage( "I", I );
//cvShowImage( "I2", I2 );
//cvShowImage( "I3", I3 );
//cvShowImage( "I4", I4 );
//cvShowImage( "L", L );
//cvShowImage( "L2", L2 );
//cvShowImage( "L3", L3 );
//cvShowImage( "L4", L4);
//cvShowImage( "SVLM", SVLM);
//cvShowImage( "S", S);
//cvShowImage( "O", O);
//cvShowImage( "E", E);
//cvShowImage( "R", R);
//cvShowImage( "G", G);
//cvShowImage( "B", B);
cvShowImage( "final", final);
cvSaveImage("E:\\matlab\\image.bmp",final);
//cvSaveImage("E:\\matlab\\Oimage.bmp",O);
//cvSaveImage("E:\\matlab\\Iimage.bmp",I);
//cvSaveImage("E:\\matlab\\SVLM.bmp",SVLM);
cvWaitKey(0);
//cvReleaseImage( &r );
//cvReleaseImage( &g );
//cvReleaseImage( & b );
//cvReleaseImage( & I );
//cvReleaseImage( & I2 );
//cvReleaseImage( & I3 );
//cvReleaseImage( & I4 );
//cvReleaseImage( & L );
//cvReleaseImage( & L2 );
//cvReleaseImage( & L3 );
//cvReleaseImage( & L4 );
//cvReleaseImage( & SVLM );
//cvReleaseImage( & S );
//cvReleaseImage( & O );
//cvReleaseImage( & E );
//cvReleaseImage( & R );
//cvReleaseImage( & G );
//cvReleaseImage( & B );
cvReleaseImage( & final );
return 0;
}
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
new.rar (88个子文件)
new
8.jpg 14KB
new.sdf 35.77MB
Release
new
Release
CL.write.1.tlog 2B
CL.read.1.tlog 2B
cl.command.1.tlog 2B
new.Build.CppClean.log 155B
new.unsuccessfulbuild 0B
new.lastbuildstate 40B
new.log 856B
vc100.pdb 36KB
new.vcxproj.filters 952B
new.vcxproj 5KB
new.cpp 7KB
new.vcxproj.user 143B
Debug
link.3064-cvtres.write.1.tlog 2B
vc100.idb 851KB
link.3880-cvtres.write.1.tlog 2B
link.4504-cvtres.read.1.tlog 2B
link.4468-cvtres.write.1.tlog 2B
link.4944-cvtres.read.1.tlog 2B
CL.write.1.tlog 176B
link.3064.write.1.tlog 2B
link.4504-cvtres.write.1.tlog 2B
CL.read.1.tlog 30KB
mt.read.1.tlog 198B
rc.write.1.tlog 178B
link.1816.read.1.tlog 2B
new.exe.intermediate.manifest 381B
rc.read.1.tlog 170B
link.724.read.1.tlog 2B
link.724.write.1.tlog 2B
link.3620-cvtres.write.1.tlog 2B
link.4504.write.1.tlog 2B
link.3880.write.1.tlog 2B
link.5424.read.1.tlog 2B
link.5424-cvtres.write.1.tlog 2B
mt.command.1.tlog 308B
link.1816.write.1.tlog 2B
link.3880.read.1.tlog 2B
link.3064.read.1.tlog 2B
cl.command.1.tlog 542B
new.Build.CppClean.log 1KB
link-cvtres.read.1.tlog 2B
link.4468.read.1.tlog 2B
link.5424-cvtres.read.1.tlog 2B
link.3620.read.1.tlog 2B
link.3620-cvtres.read.1.tlog 2B
link.4944-cvtres.write.1.tlog 2B
link.1816-cvtres.read.1.tlog 2B
link.1820.read.1.tlog 2B
link.write.1.tlog 420B
link.3620.write.1.tlog 2B
new.exe.embed.manifest.res 472B
link.4944.write.1.tlog 2B
link.6032-cvtres.write.1.tlog 2B
link.5424.write.1.tlog 2B
new.exe.embed.manifest 406B
link-cvtres.write.1.tlog 2B
link.1816-cvtres.write.1.tlog 2B
link.724-cvtres.write.1.tlog 2B
link.3880-cvtres.read.1.tlog 2B
link.command.1.tlog 2KB
link.724-cvtres.read.1.tlog 2B
rc.command.1.tlog 384B
link.read.1.tlog 4KB
new.lastbuildstate 38B
new_manifest.rc 196B
link.4468-cvtres.read.1.tlog 2B
link.1820-cvtres.write.1.tlog 2B
new.obj 107KB
mt.write.1.tlog 198B
link.1820.write.1.tlog 2B
link.3064-cvtres.read.1.tlog 2B
link.6032.read.1.tlog 2B
link.4468.write.1.tlog 2B
link.1820-cvtres.read.1.tlog 2B
new.log 3KB
link.4944.read.1.tlog 2B
vc100.pdb 884KB
link.4504.read.1.tlog 2B
link.6032-cvtres.read.1.tlog 2B
link.6032.write.1.tlog 2B
new.suo 17KB
new.sln 876B
ipch
new-7e9874d8
Debug
new.pdb 1.27MB
new.exe 34KB
after.bmp 392KB
new.ilk 355KB
共 88 条
- 1
资源评论
- qingting_1282013-07-30基于opencv的,使用颜色直方图进行图像匹配,效果还可以
- shangcanghejun2013-03-17很好,不错的,vs2010搭载opencv,很值得参考
wyql123
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功