Pythonを用いた空間統計分析のための衛星データ（LST, NDVI）の処理1

#!/usr/bin/env python
# coding: utf-8

import h5py
import pandas as pd
import tables
from pyhdf.SD import SD, SDC
import numpy as np
import pandas as pd
from scipy.interpolate import griddata
import seaborn as sns
import matplotlib.pyplot as plt
import os
import glob
import shutil


## h27v04といった場所に応じたフォルダ毎にファイルを分ける処理

#### NDVIの処理
new_dir_pathNDVI_h27v04 = './Data/
NDVI_Before_change_of_coordinate/NDVI_h27v04'
new_dir_pathNDVI_h28v04 = './Data/
NDVI_Before_change_of_coordinate/NDVI_h28v04'
new_dir_pathNDVI_h28v05 = './Data/
NDVI_Before_change_of_coordinate/NDVI_h28v05'
new_dir_pathNDVI_h29v05 = './Data/
NDVI_Before_change_of_coordinate/NDVI_h29v05'
os.makedirs(new_dir_pathNDVI_h27v04, exist_ok=True)
os.makedirs(new_dir_pathNDVI_h28v04, exist_ok=True)
os.makedirs(new_dir_pathNDVI_h28v05, exist_ok=True)
os.makedirs(new_dir_pathNDVI_h29v05, exist_ok=True)

area_folders = [new_dir_pathNDVI_h27v04, new_dir_pathNDVI_h28v04, 
new_dir_pathNDVI_h28v05, new_dir_pathNDVI_h29v05]

### 各領域のファイル名を抽出する
files_h27v04 = glob.glob('./Data/NDVI_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h27v04.csv')
files_h28v04 = glob.glob('./Data/NDVI_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h28v04.csv')
files_h28v05 = glob.glob('./Data/NDVI_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h28v05.csv')
files_h29v05 = glob.glob('./Data/NDVI_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h29v05.csv')

area_files = [files_h27v04, files_h28v04, files_h28v05, files_h29v05]

# ファイルを移動させる
for i in range(len(area_files)):
    
    for j in range(len(area_files[i])):
        new_path = shutil.move(area_files[i][j], area_folders[i])

#### LSTの処理
new_dir_pathLST_h27v04 = './Data/
LST_Before_change_of_coordinate/LST_h27v04'
new_dir_pathLST_h28v04 = './Data/
LST_Before_change_of_coordinate/LST_h28v04'
new_dir_pathLST_h28v05 = './Data/
LST_Before_change_of_coordinate/LST_h28v05'
new_dir_pathLST_h29v05 = './Data/
LST_Before_change_of_coordinate/LST_h29v05'
os.makedirs(new_dir_pathLST_h27v04, exist_ok=True)
os.makedirs(new_dir_pathLST_h28v04, exist_ok=True)
os.makedirs(new_dir_pathLST_h28v05, exist_ok=True)
os.makedirs(new_dir_pathLST_h29v05, exist_ok=True)

area_folders = [new_dir_pathLST_h27v04, new_dir_pathLST_h28v04, 
new_dir_pathLST_h28v05, new_dir_pathLST_h29v05]

### 各領域のファイル名を抽出する
files_h27v04 = glob.glob('./Data/LST_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h27v04.csv')
files_h28v04 = glob.glob('./Data/LST_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h28v04.csv')
files_h28v05 = glob.glob('./Data/LST_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h28v05.csv')
files_h29v05 = glob.glob('./Data/LST_Before_change_of_coordinate/
A[0-9][0-9][0-9][0-9][0-9][0-9][0-9].h29v05.csv')

area_files = [files_h27v04, files_h28v04, files_h28v05, files_h29v05]

# ファイルを移動させる
for i in range(len(area_files)):
    for j in range(len(area_files[i])):
        new_path = shutil.move(area_files[i][j], area_folders[i])