GeoTIFFs to NetCDF

In this example, we demonstrates the workflow for creating a GS file from the GeoTIFF (.tif/.tiff) file format. This includes adding individual TIF files as single 2-D variables, as well as how to create a 3-D variable by stacking multiple TIF files along a specified dimension.

Additionally, this example shows how to handle Raster data that have differing x-y grids. Specifically, this example creates the following Raster datasets:

1. Raster Dataset #1

   1a. 2-D magnetic grid, original x-y discretization (600 m cell size)

2. Raster Dataset #2

   2a. 2-D magnetic grid, aligned to match the x-y dimensions of the resistivity layers (1000 m cell size)

   2b. 3-D resistivity grid

Dataset References:

Minsley, B.J., James, S.R., Bedrosian, P.A., Pace, M.D., Hoogenboom, B.E., and Burton, B.L., 2021, Airborne electromagnetic, magnetic, and radiometric survey of the Mississippi Alluvial Plain, November 2019 - March 2020: U.S. Geological Survey data release, https://doi.org/10.5066/P9E44CTQ.

James, S.R., and Minsley, B.J., 2021, Combined results and derivative products of hydrogeologic structure and properties from airborne electromagnetic surveys in the Mississippi Alluvial Plain: U.S. Geological Survey data release, https://doi.org/10.5066/P9382RCI.

import matplotlib.pyplot as plt
from os.path import join
import gspy
from gspy import Survey
from pprint import pprint

Convert data from GeoTIFF to NetCDF

Initialize the Survey

# Path to example files
data_path = "..//data_files//tempest_aseg"

# Survey metadata file
metadata = join(data_path, "data//Tempest_survey_md.yml")

# Establish the Survey
survey = Survey.from_dict(metadata)


container = survey.gs.add_container('derived_products', **dict(content = "raw and processed data",
                                                        comment = "This is a test"))

Create the First Raster Dataset Import 2-D magnetic data, discretized on 600 m x 600 m grid Define input metadata file (which contains the TIF filename linked with desired variable name)

d_supp1 = join(data_path, 'data//Tempest_raster_md.yml')

# Read data and format as Raster class object
container.gs.add(key="map", metadata_file=d_supp1)

<xarray.DataTree 'map'>
Group: /survey/derived_products/map
    Dimensions:       (x: 599, nv: 2, y: 1212)
    Coordinates:
      * x             (x) float64 5kB 2.928e+05 2.934e+05 ... 6.51e+05 6.516e+05
      * nv            (nv) int64 16B 0 1
      * y             (y) float64 10kB 1.607e+06 1.606e+06 ... 8.808e+05 8.802e+05
        spatial_ref   float64 8B 0.0
    Data variables:
        x_bnds        (x, nv) float64 10kB 2.925e+05 2.931e+05 ... 6.519e+05
        y_bnds        (y, nv) float64 19kB 1.607e+06 1.606e+06 ... 8.799e+05
        magnetic_tmi  (y, x) float64 6MB nan nan nan nan nan ... nan nan nan nan nan
    Attributes:
        uuid:     d56fd019-c23a-4bff-9af4-cbaaa4f6d5dd
        comment:  <additional details or ancillary information>
        content:  gridded magnetic map

xarray.DataTree

'map'

• Dimensions:
  • x: 599
  • nv: 2
  • y: 1212
