Te dan dos imágenes en color verdadero, la Fuente y la Paleta. No necesariamente tienen las mismas dimensiones, pero se garantiza que sus áreas son las mismas, es decir, tienen el mismo número de píxeles.

Su tarea es crear un algoritmo que haga la copia más precisa de la Fuente utilizando solo los píxeles en la Paleta. Cada píxel en la paleta debe usarse exactamente una vez en una posición única en esta copia. La copia debe tener las mismas dimensiones que la Fuente.

Este script de Python se puede usar para garantizar que se cumplan estas restricciones:

from PIL import Image
def check(palette, copy):
    palette = sorted(Image.open(palette).convert('RGB').getdata())
    copy = sorted(Image.open(copy).convert('RGB').getdata())
    print 'Success' if copy == palette else 'Failed'

check('palette.png', 'copy.png')

Aquí hay varias fotos para probar. Todos tienen la misma área. Su algoritmo debería funcionar para dos imágenes con áreas iguales, no solo American Gothic y Mona Lisa. Por supuesto, debe mostrar su salida.

_{^{Gracias a Wikipedia por las imágenes de cuadros famosos.}}

Puntuación

Este es un concurso de popularidad, por lo que gana la respuesta más votada. ¡Pero estoy seguro de que hay muchas maneras de ser creativo con esto!

Animación

Millinon tuvo la idea de que sería genial ver que los píxeles se reorganizaran. Yo también lo pensé, así que escribí este script de Python que toma dos imágenes hechas de los mismos colores y dibuja las imágenes intermedias entre ellas. Actualización: acabo de revisarlo para que cada píxel se mueva la cantidad mínima que debe. Ya no es al azar.

Primero, la Mona Lisa se está convirtiendo en el gótico americano de aditsu. El siguiente es el gótico americano de bitpwner (de Mona Lisa) que se convierte en aditsu. Es sorprendente que las dos versiones compartan exactamente la misma paleta de colores.

Los resultados son realmente asombrosos. Aquí está la mona Lisa de aditsu (ralentizada para mostrar detalles).

Esta última animación no está necesariamente relacionada con el concurso. Muestra lo que sucede cuando mi script se usa para rotar una imagen 90 grados.

Java - GUI con transformación aleatoria progresiva

Probé MUCHAS cosas, algunas de ellas muy complicadas, y finalmente volví a este código relativamente simple:

import java.awt.BorderLayout;
import java.awt.event.ActionEvent;
import java.awt.event.ActionListener;
import java.awt.image.BufferedImage;
import java.io.File;
import java.io.IOException;
import java.util.Random;

import javax.imageio.ImageIO;
import javax.swing.ImageIcon;
import javax.swing.JButton;
import javax.swing.JFrame;
import javax.swing.JLabel;
import javax.swing.Timer;

@SuppressWarnings("serial")
public class CopyColors extends JFrame {
    private static final String SOURCE = "spheres";
    private static final String PALETTE = "mona";
    private static final int COUNT = 10000;
    private static final int DELAY = 20;
    private static final int LUM_WEIGHT = 10;

    private static final double[] F = {0.114, 0.587, 0.299};
    private final BufferedImage source;
    protected final BufferedImage dest;
    private final int sw;
    private final int sh;
    private final int n;
    private final Random r = new Random();
    private final JLabel l;

    public CopyColors(final String sourceName, final String paletteName) throws IOException {
        super("CopyColors by aditsu");
        source = ImageIO.read(new File(sourceName + ".png"));
        final BufferedImage palette = ImageIO.read(new File(paletteName + ".png"));
        sw = source.getWidth();
        sh = source.getHeight();
        final int pw = palette.getWidth();
        final int ph = palette.getHeight();
        n = sw * sh;
        if (n != pw * ph) {
            throw new RuntimeException();
        }
        dest = new BufferedImage(sw, sh, BufferedImage.TYPE_INT_RGB);
        for (int i = 0; i < sh; ++i) {
            for (int j = 0; j < sw; ++j) {
                final int x = i * sw + j;
                dest.setRGB(j, i, palette.getRGB(x % pw, x / pw));
            }
        }
        l = new JLabel(new ImageIcon(dest));
        add(l);
        final JButton b = new JButton("Save");
        add(b, BorderLayout.SOUTH);
        b.addActionListener(new ActionListener() {
            @Override
            public void actionPerformed(final ActionEvent e) {
                try {
                    ImageIO.write(dest, "png", new File(sourceName + "-" + paletteName + ".png"));
                } catch (IOException ex) {
                    ex.printStackTrace();
                }
            }
        });
    }

    protected double dist(final int x, final int y) {
        double t = 0;
        double lx = 0;
        double ly = 0;
        for (int i = 0; i < 3; ++i) {
            final double xi = ((x >> (i * 8)) & 255) * F[i];
            final double yi = ((y >> (i * 8)) & 255) * F[i];
            final double d = xi - yi;
            t += d * d;
            lx += xi;
            ly += yi;
        }
        double l = lx - ly;
        return t + l * l * LUM_WEIGHT;
    }

    public void improve() {
        final int x = r.nextInt(n);
        final int y = r.nextInt(n);
        final int sx = source.getRGB(x % sw, x / sw);
        final int sy = source.getRGB(y % sw, y / sw);
        final int dx = dest.getRGB(x % sw, x / sw);
        final int dy = dest.getRGB(y % sw, y / sw);
        if (dist(sx, dx) + dist(sy, dy) > dist(sx, dy) + dist(sy, dx)) {
            dest.setRGB(x % sw, x / sw, dy);
            dest.setRGB(y % sw, y / sw, dx);
        }
    }

    public void update() {
        l.repaint();
    }

    public static void main(final String... args) throws IOException {
        final CopyColors x = new CopyColors(SOURCE, PALETTE);
        x.setSize(800, 600);
        x.setLocationRelativeTo(null);
        x.setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE);
        x.setVisible(true);
        new Timer(DELAY, new ActionListener() {
            @Override
            public void actionPerformed(final ActionEvent e) {
                for (int i = 0; i < COUNT; ++i) {
                    x.improve();
                }
                x.update();
            }
        }).start();
    }
}

Todos los parámetros relevantes se definen como constantes al comienzo de la clase.

El programa primero copia la imagen de la paleta en las dimensiones fuente, luego elige repetidamente 2 píxeles aleatorios y los intercambia si eso los acerca más a la imagen fuente. "Más cerca" se define utilizando una función de distancia de color que calcula la diferencia entre los componentes r, g, b (con ponderación de luma) junto con la diferencia de luma total, con un mayor peso para la luma.

Las formas tardan unos segundos en formarse, pero un poco más de tiempo para que los colores se unan. Puede guardar la imagen actual en cualquier momento. Por lo general, esperaba entre 1 y 3 minutos antes de guardar.

Resultados:

A diferencia de otras respuestas, todas estas imágenes se generaron utilizando exactamente los mismos parámetros (distintos de los nombres de archivo).

Paleta gótica americana

Paleta Mona Lisa

Paleta de noche estrellada

La paleta Scream

Paleta de esferas

Creo que esta es la prueba más difícil y todos deberían publicar sus resultados con esta paleta:

Lo siento, la imagen del río no me pareció muy interesante, así que no la he incluido.

