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]