r.sungrass - Trực tuyến trên đám mây

CHƯƠNG TRÌNH:

TÊN

r.sun - Mô hình chiếu xạ và phát xạ mặt trời.
Tính toán bản đồ raster bức xạ mặt trời trực tiếp (chùm), khuếch tán và phản xạ cho một ngày nhất định,
vĩ độ, bề mặt và điều kiện khí quyển. Các thông số mặt trời (ví dụ: mặt trời mọc, mặt trời lặn
thời gian, độ nghiêng, bức xạ ngoài trái đất, độ dài ánh sáng ban ngày) được lưu trong bản đồ
Tập tin lịch sử. Ngoài ra, một giờ địa phương có thể được chỉ định để tính toán tỷ lệ năng lượng mặt trời
bản đồ góc và / hoặc raster bức xạ. Hiệu ứng đổ bóng của địa hình là tùy chọn
kết hợp.

TỪ KHÓA

raster, năng lượng mặt trời, năng lượng mặt trời, bóng

SYNOPSIS

r.sun
r.sun --Cứu giúp
r.sun [-pm] độ cao=chuỗi [khía cạnh=chuỗi] [giá trị khía cạnh=phao] [độ dốc=chuỗi]
[dốc_giá_trị=phao] [trái=chuỗi] [linke_value=phao] [albedo=chuỗi]
[giá trị phản chiếu=phao] [trễ=chuỗi] [Dài=chuỗi] [coeff_bh=chuỗi]
[coeff_dh=chuỗi] [Horiz_basename=tên cơ sở] [chân trời_bước=phao] [sự bù đắp=chuỗi]
[tia_rad=chuỗi] [khác_rad=chuỗi] [ref_rad=chuỗi] [global_rad=chuỗi]
[ina_time=chuỗi] ngày=số nguyên [bước=phao] [độ nghiêng=phao] [thời gian=phao]
[khoảng cách_bước=phao] [phân vùng=số nguyên] [Civil_time=phao] [-ghi đè]
[-giúp đỡ] [-dài dòng] [-yên tĩnh] [-ui]

Cờ:
-p
Không kết hợp hiệu ứng đổ bóng của địa hình

-m
Sử dụng phiên bản bộ nhớ thấp của chương trình

- 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ố:
độ cao=chuỗi [yêu cầu]
Tên của bản đồ raster độ cao đầu vào [mét]

khía cạnh=chuỗi
Tên của bản đồ khía cạnh đầu vào (khía cạnh địa hình hoặc góc phương vị của bảng điều khiển năng lượng mặt trời) [thập phân
độ]

giá trị khía cạnh=phao
Một giá trị duy nhất của hướng (khía cạnh), 270 là hướng nam
Mặc định: 270

độ dốc=chuỗi
Tên của bản đồ raster độ dốc đầu vào (độ dốc địa hình hoặc độ nghiêng của bảng điều khiển năng lượng mặt trời) [thập phân
độ]

dốc_giá_trị=phao
Một giá trị duy nhất của độ nghiêng (độ dốc)
Mặc định: 0.0

trái=chuỗi
Tên của bản đồ đầu vào hệ số đục trong khí quyển Linke [-]

linke_value=phao
Một giá trị duy nhất của hệ số độ đục trong khí quyển Linke [-]
Mặc định: 3.0

albedo=chuỗi
Tên của bản đồ raster đầu vào hệ số albedo mặt đất [-]

giá trị phản chiếu=phao
Một giá trị duy nhất của hệ số albedo trên mặt đất [-]
Mặc định: 0.2

trễ=chuỗi
Tên của bản đồ raster đầu vào có chứa vĩ độ [độ thập phân]

Dài=chuỗi
Tên của bản đồ raster đầu vào có chứa kinh độ [độ thập phân]

coeff_bh=chuỗi
Tên của hệ số bức xạ chùm tia trên bầu trời thực (đám mây dày) bản đồ raster đầu vào [0-1]

