0

我正在尝试使用 Geoviews Dataset 类绘制一个 xarray 网格。数据的形状为：(12,1300,1936)，坐标为（月、经度、纬度）。

实例化花费了太多时间（在 i5 32GB Ram 计算机中几乎需要 7 小时）。绘制一个较小的数据集是可行的（需要几秒钟）。

这是 xarray（名为 xmam）的字符串表示形式：

<xarray.DataArray (MeanTemp Month:: 12, Latitude: 1300, Longitude: 1936)>
array([[[ nan,  nan, ...,  nan,  nan],
        [ 14.,  14., ...,  nan,  nan],
        ..., 
        [ nan,  nan, ...,  nan,  nan],
        [ nan,  nan, ...,  nan,  nan]],

       [[ nan,  nan, ...,  nan,  nan],
        [ 16.,  16., ...,  nan,  nan],
        ..., 
        [ nan,  nan, ...,  nan,  nan],
        [ nan,  nan, ...,  nan,  nan]],

       ..., 
       [[ nan,  nan, ...,  nan,  nan],
        [ 17.,  17., ...,  nan,  nan],
        ..., 
        [ nan,  nan, ...,  nan,  nan],
        [ nan,  nan, ...,  nan,  nan]],

       [[ nan,  nan, ...,  nan,  nan],
        [ 14.,  14., ...,  nan,  nan],
        ..., 
        [ nan,  nan, ...,  nan,  nan],
        [ nan,  nan, ...,  nan,  nan]]], dtype=float32)
Coordinates:
  * MeanTemp Month:  (MeanTemp Month:) |S9 'January' 'February' ... 'December'
  * Latitude         (Latitude) float64 25.57 25.56 25.55 ... 14.76 14.75 14.74
  * Longitude        (Longitude) float64 -103.6 -103.6 -103.6 ... -87.49 -87.48

这就是我实例化 Geoviews 数据集的方式

gvds = gv.Dataset(xmam,kdims=['Latitude', 'Longitude'],vdims=['MeanTemp Month:'],dynamic=True)

我尝试使用和不使用dynamic给出类似结果的参数。

你认为这是什么问题？

请注意，Datashader它还没有任何作用，因为我只是在创建对象，甚至没有绘制它！

4

1 回答 1

0

问题有两个：

  1. xarray (xa) 是如何构建的。也就是说，dims参数引用坐标的名称。参数（如果coords是字典）需要与dim.

考虑以下示例：假设data是一个形状为 (12,100,100) 的 numpy narray

coords_months = range(12)
coords_lon = np.linspace(-103,87,100)
coords_lat = np.linspace(14,25,100)
dims = ['months','longitude', 'latitude']
coords = {'months':coords_months,'latitude':coords_lat,'longitude':coords_lon }
xdata = xa.DataArray(data,coords=coords,dims=dims,name='MeanTemperature')
  1. 正如philippjfr提到的，Geoviews (gv) DataSet 构造函数的声明是错误的。

相反，我使用了这个：

gvds = gv.Dataset(xdata)

绘图花费了一些时间（2分钟）。但是我可以使用定义的 regrid 方法：

from holoviews.operation.datashader import regrid
image = gvds.to(gv.Image,['Longitude','Latitude'],dynamic=True)
regrid(image)

效果很好。

其他相关问题在这里

于 2019-05-14T16:01:57.093 回答