Как разделить информацию данных (ГГГММДД) ino ГГГГ ММ ДД из файла NetCDF в python? - PullRequest
Новая
       18

Как разделить информацию данных (ГГГММДД) ino ГГГГ ММ ДД из файла NetCDF в python?

1 голос
/ 21 июня 2019

Я довольно новичок в python, и я прочитал несколько предыдущих связанных вопросов, но они не совсем точно отвечают на мой вопрос:

У меня уже есть нарезанный файл .nc для осадков для моего примера.В Python я выполняю обработку, и мне нужно время в 3 разных столбцах [2015 01 31], хотя мне дается в одном столбце ГГГГММДД [20150101 20150131 .... 20171221]

Я пробовал использоватьмодуль netCDF4 и numpy для обработки и разделения данных по мере необходимости.

#Script for reading .nc files using the netcdf4 module in Python

import numpy as np
import netCDF4 as nc


os. chdir('/mnt/lustre02/work/ch0636/g260150/sims/validation/selectedmodelRCPs/RCP2.6/Model1/')

#extract time in a single array to append it to the whole dly. weather file

data=nc.Dataset('precipitation.nc','r')
time=data.variables['time'][:]
b=np.zeros((time.size,1))
b[:,0]=time

#Extracting time sections: year, month, day

year=[]
month=[]
day=[]

for i in range(len(b)):
    year.append(b[i][0:4])
    month.append(b[i][4:6])
    day.append(b[i][6:])

print(day)  
print(month)
print(year)

Запуск этой части дает тот день, который является целой датой, включая год, месяц и день.Месяц и день - пустые массивы.Я не совсем уверен, как я мог бы на самом деле подходить к этому должным образом, это помогает мне получить то, что мне нужно для дальнейшей обработки.

1 Ответ

0 голосов
/ 26 июня 2019

Чтобы получить массивы года, месяца и дня, используя ваше решение, я предлагаю модификацию следующим образом: 

year=[]
month=[]
day=[]

for i in range(len(b)):
    yval=np.floor(b[i]/10000.)
    mval=np.floor((b[i]-yval*10000.)/100.)
    dval=np.floor(b[i]-yval*10000.-mval*100.)    
    year.append(int(yval))
    month.append(int(mval))
    day.append(int(dval))

print(day)  
print(month)
print(year)

Тем не менее, я также могу предложить такой код, где мы сначала преобразуем время взначения даты, а затем извлечь годы, месяцы и дни: 

# other way:
import datetime
def get_date(datevalin):
    fraction_of_day=datevalin-np.floor(datevalin)
    d0=datetime.datetime.strptime(str(int(np.floor(datevalin))),'%Y%m%d')
    dval=d0+datetime.timedelta(seconds=fraction_of_day*24.*3600)
    return dval

dates=[get_date(val) for val in time];
year=[val.year for val in dates];
month=[val.month for val in dates];
day=[val.day for val in dates];
...