diff --git a/pyorbital/orbital.py b/pyorbital/orbital.py index 15b6a7fd..e3399fff 100644 --- a/pyorbital/orbital.py +++ b/pyorbital/orbital.py @@ -98,6 +98,8 @@ def get_observer_look(sat_lon, sat_lat, sat_alt, utc_time, lon, lat, alt): Return: (Azimuth, Elevation) """ + + utc_time = dt2np(utc_time) (pos_x, pos_y, pos_z), (vel_x, vel_y, vel_z) = astronomy.observer_position( utc_time, sat_lon, sat_lat, sat_alt) @@ -126,22 +128,8 @@ def get_observer_look(sat_lon, sat_lat, sat_alt, utc_time, lon, lat, alt): az_ = np.arctan(-top_e / top_s) - if has_xarray and isinstance(az_, xr.DataArray): - az_data = az_.data - else: - az_data = az_ - - if has_dask and isinstance(az_data, da.Array): - az_data = da.where(top_s > 0, az_data + np.pi, az_data) - az_data = da.where(az_data < 0, az_data + 2 * np.pi, az_data) - else: - az_data[np.where(top_s > 0)] += np.pi - az_data[np.where(az_data < 0)] += 2 * np.pi - - if has_xarray and isinstance(az_, xr.DataArray): - az_.data = az_data - else: - az_ = az_data + az_ = np.where(top_s > 0, az_ + np.pi, az_) + az_ = np.where(az_ < 0, az_ + 2 * np.pi, az_) rg_ = np.sqrt(rx * rx + ry * ry + rz * rz) el_ = np.arcsin(top_z / rg_)