1

我正在使用一个新的 Python 包 (metpy)，其中包含几个子包 - 其中一个使用以下代码定义梯度和风收敛：

from __future__ import division
import numpy as np
from ..package_tools import Exporter
from ..constants import g
from ..units import atleast_2d, concatenate, units  
### the 'units' and 'constants' packages are within metpy

exporter = Exporter(globals())

def _gradient(f, *args, **kwargs):
    if len(args) < f.ndim:
        args = list(args)
        args.extend([units.Quantity(1., 'dimensionless')] * (f.ndim - len(args)))
    grad = np.gradient(f, *args, **kwargs)
    if f.ndim == 1:
        return units.Quantity(grad, f.units / args[0].units)
    return [units.Quantity(g, f.units / dx.units) for dx, g in zip(args, grad)]

def _stack(arrs):
    return concatenate([a[np.newaxis] for a in arrs], axis=0)


def _get_gradients(u, v, dx, dy):
# Helper function for getting convergence and vorticity from 2D arrays
    dudx, dudy = _gradient(u, dx, dy)
    dvdx, dvdy = _gradient(v, dx, dy)
    return dudx, dudy, dvdx, dvdy

def h_convergence(u, v, dx, dy):
    dudx, _, _, dvdy = _get_gradients(u, v, dx, dy)
return dudx + dvdy

当我尝试使用函数 h_convergence 时，我在这一行收到以下错误：

if f.ndim == 1:
    return units.Quantity(grad, f.units / args[0].units)
return [units.Quantity(g, f.units / dx.units) for dx, g in zip(args, grad)]

AttributeError: 'numpy.ndarray' object has no attribute 'units'

我是 Python 的新手，但据我所知，它试图从 metpy.units 包中调用“units”，而不是从 numpy？任何人都可以建议吗？这不是我的代码，但我不知道为什么当我调用 h_convergence 函数时它不起作用。

4

2 回答 2

3

我刚刚尝试了相同的metpy包并遇到了同样的问题：

AttributeError: 'numpy.ndarray' object has no attribute 'units'

如果您查看metpy 文档，第 1.2.1 节构造，它们解释了如何构造这些单元。

由于我也对收敛感兴趣，所以我做了以下事情：

from metpy.units import units

u1 = u*units('m/s')
v1 = v*units('m/s')

然后你应该为你的 dx 和 dy 做同样的事情：

dx = dlon * units('m')
dy = dlat * units('m')

最后的调用是：

conv = calc.kinematics.h_convergence(u1,v1,dx,dy)

这似乎对我有用。

于 2016-02-17T19:44:03.373 回答
0

您正在尝试获取( ) 的units属性，我假设它的类型为。s 没有属性，所以你得到.ff.unitsfnumpy.ndarrayndarrayunitsAttributeError

您需要检查您正在喂食的东西_gradient，f似乎是u或v来自h_convergence。它们是numpy数组还是其他类型？

于 2015-09-03T15:03:19.317 回答