Me gustaría proceso de algunos netCDF 4 archivos con arcpy (ArcGIS 10.3.1). El original proceden los datos de NCEP/NCAR proyecto y contiene vectores de viento (uwnd/vwnd) para el world-wide puntos (lon/lat) en ciertas marcas de tiempo. La dimensión "tiempo" se da como tiempo relativo en
hours since 1800-01-01 00:00:0.0
Si tengo una mirada a los datos con CDO (Clima de Datos de los Operadores), los primeros cuatro sellos de tiempo de la uwnd.10m.de gauss.2016.nc son los siguientes:
cdo sinfo uwnd.10m.gauss.2016.nc
2016-01-01 00:00:00
2016-01-01 06:00:00
2016-01-01 12:00:00
2016-01-01 18:00:00
...
Por desgracia arcpy tiene un problema de comprensión de la relación sello de tiempo y el siguiente código
inNetCDF = r"\C:\temp\uwnd.10m.gauss.2016.nc"
ncFP = arcpy.NetCDFFileProperties(inNetCDF)
ncDim = ncFP.getDimensions()
for dim in ncDim:
print "%s (%s)" % (dim, ncFP.getFieldType(dim))
top = ncFP.getDimensionSize(dim)
for i in range(0,top):
print ncFP.getDimensionValue(dim,i)
ofrece esto como salida para la dimensión "tiempo":
07.03.1800 17:50:24
07.03.1800 17:50:42
07.03.1800 17:51:00
07.03.1800 17:51:18
...
Así que en lugar de una marca de tiempo cada 6 horas ArcGIS calcula una marca de tiempo cada 18 segundos?!
No es un gran problema para volver a calcular la fecha correcta dentro de python, pero me preguntaba si hay un parámetro que puedo configurar para que arcpy ofrece directamente la fecha y la hora correctas. He tenido el mismo problema con el mal de la fecha y hora utilizando el MakeNetCDFFeatureLayer caja de herramientas de la herramienta directamente dentro de ArcMap, y no puedo hacer nada de python cálculos.
