Histogram equalization tests

[SeryZone]

Hello to all, who here!
I wrote visualization, thanks for Dunkan C for helping with this idea.

Code:

begin
        index := 0;
        inc(max);
        Histogram := nil; Percentage:=nil;
        SetLength(Histogram, (max-min));
        SetLength(Percentage, (max-min));
        NUM:=0;
        for x:=0 to (maxx*maxy-1) do
          if Iterdat[x]<max then
          inc(Histogram[IterDat[x]-min]);
        //quicks(0, Length(Histogram), Histogram);  //<----------
        for x:=1 to max-min do
          NUM:=NUM+Histogram[x];
        running_total:=0;
        for x:=1 to max-min do
          begin
            running_total := running_total + histogram[x];
            percentage[x] := running_total / NUM;
          end;
          for y := 0 to maxy-1 do
            for x := 0 to maxx-1 do
          begin
            if (IterDat[x*maxy+y] >= max) then
            begin
              buffer[index] := 0; // red
              buffer[index + 1] := 0; // green
              buffer[index + 2] := 0; // blue
            end
            else if (IterDat[x*maxy+y] > 0) and (IterDat[x*Maxy+y]<max) then
              ColorIndex := round( 4096*(percentage[IterDat[x*maxy+y]-min]) ) mod 4096;
            begin
              buffer[index] := pal[ColorIndex].b; // red
              buffer[index + 1] := pal[ColorIndex].g; // green
              buffer[index + 2] := pal[ColorIndex].r; // blue
            end;
            inc(index, 4)
          end;
          dec(max);
      end;

okay, this part have many-many-many bugs, but I don't gave up and made some experiments. How to make rank-ordered list??? If you help me, I will be very glad!

1) I test this part:

Code:

 quicks(0, Length(Histogram), Histogram);  //Sorting Histogram[] array
            ColorIndex := round( 4096*(percentage[IterDat[x*maxy+y]-min]) ) mod 4096;

Result is:



Reminds to me 'Prodigy' spider, ahaha))) 

2) Here, without sort:



3) Here, I tested squaring:

Code:

ColorIndex := round( 4096*sqr(percentage[IterDat[x*maxy+y]-min]) ) mod 4096;

Very nice result: lower density and more details!

4) And last low-density test - fourth power:

Code:

ColorIndex := round( 4096*sqr(sqr(percentage[IterDat[x*maxy+y]-min])) ) mod 4096;

So, If you can help me with solving mistakes - please, help. This equalization very important for me as designer.

P.s. Hope, you rate it...