---------------------------------------------------------------------------
UnitsError Traceback (most recent call last)
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/IPython/core/formatters.pyc in __call__(self, obj)
305 pass
306 else:
--> 307 return printer(obj)
308 # Finally look for special method names
309 method = get_real_method(obj, self.print_method)
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.pyc in <lambda>(fig)
225
226 if 'png' in formats:
--> 227 png_formatter.for_type(Figure, lambda fig: print_figure(fig, 'png', **kwargs))
228 if 'retina' in formats or 'png2x' in formats:
229 png_formatter.for_type(Figure, lambda fig: retina_figure(fig, **kwargs))
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.pyc in print_figure(fig, fmt, bbox_inches, **kwargs)
117
118 bytes_io = BytesIO()
--> 119 fig.canvas.print_figure(bytes_io, **kw)
120 data = bytes_io.getvalue()
121 if fmt == 'svg':
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/backend_bases.pyc in print_figure(self, filename, dpi, facecolor, edgecolor, orientation, format, **kwargs)
2206 orientation=orientation,
2207 dryrun=True,
-> 2208 **kwargs)
2209 renderer = self.figure._cachedRenderer
2210 bbox_inches = self.figure.get_tightbbox(renderer)
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc in print_png(self, filename_or_obj, *args, **kwargs)
505
506 def print_png(self, filename_or_obj, *args, **kwargs):
--> 507 FigureCanvasAgg.draw(self)
508 renderer = self.get_renderer()
509 original_dpi = renderer.dpi
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc in draw(self)
428 if toolbar:
429 toolbar.set_cursor(cursors.WAIT)
--> 430 self.figure.draw(self.renderer)
431 finally:
432 if toolbar:
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc in draw_wrapper(artist, renderer, *args, **kwargs)
53 renderer.start_filter()
54
---> 55 return draw(artist, renderer, *args, **kwargs)
56 finally:
57 if artist.get_agg_filter() is not None:
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/figure.pyc in draw(self, renderer)
1293
1294 mimage._draw_list_compositing_images(
-> 1295 renderer, self, artists, self.suppressComposite)
1296
1297 renderer.close_group('figure')
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/image.pyc in _draw_list_compositing_images(renderer, parent, artists, suppress_composite)
136 if not_composite or not has_images:
137 for a in artists:
--> 138 a.draw(renderer)
139 else:
140 # Composite any adjacent images together
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/core.pyc in draw(self, renderer, inframe)
338 coords.frame.update()
339 for coord in coords:
--> 340 coord._draw_grid(renderer)
341
342 for coords in self._all_coords:
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/coordinate_helpers.pyc in _draw_grid(self, renderer)
432 renderer.open_group('grid lines')
433
--> 434 self._update_ticks()
435
436 if self.grid_lines_kwargs['visible']:
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/coordinate_helpers.pyc in _update_ticks(self)
487
488 # Find the range of coordinates in all directions
--> 489 coord_range = self.parent_map.get_coord_range()
490
491 # First find the ticks we want to show
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/coordinates_map.pyc in get_coord_range(self)
161 [xmin, xmax, ymin, ymax],
162 [coord.coord_type for coord in self],
--> 163 [coord.coord_unit for coord in self])
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/coordinate_range.pyc in find_coordinate_range(transform, extent, coord_types, coord_units)
45 y = np.linspace(extent[2], extent[3], ny + 1)
46 xp, yp = np.meshgrid(x, y)
---> 47 world = transform.transform(np.vstack([xp.ravel(), yp.ravel()]).transpose())
48
49 ranges = []
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/transforms.pyc in transform(self, values)
1421
1422 # Transform the values
-> 1423 res = self.transform_affine(self.transform_non_affine(values))
1424
1425 # Convert the result back to the shape of the input values.
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/matplotlib/transforms.pyc in transform_non_affine(self, points)
2458 else:
2459 return self._b.transform_non_affine(
-> 2460 self._a.transform(points))
2461 transform_non_affine.__doc__ = Transform.transform_non_affine.__doc__
2462
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/visualization/wcsaxes/transforms.pyc in transform(self, input_coords)
249 # on all values and just ignore Numpy warnings
250 with np.errstate(all='ignore'):
--> 251 c_out = c_in.transform_to(self.output_system)
252
253 if issubclass(c_out.representation, (SphericalRepresentation, UnitSphericalRepresentation)):
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/coordinates/sky_coordinate.pyc in transform_to(self, frame, merge_attributes)
479 # Do the transformation, returning a coordinate frame of the desired
480 # final type (not generic).
--> 481 new_coord = trans(self.frame, generic_frame)
482
483 # Finally make the new SkyCoord object from the `new_coord` and
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/coordinates/transformations.pyc in __call__(self, fromcoord, toframe)
1312
1313 curr_toframe = t.tosys(**frattrs)
-> 1314 curr_coord = t(curr_coord, curr_toframe)
1315
1316 # this is safe even in the case where self.transforms is empty, because
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/coordinates/transformations.pyc in __call__(self, fromcoord, toframe)
1127 def __call__(self, fromcoord, toframe):
1128
-> 1129 M, vec = self.transform_func(fromcoord, toframe)
1130 newrep = self._apply_transform(fromcoord, M, vec)
1131
/Users/kkozarev/anaconda2/lib/python2.7/site-packages/astropy/coordinates/builtin_frames/icrs_cirs_transforms.pyc in icrs_to_hcrs(icrs_coo, hcrs_frame)
310 # this is just an origin translation so without a distance it cannot go ahead
311 if isinstance(icrs_coo.data, UnitSphericalRepresentation):
--> 312 raise u.UnitsError(_NEED_ORIGIN_HINT.format(icrs_coo.__class__.__name__))
313
314 if icrs_coo.data.differentials:
UnitsError: The input ICRS coordinates do not have length units. This probably means you created coordinates with lat/lon but no distance. Heliocentric<->ICRS transforms cannot function in this case because there is an origin shift.