[c@`sddlmZmZmZddlmZddlZddlZddlZddl Z ddl m Z ddl mZddlmZddlmZdd lmZdd lmZdd lmZdd lmZd dlmZmZd dlm Z m!Z!d dl"m#Z#ddl$m%Z%m&Z&yddl m'Z(Wne)k rbe* r_nnXdgZ+ej,ddddgZ-dddfZ.e j/dde j0Z1e2ddZ3dZ4defd YZ5de j/fd!YZ6d d"l7m8Z8m9Z9dS(#i(tabsolute_importtdivisiontprint_function(twarnNi(tunits(t constants(tQuantityInfoBase(tsix(turllib(tAstropyUserWarning(t NUMPY_LT_1_12(t broadcast_toi(t LongitudetLatitude(tCartesianRepresentationtCartesianDifferential(tUnknownSiteException(tdatat deprecated(t_erfat EarthLocationtGeodeticLocationtlontlattheighttWGS84tGRS80tWGS72g?g?cC`s@|dkr|}n|tkr<tdj|tn|S(Ns(Ellipsoid {0} not among known ones ({1})(tNonet ELLIPSOIDSt ValueErrortformat(t ellipsoidtdefault((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt_check_ellipsoid.s     cC`sLddlm}y@tjj|dtjj}tj |j j d}Wntj j k r}t|jtjr||jddq||jd|jn,tjk r||jddnX|jdg}|s||jddn|jd dd krH||jdd n|S( Ni(tNameResolveErrorttimeouttutf8tmsgsconnection timed outtresultssno results returnedtstatustOKs$unknown failure with Google maps API(t name_resolveR#RtrequestturlopenRtconftremote_timeouttjsontloadstreadtdecodeterrortURLErrort isinstancetreasontsocketR$RtgetR(turlterr_strR#trespt resp_datateR'((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt_get_json_result7s "tEarthLocationInfocB`s,eZdZdZdZdd dZRS( s Container for meta information like name, description, format. This is required when the object is used as a mixin column within a table, but can be used as a general way to store meta information. txtytzR cC`s+|jd}|j|}||_|S(NR (tpopt _parent_clsR (tselftmapR tout((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt_construct_from_dict`s Rc `s|j||d }|jd|f|jd}tjtjd|ddjddjdtfd |j D}|j |}x-|j D]\} } t |j | | qW|S( s Return a new EarthLocation instance which is consistent with the input ``cols`` and has ``length`` rows. This is intended for creating an empty column object whose elements can be set in-place for table operations like join or vstack. Parameters ---------- cols : list List of input columns length : int Length of the output column object metadata_conflicts : str ('warn'|'error'|'silent') How to handle metadata conflicts name : str Output column name Returns ------- col : EarthLocation (or subclass) Empty instance of this class consistent with ``cols`` tmetaRt descriptiontdtypetshapeitunittcopyc`s<i|]2}|dkr"|ntd||qS(txyzi(tgetattr(t.0tkey(tcolsR(s8lib/python2.7/site-packages/astropy/coordinates/earth.pys s (RIRRJ(tmerge_cols_attributesRCtutQuantitytnptzerosRKRMtFalset_represent_as_dict_attrsRHtitemstsetattrtinfo( RERStlengthtmetadata_conflictstnametattrsRLRFRGtattrtvalue((RSRs8lib/python2.7/site-packages/astropy/coordinates/earth.pytnew_likehs  % (R@RARBR N(t__name__t __module__t__doc__RZRHRRd(((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR?Xs cB`sHeZdZdZejidddgd6ejgdd6Zejejd/fZe Z dZ e d0d Ze d d0d Ze d Ze ed Ze dZe eedZedZejdZedZd0dZeeddddddZedZeeddddddZedZedZedZ d Z!d0d!Z"ee"d"d#Z#d$Z$d%Z%d&Z&ed'Z'ed(Z(ed)Z)d*Z*d+Z+d,Z,gd-Z-e.rFd.Z/nRS(1s Location on the Earth. Initialization is first attempted assuming geocentric (x, y, z) coordinates are given; if that fails, another attempt is made assuming geodetic coordinates (longitude, latitude, height above a reference ellipsoid). When using the geodetic forms, Longitudes are measured increasing to the east, so west longitudes are negative. Internally, the coordinates are stored as geocentric. To ensure a specific type of coordinates is used, use the corresponding class methods (`from_geocentric` and `from_geodetic`) or initialize the arguments with names (``x``, ``y``, ``z`` for geocentric; ``lon``, ``lat``, ``height`` for geodetic). See the class methods for details. Notes ----- This class fits into the coordinates transformation framework in that it encodes a position on the `~astropy.coordinates.ITRS` frame. To get a proper `~astropy.coordinates.ITRS` object from this object, use the ``itrs`` property. RR@RARBtnamesitformatscO`st|dkrEt|dkrEt|dtrE|djSy|j||}Wn`tjtfk r}y|j||}Wqt k r}tdj ||qXnX|S(NiislCoordinates could not be parsed as either geocentric or geodetic, with respective exceptions "{0}" and "{1}"( tlenR5RRNtfrom_geocentricRUt UnitsErrort TypeErrort from_geodetict ExceptionR(tclstargstkwargsREtexc_geocentrict exc_geodetic((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt__new__s$ cC`sV|d kr<y |j}WqKtk r8tdqKXntj|}|jdkrltjdnyLtj||dt }tj||dt }tj||dt }Wn#tjk rtjdnXt j |||\}}}t j |j |j}||||d<|d<|d`_ . If a site cannot be found in the registry (i.e., an internet connection is not available), it will fall back on a built-in list, In the future, this bundled list might include a version-controlled list of canonical observatories extracted from the online version, but it currently only contains the Greenwich Royal Observatory as an example case. Parameters ---------- site_name : str Name of the observatory (case-insensitive). Returns ------- site : This class (a `~astropy.coordinates.EarthLocation` or subclass) The location of the observatory. See Also -------- get_site_names : the list of sites that this function can access sEarthLocation.get_site_namest close_namesN( t_get_site_registryRtsiteRt __class__Rnt to_geodeticR]R`(Rpt site_nametregistrytelR=tnewel((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytof_site/s$ c C`s tjji|d6}dj|}djd|}t|d|}|ddd}|ridjd |d d |d d 6}tjj|}d j|}d jd|}t|d|} | ddtj} nd} |jd|d tjd |d tjd| S(s Return an object of this class for a given address by querying the Google maps geocoding API. This is intended as a quick convenience function to get fast access to locations. In the background, this just issues a query to the Google maps geocoding API. It is not meant to be abused! Google uses IP-based query limiting and will ban your IP if you send more than a few thousand queries per hour [1]_. .. warning:: If the query returns more than one location (e.g., searching on ``address='springfield'``), this function will use the **first** returned location. Parameters ---------- address : str The address to get the location for. As per the Google maps API, this can be a fully specified street address (e.g., 123 Main St., New York, NY) or a city name (e.g., Danbury, CT), or etc. get_height : bool (optional) Use the retrieved location to perform a second query to the Google maps elevation API to retrieve the height of the input address [2]_. Returns ------- location : This class (a `~astropy.coordinates.EarthLocation` or subclass) The location of the input address. References ---------- .. [1] https://developers.google.com/maps/documentation/geocoding/intro .. [2] https://developers.google.com/maps/documentation/elevation/intro taddresss5https://maps.googleapis.com/maps/api/geocode/json?{0}s?Unable to retrieve coordinates for address '{address}'; {{msg}}R:itgeometrytlocations{lat:.8f},{lng:.8f}Rtlngt locationss7https://maps.googleapis.com/maps/api/elevation/json?{0}s=Unable to retrieve elevation for address '{address}'; {{msg}}t elevationgRR( Rtparset urlencodeRR>RURRnR( RpRt get_heighttparstgeo_urlR:t geo_resulttloctele_urlt ele_resultR((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt of_address`s&'    cC`s |jjS(sp Get list of names of observatories for use with `~astropy.coordinates.EarthLocation.of_site`. .. note:: When this function is called, it will first attempt to download site information from the astropy data server. If it cannot (i.e., an internet connection is not available), it will fall back on the list included with astropy (which is a limited and dated set of sites). If you think a site should be added, issue a pull request to the `astropy-data repository `_ . Returns ------- names : list of str List of valid observatory names See Also -------- of_site : Gets the actual location object for one of the sites names this returns. (RRh(Rp((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytget_site_namesscC`s|r|rtdn|r4t}|_nt|dd}|sS| ry.t|tjrwt|}n t}WnZtj j j j t fk r|rnd}tt|j|jt}nX||_n|S(s Gets the site registry. The first time this either downloads or loads from the data file packaged with astropy. Subsequent calls will use the cached version unless explicitly overridden. Parameters ---------- force_download : bool or str If not False, force replacement of the cached registry with a downloaded version. If a str, that will be used as the URL to download from (if just True, the default URL will be used). force_builtin : bool If True, load from the data file bundled with astropy and set the cache to that. returns ------- reg : astropy.coordinates.sites.SiteRegistry s6Cannot have both force_builtin and force_download Truet_site_registrysCould not access the online site list. Falling back on the built-in version, which is rather limited. If you want to retry the download, do {0}._get_site_registry(force_download=True)N(Rtget_builtin_sitesRRPRR5Rt string_typestget_downloaded_sitestmovesRR3R4tIOErrorRR RRe(Rptforce_downloadt force_builtintregR&((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRs$     cC`s|jS(s>The default ellipsoid used to convert to geodetic coordinates.(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR scC`st||_dS(N(R"R(RER ((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR scC`s |jS(s:Convert to geodetic coordinates for the default ellipsoid.(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytgeodeticscC`st|d|j}|jtjj|jtj}t j t t ||\}}}t t |tjtjddtjdtt|tjtjdttj|tj|jdtS(s Convert to geodetic coordinates. Parameters ---------- ellipsoid : str, optional Reference ellipsoid to use. Default is the one the coordinates were initialized with. Available are: 'WGS84', 'GRS80', 'WGS72' Returns ------- (lon, lat, height) : tuple The tuple contains instances of `~astropy.coordinates.Longitude`, `~astropy.coordinates.Latitude`, and `~astropy.units.Quantity` Raises ------ ValueError if ``ellipsoid`` is not recognized as among the ones implemented. Notes ----- For the conversion to geodetic coordinates, the ERFA routine ``gc2gd`` is used. See https://github.com/liberfa/erfa R!R~gf@RN(R"R RRURRt _array_dtypeRWtndarrayRtgc2gdRPRR RRRYR RVRM(RER t self_arrayRRR((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRs$$s2.0t alternatives`lon`tobj_typetpropertycC`s |jdS(s5Longitude of the location, for the default ellipsoid.i(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt longitude!scC`s |jdS(s5Longitude of the location, for the default ellipsoid.i(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR'ss`lat`cC`s |jdS(s4Latitude of the location, for the default ellipsoid.i(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytlatitude,scC`s |jdS(s5Longitude of the location, for the default ellipsoid.i(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR2scC`s |jdS(s2Height of the location, for the default ellipsoid.i(R(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR7scC`s |jS(s1Convert to a tuple with X, Y, and Z as quantities(t to_geocentric(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt geocentric=scC`s|j|j|jfS(s1Convert to a tuple with X, Y, and Z as quantities(R@RARB(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRBsc C`sw|r?|jdkr?|jdkr?t||jdt}nddlm}|d|jd|jd|jd|S(s Generates an `~astropy.coordinates.ITRS` object with the location of this object at the requested ``obstime``. Parameters ---------- obstime : `~astropy.time.Time` or None The ``obstime`` to apply to the new `~astropy.coordinates.ITRS`, or if None, the default ``obstime`` will be used. Returns ------- itrs : `~astropy.coordinates.ITRS` The new object in the ITRS frame itsubok(tITRSR@RARBtobstime( tsizeR RLtTruetbuiltin_framesRR@RARB(RERR((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytget_itrsFs$tdocsAn `~astropy.coordinates.ITRS` object with for the location of this object at the default ``obstime``.cC`suddlm}|j|}tdtjtjd|jdt}t ||j j d<|j |d|S( sYGCRS position with velocity at ``obstime`` as a GCRS coordinate. Parameters ---------- obstime : `~astropy.time.Time` The ``obstime`` to calculate the GCRS position/velocity at. Returns -------- gcrs : `~astropy.coordinates.GCRS` instance With velocity included. i(tGCRSgiRtsR(i( RRRR RUtkmRRLRRRt differentialst transform_to(RERRtitrsRX((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt _get_gcrscs *cC`s;|j|j}|j}|jdj}||fS(s Calculate the GCRS position and velocity of this object at the requested ``obstime``. Parameters ---------- obstime : `~astropy.time.Time` The ``obstime`` to calculate the GCRS position/velocity at. Returns -------- obsgeoloc : `~astropy.coordinates.CartesianRepresentation` The GCRS position of the object obsgeovel : `~astropy.coordinates.CartesianRepresentation` The GCRS velocity of the object R(RRtwithout_differentialsRt to_cartesian(RERt gcrs_datat obsgeopost obsgeovel((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytget_gcrs_posvel{s c C`sddlm}d }tjdtj}tjtj|tjf}g|D]}|||^qL}g|d D]}||dj ^qr} | j t |j j gt || D]\} } | tj| ^q} t| S( sReturn the gravitational redshift at this EarthLocation. Calculates the gravitational redshift, of order 3 m/s, due to the Sun, Jupiter, the Moon, and the Earth itself. Parameters ---------- obstime : `~astropy.time.Time` The ``obstime`` to calculate the redshift at. Returns -------- redshift : `~astropy.units.Quantity` Gravitational redshift in velocity units at given obstime. i(tget_body_barycentrictsuntjupitertmoontearthgXJ\Di(RRRR(t solar_systemRtconststGRUtkgtGM_suntGM_juptGM_earthtnormtappendRRtziptctsum( RERRRhtGM_moontmassesR`t positionstpost distancestGMtdistancet redshifts((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt_gravitational_redshifts"+1cC`s|dS(s.The X component of the geocentric coordinates.R@((RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyR@scC`s|dS(s.The Y component of the geocentric coordinates.RA((RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRAscC`s|dS(s.The Z component of the geocentric coordinates.RB((RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRBscC`sNtt|j|}|j|jkr:|j|jS|jtjSdS(N(R|Rt __getitem__RKRRRURV(REtitemtresult((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRscC`s8tt|j|t|dr4|j|_ndS(NR(R|Rt__array_finalize__thasattrR(REtobj((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRscC`s5|jdkrtdntt|jSdS(Ns*0-d EarthLocation arrays cannot be indexed((RLt IndexErrorR|Rt__len__(RE((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRscC`sg|j|jtj}|gkr0|j}n|jj||d|}|j|jj|j S(s"Helper method for to and to_value.t equivalencies( RRRWRt_equivalenciesRMRRKRRL(RERMRt array_viewt new_array((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt _to_values   cC`sXddlm}d|jjd}||jtjd|}dj|||jS(Ni(t _array2stringt( trepresentationRRReRRWRRt_unitstr(RERt prefixstrtarrstr((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyt__repr__s(iN(0ReRfRgRRWRKtfloat64R{RR?R]Rut classmethodRRkRnRRYRRRRR tsetterRRRRRRRRRRRRRRRR@RARBRRRRR R(((s8lib/python2.7/site-packages/astropy/coordinates/earth.pyRsZ  031A2 "$$           (RR(:t __future__RRRtwarningsRt collectionsR7R/tnumpyRWtRRURRtunits.quantityRtexternRtextern.six.movesRtutils.exceptionsR tutils.compat.numpycompatR tutils.compat.numpyR tanglesR R RRRterrorsRtutilsRRRRt ImportErrort_ASTROPY_SETUP_t__all__t namedtupleRRRVRt OMEGA_EARTHRR"R>R?RtsitesRR(((s8lib/python2.7/site-packages/astropy/coordinates/earth.pytsB         !AP