C PROGRAM AVGN.FT C ---------------- C C PETER LEMKIN C IMAGE PROCESSING UNIT, DCBD C NATIONAL CANCER INSTITUTE C NATIONAL INSTITUTES OF HEALTH C 9000 ROCKVILLE PIKE C BETHESDA, MD. 20014 C C FEB 16, 1977 C JAN 7, 1977 C DEC 14, 1976 C C C C INTRODUCTION C ------------ C AVGN.FT IS A CHAINED PROGRAM USED WITH BMON2 C IT RECEIVES ITS ARGUMENTS FROM THE CD AREA AND THE IBM1,IHGH1 C IBM2,IHGH2, JBM, JHGH VARIABLES IN COMMON. COMMON IS RESTORED C FIRST BEFORE THE FUNCTION (TO BE INSERTED INTO THE BODY) IS C EVALUATED. AFTER THE FUNCTION IS PERFORMED, COMMON C IS SAVED AND BMON2 IS CHAINED BACK TO. C ANY COMPUTATIONS ARE DONE, THE ARGUMENTS ARE THEN CHECKED C C AVGN COMPUTES THE AVERAGE OF EACH NXN NEIGHBORHOOD of AND C STORES THAT VALUE IN . C OPDEFS C ------ C C C C S OPDEF SWAB 7431 S OPDEF SWBA 7447 C S OPDEF DAD 7443 S OPDEF DST 7445 S OPDEF DCM 7575 S OPDEF MUY 7405 S OPDEF DVI 7407 C S OPDEF DISP2 6436 S OPDEF DISP1 6435 C S OPDEF CLAMQ 7621 S OPDEF MQA 7501 S OPDEF MQL 7421 S OPDEF KRS 6034 S OPDEF BSW 7002 C C DIMENSION KP(225),ISUM(2) C C C [1] INITIALIZATION WRITE(1,900) 900 FORMAT('AVGN 2/16/77 - 4:01PM') ET=TIMER(0) C C VERIFY THAT N=0 OR N IS ODD N=ICNUM NSQUARE=N*N C IF(N)995,100,101 100 N=3 101 IF(N-15)102,102,995 102 CONTINUE C C TEST IF N IS ODD S TAD \N S AND (1 S SNA CLA S JMP \995 /N IS NOT ODD ERROR C C C [2] VERIFY BM SPECS S JMS CKIN S JMS CKOUT C@@S JMS CKIN2 C C C C [3] PROCESS IMAGE DO 200 IY1=KY1,KY2 C S CPAGE 3 S JMS TTYCTL S JMP \998 /DONE C IY=IY1-1 S TAD \IY S DISP2 C C C [3.1] PROCESS A LINE DO 200 IX1=KX1,KX2 IX=IX1-1 S TAD \IX S DISP1 C C C C [3.1.1] GET BMI NEIGHBORHOOD MEM=IBM1 IBYTE=IHGH1 IZ=N C ISUM=0 S DCA \ISUM# CALL GETNXN(KP) C C C [3.1.2] COMPUTE: AVGN S JMS KERNAL C C C C [3.1.3] IZ==>BMJ MEM=JBM IBYTE=JHGH CALL PACK2D 200 CONTINUE C C C C C [999] RETURN GOTO 998 C C N SIZE ERROR 995 WRITE(1,994)N 994 FORMAT(' ',I5,' IS EITHER NOT ODD OR IS > 15') C 999 WRITE(1,996) 996 FORMAT('BAD BM SPECIFICATION!') 998 ET=TIMER(1) CALL CHAIN('BMON2') C*************************************************** C *PROCEDURE C K O U T C****************************************************** C C C CHECK WHETHER THE OUTPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKOUT, JMP I CKOUT S CKOUT, 0 /ENTRY C C [1] TEST IF KOUTFILE="BM" S TAD \KOUTFILE S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [2] TEST IF (KOUTFILE(2) LAND '7700)=DIGIT S TAD \KOUTFILE# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKOUT /OK. C C C*************************************************** C *PROCEDURE C K I N C****************************************************** C C C CHECK WHETHER THE INPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKIN, JMP I CKIN S CKIN, 0 /ENTRY C C [1] TEST IF BMI1="BM" S TAD \SFILE S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [2] TEST IF (BMI1(2) LAND '7700)=DIGIT S TAD \SFILE# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKIN /OK. C C C*************************************************** C *PROCEDURE C K I N 2 C****************************************************** C C C CHECK WHETHER THE INPUT BM SPEC IS LEGAL ELSE GOTO 999. S CPAGE 3 SRCKIN2, JMP I CKIN2 S CKIN2, 0 /ENTRY C C [1] CHECK FIRST IBM1 SPEC S JMS CKIN C C C [2] TEST IF BMI2="BM" S TAD \SEXT S CIA S TAD BMTEXT S SZA CLA S JMP \999 /FAILED C C [3] TEST IF (BMI2((6) LAND '7700)=DIGIT S TAD \SEXT# S BSW S AND (0077 S TAD (-60 /"0" S SPA S JMP \999 /NO, FAILED S TAD (-D8 /"7" TEST S SMA CLA S JMP \999 /FAILED S JMP RCKIN2 /OK. C C C ************************************************************ C SUBROUTINE: T T Y C T L (INTERNAL) C ************************************************************ C S CPAGE 3 S RTTYC, JMP I TTYCTL S TTYCTL, 0000 /ENTRY C S KSF /ANYTHING TYPED S JMP NORMAL /NO, RETURN NORMALLY S KRB /GET TYPED CHARACTER S AND (0177 /TAKE CARE OF PARITY PROBLEMS S TAD (-17 /TEST FOR CTRL/O S SNA /SKIP IF NOT CTRL/O S JMP RTTYC /ABORT CALLING ROUTINE (ERROR RETURN) S TAD (-4 /TEST FOR CTRL/S [-17-4=-23(OCTAL)] S SZA CLA /SKIP IF CTRL/S S JMP NORMAL /NOT CTRL/O OR CTRL/S SO RETURN NORMALLY C S SLEEP,KSF /WAIT FOR CTRL/Q S JMP SLEEP /KEEP WAITING S KRB /READ CHARACTER S AND (0177 S TAD (-17 /IS IT A CTRL/O? S SNA /SKIP IF NOT S JMP RTTYC /YES, ABORT S TAD (-2 /TEST FOR CTRL/Q (-17-2=-21 OCTAL) S SZA CLA /SKIP IF SO S JMP SLEEP /NOPE, KEEP SLEEPING C S NORMAL,INC TTYCTL /INCREMENT RETURN ADDRESS FOR NORMAL RETURN S CLA /SAFETY VALVE S JMP RTTYC /RETURN C **************************************************** C *SUBROUTINE K E R N A L C ***************************************************** C C COMPUTE IZ AS A FUNTION OF KP NEIGHBORHOOD. C C S CPAGE 3 S RKERNAL, JMP I KERNAL S KERNAL, 0/ENTRY C C C GET NSQUARE VALUES DO 320 I=1,NSQUARE S CLA CMA S TAD \I S TAD PKP /KP(1) S DCA 30 S TAD I 30 S SWAB /DOES MQL S CPAGE 4 S DAD S \ISUM S DST S \ISUM S CLAMQ 320 CONTINUE C C SCALE C COMPUTE: IZ=ISUM[1:2]/NSQUARE S CPAGE 4 S DAD S \ISUM S DVI S \NSQUARE S CLA S MQA S DCA \IZ S SWBA C C RETURN S JMP RKERNAL C************** P A R A M E T E R S ************* S BMTEXT, TEXT /BM/ S PKP, \KP /POINTER END