Fix indexing/slicing into sliced CrystalMap

[hakonanes]

Description of the change

I discovered a bug where indexing/slicing into an already sliced CrystalMap assumed the initial shape (_original_shape) and not the current shape.

Progress of the PR

• Docstrings for all functions
• Unit tests with pytest for all lines
• Clean code style by running black via pre-commit
• Update release version to 0.10.1
• Merge after Deprecate minimal/brittle "3D" functionality for crystal maps #392 and Allow to pass c/color and s/sizes to StereographicPlot.scatter() #393 are merged

Minimal example of the bug fix or new feature

>>> from orix.crystal_map import CrystalMap
>>> xmap = CrystalMap.empty((4, 5))
>>> xmap.id.reshape(xmap.shape)
array([[ 0,  1,  2,  3,  4],
       [ 5,  6,  7,  8,  9],
       [10, 11, 12, 13, 14],
       [15, 16, 17, 18, 19]])
>>> xmap2 = xmap[1:4, :3]
>>> xmap2._original_shape
(4, 5)
>>> xmap2.id.reshape(xmap2.shape)
array([[ 5,  6,  7],
       [10, 11, 12],
       [15, 16, 17]])

Getting some IDs of xmap2 before this fix

>>> xmap2[1, 1].id
array([11])
>>> xmap2[:2, :2].id
array([5, 6])

and after this fix

>>> xmap2[1, 1].id
array([13])
>>> xmap2[:2, :2].id
array([ 5,  6, 10, 11])

For reviewers

• The PR title is short, concise, and will make sense 1 year later.
• New functions are imported in corresponding __init__.py.
• New features, API changes, and deprecations are mentioned in the unreleased
  section in CHANGELOG.rst.
• Contributor(s) are listed correctly in __credits__ in orix/__init__.py and in
  .zenodo.json.

[pc494]

also good, merging.