阅读提示:

    《Delphi图像处理》系列以效率为侧重点,一般代码为PASCAL,核心代码采用BASM。

    《C++图像处理》系列以代码清晰,可读性为主,全部使用C++代码。

    尽可能保持二者内容一致,可相互对照。

    本文代码必须包括文章《Delphi图像处理 -- 数据类型及公用过程》中的ImageData.pas单元。

const
_fc0: Single = 0.0;
_fc1: Single = 1.0;
_fc6: Single = 6.0;
_fc60: Single = 60.0;
_fc255: Single = 255.0;
_fc360: Single = 360.0; procedure ColorToHSV(var H, S, V: Single; Color: TARGB);
asm
push eax
push edx
push ecx
movzx ecx, Color.TARGBQuad.Blue
movzx edx, Color.TARGBQuad.Green
movzx eax, Color.TARGBQuad.Red
cmp ecx, edx // ecx = rgbMax
jge @@1 // edx = rgbMin
xchg ecx, edx
@@1:
cmp ecx, eax
jge @@2
xchg ecx, eax
@@2:
cmp edx, eax
cmova edx, eax
cvtsi2ss xmm0, ecx
movaps xmm1, xmm0
divss xmm0, _fc255
pop eax
movss [eax], xmm0 // *V = rgbMax / 255
mov eax, ecx
sub ecx, edx // delta = rgbMax - rgbmin
jnz @@3
pop edx // if (delta == 0)
pop eax // {
mov [eax], ecx // *H = *S = 0
mov [edx], ecx // return
jmp @@Exit // }
@@3:
cvtsi2ss xmm0, ecx
movaps xmm2, xmm0
divss xmm0, xmm1
pop edx
movss [edx], xmm0 // *S = delta / rgbMax
cmp al, Color.TARGBQuad.Red
jne @@5
movzx eax, Color.TARGBQuad.Green
movzx edx, Color.TARGBQuad.Blue
xor ecx, ecx // if (R == rgbMax) eax = G - B; add = 0
jmp @@7
@@5:
cmp al, Color.TARGBQuad.Green
jne @@6
movzx eax, Color.TARGBQuad.Blue
movzx edx, Color.TARGBQuad.Red
mov ecx, 120 // if (G == rgbMax) eax = B - R; add = 120
jmp @@7
@@6:
movzx eax, Color.TARGBQuad.Red
movzx edx, Color.TARGBQuad.Green
mov ecx, 240 // if (B == rgbMax) eax = R - G; add = 240
@@7:
sub eax, edx
cvtsi2ss xmm0, eax
cvtsi2ss xmm1, ecx
mulss xmm0, _fc60
divss xmm0, xmm2
addss xmm0, xmm1 // H = eax * 60 / delta + add
comiss xmm0, _fc0
jae @@8
addss xmm0, _fc360
@@8:
pop eax
movss [eax], xmm0
@@Exit:
end; function HSVToColor(H, S, V: Single): TARGB;
asm
movss xmm1, H
comiss xmm1, _fc0
jae @@1
addss xmm1, _fc360
jmp @@2
@@1:
comiss xmm1, _fc360
jb @@2
subss xmm1, _fc360
@@2:
divss xmm1, _fc60
cvtss2si edx, xmm1 // index = Round(H)
cvtsi2ss xmm2, edx
subss xmm1, xmm2 // extra = H - index
comiss xmm1, _fc0 // if (extra < 0) // 如果index发生五入
jae @@3 // {
dec edx // index --
addss xmm1, _fc1 // extra ++
@@3: // }
test edx, 1
jnz @@4
movaps xmm2, xmm1
movss xmm1, _fc1
subss xmm1, xmm2 // if (index & 1 == 0) extra = 1 - extra
@@4:
movss xmm0, V
movss xmm2, S
minss xmm0, _fc1
minss xmm2, _fc1
maxss xmm0, _fc0
maxss xmm2, _fc0
mulss xmm0, _fc255
mulss xmm2, xmm0 // V * S
mulss xmm1, xmm2 // V * S * extra
pslldq xmm1, 4
movss xmm1, xmm2 // xmm1 = 0 0 V*S*E V*S
pshufd xmm0, xmm0, 0 // xmm0 = V V V V
subps xmm0, xmm1 // xmm0 -= mm1 = V V T P
jmp @@jmpTable[edx*4].Pointer
@@jmpTable: dd offset @@H60
dd offset @@H120
dd offset @@H180
dd offset @@H240
dd offset @@H300
dd offset @@H360
@@H360: // 300 - 359 (V, P, T)
pshufd xmm0, xmm0, 11100001b
jmp @@H60
@@H300: // 240 - 299 (T, P, V)
pshufd xmm0, xmm0, 11010010b
jmp @@H60
@@H240: // 180 - 239 (P, T, V)
pshufd xmm0, xmm0, 11000110b
jmp @@H60
@@H180: // 120 - 179 (P, V, T)
pshufd xmm0, xmm0, 11001001b
jmp @@H60
@@H120: // 60 - 119 (T, V, P)
pshufd xmm0, xmm0, 11011000b
@@H60: // 0 - 59 (V, T, P)
cvtps2dq xmm0, xmm0
packssdw xmm0, xmm0
packuswb xmm0, xmm0
movd eax, xmm0
or eax, 0ff000000h
end;

《Delphi图像处理》系列使用GDI+单元下载地址和说明见文章《GDI+ for VCL基础 -- GDI+ 与 VCL》。

    因水平有限,错误在所难免,欢迎指正和指导。邮箱地址:[email protected]

这里可访问《Delphi图像处理 -- 文章索引》。

05-08 15:21