Is your feature request related to a problem? Please describe.

I am trying to drop particular indices from a dimension that doesn't have coordinates.

Following: [drop_sel() documentation](http://xarray.pydata.org/en/stable/generated/xarray.Dataset.drop_sel.html#xarray.Dataset.drop_se

Describe the bug
According to the multiscene documentation, the property all_same_area does:

Determine if all contained Scenes have the same ‘area’.

However, I have created a multiscene where all scenes have the same area (they just differ between datasets), yet the property returns Fa

https://earthenginepartners.appspot.com/science-2013-global-forest

There are various examples of how to define an AreaDefinition throughout the documentation and in test files like https://github.com/pytroll/pyresample/blob/master/pyresample/test/test_files/areas.yaml. However, some of these have projection definitions with only a defined (equatorial radius). This is technically incomplete and although it works in PROJ, rasterio/GDAL are a little more strict an

I've been using nb_black lately and it's wonderful. https://github.com/dnanhkhoa/nb_black

We can add this with %load_ext lab_black for our notebook examples and it will auto-format every cell to fit black to make it more clean and readable without us having to worry about it.

Looping over the data_vars in fit/predict/transform/... should be done with parallel.

It's really impressive how much pyomeca has progressed! Have you ever considered including some form of simple 3D marker visualization?