51 52 54 55 56 59 60 62南南南649 1083 9979南13602年5743 43444 6365 1283 9062 4026南南

cellsOfInt =数据(0:)

我为一些细胞nan值以来没有值的时间步骤。

每一行对应一个时间步长对感兴趣的细胞51,在第一次一步nan值,在第二次步骤中我们有5743等等,直到87648年。

所有的细胞都感兴趣的国家。txt文件与感兴趣的细胞950 X 1000 id或南当没有价值。

(国家id)……5152545556596062年…[时间步1]…南649年南南1083 9979南13602年…[时间步2]…5743 1283 9062 4026南6365 43444南…(时间步n)…………………………

cellsOfInt =数据(0:)

• 我为一些细胞nan值以来没有值的时间步骤。
• 每一行对应一个时间步长对感兴趣的细胞51,在第一次一步nan值,在第二次步骤中我们有5743等等,直到87648年。
• 所有的细胞都感兴趣的国家。txt文件与感兴趣的细胞950 X 1000 id或南当没有价值。

我的代码:

值= loadtxt (“countries.txt skiprows = 0) %的国家

data = loadtxt (energy.txt skiprows = 0) %能源

能量三种文件是:

51 52 54 55 56 59 60 62

649年南南南1083 9979南13602

5743 43444 6365 1283 9062 4026南南

等(直到87648步)

感兴趣的细胞在哪里第一行:

所有的细胞都感兴趣的国家。txt文件与感兴趣的细胞950 X 1000 id或南当没有价值。

y =数控。createDimension (" y ", 950) # y x =数控。createDimension (“x”, 1000) # x时间=数控。createDimension(“时间”,87648)#

纬度=数控。createVariable (“y”、“f8”(“y”))经度=数控。createVariable (“x”、“f8”(“x”))时间= nc.createVariable(“时间”、“f8”(“时间”))国家=数控。createVariable(“国家”、“f8”(" y "、" x "));能量=数控。createVariable(“能源”、“f8”(“时间”、“y”、“x”), fill_value = -9999, chunksize = (950, 1000));

时间。standard_name = '时间'时间。单位= '小时自2006-01-01 00:00:00.0 ' time.calendar =“proleptic_gregorian”

背阔肌= np.arange (5497500、747500、-5000) # y朗= np.arange (2502500, 2502500, 5000) # x

#填充变量[:]=背阔肌经度纬度[:]=经度

国家[:]=价值

xrange(0950):我的#行j xrange(0, 1000): #列在xrange n(0, 87648): #行米xrange(0,1): #如果数据列(1米)= =国家(i, j): #函件nc.variables数据(“能源”)[n] = [n m]

值= loadtxt (“countries.txt skiprows = 0) % = loadtxt国家数据(“energy.txt”, skiprows = 0) %的能量

的能源. txt文件是:

51 52 54 55 56 59 60 62 649年南南南1083 9979南13602 5743 43444 6365 1283 9062 4026南南

等(直到87648步),感兴趣的细胞在哪里吗的第一行:

所有的细胞都感兴趣的国家。txt文件与感兴趣的细胞950 X 1000 id或南当没有价值。

y =数控。createDimension (" y ", 950) # y x =数控。createDimension (“x”, 1000) # x时间=数控。createDimension(“时间”,87648)#纬度=数控。createVariable (“y”、“f8”(“y”))经度=数控。createVariable (“x”、“f8”(“x”))时间= nc.createVariable(“时间”、“f8”(“时间”))国家=数控。createVariable(“国家”、“f8”(" y "、" x "));能量=数控。createVariable(“能源”、“f8”(“时间”、“y”、“x”), fill_value = -9999, chunksize = (950, 1000));时间。standard_name = '时间'时间。units ='hours since 2006-01-01 00:00:00.0' time.calendar='proleptic_gregorian' lats = np.arange(5497500,747500,-5000) #y lons = np.arange(2502500,7502500,5000) #x # Fill variables latitude[:] = lats longitude[:] = lons countries [:]=value

xrange(0950):我的#行j xrange(0, 1000): #列在xrange n(0, 87648): #行米xrange(0,1): #如果数据列(1米)= =国家(i, j): #函件nc.variables数据(“能源”)[n] = [n m]