coeff_dh=chuỗi
Tên của bản đồ đầu vào hệ số bức xạ khuếch tán (sương mù) trên bầu trời thực [0-1]

Horiz_basename=tên cơ sở
Tên cơ sở bản đồ đầu vào thông tin đường chân trời

chân trời_bước=phao
Kích thước bước góc cho đường chân trời đa hướng [độ]

sự bù đắp=chuỗi
Bản đồ raster góc tới đầu ra (chỉ ở chế độ 1)

tia_rad=chuỗi
Bức xạ chùm tia đầu ra [Wm-2] (chế độ 1) hoặc bản đồ raster bức xạ [Wh.m-2.ngày-1] (chế độ
2)

khác_rad=chuỗi
Bức xạ khuếch tán đầu ra [Wm-2] (chế độ 1) hoặc bản đồ raster bức xạ [Wh.m-2.day-1]
(chế độ 2)

ref_rad=chuỗi
Đầu ra bức xạ phản xạ mặt đất [Wm-2] (chế độ 1) hoặc bản đồ raster bức xạ
[Wh.m-2.day-1] (chế độ 2)

global_rad=chuỗi
Đầu ra bức xạ / bức xạ toàn cầu (tổng số) [Wm-2] (chế độ 1) hoặc
Bản đồ raster bức xạ / bức xạ [Wh.m-2.day-1] (chế độ 2)

ina_time=chuỗi
Bản đồ raster thời gian cách ly đầu ra [h] (chỉ ở chế độ 2)

ngày=số nguyên [yêu cầu]
Số ngày trong năm (1-365)
Tùy chọn: 1-365

bước=phao
Bước thời gian khi tính toán tổng bức xạ cả ngày [giờ thập phân]
Mặc định: 0.5

độ nghiêng=phao
Giá trị giảm dần (ghi đè giá trị được tính toán nội bộ) [radian]

thời gian=phao
Giờ địa phương (năng lượng mặt trời) (chỉ được đặt cho chế độ 1) [giờ thập phân]
Tùy chọn: 0-24

khoảng cách_bước=phao
Hệ số bước khoảng cách lấy mẫu (0.5-1.5)
Mặc định: 1.0

phân vùng=số nguyên
Đọc các tệp đầu vào với số lượng nhiều phần này
Mặc định: 1

Civil_time=phao
Giá trị múi giờ dân sự, nếu không có, thời gian sẽ là giờ mặt trời địa phương

MÔ TẢ

r.sun tính toán bản đồ raster bức xạ mặt trời chiếu tia (trực tiếp), khuếch tán và phản xạ mặt đất
cho các điều kiện ngày, vĩ độ, bề mặt và khí quyển nhất định. Các thông số năng lượng mặt trời (ví dụ: thời gian
của mặt trời mọc và lặn, góc nghiêng, bức xạ ngoài trái đất, độ dài ánh sáng ban ngày) là
được lưu trữ trong các tệp lịch sử của bản đồ kết quả. Ngoài ra, giờ địa phương có thể là
được chỉ định để tính góc tới mặt trời và / hoặc bản đồ raster bức xạ. Bóng tối
ảnh hưởng của địa hình được kết hợp theo mặc định. Điều này có thể được thực hiện trong nội bộ bởi
tính toán hiệu ứng đổ bóng trực tiếp từ mô hình độ cao kỹ thuật số hoặc bằng
xác định bản đồ raster về độ cao đường chân trời nhanh hơn nhiều. Những đường chân trời raster
bản đồ có thể được tính toán bằng cách sử dụng r.horizon.

Đối với các tọa độ kinh độ-vĩ độ, nó yêu cầu bản đồ độ cao tính bằng mét. Các
quy tắc là:

· Tọa độ vĩ độ / kinh độ: độ cao tính bằng mét;

