Challenge to convert a simple IDL code into Python
Gary Herron
gary.herron at islandtraining.com
Sun Aug 24 13:40:27 EDT 2014
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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/python-list/attachments/20140824/724306d0/attachment.html>
More information about the Python-list
mailing list