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
Exemplo n.º 2
0
 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]