#make AVHRR svd #dated Mar 22,2009: http://faculty.washington.edu/steig/nature09data/data/cloudmaskedAVHRR.txt #unavailable MArr 22, 2009 #300 rows monthly from Jan 1982 to end 2006; 5509 columns download.file("http://faculty.washington.edu/steig/nature09data/data/cloudmaskedAVHRR.txt","temp.dat",mode="wb") grid=scan("temp.dat",n= -1) # 37800 length(grid)/5509 #300 avhrr=ts(t(array(grid,dim=c(5509,300))),start=c(1982,1),freq=12) tsp(avhrr) #[1] 1982.000 2006.917 12.000 #300 5509 avhrr=array(avhrr,dim=c(12,nrow(avhrr)/12, 5509) );dim(avhrr) avhrr=aperm(avhrr,perm=c(2,1,3)); dim(avhrr) #25 12 5509 avhrr= array(avhrr,dim=c(25, 12*5509) ) avhrr=scale(avhrr,scale=FALSE) #25 66108 avhrr= array(avhrr, dim=c( 25,12,5509) ) avhrr= aperm(avhrr, perm=c (2,1,3) ) ; dim(avhrr) #12 25 5509 avhrr=array(avhrr, dim=c( 12*25,5509) ) save(avhrr,file="d:/climate/data/steig/avhrr_anom.tab") load("d:/climate/data/steig/avhrr_anom.tab");dim(avhrr) # 300 5509 #collated in collation.steig.txt from steig/nature09data/data/cloudmaskedAVHRR.txt W=svd(t( scale(avhrr,scale=FALSE) ),nu=3,nv=3) svd.avhrr_anom=W save(svd.avhrr_anom,file="d:/climate/data/steig/svd.avhrr_anom.tab") # uploaded to www.climateaudit.org/data/steig/svd.avhrr_anom.tab U=W$v;dim(U) #300 3 V=W$u; dim(V)