netcdf4 изменяет переменные на основе других меньших значений netcdf с разными широтами и долями - PullRequest
Новая
       10

netcdf4 изменяет переменные на основе других меньших значений netcdf с разными широтами и долями

0 голосов
/ 31 августа 2018

Я пытаюсь создать пустой файл netCDF и использовать другие меньшие файлы netCDF для заполнения пустых значений netCDF. Это дает мне ошибку индекса:

логический индекс не соответствует индексированному массиву по измерению 0; измерение равно 100001, но соответствующее логическое измерение равно 3214 в sourcef [(d_AS <= 30)] = d_AS </p>

Весь код выглядит следующим образом: 

import glob

import numpy as np
from netCDF4 import Dataset

outf=Dataset('results.nc','w',format='NETCDF4_CLASSIC')
lats=np.linspace(20.0000,50.0000,100001)
lons=np.linspace(-130.0000,-100.0000,100001)


lat=outf.createDimension('lon',len(lats))
lon=outf.createDimension('lat',len(lons))
longitude=outf.createVariable('longitude',np.float64,('lon',))
latitude=outf.createVariable('latitude',np.float64,('lat',))
WFHS=outf.createVariable('AREASCORE',np.int8,('lon','lat'))

longitude[:]=lons
latitude[:]=lats
AREASCORE=np.full((100001,100001),255,dtype=np.uint8)

WFHS[:,:]=AREASCORE

outf.description="AREASCORE for 13 states"
longitude.units="degrees east"
latitude.units='degrees north'

print ("created empty array")
sourcef=outf.variables['AREASCORE'][:]

for d in glob.glob(*.nc'):
    inf=Dataset(d)

    d_AS=inf.variables['AREASCORE'][:]
    d_AS=d_AS.astype(int)
    d_AS[np.isnan(d_AS)]=255
    sourcef[(d_AS<=30)]=d_AS
    outf.variables['AREASCORE'][:]=sourcef
    d.close()



outf.close()

Это дает мне ошибку

...