高手来帮帮我》视频分析(200分)

To pengyt :
能否劳驾给我一份，让我也学习学习。 mayflower6@sina.com

由于pengyt上不了网，现按pengyt的吩咐，将他提供的算法贴出，看行不行。


刚出差回来，由于没有NOTEBOOK，在外地无法上网。
我现在正在做一个眼科的医疗影像处理系统，这些
算法都是从此系统中直接取出来的，各位可能还要
进行一些修改哟

另：现在发现我进不了大富翁论坛了，现象为在输入
名字和密码后，按<进入>，经过一轮等待之后又出现
问名字和密码的画面。怎么回事？

另：请帮忙转发给beilei一份，由于上不去大富翁，
也不知道他的邮箱。


// 计算BITMAP的平均亮度
//以上为计算BITMAP的平均亮度的函数。
//在计算中没有考虑RGB三种颜色的按照人眼的敏感度加权。
//如果只计算某种颜色通道的平均亮度，
//则将total:=total+3改为+1，
//将sum:=sum+p[j*3]+p[j*3+1]+p[j*3+2]
//改为sum:=sum+p[j*3]或p[j*3+1]或p[j*3+2]
//返回值为0-255之间的一个数值。
//要求thebitmap为24bit。
//计算某一区域的平均亮度请自行修改。

function getbmpbrightness(thebitmap: tbitmap): byte;
var
i,j: integer;
p: pbytearray;
hh,ww: integer;
sum: int64;
total: int64;
begin
sum:=0;
total:=0;
with thebitmap do
begin
hh:=thebitmap.height;
ww:=thebitmap.width;
i:=0;
while i<hh do
begin
p:=scanline;
inc(i,3);
for j:=0 to ww-1 do
begin
sum:=sum+p[j*3]+p[j*3+1]+p[j*3+2];// &amp;&amp;&amp;
total:=total+3;// &amp;&amp;&amp;
end;
end;
end;
result:=sum div total;
end;


// bitmap的平滑（模糊），其中，smoothlevel为平滑程度，
// 平滑程度的取值为1、2、3。
// x1,y1,x2,y2为此BITMAP中需要平滑的范围。
// 要求mybitmap为24bit。

procedure bmpSmooth(mybitmap: tbitmap; smoothlevel: byte; x1,y1,x2,y2: integer);
var
p,p1,p2: pbytearray;
i,j: integer;
thew,theh: integer;
s1,s2,s3: integer;
rvalue: integer;
gvalue: integer;
bvalue: integer;
tmpbitmap: tbitmap;
maxnum,minnum: integer;
begin
tmpbitmap:=tbitmap.create;
tmpbitmap.assign(mybitmap);
with mybitmap do
begin
thew:=width;
theh:=height;
if (x1=x2) or (y1=y2) then
begin
x1:=2;
y1:=2;
x2:=thew-2;
y2:=theh-2;
end
else
begin
minnum:=min(y1,y2);
maxnum:=max(y1,y2);
y1:=minnum;
y2:=maxnum;
minnum:=min(x1,x2);
maxnum:=max(x1,x2);
x1:=minnum;
x2:=maxnum;
end;
if x1<2 then x1:=2;
if y1<2 then y1:=2;
if x2>thew-2 then x2:=thew-2;
if y2>theh-2 then y2:=theh-2;
for i:=y1 to y2 do
begin
p:=ScanLine;
p1:=tmpbitmap.scanline[i-1];
p2:=tmpbitmap.scanline[i+1];
for j:=x1 to x2 do
begin
s1:=j*3;
s2:=j*3+1;
s3:=j*3+2;
case smoothlevel of
1: begin // 轻微
rvalue:=(p[s3]*3+p[s3-3]+p[s3+3]+p1[s3]+p1[s3-3]+p1[s3+3]+
p2[s3]+p2[s3-3]+p2[s3+3]) div 11;
gvalue:=(p[s2]*3+p[s2-3]+p[s2+3]+p1[s2]+p1[s2-3]+p1[s2+3]+
p2[s2]+p2[s2-3]+p2[s2+3]) div 11;
bvalue:=(p[s1]*3+p[s1-3]+p[s1+3]+p1[s1]+p1[s1-3]+p1[s1+3]+
p2[s1]+p2[s1-3]+p2[s1+3]) div 11;
end;
2: begin // 普通
rvalue:=(p[s3]*2+p[s3-3]+p[s3+3]+p1[s3]+p1[s3-3]+p1[s3+3]+
p2[s3]+p2[s3-3]+p2[s3+3]) div 10;
gvalue:=(p[s2]*2+p[s2-3]+p[s2+3]+p1[s2]+p1[s2-3]+p1[s2+3]+
p2[s2]+p2[s2-3]+p2[s2+3]) div 10;
bvalue:=(p[s1]*2+p[s1-3]+p[s1+3]+p1[s1]+p1[s1-3]+p1[s1+3]+
p2[s1]+p2[s1-3]+p2[s1+3]) div 10;
end;
3: begin // 加强
rvalue:=(p[s3]+p[s3-3]+p[s3+3]+p1[s3]+p1[s3-3]+p1[s3+3]+
p2[s3]+p2[s3-3]+p2[s3+3]) div 9;
gvalue:=(p[s2]+p[s2-3]+p[s2+3]+p1[s2]+p1[s2-3]+p1[s2+3]+
p2[s2]+p2[s2-3]+p2[s2+3]) div 9;
bvalue:=(p[s1]+p[s1-3]+p[s1+3]+p1[s1]+p1[s1-3]+p1[s1+3]+
p2[s1]+p2[s1-3]+p2[s1+3]) div 9;
end;
else
begin
rvalue:=p[s3];
gvalue:=p[s2];
bvalue:=p[s1];
end;
end;
p[s3]:=rvalue;
p[s2]:=gvalue;
p[s1]:=bvalue;
end;
end;
end;
tmpbitmap.free;
end;

