pour voir en enseble de slices d'une tomo
-----------------------------------------------
MLab.ShowTomoSlices("90s4_0002_bis.vol", range(0,2048,100) )


pour faire un simple graph
--------------------------------

 
x= Numeric.arrayrange(0,12.0,0.01)
y=Numeric.sin(x)
w=MLab.Plot_1d((y,x ) )

print w

# w est un dictionnaire qui permet une personalisation poussee


pour utiliser un filtre existant:
example avec AzimIntFilter pour transformer une image en coordonnees polaires
------------------------------------------------------------------------------

import PyDVT.EdfFileData as EdfFileData
data=EdfFileData.EdfFileData()
data.SetSource('circle.edf')
data.LoadSource()

import DataSelection
sel=DataSelection.RectSelection(data)
import ImageView
Filter=ImageView.ColormapFilter(None,sel)
import AzimIntFilter
azim=AzimIntFilter.AzimIntFilter(Filter)

import MLabData
mldata=MLabData.MLabData(azim)

mldata=MLabData.MLabData(azim.GetOutput()['data'])
MLab.newWindow(mldata)


ou bien:
en utilisant la nouvelle classe MLabEdfFileData: (18/02/03)
------------------------------------------------------------------------------

import MLabEdfFileData
data=MLabEdfFileData.MLabEdfFileData('circle.edf')
data.display()

azim=data.applyFilter('AzimIntFilter')
arr=azim.GetOutput()['data']
import MLabData
w=MLab.newWindow(MLabData.MLabData(arr))


Pour appliquer le filtre sur une partie de l'image 
et pour passer les parametres au filtre en utilisant le dictionnaire 'pars' (26/02/03)
--------------------------------------------------------------------------------------

import MLabEdfFileData
data=MLabEdfFileData.MLabEdfFileData('circle.edf')
arr=data.GetPageArray(0)[200:300,200:300]
data=MLabEdfFileData.MLabEdfFileData()
data.AppendPage({},arr)
data.display()

azim=data.applyFilter('AzimIntFilter',pars={'center':(56,56)})
arr=azim.GetOutput()['data']
import MLabData
w=MLab.newWindow(MLabData.MLabData(arr))   


Pour lire un stack d'images, appliquer un filtre sur ce stack (ici on cherche le maximum de l'image
et ses coordonnees) et ecrire les resultats dans un fichier ascii (19/03/03)
--------------------------------------------------------------------------------------

t1=time.time()
import MLabRead
dict1=MLabRead.MLabRead('/data/opid02/inhouse/pierre/beamprofile', None, 'beam.*ccd', 'EDF',None,(1050,1050),(100,100))
import MLabApplyFilter
import FindMaxFilter
dict2=MLabApplyFilter.MLabApplyFilter(FindMaxFilter.FindMaxFilter,dict1)
import MLabAsciiWrite
MLabAsciiWrite.MLabAsciiWrite('.','test',dict2)
t2=time.time()
print t2-t1 
