一、事有蹊跷
$ v/ N9 V6 v# T4 { 接篇一,前面提到在使用cube.AI生成的c语言神经网络模型API调用时,输入数据数量是24,输出数据数量是4,但上文设想采集了三轴加速度传感器的x/y/z三个各数据,按Jogging(慢跑),Walking(走了)两种态势采集了两组数据.csv,那么在实际中应该是输入数据数量是3(x/y/x-value),输出数据数量是2(Jogging,Walking两种类别)。0 n4 o# P4 M6 x- `( }0 ?
, G1 t/ c# w7 j0 n; x
由于模型API是cube.AI基于训练输出的神经网络模型而生成的,想必应该是训练神经网络模型时参数设置问题,因此我们重新回顾HAR训练项目.) | E1 p; B: v L. t l* o1 t
" ^+ f9 A9 i- \
进入“STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR”目录,运行 python .\RunMe.py -h命令,查看参数设置指令帮助,可以看到--seqLength和--stepSize参数设置都和input有关,默认数值是24,就可以笃定在API中调用时,输入数据数量是24就来自于此。
$ i- V! {6 S6 O9 r( R q$ i- PS D:\tools\arm_tool\STM32CubeIDE\STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR> python3 .\RunMe.py -h
# \9 z1 E& O6 Z) x' a* y+ K. [ - Using TensorFlow backend.
7 v7 B6 m5 ?: Q& m& u - usage: RunMe.py [-h] [--model MODEL] [--dataset DATASET] [--dataDir DATADIR]
6 y1 t( j, W% h ? - [--seqLength SEQLENGTH] [--stepSize STEPSIZE] [-m MERGE]3 }6 X2 w7 W0 x! J$ H$ f- S
- [--preprocessing PREPROCESSING] [--trainSplit TRAINSPLIT]( X" z, |3 |8 U
- [--validSplit VALIDSPLIT] [--epochs N] [--lr LR], g8 i# b% ^5 C' g
- [--decay DECAY] [--batchSize N] [--verbose N]0 M/ N+ I5 _$ c2 {. r1 R5 S9 v
- [--nrSamplesPostValid NRSAMPLESPOSTVALID]
! z, m0 x5 I8 M& ?% J, ` -
E" v, N: @4 t4 |# C2 g7 w - Human Activity Recognition (HAR) in Keras with Tensorflow as backend on WISDM @5 C& U1 \ j; _8 X
- and WISDM + self logged datasets
% \/ k, G2 B0 B8 { -
, k* \5 v" D9 _, B' V - optional arguments:2 q+ W! Z- |! X6 G( \( J. b: u; E% ?
- -h, --help show this help message and exit: d9 C, _0 s, U) O3 q5 D% f6 P
- --model MODEL choose one of the two availavle choices, IGN or GMP, (
9 r3 n6 l" [2 f3 | - default = IGN )/ p8 k0 j" J( w2 H ~
- --dataset DATASET choose a dataset to use out of two choices, WISDM or
* O7 r. k# |/ c m) N - AST, ( default = WISDM )' i B# s8 [! i# q% [
- --dataDir DATADIR path to new data collected using STM32 IoT board; }5 _ b* I( _6 a& n3 J
- recorded at 26Hz as sampling rate, (default = )! D4 Y$ D( E0 H( U# W
- --seqLength SEQLENGTH
5 z' l. S1 W/ `3 A" a3 `1 i: f - input sequence lenght (default:24)
1 k6 @; V1 r; N0 @& Q7 A. q - --stepSize STEPSIZE step size while creating segments (default:24, equal0 z9 Q m* [1 v5 T* f
- to seqLen)( l( P( G& \5 N: W4 i2 k
- -m MERGE, --merge MERGE% Q& ^0 D9 \9 l9 b, D
- if to merge activities (default: True)
' X2 y7 E8 k. R" {, ?1 u - --preprocessing PREPROCESSING9 C" E# i/ M" u0 E# c( V
- gravity rotation filter application (default = True)
' B8 \& v" ~$ ?1 f0 A% t - --trainSplit TRAINSPLIT H4 b% p7 G+ o! s3 b3 K7 H
- train and test split (default = 0.6 (60 precent for5 L9 d( x% b$ U! Y+ L
- train and 40 precent for test))2 ~; {9 a$ K; B' L# [: \
- --validSplit VALIDSPLIT
. t$ r9 P* y: H; Y - train and validation data split (default = 0.7 (70/ C, M! _/ j9 s- r. {1 C! ?0 k
- percent for train and 30 precent for validation))# @8 A$ H. s1 [' k, K- P3 f
- --epochs N number of total epochs to run (default: 20)$ N& ?. S$ U& t1 `
- --lr LR initial learning rate n5 r, P7 x# M, _' [
- --decay DECAY decay in learning rate, (default = 1e-6)
+ B3 [/ z6 F- e9 T& A! q% X+ q0 C - --batchSize N mini-batch size (default: 64)
5 t* W6 [, l- ^6 y T - --verbose N verbosity of training and test functions in keras, 0,1 X9 n. y' m1 I' m# {
- 1, or 2. Verbosity mode. 0 = silent, 1 = progress bar,
3 j9 U' _: w+ }5 {9 ~& | ?0 @2 X5 T - 2 = one line per epoch (default: 1), T/ m4 @0 f3 c, v
- --nrSamplesPostValid NRSAMPLESPOSTVALID: L. e6 L4 t1 W& U# K& g
- Number of samples to save from every class for post
% g" ?, R# s, p, p+ @: q f - training and CubeAI conversion validation. (default =8 v. Y, ~* n, [7 y
- 2)
/ O, g+ K2 b9 T" ^ - PS D:\tools\arm_tool\STM32CubeIDE\STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR>
复制代码
1 ], Z/ J1 s& z$ [2 J4 r+ \4 v 而输出数据数量是4的原因追踪源码可以看到,来自于PrepareDataset.py(数据集预处理源文件),由于在参数--dataset默认设置是WISDM,因此分类输出即为'Jogging', 'Stationary', 'Stairs', 'Walking',即输出数据数量为4:
: o! q. Z7 a" m
. B1 z0 F7 G# z# Y- _# G
: X( U9 O- L# c; i$ R/ n& N* S
, L/ P* H$ w! D3 v% k/ H8 ] 下面来深入了解HAR(Human Activity Recognition,人类行为识别)案例为何会有这样的设置。! \3 W0 ~8 @; Y, |4 e+ f1 @
6 A! B) c1 q- f) k9 N/ d( t
二、HAR训练项目分析* P# Z' j/ Q2 B
由于前文仅仅采集Jogging(慢跑),Walking(走了)两种态势数据,但通常Stationary(站立不动)态势更常见和默认姿态,因此本文按篇一方法再采集一组Stationary姿态日志数据。
: X5 S3 i. I! z& ` U' v* _8 }, N0 c% t+ K7 \ [" t6 q
4 Y" r% D: k- \
$ ^$ }2 ]) x6 O5 X$ _: {5 o; C& b
并同样拷贝到HAR/Log_data目录下,数据如下:* ]+ Z' K ^5 [) r2 m
7 X: W3 e' \- B: c
. @4 b. F6 |7 K/ X5 ^8 i' \" R
- s: c2 W+ X8 y4 H- q 在PrepareDataset.py文件中,read_dataset( self )和preprocess_data( self, data )是用来处理WISDM数据集的,而get_data_from_file( self, fileName, preparedDataFileName )和prepare_self_logged_data( self )函数是用来处理自行采集数据集的。
4 {' l- `. z+ T5 z. E* X
0 E7 X* J F$ s/ l8 D( k
; [) t* E8 L; j: m: a
) E( z( c8 t) X+ l1 c, ^; I
read_dataset会读取WISDM数据集的datasets/WISDM_ar_v1.1_raw.txt文件$ T1 M9 M4 G/ }5 G
# p6 S0 C5 L. ]2 M& k
9 ]+ O% p3 }- k* F- h, N
2 P( N* o, R- N
这些数据是已经做了转换预处理的数据集,见下图红线标注部分。2 K8 O/ E3 g5 D
. v' X4 H O z% r+ V% p3 t, `
# L9 V" p( g1 U& V" g: ?0 ? Y( a" k1 } I* n! H
而get_data_from_file( self, fileName, preparedDataFileName )函数读取自行采集数据集(.csv)文件时,做了转换预处理,确保和WISDM数据集一致。
5 {( }3 v4 b7 c( m3 i( e% N: K9 Z. q# @( A' q. `
- p, c* L, A, N' a9 r
1 `2 L h+ t" a: L 同时项目还把预处理过的自行采集数据集与训练好的模型一并以csv格式输出到目录中,以源数据目录命名,例如Log_data.csv。2 ?' u! W! ?; ^, e
# H0 Z$ r: N/ Q% [# {# o$ K
( n. b2 A) y8 M" [: ~
5 G" X5 E4 S7 K3 s+ }$ k8 M 下来就是输入数据数量问题,由于在FP-AI-SENSING1案例项目,数据采集是按一定时间间隔连续实时采集的,采集间隔也很短(毫秒级别)9 u1 E! K. y# c3 [" G0 Z9 q
, s; G0 a# N* d0 A$ @$ s! B& i- E
! \9 T [( a: c& X& T( G' f) u
. ~ l4 x# a4 n9 {1 T
因此get_segment_indices和get_data_segments函数就是将采集到的连续数据处理成一个长度为seqLength的窗口的输入数据集,即每次输入数据是一段24组(一组3个数据值,x/y/z)数据,也就对应了我们前文定义输入数据缓存是static ai_float in_data[AI_HAR_IGN_IN_1_SIZE=24*3*1];,共72个float 数据。显然HAR项目做法是通过一组连续数据集作为输入比单个态势数据更能反应人类行为姿态的持续性,更贴近实际。, _$ k6 ~ X: m m2 y+ M/ |6 _
. d9 o+ _8 A3 {6 s0 a: E/ [
$ }& a: X# O; v0 X
: ]1 T1 M* r! K N+ ~6 q 在read_dataset函数加入打印读取的WISDM数据集语句。" B4 e3 R) L' i1 C; V5 G
: h# i. @- ^% F# B8 `5 ~6 s
2 b. k/ `5 v8 H4 |: e `8 e1 i' x1 ]; ?9 S& v0 {' C
在prepare_self_logged_data函数加入打印读取自行采集数据集语句8 Z N. O$ M# ~3 ^/ `$ x5 n7 V
( i2 w8 J2 V1 ^$ d, G9 ~" I) M2 S- P* ]
% m8 ^/ w# {% c( S
+ {" m, n0 O& S! @( _+ C 因为打印输出依据后,后续会打印数据按输入数据长度来分段的信息,暂时保持seqLength和stepSize为默认的24数值,运行python3 .\RunMe.py --dataDir=Log_data命令:- x! m+ [! R H% o; _8 S
- PS D:\tools\arm_tool\STM32CubeIDE\STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR> python3 .\RunMe.py --dataDir=Log_data/ |- K5 t: p4 A, Z% Y, {
- Using TensorFlow backend.2 y+ ^& |8 ~; ] a% O1 _: B" B
- Running HAR on WISDM dataset, with following variables }5 `' L! }5 N
- merge = True
* C2 r2 k) M" r X3 k$ @; x1 C4 t$ L - modelName = IGN,$ S3 o* @- M: B5 H: W
- segmentLength = 242 v$ n" P- j y5 q
- stepSize = 24( {. J* E& v& t& ~
- preprocessing = True' ]; o1 {2 I9 e0 Q2 H
- trainTestSplit = 0.6/ A/ @8 H2 b& S+ W
- trainValidationSplit = 0.7) |! X' \, S& P' v* w* S2 e( g1 {5 X
- nEpochs = 20
# e* o+ E1 I- N; h5 H - learningRate = 0.0005# u& f% S7 W+ y6 c; t w
- decay =1e-06
0 s- D: e c- G: f" L5 u1 l - batchSize = 640 x( u& S1 H7 `; Y. g5 J2 ?
- verbosity = 1$ S1 `' @: @# F7 y5 ?3 f
- dataDir = Log_data4 F% a: v. M$ f& o
- nrSamplesPostValid = 22 i7 _; U @7 {$ j$ R
- User Activity_Label Arrival_Time x y z
m# d- M/ v V F; j - 0 33 Jogging 49105962326000 -0.694638 12.680544 0.50395286;
- B- Y; T# D6 Y. {9 z - 1 33 Jogging 49106062271000 5.012288 11.264028 0.95342433;
; {( S G l# T' j, U9 N - 2 33 Jogging 49106112167000 4.903325 10.882658 -0.08172209;
+ m' i. k, e/ z) `- v2 t - 3 33 Jogging 49106222305000 -0.612916 18.496431 3.0237172;
3 n/ X/ K. h8 t# r - 4 33 Jogging 49106332290000 -1.184970 12.108489 7.205164;
* R. Y; P$ `: F5 P% d% q - ... ... ... ... ... ... ...
7 w% H* K% ~+ h P/ o - 1098199 19 Sitting 131623331483000 9.000000 -1.570000 1.69;
/ N% Y) S' k9 b' L1 j/ M: J - 1098200 19 Sitting 131623371431000 9.040000 -1.460000 1.73;+ n0 y4 N4 ~- W1 p& X! P: S$ V# w& \
- 1098201 19 Sitting 131623411592000 9.080000 -1.380000 1.69;
! m; b3 |5 r7 h, b1 P6 W - 1098202 19 Sitting 131623491487000 9.000000 -1.460000 1.73;9 z/ w1 y" c. Q% ?2 d9 O4 D
- 1098203 19 Sitting 131623531465000 8.880000 -1.330000 1.61;
* Q& i3 r+ ?5 B% ` -
9 k3 H+ n2 j3 @* N% ] - [1098204 rows x 6 columns]
Q+ {, c4 ]" J. C5 [ - Segmenting Train data
' Q9 V2 x3 O8 N& a - Segments built : 100%|███████████████████████████████████████████████████| 27456/27456 [00:28<00:00, 954.67 segments/s]
+ }1 b+ G) ]8 `' d h& O# u( |$ H - Segmenting Test data" U0 E8 r# ]5 Z h/ f, w
- Segments built : 100%|██████████████████████████████████████████████████| 18304/18304 [00:14<00:00, 1298.22 segments/s]9 K4 q1 H; u6 L& Z. l6 V! H ?( T
- Segmentation finished!
1 D9 g. B) b% X ? - preparing data file from all the files in directory Log_data
$ |. c% R# w) {. ~ - parsing data from IoT01-MemsAnn_11_Jan_23_16h_57m_17s.csv" l* L" S8 e& l5 C. E e4 G
- parsing data from IoT01-MemsAnn_11_Jan_23_16h_57m_53s.csv
4 o6 K/ x( O$ W- B5 `7 H8 ^3 Q: c - parsing data from IoT01-MemsAnn_26_Jan_23_15h_51m_01s.csv4 G# b4 {& ?* y$ ]4 J8 J# O
- x y z Activity_Label
: p* h6 j% x/ a4 c6 n - 0 -1.965414 -0.143890 9.367359 Walking0 Q9 Y1 R( V$ p. E! C5 u% Q
- 1 -1.629783 0.664754 9.618931 Walking
7 G$ r: ^: `, B) ^4 b: {" I' a - 2 -1.720833 0.384629 9.492079 Walking* _5 c: e+ x, g! b8 Y4 c
- 3 -1.681419 0.534637 9.648173 Walking/ M& e+ z# b6 k; X5 u
- 4 -1.729849 0.421650 9.557259 Walking7 q' w v1 [4 U" w6 f( H
- ... ... ... ... ...
3 g1 K# a& i6 C3 r1 F4 ?6 \8 B5 s - 2639 -1.171046 0.033572 9.746819 Stationary% x1 D) T! {( k5 e( W4 N8 O8 v
- 2640 -1.212873 0.007256 9.759376 Stationary
$ I5 U+ B. R2 H5 { - 2641 -1.212011 0.019485 9.753982 Stationary
% Z/ c* E4 T% L2 _: m5 P; j - 2642 -1.172311 -0.004511 9.734770 Stationary7 `0 m1 W( w% T+ F# R& b- @
- 2643 -1.175431 0.035787 9.753447 Stationary% @ v U. x. Y" H5 d
- / R2 T" @9 v7 y. G9 a# l0 x
- [2644 rows x 4 columns]7 Y. p2 i5 T% z$ H& l8 w
- Segmenting the AI logged Train data: N7 @8 i' l& N D* x
- Segments built : 100%|████████████████████████████████████████████████████████| 67/67 [00:00<00:00, 2795.31 segments/s]
! g; L }0 F9 h7 |3 s: j- ? - Segmenting the AI logged Test data+ `; ]5 }; t9 [. X) T
- Segments built : 100%|████████████████████████████████████████████████████████| 45/45 [00:00<00:00, 2370.94 segments/s]
$ m: |0 v1 X( `8 \% W6 G% ` - Segmentation finished!
: Q! @2 i/ F1 N* c5 ` - _________________________________________________________________/ u3 n7 f" N; i
- Layer (type) Output Shape Param #
. X7 [- Z6 u+ }/ ^! h; E - =================================================================9 H* [- T8 f8 j3 o& B0 M
- conv2d_1 (Conv2D) (None, 9, 3, 24) 4082 \/ L3 g% x: _( g8 m! `. d
- _________________________________________________________________
6 O R$ R% R( D; P7 ? - max_pooling2d_1 (MaxPooling2 (None, 3, 3, 24) 09 x/ p) q* Y% G3 S3 V4 M4 Y! L; M
- _________________________________________________________________
: V1 X7 e0 H5 p- X6 g - flatten_1 (Flatten) (None, 216) 0* G& H' V a1 a9 t# h4 I( B* V
- _________________________________________________________________, q5 B0 _4 C2 e7 P+ a
- dense_1 (Dense) (None, 12) 26048 Y" U) U5 Y6 p+ w2 p- ?! J9 k
- _________________________________________________________________
8 Y) h6 i0 _0 C5 }" l, o( K - dropout_1 (Dropout) (None, 12) 0, D$ c! b; \4 o/ m
- _________________________________________________________________
9 ~) Q% [9 L+ Q5 B( A z9 G& _# k - dense_2 (Dense) (None, 4) 527 \, ]1 @7 n! P% } ~
- =================================================================
/ h1 t' s7 M. R - Total params: 3,0640 {7 I5 r( L% u/ A8 i* o' i
- Trainable params: 3,064* ^( E+ d& H r
- Non-trainable params: 07 {; }7 ~) @- H$ B* f* U0 q5 h- x
- _________________________________________________________________4 z/ [5 I, `7 j
- Train on 19288 samples, validate on 8233 samples! @. D! Q- J) a t Y3 U# }
- Epoch 1/20& S2 ?& \! R4 E7 e+ ?1 s
- 2023-01-26 17:12:56.882726: I tensorflow/core/platform/cpu_feature_guard.cc:142] Your CPU supports instructions that this TensorFlow binary was not compiled to use: AVX2
& d8 M3 |( P. C/ V# d9 x - 19288/19288 [==============================] - 1s 54us/step - loss: 1.2022 - acc: 0.5290 - val_loss: 0.7089 - val_acc: 0.7409
Q2 x8 N% ^, [+ ^7 g( U* ^ - Epoch 2/20
' P# t% |" Q( B/ Q6 k# { - 19288/19288 [==============================] - 1s 41us/step - loss: 0.7520 - acc: 0.7017 - val_loss: 0.5342 - val_acc: 0.7985; F% E7 |/ ~1 k# Y0 c% K" b( N
- Epoch 3/20& W+ R, ~3 @, @6 p* Y$ e& W0 Y- v7 a
- 19288/19288 [==============================] - 1s 41us/step - loss: 0.6079 - acc: 0.7571 - val_loss: 0.4573 - val_acc: 0.8153
y. u% o6 \% a - Epoch 4/20) a- @3 c# t, k% C: g
- ... ...0 R2 Y0 M. K; r7 d n% r* j7 h
- 19288/19288 [==============================] - 1s 39us/step - loss: 0.3306 - acc: 0.8899 - val_loss: 0.2669 - val_acc: 0.9113* V, H* B6 J8 s$ B9 u9 z7 C
- Epoch 20/20
: W5 W( q% E! \& B5 _) O - 19288/19288 [==============================] - 1s 40us/step - loss: 0.3194 - acc: 0.8913 - val_loss: 0.2646 - val_acc: 0.9168: I9 E9 k4 C) T# j9 b. f" H
- 12831/12831 [==============================] - 0s 22us/step" k7 t& D# A9 p
- Accuracy for each class is given below.
1 e; b0 J1 x: v* @ \; W3 Y - Jogging : 97.51 %% F0 ^7 Q4 H w
- Stationary : 98.63 %4 ]& `! Y: }6 ` ~/ y) m" q
- Stairs : 70.94 %: L1 o4 ]' G) O/ c: U+ l' V
- Walking : 81.55 %" m+ U. J5 q1 [8 G- g- Y8 l0 K
- PS D:\tools\arm_tool\STM32CubeIDE\STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR>
复制代码
* y4 |% j3 D, i 读取到的WISDM数据集大小是1098204≈24*(27456+18304),而27456+18304是按--trainSplit参数默认值0.6:0.4比例切分的。读取到的自行采集数据大小是2644≈24*(67+45)。
6 z) l. l* J7 t3 ?- N- `2 N/ c0 I8 f9 y4 |, _$ b
9 m i( Y4 Z4 p/ F3 |4 n/ x: z/ E3 m9 D! d
三、工程调整
# g% G/ l4 j) l) I2 u2 t) T 从前面分析来看,输入数据数量还是采用默认数值24问题不大(各位粉丝可以自行调整其大小来测试及观察效果),输出数据数量由于有WISDM数据集参与,保持'Jogging', 'Stationary', 'Stairs', 'Walking'是中分类也OK。现注意到,其实输入数据是经过了预处理转换实际加速度值,而非是传感器原始输出数值。: z5 ^ i4 K8 f/ J
7 k8 m, M( V( R, s0 _
由于自行采集数据增加了Stationary数据文件,因此在cubeMX上采用新的训练模型文件重新生成c语言的神经网络模型。
- s# N0 _6 u, r, S: y, M! x
) f' `# f& H/ C4 Q; A
3 ]9 Z) I7 S+ r. `( I, ~6 Z2 k/ `9 K5 \8 Y; ]& ^7 ^
然后调整输入数据及预处理函数acquire_and_process_data,采用真实数据来测试0 m* ^9 Z) c* K, P: _
2 p6 ]4 Y' E ~ o6 ?. S
7 g; m5 h/ ?3 v3 z2 m* y
D4 d1 Y8 p2 U5 J. H
该函数有原来的随机赋值,& w4 Z* E- F. Q* O4 v/ o% w; u
- int acquire_and_process_data(void *in_data,int factor)
. ?! e. g& @. N; O) D7 J. V! _ - {% R7 l6 G( {0 `3 [2 X4 L
- printf("in_data:");( l! ]; |) Y: @- K1 B+ M
- for (int i=0; i<AI_HAR_IGN_IN_1_SIZE; i++)6 u- [1 q2 @( D2 A$ V
- {1 ]1 ?. L3 A4 A
- switch(i%3){: E7 V4 ^0 S! R, c {- H
- case 0:' f% ^4 o5 u8 ~7 Y0 X
- ((ai_float*)in_data)[i] = -175+(ai_float)(i*factor*1.2)/10.0;/ {% U+ T o" g1 r
- break;
8 \2 X) K: m/ }$ s) l - case 1:
" m) ^1 [* @1 D( [' z0 `' J6 s& T - ((ai_float*)in_data)[i] = 50+(ai_float)(i*factor*0.6)/100.0;
8 v: ~" `- A/ F) Q, b9 Y3 f1 W9 L1 a - break;
' U! q' G9 v* c( g7 R0 S" ~ - case 2:
' m G9 b0 g" _ - ((ai_float*)in_data)[i] = 975-(ai_float)(i*factor*1.8)/100.0;
9 M' m0 u4 R) U - break;9 t- h' j' w2 d/ t! n
- default: f/ ^6 I9 l# p' F8 ]; U7 L% o8 m P
- break;
" D/ p8 M$ F4 B( l# ~" H% a/ g - }
# t7 z* N: `* ]3 X1 o9 m5 k# q - printf("%.4f ",((ai_float*)in_data)[i]);
4 Y& J0 O4 ^, X# Q - }& h A. N0 ^' Q. S' l
- printf("\n");
, W# P/ J# n9 D6 O; o, A" d - 6 \; T" p+ Z1 m/ o% w, A
- return 0;
5 n; E/ M+ m& q3 l0 o4 C - }
复制代码 8 {* [' z9 ]( T9 }4 H8 m1 e6 z8 S
调整为- J' D) k7 B" A! T
- ai_float in_buf[] =2 L5 @: {1 Z+ j8 w
- {9 m4 Z" Y: K# R4 |% x' h7 Y
- -1.9654135467252831,-0.14388957575400915,9.36735860765576,
) M8 i: G1 m% c" T: \ - -1.6297827960727935,0.6647544204312931,9.618930851170278,0 X/ L7 G) w1 {- U- [& M! k
- -1.7208332169161387,0.38462856310059845,9.492079217583663,6 |+ u4 Q, @2 x+ _" c
- -1.6814190066807724,0.5346365723749872,9.648173176699613,% o3 G! ]% S7 g2 x* ^5 h
- -1.7298486510592452,0.42164964981080416,9.557259289818587,
4 B/ p( D/ I @ - -1.7618787694384546,0.45864558999786653,9.653153776935605,
+ z3 g- K! B3 l& k - -1.7410197123193858,0.4236369742675384,9.55486293595946,
# H- t5 V' U. c l5 f* Q* f, D - -1.7600076822930908,0.46214612481362705,9.594426626710453,
. G4 _1 y7 ~& V( J$ o& w2 x1 G - -1.5761631958773263,0.3715109598910308,9.436853714636964,
% h2 i9 [3 D l! O0 ]% z# p( C$ p - -1.5920827364351244,0.37070313540523914,9.66189484448469,* ?; b& E$ A1 G
- -1.6178308849438598,0.37500917334673567,9.695719226290015,
h3 _! x$ I# J1 r - -1.4388296833472143,0.6108605310285585,9.464814699883437,4 w9 Z' ]1 A! f' d
- -1.5651621282887258,0.5691273914891515,9.513897717476588,6 p% `" P/ K2 K' U5 I
- -1.4637992479412343,0.5105873209777632,9.501636895304161,
2 o* _9 X. ?$ k, A" ~ - -0.6794677157685166,0.5024637601753793,8.96404801376064,0 D" P8 l2 @7 s& @( L) a2 o
- 0.2600149177042748,0.6699546179356337,8.903349009412763,
- r$ ?: f4 h$ e$ j3 n- u, _ - 1.0712686735261918,1.4889662656074603,9.520348132500752,: Q2 j' K& L8 z; l2 O; |& m' R* Z
- 0.3914123345764725,1.4210706041563634,10.557387805652848,
5 v8 ~" v) ~# n, l6 W: x% C! M$ Y' L - 1.0779003359396493,1.0582703827741018,10.454469820960814,
4 o* _$ n. ?. z - 0.12433283758079197,-0.27273511643713033,10.328552286632643,
J; d9 T/ [% l - -0.010219096051988997,0.2961821896002729,9.483084545625971,* R9 Q% Z7 w0 {: \, s/ I" \; T
- -1.6910112286007235,-0.2898761724876157,9.704755735796937,
3 B# j( J% X- ^9 B - -2.693651827312974,-0.41126025575408387,9.825328217800239,
% Z/ ^4 |0 }) u" B( r0 h" T/ a9 ] - -2.8416981790648177,-0.14586229740441406,9.880552703938179
! }( t& i. `) q8 Q2 m3 }5 m$ P - };- ^7 D3 c0 c ?, x" i& a j
-
, `' k: o Z& F- `" w8 k$ x - int acquire_and_process_data(void *in_data,int factor)
7 d9 V" I- W# U# Z0 M - {
- }2 n! J* ?: S$ Z9 |- s - printf("in_data:");' Q1 r; d d1 C: f# i1 w$ ^# V
- for (int i=0; i<AI_HAR_IGN_IN_1_SIZE; i++)
, o9 l4 `3 v1 I9 S& c - {
# e/ e' W! ~# @ - ((ai_float*)in_data)[i] =in_buf[i];
2 x" A# o) b1 m2 H% b - printf("%.4f ",((ai_float*)in_data)[i]);
) L% a+ ?! @# @$ P ?5 ?5 r - }9 ]% q T" E: q5 V) g
- printf("\n");6 b* m, o+ K+ G# W
- 4 D1 ?; `3 Q1 H+ V# M! h O# i
- return 0;1 N& _) [4 \ ~( i
- }
复制代码
# X( H# h9 j) i6 Y3 Q 重新编辑及加载到开发板,打开串口助手,输入test*,查看测试效果如下,显然和期待有很大出入:
- o, K n- D( r4 G- B; @% V+ t& p, r( X, f+ p% j6 O2 K
+ B+ e, ^0 i- `, Z4 _
/ ~' z1 E0 M7 w# f, v$ N
四、测试有误问题
3 J; T# A: ]% q M! L6 J* B: { 上述结果和HAR项目中的PrepareDataset.py分类似乎并没匹配上,那么就要分析PrepareDataset.py分类和生成后的c语言模型时如何对应的,以及用来测试的数据是否符合模型要求。+ R% {4 t( L3 }3 n/ c
& U3 q& _1 B% x2 x: R
, S2 H3 b. ]1 F6 N
7 c! A, i* q- A* `4 Z2 ~; } 为此,再次回到HAR训练模型项目,在RunMe.py文件中追加以下语句:+ e0 `; B) J% b( i& M
- print("TestX:")
7 I! ~4 Y7 a3 }' C5 f* z" V - print(TestX[0:1])+ F3 X& z0 S9 M8 J
- print("TestY:"): R$ ^( W" ^3 {& p$ J1 `% x' l, z: D
- print(TestY[0:1])
复制代码
$ C7 e7 j' c' C4 P4 x9 s 以获取用来测试的真实模拟数据,再次运行python3 .\RunMe.py --dataDir=Log_data命令,最后部分输出如下:; @8 u9 ?6 p8 L% b& e' D
- TestX:; L# a2 ?; m% X! b5 M" E
- [[[[-3.64077855e-16]
. ?! q$ M9 j' n3 z - [-4.69612372e-16]' U8 C0 Y& p6 t; o9 c' Z% w2 p) {
- [ 3.68092310e-10]]
/ }' ]. p# I" g6 V% } -
& O+ E3 z& m1 A# m5 n* t - [[ 1.02756571e+01]
7 x5 k% x9 j! l - [-1.14305305e+01]
+ G8 i( q* w9 m; B - [ 2.61872125e+01]] j! x, A& J9 }( y' D0 u) j6 |/ G
-
8 j7 N$ S' @) N% i, j9 q( B2 P6 r0 W - [[-2.84181689e+00]& e2 v- t9 U' [& ?: ^
- [-3.54747048e+00]- \0 M7 s9 d, A, N4 t6 S# L: P4 j
- [-5.51206446e+00]]
/ B ?7 S5 N ~9 e -
& F W3 u& U: O& A Y. B; i - [[-3.82102513e+00]4 ]" s$ H1 R8 ~9 l
- [-1.41233186e+01]9 E" J$ ]( W" ]! t+ G
- [-4.59900586e+00]]
! v4 N& |# L* `4 m$ n1 K - - F8 O; _* p& T) J D8 [$ O
- [[ 6.68010824e+00]
2 {( p$ H3 U5 E4 N) R, f! S) r# p9 ^ - [ 9.39457601e+00]' \& h1 ?! I( r# ~- ? @+ t1 a( f
- [-2.96397789e+00]]6 [4 u% i4 D, m1 q% [& a9 X( f, T
-
! Y5 z/ G; c/ s. J! p# j9 A7 o% ? - [[-1.71771782e+01]: f P& Z8 T" w
- [ 1.19374149e+01]
/ |7 {) k ~7 g' d+ o* f7 x* f - [ 3.05770680e+00]]
. N1 h6 m6 d2 v! i# b* ]8 y - 9 {0 k3 K* g3 A0 ^/ V8 O
- [[-6.65782005e+00]0 h' V9 l( }1 b6 T6 R( w# _8 [
- [ 2.39062819e+00]/ F6 x" |2 G8 g
- [ 3.22844912e+00]]
' X+ ^/ R# O$ K -
) n w0 o0 D9 M& j h+ i# C0 P - [[ 4.59021292e+00]/ Z2 X0 t; V" u
- [-6.27548028e+00]9 c @7 c5 Z+ i8 D3 j
- [-4.92783556e+00]]
# W0 Q1 }- g2 s6 ?$ l/ N -
' h5 }8 n5 ^" R; C a d3 ]( O* S: ~ - [[ 8.03018658e+00]6 x; O! C! T+ ^- w0 Z4 @- ]
- [-2.72208600e+00]
% O2 J5 A' _2 m m* O/ Y8 ?3 t+ _1 Z - [-6.35796053e+00]]6 M' `" {% ]7 f/ d$ r& K% ^
- ( D1 Q! w+ ]/ E
- [[ 7.73164454e+00]+ O9 I% v" q' w# b1 C$ w
- [-6.31879160e+00]1 U$ Z" J& @$ J) N; A
- [-5.90723810e+00]]1 V3 S( R/ B- ]0 n, y3 p! m
-
7 e1 i7 k& O! V) s' t3 t" k/ _ - [[ 8.53803514e-01]
. {7 R2 m/ i6 \ - [-9.75763211e+00]7 j# j1 d3 S& n5 F# V. e' P
- [ 1.02466115e+01]]# M1 j; I- B9 y+ I2 X
- {, X) O. j X2 g9 ~
- [[ 1.11299171e+01]# F5 i0 e+ N. O( ?$ {& x+ b
- [-1.70658346e+01]9 o# J9 d' y6 [/ [
- [ 2.18511283e+01]]+ U3 k5 m8 j% D9 M+ x( E
-
. U2 C& |/ j. Y- i4 A5 Q9 A - [[ 3.92044994e-01]2 A5 J/ y$ {- y0 w
- [ 5.94768181e+00]7 f2 `7 X5 p+ M3 V' {
- [ 4.30131750e+00]]
/ {+ O( g a) D: G, x - ' Q" s$ b" _# Y {4 X# j
- [[-5.61807988e+00]
. k. D; ^4 |; X6 ? - [ 1.97310400e+01]
) \, o2 C7 W- C! _4 d. R - [-2.22512540e+00]]% S8 |3 N, I4 E) N+ x0 |
-
( l3 R6 F+ M/ Q# C - [[ 3.86836548e+00]0 C- \' B# ]% @; w3 F1 U2 L1 P
- [ 1.71617325e+00]
9 Q1 e7 u) M6 u, L% E0 e - [-5.86292387e+00]]
- t' W" e% T! q7 `, N e -
+ A- w0 g- P) i9 r' p, O2 I; D+ n - [[ 7.65913325e+00]
1 ^$ t X" g) j/ Y+ N' \1 _ - [-7.19628424e+00]# J6 W5 K3 ~+ I$ w
- [ 2.01628025e+00]]
/ e, ?: [: n% K: h5 C+ g - & J P2 c& P2 P" r$ }
- [[-7.52357836e+00]
. p# P& ^/ ]" R! A' s - [ 3.68102584e+00]
) ?, ^& I1 `; [' J! V - [-1.22753233e+01]]
1 ?7 L3 o% W: p4 j2 x+ i -
6 f1 r4 v* f$ T: ^( Q& K - [[-5.12351958e+00]
; v$ U! K7 T J- p( A - [ 1.23941669e+01]
# [$ a. r* k! m3 w) A, Q4 I - [-1.77385540e+00]]
9 X, h* n- e" E: @ - 7 x8 `: W; q/ N6 ^: L
- [[-4.86155823e-01]
$ @" v. ~* r& l( _8 ?" _ - [ 1.26333902e+01]
; O6 j% n7 D: `. o0 k; y. p - [ 5.93595914e+00]]6 M$ _* t) C7 E2 v |
- ( ^0 y, \* I6 t0 b/ A
- [[-1.96569165e+01]
5 Z( i6 t% x9 ?5 b' z# R& N5 M - [ 1.00467317e+01]4 T2 V( S2 T9 }# J1 Y6 N
- [ 9.47374003e+00]]* v, B4 J4 ~4 M$ P8 Q
- ; b) b( K0 R1 g
- [[-4.34050581e+00]
s3 w9 @) g0 T# K/ M, I$ ^ - [ 5.16311148e-01]) n# o3 j( n7 l
- [-5.63004156e-01]]
3 x+ S3 s+ y4 R+ F% { p - ; Z9 u w+ e5 K1 v; a) q
- [[-3.57974669e+00]
) K) B4 G& W" [; t/ T" r, H" g - [ 4.87240857e-01]/ W2 Q. m8 t- X! y# e3 r
- [-9.38271247e-01]]
( M6 L G k q7 ~ - 9 x! f! X; s6 Q6 F
- [[ 6.11930536e+00]& b U8 I1 d$ i% z' u. I" V( G7 e
- [ 5.99067573e+00], J/ g' U0 R6 a" x, }) T/ N9 `8 z
- [-7.68834262e+00]]
" ?0 B7 b. Q/ J0 I -
9 J) l' X$ S1 H5 T - [[ 1.12153409e+01]' w% O$ l5 A) F* g& q0 U" Y
- [ 2.37168199e+00]3 M* ^; N+ w% V2 k: d
- [-7.40963357e+00]]]]
X1 v8 E2 C/ @3 n2 Y$ v - TestY:
7 n, r1 K& X" E1 Y+ S - [[1. 0. 0. 0.]]
, [* L0 _+ X% x" L3 O - PS D:\tools\arm_tool\STM32CubeIDE\STM32CubeFunctionPack_SENSING1_V4.0.3\Utilities\AI_Ressources\Training Scripts\HAR>
复制代码
/ w& P, `; \* Q M+ g/ K3 Y 显然输入24*3的一组数据,输出结果是[1. 0. 0. 0.],再将该数据用于测试 c语言模型神经网络模型API调用。) |5 u& v& v" x% A" U
- ai_float in_buf[] =
9 G b- N6 H4 e, L: d - {
/ L6 f% q0 _; Z ~) Q2 K! B F - -3.64077855e-16,-4.69612372e-16,3.68092310e-10,5 H8 Z) `0 q6 ~% K4 z, }+ V$ `' q
- 1.02756571e+01,-1.14305305e+01,2.61872125e+01,
/ t" g3 ^2 Y4 G+ } - -2.84181689e+00,-3.54747048e+00,-5.51206446e+00,
) z& r: P$ l Z% P7 Z/ [ - -3.82102513e+00,-1.41233186e+01,-4.59900586e+00,: E( H$ q+ o5 n# a" N, q
- 6.68010824e+00, 9.39457601e+00,-2.96397789e+00,
" X- H+ a$ A; F% @, \ - -1.71771782e+01,1.19374149e+01,3.05770680e+00,3 g# Q, x5 g) g% ]5 c$ n4 v- h
- -6.65782005e+00,2.39062819e+00,3.22844912e+00,) U$ `& b' Q' `1 g" O
- 4.59021292e+00,-6.27548028e+00,-4.92783556e+00,+ O# t2 k1 K" g/ u C% x* g
- 8.03018658e+00,-2.72208600e+00,-6.35796053e+00,
3 p6 o8 J' m: H E: ]2 B; V9 l) @ - 7.73164454e+00,-6.31879160e+00,-5.90723810e+00,# A7 K9 \4 G7 ]
- 8.53803514e-01,-9.75763211e+00, 1.02466115e+01,
2 t/ r; N; C! K; {* y - 1.11299171e+01,-1.70658346e+01, 2.18511283e+01,
" V& b$ B2 f6 P! |( H8 q - 3.92044994e-01, 5.94768181e+00, 4.30131750e+00,
0 E/ F& L& L, g. E" a6 s9 g - -5.61807988e+00, 1.97310400e+01,-2.22512540e+00,
5 F3 l' _6 \7 X - 3.86836548e+00, 1.71617325e+00,-5.86292387e+00,0 O$ a- N4 ^+ E3 ?
- 7.65913325e+00,-7.19628424e+00, 2.01628025e+00,
* F: ^: S" b) T - -7.52357836e+00, 3.68102584e+00,-1.22753233e+01,
+ g' b9 o+ R1 u* B g! W - -5.12351958e+00, 1.23941669e+01,-1.77385540e+00,1 e& p+ G4 l$ }
- -4.86155823e-01, 1.26333902e+01, 5.93595914e+00,7 z4 U; ~% H/ E
- -1.96569165e+01, 1.00467317e+01, 9.47374003e+00,
M$ @* i$ ^/ j+ t4 x - -4.34050581e+00, 5.16311148e-01,-5.63004156e-01,0 _4 _ U& e$ F) U+ s, ^# t
- -3.57974669e+00, 4.87240857e-01,-9.38271247e-01,6 K6 U1 u& h1 U
- 6.11930536e+00, 5.99067573e+00,-7.68834262e+00,% d, j- Q1 D0 | V# r
- 1.12153409e+01, 2.37168199e+00,-7.40963357e+00
# E) o" ?$ Y. L) B8 B6 I - };
7 z8 A& w, K( I% A& s5 M { - ; z% [/ o% F$ O7 p# {5 N8 R
- int acquire_and_process_data(void *in_data,int factor)# r0 o k n: {) `7 e( S
- {8 w9 @0 ~0 Z- ~
- printf("in_data:");
6 O- J. d2 K' d - for (int i=0; i<AI_HAR_IGN_IN_1_SIZE; i++)
% i! {9 `* E; Z) q" o, E- L - {
0 X G4 Z$ |6 x1 d! W& W - ((ai_float*)in_data)[i] =in_buf[i];% |& ]" H$ b* n) G
- printf("%.4f ",((ai_float*)in_data)[i]);7 h& H$ p7 h# f+ W4 ]* |
- }( i! c: b1 C1 H/ ~, C- w" v$ O
- printf("\n");1 P; P# ?! R9 @% O" ~/ c4 H* u, S
- 6 J6 I' k) \: c& d
- return 0;% r$ C9 |1 I! @. @# ^5 t
- }
复制代码 Z: E. F. H7 y% R
再次编译下载程序,采用串口助手测试如下,输出数据为,【0.943687 0.000000 0.000294 0.056019】,而实际数据是[1. 0. 0. 0.],显然能对应上,只是精度问题:
- c* q5 y! C+ I& W
3 W+ B- E0 N* P8 H. u
- u# Z: {2 u1 h/ \8 Q0 f
0 |' z4 C+ m& l& Q# F+ ]7 x, L0 k4 L
五、精度问题7 G6 F) L! B: Z* b( ]) I# ]2 K! w3 W
再次回到HAR训练项目上,如前文所述,训练结果按时间生成一个目录进行输出:7 _1 _" p1 Z4 Z* _% |8 L+ G! m* z: m
% O7 v5 W7 }' |) b+ K
6 j6 U6 w0 g9 E) n6 s8 Y0 H, z
- b' z" E; }1 \: \+ }
输出了训练的acc 、loss变化曲线图和混淆矩阵图,观察可以看到,用于训练数据其主要还是来自于WISDM的数据集占据大头,从混淆矩阵也可以看出,为何前面测试时输出结果0.943687 0.000000 0.000294 0.056019】,显然训练过程中,'Jogging'姿态会有一部分判断 'Stairs', 'Walking'姿态。; F1 }6 q, B* }5 N+ b% l8 E+ P
+ l O: `% |3 m5 U
2 L# u- v, W* V& J" k
2 H& j2 c+ Y4 t+ A4 [% X+ @+ X
在训练过程日志显示:
$ M( D& }1 e" l- c) _; I; J5 g: W: Z7 R; w0 O8 b
8 s) v, o) v( E" q y
1 K. r; F3 j, Q0 C 再次改写RunMe.py,加入下面打印语句,查看用于训练、验证及测试的数据集量级:3 D" p; h% |. x" v- N( B: P
- print("TrainX:"+str(len(TrainX))); m& u1 k' v8 [ S$ J( ?- y
- print("TrainY:"+str(len(TrainY)))
1 H& b- b/ W' V$ z' H/ u9 p - print("ValidationX:"+str(len(ValidationX)))& x4 }, w5 @" u: O, g6 ?+ {: e
- print("ValidationY:"+str(len(ValidationY))): b/ z, N# g2 l7 c
- print("TestX:"+str(len(TestX))) J Z( J& ?/ }0 n0 g
- print("TestY:"+str(len(TestY)))
复制代码 . s1 ?6 A, }! g) K
再次运行python3 .\RunMe.py --dataDir=Log_data命令:- v: |; r/ V* p
! X% ~, [2 C% _0 f( [8 j8 a
. h) s6 w) j6 J4 `6 U) e4 L" w
& _' s% ~* f9 O+ ~
通常在实际项目中,不断优化算法模型是主要原因,但相信ST公司及其合作团队应该就算法模型做了足够好,而本文测试验证为何还出现精度相差这么大,应该是以下几点:- H; e" v) a0 A7 J( B& P B/ S/ k
1)本文用于训练模型的数据集不是同一设备,同一批次,统一采集方法获得足够量数据集。
% q; p/ D# N; `. r {. K! j
2 ]/ ]* _* B! N3 x0 Q0 T4 ? 2)本文自行采集的数量偏小,可能预WISDM采集数据方法存在差异而出现精确性问题,并加入其中,还是会影响模型训练精度的(若仅采用WISDM采集数据集)。# c1 v/ l& E7 }7 e" O
8 H+ }" A; G8 Q3 a
3)本文仅采用了HAR训练项目建议的默认参数,实际项目中,可能需要我们根据数据集及部署环境,不断尝试不同参数设置来训练和使用神经网络模型。
) X+ J0 v2 z; p/ X
( U9 B$ X! n& d& I# a) }1 y 4)cube.AI将keras训练的神经网络模型转换为c语言的神经网络模型虽然损耗极小,但还是存在一定精度损失的。
7 p/ s# R5 S9 x8 L
. f! R+ r* U& Y4 P1 X9 C
# I7 w6 ]: P ^$ v5 `- w. Y4 |) @2 q# |! u
篇一、二,是基于FP-AI-SENSING1案例和配套的B-L475E-IOT01A开发板来使用cube.AI软件包的,那么在实际使用cube.AI时,需要与我们项目实际硬件平台,并基于该硬件平台采集的数据进行数据预处理和神经网络模型构建、训练、验证及测试,然后再通过cubeMX和cube.AI将训练好的神经网络模型转换为c语言支持的AI模型,通过嵌入式程序加载(AI模型+cube.AI)从而实现神经网络计算的前端或边缘端部署,请阅读偏三。) X& ]. N% v/ E' F
————————————————
# j0 {7 f" l7 J8 @' N* @2 a版权声明:py_free-物联智能' a5 M; A7 X5 q& |
如有侵权请联系删除; G$ U* g$ _8 v. l3 q) M0 b4 P7 l
|