utils.py lsst::skymap::detail::utils # #LSSTDataManagementSystem #Copyright2008,2009,2010LSSTCorporation. # #Thisproductincludessoftwaredevelopedbythe #LSSTProject(http://www.lsst.org/). # #Thisprogramisfreesoftware:youcanredistributeitand/ormodify #itunderthetermsoftheGNUGeneralPublicLicenseaspublishedby #theFreeSoftwareFoundation,eitherversion3oftheLicense,or #(atyouroption)anylaterversion. # #Thisprogramisdistributedinthehopethatitwillbeuseful, #butWITHOUTANYWARRANTY;withouteventheimpliedwarrantyof #MERCHANTABILITYorFITNESSFORAPARTICULARPURPOSE.Seethe #GNUGeneralPublicLicenseformoredetails. # #YoushouldhavereceivedacopyoftheLSSTLicenseStatementand #theGNUGeneralPublicLicensealongwiththisprogram.Ifnot, #see<http://www.lsstcorp.org/LegalNotices/>. # __all__=["coordFromVec"] importnumpy importlsst.sphgeom importlsst.geomasgeom _TinyFloat=numpy.finfo(float).tiny defcoordFromVec(vec,defRA=None): """ConvertanICRScartesianvectortoanICRSlsst.geom.SpherePoint Parameters ---------- vec:`list`of`float` AnICRScatesianvector. defRA:`lsst.geom.Angle`orNone TheRAtouseifthevectoristoonearapole; ignoredifnotnearapole. Raises ------ RuntimeError IfvectoonearapoleanddefRAisNone. """ ifabs(vec[0])<_TinyFloatandabs(vec[1])<_TinyFloat: ifdefRAisNone: raiseRuntimeError("AtpoleanddefRA==None") ifvec[2]>0: decDeg=90.0 else: decDeg=-90.0 returngeom.SpherePoint(defRA,decDeg*geom.degrees) returngeom.SpherePoint(lsst.sphgeom.Vector3d(*vec))