[Proj] Convert the array into required coordinates

Cleo Drakos cleo21drakos at gmail.com
Sat Aug 23 01:38:01 EST 2014


Dear guys,

I have 2d numpy array of 480 rows and 1440 columns as named by 'data' below:

The first element belongs to (49.875S,179.875W),
the second element belongs to (49.625S,179.625W),and the last element
belongs to (49.875N,179.875E).
import os, glob, gdal, numpy as np

fname = '3B42RT.2014010606.7.bin'
with open(fname, 'rb') as fi:
    fi.seek(2880,0)
    data = np.fromfile(fi,dtype=np.uint16,count=480*1440)
    data = data.byteswap()
    data = data.reshape(1440,480)

 How can I convert this numpy array so that its first element belongs to
(49.875N,179.625W), i.e., upper left latitude and longitude respectively;
and the last element belong to (49.625S,179.875E), i.e., lower right
latitute and longitude respectively.

I tried to rotate it, but I do not think it is correct.

 data = np.rot90(data,1)

 Have some of you experienced with this type of problem? The binary file I
am using is here:
ftp://trmmopen.gsfc.nasa.gov/pub/merged/3B42RT/3B42RT.2014010606.7.bin.gz
cleo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20140823/01e995c2/attachment.htm 


More information about the Proj mailing list