real from tuple ??

• To: <visad-list@xxxxxxxxxxxxx>
• Subject: real from tuple ??
• From: "Campbell, Garrett" <Campbell@xxxxxxxxxxxxxxxxxx>
• Date: Mon, 29 Oct 2001 12:02:24 -0700

I am attempting to do arithmetic with the Band4 data in the tuple
returned from
 
goese
load("adde://adde.ucar.edu/imagedata?group=rtimages&descr=ge-ir&latlon=4
0 -105&size=200 200")

the following is the display of the contents

VisAD MathType analysis
  FunctionType:
  Domain has 2 components:
   0. RealType: Latitude
   0.   Name = Latitude
   0.   Unit: deg
   1. RealType: Longitude
   1.   Name = Longitude
   1.   Unit: deg
  Range:
    RealTupleType has 1 components:
     0. RealType: Band4
     0.   Name = Band4
- - - - - - - - - - - - - - - - - - - - - - -
DataType analysis...
    FlatField of length = 357
    ((Latitude, Longitude) -> Band4)
      Domain has 2 components:
        Linear2DSet: Length = 357
          0. Linear1DSet (Latitude) Range = 39.7 to 40.3 step
0.0374999999999996
45
          1. Linear1DSet (Longitude) Range = -105.3 to -104.7 step
0.02999999999
9999714
      Range has 1 components:
        0. DoubleSet (Band4)  Dimension = 1
        0. number missing = 0

I would like to extract the values in the flatfield and compare them to
a separate array
extracted from a differet loaded object

I do not want the use the standard arithmetic functions of tuples
because it assumes
that I would want to process points at the same location.  I need to
perform the
calculation in the array space, not real space.

Following is the code, but it does not work because I can not get the
values of the data tuple.

from visad.python.JPythonMethods import *
from visad import *
from visad import Real, Unit, RealType
import graph
latc = 40.
lonc = -105
goese
load("adde://adde.ucar.edu/imagedata?group=rtimages&descr=ge-ir&latlon=4
0 -105&size=200 200")
goesw
load("adde://adde.ucar.edu/imagedata?group=rtimages&descr=gw-ir&latlon=4
0 -105&size=200 200")
print "done loading images"
graph.image(goese, colortable=None)
sampleDomain = makeDomain(RealTupleType.LatitudeLongitudeTuple, 39, 41,
61, -106, -104, 61)
d = resample(goese, sampleDomain)
print "resampled first image"
an=[]
e = goesw.resample(sampleDomain, Data.NEAREST_NEIGHBOR, Data.NO_ERRORS)
an.append(e)
f = goese.resample(sampleDomain, Data.NEAREST_NEIGHBOR, Data.NO_ERRORS)
an.append(f)
#f = goesw.resample(sampleDomain)
g = e - f
clearplot()
plot(g)
clearplot()
a=[]
dlat = .3
dlon = dlat
num=3
sampCent = makeDomain(RealTupleType.LatitudeLongitudeTuple, latc-dlat,
latc+dlat, 17, lonc-dlon,lonc+dlon, 21)
eas = goese.resample(sampCent, Data.NEAREST_NEIGHBOR, Data.NO_ERRORS)
print dumpTypes(eas)
print dumpTypes(k)
## does not work:
fe = eas.getValue()   ???
print fe
titl="east"
graph.image(eas,title=titl)
for l in range(-num,(num+1)):
   lat = latc + l*dlat/num
   latp = lat + dlat
   latm = lat - dlat
   for m in range(-num,(num+1)):
      lon = lonc + m*dlon/num
      lonp = lon + dlon
      lonm = lon - dlon
      samp = makeDomain(RealTupleType.LatitudeLongitudeTuple, latm,
latp, 17, lonm, lonp, 21)
      wes = goesw.resample(samp,Data.NEAREST_NEIGHBOR, Data.NO_ERRORS)
##   the following does not work
      for i in range(0,len(wes)):
         print wes[i].getValue()

      c = 0.
      fw = getValue(wes)
      for i in range(0,len(fw)):
         c = c + fe[i] - fw[i]
         print fw[i],fe[i],c

      print c,lat,lon
      a.append(wes)

graph,animation(a)



G. Garrett Campbell Ph. D.
Cooperative Institute for Research in the Atmosphere
Colorado State University
Ft. Collins, CO, USA

G. Garrett Campbell Ph. D.
Cooperative Institute for Research in the Atmosphere
Colorado State University
Ft. Collins, CO, USA
970 491 8497
campbell@xxxxxxxxxxxxxxxxxx 

• Follow-Ups:
  • Re: real from tuple ??
    • From: Tom Whittaker