| satellite data [message #35099] |
Tue, 13 May 2003 13:17  |
mmccabe
Messages: 4
Registered: May 2003
|
Junior Member |
|
|
Hi all.
I know this is an old question, but I'm a new poster so please be
gentle.
I muck around with all sorts of satellite data and am interested in
finding a nice simple method of plotting said data onto a map
projection or regular grid wihtout using ENVI. Some of the data I use
(TRMM-TMI) is irregular in that I just have lat/lon and the associated
data values. Others, such as when I extract info from TERRA or AQUA
files, is regular data, but only has the bounding co-ords of the
image.
Can someone, in the spirit of 'shepherding the weak through the valley
of darkness' offer me some basic tips on addressing this. It would be
nice if I didn't have to interpolate the data also - so that I could
interrogate the image for actual information rather than just
displaying it.
While I'm at it, any tips on efficiently reading binary data of
unknown size into an array - the data consists of 11 values (4byte
records) for each pixel, but you never know how many pixels there are
going to be for a particular overpass.
What I do at the moment is:
...
openr, lun, file, /get_lun, /swap_if_little_endian
fileinfo = FStat(lun)
rows = fileinfo.size/(11*4)
tmp=fltarr(11,rows)
while not eof(lun) do begin
readu,lun,tmp
endwhile
....
Then I just need to extract the particular columns I'm interested in
and make them into a 2D Array for plotting purposes (associated with
lat/lon) but I'm not sure whether this is the best way to do it...hmmm
- this probably doesn't make much sense so if your interested I could
email you with some more detailed information.
Thanks for your patience!
|
|
|
|
| Re: satellite data [message #75556 is a reply to message #35099] |
Thu, 17 March 2011 05:52  |
Haje Korth
Messages: 651
Registered: May 1997
|
Senior Member |
|
|
Try repeating the 0 longitude as 360. Haje
On Mar 16, 4:14 pm, Alessandro <alecarlo.dami...@gmail.com> wrote:
> Hi
>
> I’m trying to plot some satellite data with a polar stereographic
> projection. To do it I’m using the command GRIDDATA and then MAP_SET,
> 90,/Stereographic. It works fine but there is not good interpolation
> between the longitude 359 and 0. How can I fix this issue?
> Any suggestions?
>
> Thank you in advance,
> Alex
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
