[R-es] Unir archivos netCDF
Carlos Ortega
cof en qualityexcellence.es
Lun Ene 9 15:02:53 CET 2017
Hola,
Por si te puede ayudar:
http://stackoverflow.com/questions/29101807/how-to-merge-netcdf-file-into-one
http://stackoverflow.com/questions/18866266/combine-netcdf-files-to-average-the-values-for-variables-in-r
Saludos,
Carlos Ortega
www.qualityexcellence.es
2017-01-09 10:34 GMT+01:00 Isa García Barón <xanamb en gmail.com>:
> Buenos días lista,
>
> Lo que quiero es unir archivos netCDF. Tengo varios archivos con las mismas
> características, pero de días consecutivos y me gustaría unir todos.
>
> He probado con las funciones "transNcdfMerge" y "aggregateNcdf", ambas del
> paquete ncdf.tools. Con ninguna de las dos he conseguido nada, a
> continuación os copio el script que he escrito hasta el momento:
>
> setwd(DataDir_geo_cur)
>
> ## lista con todos los archivo netCDF contenido en el directorio:
> files = list.files(DataDir_geo_cur,pattern='*.nc',full.names=TRUE)
>
> ## leemos todos los archivos:
> for(i in seq_along(files)) {
> assign(paste0("gc",i),nc_open(files[i]))
> }
>
> # ejemplo de uno de los archivos
> gc1
>
> File C:/use/dt_global_allsat_msla_uv_20070905_20140106.nc
> (NC_FORMAT_CLASSIC):
>
> 5 variables (excluding dimension variables):
> float lat_bnds[nv,lat]
> float lon_bnds[nv,lon]
> int crs[]
> grid_mapping_name: latitude_longitude
> semi_major_axis: 6371000
> inverse_flattening: 0
> int u[lon,lat,time]
> _FillValue: -2147483647
> long_name: Geostrophic velocity anomalies: zonal component
> standard_name:
> surface_eastward_geostrophic_sea_water_velocity_assuming_
> sea_level_for_geoid
> units: m/s
> scale_factor: 1e-04
> int v[lon,lat,time]
> _FillValue: -2147483647
> long_name: Geostrophic velocity anomalies: meridian component
> standard_name:
> surface_northward_geostrophic_sea_water_velocity_assuming_
> sea_level_for_geoid
> units: m/s
> scale_factor: 1e-04
>
> 4 dimensions:
> time Size:1
> long_name: Time
> standard_name: time
> units: days since 1950-01-01 00:00:00 UTC
> calendar: julian
> axis: T
> lat Size:720
> long_name: Latitude
> standard_name: latitude
> units: degrees_north
> bounds: lat_bnds
> axis: Y
> valid_min: -90
> valid_max: 90
> lon Size:1440
> long_name: Longitude
> standard_name: longitude
> units: degrees_east
> bounds: lon_bnds
> axis: X
> valid_min: 0
> valid_max: 360
> nv Size:2
>
> 26 global attributes:
> cdm_data_type: Grid
> title: DT merged Global Ocean Gridded Geostrophic Velocities
> SSALTO/Duacs L4 product
> summary: This dataset contains Delayed Time Level-4 geostrophic
> velocities products from multi-satellite observations over Global Ocean.
> comment: Surface product; Geostrophic Velocities referenced to the
> [1993, 2012] period
> time_coverage_resolution: P1D
> product_version: 5.0
> institution: CNES, CLS
> project: SSALTO/DUACS
> references: www.aviso.altimetry.fr
> contact: aviso en altimetry.fr
> license:
> http://www.aviso.altimetry.fr/fileadmin/documents/data/License_Aviso.pdf
> platform: Envisat, Geosat Follow On, Jason-1
> date_created: 2014-02-27 20:27:23
> history: 2014-02-27 20:27:23:creation
> Conventions: CF-1.6
> standard_name_vocabulary:
> http://cf-pcmdi.llnl.gov/documents/cf-standard-names/
> standard-name-table/12/cf-standard-name-table.html
> geospatial_lat_min: -90
> geospatial_lat_max: 90
> geospatial_lon_min: 0
> geospatial_lon_max: 360
> geospatial_vertical_min: 0.0
> geospatial_vertical_max: 0.0
> geospatial_lat_units: degrees_north
> geospatial_lon_units: degrees_east
> geospatial_lat_resolution: 0.25
> geospatial_lon_resolution: 0.25
>
> ## creamos un vector con todos los archivos abiertos en el environment:
> names = ls(pattern='gc')
>
> ## Intento de unión con la función transNcdfMerge:
> Geo_cur <- transNcdfMerge(names, path.target = DataDir)
>
> Error in convertDateNcdf2R(fileT) : File gc1 is not existent!
>
> ## Intento de unión con la función aggregateNcdf:
> Geo_cur_1 <- aggregateNcdf(names, path.out = DataDir, period=1)
>
> Error in system(command, as.integer(flag), f, stdout, stderr) : character
> string expected as first argument
>
> Espero que podaís ayudarme, muchísimas gracias y feliz año!
>
> *-----------------------------------------------------------
> ------------------------------------------*
> *Isabel García Barón*
> Email: xanamb en gmail.com
> PhD Student at AZTI Foundation - AZTI Fundazioa
> Marine Ecosystems Functioning
> Herrera Kaia, Portualdea z/g 20110 - Pasaia (Gipuzkoa)
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es