comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: avhrr data extraction
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: avhrr data extraction [message #80418 is a reply to message #80412] Mon, 04 June 2012 07:33 Go to previous message
Matt[2] is currently offline  Matt[2]
Messages: 69
Registered: March 2007
Member
Hi Anil,

I think what you want is to just grab the SST data
bsloni=loni[4715:5057]
bslati=lati[1006:1142]
bsst = sst[4715:5057, 1006:1142]


That looks like Australia to me. That what you were expecting?

Hope that helps.

Matt



On Jun 4, 6:46 am, anil <akpinar.a...@gmail.com> wrote:
> Hi,
> I have downloaded some daily ,04km, night time sea surface temperature
> data from avhrr. an example:ftp://podaac-ftp.jpl.nasa.gov/allData/avhrr/L3/pathf inder_v5/daily/ni...
>
> I want to read these data and extract the values for a region (defined
> latitudes and longitudes). To read the data , I use:
>
> file=hdf_sd_start('2002252.s04d1pfv50-sst.hdf')
> indexa=hdf_sd_nametoindex(file,'lon')
> indexb=hdf_sd_nametoindex(file,'lat')
> indexc=hdf_sd_nametoindex(file,'sst')
> varida=hdf_sd_select(file,indexa)
> varidb=hdf_sd_select(file,indexb)
> varidc=hdf_sd_select(file,indexc)
> hdf_sd_getdata,varida,loni
> hdf_sd_getdata,varidb,lati
> hdf_sd_getdata,varidc,ssti
> hdf_sd_endaccess,varida
> hdf_sd_endaccess,varidb
> hdf_sd_endaccess,varidc
> hdf_sd_end,file
>
> ;bsloni=loni(4715:5057)
> ;bslati=lati(1006:1142)
> end
> I want to extract the latitudes (1006 to 1142) and longitudes (4715 to
> 5057) and the corresponding sea surface temperatures (sst's) .At the
> end I need to end up with something like:
> 40.75   36.74     sst1
> 40.75   36.88     sst2
> 40.82   37.12     sst3
> ............................... and so on. How can i do this? Defining
> a bsloni and bslati gives me the latitudes and longitudes but how do i
> get the corresponding sst's? With some kind of a where function? or a
> few for loops?or some other way?

--
Matthew Savoie - Senior Software Developer
National Snow and Ice Data Center
(303) 735-0785 http://nsidc.org
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Previous Topic: Method Overrides
Next Topic: Re: Texture Map problems when using IDLgrPolygon

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 06:56:55 PDT 2025

Total time taken to generate the page: 0.79811 seconds