rtc_util.py 文件源码

def mk_rad(jmax,imax,inc,dem,sang,tau,ref,eref):
    rad=np.zeros(imax*jmax).reshape(jmax,imax)
    #t=cv2.getTickCount() 
    for j in range(jmax):
      rad[j,:] = radiance(ref[j,:],inc[j,:],tau[j,:],dem[j,:]/1000,eref[j,:],sang[j,:])
      #if j % 100 == 0: print j,(cv2.getTickCount()-t)/cv2.getTickFrequency()
    print np.mean(rad),np.std(rad)
    return rad