• Coordinates: (4)
  • x
    (x)
    float64
    2.928e+05 2.934e+05 ... 6.516e+05

    standard_name :

    projection_x_coordinate

    long_name :

    Easting

    units :

    meter

    null_value :

    not_defined

    axis :

    X

    valid_range :

    [292800. 651600.]

    grid_mapping :

    spatial_ref

    bounds :

    x_bnds

    array([292800., 293400., 294000., ..., 650400., 651000., 651600.], shape=(599,))

  • nv
    (nv)
    int64
    0 1

    standard_name :

    nv

    long_name :

    Number of vertices for bounding variables

    units :

    not_defined

    null_value :

    not_defined

    valid_range :

    [0. 1.]

    grid_mapping :

    spatial_ref

    array([0, 1])

  • y
    (y)
    float64
    1.607e+06 1.606e+06 ... 8.802e+05

    standard_name :

    projection_y_coordinate

    long_name :

    Northing

    units :

    meter

    null_value :

    not_defined

    axis :

    Y

    valid_range :

    [ 880200. 1606800.]

    grid_mapping :

    spatial_ref

    bounds :

    y_bnds

    array([1606800., 1606200., 1605600., ...,  881400.,  880800.,  880200.],
          shape=(1212,))

  • spatial_ref
    ()
    float64
    0.0

    crs_wkt :

    PROJCRS["NAD83 / Conus Albers",BASEGEOGCRS["NAD83",DATUM["North American Datum 1983",ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4269]],CONVERSION["Conus Albers",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",23,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",-96,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",29.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",45.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["easting (X)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["northing (Y)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Data analysis and small scale data presentation for contiguous lower 48 states."],AREA["United States (USA) - CONUS onshore - Alabama; Arizona; Arkansas; California; Colorado; Connecticut; Delaware; Florida; Georgia; Idaho; Illinois; Indiana; Iowa; Kansas; Kentucky; Louisiana; Maine; Maryland; Massachusetts; Michigan; Minnesota; Mississippi; Missouri; Montana; Nebraska; Nevada; New Hampshire; New Jersey; New Mexico; New York; North Carolina; North Dakota; Ohio; Oklahoma; Oregon; Pennsylvania; Rhode Island; South Carolina; South Dakota; Tennessee; Texas; Utah; Vermont; Virginia; Washington; West Virginia; Wisconsin; Wyoming."],BBOX[24.41,-124.79,49.38,-66.91]],ID["EPSG",5070]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314140356

    inverse_flattening :

    298.257222101

    reference_ellipsoid_name :

    GRS 1980

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    NAD83

    horizontal_datum_name :

    North American Datum 1983

    projected_crs_name :

    NAD83 / Conus Albers

    grid_mapping_name :

    albers_conical_equal_area

    standard_parallel :

    (29.5, 45.5)

    latitude_of_projection_origin :

    23.0

    longitude_of_central_meridian :

    -96.0

    false_easting :

    0.0

    false_northing :

    0.0

    authority :

    EPSG

    wkid :

    5070

    GeoTransform :

    (292500.0, 600.0, 0.0, 1607100.0, 0.0, -600.0)

    array(0.)

• Data variables: (3)
  • x_bnds
    (x, nv)
    float64
    2.925e+05 2.931e+05 ... 6.519e+05

    standard_name :

    projection_x_coordinate_bounds

    long_name :

    Easting cell boundaries

    units :

    meter

    null_value :

    not_defined

    axis :

    X

    valid_range :

    [292500. 651900.]

    grid_mapping :

    spatial_ref

    array([[292500., 293100.],
           [293100., 293700.],
           [293700., 294300.],
           ...,
           [650100., 650700.],
           [650700., 651300.],
           [651300., 651900.]], shape=(599, 2))

  • y_bnds
    (y, nv)
    float64
    1.607e+06 1.606e+06 ... 8.799e+05

    standard_name :

    projection_y_coordinate_bounds

    long_name :

    Northing cell boundaries

    units :

    meter

    null_value :

    not_defined

    axis :

    Y

    valid_range :

    [ 879900. 1607100.]

    grid_mapping :

    spatial_ref

    array([[1607100., 1606500.],
           [1606500., 1605900.],
           [1605900., 1605300.],
           ...,
           [ 881700.,  881100.],
           [ 881100.,  880500.],
           [ 880500.,  879900.]], shape=(1212, 2))

  • magnetic_tmi
    (y, x)
    float64
    nan nan nan nan ... nan nan nan nan

    standard_name :

    total_magnetic_intensity

    long_name :

    Total magnetic intensity, diurnally corrected and filtered

    units :

    nT

    null_value :

    1.70141e+38

    valid_range :

    [-17504.6640625 11490.32324219]

    grid_mapping :

    spatial_ref

    array([[          nan,           nan,           nan, ..., -437.66494751,
            -438.0090332 , -435.30517578],
           [          nan,           nan,           nan, ..., -429.26696777,
            -428.53414917, -426.61795044],
           [          nan,           nan,           nan, ..., -420.28356934,
            -419.06771851, -417.57122803],
           ...,
           [          nan,           nan,           nan, ...,           nan,
                      nan,           nan],
           [          nan,           nan,           nan, ...,           nan,
                      nan,           nan],
           [          nan,           nan,           nan, ...,           nan,
                      nan,           nan]], shape=(1212, 599))

• Attributes: (3)

  uuid :

  d56fd019-c23a-4bff-9af4-cbaaa4f6d5dd

  comment :

  <additional details or ancillary information>

  content :

  gridded magnetic map

Create the Second Raster Dataset

# Import both 3-D resistivity and 2-D magnetic data, aligned onto a common 1000 m x 1000 m grid
# Define input metadata file (which contains the TIF filenames linked with desired variable names)
d_supp2 = join(data_path, 'data//Tempest_rasters_md.yml')

# Read data and format as Raster class object
container.gs.add(key="maps", metadata_file=d_supp2)

<xarray.DataTree 'maps'>
Group: /survey/derived_products/maps
    Dimensions:       (z: 5, nv: 2, x: 363, y: 770)
    Coordinates:
      * z             (z) float64 40B 0.0 5.0 10.0 15.0 20.0
      * nv            (nv) int64 16B 0 1
      * x             (x) float64 3kB 2.915e+05 2.925e+05 ... 6.525e+05 6.535e+05
      * y             (y) float64 6kB 1.648e+06 1.647e+06 ... 8.798e+05 8.788e+05
        spatial_ref   float64 8B 0.0
    Data variables:
        z_bnds        (z, nv) float64 80B -2.5 2.5 2.5 7.5 ... 12.5 17.5 17.5 22.5
        x_bnds        (x, nv) float64 6kB 2.91e+05 2.92e+05 ... 6.53e+05 6.54e+05
        y_bnds        (y, nv) float64 12kB 1.648e+06 1.647e+06 ... 8.783e+05
        resistivity   (z, y, x) float64 11MB nan nan nan nan nan ... nan nan nan nan
        magnetic_tmi  (y, x) float32 1MB nan nan nan nan nan ... nan nan nan nan nan
    Attributes:
        uuid:     fff129ce-4888-4dd9-98a6-4fb66ceb64b1
        comment:  <additional details or ancillary information>
        content:  interpolated resistivity models (3-D depth grid) and magnetic m...

xarray.DataTree

'maps'

• Dimensions:
  • z: 5
  • nv: 2
  • x: 363
  • y: 770
• Coordinates: (5)
  • z
    (z)
    float64
    0.0 5.0 10.0 15.0 20.0

    comment :

    depth is defined here inside our json file in order to stack the resistivity tifs along a third dimension. This depth information is not contained in any of the individual tifs.

    standard_name :

    z

    long_name :

    Depth below earth's surface DTM

    units :

    meters

    null_value :

    not_defined

    axis :

    Z

    positive :

    down

    datum :

    ground surface

    valid_range :

    [ 0. 20.]

    grid_mapping :

    spatial_ref

    bounds :

    z_bnds

    array([ 0.,  5., 10., 15., 20.])

  • nv
    (nv)
    int64
    0 1

    standard_name :

    nv

    long_name :

    Number of vertices for bounding variables

    units :

    not_defined

    null_value :

    not_defined

    valid_range :

    [0. 1.]

    grid_mapping :

    spatial_ref

    array([0, 1])

  • x
    (x)
    float64
    2.915e+05 2.925e+05 ... 6.535e+05

    standard_name :

    projection_x_coordinate

    long_name :

    Easting

    units :

    meter

    null_value :

    not_defined

    axis :

    X

    valid_range :

    [291455. 653455.]

    grid_mapping :

    spatial_ref

    bounds :

    x_bnds

    array([291455., 292455., 293455., ..., 651455., 652455., 653455.], shape=(363,))

  • y
    (y)
    float64
    1.648e+06 1.647e+06 ... 8.788e+05

    standard_name :

    projection_y_coordinate

    long_name :

    Northing

    units :

    meter

    null_value :

    not_defined

    axis :

    Y

    valid_range :

    [ 878785. 1647785.]

    grid_mapping :

    spatial_ref

    bounds :

    y_bnds

    array([1647785., 1646785., 1645785., ...,  880785.,  879785.,  878785.],
          shape=(770,))

  • spatial_ref
    ()
    float64
    0.0

    crs_wkt :

    PROJCRS["NAD83 / Conus Albers",BASEGEOGCRS["NAD83",DATUM["North American Datum 1983",ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4269]],CONVERSION["Conus Albers",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",23,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",-96,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",29.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",45.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["easting (X)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["northing (Y)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Data analysis and small scale data presentation for contiguous lower 48 states."],AREA["United States (USA) - CONUS onshore - Alabama; Arizona; Arkansas; California; Colorado; Connecticut; Delaware; Florida; Georgia; Idaho; Illinois; Indiana; Iowa; Kansas; Kentucky; Louisiana; Maine; Maryland; Massachusetts; Michigan; Minnesota; Mississippi; Missouri; Montana; Nebraska; Nevada; New Hampshire; New Jersey; New Mexico; New York; North Carolina; North Dakota; Ohio; Oklahoma; Oregon; Pennsylvania; Rhode Island; South Carolina; South Dakota; Tennessee; Texas; Utah; Vermont; Virginia; Washington; West Virginia; Wisconsin; Wyoming."],BBOX[24.41,-124.79,49.38,-66.91]],ID["EPSG",5070]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314140356

    inverse_flattening :

    298.257222101

    reference_ellipsoid_name :

    GRS 1980

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    NAD83

    horizontal_datum_name :

    North American Datum 1983

    projected_crs_name :

    NAD83 / Conus Albers

    grid_mapping_name :

    albers_conical_equal_area

    standard_parallel :

    (29.5, 45.5)

    latitude_of_projection_origin :

    23.0

    longitude_of_central_meridian :

    -96.0

    false_easting :

    0.0

    false_northing :

    0.0

    authority :

    EPSG

    wkid :

    5070

    GeoTransform :

    (290955.0, 1000.0, 0.0, 1648285.0, 0.0, -1000.0)

    array(0.)

• Data variables: (5)
  • z_bnds
    (z, nv)
    float64
    -2.5 2.5 2.5 7.5 ... 17.5 17.5 22.5

    comment :

    depth is defined here inside our json file in order to stack the resistivity tifs along a third dimension. This depth information is not contained in any of the individual tifs.

    standard_name :

    z_bounds

    long_name :

    Depth below earth's surface DTM cell boundaries

    units :

    meters

    null_value :

    not_defined

    axis :

    Z

    positive :

    down

    datum :

    ground surface

    valid_range :

    [-2.5 22.5]

    grid_mapping :

    spatial_ref

    array([[-2.5,  2.5],
           [ 2.5,  7.5],
           [ 7.5, 12.5],
           [12.5, 17.5],
           [17.5, 22.5]])

  • x_bnds
    (x, nv)
    float64
    2.91e+05 2.92e+05 ... 6.54e+05

    standard_name :

    projection_x_coordinate_bounds

    long_name :

    Easting cell boundaries

    units :

    meter

    null_value :

    not_defined

    axis :

    X

    valid_range :

    [290955. 653955.]

    grid_mapping :

    spatial_ref

    array([[290955., 291955.],
           [291955., 292955.],
           [292955., 293955.],
           [293955., 294955.],
           [294955., 295955.],
           [295955., 296955.],
           [296955., 297955.],
           [297955., 298955.],
           [298955., 299955.],
           [299955., 300955.],
           [300955., 301955.],
           [301955., 302955.],
           [302955., 303955.],
           [303955., 304955.],
           [304955., 305955.],
           [305955., 306955.],
           [306955., 307955.],
           [307955., 308955.],
           [308955., 309955.],
           [309955., 310955.],
    ...
           [633955., 634955.],
           [634955., 635955.],
           [635955., 636955.],
           [636955., 637955.],
           [637955., 638955.],
           [638955., 639955.],
           [639955., 640955.],
           [640955., 641955.],
           [641955., 642955.],
           [642955., 643955.],
           [643955., 644955.],
           [644955., 645955.],
           [645955., 646955.],
           [646955., 647955.],
           [647955., 648955.],
           [648955., 649955.],
           [649955., 650955.],
           [650955., 651955.],
           [651955., 652955.],
           [652955., 653955.]])

  • y_bnds
    (y, nv)
    float64
    1.648e+06 1.647e+06 ... 8.783e+05

    standard_name :

    projection_y_coordinate_bounds

    long_name :

    Northing cell boundaries

    units :

    meter

    null_value :

    not_defined

    axis :

    Y

    valid_range :

    [ 878285. 1648285.]

    grid_mapping :

    spatial_ref

    array([[1648285., 1647285.],
           [1647285., 1646285.],
           [1646285., 1645285.],
           ...,
           [ 881285.,  880285.],
           [ 880285.,  879285.],
           [ 879285.,  878285.]], shape=(770, 2))

  • resistivity
    (z, y, x)
    float64
    nan nan nan nan ... nan nan nan nan

    standard_name :

    log10_resistivity

    long_name :

    Electrical Resistivity on Logarithmic Scale

    units :

    log10_ohm_m

    null_value :

    -9999.99

    valid_range :

    [0.27706546 3.74382997]

    grid_mapping :

    spatial_ref

    array([[[nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            ...,
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan]],
    
           [[nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            ...,
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan]],
    
           [[nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            ...,
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan]],
    
           [[nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            ...,
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan]],
    
           [[nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            ...,
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan],
            [nan, nan, nan, ..., nan, nan, nan]]], shape=(5, 770, 363))

  • magnetic_tmi
    (y, x)
    float32
    nan nan nan nan ... nan nan nan nan

    standard_name :

    total_magnetic_intensity

    long_name :

    Total magnetic intensity, diurnally corrected and filtered

    units :

    nT

    null_value :

    -999999

    valid_range :

    [-12772. 11280.377]

    grid_mapping :

    spatial_ref

    array([[nan, nan, nan, ..., nan, nan, nan],
           [nan, nan, nan, ..., nan, nan, nan],
           [nan, nan, nan, ..., nan, nan, nan],
           ...,
           [nan, nan, nan, ..., nan, nan, nan],
           [nan, nan, nan, ..., nan, nan, nan],
           [nan, nan, nan, ..., nan, nan, nan]],
          shape=(770, 363), dtype=float32)

• Attributes: (3)

  uuid :

  fff129ce-4888-4dd9-98a6-4fb66ceb64b1

  comment :

  <additional details or ancillary information>

  content :

  interpolated resistivity models (3-D depth grid) and magnetic map (2-D)

Save to NetCDF file

d_out = join(data_path, 'tifs.nc')
survey.gs.to_netcdf(d_out)

uuid
title
institution
source
history
references
comment
content
conventions
created_by
type
_FillValue
contractor_project_number
contractor
client
survey_type
survey_area_name
state
country
acquisition_start
acquisition_end
dataset_created
coordinates
_FillValue
time
area
current
frequency
electromagnetic_moment
magnetometer_b_field
electromagnetic_b_field
coordinates
_FillValue
traverse_line_spacing
traverse_line_direction
nominal_terrain_clearance
final_line_kilometers
TRAVERSE LINE NUMBERS
REPEAT LINE NUMBERS
PRE ZERO LINE NUMBERS
POST ZERO LINE NUMBERS
coordinates
_FillValue
aircraft
aircraft_registration
magnetometer
magnetometer_installation
electromagnetic_system
electromagnetic_installation
electromagnetic_coil_orientations
spectrometer_system
spectrometer_installation
radar_altimeter_system
radar_altimeter_sample_rate
laser_altimeter_system
laser_altimeter_sample_rate
navigation_system
navigation_sample_rate
acquisition_system
coordinates
_FillValue
crs_wkt
semi_major_axis
semi_minor_axis
inverse_flattening
reference_ellipsoid_name
longitude_of_prime_meridian
prime_meridian_name
geographic_crs_name
horizontal_datum_name
projected_crs_name
grid_mapping_name
standard_parallel
latitude_of_projection_origin
longitude_of_central_meridian
false_easting
false_northing
authority
wkid
content
comment
type
_FillValue
crs_wkt
semi_major_axis
semi_minor_axis
inverse_flattening
reference_ellipsoid_name
longitude_of_prime_meridian
prime_meridian_name
geographic_crs_name
horizontal_datum_name
projected_crs_name
grid_mapping_name
standard_parallel
latitude_of_projection_origin
longitude_of_central_meridian
false_easting
false_northing
authority
wkid
uuid
comment
content
_FillValue
standard_name
long_name
null_value
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
null_value
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
units
null_value
valid_range
grid_mapping
coordinates
_FillValue
crs_wkt
semi_major_axis
semi_minor_axis
inverse_flattening
reference_ellipsoid_name
longitude_of_prime_meridian
prime_meridian_name
geographic_crs_name
horizontal_datum_name
projected_crs_name
grid_mapping_name
standard_parallel
latitude_of_projection_origin
longitude_of_central_meridian
false_easting
false_northing
authority
wkid
GeoTransform
_FillValue
standard_name
long_name
units
null_value
axis
valid_range
grid_mapping
bounds
standard_name
long_name
units
null_value
valid_range
grid_mapping
_FillValue
standard_name
long_name
units
null_value
axis
valid_range
grid_mapping
bounds
uuid
comment
content
_FillValue
comment
standard_name
long_name
null_value
datum
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
null_value
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
null_value
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
units
null_value
valid_range
grid_mapping
coordinates
_FillValue
standard_name
long_name
units
null_value
valid_range
grid_mapping
coordinates
_FillValue
crs_wkt
semi_major_axis
semi_minor_axis
inverse_flattening
reference_ellipsoid_name
longitude_of_prime_meridian
prime_meridian_name
geographic_crs_name
horizontal_datum_name
projected_crs_name
grid_mapping_name
standard_parallel
latitude_of_projection_origin
longitude_of_central_meridian
false_easting
false_northing
authority
wkid
GeoTransform
_FillValue
comment
standard_name
long_name
units
null_value
axis
positive
datum
valid_range
grid_mapping
bounds
standard_name
long_name
units
null_value
valid_range
grid_mapping
_FillValue
standard_name
long_name
units
null_value
axis
valid_range
grid_mapping
bounds
_FillValue
standard_name
long_name
units
null_value
axis
valid_range
grid_mapping
bounds

Reading back in the GS NetCDF file

new_survey = gspy.open_datatree(d_out)['survey']

Plotting

# Make a map-view plot of a specific data variable, using Xarray's plotter
# In this case, we slice the 3-D resistivity variable along the depth dimension
new_survey['derived_products']["maps"]['resistivity'].plot(col='z', vmax=3, cmap='jet', robust=True)

# Make a map-view plot comparing the different x-y discretization of the two magnetic variables, using Xarray's plotter
plt.figure()
ax=plt.gca()
new_survey['derived_products']["maps"]['magnetic_tmi'].plot(ax=ax, cmap='jet', robust=True)
new_survey['derived_products']["map"]['magnetic_tmi'].plot(ax=ax, cmap='Greys', cbar_kwargs={'label': ''}, robust=True)
plt.ylim([1.20556e6, 1.21476e6])
plt.xlim([3.5201e5, 3.6396e5])
plt.show()

• 
•