--- title: "TNSI: Traitement numerique du signal" date: 2021-11-10 09:00 categories: [Image S9, TNSI] tags: [Image, S9, TNSI] math: true --- Lien de la [note Hackmd](https://hackmd.io/@lemasymasa/BJxQ-Ztvt) # Presentation Parcours: - Ingenieur topographe - INSA - These en traitement du signal et de l'image - Telecom Ingenieur Chercheur a EDF R&D (Saclay) - Groupe realite virtuelle et visualisation scientifique (avec Arnaud MAS) - Numerisation et apprentissage statistique (image et 3D) ## Groupe Realite Virtuelle Visualisation Scientifique - Aider l'exploitant des sites de production nucleaire a decider lors des phases de preparation et de realisation des arretes de tranche - Aider ## OCR pour les plans techniques ![](https://i.imgur.com/rroAJpP.png) ## La reconnaissance de forme dans les images ![](https://i.imgur.com/IZye8az.jpg) Segmentation semantique: deux reseau - Segmentation pixelique - Segmentation semantique pour classifier les pixels de l'image ![](https://i.imgur.com/lq7hSnj.png) ## La reconnaissance de forme *dans les nuages de points* ![](https://i.imgur.com/kAh6IEB.jpg) > Donnees DP2D FES2 > Local R250 > 280 millions de points Actuellement: manuel - Segmentation en petit nuage de points + ajustement de forme # Traitement numerique du signal :::info **Definition** Representation de la variation d'un phenomene physique ::: > Exemples > Evolution de la temperature ou de la pression dans le temps :::info **Definition** Transcrire numeriquement (i.e. des donnees) un signal continu (monde reel) ::: On veut passer du monde continu au monde discret ![](https://i.imgur.com/cklHObs.png) *Comment on fait le passage du continu au discret ?* > On va faire un echantillonage en temps et en amplitude *Quelle est la precision de cette discretisation ?* > On utilise le theoreme de Shannon En resume: on prend notre signal, on compare dans chaque base de fonction a quel point notre signal ressemble et on va pouvoir voir a quel point notre signal est haute frequence et basse frequence. ![](https://i.imgur.com/TvaWcT2.png) :::danger Avec les transformee de Fourier, on veut passer en temporel **sans perdre d'information** ::: ## Theoreme d'interpolation :::info Decoule du theoreme de Shannon ::: On utilise un *sinus cardinal* Le theoreme de Shannon nous dit qu'on a un signal continu: $$ x(t) = \sum_{n=-\infty}^{+\infty}x[n]\underbrace{\frac{\sin\biggr(\frac{\pi(\overbrace{t-\color{red}{nT_e}}^{\text{translation}})}{\underbrace{\color{red}{T_e}}_{\text{echelle}}}\biggr)}{\frac{\pi(t-nT_e)}{T_e}}}_{\color{red}{sinc}} $$ ![](https://i.imgur.com/xPCan8j.png) ## Exercice $$ x(t) = \sin(10t) + 2\sin(2t) + \sin(5t) - 3\sin(\frac{t}{2})\\ t = [-2, 2]\to 400 $$ 1. Normaliser et centrer 2. $T_e=0.1$ 3. Interpoler avec $sinc$ On reprendre notre figure ![](https://i.imgur.com/rWmIKr9.png) On prend des echantillons a intervalles regulire ![](https://i.imgur.com/ZdBEWja.png) On va sommer les sinus cardinaux: ![](https://i.imgur.com/1FBfJc8.png) Ca va nous permettre de reconstruire notre signal: ![](https://i.imgur.com/fplEQWb.png) *Quel est l'interet du sinus cardinal ?* > On peut utiliser n'importe quelle interpolation, mais le sinus cardinal est le meilleur ## Quantification :::info **Quantification scalaire**: arrondir a l'entier le plus proche ::: ![](https://i.imgur.com/FPDCbTe.png) - $x$: amplitude des valeurs - $q(x)$: valeur de quantification ## Traitement *Comment on debruite un signal ?* > Convolution avec une fonction gaussienne ? > Convolution avec une porte ? > Non local means ? (c'est un flou gaussien ou un flou uniforme) On a un signal avec un flou gaussien: ![](https://i.imgur.com/4EWwML3.png) :::success On fait une convolution avec une porte ::: ## En resume - Comment echantilloner et reconstruire un signal - Comment analyser un signal - Comment filtrer une partie de l'information d'un signal # Transformee de Fourier 1. Rappel rapide 2. Analyser harmonique 1. Decomposition en serie de Fourier 2. Discrete Time Fourier Transform (DTFT) - Tf a temps discret 3. Discrete Fourier Transform (DFT) - TF discrete 4. Continuous Fourier Transform (CFT) 5. Fast Fourier Transform (FFT) ## Produit scalaire *Comment est-ce qu'on calcule un produit scalaire ?* $$ \langle x, y\rangle = \sum_{n=-\infty}^{+\infty}x[n]y[n]^*\\ \langle x, y\rangle = \int_{-\infty}^{+\infty}x(t)y^*(t)dt $$ :::warning Avec Fourier, on est en complexes ::: $$ \Vert x\Vert^2 = \langle x,x\rangle = \sum_{n=-\infty}^{+\infty}\vert x[n]\vert $$ ## Decomposition/reconstruction Soit $$\{f_n\}_{n\in\mathbb N}\to$$ base orthogonale. $$ \langle f_n,f_p\rangle = 0\quad n\neq p $$ $\exists$ une suite $\lambda[n]$ telle que $\lim_{N\to+\infty}\Vert x-\sum\lambda[n]f_n\Vert = 0$ $$ x=\sum_{n=0}^{+\infty}\lambda_n f_n\quad\text{avec }\lambda_n = \frac{\langle x,f_n\rangle}{\Vert f_n\Vert^2} $$ ## Exercice $$ x = \sin(2\pi t) + 2\sin(3\times 2\pi t) - 3\sin(5\times 2\pi t) + \sin(7\times 2\pi t)\\ t = [0,1] $$ 1. - Tracer $x$ avec 1000 echantillons - Decomposition de $x$ avec $$\{\sin(n\cdot 2\pi t)\}_{n\in N}\to$$ $N = [0,...,5]$ ou $N = [0,...,20]$ - Verifier l'orthogonalite - Tracer les coefficients - Reconstruction 2. Si on ajoute une phase au sinus ? - Idem