r3.out.netcdfgrass – En ligne dans le cloud

PROGRAMME:

Nom

r3.out.netcdf - Exporter une carte raster 3D sous forme de fichier netCDF.

MOTS-CLÉS

raster3d, exportation, netCDF, voxel

SYNOPSIS

r3.out.netcdf
r3.out.netcdf --Aidez-moi
r3.out.netcdf [-pm] contribution=prénom sortie=prénom [nul=flotter] [--écraser] [--vous aider]
[--verbeux] [--calme] [--ui]

Drapeaux:
-p
Exporter les informations de projection sous forme de paramètres wkt et proj4

-m
Utiliser un masque raster 3D (s'il existe) avec la carte d'entrée

--écraser
Autoriser les fichiers de sortie à écraser les fichiers existants

--Aidez-moi
Imprimer le récapitulatif d'utilisation

--verbeux
Sortie du module verbeux

--silencieux
Sortie module silencieuse

--interface utilisateur
Forcer le lancement de la boîte de dialogue GUI

Paramètres:
contribution=prénom [obligatoire]
Nom de la carte raster 3D en entrée

sortie=prénom [obligatoire]
Nom du fichier de sortie netCDF

nul=flotter
La valeur à utiliser pour les valeurs nulles, la valeur par défaut est la norme netCDF

DESCRIPTION

Le module r3.out.netcdf Exporte une carte raster 3D au format NetCDF. Les cartes sont des rasters 3D valides.
cartes dans le chemin de recherche du jeu de cartes actuel. sortie le paramètre est le nom d'un fichier netCDF
qui sera écrit dans le répertoire de travail courant. Le module est sensible au mapset
paramètres de région (définis avec g.région) et les paramètres de l'unité verticale de la carte raster 3D.
L'unité verticale peut être définie avec r3.support.

Note

Le fichier netCDF résultant aura le temps comme troisième dimension dans le cas d'une verticale temporelle
l'unité est spécifiée pour la carte d'entrée, sinon la troisième dimension est supposée être spatiale.
Si la carte raster 3D possède un horodatage avec une heure absolue, l'heure de début sera utilisée comme
Point de départ de la série chronologique netCDF. Si l'horodatage absolu est manquant, la date
La valeur par défaut sera « 01 janv. 1900 00:00:00 ». En cas d'horodatage relatif, aucun début n'est indiqué.
L'heure est définie dans le fichier netCDF. Le nombre de profondeurs de la carte raster 3D correspond au nombre
des tranches de temps dans le fichier netCDF. Vous pouvez utiliser t.rast.à.rast3 convertir un espace-temps
un ensemble de données raster dans une carte raster 3D représentant un cube voxel espace-temps, puis exportez le
séries temporelles raster sous forme de fichier netCDF.

Les coordonnées spatiales sont exportées sous forme de coordonnées centrées sur la cellule. La projection peut être
éventuellement stockés dans les métadonnées sous la forme crs attributs. Les métadonnées de projection netCDF
le stockage suit la suggestion spatial_ref GDAL/netCDF ici et le netCDF CF 1.6
Convention utilisant ici les informations de projection WKT. Une chaîne PROJ.4 est également stockée dans
le crs section d'attributs. L'exportation des paramètres de projection est supprimée lorsque
La projection XY est définie.

La portée de la carte raster 3D est définie dans le fichier de sortie netCDF. Le fichier netCDF est optionnel.
la valeur manquante et _FillValue peuvent être définies à l'aide de nul option. Par défaut, le netCDF
Les valeurs à virgule flottante _FillValues ​​sont utilisées pour remplir une cellule vide.

La bibliothèque netCDF doit être installée sur le système et activée au moment de la configuration pour
compiler ce module.

Exemples

g.region s=-90 n=90 w=-180 e=180 b=0 t=5 res=10 res3=10 tbres=1 -p3
################################################ ###################
# Nous créons une carte de volume simple avec des valeurs à virgule flottante
################################################ ###################
r3.mapcalc --o expr="volume_float = float(col() + row() + depth())"
r3.info volume_float
+---------------------------------------------------------------- ---------------------------+
| Couche : volume_float Date : jeu. 14 juin 08:40:56 2012 |
| Mapset : PERMANENT Connexion du créateur : soeren |
| Emplacement : TestLL |
| Base de données : /1/soeren/grassdata |
| Titre : volume_float |
| Unités : aucune |
| Unité verticale : unités |
| Horodatage : aucun |
| ----------------------------------------------------- --------------------------- |
| |
| Type de carte : Cellule 3d Nombre de catégories : 0 |
| Type de données : FCELL |
| Rangées : 18 |
| Colonnes : 36 |
| Profondeurs : 5 |
| Nombre total de cellules : 3240 |
| Taille totale : 10114 octets |
| Nombre de tuiles : 1 |
| Taille moyenne des tuiles : 10114 octets |
| Taille de la tuile en mémoire : 12960 octets |
| Nombre de tuiles en x, y et z : 1, 1, 1 |
| Dimension d'une tuile en x, y, z : 36, 18, 5 |
| |
| Projection : Latitude-Longitude (zone 0) |
| N: 90N S: 90S Rés: 10 |
| E : 180E O : 180W Rés : 10 |
| T : 5 B : 0 Rés : 1 |
| Plage de données : min = 3 max = 59 |
| |
| Source de données : |
| |
| |
| |
| Description des données : |
| généré par r3.mapcalc |
| |
| Commentaires : |
| r3.mapcalc expression="volume_float = float(col() + row() + depth())" |
| |
+---------------------------------------------------------------- ---------------------------+
################################################ ###################
# Nous utilisons l'outil netCDF ncdump pour examiner l'en-tête
# et coordonnées du fichier netCDF exporté
################################################ ###################
r3.out.netcdf --o entrée=volume_float sortie=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
dimensions:
longitude = 36 ;
latitude = 18 ;
z = 5 ;
variables:
longitude flottante(longitude) ;
longitude:unités = "degrés_est" ;
longitude:long_name = "Valeurs de longitude" ;
longitude:nom_standard = "longitude" ;
longitude:axe = "X" ;
latitude flottante(latitude) ;
latitude:unités = "degrés_nord" ;
latitude:long_name = "Valeurs de latitude" ;
latitude:nom_standard = "latitude" ;
latitude:axe = "Y" ;
float z(z) ;
z:unités = "mètre" ;
z:long_name = "coordonnée z de la projection" ;
z:standard_name = "coordonnée_z_projection" ;
z:positif = "haut" ;
z:axe = "Z" ;
float volume_float(z, latitude, longitude) ;
volume_float:valid_min = 3.f ;
volume_float:valid_max = 59.f ;
// attributs globaux :
:Conventions = "CF-1.5" ;
:history = "Exportation netCDF GRASS GIS 7 de r3.out.netcdf" ;
données:
longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
-65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
105, 115, 125, 135, 145, 155, 165, 175 ;
latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
-65, -75, -85 ;
z = 0.5, 1.5, 2.5, 3.5, 4.5 ;
}
################################################ ###################
# Exportation des paramètres spécifiques à la projection et d'une valeur nulle
################################################ ###################
r3.out.netcdf --o -p null=-1 entrée=volume_float sortie=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
dimensions:
longitude = 36 ;
latitude = 18 ;
z = 5 ;
variables:
char crs ;
crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degré\",0.0174532925199433]]" ;
crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degré\",0.0174532925199433]]" ;
crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
longitude flottante(longitude) ;
longitude:unités = "degrés_est" ;
longitude:long_name = "Valeurs de longitude" ;
longitude:nom_standard = "longitude" ;
longitude:axe = "X" ;
latitude flottante(latitude) ;
latitude:unités = "degrés_nord" ;
latitude:long_name = "Valeurs de latitude" ;
latitude:nom_standard = "latitude" ;
latitude:axe = "Y" ;
float z(z) ;
z:unités = "mètre" ;
z:long_name = "coordonnée z de la projection" ;
z:standard_name = "coordonnée_z_projection" ;
z:positif = "haut" ;
z:axe = "Z" ;
float volume_float(z, latitude, longitude) ;
volume_float:valid_min = 3.f ;
volume_float:valid_max = 59.f ;
volume_float:missing_value = -1.f ;
volume_float:_FillValue = -1.f ;
volume_float:grid_mapping = "crs" ;
// attributs globaux :
:Conventions = "CF-1.5" ;
:history = "Exportation netCDF GRASS GIS 7 de r3.out.netcdf" ;
données:
longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
-65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
105, 115, 125, 135, 145, 155, 165, 175 ;
latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
-65, -75, -85 ;
z = 0.5, 1.5, 2.5, 3.5, 4.5 ;
}
################################################ ###################
# Attribution du temps comme unité verticale et définition d'un horodatage absolu
################################################ ###################
r3.timestamp map=volume_float date='1er janvier 2001/5 janvier 2001'
r3.support map=volume_float vunit="jours"
r3.out.netcdf --o -p null=-1 entrée=volume_float sortie=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
dimensions:
longitude = 36 ;
latitude = 18 ;
temps = 5 ;
variables:
char crs ;
crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degré\",0.0174532925199433]]" ;
crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degré\",0.0174532925199433]]" ;
crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
longitude flottante(longitude) ;
longitude:unités = "degrés_est" ;
longitude:long_name = "Valeurs de longitude" ;
longitude:nom_standard = "longitude" ;
longitude:axe = "X" ;
latitude flottante(latitude) ;
latitude:unités = "degrés_nord" ;
latitude:long_name = "Valeurs de latitude" ;
latitude:nom_standard = "latitude" ;
latitude:axe = "Y" ;
int temps(temps) ;
time:units = "jours depuis le 2001-01-01 00:00:00" ;
time:long_name = "Durée en jours" ;
heure:calendrier = "grégorien" ;
heure:positive = "en haut" ;
temps:axe = "T" ;
float volume_float(heure, latitude, longitude) ;
volume_float:valid_min = 3.f ;
volume_float:valid_max = 59.f ;
volume_float:missing_value = -1.f ;
volume_float:_FillValue = -1.f ;
volume_float:grid_mapping = "crs" ;
// attributs globaux :
:Conventions = "CF-1.5" ;
:history = "Exportation netCDF GRASS GIS 7 de r3.out.netcdf" ;
données:
longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
-65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
105, 115, 125, 135, 145, 155, 165, 175 ;
latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
-65, -75, -85 ;
temps = 0, 1, 2, 3, 4 ;
}

Utilisez r3.out.netcdfgrass en ligne à l'aide des services onworks.net