gdptools example without the OPeNDAP catalog - gridMET data

import logging
import warnings

import hvplot.pandas
import hvplot.xarray
import numpy as np
import pandas as pd
import xarray as xr
from pynhd import NLDI, WaterData
from pygeohydro import watershed

from gdptools import AggGen, UserCatData, WeightGen

warnings.filterwarnings("ignore")

Create GeoDataFrame of Delaware River Basin HUC12s

In this cell we use clients to several https://api.water.usgs.gov tools:

• Tool: USGS NLDI Client: pynhd.NLDI

• Tool: USGS GeoServer Client: pynhd.WaterData

More examples of using the HyRiver pynhd package can be found here

# Delaware River Basin HUC12s
# --------------------------------
wbd = watershed.WBD("huc4")
delaware_basin = wbd.byids(field="huc4", fids="0204")
huc12_basins = WaterData("wbd12").bygeom(delaware_basin.iloc[0].geometry)
huc12_basins = huc12_basins[huc12_basins["huc12"].str.startswith(("020401", "020402"))]
huc12_basins.head()

	geometry	tnmid	metasourceid	sourcedatadesc	sourceoriginator	sourcefeatureid	loaddate	gnis_id	areaacres	areasqkm	states	huc12	name	hutype	humod	tohuc	noncontributingareaacres	noncontributingareasqkm	globalid
