An Evaluation of the Current State of Digital Photography

Charles Dickinson


Appendix A

PRO dc260
 

filename=PICKFILE(title='Select the TIFF image ',/R,filter='/nfs/cis/ugrad/chd4750/Senior/luminance260/*.tif')
image1=READ_TIFF(filename)
;tv, image1
xsize=N_ELEMENTS(image1(*,1))
ysize=N_ELEMENTS(image1(1,*))
mean = FLTARR(11)
;-----------------------------------------------------------------
;each of the following sections looks within the defined area (ex. ref1), takes the mean of the digital values, and also prints the standard deviation to the screen

;FOR DENSITY = .1
ref1=image1[665:685,500:525]
sum1=TOTAL(ref1)
mean[0]=sum1/N_ELEMENTS(ref1)
mom1=MOMENT(ref1,SDEV=s1)
print,s1
;------------------------------------------------------------------
 

;FOR DENSITY = .2
ref2=image1[700:712,500:525]
sum2=TOTAL(ref2)
mean[1]=sum2/N_ELEMENTS(ref2)
mom2=MOMENT(ref2, SDEV=s2)
print,s2
;------------------------------------------------------------------
 

;FOR DENSITY = .3
ref3=image1[722:735,500:525]
sum3=TOTAL(ref3)
mean[2]=sum3/N_ELEMENTS(ref3)
mom3=MOMENT(ref3, SDEV=s3)
print,s3
;------------------------------------------------------------------
 

;FOR DENSITY = .4
ref4=image1[746:757,500:525]
sum4=TOTAL(ref4)
mean[3]=sum4/N_ELEMENTS(ref4)
mom4=MOMENT(ref4,SDEV=s4)
print,s4
;------------------------------------------------------------------
 

;FOR DENSITY = .5
ref5=image1[767:780,500:525]
sum5=TOTAL(ref5)
mean[4]=sum5/N_ELEMENTS(ref5)
mom5=MOMENT(ref5,SDEV=s5)
print,s5
;------------------------------------------------------------------
 

;FOR DENSITY = .6
ref6=image1[792:803,500:525]
sum6=TOTAL(ref6)
mean[5]=sum6/N_ELEMENTS(ref6)
mom6=MOMENT(ref6,SDEV=s6)
print,s6
;------------------------------------------------------------------

;FOR DENSITY = .7
ref7=image1[815:825,500:525]
sum7=TOTAL(ref7)
mean[6]=sum7/N_ELEMENTS(ref7)
mom7=MOMENT(ref7,SDEV=s7)
print,s7
;------------------------------------------------------------------
 

;FOR DENSITY = .8
ref8=image1[836:850,500:525]
sum8=TOTAL(ref8)
mean[7]=sum8/N_ELEMENTS(ref8)
mom8=MOMENT(ref8,SDEV=s8)
print,s8
;------------------------------------------------------------------
 

;FOR DENSITY = .9
ref9=image1[860:875,500:525]
sum9=TOTAL(ref9)
mean[8]=sum9/N_ELEMENTS(ref9)
mom9=MOMENT(ref9,SDEV=s9)
print,s9
;------------------------------------------------------------------
 

;FOR DENSITY = 1.0
ref10=image1[883:895,500:525]
sum10=TOTAL(ref10)
mean[9]=sum10/N_ELEMENTS(ref10)
mom10=MOMENT(ref10,SDEV=s10)
print,s10
;------------------------------------------------------------------
 

;FOR DENSITY = 1.1
ref11=image1[905:930,500:525]
sum11=TOTAL(ref11)
mean[10]=sum11/N_ELEMENTS(ref11)
mom11=MOMENT(ref11,SDEV=s11)
print,s11
;------------------------------------------------------------------

;Here, the variables where the average digital counts are stored get read out to a text file

filename2=PICKFILE(title='Select the OECF file to write ',/W,filter= ' *.txt')
OPENW,1, filename2
PRINTF,1,mean[0]
PRINTF,1,mean[1]
PRINTF,1,mean[2]
PRINTF,1,mean[3]
PRINTF,1,mean[4]
PRINTF,1,mean[5]
PRINTF,1,mean[6]
PRINTF,1,mean[7]
PRINTF,1,mean[8]
PRINTF,1,mean[9]
PRINTF,1,mean[10]
CLOSE,1
END


Table of Contents | Thesis