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

Home » Public Forums » archive » Question regarding HDF file
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: Question regarding HDF file [message #55219 is a reply to message #55144] Tue, 07 August 2007 08:21 Go to previous messageGo to previous message
None[1] is currently offline  None[1]
Messages: 13
Registered: August 2007
Junior Member
>> HDF_SD_GETDATA, sdsIDsm, rr, start=[0,0], count=[x[j],y[j]],
>> stride=[0,0]
>
> HDF_SD_GetInfo, sdsIDsm, NAME=name, HDF_TYPE=htype, TYPE=itype
> PRINT, im[j], name, htype, itype, min(rr), max(rr)
>
>> HDF_SD_GETDATA, sdsIDsm2, r2
>
> HDF_SD_GetInfo, sdsIDsm2, NAME=name, HDF_TYPE=htype, TYPE=itype
> PRINT, bd[j], name, htype, itype, min(r2), max(r2)
>
> PRINT, min(rowp), max(rowp)
>

Hi. I did the corrections which you had mentioned, but still I am
getting the same results. I also changed the !PI and !DTOR.
I have changed it to read from the variable directly instead writing
to the text file and again reading from it. these are the information
which I got.

sdsIDsm= 9 ImageData DFNT_UINT8 BYTE 40 249
sdsIDsm2= 10 RadiometricCorrTable DFNT_FLOAT32 FLOAT
-3.57843 2.47200

Max(rowp) & Min(rowp) = 0.506396 0.000000

But if i open the hdf file, I am getting values like -3.5123 -1.5232
-3.5138 .............

code:

sdFileID2 = HDF_SD_Start(fname, /RdWr)
sdsIDsm = HDF_SD_Select(sdFileID2,im[j]) ; Image data
sdsIDsm2 = HDF_SD_Select(sdFileID2,bd[j]) ; recorrection data
HDF_SD_GETDATA, sdsIDsm, rr
HDF_SD_GETDATA, sdsIDsm2, r2
sz=size(rr, /DIMENSIONS)
x[j]=fix(strcompress(sz[0],/remove_all))
y[j]=fix(strcompress(sz[1],/remove_all))
results=fltarr(x[j],y[j])
rowp=fltarr(x[j],y[j])

; Reading D, A, G values
i=0
while ((i lt (y[j]-1)) && (n lt (x[j]-1))) do begin
d=r2[0,i]
a=r2[1,i]
g=r2[2,i]
results[n,*]=(float(a)*(rr[(n),*])/float(g))+float(d)
rowp[n,*]=(!PI*results[n,*]*sd*sd)/(esun[j]*COS((90-theta)*!
DTOR))
n=n+1
endwhile
HDF_SD_ADDDATA, sdsIDsm, rowp, start=[0,0], count=[x[j],y[j]],
stride=[0,0]
HDF_SD_EndAccess, sdsIDsm
HDF_SD_END, sdFileID2

Thanks
Rajesh
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Thinning algorithm without for loops
Next Topic: Re: Another simple one

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

Current Time: Sun Nov 30 07:36:09 PST 2025

Total time taken to generate the page: 0.00904 seconds