gsmapcsv2csv

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
# GSMaPのCSVファイルから特定の地点の日降水量のグラフを作成する
# データはJAXAのftpで公開されている地域ごとの降水量csvファイル(日データ)
# ftp://hokusai.eorc.jaxa.jp/

# 必要なモジュールのインポート
import pandas as pd
import glob
import datetime

# 定数の設定
# 対象地の緯度経度及びそれがが記録されている行
# 緯度経度から行を抽出したいが、とりあえずマニュアルで記載。要改善
terget_line = 10729

# データの読み込み
files = glob.glob('../data/*.csv')

dates_int = []
for x in files:
    dates_int.append(int(x[26:34]))

df_rain = pd.DataFrame()
for x in files:
    df_rain = df_rain.append(pd.read_csv(x)[terget_line-2:terget_line-1])

df_rain['date'] = dates_int
df_rain = df_rain.sort_values(by="date")
df_rain.to_csv("../out/rain_Kaolack.csv", sep=",")