Path: csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!feeder.erje.net!eu.feeder.erje.net!newsfeed.xs4all.nl!newsfeed2a.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail Return-Path: X-Original-To: python-list@python.org Delivered-To: python-list@mail.python.org X-Spam-Status: OK 0.001 X-Spam-Evidence: '*H*': 1.00; '*S*': 0.00; 'algorithm': 0.04; 'assuming': 0.09; 'bytes,': 0.09; 'bytes.': 0.09; 'received:80.91': 0.09; 'received:80.91.229': 0.09; 'received:gmane.org': 0.09; 'received:list': 0.09; 'thus,': 0.09; 'assume': 0.14; '"extended': 0.16; '"python': 0.16; 'int)': 0.16; 'numpy': 0.16; 'received:80.91.229.3': 0.16; 'received:dip0.t-ipconnect.de': 0.16; 'received:plane.gmane.org': 0.16; 'received:t-ipconnect.de': 0.16; 'subject:CSV': 0.16; 'do,': 0.16; 'wrote:': 0.18; 'bit': 0.19; 'seems': 0.21; '>>>': 0.22; 'header:User-Agent:1': 0.23; 'sorry,': 0.24; 'mon,': 0.24; 'looks': 0.24; 'mention': 0.26; 'values': 0.27; 'header:X -Complaints-To:1': 0.27; 'point': 0.28; "doesn't": 0.30; 'sep': 0.31; 'raw': 0.33; "can't": 0.35; 'something': 0.35; 'but': 0.35; '+0200,': 0.36; 'possible': 0.36; 'changing': 0.37; 'too': 0.37; 'to:addr:python-list': 0.38; 'specialist': 0.39; 'to:addr:python.org': 0.39; 'received:org': 0.40; 'skip:u 10': 0.60; 'read': 0.60; 'extended': 0.61; 'skip:n 10': 0.64; 'account': 0.65; 'otten': 0.84; 'skip:] 10': 0.84; 'generation.': 0.91; 'try.': 0.91 X-Injected-Via-Gmane: http://gmane.org/ To: python-list@python.org From: Peter Otten <__peter__@web.de> Subject: Re: CSV methodology Date: Tue, 16 Sep 2014 13:22:02 +0200 Organization: None References: Mime-Version: 1.0 Content-Type: text/plain; charset="ISO-8859-1" Content-Transfer-Encoding: 7Bit X-Gmane-NNTP-Posting-Host: p57bd8bbb.dip0.t-ipconnect.de User-Agent: KNode/4.13.3 X-BeenThere: python-list@python.org X-Mailman-Version: 2.1.15 Precedence: list List-Id: General discussion list for the Python programming language List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Newsgroups: comp.lang.python Message-ID: Lines: 49 NNTP-Posting-Host: 2001:888:2000:d::a6 X-Trace: 1410866541 news.xs4all.nl 2861 [2001:888:2000:d::a6]:48978 X-Complaints-To: abuse@xs4all.nl Xref: csiph.com comp.lang.python:77927 jayte wrote: > On Mon, 15 Sep 2014 09:29:02 +0200, Peter Otten <__peter__@web.de> wrote: > >>jayte wrote: >> >>> Sorry, I neglected to mention the values' significance. The MXP program >>> uses the "distance estimate" algorithm in its fractal data generation. >>> The values are thus, for each point in a 1778 x 1000 image: >>> >>> Distance, (an extended double) >>> Iterations, (a 16 bit int) >>> zc_x, (a 16 bit int) >>> zc_y (a 16 bit int) >>> >> >>Probably a bit too early in your "Python career", > > Absolutely, just thought it would be interesting to start experimenting, > while learning (plus, can't help but be anxious) > >> but you can read raw data >>with numpy. Something like >> >>with open(filename, "rb") as f: >> a = numpy.fromfile(f, dtype=[ >> ("distance", "f16"), >> ("iterations", "i2"), >> ("zc_x", "i2"), >> ("zc_y", "i2"), >> ]).reshape(1778, 1000) >> >>might do, assuming "extended double" takes 16 bytes. > > Will try. Double extended precision is ten bytes, but I assume > changing the "f16" to "f10" would account for that... Unfortunately it seems that numpy doesn't support "f10" >>> numpy.dtype("f8") dtype('float64') >>> numpy.dtype("f16") dtype('float128') >>> numpy.dtype("f10") dtype('float32') # looks strange to me But you better ask for confirmation (and possible workarounds) in a specialist forum.