statisticalTest.lyx 10 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555
  1. #LyX 2.0 created this file. For more info see http://www.lyx.org/
  2. \lyxformat 413
  3. \begin_document
  4. \begin_header
  5. \textclass book
  6. \use_default_options false
  7. \maintain_unincluded_children false
  8. \language french
  9. \language_package default
  10. \inputencoding latin9
  11. \fontencoding T1
  12. \font_roman default
  13. \font_sans helvet
  14. \font_typewriter default
  15. \font_default_family sfdefault
  16. \use_non_tex_fonts false
  17. \font_sc false
  18. \font_osf false
  19. \font_sf_scale 100
  20. \font_tt_scale 100
  21. \graphics default
  22. \default_output_format default
  23. \output_sync 0
  24. \bibtex_command default
  25. \index_command default
  26. \paperfontsize 10
  27. \spacing single
  28. \use_hyperref true
  29. \pdf_bookmarks true
  30. \pdf_bookmarksnumbered false
  31. \pdf_bookmarksopen false
  32. \pdf_bookmarksopenlevel 1
  33. \pdf_breaklinks false
  34. \pdf_pdfborder false
  35. \pdf_colorlinks true
  36. \pdf_backref false
  37. \pdf_pdfusetitle true
  38. \pdf_quoted_options "linkcolor=blue"
  39. \papersize a4paper
  40. \use_geometry true
  41. \use_amsmath 1
  42. \use_esint 1
  43. \use_mhchem 0
  44. \use_mathdots 0
  45. \cite_engine basic
  46. \use_bibtopic false
  47. \use_indices false
  48. \paperorientation portrait
  49. \suppress_date true
  50. \use_refstyle 0
  51. \index Index
  52. \shortcut idx
  53. \color #008000
  54. \end_index
  55. \leftmargin 1.5cm
  56. \topmargin 3cm
  57. \rightmargin 1.5cm
  58. \bottommargin 3cm
  59. \secnumdepth 3
  60. \tocdepth 3
  61. \paragraph_separation indent
  62. \paragraph_indentation default
  63. \quotes_language french
  64. \papercolumns 1
  65. \papersides 2
  66. \paperpagestyle fancy
  67. \tracking_changes false
  68. \output_changes false
  69. \html_math_output 0
  70. \html_css_as_file 0
  71. \html_be_strict false
  72. \end_header
  73. \begin_body
  74. \begin_layout Subsection
  75. Tools > Statistics > Statistical test
  76. \end_layout
  77. \begin_layout Standard
  78. \begin_inset CommandInset label
  79. LatexCommand label
  80. name "subsection:statisticalTest"
  81. \end_inset
  82. \end_layout
  83. \begin_layout Standard
  84. \begin_inset Float figure
  85. placement !htb
  86. wide false
  87. sideways false
  88. status open
  89. \begin_layout Plain Layout
  90. \align center
  91. \begin_inset Graphics
  92. filename ../images/Partie3_Fonctions/statisticalTestDlg.png
  93. width 35text%
  94. \end_inset
  95. \begin_inset Caption
  96. \begin_layout Plain Layout
  97. \begin_inset CommandInset label
  98. LatexCommand label
  99. name "fig:statisticalTestDlg"
  100. \end_inset
  101. Boite de dialogue pour le test statistique local par rapport à une loi (normale)
  102. \end_layout
  103. \end_inset
  104. \end_layout
  105. \end_inset
  106. \end_layout
  107. \begin_layout Standard
  108. \begin_inset Index idx
  109. status collapsed
  110. \begin_layout Plain Layout
  111. champ scalaire
  112. \end_layout
  113. \end_inset
  114. \begin_inset Index idx
  115. status collapsed
  116. \begin_layout Plain Layout
  117. statistiques!test
  118. \end_layout
  119. \end_inset
  120. \begin_inset Index idx
  121. status collapsed
  122. \begin_layout Plain Layout
  123. statistiques!paramètres
  124. \end_layout
  125. \end_inset
  126. \begin_inset Index idx
  127. status collapsed
  128. \begin_layout Plain Layout
  129. Gauss|see
  130. \begin_inset ERT
  131. status collapsed
  132. \begin_layout Plain Layout
  133. {
  134. \end_layout
  135. \end_inset
  136. loi normale
  137. \begin_inset ERT
  138. status collapsed
  139. \begin_layout Plain Layout
  140. }
  141. \end_layout
  142. \end_inset
  143. \end_layout
  144. \end_inset
  145. Cette fonction permet d'appliquer un test du
  146. \begin_inset Formula $\chi^{2}$
  147. \end_inset
  148. \shape italic
  149. local
  150. \shape default
  151. sur un nuage de point muni d'un champ scalaire : le test du
  152. \begin_inset Formula $\chi^{2}$
  153. \end_inset
  154. est appliqué à chaque point, en confrontant l'histogramme
  155. \begin_inset Index idx
  156. status collapsed
  157. \begin_layout Plain Layout
  158. histogramme
  159. \end_layout
  160. \end_inset
  161. des valeurs scalaires du point et de ses
  162. \begin_inset Formula $n$
  163. \end_inset
  164. voisins (
  165. \begin_inset Formula $n$
  166. \end_inset
  167. étant un des paramètres de l'algorithme) avec la distribution spécifiée
  168. par l’utilisateur.
  169. \begin_inset Newline newline
  170. \end_inset
  171. \end_layout
  172. \begin_layout Standard
  173. Avant de spécifier les paramètres, l'utilisateur doit choisir le type de
  174. distribution (il a le choix actuellement entre
  175. \emph on
  176. Gauss
  177. \emph default
  178. et
  179. \emph on
  180. Weibull
  181. \emph default
  182. , les deux distributions étant aussi proposées dans l'outil
  183. \shape italic
  184. Compute stat.
  185. params
  186. \shape default
  187. - section
  188. \begin_inset space ~
  189. \end_inset
  190. \begin_inset CommandInset ref
  191. LatexCommand ref
  192. reference "subsection:computeStatParams"
  193. \end_inset
  194. ).
  195. Le résultat est un nouveau champ scalaire (une valeur pour chaque point
  196. - la métrique du
  197. \begin_inset Formula $\chi^{2}$
  198. \end_inset
  199. - qui donne une information sur la concordance locale entre la valeur scalaire
  200. active lors du test et la distribution testée).
  201. La théorie du test du
  202. \begin_inset Formula $\chi^{2}$
  203. \end_inset
  204. nous fournit un seuil (calculé à partir de la marge d'erreur
  205. \begin_inset Formula $p(\chi^{2})$
  206. \end_inset
  207. , dernier paramètre de l'algorithme) qui permet de classer les points en
  208. fonction de leur non-appartenance à la loi testée.
  209. Cette loi représentera typiquement le bruit de mesure, et on obtiendra
  210. ainsi l'ensemble des points dont la distance (à l'autre nuage/maillage)
  211. ne fait pas partie du bruit de mesure.
  212. Autrement dit, on pourra considérer que ces points correspondent effectivement
  213. à une différence, un changement, et on évitera de prendre en compte les
  214. autres points qui sont potentiellement des
  215. \shape italic
  216. faux positifs
  217. \shape default
  218. (immobiles mais dont la distance n'est pas nulle car elle est bruitée).
  219. Une fois le nuage séparé en deux classes, on peut garder le groupe des
  220. points
  221. \emph on
  222. hors distribution
  223. \emph default
  224. (voir figure
  225. \begin_inset space ~
  226. \end_inset
  227. \begin_inset CommandInset ref
  228. LatexCommand ref
  229. reference "fig:statisticalTestExample"
  230. \end_inset
  231. , en rouge) et les segmenter
  232. \begin_inset Index idx
  233. status collapsed
  234. \begin_layout Plain Layout
  235. segmentation
  236. \end_layout
  237. \end_inset
  238. par exemple en fonction de la proximité relative des points (par une extraction
  239. des composantes connexes - Cf.
  240. section
  241. \begin_inset space ~
  242. \end_inset
  243. \begin_inset CommandInset ref
  244. LatexCommand ref
  245. reference "subsection:labelConnectedComponents"
  246. \end_inset
  247. ).
  248. \begin_inset Index idx
  249. status collapsed
  250. \begin_layout Plain Layout
  251. composantes connexes
  252. \end_layout
  253. \end_inset
  254. \end_layout
  255. \begin_layout Standard
  256. \begin_inset Float figure
  257. placement !htb
  258. wide false
  259. sideways false
  260. status open
  261. \begin_layout Plain Layout
  262. \align center
  263. \begin_inset Graphics
  264. filename ../images/Partie3_Fonctions/statisticalTestExample.jpg
  265. width 90text%
  266. \end_inset
  267. \begin_inset Caption
  268. \begin_layout Plain Layout
  269. \begin_inset CommandInset label
  270. LatexCommand label
  271. name "fig:statisticalTestExample"
  272. \end_inset
  273. Champ des écarts initial (en haut à gauche), filtrage statistique (en haut
  274. à droite), puis extraction des points
  275. \emph on
  276. hors distribution théorique
  277. \emph default
  278. (en bas à gauche) et enfin extraction des composantes connexes (en bas
  279. à droite).
  280. \end_layout
  281. \end_inset
  282. \end_layout
  283. \end_inset
  284. \end_layout
  285. \begin_layout Standard
  286. Remarques :
  287. \end_layout
  288. \begin_layout Itemize
  289. Pour appeler cette fonction, il faut sélectionner une seule entité 3D, munie
  290. d'un champ scalaire actif.
  291. \end_layout
  292. \begin_layout Itemize
  293. Pour régler le paramètre
  294. \begin_inset Formula $p(\chi^{2})$
  295. \end_inset
  296. , il est important de comprendre que le test du
  297. \begin_inset Formula $\chi^{2}$
  298. \end_inset
  299. permet uniquement de rejeter l'hypothèse selon laquelle
  300. \emph on
  301. \begin_inset Quotes fld
  302. \end_inset
  303. les valeurs du champ scalaire prises sur le voisinage de chaque point suivent
  304. la loi testée
  305. \emph default
  306. \begin_inset Quotes frd
  307. \end_inset
  308. , mais pas l'inverse.
  309. Ainsi, plus la marge d'erreur est faible, et plus le seuil du
  310. \begin_inset Formula $\chi^{2}$
  311. \end_inset
  312. sera grand (on rejette moins souvent l'hypothèse citée précédemment, et
  313. on classe donc moins de points comme
  314. \emph on
  315. ne suivant pas la loi testée
  316. \emph default
  317. ).
  318. \end_layout
  319. \begin_layout Itemize
  320. \begin_inset ERT
  321. status open
  322. \begin_layout Plain Layout
  323. \backslash
  324. textcolor{red}
  325. \end_layout
  326. \end_inset
  327. \begin_inset ERT
  328. status collapsed
  329. \begin_layout Plain Layout
  330. {
  331. \end_layout
  332. \end_inset
  333. Inversement, plus
  334. \begin_inset Formula $p(\chi^{2})$
  335. \end_inset
  336. est grand, plus on aura de points "hors la loi" (colorés en rouge dans
  337. la figure
  338. \begin_inset space ~
  339. \end_inset
  340. \begin_inset CommandInset ref
  341. LatexCommand ref
  342. reference "fig:statisticalTestExample"
  343. \end_inset
  344. ).
  345. \begin_inset ERT
  346. status collapsed
  347. \begin_layout Plain Layout
  348. }
  349. \end_layout
  350. \end_inset
  351. Notez que ce paramètre sert uniquement à pré-positionner les potentiomètres
  352. de réglage des fausses couleurs (seuils de coupure et de saturation des
  353. valeurs du champ scalaire) pour l'affichage du résultat à l'écran (Cf.
  354. section
  355. \begin_inset space ~
  356. \end_inset
  357. \begin_inset CommandInset ref
  358. LatexCommand ref
  359. reference "Champs-scalaires"
  360. \end_inset
  361. ).
  362. Ces potentiomètres peuvent être ensuite déplacés par l'utilisateur avant
  363. extraction effective des points (par appel de la fonction
  364. \emph on
  365. Scalar Fields > Filter by Value
  366. \emph default
  367. , qui va créer un nouveau nuage de points ne comportant que les points présentem
  368. ent affichés à l'écran, c.à.d.
  369. les points ne suivant pas la distribution théorique).
  370. De plus, la distance du
  371. \begin_inset Formula $\chi^{2}$
  372. \end_inset
  373. est extrêmement divergente et ceci donne une grande marge de manœuvre à
  374. l'algorithme.
  375. Ainsi, une modification relativement grande du seuil de coupure n'aura
  376. que peut d'effet sur la classification.
  377. Au pire, on risque de rater un tout petit nombre de points (au niveau des
  378. bordures des zones limites).
  379. \end_layout
  380. \begin_layout Itemize
  381. Pour obtenir des résultats précis, il faut par contre connaitre ou mesurer
  382. la distribution du bruit de mesure (une sorte de bruit moyen, en première
  383. approximation, comprenant les erreurs de mesure due au capteur lui-même
  384. et aux conditions de mesure - nature de la surface, lumière, température,
  385. etc.
  386. - ou aussi les erreurs provenant de la création du maillage dans le cas
  387. d'une comparaison nuage/maillage par exemple, etc.).
  388. Les paramètres de la distribution statistique correspondante peuvent donc
  389. être définis à partir de connaissances a priori mais peuvent aussi être
  390. déterminés à partir des valeurs scalaires associées à une portion du nuage
  391. typiquement (voir la fonction
  392. \shape italic
  393. Compute stat.
  394. params
  395. \shape default
  396. - section
  397. \begin_inset space ~
  398. \end_inset
  399. \begin_inset CommandInset ref
  400. LatexCommand ref
  401. reference "subsection:computeStatParams"
  402. \end_inset
  403. ).
  404. \end_layout
  405. \begin_layout Itemize
  406. L'algorithme crée un nouveau champ scalaire nomme (
  407. \emph on
  408. Chi2 Distances
  409. \emph default
  410. ).
  411. Ce champ est ajouté au nuage de points courant.
  412. Si celui-ci existe déjà il sera écrasé.
  413. \end_layout
  414. \end_body
  415. \end_document