· Các tọa độ khác: độ cao trong cùng một đơn vị với tọa độ đông bắc.
Hình dạng mặt trời của mô hình dựa trên các công trình của Krcho (1990), sau đó được cải tiến bởi
Jenco (1992). Các phương trình mô tả vị trí Mặt trời - Trái đất cũng như sự tương tác của
bức xạ mặt trời với khí quyển ban đầu dựa trên các công thức được đề xuất bởi
Kitler và Mikler (1986). Thành phần này đã được cập nhật đáng kể bởi các kết quả và
đề xuất của nhóm làm việc do Scharmer và Greif (2000) điều phối (thuật toán này
có thể được thay thế bằng thư viện thuật toán SOLPOS có trong GRASS trong lệnh r.sunmask).
Mô hình tính toán tất cả ba thành phần của bức xạ toàn cầu (chùm, khuếch tán và phản xạ)
đối với điều kiện bầu trời quang đãng, tức là không tính đến không gian và thời gian
sự biến đổi của các đám mây. Phạm vi và độ phân giải không gian của khu vực được mô hình hóa, cũng như
tích hợp theo thời gian, chỉ bị giới hạn bởi bộ nhớ và tài nguyên lưu trữ dữ liệu. Các
mô hình được xây dựng để đáp ứng nhu cầu của người dùng trong các lĩnh vực khoa học khác nhau (thủy văn, khí hậu học,
sinh thái học và khoa học môi trường, quang điện, kỹ thuật, v.v.) cho lục địa,
khu vực theo quy mô cảnh quan.

Mô hình coi là hiệu ứng đổ bóng của địa hình cục bộ trừ khi được tắt bằng
các -p cờ. r.sun hoạt động ở hai chế độ: Ở chế độ đầu tiên, nó sẽ tính toán cho địa phương đã đặt
thời gian góc tới mặt trời [độ] và giá trị bức xạ mặt trời [Wm-2]. Trong lần thứ hai
chế độ hàng ngày tổng bức xạ mặt trời [Wh.m-2.ngày-1] được tính trong một ngày nhất định. Bởi một
viết kịch bản cho hai chế độ có thể được sử dụng riêng biệt hoặc kết hợp để cung cấp các ước tính
cho bất kỳ khoảng thời gian mong muốn nào. Mô hình giải quyết sự cản trở bầu trời bằng cách cứu trợ địa phương
Tính năng, đặc điểm. Một số thông số năng lượng mặt trời được lưu trong các tệp lịch sử của bản đồ kết quả,
có thể được xem bằng lệnh r.info.

Bản đồ raster góc tới mặt trời sự bù đắp được tính toán để chỉ định bản đồ độ cao
độ cao, bản đồ raster khía cạnh khía cạnh, bản đồ raster độ dốc dốc, cho ngày ngày
và giờ địa phương thời gian. Không cần xác định vĩ độ cho các vị trí đã biết và
hệ tọa độ / phép chiếu xác định (kiểm tra bằng lệnh g.proj). Nếu bạn có
phép chiếu không xác định, hệ thống (x, y), v.v. thì vĩ độ có thể được xác định rõ ràng cho
các khu vực rộng lớn bằng bản đồ raster đầu vào Latin với các giá trị vĩ độ nội suy. Tất cả raster đầu vào
bản đồ phải là bản đồ raster dấu chấm động (FCELL). Dữ liệu rỗng trong bản đồ bị loại trừ khỏi
tính toán (và cũng tăng tốc độ tính toán), vì vậy mỗi bản đồ raster đầu ra sẽ chứa
dữ liệu rỗng trong các ô theo tất cả các bản đồ raster đầu vào. Người dùng có thể sử dụng lệnh r.null để
tạo / đặt lại tệp null cho bản đồ raster đầu vào của bạn.
Ngày được chỉ định ngày là số của ngày trong năm, trong đó ngày 1 tháng XNUMX là ngày
số 1 ​​và ngày 31 tháng 365 là XNUMX. Thời gian thời gian phải là giờ địa phương (mặt trời) (tức là KHÔNG phải là múi giờ,
ví dụ như GMT, CET) trong hệ thập phân, ví dụ 7.5 (= 7 giờ 30 phút sáng), 16.1 = 4 giờ 6 giờ chiều.

