ANTICHAT - Показать сообщение отдельно - [PHP][HELP] Скелетолизация изображения

draliokero Изображение мы должны изменять только после завершения итерации. Из-за того что мы насилуем imageGD постоянно - код получается очень медленным, по хорошему - лучше при бинаризации создавать двумерный массив, а потом уже раскрашивать обратно. Работающий код, без оптимизации:

PHP код:


		
			
<?php



set_time_limit(0);





$im = imagecreatefrompng("skeleton.png");

$width = imagesx($im);

$height = imagesy($im);

for ($y = 0; $y < $height; $y++) {

    for ($x = 0; $x < $width; $x++) {

        if (imagecolorat($im, $x , $y) <= 0x7FFFFF) imagesetpixel($im, $x, $y, 0x000000);

        else imagesetpixel($im, $x, $y, 0xFFFFFF);

    } 

} 

// Функция скелетонизации бинарного изображения

function zhang_suen_thinning_iteration($condition_switch)

{

    global $width, $height, $im;

    

    $im2 = imagecreatetruecolor($width, $height);

    

    $running = false;

    for($y = 1; $y < $height - 1; $y++) {

        for($x = 1; $x < $width - 1; $x++) {

            $p1 = imagecolorat($im, $x, $y) > 0 ? 1 : 0;

        

            if ($p1 > 0) {

            imagesetpixel($im2,$x,$y,0xFFFFFF);

                $p2 = imagecolorat($im, $x, $y - 1) > 0 ? 1 : 0;

                $p3 = imagecolorat($im, $x + 1, $y - 1) > 0 ? 1 : 0;

                $p4 = imagecolorat($im, $x + 1, $y) > 0 ? 1 : 0;

                $p5 = imagecolorat($im, $x + 1, $y + 1) > 0 ? 1 : 0;

                $p6 = imagecolorat($im, $x, $y + 1) > 0 ? 1 : 0;

                $p7 = imagecolorat($im, $x - 1, $y + 1) > 0 ? 1 : 0;

                $p8 = imagecolorat($im, $x - 1, $y) > 0 ? 1 : 0;

                $p9 = imagecolorat($im, $x - 1, $y - 1) > 0 ? 1 : 0;

                $connectivity =

                ($p2 == 0 && $p3 == 1 ? 1 : 0) +

                ($p3 == 0 && $p4 == 1 ? 1 : 0) +

                ($p4 == 0 && $p5 == 1 ? 1 : 0) +

                ($p5 == 0 && $p6 == 1 ? 1 : 0) +

                ($p6 == 0 && $p7 == 1 ? 1 : 0) +

                ($p7 == 0 && $p8 == 1 ? 1 : 0) +

                ($p8 == 0 && $p9 == 1 ? 1 : 0) +

                ($p9 == 0 && $p2 == 1 ? 1 : 0);

                $non_zero_neighbors = $p2 + $p3 + $p4 + $p5 + $p6 + $p7 + $p8 + $p9;



                if (2 <= $non_zero_neighbors && $non_zero_neighbors <= 6 && $connectivity == 1) {

                    if ($condition_switch == true) {

                        if (($p2 * $p4 * $p6 == 0) && ($p4 * $p6 * $p8 == 0)) {

                            imagesetpixel($im2, $x, $y, 0x000000);

                            $running = true;

                        } 

                    } else {

                        if (($p2 * $p4 * $p8 == 0) && ($p2 * $p6 * $p8 == 0)) {

                            imagesetpixel($im2, $x, $y, 0x000000);

                            $running = true;

                        } 

                    } 

                } 

            } else imagesetpixel($im2, $x, $y, 0x000000);

        } 

    } 

    imagecopy($im,$im2,0,0,0,0,$width,$height);

    return $running;

} 

$flag = true;$count = 0;



//zhang_suen_thinning_iteration($im, $flag);

while (zhang_suen_thinning_iteration($im, $flag = !$flag) && $count++ < 60) {}



header('Content-type: image/png');

imagepng($im);

За 60 итераций получаем такое изображение из левой верхней картинки первого поста (С увеличением чиасла итераций треугольник конечно пропадет, но достаточно не быстро):

З.Ы. получить ЧБ изображение можно через среднее арифметическое трех каналов цветного (что, кстати, и будет яркостью пиксела)