1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:
80:
81:
82:
83:

procedure ResampleBilinear(Src, Dst: TBitmap);
var
  X10, Y10, X11, Y11: Integer;
  X1, Y1: Single;
  Fx, Fy: Single;
  F: Single;
  I: Integer;
  X2, Y2: Integer;
  W1, H1, W2, H2: Integer;
  SrcRow0, SrcRow1, DstRow: PByteArray;
begin
  // Bitmaps auf 24-Bit RGB setzen
  Src.PixelFormat := pf24Bit;
  Dst.PixelFormat := pf24Bit;

  // Größe der Bitmaps speichern
  W1 := Src.Width  - 1;
  H1 := Src.Height - 1;
  W2 := Dst.Width  - 1;
  H2 := Dst.Height - 1;

  // Umrechnungsfaktoren berechnen
  Fx := W1 / W2;
  Fy := H1 / H2;

  // Zeilen durchgehen
  for Y2 := 0 to H2 do
  begin
    { Y1   Position des Pixels im Bereich 0..1
      Y10  Erste Scanline
      Y11  Zweite Scanline }
    Y1 := Y2 * Fy;
    Y10 := Trunc(Y1);
    Y11 := Y10;
    Y1 := Y1 - Y10;
    if (Y1 > 0) then
      Y11 := Y11 + 1;
    { Entspricht folgendem Code:
      Y1  :=  Frac(Y2 * Fy);
      Y10 := Floor(Y2 * Fy);
      Y11 :=  Ceil(Y2 * Fy); }

    // Scanlines laden
    SrcRow0 := Src.ScanLine[Y10];
    SrcRow1 := Src.ScanLine[Y11];
    DstRow  := Dst.ScanLine[Y2];

    for X2 := 0 to W2 do
    begin
      { X1   Position des Pixels im Bereich 0..1
        X10  Erster Pixel
        X11  Zweiter Pixel }
      X1 := X2 * Fx;
      X10 := Trunc(X1);
      X1 := X1 - X10;
      X10 := X10 * 3;
      X11 := X10;
      if (X1 > 0) then
        X11 := X11 + 3;
      { Entspricht folgendem Code:
        X1  :=  Frac(X2 * Fy);
        X10 := Floor(X2 * Fy) * 3;
        X11 :=  Ceil(X2 * Fy) * 3; }

      // Pixel setzen (R/G/B)
      DstRow[X2 * 3    ] := Round(
        SrcRow0[X10    ] * (1 - X1) * (1 - Y1) +
        SrcRow0[X11    ] *      X1  * (1 - Y1) +
        SrcRow1[X10    ] * (1 - X1) *      Y1  +
        SrcRow1[X11    ] *      X1  *      Y1);
      DstRow[X2 * 3 + 1] := Round(
        SrcRow0[X10 + 1] * (1 - X1) * (1 - Y1) +
        SrcRow0[X11 + 1] *      X1  * (1 - Y1) +
        SrcRow1[X10 + 1] * (1 - X1) *      Y1  +
        SrcRow1[X11 + 1] *      X1  *      Y1);
      DstRow[X2 * 3 + 2] := Round(
        SrcRow0[X10 + 2] * (1 - X1) * (1 - Y1) +
        SrcRow0[X11 + 2] *      X1  * (1 - Y1) +
        SrcRow1[X10 + 2] * (1 - X1) *      Y1  +
        SrcRow1[X11 + 2] *      X1  *      Y1);
    end;
  end;
end;