También agregué un video en https://www.youtube.com/watch?v=_-w3cKL5teM , muestra lo que hace el programa (no exactamente en tiempo real sino similar) y luego muestra el movimiento gradual de píxeles usando la pitón de Calvin guión. Lamentablemente, la calidad del video está significativamente dañada por la codificación / compresión de YouTube.

Java

import java.awt.Point;
import java.awt.image.BufferedImage;
import java.io.File;
import java.io.IOException;
import java.util.ArrayList;
import java.util.Collections;
import java.util.Comparator;
import java.util.HashSet;
import java.util.Iterator;
import java.util.LinkedList;
import java.util.List;
import java.util.Random;
import javax.imageio.ImageIO;

/**
 *
 * @author Quincunx
 */
public class PixelRearranger {

    public static void main(String[] args) throws IOException {
        BufferedImage source = ImageIO.read(resource("American Gothic.png"));
        BufferedImage palette = ImageIO.read(resource("Mona Lisa.png"));
        BufferedImage result = rearrange(source, palette);
        ImageIO.write(result, "png", resource("result.png"));
        validate(palette, result);
    }

    public static class MInteger {
        int val;

        public MInteger(int i) {
            val = i;
        }
    }

    public static BufferedImage rearrange(BufferedImage source, BufferedImage palette) {
        BufferedImage result = new BufferedImage(source.getWidth(),
                source.getHeight(), BufferedImage.TYPE_INT_RGB);

        //This creates a list of points in the Source image.
        //Then, we shuffle it and will draw points in that order.
        List<Point> samples = getPoints(source.getWidth(), source.getHeight());
        System.out.println("gotPoints");

        //Create a list of colors in the palette.
        rgbList = getColors(palette);
        Collections.sort(rgbList, rgb);
        rbgList = new ArrayList<>(rgbList);
        Collections.sort(rbgList, rbg);
        grbList = new ArrayList<>(rgbList);
        Collections.sort(grbList, grb);
        gbrList = new ArrayList<>(rgbList);
        Collections.sort(gbrList, gbr);
        brgList = new ArrayList<>(rgbList);
        Collections.sort(brgList, brg);
        bgrList = new ArrayList<>(rgbList);
        Collections.sort(bgrList, bgr);

        while (!samples.isEmpty()) {
            Point currentPoint = samples.remove(0);
            int sourceAtPoint = source.getRGB(currentPoint.x, currentPoint.y);
            int bestColor = search(new MInteger(sourceAtPoint));
            result.setRGB(currentPoint.x, currentPoint.y, bestColor);
        }
        return result;
    }

    public static List<Point> getPoints(int width, int height) {
        HashSet<Point> points = new HashSet<>(width * height);
        for (int x = 0; x < width; x++) {
            for (int y = 0; y < height; y++) {
                points.add(new Point(x, y));
            }
        }
        List<Point> newList = new ArrayList<>();
        List<Point> corner1 = new LinkedList<>();
        List<Point> corner2 = new LinkedList<>();
        List<Point> corner3 = new LinkedList<>();
        List<Point> corner4 = new LinkedList<>();

        Point p1 = new Point(width / 3, height / 3);
        Point p2 = new Point(width * 2 / 3, height / 3);
        Point p3 = new Point(width / 3, height * 2 / 3);
        Point p4 = new Point(width * 2 / 3, height * 2 / 3);

        newList.add(p1);
        newList.add(p2);
        newList.add(p3);
        newList.add(p4);
        corner1.add(p1);
        corner2.add(p2);
        corner3.add(p3);
        corner4.add(p4);
        points.remove(p1);
        points.remove(p2);
        points.remove(p3);
        points.remove(p4);

        long seed = System.currentTimeMillis();
        Random c1Random = new Random(seed += 179426549); //The prime number pushes the first numbers apart
        Random c2Random = new Random(seed += 179426549); //Or at least I think it does.
        Random c3Random = new Random(seed += 179426549);
        Random c4Random = new Random(seed += 179426549);

        Dir NW = Dir.NW;
        Dir N = Dir.N;
        Dir NE = Dir.NE;
        Dir W = Dir.W;
        Dir E = Dir.E;
        Dir SW = Dir.SW;
        Dir S = Dir.S;
        Dir SE = Dir.SE;
        while (!points.isEmpty()) {
            putPoints(newList, corner1, c1Random, points, NW, N, NE, W, E, SW, S, SE);
            putPoints(newList, corner2, c2Random, points, NE, N, NW, E, W, SE, S, SW);
            putPoints(newList, corner3, c3Random, points, SW, S, SE, W, E, NW, N, NE);
            putPoints(newList, corner4, c4Random, points, SE, S, SW, E, W, NE, N, NW);
        }
        return newList;
    }

    public static enum Dir {
        NW(-1, -1), N(0, -1), NE(1, -1), W(-1, 0), E(1, 0), SW(-1, 1), S(0, 1), SE(1, 1);
        final int dx, dy;

        private Dir(int dx, int dy) {
            this.dx = dx;
            this.dy = dy;
        }

        public Point add(Point p) {
            return new Point(p.x + dx, p.y + dy);
        }
    }

    public static void putPoints(List<Point> newList, List<Point> listToAddTo, Random rand,
                                 HashSet<Point> points, Dir... adj) {
        List<Point> newPoints = new LinkedList<>();
        for (Iterator<Point> iter = listToAddTo.iterator(); iter.hasNext();) {
            Point p = iter.next();
            Point pul = adj[0].add(p);
            Point pu = adj[1].add(p);
            Point pur = adj[2].add(p);
            Point pl = adj[3].add(p);
            Point pr = adj[4].add(p);
            Point pbl = adj[5].add(p);
            Point pb = adj[6].add(p);
            Point pbr = adj[7].add(p);
            int allChosen = 0;
            if (points.contains(pul)) {
                if (rand.nextInt(5) == 0) {
                    allChosen++;
                    newPoints.add(pul);
                    newList.add(pul);
                    points.remove(pul);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pu)) {
                if (rand.nextInt(5) == 0) {
                    allChosen++;
                    newPoints.add(pu);
                    newList.add(pu);
                    points.remove(pu);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pur)) {
                if (rand.nextInt(3) == 0) {
                    allChosen++;
                    newPoints.add(pur);
                    newList.add(pur);
                    points.remove(pur);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pl)) {
                if (rand.nextInt(5) == 0) {
                    allChosen++;
                    newPoints.add(pl);
                    newList.add(pl);
                    points.remove(pl);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pr)) {
                if (rand.nextInt(2) == 0) {
                    allChosen++;
                    newPoints.add(pr);
                    newList.add(pr);
                    points.remove(pr);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pbl)) {
                if (rand.nextInt(5) == 0) {
                    allChosen++;
                    newPoints.add(pbl);
                    newList.add(pbl);
                    points.remove(pbl);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pb)) {
                if (rand.nextInt(3) == 0) {
                    allChosen++;
                    newPoints.add(pb);
                    newList.add(pb);
                    points.remove(pb);
                }
            } else {
                allChosen++;
            }
            if (points.contains(pbr)) {
                newPoints.add(pbr);
                newList.add(pbr);
                points.remove(pbr);
            }
            if (allChosen == 7) {
                iter.remove();
            }
        }
        listToAddTo.addAll(newPoints);
    }

