Re: avhrr data extraction [message #80418 is a reply to message #80412] |
Mon, 04 June 2012 07:33  |
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
|
|
|