| doubt in EOS routines! [message #71735] |
Fri, 16 July 2010 17:48 |
gtbargavi
Messages: 2
Registered: March 2007
|
Junior Member |
|
|
Hi all,
I am using AMSR-E Level 3 data which is in HDF-EOS V2 format. I am
trying to subset a region with specific lat/lon values using
EOS_GD_DEFBOXREGION and EOS_GD_EXTRACTREGION routines. I would like to
know if there is any routine to check the lat/lon values obtained for
the subset. i.e. say, I am subsetting a region for lat: 30N to 40 N,
lon: 50 E to 60 E, and I would like to know if the result (subset)
obtained has correct values for this region. how do i check it?
here is my piece of code..
file = 'C:\AMSR-E_2003-04\AMSR-E_2003-04_Data
\AMSR_E_L3_DailyOcean_V04_20031106.hdf'
grid = 'GlobalGrid'
dataf = 'Med_res_vapor'
fid = eos_gd_open(file)
gid = eos_gd_attach(fid, grid)
r1 = eos_gd_pixreginfo(gid, pixregcode)
print, "px", pixregcode
s1 = eos_gd_readfield(gid, dataf, wv)
print, "s1", s1
help, wv
clon = fltarr(2)
clat = fltarr(2)
clon[0] = 94.d
clon[1] = 96.d
clat[0] = 2.d
clat[1] = 4.d
regionid = eos_gd_defboxregion(gid, clon, clat)
help, regionid
s2 = eos_gd_extractregion(gid, regionid, dataf, rdata)
print, "s2", s2
help, rdata
end
also, is 'pixel registration code' returned by 'eos_gd_pixreginfo' the
one used as 'X or Y offset' value used in the below mentioned formula
to find out 'i-th longitude or latitude value'?
(i + offsetX) * scaleX + leftX
offsetX - pixel registration
scaleX and leftX - obtained from the range of longitude
Thanks,
Maha!
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
