C++*******************************************************************
C
C ROT32.F                          USED ALLOCATE AUGUST 00 ARDEAN LEITH
C
C **********************************************************************
C=* FROM: SPIDER - MODULAR IMAGE PROCESSING SYSTEM.   AUTHOR: J.FRANK  *
C=* Copyright (C) 1985-2005  Health Research Inc.                      *
C=*                                                                    *
C=* HEALTH RESEARCH INCORPORATED (HRI),                                *   
C=* ONE UNIVERSITY PLACE, RENSSELAER, NY 12144-3455.                   *
C=*                                                                    *
C=* Email:  spider@wadsworth.org                                       *
C=*                                                                    *
C=* This program is free software; you can redistribute it and/or      *
C=* modify it under the terms of the GNU General Public License as     *
C=* published by the Free Software Foundation; either version 2 of the *
C=* License, or (at your option) any later version.                    *
C=*                                                                    *
C=* This program is distributed in the hope that it will be useful,    *
C=* but WITHOUT ANY WARRANTY; without even the implied warranty of     *
C=* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU  *
C=* General Public License for more details.                           *
C=*                                                                    *
C=* You should have received a copy of the GNU General Public License  *
C=* along with this program; if not, write to the                      *
C=* Free Software Foundation, Inc.,                                    *
C=* 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.      *
C=*                                                                    *
C **********************************************************************
C
C    PURPOSE: 
C     ROTATE AN IMAGE BY AN ARBITRARY ANGLE OF
C     DEGREE BETWEEN 0 AND 360.  IT IS CALLED WHEN THE IMAGE IS SMALL
C     ENOUGH TO FIT IN THE BUFFER.
C
C    ROT32(LUNI,LUNO,NSAM,NROWS,NROWE,NROWSK,BUF,THETA,BACK,SHX,SHY)
C
C    PARAMETERS:
C
C       LUNI         LOGICAL UNIT NUMBER OF INPUT IMAGE
C       LUNO         LOGICAL UNIT NUMBER OF OUTPUT IMAGE
C	NSAM	     ROW LENGTH
C	NROWS,NROWE  STARTING AND ENDING ROW
C       NROWSK       SKIPPING FACTOR FOR ROWS
C                    (A NEGATIVE VALUE MEANS NONSEQUENTIAL INPUT AND OUTPUT) 
C       THETA        ROTATION ANGLE IN RADIANS
C       BACK         AVERAGE OF INPUT IMAGE
C       SHX,SHY      ORIGIN SHIFT
C
C       POSITIVE THETA: COUNTER-CLOCKWISE ROTATION
C                       (MATHEMATICALLY POSITIVE DIRECTION)
C
C        AN IMAGE CAN BE A SLICE OF A THREE-DIMENSIONAL DENSITY
C        DISTRIBUTION. FOR THIS REASON, A STARTING ROW UNEQUAL
C        TO 1 AND AN ENDING ROW UNEQUAL TO NROW IS PERMITTED.
C
C--*******************************************************************

	SUBROUTINE ROT32
     &     (LUNI,LUNO,NSAM,NROWS,NROWE,NROWSK,THETA,BACK,SHX,SHY)

        INCLUDE 'CMBLOCK.INC'

        INCLUDE 'CMLIMIT.INC'
        COMMON /IOBUF/ RBUF(NBUFSIZ)
 
C       ACTUAL DIMENSIONS OF BUF ARE NSAM + NSAM * NROW,  WHERE
C       NROW = ((NROWE-NROWS) / NROWSK) + 1
        REAL, ALLOCATABLE, DIMENSION(:) :: BUF
  
	DATA PI/3.14159/

        IFLAG1 = 0
        IF (NROWSK.LT.0) IFLAG1 = 1
        IF (NROWSK.LT.0) NROWSK = -NROWSK

        IF (THETA .EQ. 0.) THEN
C          ZERO DEGREE ROTATION
           II = 0
	   DO  I=NROWS,NROWE,NROWSK
              IF (IFLAG1 .EQ. 0) II = II + 1
              IF (IFLAG1 .EQ. 1) II = I
              CALL REDLIN(LUNI,RBUF,NSAM,I)
	      CALL WRTLIN(LUNO,RBUF,NSAM,II)
	   ENDDO
	   RETURN
        ENDIF

 	NSAMH = NSAM/2
	NROW  = ((NROWE-NROWS) / NROWSK) + 1
	NROWH = NROW/2
	KCENT = NSAMH+1
	ICENT = NROWH+1

        ALLOCATE (BUF(NSAM*NROW), STAT=IRTFLG)
        IF (IRTFLG .NE. 0) THEN 
            CALL ERRT(46,'ROT32, BUF',IER)
            RETURN
        ENDIF     

10	IF (THETA .GT. PI)  THETA = -2.*PI+THETA
	IF (THETA .LT. -PI) THETA = 2.*PI+THETA
	COD = COS(THETA)
	SID = SIN(THETA)

C       READ IN WHOLE INPUT IMAGE
        J = 0
	DO  I = NROWS,NROWE,NROWSK
           J = J + 1
           L = (J-1)*(NSAM)+1
           CALL REDLIN(LUNI,BUF(L),NSAM,I)        
	ENDDO

C       NOW GO THROUGH OUTPUT COO SYSTEM; COMPUTE, FOR EACH POINT IN
C       ROW, THE POSITION IN THE OLD COO SYSTEM.
C       THEN CALCULATE POINT FROM FOUR SURROUNDING POINTS USING BILINEAR
C       INTERPOLATION. WRITE OUT EACH LINE AS YOU GO ALONG.

C       JUST TO ALLOW A CHANGE IN THE ROTATIONAL CENTER
C       (RELATIVE TO THE CENTRAL PIXEL)

        RICENT = ICENT + SHY
        RKCENT = KCENT + SHX  

        JJ     = 0
	DO  I = 1,NROW
           JJ = JJ+1
           IF (IFLAG1 .EQ. 1) II = NROWS + (I-1) * NROWSK
           IF (IFLAG1 .EQ. 0) II = I
           Y    = I - RICENT
           YCOD = Y * COD + RICENT
           YSID = -Y * SID + RKCENT
           DO K = 1,NSAM
              RBUF(K) = BACK
              X       = K - RKCENT
              XOLD    = X * COD + YSID
              YOLD    = X * SID + YCOD
              IYOLD   = YOLD
              YDIF    = YOLD - IYOLD
              YREM    = 1.   - YDIF
              IXOLD   = XOLD
              IF ((IYOLD .GE. 1 .AND. IYOLD .LE. NROW-1) .AND.
     &            (IXOLD .GE. 1 .AND. IXOLD .LE. NSAM-1)) THEN
c                INSIDE BOUNDARIES OF OUTPUT IMAGE
                 XDIF    = XOLD - IXOLD
                 XREM    = 1. - XDIF
                 NADDR   = (IYOLD-1) * NSAM + IXOLD 
                 RBUF(K) = YDIF*(BUF(NADDR+NSAM)*XREM
     &                    +BUF(NADDR+NSAM+1)*XDIF)
     &                    +YREM*(BUF(NADDR)*XREM + BUF(NADDR+1)*XDIF)
              ENDIF
           ENDDO

           CALL WRTLIN(LUNO,RBUF,NSAM,II)
	ENDDO

        IF (ALLOCATED(BUF)) DEALLOCATE(BUF)

	RETURN
	END