! rd_i12 read images in raw format and determine if it''s 8 or 12 bits coded
! Copyright (C) 1999,2010 Gauthier Delerce
!
! This program is free software; you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation; either version 2, or (at your option)
! any later version.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program; if not, write to the Free Software Foundation,
! Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
!
! Gauthier Delerce <gauthier@delerce.fr>
subroutine rd_i12(filename,image,m,ierr)
!****************************************************
! PORTABILITY PROBLEMS
! UNIX AND WINDOWS SYNTAXES USED ARE NOT COMPATIBLES
! SPEED UP , CHECK 16 BITS PLEASE
!****************************************************
implicit none
integer::m
!integer(2),dimension(m)::img2
integer,dimension(m)::image,image2
integer(4),dimension(m)::img
integer::ierr,i,j,bpp,seed,inde
integer(1),dimension(2*m)::in
character(256)::filename
real::rnunf
!external rnunf
inde=scan(filename,' ')
inde=inde-1
!image=0.
!image2=0.
!bpp=1
ierr=0
!open(unit=10,file=filename,access='direct',recl=m,form='binary',status='old',ACTION='READ',err=10)
open(unit=10,file=filename,access='direct',recl=m,status='old',ACTION='READ',err=10)
goto 20
10 ierr=1
write(*,*) 'Error reading image'
return
20 continue
print*, 'Reading image ',filename(1:inde)
read(10,rec=1,err=10) in(1:m)
read(10,rec=2,err=11) in(m+1:2*m)
goto 12
11 print*,'8 bits image '
image(1:m)=in(1:m)
forall(i=1:m,image(i)<0.)
image(i)=image(i)+256
endforall
bpp=1
!print*,image
goto 155
12 bpp=2
print*,' Image 12 bits'
image(1:m)=in(1:2*m:2)
forall(i=1:m,image(i)<0.)
image(i)=image(i)+256
endforall
image2(1:m)=in(2:2*m:2)
forall(i=1:m,image2(i)<0.)
image2(i)=image2(i)+256
endforall
image=image+image2*256
155 close(10)
print*,'Intensite min/max: ',minval(image),maxval(image)
print*,''
return
end