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

[<ni...@ma...>]

Greetings,

####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...

nap "red = get_gridascii('red.ASC')";
nap "near = get_gridascii('near.ASC')";
nap [result = (near-red)/(near+red)]
nap "result = (near-red)/(near+red)"
$result
plot_nao result -print 0 -filename result.jpeg

####The resulting plot is correct....

####Trying the same, but with using hdf files,


nap "red = [nap_get hdf red.hdf "3-dimensional Scientific Dataset"]"
nap "near = [nap_get hdf near.hdf "3-dimensional Scientific Dataset"]"


nap "result = (near-red)/(near+red)"

plot_nao result -print 0 -filename result1.jpeg

####The resulting plot is only the first row of the hdf SDS...
####Do we need to loop in tcl and read all the rows?
#### or is there an "append nao" type of method?

regards,

nikos

thanx,