0	MULTIPOLYGON (((-75.02421 39.44397, -75.02353 ...	{001055CD-7B5E-4F96-ABB1-22923FD2B454}	NaN	None	None	None	2013-01-18T07:08:41Z	None	16178.10	65.47	NJ	020402060505	White Marsh Run-Maurice River	S	NM	020402060507	0	0	{A70A97FA-E29C-11E2-8094-0021280458E6}
1	MULTIPOLYGON (((-74.62113 41.00565, -74.62139 ...	{009DA440-393A-46A0-BBA3-C38139823414}	{ED602145-9201-4827-9CE1-05D252484579}	None	None	None	2017-10-03T20:11:06Z	None	29107.93	117.80	NJ	020401050502	Lubbers Run-Musconetcong River	S	NM	020401050503	0	0	{B830B4F8-E29C-11E2-8094-0021280458E6}
2	MULTIPOLYGON (((-75.75805 40.65996, -75.75783 ...	{00F73634-FA00-464F-93F0-F0BCCB464567}	{ED602145-9201-4827-9CE1-05D252484579}	None	None	None	2017-10-03T20:10:58Z	None	35229.84	142.57	PA	020402030305	Upper Maiden Creek	S	NM	020402030307	0	0	{BAF5E77A-E29C-11E2-8094-0021280458E6}
3	MULTIPOLYGON (((-75.13279 39.88601, -75.13439 ...	{017B8A58-E959-4600-B212-8EB921F0F9F7}	{E9F5C988-2313-440E-A05E-C5871E2773A6}	None	None	None	2017-10-03T20:11:04Z	None	24920.90	100.85	NJ,PA	020402020507	Woodbury Creek-Delaware River	S	TF	020402020607	0	0	{B122D6B7-E29C-11E2-8094-0021280458E6}
4	MULTIPOLYGON (((-74.92817 40.55256, -74.92833 ...	{01B4598B-0623-407F-9133-3040D99A8D1A}	{ED602145-9201-4827-9CE1-05D252484579}	None	None	None	2017-10-03T20:10:52Z	None	15086.80	61.05	NJ	020401050905	Lockatong Creek	S	NM	020401050908	0	0	{A6AB04BE-E29C-11E2-8094-0021280458E6}

from holoviews.element.tiles import EsriTerrain

drb = huc12_basins.hvplot(
    geo=True, coastline='50m', alpha=0.2, c='r', frame_width=300,
    xlabel="longitude", ylabel="latitude", clabel="monthly water evaporation (mm)",
    title="Delaware River HUC12 basins", xlim=(-76.8, -73.8), aspect='equal'
)
EsriTerrain() * drb

Load gridMET OPeNDAP endpoint with xarray

• Note currently the gridMET data require ‘#fillmismatch’ appended to the end of the OPeNDAP endpoint. For example

gridmet_URL = [
    'http://thredds.northwestknowledge.net:8080/thredds/dodsC/agg_met_tmmx_1979_CurrentYear_CONUS.nc#fillmismatch',
    'http://thredds.northwestknowledge.net:8080/thredds/dodsC/agg_met_tmmn_1979_CurrentYear_CONUS.nc#fillmismatch'
]
data = xr.open_mfdataset(gridmet_URL)
data

<xarray.Dataset> Size: 224GB
Dimensions:                    (day: 17273, lat: 585, lon: 1386, crs: 1)
Coordinates:
  * day                        (day) datetime64[ns] 138kB 1979-01-01 ... 2026...
  * lat                        (lat) float64 5kB 49.4 49.36 ... 25.11 25.07
  * lon                        (lon) float64 11kB -124.8 -124.7 ... -67.1 -67.06
  * crs                        (crs) float32 4B 3.0
Data variables:
    daily_maximum_temperature  (day, lat, lon) float64 112GB dask.array<chunksize=(17273, 585, 1386), meta=np.ndarray>
    daily_minimum_temperature  (day, lat, lon) float64 112GB dask.array<chunksize=(17273, 585, 1386), meta=np.ndarray>
Attributes: (12/19)
    geospatial_bounds_crs:      EPSG:4326
    Conventions:                CF-1.0
    geospatial_bounds:          POLYGON((-124.7666666333333 49.40000000000000...
    geospatial_lat_min:         25.066666666666666
    geospatial_lat_max:         49.40000000000000
    geospatial_lon_min:         -124.7666666333333
    ...                         ...
    date:                       03 March 2026
    note1:                      The projection information for this file is: ...
    note2:                      Citation: Abatzoglou, J.T., 2013, Development...
    note3:                      Data in slices after last_permanent_slice (1-...
    note4:                      Data in slices after last_provisional_slice (...
    note5:                      Days correspond approximately to calendar day...

xarray.Dataset

• Dimensions:
  • day: 17273
  • lat: 585
  • lon: 1386
  • crs: 1
• Coordinates: (4)
  • day
    (day)
    datetime64[ns]
    1979-01-01 ... 2026-04-16

    description :

    days since 1900-01-01

    long_name :

    time

    standard_name :

    time

    _ChunkSizes :

    365

    array(['1979-01-01T00:00:00.000000000', '1979-01-02T00:00:00.000000000',
           '1979-01-03T00:00:00.000000000', ..., '2026-04-14T00:00:00.000000000',
           '2026-04-15T00:00:00.000000000', '2026-04-16T00:00:00.000000000'],
          shape=(17273,), dtype='datetime64[ns]')

  • lat
    (lat)
    float64
    49.4 49.36 49.32 ... 25.11 25.07

    units :

    degrees_north

    description :

    latitude

    long_name :

    latitude

    standard_name :

    latitude

    axis :

    Y

    _ChunkSizes :

    585

    array([49.4     , 49.358333, 49.316667, ..., 25.15    , 25.108333, 25.066667],
          shape=(585,))

  • lon
    (lon)
    float64
    -124.8 -124.7 ... -67.1 -67.06

    units :

    degrees_east

    description :

    longitude

    long_name :

    longitude

    standard_name :

    longitude

    axis :

    X

    _ChunkSizes :

    1386

    array([-124.766667, -124.725   , -124.683333, ...,  -67.141667,  -67.1     ,
            -67.058333], shape=(1386,))

  • crs
    (crs)
    float32
    3.0

    grid_mapping_name :

    latitude_longitude

    longitude_of_prime_meridian :

    0.0

    semi_major_axis :

    6378137.0

    long_name :

    WGS 84

    inverse_flattening :

    298.257223563

    GeoTransform :

    -124.7666666333333 0.041666666666666 0 49.400000000000000 -0.041666666666666

    spatial_ref :

    GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]

    _ChunkSizes :

    1

    array([3.], dtype=float32)

• Data variables: (2)
  • daily_maximum_temperature
    (day, lat, lon)
    float64
    dask.array<chunksize=(17273, 585, 1386), meta=np.ndarray>

    units :

    K

    description :

    Daily Maximum Temperature (2m)

    long_name :

    tmmx

    standard_name :

    tmmx

    dimensions :

    lon lat time

    grid_mapping :

    crs

    coordinate_system :

    WGS84,EPSG:4326

    _ChunkSizes :

    [ 61 98 231]

    Array Chunk Bytes 104.35 GiB 104.35 GiB Shape (17273, 585, 1386) (17273, 585, 1386) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

  • daily_minimum_temperature
    (day, lat, lon)
    float64
    dask.array<chunksize=(17273, 585, 1386), meta=np.ndarray>

    units :

    K

    description :

    Daily Minimum Temperature (2m)

    long_name :

    tmmn

    standard_name :

    tmmn

    dimensions :

    lon lat time

    grid_mapping :

    crs

    coordinate_system :

    WGS84,EPSG:4326

    _ChunkSizes :

    [ 61 98 231]

    Array Chunk Bytes 104.35 GiB 104.35 GiB Shape (17273, 585, 1386) (17273, 585, 1386) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

• Attributes: (19)

  geospatial_bounds_crs :

  EPSG:4326

  Conventions :

  CF-1.0

  geospatial_bounds :

  POLYGON((-124.7666666333333 49.400000000000000, -124.7666666333333 25.066666666666666, -67.058333300000015 25.066666666666666, -67.058333300000015 49.400000000000000, -124.7666666333333 49.400000000000000))

  geospatial_lat_min :

  25.066666666666666

  geospatial_lat_max :

  49.40000000000000

  geospatial_lon_min :

  -124.7666666333333

  geospatial_lon_max :

  -67.058333300000015

  geospatial_lon_resolution :

  0.041666666666666

  geospatial_lat_resolution :

  0.041666666666666

  geospatial_lat_units :

  decimal_degrees north

  geospatial_lon_units :

  decimal_degrees east

  coordinate_system :

  EPSG:4326

  author :

  John Abatzoglou - University of California Merced, jabatzoglou@ucmerced.edu

  date :

  03 March 2026

  note1 :

  The projection information for this file is: GCS WGS 1984.

  note2 :

  Citation: Abatzoglou, J.T., 2013, Development of gridded surface meteorological data for ecological applications and modeling, International Journal of Climatology, DOI: 10.1002/joc.3413

  note3 :

  Data in slices after last_permanent_slice (1-based) are considered provisional and subject to change with subsequent updates

  note4 :

  Data in slices after last_provisional_slice (1-based) are considered early and subject to change with subsequent updates

  note5 :

  Days correspond approximately to calendar days ending at midnight, Mountain Standard Time (7 UTC the next calendar day)

Calculate weights

• To determine an area-weighted statistic for each HUC12 polygon, the area of all grid cells (a Shapely Polygon object representing each cell of the grid).

• The weight of each cell contributing to a HUC12 polygon is the cell-intersected-area / HUC12 polygon-area.

• Not all all data readable by xarray are made the same. Therefor, when using non-OPeNDAP catalog endpoints, the user must supply some metadata about the xarray data and GeoDataFrame data.

Enable logging

gdptools uses Python’s standard logging module. By default it emits no output. Configuring the gdptools logger before running weight generation and aggregation lets you monitor progress, timing, and data dimensions throughout the workflow. Here we set gdptools to INFO while keeping third-party libraries at WARNING to avoid noisy output from rasterio, fiona, etc.

# Configure logging: gdptools at INFO, everything else at WARNING
logging.basicConfig(
    level=logging.WARNING,
    format="%(name)s | %(levelname)s | %(message)s",
    force=True,
)
logging.getLogger("gdptools").setLevel(logging.INFO)

data['crs'].attrs.get('spatial_ref')

'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'

# metadata for calculating weights

data_crs = data['crs'].attrs.get('spatial_ref')
x_coord = "lon"
y_coord = "lat"
t_coord = "day"
sdate = "1979-01-01"
edate = "1979-01-07"
var = ["daily_maximum_temperature", "daily_minimum_temperature"]
shp_crs = 4326
shp_poly_idx = "huc12"
wght_gen_crs = 6931

user_data = UserCatData(
    source_ds=data,
    source_crs=data_crs,
    source_x_coord=x_coord,
    source_y_coord=y_coord,
    source_t_coord=t_coord,
    source_var=var,
    target_gdf=huc12_basins,
    target_crs=shp_crs,
    target_id=shp_poly_idx,
    source_time_period=[sdate, edate]
)

wght_gen = WeightGen(
    user_data=user_data,
    method="serial",
    output_file='example_weights.csv',
    weight_gen_crs=6931,
)

wdf = wght_gen.calculate_weights()

To illustrate the format of the saved weights file it’s opened here, however in the next step we use the returned datafame.

pd.read_csv('example_weights.csv')

	huc12	i	j	wght
0	20402060505	74	35	0.081887
1	20402060505	77	35	0.002041
2	20402060505	76	35	0.168929
3	20402060505	75	35	0.201884
4	20402060505	76	34	0.119210
...	...	...	...	...
5677	20401050103	37	36	0.032199
5678	20401050103	36	37	0.184037
5679	20401050103	35	37	0.100861
5680	20401050103	34	37	0.000047
5681	20401050103	36	36	0.023088

5682 rows × 4 columns

Run the grid-to-polygon area weighted average interpolation using the calculated weights.

• Future releases will add additional area-weighted statistics. Currently only the mean is returned.

• If more that one variable were available in the dataset then the run_weights function could be embedded in a loop that indexes each variable.

run_weights Output

Run weighs returns:

• ngdf is the GeoDataFrame passed to run_weights. huc12_basins maybe reordered when running weights so the re-ordered GeoDataFrame is returned.

• vals are the area-weighted interpolated values ordered by the ngdf.index

agg_gen = AggGen(
    user_data=user_data,
    stat_method="masked_mean",
    agg_engine="serial",
    agg_writer="csv",
    weights='example_weights.csv',
    out_path='.',
    file_prefix="test_gm_nc"
)
ngdf, ds_out = agg_gen.calculate_agg()

ds_out

<xarray.Dataset> Size: 58kB
Dimensions:                    (time: 7, huc12: 427)
Coordinates:
  * time                       (time) datetime64[ns] 56B 1979-01-01 ... 1979-...
  * huc12                      (huc12) object 3kB '020401010101' ... '0204020...
    lat                        (huc12) float64 3kB 42.39 42.38 ... 38.76 38.83
    lon                        (huc12) float64 3kB -74.62 -74.71 ... -75.18
Data variables:
    daily_maximum_temperature  (time, huc12) float64 24kB 283.3 283.8 ... 283.8
    crs                        float64 8B 1.0
    daily_minimum_temperature  (time, huc12) float64 24kB 276.9 277.3 ... 273.5
Attributes:
    Conventions:  CF-1.8
    featureType:  timeSeries
    history:      2026_04_18_03_42_45 Original file created by gdptools packa...

xarray.Dataset

• Dimensions:
  • time: 7
  • huc12: 427
• Coordinates: (4)
  • time
    (time)
    datetime64[ns]
    1979-01-01 ... 1979-01-07

    standard_name :

    time

    axis :

    T

    description :

    days since 1900-01-01

    long_name :

    time

    _ChunkSizes :

    365

    array(['1979-01-01T00:00:00.000000000', '1979-01-02T00:00:00.000000000',
           '1979-01-03T00:00:00.000000000', '1979-01-04T00:00:00.000000000',
           '1979-01-05T00:00:00.000000000', '1979-01-06T00:00:00.000000000',
           '1979-01-07T00:00:00.000000000'], dtype='datetime64[ns]')

  • huc12
    (huc12)
    object
    '020401010101' ... '020402070605'

    feature_id :

    huc12

    array(['020401010101', '020401010102', '020401010103', ..., '020402070603',
           '020402070604', '020402070605'], shape=(427,), dtype=object)

  • lat
    (huc12)
    float64
    42.39 42.38 42.33 ... 38.76 38.83

    long_name :

    Latitude of HRU centroid

    units :

    degrees_north

    standard_name :

    latitude

    axis :

    Y

    array([42.38709124, 42.38194564, 42.33066568, 42.34506257, 42.27096486,
           42.33286078, 42.27268526, 42.21477835, 42.23599231, 42.16431883,
           42.19773905, 42.19367479, 42.10830645, 42.15231851, 42.14804815,
           42.0610358 , 42.06494591, 42.00364181, 42.01807534, 41.92807218,
           41.82020995, 41.8177342 , 41.88759791, 41.89359282, 41.79416696,
           41.88143429, 41.81074205, 41.8325948 , 41.79741774, 41.7579096 ,
           41.75263053, 41.72747402, 41.70039575, 41.64945824, 41.63989504,
           41.63528063, 41.56591882, 41.52355406, 41.8601189 , 41.94431143,
           41.93807835, 41.91149244, 42.02635622, 41.99380184, 41.93423675,
           41.97821095, 42.24438706, 42.28232371, 42.19237914, 42.09355086,
           42.2066441 , 42.14669668, 42.08038278, 42.17847768, 42.10887865,
           42.15423287, 42.07162316, 42.12396886, 42.06025285, 42.04818976,
           42.04152554, 41.96271441, 42.02238788, 41.96134086, 41.7444889 ,
           41.60007067, 41.66894526, 41.72628504, 41.72767136, 41.65829199,
           41.51024501, 41.51510833, 41.44041283, 41.37294254, 41.28511086,
           41.31913834, 41.28472112, 41.42847975, 41.40974851, 41.5768005 ,
           41.39814576, 41.48152051, 41.79106472, 41.75846873, 41.76022896,
           41.66332337, 41.63746458, 41.58696315, 41.50037311, 41.61167908,
           41.55134211, 41.52414727, 41.98721994, 41.95663128, 41.87401815,
           41.70078989, 41.74135134, 41.55495864, 41.56037636, 41.38977945,
    ...
           40.01258757, 40.07769698, 40.01770113, 39.97326735, 39.97273175,
           40.03096053, 39.89589448, 39.94976482, 39.93668621, 39.81751694,
           39.77981606, 39.82517868, 39.84958314, 39.85392694, 39.71559434,
           39.79130466, 39.73694403, 39.90072259, 39.83987768, 39.77262659,
           39.6134099 , 39.67609231, 39.66175928, 39.74792837, 39.71299305,
           39.78513135, 39.77710259, 39.68562058, 39.64478301, 39.60199585,
           39.57457983, 39.54110027, 39.49741806, 39.44699068, 39.3889688 ,
           39.70274646, 39.64665701, 39.72981077, 39.6086974 , 39.62224572,
           39.60094347, 39.56590078, 39.5273166 , 39.5073717 , 39.49510458,
           39.40954598, 39.39172393, 39.63264235, 39.64426141, 39.59733826,
           39.5171156 , 39.52402719, 39.441139  , 39.44956015, 39.44380073,
           39.3467356 , 39.38872798, 39.40158597, 39.28903075, 39.47654928,
           39.45578049, 39.40144456, 39.3133119 , 39.21451082, 39.03033286,
           39.16548315, 39.16780381, 39.29300703, 39.34802074, 39.38971414,
           39.2371201 , 39.24609035, 39.29600579, 39.20335525, 39.14166802,
           39.17037489, 39.11958873, 39.08468819, 39.08814104, 38.98365802,
           38.93646378, 39.03847816, 38.98125807, 39.02215901, 39.04138225,
           38.8888352 , 38.83920929, 38.92373079, 38.87063938, 38.89500112,
           38.95868477, 38.95182021, 38.74409774, 38.8183628 , 38.76953382,
           38.76308549, 38.82734413])

  • lon
    (huc12)
    float64
    -74.62 -74.71 ... -75.2 -75.18

    long_name :

    Longitude of HRU centroid

    units :

    degrees_east

    standard_name :

    longitude

    axis :

    X

    array([-74.62155934, -74.71393632, -74.71097823, -74.82334628,
           -74.78436641, -74.9050593 , -74.982048  , -74.94972416,
           -75.07046025, -75.05487613, -75.14871651, -75.27221166,
           -75.25187928, -75.37734435, -75.46499607, -75.50578141,
           -75.41324413, -75.39714836, -75.29989383, -75.35319139,
           -75.35504196, -75.28379964, -75.27209008, -75.06522679,
           -75.18527034, -75.16809027, -75.08286154, -74.97627653,
           -74.87211116, -74.92359403, -74.9951963 , -75.08496109,
           -75.14992347, -75.14852115, -74.96700246, -75.05089549,
           -75.11226437, -75.00260609, -74.76993575, -74.64269853,
           -74.74286003, -74.86972929, -74.66558223, -74.83597523,
           -75.03836962, -74.99689157, -74.53846443, -74.58432546,
           -74.50728543, -74.55079353, -74.69155338, -74.65326962,
           -74.65238065, -74.80481052, -74.76679719, -74.89627206,
           -74.87883785, -74.99265226, -75.02212726, -75.11164318,
           -75.18426846, -75.14992755, -75.25357205, -75.22204532,
           -75.39204394, -75.38822378, -75.37595071, -75.26653723,
           -75.33089016, -75.25811999, -75.40956468, -75.2948941 ,
           -75.43747734, -75.43356639, -75.25358723, -75.22234631,
           -75.3372204 , -75.34772041, -75.26205343, -75.21616626,
    ...
           -75.58039055, -75.56309396, -75.73818909, -75.77393422,
           -75.66475678, -75.59799505, -75.57469113, -75.48695165,
           -75.52264116, -75.55529616, -75.61375393, -75.65002735,
           -75.65619833, -75.64664153, -75.61775882, -75.6773595 ,
           -75.640833  , -75.41428481, -75.33325574, -75.46383287,
           -75.38803708, -75.53296685, -75.468395  , -75.32641408,
           -75.41773566, -75.48197819, -75.26006972, -75.20690026,
           -75.32115836, -75.03572653, -75.10827543, -75.18014582,
           -75.1648013 , -75.05421537, -75.13225795, -75.04151253,
           -74.96886142, -75.08947663, -75.0538509 , -74.92241422,
           -74.97670576, -75.34450956, -75.45149103, -75.391126  ,
           -75.24124881, -75.11145662, -74.92115948, -74.85756398,
           -74.97804183, -75.63691036, -75.58186614, -75.54043255,
           -75.61363294, -75.50650617, -75.45470277, -75.43581167,
           -75.41557127, -75.59035423, -75.58455701, -75.55923319,
           -75.43367018, -75.5768933 , -75.5411715 , -75.52643162,
           -75.48143643, -75.42474156, -75.34191737, -75.48741851,
           -75.41582301, -75.42282182, -75.30484825, -75.35384928,
           -75.36063472, -75.28481889, -75.34372597, -75.31304137,
           -75.2600868 , -75.19528858, -75.17671553])

• Data variables: (3)
  • daily_maximum_temperature
    (time, huc12)
    float64
    283.3 283.8 283.9 ... 284.2 283.8

    units :

    K

    long_name :

    tmmx

    coordinates :

    time lat lon

    grid_mapping :

    crs

    array([[283.34210068, 283.78952008, 283.89178203, ..., 287.57650179,
            286.57301094, 285.8546885 ],
           [282.64553906, 282.96178847, 283.04012477, ..., 286.07571275,
            285.56497204, 285.07155596],
           [270.35466739, 270.31726438, 270.10214523, ..., 278.55367893,
            278.52993597, 278.16978482],
           ...,
           [264.06583868, 264.05605069, 263.84202038, ..., 277.05952989,
            277.26357088, 277.05864202],
           [269.96373084, 270.23207355, 270.22019487, ..., 278.40812795,
            278.69958657, 278.54754627],
           [275.40301559, 275.93843299, 276.21703238, ..., 284.42459739,
            284.17531281, 283.79556951]], shape=(7, 427))

  • crs
    ()
    float64
    1.0

    crs_wkt :

    GEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],USAGE[SCOPE["Horizontal component of 3D system."],AREA["World."],BBOX[-90,-180,90,180]],ID["EPSG",4326]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984 ensemble

    grid_mapping_name :

    latitude_longitude

    array(1.)

  • daily_minimum_temperature
    (time, huc12)
    float64
    276.9 277.3 277.2 ... 273.2 273.5

    units :

    K

    long_name :

    tmmn

    coordinates :

    time lat lon

    grid_mapping :

    crs

    array([[276.87512065, 277.31779651, 277.20568237, ..., 280.42187048,
            280.03005589, 280.07989482],
           [269.25666708, 269.413816  , 269.07031036, ..., 275.42378863,
            275.36636168, 275.48730845],
           [257.50149145, 257.83603238, 257.34029282, ..., 263.49835217,
            263.64024802, 263.86974059],
           ...,
           [259.29726844, 259.59266588, 259.18771484, ..., 268.4642002 ,
            269.10390725, 269.62926068],
           [261.07456999, 261.39755945, 261.16880077, ..., 272.1326193 ,
            272.37994001, 272.60068553],
           [262.67284447, 262.95187546, 262.71591502, ..., 273.02222678,
            273.22654557, 273.48730845]], shape=(7, 427))

• Attributes: (3)

  Conventions :

  CF-1.8

  featureType :

  timeSeries

  history :

  2026_04_18_03_42_45 Original file created by gdptools package: https://code.usgs.gov/wma/nhgf/toolsteam/gdptools

pd.read_csv('test_gm_nc.csv', index_col=0)

	time	varname	units	020401010101	020401010102	020401010103	020401010104	020401010105	020401010106	020401010201	...	020402070503	020402070504	020402070505	020402070506	020402070507	020402070601	020402070602	020402070603	020402070604	020402070605
index
0	1979-01-01	daily_maximum_temperature	K	283.342101	283.789520	283.891782	284.313867	284.279681	283.724576	283.688328	...	288.614109	287.102319	287.257388	287.192510	285.911595	288.647165	287.673752	287.576502	286.573011	285.854688
1	1979-01-02	daily_maximum_temperature	K	282.645539	282.961788	283.040125	283.368031	283.418329	282.748207	282.790915	...	287.766527	285.970048	286.572958	286.722795	285.626663	286.617455	286.196586	286.075713	285.564972	285.071556
2	1979-01-03	daily_maximum_temperature	K	270.354667	270.317264	270.102145	270.433602	270.454707	270.254940	270.760374	...	278.432881	278.402124	278.703507	278.566303	278.221036	278.464381	278.479914	278.553679	278.529936	278.169785
3	1979-01-04	daily_maximum_temperature	K	263.918458	263.961688	263.736135	264.144980	264.093689	264.068521	264.586014	...	274.790052	274.646453	275.191191	275.202789	274.948669	274.324795	274.502752	274.522852	274.589280	274.395581
4	1979-01-05	daily_maximum_temperature	K	264.065839	264.056051	263.842020	264.165891	264.218206	264.084139	264.706366	...	277.148031	277.046858	277.584110	277.602286	277.298743	276.941864	276.956516	277.059530	277.263571	277.058642
5	1979-01-06	daily_maximum_temperature	K	269.963731	270.232074	270.220195	270.616740	270.654793	270.416690	270.891040	...	278.232881	278.396659	278.907483	278.801530	278.666241	278.123928	278.283779	278.408128	278.699587	278.547546
6	1979-01-07	daily_maximum_temperature	K	275.403016	275.938433	276.217032	276.513022	276.740472	275.976119	276.418405	...	284.271129	284.194942	284.334873	284.167058	283.541393	284.590986	284.436977	284.424597	284.175313	283.795570
0	1979-01-01	daily_minimum_temperature	K	276.875121	277.317797	277.205682	277.407442	277.398712	277.034636	276.912089	...	281.526809	280.340967	280.397078	280.529055	280.087390	280.910730	280.530568	280.421870	280.030056	280.079895
1	1979-01-02	daily_minimum_temperature	K	269.256667	269.413816	269.070310	269.239835	269.257574	269.118803	269.312846	...	274.801765	275.134578	274.757241	274.653814	274.893253	275.437823	275.296745	275.423789	275.366362	275.487308
2	1979-01-03	daily_minimum_temperature	K	257.501491	257.836032	257.340293	257.918269	257.588494	258.058366	258.241108	...	263.393736	263.564219	263.443779	263.438990	263.702673	263.170775	263.485903	263.498352	263.640248	263.869741
3	1979-01-04	daily_minimum_temperature	K	258.095835	258.408270	257.959467	258.499743	258.186048	258.647050	258.789096	...	263.532182	263.871596	263.858109	263.911718	264.327853	263.159602	263.670639	263.696368	264.030124	264.398075
4	1979-01-05	daily_minimum_temperature	K	259.297268	259.592666	259.187715	259.580379	259.406128	259.667843	259.853056	...	268.417116	268.802196	269.038191	269.138270	269.711928	267.564882	268.407991	268.464200	269.103907	269.629261
5	1979-01-06	daily_minimum_temperature	K	261.074570	261.397559	261.168801	261.430951	261.427480	261.406271	261.782945	...	271.789742	272.279462	272.127984	272.111718	272.457381	271.686187	272.149811	272.132619	272.379940	272.600686
6	1979-01-07	daily_minimum_temperature	K	262.672844	262.951875	262.715915	262.921207	263.001742	262.820974	263.269528	...	272.644310	273.099901	273.012330	272.954145	273.317758	272.646589	272.972866	273.022227	273.226546	273.487308

14 rows × 430 columns

aggdata = agg_gen.agg_data
ds_ss = aggdata.get("daily_maximum_temperature").da

Plot the gridded data and interpolated data

var = "daily_maximum_temperature"
time_step = "1979-01-07"
tmin = np.nanmin(ds_ss.sel(day=time_step).values)
tmax = np.nanmax(ds_ss.sel(day=time_step).values)
print(tmin,tmax)

273.1 286.1

ds_out[var].shape

(7, 427)

ngdf[var] = ds_out[var].sel(time=time_step)
gridmet_agg = ngdf.to_crs(4326).hvplot(
    c='daily_maximum_temperature', geo=True, coastline='50m', frame_width=300, alpha=1, cmap='viridis', clim=(tmin, tmax),
    xlabel="longitude", ylabel="latitude", clabel="daily_maximum_temperature (K)", xlim=(-76.8, -73.8),
    title="HUC12 area-weighted average", aspect='equal'
)
gridmet_raw = ds_ss.isel(day=6).hvplot.quadmesh(
    'lon', 'lat', 'daily_maximum_temperature', cmap='viridis', alpha=1, grid=True, geo=True, coastline='50m', frame_width=300, clim=(tmin, tmax),
    xlabel="longitude", ylabel="latitude", clabel="daily_maximum_temperature (K)", xlim=(-76.8, -73.8),
    title="Gridmet daily max temperature", aspect='equal'
)

(EsriTerrain() * gridmet_raw) + (EsriTerrain() * gridmet_agg) + (EsriTerrain() * gridmet_raw * gridmet_agg)