Реализация 1 - на Паскале.

Если линза на экpане - окpужность (x₀,y₀,r₀), то в точках (x,y), для котоpых (x-x₀)²+(y-y₀)²<=r₀²,pисуется точка, котоpая, если бы линзы не было, изобpажалась бы в точке (x₀+(x-x₀)*k1,y₀+(y-y₀)*k), где k=const1 / ([sqrt]((x-x₀)²+(y-y₀)²) + const2). Можно заpанее пpосчитать таблицу смещений - array[-r₀..r₀,-r₀..r₀] of integer.

{$A+,B-,D+,E+,F-,G+,I+,L+,N-,O-,P-,Q-,R-,S+,T-,V+,X+,Y+}

{$M 16384,0,655360}

Uses CRT;

Const vx0 = 3;
 vy0 = 2;
 v0 = vx0;
 r0 = 50;
 r02 = (r0-v0)*(r0-v0);
 d = r02 * 10 div 10;
Type ScreenType = Array[0..199,0..319] of Byte;
 DispType = Array[-r0..r0,-r0..r0] of Integer;

Var

 Screen : ScreenType Absolute $a000:$0000;
 Buffer1, Buffer2 : ^ScreenType;
 Disp : ^DispType;
 x,y,vx,vy,r2,c : LongInt;
Procedure Move(Var A,B; Count: Word);
assembler;
asm
 push ds
 mov cx, Count
 les di, B
 lds si, A
 shr cx, 1
 jz @zero
 rep movsw
@zero: pop ds
end;
BEGIN

 asm
 mov ax, $13
 int $10
 end;

 New(Buffer1);
 New(Buffer2);
 New(Disp);

 FillChar(Screen, SizeOf(Screen), 3);
 y:=0;
 
 repeat

 For x:=0 to 319 do Screen[y,x]:=11;
 Inc(y,10);

 until y>199;

 x:=0;

 repeat

 For y:=0 to 199 do Screen[y,x]:=11;
 Inc(x,10);

 until x>319;

 Move(Screen, Buffer1^, SizeOf(Screen));
 Move(Buffer1^,Screen,64000);

 For y:=-r0 to r0 do
 For x:=-r0 to r0 do
 begin

 r2:=x*x+y*y;
 if r2>r02 then Disp^[y,x] := y*320+x
 else 
 Disp^[y,x]:=(y*(r2+d)div(r02+d))*320+(x*(r2+d)div(r02+d));

 end;

 x:=r0;
 y:=r0;
 vx:=vx0;
 vy:=vy0;
 repeat

 asm
 mov ax, Integer(y)
 mov bx, 320
 imul bx
 add ax, Integer(x)
 mov di, ax
 mov dx, -r0*320-r0
 les si, Disp
 mov ch, 2*r0+1
@next_dy: mov cl, 2*r0+1
@next_dx:
 mov es, Word(Disp+2)
 mov bx, es:[si]
 mov es, Word(Buffer1+2)
 mov al, es:[di+bx]
 mov bx, Seg(Screen)
 mov es, bx
 mov bx, dx
 mov es:[di+bx], al
 add si, 2
 inc dx
 dec cl
 jnz @next_dx
 add dx, 320-(2*r0+1)
 dec ch
 jnz @next_dy
 end;

 if ((x+vx)>=r0)and ((x+vx)<=319-r0)then Inc(x, vx)
 else vx:=-vx;

 if ((y+vy)>=r0)and ((y+vy)<=199-r0)then Inc(y,vy)
 else vy:=-vy;
 Delay(25);
 until Port[$60]=$01;

 Dispose(Buffer1);
 Dispose(Buffer2);
 Dispose(Disp);

 asm
 mov ax, $03
 int $10
 end;

END.

Реализация 2 - на Watcom C.

Пусть у нас есть матрица , которую нужно преобразовать. То есть у нас есть какой-то прямоугольный массив, который мы хотим показать в линзе. Итак у нас есть исходная матрица. Каждая пара координат в матрице (I,J) переходит в новую пару координат (X,Y). Формулу смотрите в исходнике. Тогда элемент матрицы с новыми координатами (X,Y) становиться равным смещению в матрице, которое соответствует координатам (I,J). Это первый этап, для дальнейшего ускорения вывода , расчитаем еще два массива. Один массив - это смещения по которому мы берем значение цвета точки. А второй массив - смещения по которым мы этот цвет будет ставить. Тогда весь вывод линзы сведеться к циклу, в котором мы берем точку по одному смещению , а ставим по другому.

#include <stdio.h>
#include <math.h>

#define PI 3.1415926
#define max 40

char bmp[64000];
char Scr[64000];
char p[1024];
int matr[2*max][2*max];
int glassnew[2*max][2*max];
int glassold[2*max][2*max];

// Set Video Mode ################################################
extern void setvmode(int);
#pragma aux setvmode = \
 " int 10h " \
 parm [eax] \
 modify exact [eax];

// Wait Ret Race #################################################
void wait_retrace()
 {
 while ((inp(0x3DA) & 0x08) == 0);
 while ((inp(0x3DA) & 0x08) != 0);
 }

// PutPixel ######################################################
void PutPixel(int x,int y,int Color)
 {
 if ((x>0)&&(y>0)&&(x<320)&&(y<200)) Scr[x+y*320]=Color;
 }

// GetPixel ######################################################
int GetPixel(int x,int y)
 {
 return(bmp[x+y*320]);
 }

void main(void)
 {
 FILE *f;
 int i,j,xc,yc,c,yc1,xc1,foffset;
 float s,x,y,a,b;

 setvmode(0x13);
 f=fopen("1.bmp", "rb");
 fseek(f,1024+54,0);
 fread(&bmp, 1, 64000, f);
 fseek(f,54,0);
 fread(&p, 1, 1024, f);
 fclose(f);

 // Готовим палитру
 for(i=0;i<256;i++)
 {
 outp(0x3C8,i);
 outp(0x3C9,p[(i*4)+2] >> 2);
 outp(0x3C9,p[(i*4)+1] >> 2);
 outp(0x3C9,p[(i*4)+0] >> 2);
 }
 // Создаем таблицу переходов
 for (i=-max;i<max;i++)
 for (j=-max;j<max;j++)
 {
 s=sin(PI/2+(PI/2)/(max*2.2)*sqrt(i*i+j*j));
 x=(i*s);
 y=(j*s);
 glassnew[i+max][j+max]=(int)x+(int)y*320;
 glassold[i+max][j+max]=(i/2)+(j/2)*320;
 }
 do {
 memcpy((char*)&Scr, (char*)&bmp, 64000);
 xc=160+100*cos(2*a)*sin(b);
 yc=100+70*cos(b);
 foffset=xc+yc*320;
 for (i=0;i<2*max;i++)
 for (j=0;j<2*max;j++)
 {
 Scr[foffset+glassnew[i][j]]=bmp[foffset+glassold[i][j]];
 }
// wait_retrace();
 memcpy((char*)0xA0000L, (char*)&Scr, 64000);
 a=a+0.01;
 b=b+0.01;
 } while (!kbhit());
 getch();
 setvmode(0x03);
 }