Маленькая планета с помозью Perlin Noise

[ZBEP]

Здравствуйте!
Необходимо сгенерировать 2D изображение планеты 64х64
Наткнулся на Perlin Noise:
http://freespace.virgin.net/hugo.eli...s/m_perlin.htm
http://wiki.delphigl.com/index.php/Perlin_Noise

Попытался реализовать, но в итоге что-то не то получается..

Вот полный код проекта, рисую прямо на форме:

unit Unit1;

interface

uses
  Winapi.Windows, Winapi.Messages, System.SysUtils, System.Variants, System.Classes, Vcl.Graphics,
  Vcl.Controls, Vcl.Forms, Vcl.Dialogs, Vcl.ExtCtrls;

type
  TForm1 = class(TForm)
    procedure FormClick(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;
  Seed,Number_Of_Octaves,persistence:integer;
  tXY:array [0..256,0..256] of Real;

implementation

{$R *.dfm}

//Создание шумовой функции
function Noise(x,y:integer):Real;
var n:integer;
begin
  n := x + y * 57 + Seed;
  n := ( n shl 13 ) xor n;
  Result := 1 - ( (n * (n * n * 15731 + 789221) + 1376312589) and $7fffffff) / 1073741824;
end;

//Косинусная интерполяция
function InterpoliereCos(a, b, x: Real): Real;
var f:Real;
begin
  result := a*((cos(pi/2*x)+1)/2) + b*(1-(cos(pi/2*x)+1)/2)
end;

//Сглажевание шума
function SmoothedNoise(x,y:integer):Real;
var corners,sides,center:Real;
begin
  corners := ( Noise(x-1, y-1)+Noise(x+1, y-1)+Noise(x-1, y+1)+Noise(x+1, y+1) ) / 16;
  sides   := ( Noise(x-1, y)  +Noise(x+1, y)  +Noise(x, y-1)  +Noise(x, y+1) ) /  8;
  center  :=  Noise(x, y) / 4;
  Result  := corners + sides + center;
end;


function InterpolatedNoise(x,y:Real):Real;
var integer_X,integer_Y:integer;
fractional_X,fractional_Y:Real;
v1,v2,v3,v4,i1,i2:Real;
begin
  integer_X    := Trunc(x);
  fractional_X := x - integer_X;
  integer_Y    := round(y);
  fractional_Y := y - integer_Y;

  v1 := SmoothedNoise(integer_X,     integer_Y);
  v2 := SmoothedNoise(integer_X + 1, integer_Y);
  v3 := SmoothedNoise(integer_X,     integer_Y + 1);
  v4 := SmoothedNoise(integer_X + 1, integer_Y + 1);
  i1 := InterpoliereCos(v1 , v2 , fractional_X);
  i2 := InterpoliereCos(v3 , v4 , fractional_X);
  Result := InterpoliereCos(i1 , i2 , fractional_Y);
end;

function PerlinNoise_2D(x,y:Real):Real;
var i,n,p:integer;
total,frequency,amplitude:Real;
begin
  total := 0;
  p := persistence;
  n := Number_Of_Octaves - 1;
  for i:=0 to n do
  begin
    frequency := 2 xor i;
    amplitude := p xor i;
    total := total + InterpolatedNoise(x * frequency, y * frequency) * amplitude;
  end;
  Result := total
end;

procedure TForm1.FormClick(Sender: TObject);
var x,y:integer;
begin
  Seed:=128;
  Number_Of_Octaves:=4;
  persistence:=25;

  Form1.Canvas.FillRect(rect(0,0,Form1.Width,Form1.Height));
  for x:=0 to 256 do
  for y:=0 to 256 do
  tXY[x,y]:=PerlinNoise_2D(x,y);

  for x:=0 to 256 do
  for y:=0 to 256 do
  if tXY[x,y]>0 then
  Form1.Canvas.Pixels[x,y]:=RGB(round(255*    tXY[x,y] ),round(255*    tXY[x,y] ),round(255*    tXY[x,y])) else
  Form1.Canvas.Pixels[x,y]:=RGB(round(255*abs(tXY[x,y])),round(255*abs(tXY[x,y])),round(255*abs(tXY[x,y])));
end;

end.

Подскажите пожалуйста в чем ошибка.

Прикрепил изображение на котором показал, как сейчас отрабатывает алгоритм.

[ZBEP]

И да, если вдруг есть более подходящие методы, чем Perlin Noise, то подскажите куда копать.

[ZBEP]

Наткнулся на пример когерентного шума и понял что не верно определял цвет пикселя.

[ProQsy]

Спасибо за пример)