数学度坐标到气象度坐标

Question

我有一组激光雷达数据，记录在数学度坐标（即90度是北），但我试图使这在气象坐标（即0，360度是北）绘图

我知道需要执行的数学，以便转换为工作，但实施是困难的

我不确定将数据作为list或np.array处理是最好的。我目前正在尝试以数组的形式使用它

if np.where((wdir1>0)&(wdir1<90)):
    wdir_1 = abs(wdir1-90)

elif np.where((wdir1>90)&(wdir1<360)):
    wdir_1 = abs (wdir1-450)'

但是，在执行此操作时，数组中的所有值都要减去90。而我只希望0和90之间的值减去90。大于90和小于360的值要减去450

Answer 1

在应用该操作之前，需要对数组进行切片。 这样，数组只有在条件为True时才会改变

wdir_1[np.logical_and(0 <= wdir_1, wdir_1 < 90)] -= 90
wdir_1[np.logical_and(90 <= wdir_1, wdir_1 < 360)] -= 450
wdir_1 = np.abs(wdir_1)