    public static List<MInteger> getColors(BufferedImage img) {
        int width = img.getWidth();
        int height = img.getHeight();
        List<MInteger> colors = new ArrayList<>(width * height);
        for (int x = 0; x < width; x++) {
            for (int y = 0; y < height; y++) {
                colors.add(new MInteger(img.getRGB(x, y)));
            }
        }
        return colors;
    }

    public static int search(MInteger color) {
        int rgbIndex = binarySearch(rgbList, color, rgb);
        int rbgIndex = binarySearch(rbgList, color, rbg);
        int grbIndex = binarySearch(grbList, color, grb);
        int gbrIndex = binarySearch(gbrList, color, gbr);
        int brgIndex = binarySearch(brgList, color, brg);
        int bgrIndex = binarySearch(bgrList, color, bgr);

        double distRgb = dist(rgbList.get(rgbIndex), color);
        double distRbg = dist(rbgList.get(rbgIndex), color);
        double distGrb = dist(grbList.get(grbIndex), color);
        double distGbr = dist(gbrList.get(gbrIndex), color);
        double distBrg = dist(brgList.get(brgIndex), color);
        double distBgr = dist(bgrList.get(bgrIndex), color);

        double minDist = Math.min(Math.min(Math.min(Math.min(Math.min(
                distRgb, distRbg), distGrb), distGbr), distBrg), distBgr);

        MInteger ans;
        if (minDist == distRgb) {
            ans = rgbList.get(rgbIndex);
        } else if (minDist == distRbg) {
            ans = rbgList.get(rbgIndex);
        } else if (minDist == distGrb) {
            ans = grbList.get(grbIndex);
        } else if (minDist == distGbr) {
            ans = grbList.get(grbIndex);
        } else if (minDist == distBrg) {
            ans = grbList.get(rgbIndex);
        } else {
            ans = grbList.get(grbIndex);
        }
        rgbList.remove(ans);
        rbgList.remove(ans);
        grbList.remove(ans);
        gbrList.remove(ans);
        brgList.remove(ans);
        bgrList.remove(ans);
        return ans.val;
    }

    public static int binarySearch(List<MInteger> list, MInteger val, Comparator<MInteger> cmp){
        int index = Collections.binarySearch(list, val, cmp);
        if (index < 0) {
            index = ~index;
            if (index >= list.size()) {
                index = list.size() - 1;
            }
        }
        return index;
    }

    public static double dist(MInteger color1, MInteger color2) {
        int c1 = color1.val;
        int r1 = (c1 & 0xFF0000) >> 16;
        int g1 = (c1 & 0x00FF00) >> 8;
        int b1 = (c1 & 0x0000FF);

        int c2 = color2.val;
        int r2 = (c2 & 0xFF0000) >> 16;
        int g2 = (c2 & 0x00FF00) >> 8;
        int b2 = (c2 & 0x0000FF);

        int dr = r1 - r2;
        int dg = g1 - g2;
        int db = b1 - b2;
        return Math.sqrt(dr * dr + dg * dg + db * db);
    }

    //This method is here solely for my ease of use (I put the files under <Project Name>/Resources/ )
    public static File resource(String fileName) {
        return new File(System.getProperty("user.dir") + "/Resources/" + fileName);
    }

    static List<MInteger> rgbList;
    static List<MInteger> rbgList;
    static List<MInteger> grbList;
    static List<MInteger> gbrList;
    static List<MInteger> brgList;
    static List<MInteger> bgrList;
    static Comparator<MInteger> rgb = (color1, color2) -> color1.val - color2.val;
    static Comparator<MInteger> rbg = (color1, color2) -> {
        int c1 = color1.val;
        int c2 = color2.val;
        c1 = ((c1 & 0xFF0000)) | ((c1 & 0x00FF00) >> 8) | ((c1 & 0x0000FF) << 8);
        c2 = ((c2 & 0xFF0000)) | ((c2 & 0x00FF00) >> 8) | ((c2 & 0x0000FF) << 8);
        return c1 - c2;
    };
    static Comparator<MInteger> grb = (color1, color2) -> {
        int c1 = color1.val;
        int c2 = color2.val;
        c1 = ((c1 & 0xFF0000) >> 8) | ((c1 & 0x00FF00) << 8) | ((c1 & 0x0000FF));
        c2 = ((c2 & 0xFF0000) >> 8) | ((c2 & 0x00FF00) << 8) | ((c2 & 0x0000FF));
        return c1 - c2;
    };

    static Comparator<MInteger> gbr = (color1, color2) -> {
        int c1 = color1.val;
        int c2 = color2.val;
        c1 = ((c1 & 0xFF0000) >> 16) | ((c1 & 0x00FF00) << 8) | ((c1 & 0x0000FF) << 8);
        c2 = ((c2 & 0xFF0000) >> 16) | ((c2 & 0x00FF00) << 8) | ((c2 & 0x0000FF) << 8);
        return c1 - c2;
    };

    static Comparator<MInteger> brg = (color1, color2) -> {
        int c1 = color1.val;
        int c2 = color2.val;
        c1 = ((c1 & 0xFF0000) >> 8) | ((c1 & 0x00FF00) >> 8) | ((c1 & 0x0000FF) << 16);
        c2 = ((c2 & 0xFF0000) >> 8) | ((c2 & 0x00FF00) >> 8) | ((c2 & 0x0000FF) << 16);
        return c1 - c2;
    };

    static Comparator<MInteger> bgr = (color1, color2) -> {
        int c1 = color1.val;
        int c2 = color2.val;
        c1 = ((c1 & 0xFF0000) >> 16) | ((c1 & 0x00FF00)) | ((c1 & 0x0000FF) << 16);
        c2 = ((c2 & 0xFF0000) >> 16) | ((c2 & 0x00FF00)) | ((c2 & 0x0000FF) << 16);
        return c1 - c2;
    };

    public static void validate(BufferedImage palette, BufferedImage result) {
        List<Integer> paletteColors = getTrueColors(palette);
        List<Integer> resultColors = getTrueColors(result);
        Collections.sort(paletteColors);
        Collections.sort(resultColors);
        System.out.println(paletteColors.equals(resultColors));
    }

    public static List<Integer> getTrueColors(BufferedImage img) {
        int width = img.getWidth();
        int height = img.getHeight();
        List<Integer> colors = new ArrayList<>(width * height);
        for (int x = 0; x < width; x++) {
            for (int y = 0; y < height; y++) {
                colors.add(img.getRGB(x, y));
            }
        }
        Collections.sort(colors);
        return colors;
    }
}

Mi enfoque funciona al encontrar el color más cercano a cada píxel (bueno, probablemente el más cercano), en 3 espacios, ya que los colores son 3D.

Esto funciona creando una lista de todos los puntos que necesitamos llenar y una lista de todos los colores posibles que podemos usar. Aleatorizamos la lista de puntos (para que la imagen resulte mejor), luego revisamos cada punto y obtenemos el color de la imagen de origen.

Actualización: Solía ​​simplemente buscar binarios, por lo que el rojo coincidía mejor que el verde, que coincidía mejor que el azul. Ahora lo cambié para hacer seis búsquedas binarias (todas las permutaciones posibles), luego elijo el color más cercano. Solo toma ~ 6 veces más tiempo (es decir, 1 minuto). Si bien las imágenes aún son granuladas, los colores combinan mejor.