Mặt trời độ nghiêng tham số là một tùy chọn để ghi đè giá trị được tính toán bởi
thông lệ nội bộ cho các ngày trong năm. Giá trị của vĩ độ địa lý có thể được đặt là
một hằng số cho toàn bộ vùng được tính toán hoặc, dưới dạng một tùy chọn, một lưới biểu diễn theo không gian
giá trị phân phối trên một khu vực rộng lớn. Vĩ độ địa lý cũng phải ở dạng thập phân
hệ thống có giá trị dương đối với bán cầu bắc và âm đối với bán cầu nam. Trong
nguyên tắc tương tự yếu tố độ đục Linke (trái, lin ) và albedo mặt đất (albedo, trắng)
có thể thiết lập.

Bên cạnh bức xạ bầu trời trong, người dùng có thể tính toán bức xạ bầu trời thực (chùm, khuếch tán)
sử dụng coeff_bh và coeff_dh bản đồ raster đầu vào xác định phần của
bức xạ bầu trời trong giảm do các yếu tố khí quyển (ví dụ như mây). Giá trị là
giữa 0-1. Thông thường những hệ số này có thể thu được từ một khí tượng dài hạn
các phép đo được cung cấp dưới dạng bản đồ raster với sự phân bố không gian của các hệ số này
riêng biệt cho chùm tia và bức xạ khuếch tán (xem Suri và Hofierka, 2004, phần 3.2).

Bản đồ chiếu xạ mặt trời hoặc raster bức xạ tia_rad, khác_rad, ref_rad được tính toán
trong một ngày nhất định ngày, vĩ độ Latin, độ cao độ cao, dốc độ dốc và khía cạnh khía cạnh
bản đồ raster. Để thuận tiện, đầu ra raster được đưa ra dưới dạng global_rad sẽ xuất ra tổng của
ba thành phần bức xạ. Chương trình sử dụng hệ số độ đục của khí quyển Linke và
hệ số albedo mặt đất. Một giá trị mặc định, đơn lẻ của yếu tố Linke là lin= 3.0 và gần
mức trung bình hàng năm của khu vực nông thôn - thành phố. Yếu tố Linke hoàn toàn rõ ràng
bầu không khí là lin= 1.0. Xem ghi chú bên dưới để tìm hiểu thêm về yếu tố này. Tai nạn
góc mặt trời là góc giữa đường chân trời và vectơ chùm tia mặt trời.

Các bản đồ bức xạ mặt trời cho một ngày nhất định được tính toán bằng cách tích hợp các
bức xạ giữa thời gian mặt trời mọc và mặt trời lặn cho ngày đó. Người dùng có thể đặt mức tốt hơn hoặc
bước thời gian thô hơn được sử dụng để tính toán bức xạ cả ngày với bước Lựa chọn. Các
giá trị mặc định của bước là 0.5 giờ. Các bước lớn hơn (ví dụ: 1.0-2.0) có thể tăng tốc độ tính toán
nhưng tạo ra kết quả kém tin cậy hơn (và nhiều răng cưa hơn). Khi mặt trời di chuyển qua khoảng. 15 °
bầu trời trong một giờ, mặc định bước trong nửa giờ sẽ tạo ra các bước 7.5 ° trong
dữ liệu. Cho đầu ra tương đối mượt mà với mặt trời được đặt ở mọi mức độ chuyển động trong
bầu trời bạn nên thiết lập bước đến 4 phút hoặc ít hơn. bước= 0.05 tương đương với mỗi 3
phút. Tất nhiên, việc đặt bước thời gian rất tốt sẽ làm tăng
thời gian chạy của mô-đun.

