※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

#hdpで
  • a:全てのファイルを、-d:データだけ、-b:バイナリ -x:テキスト
#hdp dumpsds -a -d -b A20132132013243.L3m_MO_CHL_chlor_a_9km.hdf > data.bin
#hdp dumpsds -a -h A20132132013243.L3m_MO_CHL_chlor_a_9km.hdf > read.txt
hdp dumpsds -a -d -x A20132132013243.L3m_MO_CHL_chlor_a_9km.hdf > data.txt

#R
テキストファイルでやる
d <- scan("data.txt", ncol=4320)
dm <- matrix(d, nrow=2160, ncol=4320, byrow=TRUE)
dm4 <- dm[517:600, 3805:3936] #北海道周辺の切り出しN40-47 E137-148
dm5 <- replace(dm4,which(dm4==-32767),NA) #-32767を0に置換

#1度グリッド #抽出
for loop 横に進んで一段下がりまた横に進む
dt <- c()
for(j in seq(from=1,to=84,by=12)){
	for(i in seq(from=1,to=132,by=12)){
		dx <- dm5[j:(j+11),i:(i+11)]
		dt <- c(dt,mean(dx, na.rm=TRUE))
		}
	}

#1
/2度グリッド
#抽出
for loop 横に進んで一段下がりまた横に進む
dt <- c()
for(j in seq(from=1,to=84,by=6)){
	for(i in seq(from=1,to=132,by=6)){
		dx <- dm5[j:(j+5),i:(i+5)]
		dt <- c(dt,mean(dx, na.rm=TRUE))
		}
	}
|