[R-es] Unir archivos netCDF

[Carlos Ortega]

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]]