Đơn vị đầu ra tính bằng Wh trên mét bình phương mỗi ngày [Wh / (m * m) / ngày]. Tai nạn
bản đồ góc và bức xạ / bức xạ được tính toán với ảnh hưởng che bóng của sự giảm nhẹ
theo mặc định. Cũng có thể tính toán chúng mà không có ảnh hưởng này bằng cách sử dụng
cờ phẳng (-p). Ở những vùng miền núi, điều này có thể dẫn đến những kết quả rất khác! Người dùng
nên lưu ý rằng việc tính đến hiệu ứng bóng mờ của việc giảm nhẹ có thể làm chậm
tốc độ tính toán, đặc biệt là khi độ cao mặt trời thấp.

Khi xem xét hiệu ứng đổ bóng, tốc độ và độ chính xác của tính toán có thể
được kiểm soát bởi khoảng cách_bước tham số xác định mật độ lấy mẫu tại đó
khả năng hiển thị của một tế bào lưới được tính theo hướng của đường đi của dòng năng lượng mặt trời. Nó
cũng xác định phương pháp tính độ cao của chướng ngại vật. Khi chọn một
khoảng cách_bước nhỏ hơn 1.0 (nghĩa là các điểm lấy mẫu sẽ được tính tại khoảng cách_bước *
khoảng cách kích thước ô), r.sun lấy độ cao từ điểm lưới gần nhất. Giá trị trên 1.0
sẽ sử dụng giá trị độ cao tối đa được tìm thấy trong 4 điểm lưới xung quanh gần nhất. Các
giá trị mặc định khoảng cách_bước= 1.0 sẽ cho kết quả hợp lý cho hầu hết các trường hợp (ví dụ: trên
DEM). Các khoảng cách_bước giá trị xác định một hệ số nhân cho khoảng cách lấy mẫu.
Khoảng cách lấy mẫu cơ bản này bằng trung bình cộng của cả hai kích thước ô. Các
giá trị hợp lý nằm trong khoảng 0.5-1.5. Các giá trị dưới 0.5 sẽ giảm và các giá trị
trên 1.0 sẽ làm tăng tốc độ tính toán. Các giá trị lớn hơn 2.0 có thể tạo ra các ước tính
với độ chính xác thấp hơn trong cứu trợ được mổ xẻ cao. Các khu vực bị che khuất hoàn toàn được ghi vào
bản đồ đầu ra dưới dạng giá trị XNUMX. Các khu vực có dữ liệu NULL được coi là không có rào cản với
hiệu ứng đổ bóng.

Các tệp lịch sử của bản đồ được tạo có chứa các thông số được liệt kê sau đây được sử dụng trong
tính toán:
- Hằng số mặt trời 1367 Wm-2
- Bức xạ ngoài Trái đất trên mặt phẳng vuông góc với chùm tia mặt trời [Wm-2]
- Ngày trong năm
- Declination [radian]
- Giờ thập phân (Chỉ thay thế 1)
- Bình minh và hoàng hôn (tối thiểu-tối đa) trên một mặt phẳng nằm ngang
- Độ dài ánh sáng ban ngày
- Vĩ độ địa lý (tối thiểu-tối đa)
- Hệ số độ đục liên kết (tối thiểu-tối đa)
- Albedo mặt đất (tối thiểu-tối đa)

Người dùng có thể sử dụng một shellcript đẹp với ngày thay đổi để tính toán bức xạ trong một thời gian
khoảng thời gian trong năm (ví dụ như thảm thực vật hoặc thời kỳ mùa đông). Độ cao, khía cạnh và độ dốc
các giá trị đầu vào không được phân loại lại thành các loại thô hơn. Điều này có thể dẫn đến
kết quả không chính xác.

LỰA CHỌN

Hiện tại, có hai chế độ r.sun. Trong chế độ đầu tiên, nó tính toán tỷ lệ năng lượng mặt trời
bản đồ raster góc và bức xạ mặt trời sử dụng giờ địa phương đã đặt. Ở chế độ thứ hai hàng ngày
tổng lượng bức xạ mặt trời [Wh.m-2.day-1] được tính cho một ngày cụ thể.

