这是 r3.out.netcdfgrass 命令，可以使用我们的多个免费在线工作站之一在 OnWorks 免费托管服务提供商中运行，例如 Ubuntu Online、Fedora Online、Windows 在线模拟器或 MAC OS 在线模拟器
程序：
您的姓名
r3.out.netcdf文件 - 将 3D 光栅图导出为 netCDF 文件。
关键词
raster3d、导出、netCDF、体素
概要
r3.out.netcdf文件
r3.out.netcdf文件 - 帮帮我
r3.out.netcdf文件 [-pm] 输入=姓名 产量=姓名 [空=浮动] [--覆盖] [--帮助]
[--详细] [--安静] [--ui]
标志：
-p
将投影信息导出为 wkt 和 proj4 参数
-m
对输入地图使用 3D 光栅掩码（如果存在）
--覆盖
允许输出文件覆盖现有文件
- 帮帮我
打印使用摘要
--详细
详细模块输出
- 安静的
静音模块输出
--用户界面
强制启动 GUI 对话框
参数：
输入=姓名 [必需的]
输入 3D 栅格地图的名称
产量=姓名 [必需的]
netCDF 输出文件的名称
空=浮动
用于空值的值，默认为 netCDF 标准
商品描述
模块 r3.out.netcdf文件 将 3D 光栅图导出为 netCDF 文件。 地图是有效的 3D 栅格
当前地图集搜索路径中的地图。 这 产量 参数是 netCDF 文件的名称
将写入当前工作目录中。 模块对mapset敏感
区域设置（设置为 区域) 和 3D 栅格地图的垂直单位设置。
可以使用 r3.support 设置垂直单位。
备注
生成的 netCDF 文件将有时间作为第三维，以防时间垂直
单位为输入地图指定，否则第三维被假定为空间。
如果 3D 栅格地图具有带绝对时间的时间戳，则开始时间将用作
netCDF 时间序列的起点。 如果绝对时间戳缺少日期
“Jan. 01. 1900 00:00:00”将用作默认值。 在相对时间戳没有开始的情况下
时间在 netCDF 文件中设置。 3D光栅图的深度数为数
netCDF 文件中的时间片。 您可以使用 t.rast.to.rast3 转换时空
将栅格数据集转换为表示时空体素立方体的 3D 栅格地图，然后导出
栅格时间序列作为 netCDF 文件。
空间坐标导出为以单元格为中心的坐标。 投影可以是
可选地存储在元数据中作为 CRS 属性 。 netCDF 投影元数据
存储遵循此处的 spatial_ref GDAL/netCDF 建议和 netCDF CF 1.6
这里使用 WKT 投影信息的约定。 另外一个 PROJ.4 字符串存储在
这些因素包括原料奶的可用性以及达到必要粉末质量水平所需的工艺。 CRS 属性部分。 投影参数的导出被抑制时
XY 投影已设置。
3D 光栅图的范围在 netCDF 输出文件中设置。 可选的 netCDF
缺失值和 _FillValue 可以使用 空 选项。 默认情况下，netCDF
浮点 _FillValues 用于填充空单元格。
netCDF 库必须安装在系统上并在配置时激活
编译这个模块。
例子
g.区域 s=-90 n=90 w=-180 e=180 b=0 t=5 res=10 res3=10 tbres=1 -p3
############################################### ###################
# 我们创建一个带有浮点值的简单体积图
############################################### ###################
r3.mapcalc --o expr="volume_float = float(col() + row() + depth())"
r3.info 卷_浮点
+------------------------------------------------- ---------------------------+
| 层：volume_float 日期：Thu Jun 14 08:40:56 2012 |
| 地图集：永久登录创建者：soeren |
| 地点：TestLL |
| 数据库：/1/soeren/grassdata |
| 标题：volume_float |
| 单位：无 |
| 垂直单位：单位|
| 时间戳：无 |
|------------------------------------------------- ---------------------------|
| |
| 地图类型：3d 单元格类别数量：0 |
| 数据类型：FCELL |
| 行数：18 |
| 列：36 |
| 深度：5 |
| 总细胞数：3240 |
| 总大小：10114 字节 |
| 瓷砖数量：1 |
| 平均分块大小：10114 字节 |
| 内存中的切片大小：12960 字节 |
| x、y 和 z 中的图块数：1, 1, 1 |
| 瓷砖在 x、y、z 方向的尺寸：36、18、5 |
| |
| 投影：经纬度（区域 0）|
| N：90N S：90S 分辨率：10 |
| E：180E W：180W 分辨率：10 |
| T：5 B：0 分辨率：1 |
| 数据范围：min = 3 max = 59 |
| |
| 数据来源：|
| |
| |
| |
| 数据说明： |
| 由 r3.mapcalc 生成 |
| |
| 评论： |
| r3.mapcalc 表达式="volume_float = float(col() + row() + depth())" |
| |
+------------------------------------------------- ---------------------------+
############################################### ###################
# 我们使用netCDF ncdump工具来查看header
# 和导出的 netCDF 文件的坐标
############################################### ###################
r3.out.netcdf --o 输入=volume_float 输出=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
尺寸：
经度 = 36 ;
纬度 = 18 ;
z = 5 ；
变量：
浮动经度（经度）；
经度：单位=“度数东”；
longitude:long_name = "经度值";
经度：标准名称=“经度”；
经度：轴=“X”；
浮动纬度（纬度）；
纬度：单位=“度数_北”；
latitude:long_name = "纬度值" ;
纬度：标准名称=“纬度”；
纬度：轴=“Y”；
浮动 z(z) ;
z：单位=“米”；
z:long_name = "投影的z坐标";
z:standard_name = "projection_z_coordinate" ;
z：正=“向上”；
z:axis = "Z" ;
float volume_float(z, 纬度, 经度) ;
Volume_float:valid_min = 3.f ;
volume_float:valid_max = 59.f;
// 全局属性：
:Conventions = "CF-1.5" ;
:history = "GRASS GIS 7 netCDF 导出 r3.out.netcdf";
数据：
经度 = -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；
纬度 = 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；
}
############################################### ###################
# 导出投影特定设置和空值
############################################### ###################
r3.out.netcdf --o -p null=-1 输入=volume_float 输出=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
尺寸：
经度 = 36 ;
纬度 = 18 ;
z = 5 ；
变量：
字符 crs ;
crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0",UNIT度\",0.0174532925199433]]" ;
crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[度\",0.0174532925199433]]" ;
crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000";
浮动经度（经度）；
经度：单位=“度数东”；
longitude:long_name = "经度值";
经度：标准名称=“经度”；
经度：轴=“X”；
浮动纬度（纬度）；
纬度：单位=“度数_北”；
latitude:long_name = "纬度值" ;
纬度：标准名称=“纬度”；
纬度：轴=“Y”；
浮动 z(z) ;
z：单位=“米”；
z:long_name = "投影的z坐标";
z:standard_name = "projection_z_coordinate" ;
z：正=“向上”；
z:axis = "Z" ;
float volume_float(z, 纬度, 经度) ;
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";
// 全局属性：
:Conventions = "CF-1.5" ;
:history = "GRASS GIS 7 netCDF 导出 r3.out.netcdf";
数据：
经度 = -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；
纬度 = 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；
}
############################################### ###################
# 指定时间为垂直单位并设置绝对时间戳
############################################### ###################
r3.timestamp map=volume_float date='1 Jan 2001/5 Jan 2001'
r3.support map=volume_float vunit="days"
r3.out.netcdf --o -p null=-1 输入=volume_float 输出=test_float.nc
ncdump -c test_float.nc
netcdf test_float {
尺寸：
经度 = 36 ;
纬度 = 18 ;
时间 = 5 ;
变量：
字符 crs ;
crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0",UNIT度\",0.0174532925199433]]" ;
crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[度\",0.0174532925199433]]" ;
crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000";
浮动经度（经度）；
经度：单位=“度数东”；
longitude:long_name = "经度值";
经度：标准名称=“经度”；
经度：轴=“X”；
浮动纬度（纬度）；
纬度：单位=“度数_北”；
latitude:long_name = "纬度值" ;
纬度：标准名称=“纬度”；
纬度：轴=“Y”；
整数时间（时间）；
time:units = "自 2001-01-01 00:00:00 以来的天数";
time:long_name = "时间天数";
时间：日历=“公历”；
时间：正=“向上”；
时间：轴=“T”；
浮动volume_float（时间，纬度，经度）；
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";
// 全局属性：
:Conventions = "CF-1.5" ;
:history = "GRASS GIS 7 netCDF 导出 r3.out.netcdf";
数据：
经度 = -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；
纬度 = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
-65, -75, -85 ;
时间 = 0, 1, 2, 3, 4 ;
}
使用 onworks.net 服务在线使用 r3.out.netcdfgrass
