恳切希望大姐大哥解释一下这个有关图像处理的程序，敝人再次向您致敬。非常着急。

#include <math.h>
#include <vcl.h>
//------------------------------亅亅亅亅亅亅亅亅亅亅rgb to XYZ (sRGB婯奿偱)

void sRGBtoLAB(int R1,int G1,int B1,float &L,float &A,float &B)
{
  double R3,G3,B3;   //僨僶僀僗怣崋
  double X,Y,Z;   //嶰巋寖抣
  double Xw,Yw,Zw;//敀怓偺嶰巋寖抣
  double XXw,YYw,ZZw;    //搑拞寁嶼
  double fXXw,fYYw,fZZw; //搑拞寁嶼

  Xw=0.95045;//D65
  Yw=1.0;
  Zw=1.08905;

  R3=(double)R1/255.;//255傪1偵惓婯壔
  G3=(double)G1/255.;
  B3=(double)B1/255.;

  if(R3<=0.03928)R3/=12.92; else R3=pow((R3+0.055)/1.055,2.4);
  if(G3&lt;=0.03928)G3/=12.92; else G3=pow((G3+0.055)/1.055,2.4);
  if(B3&lt;=0.03928)B3/=12.92; else B3=pow((B3+0.055)/1.055,2.4);

  X=0.4124*R3+0.3576*G3+0.1805*B3;
  Y=0.2126*R3+0.7152*G3+0.0722*B3;
  Z=0.0193*R3+0.1192*G3+0.9505*B3;


//xreal=X*100.0; yreal=Y*100.0;  zreal=Z*100.0;



//------------------------------亅亅亅亅亅亅亅亅亅亅XYZ to LAB
  XXw=X/Xw;
  YYw=Y/Yw;
  ZZw=Z/Zw;

  if(XXw&lt;=0.008856)fXXw=7.787*XXw+16.0/116.0; else fXXw=pow(XXw,1.0/3.0);
  if(YYw&lt;=0.008856)fYYw=7.787*YYw+16.0/116.0; else fYYw=pow(YYw,1.0/3.0);
  if(ZZw&lt;=0.008856)fZZw=7.787*ZZw+16.0/116.0; else fZZw=pow(ZZw,1.0/3.0);

  L=116.0*fYYw-16.0;
  A=500.0*(fXXw-fYYw);
  B=200.0*(fYYw-fZZw);
  if(L>100.0)L=100.0;
  if(L<0.0)L=0.0;
}
//------------------------------亅亅亅亅亅亅亅亅亅亅

void LABtosRGB(float Ll,float Aa,float Bb,int &dR,int &dG,int &dB)
{
  double r,g,b;
  double X,Y,Z;   //嶰巋寖抣
  double Xw,Yw,Zw;//敀怓偺嶰巋寖抣
  double XXw,YYw,ZZw;    //搑拞寁嶼
  double fXXw,fYYw,fZZw; //搑拞寁嶼

  Xw=0.95045;//D65
  Yw=1.;
  Zw=1.08905;

  fYYw=(Ll+16.0)/116.0;
  fXXw=Aa/500.0+fYYw;
  fZZw=fYYw-Bb/200.0;

  if(fXXw>0.206896)XXw=pow(fXXw,3.0); else XXw=(fXXw-16.0/116)/7.787;
  if(fYYw&gt;0.206896)YYw=pow(fYYw,3.0); else YYw=(fYYw-16.0/116)/7.787;
  if(fZZw&gt;0.206896)ZZw=pow(fZZw,3.0); else ZZw=(fZZw-16.0/116)/7.787;

  X=XXw*Xw;
  Y=YYw*Yw;
  Z=ZZw*Zw;

// xreal=X*100.0; yreal=Y*100.0;  zreal=Z*100.0;



//------------------------------亅亅亅亅亅亅亅亅亅亅XYZ to rgb (sRGB婯奿偱)
  r= 3.2410*X-1.5374*Y-0.4986*Z;
  g=-0.9692*X+1.8760*Y+0.0416*Z;
  b= 0.0556*X-0.2040*Y+1.0570*Z;

  if(r<=0.00304)r*=12.92; else r=1.055*pow(r,1.0/2.4)-0.055;
  if(g&lt;=0.00304)g*=12.92; else g=1.055*pow(g,1.0/2.4)-0.055;
  if(b&lt;=0.00304)b*=12.92; else b=1.055*pow(b,1.0/2.4)-0.055;

  r*=255.;//1傪255偵惓婯壔
  g*=255.;
  b*=255.;
  dR=(int)r+0.5;
  dG=(int)g+0.5;
  dB=(int)b+0.5;
  if(dR>255)dR=255;
  if(dG&gt;255)dG=255;
  if(dB&gt;255)dB=255;
  if(dR<0)dR=0;
  if(dG&lt;0)dG=0;
  if(dB&lt;0)dB=0;

}
//------------------------------亅亅亅亅亅亅亅亅亅亅r,g,b傪24bit偵偡傞

void makeColor1(TColor color1,int &lutR1,int &lutG1,int &lutB1)
{
    TColor temp_color = color1;
    lutR1=(temp_color >&gt; (0))  & ~(~0 <&lt; 8);
    lutG1=(temp_color >&gt; (8))  & ~(~0 <&lt; 8);
    lutB1=(temp_color >&gt; (16)) & ~(~0 <&lt; 8);
}
//------------------------------亅亅亅亅亅亅亅亅亅亅24bit傪Tcolor宆偵偡傞

void makeRGB1(int lutR,int lutG,int lutB,TColor &color)
{
    int c = 0L;
    c = lutB &lt;&lt; 16; c += lutG &lt;&lt; 8; c += lutR ;
    color = (TColor)c;
}
//------------------------------亅亅亅亅亅亅亅亅亅亅留言


   
TOP
我的空间   联系我们   Powered by UCenter Home 1.2 © 2001-2008 Comsenz Inc.