Actualización 2: ya no aleatorizo ​​la lista. En cambio, elijo 4 puntos siguiendo la regla de los tercios, luego organizo los puntos al azar, con preferencia para completar el centro.

Nota: Consulte el historial de revisiones para las imágenes antiguas.

Mona Lisa -> Río:

Mona Lisa -> Gótico Americano:

Mona Lisa -> Esferas Raytraced:

Noche estrellada -> Mona Lisa:

Aquí hay un GIF animado que muestra cómo se construyó la imagen:

Y mostrando los píxeles tomados de la Mona Lisa:

Perl, con espacio de color Lab y tramado

Nota: Ahora tengo una solución C también.

Utiliza un enfoque similar al de aditsu (elija dos posiciones aleatorias e intercambie los píxeles en esas posiciones si haría que la imagen se pareciera más a la imagen objetivo), con dos mejoras importantes:

1. Utiliza el CIE L espacio de color a b * para comparar colores: la métrica euclidiana en este espacio es una muy buena aproximación a la diferencia perceptiva entre dos colores, por lo que las asignaciones de color deberían ser más precisas que RGB o incluso HSV / HSL.
2. Después de un pase inicial que coloca los píxeles en la mejor posición individual posible, realiza un pase adicional con un tramado aleatorio. En lugar de comparar los valores de píxel en las dos posiciones de intercambio, calcula el valor de píxel promedio de un vecindario de 3x3 centrado en las posiciones de intercambio. Si un intercambio mejora los colores promedio de los vecindarios, está permitido, incluso si hace que los píxeles individuales sean menos precisos. Para algunos pares de imágenes, esto tiene un efecto dudoso sobre la calidad (y hace que el efecto de la paleta sea menos llamativo), pero para algunos (como esferas -> cualquier cosa) ayuda bastante. El factor "detalle" enfatiza el píxel central en un grado variable. Al aumentarlo, disminuye la cantidad total de interpolación, pero conserva más detalles finos de la imagen de destino. La optimización difuminada es más lenta,

Promediar valores de laboratorio, como lo hace el tramado, no es realmente justificado (deben convertirse a XYZ, promediar y volver a convertir), pero funciona bien para estos fines.

Estas imágenes tienen límites de terminación de 100 y 200 (finalizan la primera fase cuando se aceptan menos de 1 en 5000 swaps, y la segunda fase en 1 en 2500), y un factor de detalle de tramado de 12 (un tramado más ajustado que el conjunto anterior ) Con esta configuración de súper alta calidad, las imágenes tardan mucho tiempo en generarse, pero con la paralelización, todo el trabajo aún finaliza en una hora en mi caja de 6 núcleos. Al subir los valores hasta 500 o más, termina las imágenes en pocos minutos, solo se ven un poco menos pulidas. Quería mostrar el algoritmo a los mejores aquí.

El código no es de ninguna manera bonito:

#!/usr/bin/perl
use strict;
use warnings;
use Image::Magick;
use Graphics::ColorObject 'RGB_to_Lab';
use List::Util qw(sum max);

my $source = Image::Magick->new;
$source->Read($ARGV[0]);
my $target = Image::Magick->new;
$target->Read($ARGV[1]);
my ($limit1, $limit2, $detail) = @ARGV[2,3,4];

my ($width, $height) = ($target->Get('width'), $target->Get('height'));

# Transfer the pixels of the $source onto a new canvas with the diemnsions of $target
$source->Set(magick => 'RGB');
my $img = Image::Magick->new(size => "${width}x${height}", magick => 'RGB', depth => 8);
$img->BlobToImage($source->ImageToBlob);

my ($made, $rejected) = (0,0);

system("rm anim/*.png");

my (@img_lab, @target_lab);
for my $x (0 .. $width) {
  for my $y (0 .. $height) {
    $img_lab[$x][$y] = RGB_to_Lab([$img->getPixel(x => $x, y => $y)], 'sRGB');
    $target_lab[$x][$y] = RGB_to_Lab([$target->getPixel(x => $x, y => $y)], 'sRGB');
  }
}

my $n = 0;
my $frame = 0;
my $mode = 1;

while (1) {
  $n++;

  my $swap = 0;
  my ($x1, $x2, $y1, $y2) = (int rand $width, int rand $width, int rand $height, int rand $height);
  my ($dist, $dist_swapped);

  if ($mode == 1) {
    $dist = (sum map { ($img_lab[$x1][$y1][$_] - $target_lab[$x1][$y1][$_])**2 } 0..2)
          + (sum map { ($img_lab[$x2][$y2][$_] - $target_lab[$x2][$y2][$_])**2 } 0..2);

    $dist_swapped = (sum map { ($img_lab[$x2][$y2][$_] - $target_lab[$x1][$y1][$_])**2 } 0..2)
                  + (sum map { ($img_lab[$x1][$y1][$_] - $target_lab[$x2][$y2][$_])**2 } 0..2);

  } else { # dither mode
    my $xoffmin = ($x1 == 0 || $x2 == 0 ? 0 : -1);
    my $xoffmax = ($x1 == $width - 1 || $x2 == $width - 1 ? 0 : 1);
    my $yoffmin = ($y1 == 0 || $y2 == 0 ? 0 : -1);
    my $yoffmax = ($y1 == $height - 1 || $y2 == $height - 1 ? 0 : 1);

    my (@img1, @img2, @target1, @target2, $points);
    for my $xoff ($xoffmin .. $xoffmax) {
      for my $yoff ($yoffmin .. $yoffmax) {
        $points++;
        for my $chan (0 .. 2) {
          $img1[$chan] += $img_lab[$x1+$xoff][$y1+$yoff][$chan];
          $img2[$chan] += $img_lab[$x2+$xoff][$y2+$yoff][$chan];
          $target1[$chan] += $target_lab[$x1+$xoff][$y1+$yoff][$chan];
          $target2[$chan] += $target_lab[$x2+$xoff][$y2+$yoff][$chan];
        }
      }
    }

    my @img1s = @img1;
    my @img2s = @img2;
    for my $chan (0 .. 2) {
      $img1[$chan] += $img_lab[$x1][$y1][$chan] * ($detail - 1);
      $img2[$chan] += $img_lab[$x2][$y2][$chan] * ($detail - 1);

      $target1[$chan] += $target_lab[$x1][$y1][$chan] * ($detail - 1);
      $target2[$chan] += $target_lab[$x2][$y2][$chan] * ($detail - 1);

      $img1s[$chan] += $img_lab[$x2][$y2][$chan] * $detail - $img_lab[$x1][$y1][$chan];
      $img2s[$chan] += $img_lab[$x1][$y1][$chan] * $detail - $img_lab[$x2][$y2][$chan];
    }

    $dist = (sum map { ($img1[$_] - $target1[$_])**2 } 0..2)
          + (sum map { ($img2[$_] - $target2[$_])**2 } 0..2);

    $dist_swapped = (sum map { ($img1s[$_] - $target1[$_])**2 } 0..2)
                  + (sum map { ($img2s[$_] - $target2[$_])**2 } 0..2);

  }

  if ($dist_swapped < $dist) {
    my @pix1 = $img->GetPixel(x => $x1, y => $y1);
    my @pix2 = $img->GetPixel(x => $x2, y => $y2);
    $img->SetPixel(x => $x1, y => $y1, color => \@pix2);
    $img->SetPixel(x => $x2, y => $y2, color => \@pix1);
    ($img_lab[$x1][$y1], $img_lab[$x2][$y2]) = ($img_lab[$x2][$y2], $img_lab[$x1][$y1]);
    $made ++;
  } else {
    $rejected ++;
  }

  if ($n % 50000 == 0) {
#    print "Made: $made Rejected: $rejected\n";
    $img->Write('png:out.png');
    system("cp", "out.png", sprintf("anim/frame%05d.png", $frame++));
    if ($mode == 1 and $made < $limit1) {
      $mode = 2;
      system("cp", "out.png", "nodither.png");
    } elsif ($mode == 2 and $made < $limit2) {
      last;
    }
    ($made, $rejected) = (0, 0);
  }
}

