Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.lang.python > #76934

Re: Challenge to convert a simple IDL code into Python

Show all headers | View raw

[Multipart message — attachments visible in raw view] - view raw

On 08/24/2014 08:38 AM, Cleo Drakos wrote:
>
> Here is IDL code:
>
> |pro read_binary_file
>
> file=  "3B42RT.2014010318.7.bin"
>
> num_lon=  1440
> num_lat=  480
>
> data=  {header:  bytarr(num_lon*2),  precip:  intarr(num_lon,num_lat),  precip_error:  intarr(num_lon,num_lat),  $
>        source_of_estimate:  bytarr(num_lon,num_lat),  precip_uncal:  intarr(num_lon,num_lat)}
>
> close,  1
> openr,  1,  file
> readu,  1,  data
> close,  1     
> precip=  swap_endian(data.precip)
>
> print,  precip
> end|
>
> My attempt in Python is here:
>
> |fname=  '3B42RT.2014010318.7.bin'
>
> num_lon=  1440
> num_lat=  480
>
> with  open(fname,  'rb')  as  fi:
>      contents=  fi.read()
>      precip=  struct.unpack_from('>'+str(num_lon*num_lat)+'I',  contents,  offset=  num_lon*2)
>
>      precip=  np.array(data,dtype=int)
>      precip data|
>
> But, the results obtained from two codes above are not same. The IDL 
> code is correct, but what is wrong with my python code?
>

You are not giving us much to go on here.

First that's not legal Python .  The |last line was probably meant to be 
a print|||, but if you expect us to examine your code, you absolutely 
*must* give us the actual code.  Please cut and paste.

Second, saying it doesn't work us useless.  Please tell us what the code 
is supposed to do, and what it actually does, what you expected...  If 
there is output, cut and paste it.  If there is an error traceback, cut 
and paste it.

While you are at it, tell us what the IDL code does and cut and paste 
it's output.

Gary Herron


> Here is binary file I was working 
> with:ftp://trmmopen.gsfc.nasa.gov/pub/merged/3B42RT/3B42RT.2014010318.7.bin.gz
>
>
> cleo
>
>

Back to comp.lang.python | Previous | Next | Find similar | Unroll thread

Thread

Re: Challenge to convert a simple IDL code into Python Gary Herron <gary.herron@islandtraining.com> - 2014-08-24 10:40 -0700

csiph-web