# Copyright (c) 2023, TU Wien, Department of Geodesy and Geoinformation
# All rights reserved.
# Redistribution and use in source and binary forms, with or without
# modification, are permitted provided that the following conditions are met:
# * Redistributions of source code must retain the above copyright
# notice, this list of conditions and the following disclaimer.
# * Redistributions in binary form must reproduce the above copyright
# notice, this list of conditions and the following disclaimer in the
# documentation and/or other materials provided with the distribution.
# * Neither the name of TU Wien, Department of Geodesy and Geoinformation
# nor the names of its contributors may be used to endorse or promote
# products derived from this software without specific prior written
# permission.
# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
# ARE DISCLAIMED. IN NO EVENT SHALL TU WIEN, DEPARTMENT OF GEODESY AND
# GEOINFORMATION BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
# EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
# PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
# OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
# WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
# OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
# ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
import os
import abc
import glob
import copy
import warnings
from datetime import datetime
import numpy as np
from pygeobase.utils import split_daterange_in_intervals
from pygeobase.object_base import Image
[docs]
class StaticBase:
"""
The StaticBase class serves as a template for i/o objects used in
GriddedStaticBase.
Parameters
----------
filename : str
File name.
mode : str, optional
Opening mode. Default: r
"""
__metaclass__ = abc.ABCMeta
def __init__(self, filename, mode='r', **kwargs):
"""
Initialization of i/o object.
"""
self.filename = filename
self.mode = mode
self.kwargs = kwargs
[docs]
@abc.abstractmethod
def read(self, gpi):
"""
Read data for given grid point.
Parameters
----------
gpi : int
Grid point index.
Returns
-------
data : numpy.ndarray
Data set.
"""
return
[docs]
@abc.abstractmethod
def write(self, data):
"""
Write data.
Parameters
----------
data : numpy.ndarray
Data records.
"""
return
[docs]
@abc.abstractmethod
def flush(self):
"""
Flush data.
"""
return
[docs]
@abc.abstractmethod
def close(self):
"""
Close file.
"""
return
[docs]
class TsBase:
"""
The TsBase class serves as a template for i/o objects used in
GriddedTsBase.
Parameters
----------
filename : str
File name.
mode : str, optional
Opening mode. Default: r
"""
__metaclass__ = abc.ABCMeta
def __init__(self, filename, mode='r', **kwargs):
"""
Initialization of i/o object.
"""
self.filename = filename
self.mode = mode
self.kwargs = kwargs
[docs]
@abc.abstractmethod
def read(self, gpi, **kwargs):
"""
Read time series data for given grid point.
Parameters
----------
gpi : int
Grid point index.
Returns
-------
data : object
pygeobase.object_base.TS object.
"""
return
[docs]
@abc.abstractmethod
def write(self, gpi, data, **kwargs):
"""
Write data.
Parameters
----------
gpi : int
Grid point index.
data : object
pygeobase.object_base.TS object.
"""
return
[docs]
def flush(self):
"""
Flush data.
"""
return
[docs]
def close(self):
"""
Close file.
"""
return
[docs]
class ImageBase:
"""
ImageBase class serves as a template for i/o objects used for reading
and writing image data.
Parameters
----------
filename : str
Filename path.
mode : str, optional
Opening mode. Default: r
"""
__metaclass__ = abc.ABCMeta
def __init__(self, filename, mode='r', **kwargs):
"""
Initialization of i/o object.
"""
self.filename = filename
self.mode = mode
self.kwargs = kwargs
[docs]
@abc.abstractmethod
def read(self, **kwargs):
"""
Read data of an image file.
Returns
-------
image : object
pygeobase.object_base.Image object
"""
return
[docs]
def read_masked_data(self, **kwargs):
"""
Read data of an image file and mask the data according to
specifications.
Returns
-------
image : object
pygeobase.object_base.Image object
"""
raise NotImplementedError('Please implement to enable.')
[docs]
def resample_data(self, image, index, distance, weights, **kwargs):
"""
Takes an image and resample (interpolate) the image data to
arbitrary defined locations given by index and distance.
The default implementation just takes the weighted mean of
all defined distances.
Parameters
----------
image : :py:class`pygeobase.object_base.Image` or numpy.recarray
Image or numpy.recarray like object with shape = (x, )
index : np.array
Index into image data defining a look-up table for data elements
used in the interpolation process for each defined target
location. For each point in image the neighbors in the targed
grid are in the index array. This array is of shape (x, max_neighbors)
distance : np.array
Array representing the distances of the image data to the
arbitrary defined locations.
The distances of points not to use are set to np.inf
This array is of shape (x, max_neighbors)
weights : np.array
Array representing the weights of the image data that should be
used during resampling.
The weights of points not to use are set to np.nan
This array is of shape (x, max_neighbors)
Returns
-------
target : dict
dictionary with a numpy.ndarray for each field in
the input image. We can not return a image here
since we do not know the target latitudes and longitudes.
"""
total_weights = np.nansum(weights, axis=1)
target = {}
for name in image.dtype.names:
target[name] = np.nansum(image[name][index] * weights,
axis=1) / total_weights
return target
[docs]
@abc.abstractmethod
def write(self, image, **kwargs):
"""
Write data to an image file.
Parameters
----------
image : object
pygeobase.object_base.Image object
"""
return
[docs]
@abc.abstractmethod
def flush(self):
"""
Flush data.
"""
return
[docs]
@abc.abstractmethod
def close(self):
"""
Close file.
"""
return
[docs]
class GriddedBase:
"""
The GriddedBase class uses another IO class together with a grid
object to read/write a dataset under the given path.
Parameters
----------
path : string
Path to dataset.
grid : pygeogrids.BasicGrid of CellGrid instance
Grid on which the time series data is stored.
ioclass : class
IO class.
mode : str, optional
File mode and can be read 'r', write 'w' or append 'a'. Default: 'r'
fn_format : str, optional
The string format of the cell files. Default: '{:04d}'
ioclass_kws : dict, optional
Additional keyword arguments for the ioclass. Default: None
"""
__metaclass__ = abc.ABCMeta
def __init__(self,
path,
grid,
ioclass,
mode='r',
fn_format='{:04d}',
ioclass_kws=None):
self.path = path
self.grid = grid
self.ioclass = ioclass
self.mode = mode
self.fn_format = fn_format
self.previous_cell = None
self.fid = None
if ioclass_kws is None:
self.ioclass_kws = {}
else:
self.ioclass_kws = ioclass_kws
def __enter__(self):
"""
Context manager initialization.
Returns
-------
self : GriddedStaticBase object
self
"""
return self
def __exit__(self, exc_type, exc_value, traceback):
"""
Exit the runtime context related to this object. The file will be
closed. The parameters describe the exception that caused the
context to be exited.
exc_type :
exc_value :
traceback :
"""
self.close()
def _open(self, gp):
"""
Open file.
Parameters
----------
gp : int
Grid point.
Returns
-------
success : boolean
Flag if opening the file was successful.
"""
success = True
cell = self.grid.gpi2cell(gp)
filename = os.path.join(self.path, self.fn_format.format(cell))
if self.mode == 'r':
if self.previous_cell != cell:
self.close()
try:
self.fid = self.ioclass(filename,
mode=self.mode,
**self.ioclass_kws)
except IOError as e:
success = False
self.fid = None
msg = f"I/O error({e.errno}): {e.strerror}, {filename}"
warnings.warn(msg, RuntimeWarning)
self.previous_cell = None
else:
self.previous_cell = cell
if self.mode in ['w', 'a']:
if self.previous_cell != cell:
self.flush()
self.close()
try:
self.fid = self.ioclass(filename,
mode=self.mode,
**self.ioclass_kws)
except IOError as e:
success = False
self.fid = None
msg = f"I/O error({e.errno}): {e.strerror}, {filename}"
warnings.warn(msg, RuntimeWarning)
self.previous_cell = None
else:
self.previous_cell = cell
return success
def _read_lonlat(self, lon, lat, **kwargs):
"""
Reading data for given longitude and latitude coordinate.
Parameters
----------
lon : float
Longitude coordinate.
lat : float
Latitude coordinate.
Returns
-------
data : dict of values
data record.
"""
max_dist = kwargs.pop("max_dist", np.Inf)
gp, _ = self.grid.find_nearest_gpi(lon, lat, max_dist)
return self._read_gp(gp, **kwargs)
def _read_gp(self, gp, **kwargs):
"""
Read data for given grid point.
Parameters
----------
gp : int
Grid point.
Returns
-------
data : numpy.ndarray
Data set or None in case of an error.
"""
if self.mode in ['w']:
raise IOError("File is not open in read mode")
data = None
if self._open(gp):
data = self.fid.read(gp, **kwargs)
return data
[docs]
def read(self, *args, **kwargs):
"""
Takes either 1 or 2 arguments and calls the correct function
which is either reading the gpi directly or finding
the nearest gpi from given lat,lon coordinates and then reading it
"""
if len(args) == 1:
data = self._read_gp(args[0], **kwargs)
if len(args) == 2:
data = self._read_lonlat(args[0], args[1], **kwargs)
if len(args) < 1 or len(args) > 2:
raise ValueError("Wrong number of arguments")
return data
def _write_lonlat(self, lon, lat, data, **kwargs):
"""
Write time series for given longitude and latitude coordinate.
Parameters
----------
lon : float
Longitude coordinate.
lat : float
Latitude coordinate.
data : numpy.ndarray
Data records.
"""
gp, _ = self.grid.find_nearest_gpi(lon, lat)
return self._write_gp(gp, data, **kwargs)
[docs]
def write(self, *args, **kwargs):
"""
Takes either 1 or 2 arguments and calls the correct function
which is either writing the gpi directly or finding
the nearest gpi from given lat,lon coordinates and then writing it.
"""
if len(args) == 1:
# args: data
self._write_gp(args[0]['gpi'], args[0], **kwargs)
if len(args) == 2:
# args: gp, data
self._write_gp(args[0], args[1], **kwargs)
if len(args) == 3:
# args: lon, lat, data
self._write_lonlat(args[0], args[1], args[2], **kwargs)
if len(args) < 1 or len(args) > 3:
raise ValueError("Wrong number of arguments")
def _write_gp(self, gp, data, **kwargs):
"""
Write data for given grid point.
Parameters
----------
gp : int
Grid point.
data : numpy.ndarray
Data
"""
if self.mode in ['r']:
raise IOError("File is not open in write/append mode")
if self._open(gp):
lon, lat = self.grid.gpi2lonlat(gp)
self.fid.write(gp, data, lon=lon, lat=lat, **kwargs)
[docs]
def iter_gp(self, **kwargs):
"""
Yield all values for all grid points.
Yields
------
data : pandas.DataFrame
Data set.
gp : int
Grid point.
"""
if 'll_bbox' in kwargs:
latmin, latmax, lonmin, lonmax = kwargs['ll_bbox']
gps = self.grid.get_bbox_grid_points(latmin, latmax, lonmin,
lonmax)
kwargs.pop('ll_bbox', None)
elif 'gpis' in kwargs:
subgrid = self.grid.subgrid_from_gpis(kwargs['gpis'])
gp_info = list(subgrid.grid_points())
gps = np.array(gp_info, dtype=np.int32)[:, 0]
kwargs.pop('gpis', None)
else:
gp_info = list(self.grid.grid_points())
gps = np.array(gp_info, dtype=np.int32)[:, 0]
for gp in gps:
try:
data = self._read_gp(gp, **kwargs)
except IOError as e:
msg = f"I/O error({e.errno}): {e.strerror}, {gp}"
warnings.warn(msg, RuntimeWarning)
data = None
yield data, gp
[docs]
def flush(self):
"""
Flush data.
"""
if self.fid is not None:
self.fid.flush()
[docs]
def close(self):
"""
Close file.
"""
if self.fid is not None:
self.fid.close()
self.fid = None
[docs]
def get_spatial_subset(self,
gpis=None,
cells=None,
ll_bbox=None,
grid=None):
"""
Select spatial subset and return data set with new grid.
Parameters
----------
gpis : numpy.ndarray
Grid point indices.
cells : numpy.ndarray
Cell number.
ll_bbox : tuple (latmin, latmax, lonmin, lonmax)
Lat/Lon bounding box
grid : pygeogrids.CellGrid
Grid object.
Returns
-------
dataset : GriddedBase or child
New data set with for spatial subset.
"""
if gpis:
new_grid = self.grid.subgrid_from_gpis(gpis)
if cells:
new_grid = self.grid.subgrid_from_cells(cells)
if ll_bbox:
latmin, latmax, lonmin, lonmax = ll_bbox
gps = self.grid.get_bbox_grid_points(latmin, latmax, lonmin,
lonmax)
new_grid = self.grid.subgrid_from_gpis(gps)
if grid:
new_grid = grid
dataset = copy.deepcopy(self)
dataset.grid = new_grid
return dataset
[docs]
class GriddedTsBase(GriddedBase):
"""
The GriddedTsBase class uses another IO class together with a grid object
to read/write a time series dataset under the given path.
"""
def _read_gp(self, gp, **kwargs):
"""
Reads time series for a given grid point index.
Parameters
----------
gp : int
Grid point.
Returns
-------
data : object
pygeobase.object_base.TS object
"""
if self.mode in ['w']:
raise IOError("File is not open in read mode")
data = None
if self._open(gp):
data = self.fid.read(gp, **kwargs)
return data
def _write_gp(self, gp, data, **kwargs):
"""
Write data for given grid point.
Parameters
----------
gp : int
Grid point.
data : object
pygeobase.object_base.TS object
"""
if self.mode in ['r']:
raise IOError("File is not open in write/append mode")
if self._open(gp):
lon, lat = self.grid.gpi2lonlat(gp)
self.fid.write(gp, data, lon=lon, lat=lat, **kwargs)
[docs]
class MultiTemporalImageBase:
"""
The MultiTemporalImageBase class make use of an ImageBase object to
read/write a sequence of multi temporal images under a given path.
Parameters
----------
path : string
Path to dataset.
ioclass : class
IO class.
mode : str, optional
File mode and can be read 'r', write 'w' or append 'a'. Default: 'r'
fname_templ : str
Filename template of the data to read. Default placeholder for
parsing datetime information into the fname_templ is "{datetime}".
e.g. "ASCAT_{datetime}_image.nc" will be translated into the filename
ASCAT_20070101_image.nc for the date 2007-01-01.
datetime_format : str
String specifying the format of the datetime object to be parsed
into the fname_template.
e.g. "%Y/%m" will result in 2007/01 for datetime 2007-01-01 12:15:00
subpath_templ : list, optional
If given it is used to generate a sub-paths from the given timestamp.
Each item in the list represents one folder level. This can be used
if the files for May 2007 are e.g. in folders 2007/05/ then the
files can be accessed via the list ['%Y', '%m'].
ioclass_kws : dict
Additional keyword arguments for the ioclass.
exact_templ : boolean, optional
If True then the fname_templ matches the filename exactly.
If False then the fname_templ will be used in glob to find the file.
dtime_placeholder : str
String used in fname_templ as placeholder for datetime.
Default value is "datetime".
"""
__metaclass__ = abc.ABCMeta
def __init__(self,
path,
ioclass,
mode='r',
fname_templ="",
datetime_format="",
subpath_templ=None,
ioclass_kws=None,
exact_templ=True,
dtime_placeholder="datetime"):
self.path = path
self.ioclass = ioclass
self.mode = mode
self.fname_templ = fname_templ
self.datetime_format = datetime_format
self.subpath_templ = subpath_templ
self.exact_templ = exact_templ
self.dtime_placeholder = dtime_placeholder
self.fid = None
if ioclass_kws is None:
self.ioclass_kws = {}
else:
self.ioclass_kws = ioclass_kws
def __enter__(self):
"""
Context manager initialization.
Returns
-------
self : GriddedBaseTs object
self
"""
return self
def __exit__(self, exc_type, exc_value, traceback):
"""
Exit the runtime context related to this object. The file will be
closed. The parameters describe the exception that caused the
context to be exited.
exc_type :
exc_value :
traceback :
"""
self.close()
[docs]
def flush(self):
"""
Flush data.
"""
if self.fid is not None:
self.fid.flush()
[docs]
def close(self):
"""
Close file.
"""
if self.fid is not None:
self.fid.close()
self.fid = None
def _open(self, filepath):
"""
Open file.
Parameters
----------
filepath : str
Path to file.
Returns
-------
success : boolean
Flag if opening the file was successful.
"""
success = True
self.close()
try:
self.fid = self.ioclass(filepath,
mode=self.mode,
**self.ioclass_kws)
except IOError as e:
self.fid = None
success = False
warnings.warn(f"I/O error({e.errno}): {e.strerror}", RuntimeWarning)
return success
def _search_files(self,
timestamp,
custom_templ=None,
str_param=None,
custom_datetime_format=None):
"""
searches for filenames with the given timestamp. This function is
used by _build_filename which then checks if a unique filename was
found.
Parameters
----------
timestamp: datetime
Datetime for given filename
custom_tmpl : string, optional
If given the custom_templ is used instead of the fname_templ. This
is convenient for some datasets where not all file names follow
the same convention and where the read_image function can choose
between templates based on some condition.
custom_datetime_format: string, optional
If given the custom_datetime_format will be used instead of the
datetime_format. This adds support to search for multiple files
for example for a given day, a given month or a specific year.
str_param : dict, optional
If given then this dict will be applied to the fname_templ using
the fname_templ.format(**str_param) notation before the resulting
string is put into datetime.strftime.
- example from python documentation:
coord = {'latitude': '37.24N', 'longitude': '-115.81W'}
'Coordinates: {latitude}, {longitude}'.format(**coord)
'Coordinates: 37.24N, -115.81W'
"""
if custom_templ is not None:
fname_templ = custom_templ
else:
fname_templ = self.fname_templ
if custom_datetime_format is not None:
dFormat = {self.dtime_placeholder: custom_datetime_format}
else:
dFormat = {self.dtime_placeholder: self.datetime_format}
fname_templ = fname_templ.format(**dFormat)
if str_param is not None:
fname_templ = fname_templ.format(**str_param)
sub_path = ''
if self.subpath_templ is not None:
for s in self.subpath_templ:
sub_path = os.path.join(sub_path, timestamp.strftime(s))
search_file = os.path.join(self.path, sub_path,
timestamp.strftime(fname_templ))
if self.exact_templ:
return [search_file]
else:
filename = glob.glob(search_file)
if not filename:
filename = []
return filename
def _build_filename(self, timestamp, custom_templ=None, str_param=None):
"""
This function uses _search_files to find the correct
filename and checks if the search was unambiguous
Parameters
----------
timestamp: datetime
datetime for given filename
custom_tmpl : string, optional
If given the fname_templ is not used but the custom_templ. This
is convenient for some datasets where not all file names follow
the same convention and where the read_image function can choose
between templates based on some condition.
str_param : dict, optional
If given then this dict will be applied to the fname_templ using
the fname_templ.format(**str_param) notation before the resulting
string is put into datetime.strftime.
example from python documentation
>>> coord = {'latitude': '37.24N', 'longitude': '-115.81W'}
>>> 'Coordinates: {latitude}, {longitude}'.format(**coord)
'Coordinates: 37.24N, -115.81W'
"""
filename = self._search_files(timestamp,
custom_templ=custom_templ,
str_param=str_param)
if len(filename) == 0:
raise IOError(f"No file found for {timestamp.ctime()}")
if len(filename) > 1:
raise IOError(f"File search is ambiguous {filename}")
return filename[0]
def _assemble_img(self, timestamp, mask=False, **kwargs):
"""
Function between read_img and _build_filename that can
be used to read a different file for each parameter in a image
dataset. In the standard implementation it is assumed
that all necessary information of an image is stored in the
one file whose filename is built by the _build_filname function.
Parameters
----------
timestamp : datetime
timestamp of the image to assemble
mask : optional, boolean
Switch to read already masked data which requires the
implementation of an read_mask_data() in the ioclass
Returns
-------
img: object
pygeobase.object_base.Image object
"""
filepath = self._build_filename(timestamp)
img = None
if self._open(filepath):
kwargs['timestamp'] = timestamp
if mask is False:
img = self.fid.read(**kwargs)
else:
img = self.fid.read_masked_data(**kwargs)
return img
[docs]
def read(self, timestamp, **kwargs):
"""
Return an image for a specific timestamp.
Parameters
----------
timestamp : datetime.datetime
Time stamp.
Returns
-------
image : object
pygeobase.object_base.Image object
"""
return self._assemble_img(timestamp, **kwargs)
[docs]
def write(self, timestamp, data, **kwargs):
"""
Write image data for a given timestamp.
Parameters
----------
timestamp : datetime.datetime
exact timestamp of the image
data : object
pygeobase.object_base.Image object
"""
if self.mode in ['r']:
raise IOError("File is not open in write/append mode")
filename = self._build_filename(timestamp)
self.fid.write(filename, data, **kwargs)
[docs]
def get_tstamp_from_filename(self, filename):
"""
Return the timestamp contained in a given file name in accordance to
the defined fname_templ.
Parameters
----------
filename : string
File name.
Returns
-------
tstamp : datetime.dateime
Time stamp according to fname_templ as datetime object.
"""
StartPos = self.fname_templ.find(self.dtime_placeholder) - 1
EndPos = StartPos + len(datetime.now().strftime(self.datetime_format))
StringDate = filename[StartPos:EndPos]
return datetime.strptime(StringDate, self.datetime_format)
[docs]
def tstamps_for_daterange(self, start_date, end_date):
"""
Return all valid timestamps in a given date range.
This method must be implemented if iteration over
images should be possible.
Parameters
----------
start_date : datetime.date or datetime.datetime
start date
end_date : datetime.date or datetime.datetime
end date
Returns
-------
dates : list
list of datetimes
"""
raise NotImplementedError(
"Please implement to enable iteration over date ranges.")
[docs]
def iter_images(self, start_date, end_date, **kwargs):
"""
Yield all images for a given date range.
Parameters
----------
start_date : datetime.date or datetime.datetime
start date
end_date : datetime.date or datetime.datetime
end date
Returns
-------
image : object
pygeobase.object_base.Image object
"""
timestamps = self.tstamps_for_daterange(start_date, end_date)
if timestamps:
for timestamp in timestamps:
yield_img = self.read(timestamp, **kwargs)
yield yield_img
else:
raise IOError("no files found for given date range")
[docs]
def daily_images(self, day, **kwargs):
"""
Yield all images for a day.
Parameters
----------
day : datetime.date
Returns
-------
img : object
pygeobase.object_base.Image object
"""
start = datetime(day.year, day.month, day.day)
end = datetime(day.year, day.month, day.day, 23, 59, 59, 999999)
for img in self.iter_images(start, end, **kwargs):
yield img
[docs]
def resample_image(self, *args, **kwargs):
return self.fid.resample_data(*args, **kwargs)
[docs]
class IntervalReadingMixin:
"""
Class overwrites functions to enable reading of
multiple images in a time interval as one chunk.
E.g. reading 3 minute files in 50 minute half-orbit chunks.
"""
def __init__(self, *args, **kwargs):
if 'chunk_minutes' in kwargs:
self.chunk_minutes = kwargs.pop('chunk_minutes')
else:
self.chunk_minutes = 50
super(IntervalReadingMixin, self).__init__(*args, **kwargs)
[docs]
def tstamps_for_daterange(self, startdate, enddate):
"""
Here we split the period between startdate and enddate into
intervals of size self.chunk_minutes.
These interval reference dates are then translated to
the actual file dates during reading of the chunks.
Returns
-------
intervals: list of tuples
list of (start, end) of intervals
"""
intervals = split_daterange_in_intervals(startdate, enddate,
self.chunk_minutes)
return intervals
[docs]
def read(self, interval, **kwargs):
"""
Return an image for a specific interval.
Parameters
----------
interval : tuple
(start, end)
Returns
-------
image : object
pygeobase.object_base.Image object
"""
start, end = interval
timestamps = super(IntervalReadingMixin,
self).tstamps_for_daterange(start, end)
if len(timestamps) == 0:
return None
dataset = {}
metadataset = {}
lons = []
lats = []
for timestamp in timestamps:
img = super(IntervalReadingMixin, self).read(timestamp)
for key in img.data:
if key not in dataset:
dataset[key] = []
dataset[key].append(img.data[key])
metadataset[timestamp] = img.metadata
lons.append(img.lon)
lats.append(img.lat)
for key in dataset:
dataset[key] = np.concatenate(dataset[key])
lons = np.concatenate(lons)
lats = np.concatenate(lats)
return Image(lons,
lats,
dataset,
metadataset,
interval[0],
timekey=img.timekey)