Đây là lệnh r3.out.netcdfgrass có thể được chạy trong nhà cung cấp dịch vụ lưu trữ miễn phí OnWorks bằng cách sử dụng một trong nhiều máy trạm trực tuyến miễn phí của chúng tôi như Ubuntu Online, Fedora Online, trình giả lập trực tuyến Windows hoặc trình mô phỏng trực tuyến MAC OS
CHƯƠNG TRÌNH:
TÊN
r3.out.netcdf - Xuất bản đồ raster 3D dưới dạng tệp netCDF.
TỪ KHÓA
raster3d, xuất, netCDF, voxel
SYNOPSIS
r3.out.netcdf
r3.out.netcdf --Cứu giúp
r3.out.netcdf [-pm] đầu vào=tên đầu ra=tên [vô giá trị=phao] [-ghi đè] [-giúp đỡ]
[-dài dòng] [-yên tĩnh] [-ui]
Cờ:
-p
Xuất thông tin chiếu dưới dạng tham số wkt và proj4
-m
Sử dụng mặt nạ raster 3D (nếu tồn tại) với bản đồ đầu vào
- ghi đè
Cho phép các tệp đầu ra ghi đè lên các tệp hiện có
--Cứu giúp
In tóm tắt sử dụng
--dài dòng
Đầu ra mô-đun dài dòng
--Yên lặng
Đầu ra mô-đun yên tĩnh
--ui
Buộc khởi chạy hộp thoại GUI
Tham số:
đầu vào=tên [yêu cầu]
Tên của bản đồ raster 3D đầu vào
đầu ra=tên [yêu cầu]
Tên cho tệp đầu ra netCDF
vô giá trị=phao
Giá trị được sử dụng cho các giá trị null, mặc định là tiêu chuẩn netCDF
MÔ TẢ
Mô-đun r3.out.netcdf xuất bản đồ raster 3D dưới dạng tệp netCDF. Bản đồ là kẻ vạch 3D hợp lệ
bản đồ trong đường dẫn tìm kiếm tập bản đồ hiện tại. Các đầu ra tham số là tên của tệp netCDF
sẽ được viết trong thư mục làm việc hiện tại. Mô-đun nhạy cảm với tập bản đồ
cài đặt khu vực (đặt bằng g. khu vực) và cài đặt đơn vị dọc của bản đồ 3D raster.
Đơn vị dọc có thể được đặt bằng r3.support.
Chú thích
Tệp netCDF kết quả sẽ có thời gian là chiều thứ ba trong trường hợp là chiều dọc tạm thời
đơn vị được chỉ định cho bản đồ đầu vào, nếu không thì kích thước thứ ba được giả định là không gian.
Nếu bản đồ raster 3D có dấu thời gian với thời gian tuyệt đối, thì thời gian bắt đầu sẽ được sử dụng như
điểm bắt đầu cho chuỗi thời gian netCDF. Nếu dấu thời gian tuyệt đối bị thiếu ngày
"Ngày 01 tháng 1900. 00 00:00:XNUMX" sẽ được sử dụng làm mặc định. Trong trường hợp tem thời gian tương đối không bắt đầu
thời gian được đặt trong tệp netCDF. Số độ sâu của bản đồ 3D raster là số
thời gian trong tệp netCDF. Bạn có thể dùng t.rast.to.rast3 chuyển đổi thời gian không gian
tập dữ liệu raster thành một bản đồ raster 3D đại diện cho một khối voxel thời gian không gian, sau đó xuất
chuỗi thời gian raster dưới dạng tệp netCDF.
Tọa độ không gian được xuất dưới dạng tọa độ căn giữa ô. Phép chiếu có thể được
được lưu trữ tùy chọn trong siêu dữ liệu như crs thuộc tính . Siêu dữ liệu chiếu netCDF
lưu trữ tuân theo gợi ý spatial_ref GDAL / netCDF ở đây và netCDF CF 1.6
quy ước ở đây bằng cách sử dụng thông tin chiếu WKT. Bổ sung một chuỗi PROJ.4 được lưu trữ trong
các crs phần thuộc tính. Việc xuất các tham số chiếu bị chặn khi
Phép chiếu XY được thiết lập.
Phạm vi của bản đồ raster 3D được đặt trong tệp đầu ra netCDF. Tùy chọn là netCDF
giá trị bị thiếu và _FillValue có thể được đặt bằng cách sử dụng vô giá trị Lựa chọn. Theo mặc định, netCDF
dấu chấm động _FillValues được sử dụng để điền vào ô trống.
Thư viện netCDF phải được cài đặt trên hệ thống và được kích hoạt tại thời điểm cấu hình
biên dịch mô-đun này.
Các ví dụ
g. vùng s = -90 n = 90 w = -180 e = 180 b = 0 t = 5 res = 10 res3 = 10 tbres = 1 -p3
###################################################### ####################
# Chúng tôi tạo một bản đồ khối lượng đơn giản với các giá trị dấu chấm động
###################################################### ####################
r3.mapcalc --o expr = "volume_float = float (col () + row () + depth ())"
r3.info Volume_float
+ ------------------------------------------------- --------------------------- +
| Layer: volume_float Ngày: Thứ Sáu 14 tháng 08 40:56:2012 XNUMX |
| Mapset: PERMANENT Đăng nhập của Người tạo: soeren |
| Vị trí: TestLL |
| DataBase: / 1 / soeren / grassdata |
| Tiêu đề: volume_float |
| Đơn vị: không có |
| Đơn vị dọc: đơn vị |
| Dấu thời gian: không có |
| ------------------------------------------------- --------------------------- |
| |
| Loại bản đồ: ô 3d Số danh mục: 0 |
| Loại dữ liệu: FCELL |
| Hàng: 18 |
| Các cột: 36 |
| Độ sâu: 5 |
| Tổng số ô: 3240 |
| Tổng kích thước: 10114 Byte |
| Số lượng gạch: 1 |
| Kích thước gạch trung bình: 10114 Byte |
| Kích thước ô trong bộ nhớ: 12960 Byte |
| Số ô x, y và z: 1, 1, 1 |
| Kích thước của ô theo x, y, z: 36, 18, 5 |
| |
| Phép chiếu: Vĩ độ-Kinh độ (vùng 0) |
| N: 90N S: 90S Độ phân giải: 10 |
| E: 180E W: 180W Độ phân giải: 10 |
| T: 5 B: 0 Độ phân giải: 1 |
| Dải dữ liệu: min = 3 max = 59 |
| |
| Nguồn dữ liệu: |
| |
| |
| |
| Mô tả dữ liệu: |
| được tạo bởi r3.mapcalc |
| |
| Bình luận: |
| r3.mapcalc expression = "volume_float = float (col () + row () + depth ())" |
| |
+ ------------------------------------------------- --------------------------- +
###################################################### ####################
# Chúng tôi sử dụng công cụ ncdump netCDF để xem tiêu đề
# và tọa độ của tệp netCDF đã xuất
###################################################### ####################
r3.out.netcdf --o input = volume_float output = test_float.nc
ncdump -c test_float.nc
netcdf test_float {
kích thước:
kinh độ = 36;
vĩ độ = 18;
z = 5;
biến:
float longitude (kinh độ);
kinh độ: đơn vị = "độ_độ";
longitude: long_name = "Giá trị kinh độ";
kinh độ: standard_name = "longitude";
kinh độ: axis = "X";
vĩ độ nổi (vĩ độ);
vĩ độ: đơn vị = "độ_ chu vi";
latitude: long_name = "Giá trị vĩ độ";
vĩ độ: standard_name = "latitude";
vĩ độ: axis = "Y";
float z (z);
z: đơn vị = "mét";
z: long_name = "z tọa độ của hình chiếu";
z: standard_name = "chiếu_z_cấp";
z: positive = "up";
z: axis = "Z";
float volume_float (z, vĩ độ, kinh độ);
volume_float: valid_min = 3.f;
volume_float: valid_max = 59.f;
// thuộc tính toàn cục:
: Quy ước = "CF-1.5";
: history = "GRASS GIS 7 netCDF export of r3.out.netcdf";
dữ liệu:
kinh độ = -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;
vĩ độ = 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;
}
###################################################### ####################
# Xuất cài đặt cụ thể về phép chiếu và giá trị null
###################################################### ####################
r3.out.netcdf --o -p null = -1 input = volume_float output = test_float.nc
ncdump -c test_float.nc
netcdf test_float {
kích thước:
kinh độ = 36;
vĩ độ = 18;
z = 5;
biến:
ký tự;
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";
float longitude (kinh độ);
kinh độ: đơn vị = "độ_độ";
longitude: long_name = "Giá trị kinh độ";
kinh độ: standard_name = "longitude";
kinh độ: axis = "X";
vĩ độ nổi (vĩ độ);
vĩ độ: đơn vị = "độ_ chu vi";
latitude: long_name = "Giá trị vĩ độ";
vĩ độ: standard_name = "latitude";
vĩ độ: axis = "Y";
float z (z);
z: đơn vị = "mét";
z: long_name = "z tọa độ của hình chiếu";
z: standard_name = "chiếu_z_cấp";
z: positive = "up";
z: axis = "Z";
float volume_float (z, vĩ độ, kinh độ);
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";
// thuộc tính toàn cục:
: Quy ước = "CF-1.5";
: history = "GRASS GIS 7 netCDF export of r3.out.netcdf";
dữ liệu:
kinh độ = -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;
vĩ độ = 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;
}
###################################################### ####################
# Gán thời gian dưới dạng đơn vị dọc và đặt dấu thời gian tuyệt đối
###################################################### ####################
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 input = volume_float output = test_float.nc
ncdump -c test_float.nc
netcdf test_float {
kích thước:
kinh độ = 36;
vĩ độ = 18;
thời gian = 5;
biến:
ký tự;
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";
float longitude (kinh độ);
kinh độ: đơn vị = "độ_độ";
longitude: long_name = "Giá trị kinh độ";
kinh độ: standard_name = "longitude";
kinh độ: axis = "X";
vĩ độ nổi (vĩ độ);
vĩ độ: đơn vị = "độ_ chu vi";
latitude: long_name = "Giá trị vĩ độ";
vĩ độ: standard_name = "latitude";
vĩ độ: axis = "Y";
int time (thời gian);
time: units = "days kể từ 2001-01-01 00:00:00";
time: long_name = "Thời gian tính bằng ngày";
thời gian: lịch = "gregorian";
time: positive = "up";
time: axis = "T";
float volume_float (thời gian, vĩ độ, kinh độ);
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";
// thuộc tính toàn cục:
: Quy ước = "CF-1.5";
: history = "GRASS GIS 7 netCDF export of r3.out.netcdf";
dữ liệu:
kinh độ = -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;
vĩ độ = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
-65, -75, -85;
thời gian = 0, 1, 2, 3, 4;
}
Sử dụng r3.out.netcdfgrass trực tuyến bằng các dịch vụ onworks.net