Resultados

Paleta gótica americana

Poca diferencia aquí con vacilaciones o no.

Paleta Mona Lisa

El tramado reduce las bandas en las esferas, pero no es especialmente bonito.

Paleta de noche estrellada

Mona Lisa conserva un poco más de detalle con el tramado. Esferas es casi la misma situación que la última vez.

Paleta de gritos

La noche estrellada sin vacilaciones es lo más asombroso de la historia. El tramado hace que sea más fotocorrecto, pero mucho menos interesante.

Paleta de esferas

Como dice aditsu, la verdadera prueba. Creo que paso.

El tramado ayuda inmensamente con American Gothic y Mona Lisa, mezclando algunos grises y otros colores con los píxeles más intensos para producir tonos de piel semi-precisos en lugar de manchas horribles. El grito se ve afectado mucho menos.

Camaro - Mustang

Fuente de imágenes de la publicación de flawr.

Camaro

Mustango:

Paleta Camaro

Se ve bastante bien sin vacilaciones.

Un tramado "apretado" (mismo factor de detalle que el anterior) no cambia mucho, solo agrega un pequeño detalle en los reflejos del capó y el techo.

Un tramado "flojo" (factor de detalle reducido a 6) realmente suaviza la tonalidad, y se ven muchos más detalles a través del parabrisas, pero los patrones de tramado son más obvios en todas partes.

Paleta Mustang

Las partes del automóvil se ven geniales, pero los píxeles grises se ven con fallas. Lo que es peor, todos los píxeles amarillos más oscuros se distribuyeron sobre el cuerpo rojo del Camaro, y el algoritmo de no difuminado no puede encontrar nada que ver con los más claros (moverlos al automóvil empeoraría la combinación y trasladarlos a otro punto en el fondo no hace ninguna diferencia neta), por lo que hay un Mustang fantasma en el fondo.

El tramado es capaz de extender esos píxeles amarillos adicionales para que no se toquen, dispersándolos más o menos uniformemente sobre el fondo en el proceso. Los reflejos y las sombras en el automóvil se ven un poco mejor.

Una vez más, el tramado suelto tiene la tonalidad más grave, revela más detalles sobre los faros y el parabrisas. El auto casi se ve rojo nuevamente. el fondo es más aglomerado por alguna razón. No estoy seguro si me gusta.

Vídeo

( Enlace HQ )

Pitón

La idea es simple: cada píxel tiene un punto en el espacio 3D RGB. El objetivo es hacer coincidir cada píxel de la imagen de origen y uno de la imagen de destino, preferiblemente deben estar 'cerca' (representan el 'mismo' color). Dado que se pueden distribuir de maneras bastante diferentes, no podemos simplemente hacer coincidir el vecino más cercano.

Estrategia

Dejar nser un número entero (pequeño, 3-255 más o menos). Ahora la nube de píxeles en el espacio RGB se ordena por el primer eje (R). Este conjunto de píxeles ahora está dividido en n particiones. Cada una de las particiones ahora se ordena a lo largo del segundo eje (B), que nuevamente se ordena de la misma manera. Hacemos esto con ambas imágenes, y ahora tenemos para ambos una matriz de puntos. Ahora solo podemos hacer coincidir los píxeles por su posición en la matriz, dividir un píxel en la misma posición en cada matriz tiene una posición similar con respecto a cada nube de píxeles en el espacio RGB.

Si la distribución de los píxeles en el espacio RGB de las dos imágenes es similar (significa solo desplazado y / o estirado a lo largo del eje 3), el resultado será bastante predecible. Si las distribuciones se ven completamente diferentes, este algoritmo no producirá resultados tan buenos (como se ve en el último ejemplo), pero este es también uno de los casos más difíciles, creo. Lo que no hace es usar efectos de interacción de píxeles vecinos en la percepción.

Código

Descargo de responsabilidad: soy un novato absoluto en python.

from PIL import Image

n = 5 #number of partitions per channel.

src_index = 3 #index of source image
dst_index = 2 #index of destination image

images =  ["img0.bmp","img1.bmp","img2.bmp","img3.bmp"];
src_handle = Image.open(images[src_index])
dst_handle = Image.open(images[dst_index])
src = src_handle.load()
dst = dst_handle.load()
assert src_handle.size[0]*src_handle.size[1] == dst_handle.size[0]*dst_handle.size[1],"images must be same size"

def makePixelList(img):
    l = []
    for x in range(img.size[0]):
        for y in range(img.size[1]):
            l.append((x,y))
    return l

lsrc = makePixelList(src_handle)
ldst = makePixelList(dst_handle)

def sortAndDivide(coordlist,pixelimage,channel): #core
    global src,dst,n
    retlist = []
    #sort
    coordlist.sort(key=lambda t: pixelimage[t][channel])
    #divide
    partitionLength = int(len(coordlist)/n)
    if partitionLength <= 0:
        partitionLength = 1
    if channel < 2:
        for i in range(0,len(coordlist),partitionLength):
            retlist += sortAndDivide(coordlist[i:i+partitionLength],pixelimage,channel+1)
    else:
        retlist += coordlist
    return retlist

print(src[lsrc[0]])

lsrc = sortAndDivide(lsrc,src,0)
ldst = sortAndDivide(ldst,dst,0)

for i in range(len(ldst)):
    dst[ldst[i]] = src[lsrc[i]]

dst_handle.save("exchange"+str(src_index)+str(dst_index)+".png")

Resultado

Creo que no estuvo mal teniendo en cuenta la solución simple. Por supuesto, puede obtener mejores resultados al jugar con el parámetro, o al transformar primero los colores en otro espacio de color, o incluso al optimizar la partición.

Galería completa aquí: https://imgur.com/a/hzaAm#6

Detallado para River

monalisa> río

gente> río

bolas> río

noche estrellada> río

el llanto> río

bolas> MonaLisa, variando n = 2,4,6, ..., 20

Creo que esta fue la tarea más desafiante, lejos de buenas imágenes, aquí un gif (tuvo que reducirse a 256 colores) de los valores de los parámetros de diferencia n = 2,4,6, ..., 20. Para mí fue sorprendente que valores muy bajos produjeran mejores imágenes (al mirar a la cara de la Sra. Lisa):

Lo siento no puedo parar

¿Cuál te gusta más? ¿Chevy Camaro o Ford Mustang? Quizás esta técnica podría mejorarse y usarse para colorear imágenes bw. Ahora aquí: primero corté los autos del fondo pintando de blanco (en pintura, no muy profesional ...) y luego utilicé el programa Python en cada dirección.

