Path: csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!news.albasani.net!newsfeed.freenet.ag!news2.euro.net!newsfeed.xs4all.nl!newsfeed1.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail
MIME-Version: 1.0
In-Reply-To: <c4a4e4ee-11c9-4b5e-a24a-10be49a4ce5d@googlegroups.com>
References: <c4a4e4ee-11c9-4b5e-a24a-10be49a4ce5d@googlegroups.com>
From: Joshua Landau <joshua.landau.ws@gmail.com>
Date: Mon, 6 May 2013 17:28:34 +0100
Subject: Re: Help with loading file into an array
To: peter berrett <pwberrett@gmail.com>
Content-Type: multipart/alternative; boundary=001a11c2a9b89207fd04dc0f33f7
Cc: python-list <python-list@python.org>
Precedence: list
Newsgroups: comp.lang.python
Message-ID: <mailman.1366.1367857761.3114.python-list@python.org>
Lines: 164
NNTP-Posting-Host: 2001:888:2000:d::a6
Xref: csiph.com comp.lang.python:44836

--001a11c2a9b89207fd04dc0f33f7
Content-Type: text/plain; charset=ISO-8859-1

On 5 May 2013 07:06, peter berrett <pwberrett@gmail.com> wrote:

> I am trying to build a program that can find comets in a series of
> astronomical images. I have already written functions to find the comet in
> a series of images, the data of which is stored in embedded lists.
>
> The area I am having difficulty with is take a standard gif file (1024 x
> 1024) and reading it into an array or embedded lists.
>

This is not what you should really be doing. There are specialised things
for this.

Take a look at PIL (http://www.pythonware.com/products/pil/).
It will install as "Image" and from there you can do:

import Image
image = Image.open("FILE")
pixels = image.load()

and access the pixels like:

pixels[x, y] # You get (r, g, b) tuple of integers from 0 to 255


> In a nutshell here is an example of what I want to do
>
> Let's say I have a gif file called 20130428_0000_c2_1024.gif in a folder
> called c:\comets
>
> I want to read the data from that gif file taking the red data (excluding
> the green and blue data) and store that in an array called Image[][] which
> is a nested array length 1024 with a list in each item of 1024 length (ie
> 1024 x 1024)
>

*Ahem*:
1) Python standard is to only uppercase class (and sometimes module) names,
so it should just be "image".
2) Technically, you mean list.
3) Chose the above method instead, so instead of indexing it with
"image[x][y]" you use "image[x, y]".

Now, you only want the red channel. What this means is that you want to
"split" the image into channels before using image.load():

import Image

image = Image.open("FILE")
image.load() # .open is lazy, so you have to load now

red, green, blue = image.split()
red_pixels = red.load()

And you access as before:

red_pixels[x, y] # You get integer from 0 to 255

Could someone please provide a piece of code to do the above so I can then
> go on to modify it to pick up different files from different folders? In
> particular I am keen to seen how you read in the data and also how you
> change the directory from which you are reading the image.
>

Changing directories is done with "os.chdir" (
http://docs.python.org/3/library/os.html#os.chdir).


> For the record this is not for homework but is a pet project of mine. I
> have already written a version of the program in Justbasic but python is
> faster. I am also interested in readers views as to which is the simplest
> and best way to achieve what I am trying to do.
>

Thanks for the clarification, btw. Good luck, as well.

--001a11c2a9b89207fd04dc0f33f7
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">On 5 May 2013 07:06, peter berrett <span dir=3D"ltr">&lt;<=
a href=3D"mailto:pwberrett@gmail.com" target=3D"_blank">pwberrett@gmail.com=
</a>&gt;</span> wrote:<br><div class=3D"gmail_extra"><div class=3D"gmail_qu=
ote"><blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px 0.8ex;bo=
rder-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:so=
lid;padding-left:1ex">

I am trying to build a program that can find comets in a series of astronom=
ical images. I have already written functions to find the comet in a series=
 of images, the data of which is stored in embedded lists.<br>
<br>
The area I am having difficulty with is take a standard gif file (1024 x 10=
24) and reading it into an array or embedded lists.<br></blockquote><div><b=
r></div><div>This is not what you should really be doing. There are special=
ised things for this.</div>

<div><br></div><div>Take a look at PIL (<a href=3D"http://www.pythonware.co=
m/products/pil/">http://www.pythonware.com/products/pil/</a>).</div><div>It=
 will install as &quot;Image&quot; and from there you can do:</div><div>
<br>
</div><div>import Image</div><div>image =3D Image.open(&quot;FILE&quot;)</d=
iv><div>pixels =3D image.load()</div><div><br></div><div>and access the pix=
els like:</div><div><br></div><div>pixels[x, y] # You get (r, g, b) tuple o=
f integers from 0 to 255</div>

<div>=A0</div><blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px=
 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left=
-style:solid;padding-left:1ex">
In a nutshell here is an example of what I want to do<br>
<br>
Let&#39;s say I have a gif file called 20130428_0000_c2_1024.gif in a folde=
r called c:\comets<br>
<br>
I want to read the data from that gif file taking the red data (excluding t=
he green and blue data) and store that in an array called Image[][] which i=
s a nested array length 1024 with a list in each item of 1024 length (ie 10=
24 x 1024)<br>

</blockquote><div>=A0</div><div>*Ahem*:</div><div>1) Python standard is to =
only uppercase class (and sometimes module) names, so it should just be &qu=
ot;image&quot;.</div><div>2) Technically, you mean list.</div><div>3) Chose=
 the above method instead, so instead of indexing it with &quot;image[x][y]=
&quot; you use &quot;image[x, y]&quot;.</div>

<div><br></div><div>Now, you only want the red channel. What this means is =
that you want to &quot;split&quot; the image into channels before using ima=
ge.load():</div><div><br></div><div><div>import Image</div><div><br></div>

<div>image =3D Image.open(&quot;FILE&quot;)</div><div>image.load() # .open =
is lazy, so you have to load now</div></div><div><br></div><div>red, green,=
 blue =3D image.split()</div><div>red_pixels =3D red.load()</div><div><br><=
/div>

<div>And you access as before:</div><div><br></div><div>red_pixels[x, y] # =
You get integer from 0 to 255</div><div><br></div><blockquote class=3D"gmai=
l_quote" style=3D"margin:0px 0px 0px 0.8ex;border-left-width:1px;border-lef=
t-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">


Could someone please provide a piece of code to do the above so I can then =
go on to modify it to pick up different files from different folders? In pa=
rticular I am keen to seen how you read in the data and also how you change=
 the directory from which you are reading the image.<br>

</blockquote><div><br></div><div>Changing directories is done with &quot;os=
.chdir&quot; (<a href=3D"http://docs.python.org/3/library/os.html#os.chdir"=
>http://docs.python.org/3/library/os.html#os.chdir</a>).</div><div>=A0</div=
>

<blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px 0.8ex;border-=
left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;p=
adding-left:1ex">
For the record this is not for homework but is a pet project of mine. I hav=
e already written a version of the program in Justbasic but python is faste=
r. I am also interested in readers views as to which is the simplest and be=
st way to achieve what I am trying to do.<br>

</blockquote><div><br></div><div>Thanks for the clarification, btw. Good lu=
ck, as well.=A0</div></div></div></div>

--001a11c2a9b89207fd04dc0f33f7--