ó ›ßÈ[c@`sšddlmZmZmZddlZddlmZddlm Z ddl m Z ddl mZmZdgZd „Zdefd „ƒYZdS( i(tprint_functiontdivisiontabsolute_importN(tPolygoni(tunits(tUnitSphericalRepresentation(trotation_matrixtmatrix_producttSphericalCirclecC`s–td|d|ƒ}tdtjtj| ddƒ}t| ddƒ}t||ƒ}|jƒ}|j|ƒ}tj |ƒ}|j |j fS(s& Given a polygon with vertices defined by (lon, lat), rotate the polygon such that the North pole of the spherical coordinates is now at (lon0, lat0). Therefore, to end up with a polygon centered on (lon0, lat0), the polygon should initially be drawn around the North pole. tlontlatgà?taxistytz( RRtnptpitutradianRt to_cartesiant transformtfrom_cartesianR R (R R tlon0tlat0tpolygontm1tm2ttransform_matrix((sDlib/python2.7/site-packages/astropy/visualization/wcsaxes/patches.pyt_rotate_polygons % cB`s eZdZdejd„ZRS(sï Create a patch representing a spherical circle - that is, a circle that is formed of all the points that are within a certain angle of the central coordinates on a sphere. Here we assume that latitude goes from -90 to +90 This class is needed in cases where the user wants to add a circular patch to a celestial image, since otherwise the circle will be distorted, because a fixed interval in longitude corresponds to a different angle on the sky depending on the latitude. Parameters ---------- center : tuple or `~astropy.units.Quantity` This can be either a tuple of two `~astropy.units.Quantity` objects, or a single `~astropy.units.Quantity` array with two elements. radius : `~astropy.units.Quantity` The radius of the circle resolution : int, optional The number of points that make up the circle - increase this to get a smoother circle. vertex_unit : `~astropy.units.Unit` The units in which the resulting polygon should be defined - this should match the unit that the transformation (e.g. the WCS transformation) expects as input. Notes ----- Additional keyword arguments are passed to `~matplotlib.patches.Polygon` idc K`sØ|\}}tjddtj|dƒd tj}tjdtj|jtjƒ|ƒtj} t|| ||ƒ\}} |j|ƒ}| j|ƒ} tj|| gƒj ƒ} t t |ƒj | |dS(Ngiiiÿÿÿÿgà?( RtlinspaceRRRtrepeattto_valueRtarrayt transposetsuperRt__init__( tselftcentertradiust resolutiont vertex_unittkwargst longitudetlatitudeR R tvertices((sDlib/python2.7/site-packages/astropy/visualization/wcsaxes/patches.pyR"Hs +0(t__name__t __module__t__doc__RtdegreeR"(((sDlib/python2.7/site-packages/astropy/visualization/wcsaxes/patches.pyR)s(t __future__RRRtnumpyRtmatplotlib.patchesRtRRtcoordinates.representationRtcoordinates.matrix_utilitiesRRt__all__RR(((sDlib/python2.7/site-packages/astropy/visualization/wcsaxes/patches.pyts