東京の地価変動グラフ

前回に続き、東京都内の地価公示データを時系列のグラフにしてみました。
1984年から2011年までの27年間で7万1千件の地価のデータがあります。
二つの灰色の点があるところが日経平均株価の最高値を記録した1989年12月29日と
リーマン・ブラザーズが連邦倒産法を申請した2008年9月15日です。





Land price change in Tokyo Flickr Original (2942 x 1674)

今回の作り方

  • 国土交通省のXMLデータをリバースジオコーディングして市町村ごとにまとめる
  • 任意の市町村ごとに1平方メートルあたりの平均価格を求める
  • 時系列のcsvテキストにまとめたものを、javascriptで読み込みillustratorにチャートのパスとして描画
  • カラー調整など仕上げ
  • resouce
    国土数値情報 地価公示データ

    Posted: January 30th, 2012 | Author: | Filed under: 技術 | Tags: , , , | No Comments »

    首都圏の地価の変動を視覚化してみた

    国土交通省が公開している地価公示データを使って、首都圏の土地の値段が
    前年と比べてどのように変動しているかを出してみました。

    2008年は前年と比べて、特に東京周辺が上昇していましたが、
    2009, 2010年は下落が目立ちます。これがいわゆるリーマン・ショックなんでしょうか。
    2012年のデータが公開されれば、地震の影響が見れるかも知れません。




    Visualizing land prices around Tokyo (Tokyo, Saitama, Chiba and Kanagawa).
    This represents the impact of Lehman ??

    resouce
    国土数値情報 地価公示データ

    Posted: January 22nd, 2012 | Author: | Filed under: 技術 | Tags: , , , , | No Comments »

    rpy2 Scatter plot

    rpy2で散布図を書く

    import rpy2.robjects as robjects
    from rpy2.robjects.packages import importr
    
    grdevices = importr('grDevices')
    grdevices.png(file="file.png", width=512, height=512)
    # plotting code here
    
    # 月齢
    lstx = [0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0,
            10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0,
            20.0, 21.0, 22.0, 23.0, 24.0, 25.0, 26.0, 27.0, 28.0, 29.0]
    
    # 漁獲高
    lsty = [190.0, 200.0, 50.0, 1700.0, 360.0, 620.0, 450.0, 2100.0,
            120.0, 5000.0, 1900.0, 780.0, 880.0, 1500.0, 1500.0, 1900.0,
            200.0, 270.0, 2100.0, 7000.0, 900.0, 1200.0, 3700.0, 2000.0,
            2900.0, 1300.0, 140.0, 2250.0, 120.0, 1000.0]
    
    rx = robjects.FloatVector(lstx)
    ry = robjects.FloatVector(lsty)
    
    robjects.r.plot(x=rx, y=ry, xlab="moon age", ylab="fish", col="blue")
    
    grdevices.dev_off()
    

    file
    csvファイルから読み込む場合は、Rの命令文をそのまま書いてもよい

    import rpy2.robjects as robjects
    from rpy2.robjects.packages import importr
    grdevices = importr('grDevices')
    grdevices.png(file="file.png", width=512, height=512)
    
    robjects.r('''
    data <-read.csv('fish_moonage.csv', header=TRUE)
    plot(fish~moonage, xlab="foo", data=data, col="purple")
    ''')
    
    grdevices.dev_off()
    
    Posted: August 27th, 2010 | Author: | Filed under: 技術 | Tags: , , , , | No Comments »