GHI CHÚ

Năng lượng mặt trời là một tham số đầu vào quan trọng trong các mô hình khác nhau liên quan đến năng lượng
công nghiệp, cảnh quan, thảm thực vật, thoát hơi nước, tuyết tan hoặc viễn thám. Mặt trời
bản đồ góc tới tia có thể được sử dụng hiệu quả trong đo bức xạ và địa hình
sửa chữa ở địa hình đồi núi, nơi cần điều tra rất chính xác
đã thực hiện.

Mô hình bức xạ mặt trời trong bầu trời quang đãng được áp dụng trong r.sun dựa trên công việc đã thực hiện
để phát triển Bản đồ Bức xạ Mặt trời Châu Âu (Scharmer và Greif 2000, Page et al.
2001, Rigollier 2001). Mô hình bầu trời quang đãng ước tính bức xạ toàn cầu từ tổng
các thành phần chùm tia, khuếch tán và phản xạ của nó. Sự khác biệt chính giữa bức xạ mặt trời
các mô hình cho bề mặt nghiêng ở Châu Âu là xử lý thành phần khuếch tán. bên trong
Khí hậu châu Âu thành phần này thường là nguồn sai số ước tính lớn nhất. Đang lấy
xem xét các mô hình hiện có và hạn chế của chúng đối với Bức xạ Mặt trời Châu Âu
Nhóm Atlas đã chọn mô hình Muneer (1990) vì nó có cơ sở lý thuyết đúng đắn và do đó
nhiều tiềm năng để cải thiện sau này.

Chi tiết về các phương trình cơ bản được sử dụng trong chương trình này có thể được tìm thấy trong tài liệu tham khảo
tài liệu được trích dẫn bên dưới hoặc sách được xuất bản bởi Neteler và Mitasova: Nguồn mở GIS: A GRASS
Phương pháp tiếp cận GIS (xuất bản trên Kluwer Academic Publishers năm 2002).

Giá trị trung bình hàng tháng của hệ số độ đục Linke đối với khí hậu ôn hòa ở
bán cầu bắc (xem tài liệu tham khảo cho khu vực nghiên cứu của bạn):

Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX Tháng XNUMX hàng năm

núi 1.5 1.6 1.8 1.9 2.0 2.3 2.3 2.3 2.1 1.8 1.6 1.5 1.90

nông thôn 2.1 2.2 2.5 2.9 3.2 3.4 3.5 3.3 2.9 2.6 2.3 2.2 2.75

thành phố 3.1 3.2 3.5 4.0 4.2 4.3 4.4 4.3 4.0 3.6 3.3 3.1 3.75

công nghiệp 4.1 4.3 4.7 5.3 5.5 5.7 5.8 5.7 5.3 4.9 4.5 4.2 5.00 XNUMX XNUMX XNUMX XNUMX XNUMX XNUMX

Các cải tiến theo kế hoạch bao gồm việc sử dụng thuật toán SOLPOS cho hình học năng lượng mặt trời
tính toán và tính toán nội bộ của khía cạnh và độ dốc.

Năng lượng mặt trời thời gian
Theo mặc định, r.sun tính thời gian là giờ mặt trời thực, theo đó giữa trưa mặt trời luôn chính xác
12 giờ ở khắp mọi nơi trong khu vực hiện tại. Tùy thuộc vào khu vực quan tâm ở đâu
nằm trong múi giờ liên quan, điều này có thể gây ra chênh lệch lên đến một giờ, trong một số
các trường hợp (như Tây Ban Nha) thậm chí nhiều hơn. Ngoài ra, phần bù thay đổi trong năm
theo phương trình thời gian.

