code: def updateWCS(drizwcs,inwcs):
      
      crpix = np.array([drizwcs[0],drizwcs[2]], dtype=np.float64)
      crval = np.array([drizwcs[1],drizwcs[3]], dtype=np.float64)
      cd = np.array([[drizwcs[4],drizwcs[6]],[drizwcs[5],drizwcs[7]]], dtype=np.float64)
      inwcs.cd = cd
      inwcs.crval = crval
      inwc.crpix = crpix
      inwcs.pscale = N.sqrt(N.power(inwcs.cd[0][0],2)+N.power(inwcs.cd[1][0],2)) * 3600.
      inwcs.orient = N.arctan2(inwcs.cd[0][1],inwcs.cd[1][1]) * 180./N.pi
docs: Copy output WCS array from Drizzle into WCSObject.