Simulation radiographique

ça c'est vrai pour améliorer la qualité des détails...

Mais moi j'ai des énormes bandes ignobles qui dépendent du point de vue selon lequel on se place (les artefacts sur une sphère ne sont pas les mêmes selon les points de vue alors que c'est isotrope : ça varie de lignes à cercles concentriques en passant par des trucs indesciptibles) !

plus étrange encore, j'ai tout refait car je trouvais pas le bug, et les artefacts sont à présent les même mais que sur une moitié de la sphère !!!
Truc incompréhensible... L'autre moitié est parfaite et totalement intacte...

Sûr c'est un bug de mon programme, mais où, mystère.... Je trouve pas !!!!!
la seule piste que j'ai, c'est que si au lieu de faire:

Code:

---

1 2 3

Pour tout échantillon de la droite Ajouter à l'absorption (IntensitéImage par interpolation trilinéaire)*pas FinPour

---

je fais:

Code:

---

1 2 3

Pour tout échantillon de la droite Ajouter à l'absorption "pas" FinPou

---

je n'ai plus les artefacts (je commence à échantillonner que sur mon modèle 3D par l'étude de l'intersection du rayon avec le modèle)..

Donc ça doit venir de l'interpolation sauf que je l'ai reprise d'un de mes autres programmes qui fonctionne parfaitement!!!!

Non de D...

Cela ressemble à 1 problème de mise en forme du soft.
Lorsque l'on considère une droite joignant le point dont on veut calculer la couleur ( Intensité transmise des RX) à la source, on doit intégrer en faite un ensemble de droite allant de ce point et scannant à " l'autre bout" la source RX
On peut pondérer l''influence de chaque droite si la source n'émet pas autant de photon X sur chacun de ses élément de surface et si l'émission n'est pas isotrope.
Chaque droite "entre" dans le modèle " en un point P1 et "en sort" en P2. Elle le traverese sur une longure L qui varie avce le scan de la source et le scan de l'image à créer. Losque ses droites croisent le modèle la densité de points calculés peut changer et il faut pondérer correctement l'influence de l'attenuation avec la "porté" d'un élément. Cela peut conduire à des sauts dans le résultat. Je préfèrerais avoir une densité de point constante sur la droite et, à chaque noeud sur la droite, calculer par interpolation la valeur de mu en ce point (si la source RX n'est pas monochromatique il faut encre intégrer sur lambda car Mu = Mu(lambda) ). + simplement on peut aussi prendre la valeur d'absobtion la + proche dans le modèle; tout dépand de la densité de points sur la droite en regard de la taille des objets et de la qualité du maillage du modèle.
De toute facon dans ce deuxieme cas, ( densité linéaire de point = constante), Il y a un bruit de quatification de +- 1 point. ( il est en général + important dans la 1ere approche ).
Ce bruit peu être parfois genant et necessite parfoi un filtrage du résultat.
Pour accroitre la qualité on peut aussi
- faire des corrections d'intercept
- l'utilisation des ondelettes permet d'ameliorer les contours
- ....

Bonjour,
Je remets au goût du jour simplement pour vous signaler que j'ai enfin résolu le problème : ça n'était pas dû à un aliasing (que j'évite en faisant un échantillonnage suffisamment fin de ma droite) mais à un simple bug de programmation : j'avais une macro IND sensée me calculer les indices de mes tableaux qui fonctionnait mal...

Vous saviez vous que

Code:

---

IND((int)x+1,y,wx)

---

n'est pas pareil que

Code:

---

IND((int)(x+1),y,wx)

---

...
Halluciant quand même !! Tant d'heures de recherche du pb pour une simple macro qui prend mal les paramètres !!!

Aaargh!!!

Sauf si tu mets adroitement tes parenthéses à la definition de ta macro

Code:

---

1 2	#define TOTO(a) f( (a) )

---

Je me disais bien que ce genre d'images, je les connaissais ;)
En plus avec la connerie que j'ai dite dans l'autre topic, j'ai honte :D
J'ai peut-être surveillé certains de tes exams, Mathieu...

???
Mais qui es-tu ô étrange créature ??? :lol:
???

Sinon ok pour les macros je savais pas...

A+

Bonne question...
Tiens, j'ai appris qu'un 3A de l'ENSPS allait faire une thèse au LSIIT et à l'IPB, mais je ne sais pas qui exactement, ça s'est fait hier, en fait, t'es au courant - question sous-jacente : est-ce toi ? -

Alors là t vraiment très fort !!!
Oui c'est moi !!!

Allez, dis moi qui tu es, trop de suspense tue le suspense !!!! :lol:

Je suis ton futur collègue de bureau ;)
Trop trop fort :D