.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        Click :ref:`here <sphx_glr_download_auto_examples_quickflat_plot_zoom_to_roi.py>`     to download the full example code
    .. rst-class:: sphx-glr-example-title

    .. _sphx_glr_auto_examples_quickflat_plot_zoom_to_roi.py:


=================================
Plotting a zoomed view of one ROI
=================================

Sometimes it is useful to create a flatmap that is zoomed in one just a single
ROI. The location of the ROI can be inferred automatically, making it easy to
show just the region around that ROI.



.. rst-class:: sphx-glr-horizontal


    *

      .. image:: /auto_examples/quickflat/images/sphx_glr_plot_zoom_to_roi_001.png
          :alt: plot zoom to roi
          :class: sphx-glr-multi-img

    *

      .. image:: /auto_examples/quickflat/images/sphx_glr_plot_zoom_to_roi_002.png
          :alt: plot zoom to roi
          :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 Out:

 .. code-block:: none

    Generating a flatmap cache






|


.. code-block:: default

    # sphinx_gallery_thumbnail_number = 2

    import cortex
    import numpy as np
    np.random.seed(1234)
    from matplotlib import pyplot as plt

    def zoom_to_roi(subject, roi, hem, margin=10.0):
        roi_verts = cortex.get_roi_verts(subject, roi)[roi]
        roi_map = cortex.Vertex.empty(subject)
        roi_map.data[roi_verts] = 1
    
        (lflatpts, lpolys), (rflatpts, rpolys) = cortex.db.get_surf(subject, "flat",
                                                                    nudge=True)
        sel_pts = dict(left=lflatpts, right=rflatpts)[hem]
        roi_pts = sel_pts[np.nonzero(getattr(roi_map, hem))[0],:2]
    
        xmin, ymin = roi_pts.min(0) - margin
        xmax, ymax = roi_pts.max(0) + margin
        plt.axis([xmin, xmax, ymin, ymax])

    # Create dataset
    data = cortex.Volume.random('S1', 'fullhead')

    # Plot it using quickflat
    cortex.quickshow(data)

    # Zoom on just one region
    zoom_to_roi('S1', 'AC', 'left')

    # notice that the quality of this figure is now quite poor/grainy
    # we can improve this by changing the 'height' argument to quickflat

    cortex.quickshow(data, height=2048)
    zoom_to_roi('S1', 'AC', 'left')

.. rst-class:: sphx-glr-timing

   **Total running time of the script:** ( 0 minutes  43.348 seconds)


.. _sphx_glr_download_auto_examples_quickflat_plot_zoom_to_roi.py:


.. only :: html

 .. container:: sphx-glr-footer
    :class: sphx-glr-footer-example



  .. container:: sphx-glr-download sphx-glr-download-python

     :download:`Download Python source code: plot_zoom_to_roi.py <plot_zoom_to_roi.py>`



  .. container:: sphx-glr-download sphx-glr-download-jupyter

     :download:`Download Jupyter notebook: plot_zoom_to_roi.ipynb <plot_zoom_to_roi.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_