RE: [Tcl-nap-users] Re: Welcome to the 'Tcl-nap-users' mailing list

[<Har...@cs...>]

  > -----Original Message-----
  > From: tcl...@li...
[mailto:tcl-nap-users-
  > ad...@li...] On Behalf Of ni...@ma...
  > Sent: Sunday, 12 February 2006 21:13
  > To: tcl...@li...
  > Subject: [Tcl-nap-users] Re: Welcome to the 'Tcl-nap-users' mailing
list
  >=20
  > Greetings,
  >=20
  > ####Below we open ARC info ascii grids using nap
  > ####the grids were created using gdal_translate
  > ####iirc, WE NEED TO ADD "NODATA_value  -9999" line to the grid
files as
  > gdal doe
  > ####not seem to do this automatically...
  >=20
  > nap "red =3D get_gridascii('red.ASC')";
  > nap "near =3D get_gridascii('near.ASC')";
  > nap [result =3D (near-red)/(near+red)]
  > nap "result =3D (near-red)/(near+red)"
  > $result
  > plot_nao result -print 0 -filename result.jpeg
  >=20
  > ####The resulting plot is correct....
  >=20
  > ####Trying the same, but with using hdf files,
  >=20
  >=20
  > nap "red =3D [nap_get hdf red.hdf "3-dimensional Scientific =
Dataset"]"
  > nap "near =3D [nap_get hdf near.hdf "3-dimensional Scientific
Dataset"]"
[>>>] I gather the SDS name is "3-dimensional Scientific Dataset". I
think the spaces should be OK, but I could be wrong. It is unusual to
have such an SDS name.
  =20
  > nap "result =3D (near-red)/(near+red)"
  >=20
  > plot_nao result -print 0 -filename result1.jpeg
  >=20
  > ####The resulting plot is only the first row of the hdf SDS...
[>>>] This is strange.

  > ####Do we need to loop in tcl and read all the rows?
[>>>] No. If you do not specify an index then the whole SDS is read.

  > #### or is there an "append nao" type of method?
  >=20
  > regards,
  >=20
  > nikos
  >=20
  > thanx,
[>>>]=20
Can you send me your HDF files and your script?

Harvey