[R] Min , Max

Jonsson amen.alyaari at Bordeaux.inra.fr
Fri Apr 27 17:46:59 CEST 2012


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.



More information about the R-help mailing list