#  This software and supporting documentation are distributed by
#      Institut Federatif de Recherche 49
#      CEA/NeuroSpin, Batiment 145,
#      91191 Gif-sur-Yvette cedex
#      France
#
# This software is governed by the CeCILL-B license under
# French law and abiding by the rules of distribution of free software.
# You can  use, modify and/or redistribute the software under the 
# terms of the CeCILL-B license as circulated by CEA, CNRS
# and INRIA at the following URL "http://www.cecill.info". 
# 
# As a counterpart to the access to the source code and  rights to copy,
# modify and redistribute granted by the license, users are provided only
# with a limited warranty  and the software's author,  the holder of the
# economic rights,  and the successive licensors  have only  limited
# liability.
#
# In this respect, the user's attention is drawn to the risks associated
# with loading,  using,  modifying and/or developing or reproducing the
# software by the user in light of its specific status of free software,
# that may mean  that it is complicated to manipulate,  and  that  also
# therefore means  that it is reserved for developers  and  experienced
# professionals having in-depth computer knowledge. Users are therefore
# encouraged to load and test the software's suitability as regards their
# requirements in conditions enabling the security of their systems and/or 
# data to be ensured and,  more generally, to use and operate it in the 
# same conditions as regards security.
#
# The fact that you are presently reading this means that you have had
# knowledge of the CeCILL-B license and that you accept its terms.
from soma import aims
import sys

app = aims.AimsApplication( sys.argv, 'AimsData test in python' )

app.initialize()

reader = aims.Reader()
# read volumes as AimsData
reader.mapType( 'Volume', 'AimsData' )
data = reader.read( '/home/riviere/data/irm.ima' )
print 'data:', data
h = data.header()
print 'header:', h

dx = data.dimX()
dy = data.dimY()
dz = data.dimZ()
dt = data.dimT()

mi = data.value( 0 )
ma = data.value( 0 )

dims = [ dx, dy, dz, dt ]
print 'volume_dimensions:', dims
print 'voxel_size:', [ data.sizeX(), data.sizeY(), data.sizeZ(), data.sizeT() ]

print 'C++ min/max:', data.minimum(), '-', data.maximum()
print 'setting value 1618 at pos (125, 63, 12 )'
data.setValue( 1618, 125, 63, 12 )
print 'C++ min/max:', data.minimum(), '-', data.maximum()

print 'minIndex:', data.minIndex()
print 'maxIndex:', data.maxIndex()

print 'python min/max...'
for t in xrange( dt ):
    for z in xrange( dz ):
        print 't:', t, 'z:', z
        for y in xrange( dy ):
            for x in xrange( dx ):
                v = data.value( x, y, z, t )
                if mi > v:
                    mi = v
                if ma < v:
                    ma = v

print 'min:', mi
print 'max:', ma

fileout = '/tmp/toto.vimg'
print 'writing volume to', fileout
w = aims.Writer()
w.write( data, fileout )
print 'object type:', w.writtenObjectType()
print 'data type:', w.writtenObjectDataType()
print 'full type:', w.writtenObjectFullType()