Hay algunos artefactos, creo que porque el área de un automóvil era un poco más grande que el otro y porque mis habilidades artísticas son bastante malas =)

Python: una solución teóricamente óptima

Digo teóricamente óptimo porque la solución verdaderamente óptima no es del todo factible de calcular. Comienzo describiendo la solución teórica, y luego explico cómo la modifiqué para que sea computacionalmente factible tanto en espacio como en tiempo.

Considero la solución más óptima como la que produce el error total más bajo en todos los píxeles entre las imágenes antiguas y las nuevas. El error entre dos píxeles se define como la distancia euclidiana entre los puntos en el espacio 3D donde cada valor de color (R, G, B) es una coordenada. En la práctica, debido a la forma en que los humanos ven las cosas, la solución óptima puede no ser la mejor solución. Sin embargo, parece hacerlo bastante bien en todos los casos.

Para calcular el mapeo, consideré esto como un problema de coincidencia bipartita de peso mínimo . En otras palabras, hay dos conjuntos de nodos: los píxeles originales y los píxeles de la paleta. Se crea un borde entre cada píxel en los dos conjuntos (pero no se crean bordes dentro de un conjunto). El costo, o peso, de un borde es la distancia euclidiana entre los dos píxeles, como se describió anteriormente. Cuanto más cerca estén dos colores visualmente, menor será el costo entre los píxeles.

Esto crea una matriz de costos de tamaño N ² . Para estas imágenes donde N = 123520, se requieren aproximadamente 40 GB de memoria para representar los costos como enteros, y la mitad como enteros cortos. De cualquier manera, no tenía suficiente memoria en mi máquina para intentarlo. Otro problema es que el algoritmo húngaro , o el algoritmo Jonker-Volgenant , que se puede usar para resolver este problema, se ejecuta en tiempo N ³ . Si bien es definitivamente computable, generar una solución por imagen probablemente habría llevado horas o días.

Para solucionar este problema, clasifico aleatoriamente ambas listas de píxeles, divido las listas en fragmentos C, ejecuto una implementación en C ++ del algoritmo Jonker-Volgenant en cada par de sublistas y luego vuelvo a unir las listas para crear el mapeo final. Por lo tanto, el siguiente código permitiría encontrar la solución verdaderamente óptima siempre que establezca el tamaño de fragmento C en 1 (sin fragmentación) y tenga suficiente memoria. Para estas imágenes, configuro C en 16, de modo que N se convierta en 7720, tomando solo unos minutos por imagen.

Una manera simple de pensar por qué esto funciona es que ordenar aleatoriamente la lista de píxeles y luego tomar un subconjunto es como muestrear la imagen. Entonces, al configurar C = 16, es como tomar 16 muestras aleatorias distintas de tamaño N / C tanto del original como de la paleta. De acuerdo, probablemente haya mejores formas de dividir las listas, pero un enfoque aleatorio proporciona resultados decentes.

import subprocess
import multiprocessing as mp
import sys
import os
import sge
from random import shuffle
from PIL import Image
import numpy as np
import LAPJV
import pdb

def getError(p1, p2):
    return (p1[0]-p2[0])**2 + (p1[1]-p2[1])**2 + (p1[2]-p2[2])**2

def getCostMatrix(pallete_list, source_list):
    num_pixels = len(pallete_list)
    matrix = np.zeros((num_pixels, num_pixels))

    for i in range(num_pixels):
        for j in range(num_pixels):
            matrix[i][j] = getError(pallete_list[i], source_list[j])

    return matrix

def chunks(l, n):
    if n < 1:
        n = 1
    return [l[i:i + n] for i in range(0, len(l), n)]

def imageToColorList(img_file):
    i = Image.open(img_file)

    pixels = i.load()
    width, height = i.size

    all_pixels = []
    for x in range(width):
        for y in range(height):
            pixel = pixels[x, y]
            all_pixels.append(pixel)

    return all_pixels

def colorListToImage(color_list, old_img_file, new_img_file, mapping):
    i = Image.open(old_img_file)

    pixels = i.load()
    width, height = i.size
    idx = 0

    for x in range(width):
        for y in range(height):
            pixels[x, y] = color_list[mapping[idx]]
            idx += 1

    i.save(new_img_file)

def getMapping(pallete_list, source_list):
    matrix = getCostMatrix(source_list, pallete_list)
    result = LAPJV.lap(matrix)[1]
    ret = []
    for i in range(len(pallete_list)):
        ret.append(result[i])
    return ret

def randomizeList(l):
    rdm_l = list(l)
    shuffle(rdm_l)
    return rdm_l

def getPartialMapping(zipped_chunk):
    pallete_chunk = zipped_chunk[0]
    source_chunk = zipped_chunk[1]
    subl_pallete = map(lambda v: v[1], pallete_chunk)
    subl_source = map(lambda v: v[1], source_chunk)
    mapping = getMapping(subl_pallete, subl_source)
    return mapping

def getMappingWithPartitions(pallete_list, source_list, C = 1):
    rdm_pallete = randomizeList(enumerate(pallete_list))
    rdm_source = randomizeList(enumerate(source_list))
    num_pixels = len(rdm_pallete)
    real_mapping = [0] * num_pixels

    chunk_size = int(num_pixels / C)

    chunked_rdm_pallete = chunks(rdm_pallete, chunk_size)
    chunked_rdm_source = chunks(rdm_source, chunk_size)
    zipped_chunks = zip(chunked_rdm_pallete, chunked_rdm_source)

    pool = mp.Pool(2)
    mappings = pool.map(getPartialMapping, zipped_chunks)

    for mapping, zipped_chunk in zip(mappings, zipped_chunks):
        pallete_chunk = zipped_chunk[0]
        source_chunk = zipped_chunk[1]
        for idx1,idx2 in enumerate(mapping):
            src_px = source_chunk[idx1]
            pal_px = pallete_chunk[idx2]
            real_mapping[src_px[0]] = pal_px[0]

    return real_mapping

def run(pallete_name, source_name, output_name):
    print("Getting Colors...")
    pallete_list = imageToColorList(pallete_name)
    source_list = imageToColorList(source_name)

    print("Creating Mapping...")
    mapping = getMappingWithPartitions(pallete_list, source_list, C = 16)

    print("Generating Image...");
    colorListToImage(pallete_list, source_name, output_name, mapping)

if __name__ == '__main__':
    pallete_name = sys.argv[1]
    source_name = sys.argv[2]
    output_name = sys.argv[3]
    run(pallete_name, source_name, output_name)

Resultados:

Al igual que con la solución de aditsu, todas estas imágenes se generaron utilizando exactamente los mismos parámetros. El único parámetro aquí es C, que debe establecerse lo más bajo posible. Para mí, C = 16 fue un buen equilibrio entre velocidad y calidad.

Todas las imágenes: http://imgur.com/a/RCZiX#0

Paleta gótica americana

Paleta Mona Lisa

Paleta de noche estrellada

Paleta de gritos

Paleta de río

Paleta de esferas

Pitón

Editar: Acabo de darme cuenta de que en realidad puedes agudizar la fuente con ImageFilter para hacer que los resultados estén mejor definidos.

Rainbow -> Mona Lisa (fuente de Mona Lisa afilada, solo Luminance)

Arco iris -> Mona Lisa (fuente no afilada, ponderada con Y = 10, I = 10, Q = 0)

