eo2cube / odcr Goto Github PK

band_ratio= data.VV- data.VH
data["polRatio"]= band_ratio
print(data)

<xarray.Dataset>
Dimensions:      (time: 59, y: 1538, x: 1718)
Coordinates:
  * time         (time) datetime64[ns] 2019-01-10T16:52:32 ... 2019-12-30T16:...
  * y            (y) float64 5.96e+06 5.96e+06 5.96e+06 ... 5.991e+06 5.991e+06
  * x            (x) float64 3.599e+05 3.599e+05 ... 3.942e+05 3.942e+05
    spatial_ref  int32 32633
Data variables:
    VH           (time, y, x) float32 -19.91 -18.74 -17.19 ... -19.88 -20.44
    VV           (time, y, x) float32 -12.08 -12.63 -10.78 ... -10.14 -11.54
    polRatio     (time, y, x) float32 7.832 6.111 6.404 ... 12.11 9.747 8.906
Attributes:
    crs:           EPSG:32633
    grid_mapping:  spatial_ref
    

polRatio <- ds$VV - ds$VH
polRatio$attrs$units= "VV - VH"

ds[["CR"]] <- polRatio

Error in py_set_attr_impl(x, name, value): AttributeError: cannot set attribute 'CR' on a 'Dataset' object. Use __setitem__ styleassignment (e.g., `ds['name'] = ...`) instead of assigning variables.

Detailed traceback:
  File "/home/datacube/anaconda3/envs/jupyterhub/lib/python3.9/site-packages/xarray/core/common.py", line 275, in __setattr__
    raise AttributeError(

Traceback:

1. `[[<-`(`*tmp*`, "CR", value = <environment>)
2. `[[<-.python.builtin.object`(`*tmp*`, "CR", value = <environment>)
3. py_set_attr(x, name, value)
4. py_set_attr_impl(x, name, value)

polRatio <- ds$VV/ds$VH
polRatio$attrs$units= "VV/VH"

load_sq= seq.Date(from = as.Date("01/10/2019", format= "%m/%d/%Y"), to= as.Date("12/31/2019", format= "%m/%d/%Y"), by= "6 day")
intp_sq= seq.POSIXt(from = as.POSIXct("01/10/2019", format= "%m/%d/%Y"), to= as.POSIXct("12/31/2019", format= "%m/%d/%Y"), by= "day")

intp_sq= intp_sq[-which(intp_sq %in% load_sq)]

CR_intp= polRatio$interp(time= intp_sq[1:length(intp_sq)])

Error in py_call_impl(callable, dots$args, dots$keywords): KeyError: numpy.datetime64('NaT')

Detailed traceback:
  File "/home/datacube/anaconda3/envs/jupyterhub/lib/python3.9/site-packages/xarray/core/dataarray.py", line 1725, in interp
    ds = self._to_temp_dataset().interp(
  File "/home/datacube/anaconda3/envs/jupyterhub/lib/python3.9/site-packages/xarray/core/dataset.py", line 3131, in interp
    obj, newidx = missing._localize(obj, {k: v})
  File "/home/datacube/anaconda3/envs/jupyterhub/lib/python3.9/site-packages/xarray/core/missing.py", line 559, in _localize
    imin = index.get_loc(minval, method="nearest")
  File "/home/datacube/anaconda3/envs/jupyterhub/lib/python3.9/site-packages/pandas/core/indexes/datetimes.py", line 705, in get_loc
    raise KeyError(orig_key) from err

Traceback:

1. polRatio$interp(time = intp_sq[1:length(intp_sq)])
2. py_call_impl(callable, dots$args, dots$keywords)