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:

procedure TForm1.Button1Click(Sender: TObject);
const bb=960;
var i,j,waagerecht,senkrecht:integer;
    xk,yk,zk,q,xi,yi,sq :double;
    xm,ym:integer;
    bitmap:tbitmap;
    rect:trect;
    P : PINTEGERArray;
begin
   Button1.Caption := 'Warte';
   //Verschiebung des Spiegelsmittelpunkts
   waagerecht:=bb;//spinedit1.value;
   senkrecht:=bb;//spinedit2.value;

   bitmap:=tbitmap.Create;
   bitmap.width:=bb;
   bitmap.height:=bb;
   //Abbildung auf doppelte Größe strecken
   rect.left:=0;
   rect.Top:=0;
   rect.right:=bb;
   rect.Bottom:=bb;
   bitmap.canvas.copyrect(rect,image1.picture.bitmap.Canvas,image1.clientrect);
   //Mittelpunkt
   xm:=paintbox1.Width div 2;
   ym:=paintbox1.height div 2;
   paintbox1.Canvas.Brush.Color:=clwhite;
   paintbox1.Canvas.Rectangle(-1,-1,961,961);
   //zeilenweise pixel ermitteln
   for i:=0 to bb-1 do begin
     P := BitMap.ScanLine[i];
     yi:=(2*i-senkrecht)/bb;
     for j:=0 to bb-1 do begin
       //Umwandeln in komplexe Koordinaten
       xi:=(2*j-waagerecht)/bb;
       //Transformation auf die Riemann-Kugel
       sq := sqr(xi)+sqr(yi);
       IF sq < 1.0 then
       begin
         q:=1/(1+sq);
         zk:=sq*q;
         xk:=xi*q;
         yk:=yi*q;
         //nur zeichnen, wenn in unterer Hälfte der Riemann-Kugel
//        if zk<0.5 then
//         paintbox1.Canvas.pixels[trunc(xm+480*xk),trunc(ym+480*yk)]:=rgb(p[4*j+2],p[4*j+1],p[4*j]);
         paintbox1.Canvas.pixels[trunc(xm+480*xk),trunc(ym+480*yk)]:=p[j];
       end;
     end;
   end;
   Button1.Caption := 'Start';
   bitmap.free;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:

   bitmap2:=tbitmap.create;
   bitmap2.width:=image1.Width;
   bitmap2.height:=image1.height;
  
   bitmap2.canvas.draw(0,0,image1.picture.Bitmap);
oder
   rect2.left:=0;
   rect2.Top:=0;
   rect2.right:=bitmap2.width;
   rect2.Bottom:=bitmap2.height; 
   bitmap2.canvas.copyrect(rect2,image1.picture.Bitmap.canvas,rect2);					

1:

   bitmap2.canvas.copyrect(image1.clientrect,image1.canvas,image1.clientrect);					

1:
2:
3:
4:
5:
6:
7:
8:

jpeg: TJPEGImage;
bmp: TBitmap;

  jpeg := TJPEGImage.Create;
  jpeg.LoadFromFile(Filename);
  bmp := TBitmap.Create;
  bmp.Assign(jpeg);
  bmp.PixelFormat := pf32bit; // sicherheitshalber erzwingen (default für jpg ist vermutlich pf24bit)					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:

d = 70; % Distanz der Kugel
R = 30; % Kugelradius
alpha = 6*pi/180; % Bogenmaß

% Sinussatz
beta =  asin(sin(alpha)*(d+R)/R); % Winkel oben im Dreieck zwischen Kamera,
                  % Kugelmittelpunkt und Reflexpunkt                 
% asin gibt hier das "kleine" Ergebnis aus, wir brauchen aber das Große...
beta = pi-beta;

gamma = pi - alpha - beta; % Innenwinkelsumme Dreieck

% pi - beta ist der Winkel zwischen dem Lot und dem Eingangssttrahl.
% Eingangswinkel = Ausgangwinkel, daher verdoppeln
phi = 2*(pi - beta);

% Ausgangsstrahllänge bis zum Schnittpunkt mit der Kugel
c = R * sin(gamma)/sin(alpha);

theta = pi/2 - alpha; % Komplementärwinkel
psi = pi - theta - phi; % Wenn negativ, Spiegelung nicht mehr in der Ebene

x = c * sin(phi)/sin(psi) % Distanz an der "Wand"