Hi Joleen -
I'm having trouble replicating your exact problems on either Windows or Mac. Here's what I'm doing:
1. Through the GUI, individually add the following 3 files:
SVI01_npp_d20140122_t1843063_e1844304_b11597_c20140123010724640954_noaa_ops.h5
SVI02_npp_d20140122_t1843063_e1844304_b11597_c20140123010718218917_noaa_ops.h5
SVM15_npp_d20140122_t1843063_e1844304_b11597_c20140123010700451913_noaa_ops.h5
2. Open the Jython Shell, and run the following:
Code: Select all
I1=selectData('I1') #SVI01 Radiance
I2=selectData('I2') #SVI02 Radiance
T15=selectData('T15') #SVM15 Radiance
ds = I1.getDomainSet()
rI2=resample(I2,ds)
rT15=resample(T15,ds)
3. In the first Field Selector window, I chose SVI01 Radiance and subsetted out an area near the center of the granule. In the next 2 Field Selector windows, I selected SVI02 radiance and SVM15 radiance without subsetting a region.
4. I then created a display window with 3 panels, and ran the following commands in the Jython Shell for each panel in the display:
Code: Select all
layer1=activeDisplay().createLayer('Image Display',I1)
layer2=activeDisplay().createLayer('Image Display',rI2)
layer3=activeDisplay().createLayer('Image Display',rT15)
- Before running each line, I just clicked on my desired panel in the display so activeDisplay() would place the display in the correct panel
When I do this, I get results like the display attached to my first post, where everything appears to be valid data that can be probed... just the SVM15 doesn't correctly pick up on the domain from SVI01. So it seems like I'm not getting your results of a display with all missing values.
If I try multiplying the M band by the I band without any subsetting, I'm seeing a strange results, but I'm not seeing an image of 1 to 1. In the Jython Shell, I ran the following:
Code: Select all
I1=selectData('I1') #SVI01 Radiance
T15=selectData('T15') #SVM15 Radiance
IbyM = mul(I1,T15)
MbyI = mul(T15,I1)
panel1=buildWindow()
panel1[0].createLayer('Image Display',IbyM)
panel2=buildWindow()
panel2[0].createLayer('Image Display',MbyI)
For the display of 'IbyM', the data displays, but not correctly... it is geographically shifted to the east of where it should be and the bowtie takes up more of the granule than it should. The data can be probed, however. For the display of 'MbyI', the data displays at the correct location, but again the bowtie extends further into the granule than it should, and the data can be probed fine. This problem may require a new inquiry. I get the same results if I do this w/o scripting with a simple formula: "mul(a,b)".
Since we are seeing different results, I'm wondering if it is because we are using different data. Are you using single-banded data from peate, or multi-banded data from CLASS? Can you please post your data on ftp for me to try, myself? Also, when you subset your data for SVI01, are you choosing an area that goes right up to the edge of the granule?
Thanks again -
Bob