Mona Lisa -> American Gothic (fuente no afilada, solo Luminance)

Mona Lisa -> American Gothic (fuente no afilada, ponderada con Y = 1, I = 10, Q = 1)

River -> Rainbow (fuente no afilada, solo luminancia)

Básicamente, obtiene todos los píxeles de las dos imágenes en dos listas.

Ordénelos con la luminancia como clave. Y en YIQ representa la luminancia.

Luego, para cada píxel en la fuente (que está en orden ascendente de luminancia) obtenga el valor RGB del píxel del mismo índice en la lista de paletas.

import Image, ImageFilter, colorsys

def getPixels(image):
    width, height = image.size
    pixels = []
    for x in range(width):
        for y in range(height):
            pixels.append([(x,y), image.getpixel((x,y))])
    return pixels

def yiq(pixel):
    # y is the luminance
    y,i,q = colorsys.rgb_to_yiq(pixel[1][0], pixel[1][6], pixel[1][7])
    # Change the weights accordingly to get different results
    return 10*y + 0*i + 0*q

# Open the images
source  = Image.open('ml.jpg')
pallete = Image.open('rainbow.png')

# Sharpen the source... It won't affect the palette anyway =D
source = source.filter(ImageFilter.SHARPEN)

# Sort the two lists by luminance
sourcePixels  = sorted(getPixels(source),  key=yiq)
palletePixels = sorted(getPixels(pallete), key=yiq)

copy = Image.new('RGB', source.size)

# Iterate through all the coordinates of source
# And set the new color
index = 0
for sourcePixel in sourcePixels:
    copy.putpixel(sourcePixel[0], palletePixels[index][8])
    index += 1

# Save the result
copy.save('copy.png')

Para mantenerse al día con la tendencia de las animaciones ...

Los píxeles en el grito se clasifican rápidamente en la noche estrellada y viceversa

C # Winform - Visual Studio 2010

Editar tramado añadido.

Esa es mi versión del algoritmo de intercambio aleatorio: sabor @hobbs. Todavía siento que algún tipo de vacilación no aleatoria puede mejorar ...

Elaboración del color en el espacio Y-Cb-Cr (como en la compresión jpeg)

Elaboración en dos fases:

1. Copia del píxel de la fuente en orden de luminancia. Esto ya da una buena imagen, pero desaturada, casi escala de grises, en un tiempo cercano a 0
2. Intercambio aleatorio repetido de píxeles. El intercambio se realiza si esto da un mejor delta (respecto a la fuente) en la celda de 3x3 que contiene el píxel. Entonces es un efecto de vacilación. El delta se calcula en el espacio Y-Cr-Cb sin ponderar los diferentes componentes.

Este es esencialmente el mismo método utilizado por @hobbs, sin el primer intercambio aleatorio sin vacilar. Simplemente, mis tiempos son más cortos (¿el idioma cuenta?) Y creo que mis imágenes son mejores (probablemente el espacio de color utilizado sea más preciso).

Uso del programa: coloque imágenes .png en su carpeta c: \ temp, verifique el elemento en la lista para elegir la imagen de la paleta, seleccione el elemento en la lista para elegir la imagen de origen (no tan fácil de usar). Haga clic en el botón de inicio para comenzar la elaboración, el guardado es automático (incluso cuando prefiera no hacerlo, tenga cuidado).

Tiempo de elaboración inferior a 90 segundos.

Resultados actualizados

Paleta: gótica americana

Paleta: Monna Lisa

Paleta: arcoiris

Paleta: Río

Paleta: Grito

Paleta: Noche estrellada

Form1.cs

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
using System.Drawing.Imaging;
using System.IO;

namespace Palette
{
    public struct YRB
    {
        public int y, cb, cr;

        public YRB(int r, int g, int b)
        {
            y = (int)(0.299 * r + 0.587 * g + 0.114 * b);
            cb = (int)(128 - 0.168736 * r - 0.331264 * g + 0.5 * b);
            cr = (int)(128 + 0.5 * r - 0.418688 * g - 0.081312 * b);
        }
    }

    public struct Pixel
    {
        private const int ARGBAlphaShift = 24;
        private const int ARGBRedShift = 16;
        private const int ARGBGreenShift = 8;
        private const int ARGBBlueShift = 0;

        public int px, py;
        private uint _color;
        public YRB yrb;

        public Pixel(uint col, int px = 0, int py = 0)
        {
            this.px = px;
            this.py = py;
            this._color = col;
            yrb = new YRB((int)(col >> ARGBRedShift) & 255, (int)(col >> ARGBGreenShift) & 255, (int)(col >> ARGBBlueShift) & 255); 
        }

        public uint color
        {
            get { 
                return _color; 
            }
            set {
                _color = color;
                yrb = new YRB((int)(color >> ARGBRedShift) & 255, (int)(color >> ARGBGreenShift) & 255, (int)(color >> ARGBBlueShift) & 255);
            }
        }