Để khắc phục sự cố này, người dùng có thể sử dụng tùy chọn Civil_time = in
r. chạy để làm cho nó sử dụng thời gian trong thế giới thực (đồng hồ treo tường). Ví dụ: đối với Trung Âu,
Độ lệch múi giờ là +1, +2 khi thời gian tiết kiệm ánh sáng ban ngày có hiệu lực.

Khai thác of bóng tối bản đồ
Bản đồ bóng đổ có thể được trích xuất từ ​​bản đồ góc tới của mặt trời (incidout). Khu vực
với các giá trị XNUMX sẽ bị che khuất. Điều này sẽ không hoạt động nếu -p cờ đã được sử dụng.

Chó cái bản đồ và ra of trí nhớ vấn đề
Với một số lượng lớn hoặc các cột và hàng, r.sun có thể tiêu tốn lượng bộ nhớ đáng kể.
Mặc dù bản đồ raster đầu ra không thể phân vùng, nhưng bản đồ raster đầu vào đang sử dụng
phân vùng tham số. Trong trường hợp lỗi hết bộ nhớ (ERROR: G_malloc: hết bộ nhớ),
các phân vùng tham số có thể được sử dụng để chạy một phép tính được phân đoạn tiêu thụ ít hơn
bộ nhớ trong quá trình tính toán. Dung lượng bộ nhớ bằng r.sun được ước tính như sau:
# mà không cần phân vùng bản đồ raster đầu vào:
# yêu cầu bộ nhớ: 4 byte cho mỗi ô raster
# row, cols: hàng và cột của khu vực hiện tại (tìm hiểu với g.region)
# IR: số lượng bản đồ raster đầu vào không có bản đồ đường chân trời
# HOẶC: số lượng bản đồ raster đầu ra
memory_bytes = row * cols * (IR * 4 + Horizon_steps + OR * 4)
# với phân vùng bản đồ raster đầu vào:
memory_bytes = row * cols * ((IR * 4 + chân trời_ bước) / npartitions + OR * 4)

VÍ DỤ

Ví dụ về Bắc Carolina (xem xét cả bóng đổ):
g.region raster = elevation -p
# tính toán góc đường chân trời (để tăng tốc độ tính toán r.sun tiếp theo)
r.horizon elevation = elevation step = 30 bufferzone = 200 basename = horangle \
maxdistance = 5000
# độ dốc + khía cạnh
r.slope.aspect elevation = elevationpect = khía cạnh.dem độ dốc = độ dốc.dem
# tính toán bức xạ toàn cầu cho ngày 180 lúc 2 giờ chiều, sử dụng đầu ra r.horizon
r.sun elevation = elevation height_basename = horangle Horizon_step = 30 \
khía cạnh = khía cạnh.dem độ dốc = độ dốc.dem global_rad = ngày toàn cầu = 180 thời gian = 14
# kết quả: đầu ra bức xạ / bức xạ toàn cầu (tổng cộng) [Wm-2] cho ngày / giờ nhất định
r.univar toàn cầu_rad

Tính toán lượng chiếu xạ hàng ngày tích hợp cho một khu vực ở Bắc Carolina cho một
ngày trong năm ở độ phân giải 30m. Đây là ngày 172 (tức là ngày 21 tháng XNUMX năm không nhuận):
g.region raster = elev_ned_30m -p
# xem xét đổ bóng
r.sun elevation = elev_ned_30m linke_value = 2.5 albedo_value = 0.2 ngày = 172 \
beam_rad = b172 diff_rad = d172 \
refl_rad = r172 secure_time = it172
d.mon wx0.
# hiển thị bản đồ raster chiếu xạ [Wh.m-2.day-1]
d.rast.leg b172
# hiển thị bản đồ thời gian cách ly [h]
d.rast.leg it172
Chúng ta có thể tính ngày trong năm từ một ngày cụ thể trong Python shell:
>>> nhập ngày giờ
>>> datetime.datetime (2014, 6, 21) .timetuple (). tm_yday
172

Sử dụng r.sungrass trực tuyến bằng các dịch vụ onworks.net