[R] Min , Max

[Jonsson]

Hellow everyone,
This code bellow will calculate average daily wind speed(measurements are
taken every three hours).Any ideas how to take the Min and Max instead of
average.

library(Matrix)
setwd("C:\\Users\\aalyaari\\Desktop\\img")
listfile<-dir()
long <- file("C:\\Users\\aalyaari\\Desktop\\New folder (5)\\inra.bin", "rb")
 A=readBin(long, integer(), size=2,n=67420*1, signed=F)
 ta<-t(A)
 lot <- file("C:\\Users\\aalyaari\\Desktop\\New folder (5)\\lat.img", "rb")
 B=readBin(lot, integer(), size=2,n=67420*1, signed=F)
 tb<-t(B)

for (n in 1:length(listfile))
{
    
	#h[n]=listfile[n]
	h=listfile[n]
    	#b[n]=file.info(h[n])$size/67420/4
	b=file.info(h[n])$size/67420/4

 	wind <- file(h, "rb")
 	C=readBin(wind, double(), size=4,n=67420*b, signed=TRUE)
 
 	D<-matrix(C,nrow=b,ncol=67420)
 
	for(d in 1:b)
 	{
	 	M <- Matrix(-9999, 360, 720)
       	tm<-t(M)
	 	for(i in 1:67420)
	 	{
			tm[ta[i],tb[i]]= round(10 * ((D[(d-1)*8+1,i] + D[(d-1)*8+2,i]
+D[(d-1)*8+3,i] +D[(d-1)*8+4,i] +D[(d-1)*8+5,i] +D[(d-1)*8+6,i]
+D[(d-1)*8+7,i] +D[(d-1)*8+8,i] ) / 8))
	
	 	}
	 	to.write <- sprintf("C:\\Users\\aalyaari\\Desktop\\New folder
(6)\\Yar_%00d.bin", d)
	 	writeBin(as.integer(tm at x), size=2,to.write)
 	}
}


--
View this message in context: http://r.789695.n4.nabble.com/Min-Max-tp4593065p4593065.html
Sent from the R help mailing list archive at Nabble.com.