        public int y
        {
            get { return yrb.y; }
        }
        public int cr
        {
            get { return yrb.cr; }
        }
        public int cb
        {
            get { return yrb.cb; }
        }
    }

    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }

        private void Form1_Load(object sender, EventArgs e)
        {
            DirectoryInfo di = new System.IO.DirectoryInfo(@"c:\temp\");
            foreach (FileInfo file in di.GetFiles("*.png"))
            {
                ListViewItem item = new ListViewItem(file.Name);
                item.SubItems.Add(file.FullName);
                lvFiles.Items.Add(item);
            }
        }

        private void lvFiles_ItemSelectionChanged(object sender, ListViewItemSelectionChangedEventArgs e)
        {
            if (e.IsSelected)
            {
                string file = e.Item.SubItems[1].Text;
                GetImagePB(pbSource, file);
                pbSource.Tag = file; 
                DupImage(pbSource, pbOutput);

                this.Width = pbOutput.Width + pbOutput.Left + 20;
                this.Height = Math.Max(pbOutput.Height, pbPalette.Height)+lvFiles.Height*2;   
            }
        }

        private void lvFiles_ItemCheck(object sender, ItemCheckEventArgs e)
        {
            foreach (ListViewItem item in lvFiles.CheckedItems)
            {
                if (item.Index != e.Index) item.Checked = false;
            }
            string file = lvFiles.Items[e.Index].SubItems[1].Text;
            GetImagePB(pbPalette, file);
            pbPalette.Tag = lvFiles.Items[e.Index].SubItems[0].Text; 

            this.Width = pbOutput.Width + pbOutput.Left + 20;
            this.Height = Math.Max(pbOutput.Height, pbPalette.Height) + lvFiles.Height * 2;   
        }

        Pixel[] Palette;
        Pixel[] Source;

        private void BtnStart_Click(object sender, EventArgs e)
        {
            lvFiles.Enabled = false;
            btnStart.Visible = false;
            progressBar.Visible = true; 
            DupImage(pbSource, pbOutput);

            Work(pbSource.Image as Bitmap, pbPalette.Image as Bitmap, pbOutput.Image as Bitmap);

            string newfile = (string)pbSource.Tag +"-"+ (string)pbPalette.Tag;
            pbOutput.Image.Save(newfile, ImageFormat.Png);   

            lvFiles.Enabled = true;
            btnStart.Visible = true;
            progressBar.Visible = false;
        }

        private void Work(Bitmap srcb, Bitmap palb, Bitmap outb)
        {
            GetData(srcb, out Source);
            GetData(palb, out Palette);

            FastBitmap fout = new FastBitmap(outb);
            FastBitmap fsrc = new FastBitmap(srcb);
            int pm = Source.Length;
            int w = outb.Width;
            int h = outb.Height;
            progressBar.Maximum = pm;

            fout.LockImage();
            for (int p = 0; p < pm; p++)
            {
                fout.SetPixel(Source[p].px, Source[p].py, Palette[p].color);
            }
            fout.UnlockImage();

            pbOutput.Refresh();

            var rnd = new Random();
            int totsw = 0;
            progressBar.Maximum = 200;
            for (int i = 0; i < 200; i++)
            {
                int nsw = 0;
                progressBar.Value = i;
                fout.LockImage();
                fsrc.LockImage();
                for (int j = 0; j < 200000; j++)
                {
                    nsw += CheckSwap(fsrc, fout, 1 + rnd.Next(w - 2), 1 + rnd.Next(h - 2), 1 + rnd.Next(w - 2), 1 + rnd.Next(h - 2));
                }
                totsw += nsw;
                lnCurSwap.Text = nsw.ToString();
                lnTotSwap.Text = totsw.ToString();
                fout.UnlockImage();
                fsrc.UnlockImage();
                pbOutput.Refresh();
                Application.DoEvents();
                if (nsw == 0)
                {
                    break;
                }
            }            
        }

        int CheckSwap(FastBitmap fsrc, FastBitmap fout, int x1, int y1, int x2, int y2)
        {
            const int fmax = 3;
            YRB ov1 = new YRB();
            YRB sv1 = new YRB();
            YRB ov2 = new YRB();
            YRB sv2 = new YRB();

            int f;
            for (int dx = -1; dx <= 1; dx++)
            {
                for (int dy = -1; dy <= 1; dy++)
                {
                    f = (fmax - Math.Abs(dx) - Math.Abs(dy));
                    {
                        Pixel o1 = new Pixel(fout.GetPixel(x1 + dx, y1 + dy));
                        ov1.y += o1.y * f;
                        ov1.cb += o1.cr * f;
                        ov1.cr += o1.cb * f;

                        Pixel s1 = new Pixel(fsrc.GetPixel(x1 + dx, y1 + dy));
                        sv1.y += s1.y * f;
                        sv1.cb += s1.cr * f;
                        sv1.cr += s1.cb * f;

                        Pixel o2 = new Pixel(fout.GetPixel(x2 + dx, y2 + dy));
                        ov2.y += o2.y * f;
                        ov2.cb += o2.cr * f;
                        ov2.cr += o2.cb * f;

                        Pixel s2 = new Pixel(fsrc.GetPixel(x2 + dx, y2 + dy));
                        sv2.y += s2.y * f;
                        sv2.cb += s2.cr * f;
                        sv2.cr += s2.cb * f;
                    }
                }
            }
            YRB ox1 = ov1;
            YRB ox2 = ov2;
            Pixel oc1 = new Pixel(fout.GetPixel(x1, y1));
            Pixel oc2 = new Pixel(fout.GetPixel(x2, y2));
            ox1.y += fmax * oc2.y - fmax * oc1.y;
            ox1.cb += fmax * oc2.cr - fmax * oc1.cr;
            ox1.cr += fmax * oc2.cb - fmax * oc1.cb;
            ox2.y += fmax * oc1.y - fmax * oc2.y;
            ox2.cb += fmax  * oc1.cr - fmax * oc2.cr;
            ox2.cr += fmax * oc1.cb - fmax * oc2.cb;

            int curd = Delta(ov1, sv1, 1) + Delta(ov2, sv2, 1);
            int newd = Delta(ox1, sv1, 1) + Delta(ox2, sv2, 1);
            if (newd < curd)
            {
                fout.SetPixel(x1, y1, oc2.color);
                fout.SetPixel(x2, y2, oc1.color);
                return 1;
            }
            return 0;
        }

        int Delta(YRB p1, YRB p2, int sf)
        {
            int dy = (p1.y - p2.y);
            int dr = (p1.cr - p2.cr);
            int db = (p1.cb - p2.cb);

            return dy * dy * sf + dr * dr + db * db;
        }

        Bitmap GetData(Bitmap bmp, out Pixel[] Output)
        {
            FastBitmap fb = new FastBitmap(bmp);
            BitmapData bmpData = fb.LockImage(); 

            Output = new Pixel[bmp.Width * bmp.Height];

            int p = 0;
            for (int y = 0; y < bmp.Height; y++)
            {
                uint col = fb.GetPixel(0, y);
                Output[p++] = new Pixel(col, 0, y);

                for (int x = 1; x < bmp.Width; x++)
                {
                    col = fb.GetNextPixel();
                    Output[p++] = new Pixel(col, x, y);
                }
            }
            fb.UnlockImage(); // Unlock the bits.

            Array.Sort(Output, (a, b) => a.y - b.y);

            return bmp;
        }

        void DupImage(PictureBox s, PictureBox d)
        {
            if (d.Image != null)
                d.Image.Dispose();
            d.Image = new Bitmap(s.Image.Width, s.Image.Height);  
        }

        void GetImagePB(PictureBox pb, string file)
        {
            Bitmap bms = new Bitmap(file, false);
            Bitmap bmp = bms.Clone(new Rectangle(0, 0, bms.Width, bms.Height), PixelFormat.Format32bppArgb);
            bms.Dispose(); 
            if (pb.Image != null)
                pb.Image.Dispose();
            pb.Image = bmp;
        }
    }

    //Adapted from Visual C# Kicks - http://www.vcskicks.com/
    unsafe public class FastBitmap
    {
        private Bitmap workingBitmap = null;
        private int width = 0;
        private BitmapData bitmapData = null;
        private Byte* pBase = null;

        public FastBitmap(Bitmap inputBitmap)
        {
            workingBitmap = inputBitmap;
        }

        public BitmapData LockImage()
        {
            Rectangle bounds = new Rectangle(Point.Empty, workingBitmap.Size);

            width = (int)(bounds.Width * 4 + 3) & ~3;

            //Lock Image
            bitmapData = workingBitmap.LockBits(bounds, ImageLockMode.ReadWrite, PixelFormat.Format32bppArgb);
            pBase = (Byte*)bitmapData.Scan0.ToPointer();
            return bitmapData;
        }

        private uint* pixelData = null;

        public uint GetPixel(int x, int y)
        {
            pixelData = (uint*)(pBase + y * width + x * 4);
            return *pixelData;
        }

        public uint GetNextPixel()
        {
            return *++pixelData;
        }

        public void GetPixelArray(int x, int y, uint[] Values, int offset, int count)
        {
            pixelData = (uint*)(pBase + y * width + x * 4);
            while (count-- > 0)
            {
                Values[offset++] = *pixelData++;
            }
        }

        public void SetPixel(int x, int y, uint color)
        {
            pixelData = (uint*)(pBase + y * width + x * 4);
            *pixelData = color;
        }

        public void SetNextPixel(uint color)
        {
            *++pixelData = color;
        }

        public void UnlockImage()
        {
            workingBitmap.UnlockBits(bitmapData);
            bitmapData = null;
            pBase = null;
        }
    }

}