Services & ResourcesWolfram ForumsMathGroup Archive

[Date Index] [Thread Index] [Author Index]

Re: Memory when reding in data

• To: mathgroup at smc.vnet.net
• Subject: [mg120990] Re: Memory when reding in data
• From: "Scot T. Martin" <smartin at seas.harvard.edu>
• Date: Sat, 20 Aug 2011 06:18:31 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• References: <201108191034.GAA24015@smc.vnet.net>

When you don't need Mathematica to remember anything for In[..] and Out[..] functionality, set $HistoryLength, e.g., default value is Infinity. If you really don't need the In/Out functionality, use $HistoryLength = 1.
________________________________________
From: Peter Rodenbach [peter.rodenbach at googlemail.com]
Sent: Friday, August 19, 2011 06:34
To: mathgroup at smc.vnet.net
Subject: [mg120990] Memory when reding in data

Hi everyone.

I am analysing image data and I have the problem that after about 50
images or so Mathematica quits the calculation with the hint out of
memory. The interesting part of the code i wrote is :

(*Import the first image into Maxvals, which will be compared to all \
other images and overwritten if any point is found to be larger*)

SetDirectory[
  "E:\\...];

Initialise = OpenRead["s2_260_0000.edf", BinaryFormat -> True];
MaxValues =
  BinaryRead["s2_260_0000.edf",
   Table["UnsignedInteger16", {512}, {512}]];
counter = Dimensions[MaxValues];

(*Load all experimental images and perform maximum projection*)

MaxProjection[filename_] := {OpenRead[filename, BinaryFormat -> True],
   Imagevalues =
   BinaryRead[filename, Table["UnsignedInteger16", {512}, {512}]],
  For[i = 1, i <= counter[[1]], i++,
        For [j = 1, j <= counter[[2]], j++,
    If[MaxValues[[i, j]] >= Imagevalues[[i, j]],
     MaxValues[[i, j]] = MaxValues[[i, j]],
     MaxValues[[i, j]] = Imagevalues[[i, j]]]
        ]
   ], Close[Imagevalues]}


MaxProjection[#] & /@ FileNames@"*edf"

So I guess Mathm saves all the data of all the images after each
other. But once the comparision is done with the one MaxProjection
Image I want to have, all the data can be dumped from the memory and
then the next image can be compared with this MaxProjection image.

Does anyone have an idea? Thx.

• References:
  • Memory when reding in data
    • From: Peter Rodenbach <peter.rodenbach@googlemail.com>