print 'Vis nominal sample resolution:',num.min(sampResolutionsVis[W])
resRat=num.min(sampResolutionsVis[W])/sampResolutionsIR[argI]
print resRat
###inot=int(CassIR[0,5])
###jnot=int(CassIR[0,6])
###Tnot=CassIR[0,3]
w=[[0]]
#mid_i=int(num.median(CassIR[:,5]))
#w=num.where(CassIR[:,5]==mid_i)
#mid_j=int(num.median(CassIR[w][:,6]))
#w=num.where((CassIR[:,5]==mid_i)&(CassIR[:,6]==mid_j))
inot=int(CassIR[w][:,5])#same as mid_i
jnot=int(CassIR[w][:,6])#same as mid_j
Tnot=CassIR[w][0][3]
pixelTimesVis=getPixTimes(CassVis,A,B,Tnot,inot,jnot)
pixelTimesIR=getPixTimes(CassIR,A,B,Tnot,inot,jnot)
#distance units are in km, resolution units in km/pix
print CassIR[w][0][:3]
[Xnot,Ynot,Znot]=CassIR[w][0][:3]#latLongObjIR[inot][jnot]['SpacecraftPosition']/2
distancenot=CassIR[w][0][0]#latLongObjIR[inot][jnot]['TargetCenterDistance']/2
#sampleResolutionnot=latLongObjIR[inot][jnot]['SampleResolution']/1000.
#nominal values from the VIMS manual are 0.17*3 and 0.5 mrad
if resRat<0.5:
pixScaleVis=0.17
else:
pixScaleVis=0.17*3
argI=num.argmin(sampResolutionsIR)
timeDelt=num.abs(CassVis[:,3]-CassIR[argI][3])
W=num.where(timeDelt==num.min(timeDelt))
resolutions.append(sampResolutionsIR[argI])
w=[[0]]
inot=int(CassIR[w][:,5])#same as mid_i
jnot=int(CassIR[w][:,6])#same as mid_j
Tnot=CassIR[w][0][3]
pixelTimesVis=getPixTimes(CassVis,imDataVis.shape[0],imDataVis.shape[1],Tnot,inot,jnot)
pixelTimesIR=getPixTimes(CassIR,imDataIR.shape[0],imDataIR.shape[1],Tnot,inot,jnot)
#distance units are in km, resolution units in km/pix
[Xnot,Ynot,Znot]=CassIR[w][0][:3]
if distancenot==None:
distancenot=CassIR[w][0][4]
else:
Xnot*=distancenot/CassIR[w][0][0]
Ynot*=distancenot/CassIR[w][0][0]
Znot*=distancenot/CassIR[w][0][0]
#nominal values from the VIMS manual are 0.17*3 and 0.5 mrad
resRat=num.min(sampResolutionsVis[W])/sampResolutionsIR[argI]