#!/usr/bin/env python

#

##############################################################################

#

# MODULE: r.skyview

#

# AUTHOR(S): Anna Petrasova (kratochanna gmail.com)

# Vaclav Petras <wenzeslaus gmail com> (colorize addition)

# Steven Pawley <steven.pawley gmail.com> (parallelization)

#

# PURPOSE: Implementation of Sky-View Factor visualization technique

#

# COPYRIGHT: (C) 2013-2021 by the GRASS Development Team

#

# This program is free software under the GNU General Public

# License (version 2). Read the file COPYING that comes with GRASS

# for details.

#

##############################################################################

#%module

#% description: Computes Sky-View Factor visualization technique

#% keyword: raster

#% keyword: visualization

#%end

#%option G_OPT_R_INPUT

#%end

#%option G_OPT_R_OUTPUT

#%end

#%option

#% key: ndir

#% description: Number of directions (8 to 32 recommended)

#% type: integer

#% required: yes

#% answer: 16

#% options: 2-360

#%end

#%option

#% key: maxdistance

#% description: The maximum distance to consider when finding the horizon height

#% type: double

#% required: no

#%end

#%option

#% key: color_source

#% type: string

#% label: Source raster for colorization

#% description: Input and color_input are taken from input and color_input options respectively. The rest is computed using r.slope.aspect

#% descriptions: input; use the raster from the input option;color_input;use the raster from the color_input option;slope;compute and use slope;aspect;compute and use aspect;dxy;compute and use second order partial derivative dxy

#% multiple: no

#% required: no

#% options: input, color_input, slope, aspect, dxy

#% answer: input

#% guisection: Colorize

#%end

#%option G_OPT_R_INPUT

#% key: color_input

#% required: no

#% description: Custom raster map to be used for colorization

#% guisection: Colorize

#%end

#%option

#% key: color_table

#% type: string

#% label: Color table for colorization raster (preset color table by default)

#% description: If empty, the color table of the created raster is used (not used at all for input and color_input)

#% multiple: no

#% required: no

#% options: reds, blues, greens, oranges, sepia, aspectcolr

#% guisection: Colorize

#%end

#%option G_OPT_R_OUTPUT

#% key: colorized_output

#% required: no

#% description: Colorized sky-view factor

#% guisection: Colorize

#%end

#%option

#% key: basename

#% type: string

#% multiple: no

#% description: Set the basename for the intermediate maps

#%end

#%option

#% key: processes

#% type: integer

#% multiple: no

#% description: Number of processes to run in parallel

#% answer: 1

#%end

#%flag

#% key: o

#% label: Compute openness instead of skyview factor

#% description: Openness considers zenith angles > 90 degrees

#%end

#%flag

#% key: n

#% label: Invert color table for colorization raster

#% description: Ignored for input and color_input

#% guisection: Colorize

#%end

#%flag

#% key: p

#% label: Do not process in parallel

#% description: Do not process horizon rasters in parallel

#%end

import sys

import os

import atexit

import numpy as np

from grass.exceptions import CalledModuleError

import grass.script.core as gcore

import grass.script.raster as grast

from grass.pygrass.messages import get_msgr

from grass.pygrass.modules import Module, ParallelModuleQueue

from copy import deepcopy

# TODO: also used for r.slope.aspect result

TMP_NAME = "tmp_horizon_" + str(os.getpid())

CLEANUP = True

def cleanup():

)

def main():

return 0

def _get_horizon_maps():

return gcore.list_grouped("rast", pattern=TMP_NAME + "*")[gcore.gisenv()["MAPSET"]]

def check_map_name(name):

)

if __name__ == "__main__":

options, flags = gcore.parser()

atexit.register(cleanup)

sys.exit(main())