Question

我正在使用一个新的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的新手，但据我所知，它试图打电话给＆＃39; unit＆＃39;来自metpy.units包，不是来自numpy？任何人都可以建议吗？这不是我的代码，但是当我调用h_convergence函数时，我不知道为什么它不起作用。

Answer 1

我刚试过同样的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)

这似乎对我有用。

Answer 2

您正在尝试获取units f的{{1}}属性，我假设f.units的类型为f。 numpy.ndarray没有ndarray属性，因此您获得units。

您需要检查AttributeError为_gradient提供的内容，f，u或v来自h_convergence。它们是numpy数组还是其他类型？

AttributeError：＆＃39; numpy.ndarray＆＃39;对象没有属性＆＃39;单位＆＃39;

2 个答案: