R: 東京電力の電力使用状況グラフ (2011年7月版)

R: 東京電力の電力使用状況グラフをえがくRコードを更新。



コード

## TEPCO site
csvurl <- "http://www.tepco.co.jp/forecast/html/images/juyo-j.csv"

## データ開始行-1
l60 <- 8
l05 <- 44

## peak
con <- url(csvurl, encoding = "Shift_JIS")
lines3 <- readLines(con, 3)
maxpow <- as.numeric(strsplit(lines3[3], ",")[[1]][1])
maxpowctime <- strsplit(lines3[3], ",")[[1]][2]
maxpowctime <- strsplit(maxpowctime, "〜")[[1]]

## 1時間ごとのデータ
if (FALSE) {
pow <- read.csv(con, skip = l60, nrows = 24, header = FALSE,
                col.names = c("cdate", "ctime",
                              "today", "predict"))
pow$date <- strptime(paste(pow$cdate, pow$ctime),
                     "%Y/%m/%d %H:%M", tz = "Asia/Tokyo")
pow$today <- ifelse(pow$today == 0, NA, pow$today)
}

## 5分ごとのデータ
pow5 <- read.csv(con, skip = l05, nrows = 24*12, header = FALSE,
                col.names = c("cdate", "ctime",
                              "today"))
pow5$date <- strptime(paste(pow5$cdate, pow5$ctime),
                     "%Y/%m/%d %H:%M", tz = "Asia/Tokyo")
pow5$today <- ifelse(pow5$today == 0, NA, pow5$today)

## 予想最大電力時間帯
maxpowtime1 <- strptime(paste(pow5$cdate[1], maxpowctime[1]),
                       "%Y/%m/%d %H:%M", tz = "Asia/Tokyo")
maxpowtime2 <- strptime(paste(pow5$cdate[1], maxpowctime[2]),
                       "%Y/%m/%d %H:%M", tz = "Asia/Tokyo")

## plot
colors <- c(3, 1, 2)
par(mar = c(5, 5, 4, 1), mgp = c(4, 1, 0))
plot(pow5$date, rep(0, length(pow5$date)), type = "n",
     ylim = c(0, floor(maxpow / 500 + 1) * 500),
     xlab = "時", ylab = "消費電力(万kW)",
     main = pow5$cdate[1], las = 1, xaxt = "n")
axis.POSIXct(1, at = pow5$date[seq(1, 24*12, 3*12)])

lines(c(maxpowtime1, maxpowtime2),
      rep(maxpow, 2),
      type = "l", col = colors[1])
lines(pow5$date, pow5$today, col = colors[2])
#lines(pow$date, pow$predict, col = colors[3])
legend("bottomright",
       c("ピーク時供給力", "当日実績"),
       col = colors, lty = 1)
#legend("bottomright",
#       c("ピーク時供給力", "当日実績", "予測値"),
#       col = colors, lty = 1)

タグ：R