# -*- coding: utf-8 -*-
"""
Created on Sun May 17 17:07:43 2015

@author: JC
"""

import os
from PIL import Image
from pylab import *

path = r'Your path here'
jpglist = os.listdir(path)

outfile = os.path.join(path,'Milow_Star_Trail.jpg')
# check size of one of the photos first
om = zeros((683, 1024, 3),dtype=uint8)

for jpg in [os.path.join(path,f) for f in os.listdir(path) if f.endswith('Star_Field.JPG')]:
    pil_im = Image.open(jpg)
    im = array(pil_im)

    for row in range(im.shape[0]):
        if row%250 == 0:
            print row
        for col in range(im.shape[1]):
            r,g,b = im[row,col,:]
            r0,g0,b0 = om[row,col,:]

            if long(r)*long(r)+long(g)*long(g)+long(b)*long(b) > long(r0)*long(r0)+long(g0)*long(g0)+long(b0)*long(b0):
                om[row,col,:] = array([r,g,b])
    print jpg
    print im.shape


Image.fromarray(om).save(outfile)