.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "gallery/lightning-creek-magnetic.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_gallery_lightning-creek-magnetic.py: Magnetic anomaly grid of the Lightning Creek Sill Complex, Australia -------------------------------------------------------------------- This is a section of an airborne survey acquired in 1990 by the Queensland Government, Australia. The grid has 50 m resolution (UTM coordinates) and is at a uniform orthometric height of 500 m. **Original source:** `Geophysical Acquisition & Processing Section 2019. MIM Data from Mt Isa Inlier, QLD (P1029), magnetic line data, AWAGS levelled. Geoscience Australia, Canberra `__ **Pre-processing:** `Source code for preparation of the original dataset for redistribution in Ensaio `__ .. GENERATED FROM PYTHON SOURCE LINES 25-31 .. code-block:: Python import pygmt import xarray as xr import ensaio .. GENERATED FROM PYTHON SOURCE LINES 32-33 Download and cache the data and return the path to it on disk .. GENERATED FROM PYTHON SOURCE LINES 33-36 .. code-block:: Python fname = ensaio.fetch_lightning_creek_magnetic(version=1) print(fname) .. rst-class:: sphx-glr-script-out .. code-block:: none /home/runner/work/_temp/cache/ensaio/v1/lightning-creek-magnetic-grid.nc .. GENERATED FROM PYTHON SOURCE LINES 37-38 Load the netCDF grid with xarray .. GENERATED FROM PYTHON SOURCE LINES 38-41 .. code-block:: Python data = xr.load_dataarray(fname) data .. raw:: html 
<xarray.DataArray 'total_field_anomaly' (northing: 370, easting: 346)> Size: 1MB
    array([[  34.99995117,   36.19995117,   36.69995117, ..., -101.10004883,
            -100.40004883,  -99.60004883],
           [  36.49995117,   37.59995117,   37.99995117, ..., -102.20004883,
            -101.50004883, -100.70004883],
           [  37.09995117,   38.19995117,   38.59995117, ..., -103.30004883,
            -102.60004883, -101.90004883],
           ...,
           [ 182.79995117,  172.39995117,  160.79995117, ...,    0.79995117,
             -24.20004883,  -41.80004883],
           [ 182.09995117,  172.59995117,  161.39995117, ...,    5.99995117,
             -21.50004883,  -41.00004883],
           [ 178.79995117,  170.39995117,  160.29995117, ...,   11.39995117,
             -16.00004883,  -35.80004883]], shape=(370, 346))
    Coordinates:
      * northing  (northing) float64 3kB 7.576e+06 7.576e+06 ... 7.595e+06 7.595e+06
      * easting   (easting) float64 3kB 4.655e+05 4.656e+05 ... 4.827e+05 4.828e+05
        height    (northing, easting) float64 1MB 500.0 500.0 500.0 ... 500.0 500.0
    Attributes:
        Conventions:   CF-1.8
        title:         Magnetic total-field anomaly of the Lightning Creek sill c...
        crs:           proj=utm zone=54 south datum=WGS84 units=m no_defs ellps=W...
        source:        Interpolated from airborne magnetic line data using gradie...
        license:       Creative Commons Attribution 4.0 International Licence
        references:    Geophysical Acquisition & Processing Section 2019. MIM Dat...
        long_name:     total-field magnetic anomaly
        units:         nT
        actual_range:  [-1785.  3798.]


.. GENERATED FROM PYTHON SOURCE LINES 42-43 Make a PyGMT pseudo-color map of the total field magnetic anomaly. .. GENERATED FROM PYTHON SOURCE LINES 43-59 .. code-block:: Python fig = pygmt.Figure() scale = 2500 pygmt.makecpt(cmap="polar+h", series=[-scale, scale], background=True) fig.grdimage( grid=data, cmap=True, shading="+a45+nt0.1", projection="X15c/17c", frame="af", ) fig.colorbar( frame='af+l"total field magnetic anomaly [nT]"', position="JBC+h+o0/1c+e", ) fig.show() .. image-sg:: /gallery/images/sphx_glr_lightning-creek-magnetic_001.png :alt: lightning creek magnetic :srcset: /gallery/images/sphx_glr_lightning-creek-magnetic_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 60-61 The anomaly at the top right is the Lightning Creek sill complex. .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 3.228 seconds) .. _sphx_glr_download_gallery_lightning-creek-magnetic.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: lightning-creek-magnetic.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: lightning-creek-magnetic.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_