// 直方图均衡
// darkness为黑度，即均衡输出亮度的下限。
// brightness为白度，为均衡输出亮度的上限。
// 可以参考PHOTOSHOP类似功能的概念。
// x1,y1,x2,y2为处理的区域。
procedure grayextend(mybitmap: tbitmap; darkness,brightness: byte;
x1,y1,x2,y2: integer);
var
his: array [0..255,1..3] of real;
ihis: array [0..255,1..3] of extended;
hiseq: array [0..255,1..3] of integer;
sum: extended;
i,j: integer;
thew,theh: integer;
p: pbytearray;
sumreal: array [1..3] of real;
rate: real;
maxnum,minnum: integer;
begin
for i:=0 to 255 do
for j:=1 to 3 do
ihis[i,j]:=0;
sum:=0;
with mybitmap do
begin
thew:=width;
theh:=height;
if (x1=x2) or (y1=y2) then
begin
x1:=1;
y1:=1;
x2:=thew-1;
y2:=theh-1;
end
else
begin
minnum:=min(y1,y2);
maxnum:=max(y1,y2);
y1:=minnum;
y2:=maxnum;
minnum:=min(x1,x2);
maxnum:=max(x1,x2);
x1:=minnum;
x2:=maxnum;
end;
if x1<1 then x1:=1;
if y1<1 then y1:=1;
if x2>thew-1 then x2:=thew-1;
if y2>theh-1 then y2:=theh-1;
for i:=y1 to y2 do
begin
p:=scanline;
for j:=x1 to x2 do
begin
ihis[p[j*3],1]:=ihis[p[j*3],1]+i*i;
ihis[p[j*3+1],2]:=ihis[p[j*3+1],2]+i*i;
ihis[p[j*3+2],3]:=ihis[p[j*3+2],3]+i*i;
sum:=sum+i*i;
end;
end;
for i:=0 to 255 do
for j:=1 to 3 do his[i,j]:=ihis[i,j]/sum;
for i:=0 to 255 do
begin
sumreal[1]:=0;
sumreal[2]:=0;
sumreal[3]:=0;
for j:=0 to i do
begin
sumreal[1]:=sumreal[1]+his[j,1];
sumreal[2]:=sumreal[2]+his[j,2];
sumreal[3]:=sumreal[3]+his[j,3];
end;
rate:=1;
hiseq[i,1]:=round((brightness-darkness)*sumreal[1]*rate+darkness);
hiseq[i,2]:=round((brightness-darkness)*sumreal[2]*rate+darkness);
hiseq[i,3]:=round((brightness-darkness)*sumreal[3]*rate+darkness);
end;
for i:=y1 to y2 do
begin
p:=scanline;
for j:=x1 to x2 do
begin
p[j*3]:=hiseq[p[j*3],1];
p[j*3+1]:=hiseq[p[j*3+1],2];
p[j*3+2]:=hiseq[p[j*3+2],3];
end;
end;
end; // end of with
end;

我知道为何上不了大富翁了：如果你禁止COOKIE，就登录不了；反之则可以。

不行，我还是帮pengyt发给大家，OK

to pengyt ：
算法已收到，多谢。

用红外检测不行吗？

to chaoro:
你说的是用视频采集卡上的红外检测吗？？？
如果不是的，请你说来听听，也可发邮件至
yeiping@yeah.net
谢谢你的提醒！！！！

接受答案了.

to pengyt:
我现在正在做类似的程序,可否EMAIL一份资料给我,谢谢你

t163t163@163.net

to pengyt及楼上各位:
我也在类似程序,哪位看到麻烦EMAIL一份资料给我,谢谢先。：）
Micro_Knife@sina.com