Hi Sandrine,

Thanks for the info! That sounds like a good approach.

I’ll read a bit more into ZPL. I think your macro uses most of the functions and syntax I would need to learn to make the Bayer code.

It will probably end up being very similar to this Matlab code I just wrote:

%import image
colorimg=imread('park.jpg');
dim=size(colorimg);
xpix=dim(2);
ypix=dim(1);

%These are placeholders for now, I would update them for color filter transmission + QE estimates

rimg=colorimg(:,:,1)*0.2126; 
gimg=colorimg(:,:,2)*0.7152;
bimg=colorimg(:,:,3)*0.0722;

%rggb
for n=1:ypix
    for m=1:xpix
        
        if mod(n,2)==1
            
                    %r
                    if mod(m,2)==1;
                    raw_mosaic(n,m)=rimg(n,m);

                    %g
                    else
                    raw_mosaic(n,m)=gimg(n,m);
                    end
        else
                    %g
                    if mod(m,2)==1;
                    raw_mosaic(n,m)=gimg(n,m);

                    %b
                    else
                    raw_mosaic(n,m)=bimg(n,m);
                    end
                    
        end
        
    end
end

subplot(1,2,1)
imagesc(colorimg), axis tight, axis equal, axis off, title('Unfiltered Color Image')
subplot(1,2,2)
imagesc(raw_mosaic), axis tight, axis equal, axis off, colormap gray, title('Bayer Filter Sim')