WBT_tutorial_geomorfometria.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "marine-gregory",
   "metadata": {},
   "source": [
    "# WhiteboxTools no Python: Geomorfometria\n",
    "\n",
    "### Prof. Carlos H. Grohmann - IEE-USP  \n",
    "\n",
    "### Exemplo de uso do WhiteboxTools com análise de um Modelo de Elevação  \n",
    "\n",
    "adaptado de: https://www.whiteboxgeo.com/manual/wbt_book/python_scripting/scripting.html\n",
    "\n",
    "veja também: https://github.com/opengeos/whitebox-python"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "arbitrary-kelly",
   "metadata": {},
   "source": [
    "## Instalação WBT e setup"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "seasonal-driver",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Collecting whitebox\n",
      "  Downloading whitebox-2.3.1-py2.py3-none-any.whl (72 kB)\n",
      "\u001b[2K     \u001b[38;2;114;156;31m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m72.1/72.1 kB\u001b[0m \u001b[31m3.4 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
      "\u001b[?25hRequirement already satisfied: Click>=6.0 in /home/guano/Documents/mcenv/lib/python3.11/site-packages (from whitebox) (8.1.6)\n",
      "Installing collected packages: whitebox\n",
      "Successfully installed whitebox-2.3.1\n"
     ]
    }
   ],
   "source": [
    "# Lembre-se que você precisa ter o WhiteboxTools instalado. Se não tiver, pode fazer o download na próxima célula\n",
    "!pip install whitebox"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "searching-region",
   "metadata": {},
   "outputs": [],
   "source": [
    "# setup inicial - importe o módulo e defina a localização do executavel (ou faça o download)\n",
    "import whitebox\n",
    "# whitebox.download_wbt(linux_musl=True, reset=True)\n",
    "wbt = whitebox.WhiteboxTools()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "arbitrary-simon",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "wbt.set_whitebox_dir('/mnt/c/Users/PC/Documents/py3/lib/python3.8/site-packages/whitebox/WBT')\n",
    "wbt.set_working_dir('/mnt/d/Dropbox/USP/extensao/videos_whitebox')\n",
    "\n",
    "# wbt.set_whitebox_dir('/home/guano/Documents/mcenv/lib/python3.11/site-packages/whitebox')\n",
    "# wbt.set_working_dir('/home/guano/Dropbox/USP/extensao/videos_whitebox')\n",
    "\n",
    "# Setting the following to True enables tools to output DEFLATE compressed GeoTIFFs.\n",
    "# You only need to do this once, if you wish all tools to compress their raster\n",
    "# outputs.\n",
    "wbt.set_compress_rasters(True) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "conceptual-ground",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'absolute_value': 'Calculates the absolute value of every cell in a raster.',\n",
       " 'adaptive_filter': 'Performs an adaptive filter on an image.',\n",
       " 'add': 'Performs an addition operation on two rasters or a raster and a constant value.',\n",
       " 'add_point_coordinates_to_table': \"Modifies the attribute table of a point vector by adding fields containing each point's X and Y coordinates.\",\n",
       " 'aggregate_raster': 'Aggregates a raster to a lower resolution.',\n",
       " 'and': 'Performs a logical AND operator on two Boolean raster images.',\n",
       " 'anova': 'Performs an analysis of variance (ANOVA) test on a raster dataset.',\n",
       " 'arc_cos': 'Returns the inverse cosine (arccos) of each values in a raster.',\n",
       " 'arc_sin': 'Returns the inverse sine (arcsin) of each values in a raster.',\n",
       " 'arc_tan': 'Returns the inverse tangent (arctan) of each values in a raster.',\n",
       " 'arcosh': 'Returns the inverse hyperbolic cosine (arcosh) of each values in a raster.',\n",
       " 'arsinh': 'Returns the inverse hyperbolic sine (arsinh) of each values in a raster.',\n",
       " 'artanh': 'Returns the inverse hyperbolic tangent (arctanh) of each values in a raster.',\n",
       " 'ascii_to_las': 'Converts one or more ASCII files containing LiDAR points into LAS files.',\n",
       " 'aspect': 'Calculates an aspect raster from an input DEM.',\n",
       " 'atan2': 'Returns the 2-argument inverse tangent (atan2).',\n",
       " 'attribute_correlation': 'Performs a correlation analysis on attribute fields from a vector database.',\n",
       " 'attribute_correlation_neighbourhood_analysis': 'Performs a correlation on two input vector attributes within a neighbourhood search windows.',\n",
       " 'attribute_histogram': \"Creates a histogram for the field values of a vector's attribute table.\",\n",
       " 'attribute_scattergram': \"Creates a scattergram for two field values of a vector's attribute table.\",\n",
       " 'average_flowpath_slope': 'Measures the average slope gradient from each grid cell to all upslope divide cells.',\n",
       " 'average_normal_vector_angular_deviation': 'Calculates the circular variance of aspect at a scale for a DEM.',\n",
       " 'average_overlay': 'Calculates the average for each grid cell from a group of raster images.',\n",
       " 'average_upslope_flowpath_length': 'Measures the average length of all upslope flowpaths draining each grid cell.',\n",
       " 'balance_contrast_enhancement': 'Performs a balance contrast enhancement on a colour-composite image of multispectral data.',\n",
       " 'basins': 'Identifies drainage basins that drain to the DEM edge.',\n",
       " 'bilateral_filter': 'A bilateral filter is an edge-preserving smoothing filter introduced by Tomasi and Manduchi (1998).',\n",
       " 'block_maximum_gridding': 'Creates a raster grid based on a set of vector points and assigns grid values using a block maximum scheme.',\n",
       " 'block_minimum_gridding': 'Creates a raster grid based on a set of vector points and assigns grid values using a block minimum scheme.',\n",
       " 'boundary_shape_complexity': 'Calculates the complexity of the boundaries of raster polygons.',\n",
       " 'breach_depressions': \"Breaches all of the depressions in a DEM using Lindsay's (2016) algorithm. This should be preferred over depression filling in most cases.\",\n",
       " 'breach_depressions_least_cost': 'Breaches the depressions in a DEM using a least-cost pathway method.',\n",
       " 'breach_single_cell_pits': 'Removes single-cell pits from an input DEM by breaching.',\n",
       " 'buffer_raster': 'Maps a distance-based buffer around each non-background (non-zero/non-nodata) grid cell in an input image.',\n",
       " 'burn_streams_at_roads': 'Burns-in streams at the sites of road embankments.',\n",
       " 'ceil': 'Returns the smallest (closest to negative infinity) value that is greater than or equal to the values in a raster.',\n",
       " 'centroid': 'Calculates the centroid, or average location, of raster polygon objects.',\n",
       " 'centroid_vector': 'Identifies the centroid point of a vector polyline or polygon feature or a group of vector points.',\n",
       " 'change_vector_analysis': 'Performs a change vector analysis on a two-date multi-spectral dataset.',\n",
       " 'circular_variance_of_aspect': 'Calculates the circular variance of aspect at a scale for a DEM.',\n",
       " 'classify_buildings_in_lidar': 'Reclassifies a LiDAR points that lie within vector building footprints.',\n",
       " 'classify_overlap_points': 'Classifies or filters LAS points in regions of overlapping flight lines.',\n",
       " 'clean_vector': 'Removes null features and lines/polygons with fewer than the required number of vertices.',\n",
       " 'clip': 'Extract all the features, or parts of features, that overlap with the features of the clip vector.',\n",
       " 'clip_lidar_to_polygon': 'Clips a LiDAR point cloud to a vector polygon or polygons.',\n",
       " 'clip_raster_to_polygon': 'Clips a raster to a vector polygon.',\n",
       " 'closing': 'A closing is a mathematical morphology operation involving an erosion (min filter) of a dilation (max filter) set.',\n",
       " 'clump': 'Groups cells that form discrete areas, assigning them unique identifiers.',\n",
       " 'compactness_ratio': 'Calculates the compactness ratio (A/P), a measure of shape complexity, for vector polygons.',\n",
       " 'conditional_evaluation': 'Performs a conditional evaluation (if-then-else) operation on a raster.',\n",
       " 'conditioned_latin_hypercube': 'Implements conditioned Latin Hypercube sampling.',\n",
       " 'conservative_smoothing_filter': 'Performs a conservative-smoothing filter on an image.',\n",
       " 'construct_vector_tin': 'Creates a vector triangular irregular network (TIN) for a set of vector points.',\n",
       " 'contours_from_points': 'Creates a contour coverage from a set of input points.',\n",
       " 'contours_from_raster': 'Derives a vector contour coverage from a raster surface.',\n",
       " 'convert_nodata_to_zero': 'Converts nodata values in a raster to zero.',\n",
       " 'convert_raster_format': 'Converts raster data from one format to another.',\n",
       " 'corner_detection': 'Identifies corner patterns in boolean images using hit-and-miss pattern matching.',\n",
       " 'correct_vignetting': 'Corrects the darkening of images towards corners.',\n",
       " 'cos': 'Returns the cosine (cos) of each values in a raster.',\n",
       " 'cosh': 'Returns the hyperbolic cosine (cosh) of each values in a raster.',\n",
       " 'cost_allocation': 'Identifies the source cell to which each grid cell is connected by a least-cost pathway in a cost-distance analysis.',\n",
       " 'cost_distance': 'Performs cost-distance accumulation on a cost surface and a group of source cells.',\n",
       " 'cost_pathway': 'Performs cost-distance pathway analysis using a series of destination grid cells.',\n",
       " 'count_if': 'Counts the number of occurrences of a specified value in a cell-stack of rasters.',\n",
       " 'create_colour_composite': 'Creates a colour-composite image from three bands of multispectral imagery.',\n",
       " 'create_hexagonal_vector_grid': 'Creates a hexagonal vector grid.',\n",
       " 'create_plane': 'Creates a raster image based on the equation for a simple plane.',\n",
       " 'create_rectangular_vector_grid': 'Creates a rectangular vector grid.',\n",
       " 'crispness_index': 'Calculates the Crispness Index, which is used to quantify how crisp (or conversely how fuzzy) a probability image is.',\n",
       " 'cross_tabulation': 'Performs a cross-tabulation on two categorical images.',\n",
       " 'csv_points_to_vector': 'Converts a CSV text file to vector points.',\n",
       " 'cumulative_distribution': 'Converts a raster image to its cumulative distribution function.',\n",
       " 'd8_flow_accumulation': 'Calculates a D8 flow accumulation raster from an input DEM or flow pointer.',\n",
       " 'd8_mass_flux': 'Performs a D8 mass flux calculation.',\n",
       " 'd8_pointer': 'Calculates a D8 flow pointer raster from an input DEM.',\n",
       " 'd_inf_flow_accumulation': 'Calculates a D-infinity flow accumulation raster from an input DEM.',\n",
       " 'd_inf_mass_flux': 'Performs a D-infinity mass flux calculation.',\n",
       " 'd_inf_pointer': 'Calculates a D-infinity flow pointer (flow direction) raster from an input DEM.',\n",
       " 'decrement': 'Decreases the values of each grid cell in an input raster by 1.0 (see also InPlaceSubtract).',\n",
       " 'depth_in_sink': 'Measures the depth of sinks (depressions) in a DEM.',\n",
       " 'dev_from_mean_elev': 'Calculates deviation from mean elevation.',\n",
       " 'diff_from_mean_elev': 'Calculates difference from mean elevation (equivalent to a high-pass filter).',\n",
       " 'diff_of_gaussian_filter': 'Performs a Difference of Gaussian (DoG) filter on an image.',\n",
       " 'difference': 'Outputs the features that occur in one of the two vector inputs but not both, i.e. no overlapping features.',\n",
       " 'direct_decorrelation_stretch': 'Performs a direct decorrelation stretch enhancement on a colour-composite image of multispectral data.',\n",
       " 'directional_relief': 'Calculates relief for cells in an input DEM for a specified direction.',\n",
       " 'dissolve': 'Removes the interior, or shared, boundaries within a vector polygon coverage.',\n",
       " 'distance_to_outlet': 'Calculates the distance of stream grid cells to the channel network outlet cell.',\n",
       " 'diversity_filter': 'Assigns each cell in the output grid the number of different values in a moving window centred on each grid cell in the input raster.',\n",
       " 'divide': 'Performs a division operation on two rasters or a raster and a constant value.',\n",
       " 'downslope_distance_to_stream': 'Measures distance to the nearest downslope stream cell.',\n",
       " 'downslope_flowpath_length': 'Calculates the downslope flowpath length from each cell to basin outlet.',\n",
       " 'downslope_index': 'Calculates the Hjerdt et al. (2004) downslope index.',\n",
       " 'edge_contamination': 'Identifies grid cells within an input DEM that may be impacted by edge contamination for hydrological applications.',\n",
       " 'edge_density': 'Calculates the density of edges, or breaks-in-slope within DEMs.',\n",
       " 'edge_preserving_mean_filter': 'Performs a simple edge-preserving mean filter on an input image.',\n",
       " 'edge_proportion': 'Calculate the proportion of cells in a raster polygon that are edge cells.',\n",
       " 'elev_above_pit': 'Calculate the elevation of each grid cell above the nearest downstream pit cell or grid edge cell.',\n",
       " 'elev_percentile': 'Calculates the elevation percentile raster from a DEM.',\n",
       " 'elev_relative_to_min_max': 'Calculates the elevation of a location relative to the minimum and maximum elevations in a DEM.',\n",
       " 'elev_relative_to_watershed_min_max': 'Calculates the elevation of a location relative to the minimum and maximum elevations in a watershed.',\n",
       " 'elevation_above_stream': 'Calculates the elevation of cells above the nearest downslope stream cell.',\n",
       " 'elevation_above_stream_euclidean': 'Calculates the elevation of cells above the nearest (Euclidean distance) stream cell.',\n",
       " 'eliminate_coincident_points': 'Removes any coincident, or nearly coincident, points from a vector points file.',\n",
       " 'elongation_ratio': 'Calculates the elongation ratio for vector polygons.',\n",
       " 'embankment_mapping': 'Maps and/or removes road embankments from an input fine-resolution DEM.',\n",
       " 'emboss_filter': 'Performs an emboss filter on an image, similar to a hillshade operation.',\n",
       " 'equal_to': 'Performs a equal-to comparison operation on two rasters or a raster and a constant value.',\n",
       " 'erase': 'Removes all the features, or parts of features, that overlap with the features of the erase vector polygon.',\n",
       " 'erase_polygon_from_lidar': 'Erases (cuts out) a vector polygon or polygons from a LiDAR point cloud.',\n",
       " 'erase_polygon_from_raster': 'Erases (cuts out) a vector polygon from a raster.',\n",
       " 'euclidean_allocation': 'Assigns grid cells in the output raster the value of the nearest target cell in the input image, measured by the Shih and Wu (2004) Euclidean distance transform.',\n",
       " 'euclidean_distance': 'Calculates the Shih and Wu (2004) Euclidean distance transform.',\n",
       " 'exp': 'Returns the exponential (base e) of values in a raster.',\n",
       " 'exp2': 'Returns the exponential (base 2) of values in a raster.',\n",
       " 'export_table_to_csv': 'Exports an attribute table to a CSV text file.',\n",
       " 'exposure_towards_wind_flux': 'Evaluates hydrologic connectivity within a DEM',\n",
       " 'extend_vector_lines': 'Extends vector lines by a specified distance.',\n",
       " 'extract_nodes': 'Converts vector lines or polygons into vertex points.',\n",
       " 'extract_raster_values_at_points': 'Extracts the values of raster(s) at vector point locations.',\n",
       " 'extract_streams': 'Extracts stream grid cells from a flow accumulation raster.',\n",
       " 'extract_valleys': 'Identifies potential valley bottom grid cells based on local topolography alone.',\n",
       " 'fd8_flow_accumulation': 'Calculates an FD8 flow accumulation raster from an input DEM.',\n",
       " 'fd8_pointer': 'Calculates an FD8 flow pointer raster from an input DEM.',\n",
       " 'farthest_channel_head': 'Calculates the distance to the furthest upstream channel head for each stream cell.',\n",
       " 'fast_almost_gaussian_filter': 'Performs a fast approximate Gaussian filter on an image.',\n",
       " 'feature_preserving_smoothing': 'Reduces short-scale variation in an input DEM using a modified Sun et al. (2007) algorithm.',\n",
       " 'fetch_analysis': 'Performs an analysis of fetch or upwind distance to an obstacle.',\n",
       " 'fill_burn': 'Burns streams into a DEM using the FillBurn (Saunders, 1999) method.',\n",
       " 'fill_depressions': 'Fills all of the depressions in a DEM. Depression breaching should be preferred in most cases.',\n",
       " 'fill_depressions_planchon_and_darboux': 'Fills all of the depressions in a DEM using the Planchon and Darboux (2002) method.',\n",
       " 'fill_depressions_wang_and_liu': 'Fills all of the depressions in a DEM using the Wang and Liu (2006) method. Depression breaching should be preferred in most cases.',\n",
       " 'fill_missing_data': 'Fills NoData holes in a DEM.',\n",
       " 'fill_single_cell_pits': 'Raises pit cells to the elevation of their lowest neighbour.',\n",
       " 'filter_lidar_classes': 'Removes points in a LAS file with certain specified class values.',\n",
       " 'filter_lidar_scan_angles': 'Removes points in a LAS file with scan angles greater than a threshold.',\n",
       " 'filter_raster_features_by_area': 'Removes small-area features from a raster.',\n",
       " 'find_flightline_edge_points': \"Identifies points along a flightline's edge in a LAS file.\",\n",
       " 'find_lowest_or_highest_points': 'Locates the lowest and/or highest valued cells in a raster.',\n",
       " 'find_main_stem': 'Finds the main stem, based on stream lengths, of each stream network.',\n",
       " 'find_no_flow_cells': 'Finds grid cells with no downslope neighbours.',\n",
       " 'find_parallel_flow': 'Finds areas of parallel flow in D8 flow direction rasters.',\n",
       " 'find_patch_or_class_edge_cells': 'Finds all cells located on the edge of patch or class features.',\n",
       " 'find_ridges': 'Identifies potential ridge and peak grid cells.',\n",
       " 'flatten_lakes': 'Flattens lake polygons in a raster DEM.',\n",
       " 'flightline_overlap': 'Reads a LiDAR (LAS) point file and outputs a raster containing the number of overlapping flight-lines in each grid cell.',\n",
       " 'flip_image': 'Reflects an image in the vertical or horizontal axis.',\n",
       " 'flood_order': 'Assigns each DEM grid cell its order in the sequence of inundations that are encountered during a search starting from the edges, moving inward at increasing elevations.',\n",
       " 'floor': 'Returns the largest (closest to positive infinity) value that is less than or equal to the values in a raster.',\n",
       " 'flow_accumulation_full_workflow': 'Resolves all of the depressions in a DEM, outputting a breached DEM, an aspect-aligned non-divergent flow pointer, and a flow accumulation raster.',\n",
       " 'flow_length_diff': 'Calculates the local maximum absolute difference in downslope flowpath length, useful in mapping drainage divides and ridges.',\n",
       " 'gamma_correction': 'Performs a gamma correction on an input images.',\n",
       " 'gaussian_contrast_stretch': 'Performs a Gaussian contrast stretch on input images.',\n",
       " 'gaussian_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'gaussian_filter': 'Performs a Gaussian filter on an image.',\n",
       " 'gaussian_scale_space': 'Uses the fast Gaussian approximation algorithm to produce scaled land-surface parameter measurements from an input DEM.',\n",
       " 'geomorphons': 'Computes geomorphon patterns.',\n",
       " 'greater_than': 'Performs a greater-than comparison operation on two rasters or a raster and a constant value.',\n",
       " 'hack_stream_order': 'Assigns the Hack stream order to each tributary in a stream network.',\n",
       " 'heat_map': 'Calculates a heat map, or kernel density estimation (KDE), for an input point set.',\n",
       " 'height_above_ground': 'Normalizes a LiDAR point cloud, providing the height above the nearest ground-classified point.',\n",
       " 'high_pass_bilateral_filter': 'Performs a high-pass bilateral filter, by differencing an input image by the bilateral filter by Tomasi and Manduchi (1998).',\n",
       " 'high_pass_filter': 'Performs a high-pass filter on an input image.',\n",
       " 'high_pass_median_filter': 'Performs a high pass median filter on an input image.',\n",
       " 'highest_position': 'Identifies the stack position of the maximum value within a raster stack on a cell-by-cell basis.',\n",
       " 'hillshade': 'Calculates a hillshade raster from an input DEM.',\n",
       " 'hillslopes': 'Identifies the individual hillslopes draining to each link in a stream network.',\n",
       " 'histogram_equalization': 'Performs a histogram equalization contrast enhancement on an image.',\n",
       " 'histogram_matching': 'Alters the statistical distribution of a raster image matching it to a specified PDF.',\n",
       " 'histogram_matching_two_images': 'Alters the cumulative distribution function of a raster image to that of another image.',\n",
       " 'hole_proportion': \"Calculates the proportion of the total area of a polygon's holes relative to the area of the polygon's hull.\",\n",
       " 'horizon_angle': 'Calculates horizon angle (maximum upwind slope) for each grid cell in an input DEM.',\n",
       " 'horton_stream_order': 'Assigns the Horton stream order to each tributary in a stream network.',\n",
       " 'hypsometric_analysis': 'Calculates a hypsometric curve for one or more DEMs.',\n",
       " 'hypsometrically_tinted_hillshade': 'Creates an colour shaded relief image from an input DEM.',\n",
       " 'idw_interpolation': 'Interpolates vector points into a raster surface using an inverse-distance weighted scheme.',\n",
       " 'ihs_to_rgb': 'Converts intensity, hue, and saturation (IHS) images into red, green, and blue (RGB) images.',\n",
       " 'image_autocorrelation': \"Performs Moran's I analysis on two or more input images.\",\n",
       " 'image_correlation': 'Performs image correlation on two or more input images.',\n",
       " 'image_correlation_neighbourhood_analysis': 'Performs image correlation on two input images neighbourhood search windows.',\n",
       " 'image_regression': 'Performs image regression analysis on two input images.',\n",
       " 'image_stack_profile': 'Plots an image stack profile (i.e. signature) for a set of points and multispectral images.',\n",
       " 'impoundment_size_index': 'Calculates the impoundment size resulting from damming a DEM.',\n",
       " 'in_place_add': 'Performs an in-place addition operation (input1 += input2).',\n",
       " 'in_place_divide': 'Performs an in-place division operation (input1 /= input2).',\n",
       " 'in_place_multiply': 'Performs an in-place multiplication operation (input1 *= input2).',\n",
       " 'in_place_subtract': 'Performs an in-place subtraction operation (input1 -= input2).',\n",
       " 'increment': 'Increases the values of each grid cell in an input raster by 1.0. (see also InPlaceAdd)',\n",
       " 'individual_tree_detection': 'Identifies points in a LiDAR point cloud that are associated with the tops of individual trees.',\n",
       " 'insert_dams': 'Calculates the impoundment size resulting from damming a DEM.',\n",
       " 'install_wb_extension': 'Use to install a Whitebox extension product.',\n",
       " 'integer_division': 'Performs an integer division operation on two rasters or a raster and a constant value.',\n",
       " 'integral_image': 'Transforms an input image (summed area table) into its integral image equivalent.',\n",
       " 'intersect': 'Identifies the parts of features in common between two input vector layers.',\n",
       " 'is_no_data': 'Identifies NoData valued pixels in an image.',\n",
       " 'isobasins': 'Divides a landscape into nearly equal sized drainage basins (i.e. watersheds).',\n",
       " 'jenson_snap_pour_points': 'Moves outlet points used to specify points of interest in a watershedding operation to the nearest stream cell.',\n",
       " 'join_tables': \"Merge a vector's attribute table with another table based on a common field.\",\n",
       " 'k_means_clustering': 'Performs a k-means clustering operation on a multi-spectral dataset.',\n",
       " 'k_nearest_mean_filter': 'A k-nearest mean filter is a type of edge-preserving smoothing filter.',\n",
       " 'kappa_index': 'Performs a kappa index of agreement (KIA) analysis on two categorical raster files.',\n",
       " 'ks_test_for_normality': 'Evaluates whether the values in a raster are normally distributed.',\n",
       " 'laplacian_filter': 'Performs a Laplacian filter on an image.',\n",
       " 'laplacian_of_gaussian_filter': 'Performs a Laplacian-of-Gaussian (LoG) filter on an image.',\n",
       " 'las_to_ascii': 'Converts one or more LAS files into ASCII text files.',\n",
       " 'las_to_multipoint_shapefile': 'Converts one or more LAS files into MultipointZ vector Shapefiles. When the input parameter is not specified, the tool grids all LAS files contained within the working directory.',\n",
       " 'las_to_shapefile': 'Converts one or more LAS files into a vector Shapefile of POINT ShapeType.',\n",
       " 'las_to_zlidar': 'Converts one or more LAS files into the zlidar compressed LiDAR data format.',\n",
       " 'launch_wb_runner': 'Opens the Whitebox Runner application.',\n",
       " 'layer_footprint': 'Creates a vector polygon footprint of the area covered by a raster grid or vector layer.',\n",
       " 'lee_sigma_filter': 'Performs a Lee (Sigma) smoothing filter on an image.',\n",
       " 'length_of_upstream_channels': 'Calculates the total length of channels upstream.',\n",
       " 'less_than': 'Performs a less-than comparison operation on two rasters or a raster and a constant value.',\n",
       " 'lidar_block_maximum': 'Creates a block-maximum raster from an input LAS file. When the input/output parameters are not specified, the tool grids all LAS files contained within the working directory.',\n",
       " 'lidar_block_minimum': 'Creates a block-minimum raster from an input LAS file. When the input/output parameters are not specified, the tool grids all LAS files contained within the working directory.',\n",
       " 'lidar_classify_subset': 'Classifies the values in one LiDAR point cloud that correspond with points in a subset cloud.',\n",
       " 'lidar_colourize': 'Adds the red-green-blue colour fields of a LiDAR (LAS) file based on an input image.',\n",
       " 'lidar_digital_surface_model': 'Creates a top-surface digital surface model (DSM) from a LiDAR point cloud.',\n",
       " 'lidar_elevation_slice': 'Outputs all of the points within a LiDAR (LAS) point file that lie between a specified elevation range.',\n",
       " 'lidar_ground_point_filter': 'Identifies ground points within LiDAR dataset using a slope-based method.',\n",
       " 'lidar_hex_binning': 'Hex-bins a set of LiDAR points.',\n",
       " 'lidar_hillshade': 'Calculates a hillshade value for points within a LAS file and stores these data in the RGB field.',\n",
       " 'lidar_histogram': 'Creates a histogram of LiDAR data.',\n",
       " 'lidar_idw_interpolation': 'Interpolates LAS files using an inverse-distance weighted (IDW) scheme. When the input/output parameters are not specified, the tool interpolates all LAS files contained within the working directory.',\n",
       " 'lidar_info': 'Prints information about a LiDAR (LAS) dataset, including header, point return frequency, and classification data and information about the variable length records (VLRs) and geokeys.',\n",
       " 'lidar_join': 'Joins multiple LiDAR (LAS) files into a single LAS file.',\n",
       " 'lidar_kappa_index': 'Performs a kappa index of agreement (KIA) analysis on the classifications of two LAS files.',\n",
       " 'lidar_nearest_neighbour_gridding': 'Grids LiDAR files using nearest-neighbour scheme. When the input/output parameters are not specified, the tool grids all LAS files contained within the working directory.',\n",
       " 'lidar_point_density': 'Calculates the spatial pattern of point density for a LiDAR data set. When the input/output parameters are not specified, the tool grids all LAS files contained within the working directory.',\n",
       " 'lidar_point_stats': 'Creates several rasters summarizing the distribution of LAS point data. When the input/output parameters are not specified, the tool works on all LAS files contained within the working directory.',\n",
       " 'lidar_ransac_planes': 'Performs a RANSAC analysis to identify points within a LiDAR point cloud that belong to linear planes.',\n",
       " 'lidar_rbf_interpolation': 'Interpolates LAS files using a radial basis function (RBF) scheme. When the input/output parameters are not specified, the tool interpolates all LAS files contained within the working directory.',\n",
       " 'lidar_remove_duplicates': 'Removes duplicate points from a LiDAR data set.',\n",
       " 'lidar_remove_outliers': 'Removes outliers (high and low points) in a LiDAR point cloud.',\n",
       " 'lidar_rooftop_analysis': 'Identifies roof segments in a LiDAR point cloud.',\n",
       " 'lidar_segmentation': 'Segments a LiDAR point cloud based on differences in the orientation of fitted planar surfaces and point proximity.',\n",
       " 'lidar_segmentation_based_filter': 'Identifies ground points within LiDAR point clouds using a segmentation based approach.',\n",
       " 'lidar_shift': 'Shifts the x,y,z coordinates of a LiDAR file.',\n",
       " 'lidar_tin_gridding': 'Creates a raster grid based on a Delaunay triangular irregular network (TIN) fitted to LiDAR points.',\n",
       " 'lidar_thin': 'Thins a LiDAR point cloud, reducing point density.',\n",
       " 'lidar_thin_high_density': 'Thins points from high density areas within a LiDAR point cloud.',\n",
       " 'lidar_tile': 'Tiles a LiDAR LAS file into multiple LAS files.',\n",
       " 'lidar_tile_footprint': 'Creates a vector polygon of the convex hull of a LiDAR point cloud. When the input/output parameters are not specified, the tool works with all LAS files contained within the working directory.',\n",
       " 'lidar_tophat_transform': 'Performs a white top-hat transform on a Lidar dataset; as an estimate of height above ground, this is useful for modelling the vegetation canopy.',\n",
       " 'line_detection_filter': 'Performs a line-detection filter on an image.',\n",
       " 'line_intersections': 'Identifies points where the features of two vector line layers intersect.',\n",
       " 'line_thinning': 'Performs line thinning a on Boolean raster image; intended to be used with the RemoveSpurs tool.',\n",
       " 'linearity_index': 'Calculates the linearity index for vector polygons.',\n",
       " 'lines_to_polygons': 'Converts vector polylines to polygons.',\n",
       " 'list_unique_values': \"Lists the unique values contained in a field within a vector's attribute table.\",\n",
       " 'list_unique_values_raster': \"Lists the unique values contained in a field within a vector's attribute table.\",\n",
       " 'ln': 'Returns the natural logarithm of values in a raster.',\n",
       " 'local_quadratic_regression': 'An implementation of the constrained quadratic regression algorithm using a flexible window size described in Wood (1996).',\n",
       " 'log10': 'Returns the base-10 logarithm of values in a raster.',\n",
       " 'log2': 'Returns the base-2 logarithm of values in a raster.',\n",
       " 'long_profile': 'Plots the stream longitudinal profiles for one or more rivers.',\n",
       " 'long_profile_from_points': 'Plots the longitudinal profiles from flow-paths initiating from a set of vector points.',\n",
       " 'longest_flowpath': 'Delineates the longest flowpaths for a group of subbasins or watersheds.',\n",
       " 'lowest_position': 'Identifies the stack position of the minimum value within a raster stack on a cell-by-cell basis.',\n",
       " 'md_inf_flow_accumulation': 'Calculates an FD8 flow accumulation raster from an input DEM.',\n",
       " 'majority_filter': 'Assigns each cell in the output grid the most frequently occurring value (mode) in a moving window centred on each grid cell in the input raster.',\n",
       " 'map_off_terrain_objects': 'Maps off-terrain objects in a digital elevation model (DEM).',\n",
       " 'max': 'Performs a MAX operation on two rasters or a raster and a constant value.',\n",
       " 'max_absolute_overlay': 'Evaluates the maximum absolute value for each grid cell from a stack of input rasters.',\n",
       " 'max_anisotropy_dev': 'Calculates the maximum anisotropy (directionality) in elevation deviation over a range of spatial scales.',\n",
       " 'max_anisotropy_dev_signature': 'Calculates the anisotropy in deviation from mean for points over a range of spatial scales.',\n",
       " 'max_branch_length': \"Lindsay and Seibert's (2013) branch length index is used to map drainage divides or ridge lines.\",\n",
       " 'max_difference_from_mean': 'Calculates the maximum difference from mean elevation over a range of spatial scales.',\n",
       " 'max_downslope_elev_change': 'Calculates the maximum downslope change in elevation between a grid cell and its eight downslope neighbors.',\n",
       " 'max_elev_dev_signature': 'Calculates the maximum elevation deviation over a range of spatial scales and for a set of points.',\n",
       " 'max_elevation_deviation': 'Calculates the maximum elevation deviation over a range of spatial scales.',\n",
       " 'max_overlay': 'Evaluates the maximum value for each grid cell from a stack of input rasters.',\n",
       " 'max_upslope_elev_change': 'Calculates the maximum upslope change in elevation between a grid cell and its eight downslope neighbors.',\n",
       " 'max_upslope_flowpath_length': 'Measures the maximum length of all upslope flowpaths draining each grid cell.',\n",
       " 'max_upslope_value': 'Calculates the maximum upslope value from an input values raster along flowpaths.',\n",
       " 'maximal_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'maximum_filter': 'Assigns each cell in the output grid the maximum value in a moving window centred on each grid cell in the input raster.',\n",
       " 'mean_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'mean_filter': 'Performs a mean filter (low-pass filter) on an input image.',\n",
       " 'median_filter': 'Performs a median filter on an input image.',\n",
       " 'medoid': 'Calculates the medoid for a series of vector features contained in a shapefile.',\n",
       " 'merge_line_segments': 'Merges vector line segments into larger features.',\n",
       " 'merge_table_with_csv': \"Merge a vector's attribute table with a table contained within a CSV text file.\",\n",
       " 'merge_vectors': 'Combines two or more input vectors of the same ShapeType creating a single, new output vector.',\n",
       " 'min': 'Performs a MIN operation on two rasters or a raster and a constant value.',\n",
       " 'min_absolute_overlay': 'Evaluates the minimum absolute value for each grid cell from a stack of input rasters.',\n",
       " 'min_downslope_elev_change': 'Calculates the minimum downslope change in elevation between a grid cell and its eight downslope neighbors.',\n",
       " 'min_max_contrast_stretch': 'Performs a min-max contrast stretch on an input greytone image.',\n",
       " 'min_overlay': 'Evaluates the minimum value for each grid cell from a stack of input rasters.',\n",
       " 'minimal_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'minimum_bounding_box': 'Creates a vector minimum bounding rectangle around vector features.',\n",
       " 'minimum_bounding_circle': 'Delineates the minimum bounding circle (i.e. smallest enclosing circle) for a group of vectors.',\n",
       " 'minimum_bounding_envelope': 'Creates a vector axis-aligned minimum bounding rectangle (envelope) around vector features.',\n",
       " 'minimum_convex_hull': 'Creates a vector convex polygon around vector features.',\n",
       " 'minimum_filter': 'Assigns each cell in the output grid the minimum value in a moving window centred on each grid cell in the input raster.',\n",
       " 'modified_k_means_clustering': 'Performs a modified k-means clustering operation on a multi-spectral dataset.',\n",
       " 'modify_no_data_value': 'Modifies nodata values in a raster.',\n",
       " 'modulo': 'Performs a modulo operation on two rasters or a raster and a constant value.',\n",
       " 'mosaic': 'Mosaics two or more images together.',\n",
       " 'mosaic_with_feathering': 'Mosaics two images together using a feathering technique in overlapping areas to reduce edge-effects.',\n",
       " 'multi_part_to_single_part': 'Converts a vector file containing multi-part features into a vector containing only single-part features.',\n",
       " 'multidirectional_hillshade': 'Calculates a multi-direction hillshade raster from an input DEM.',\n",
       " 'multiply': 'Performs a multiplication operation on two rasters or a raster and a constant value.',\n",
       " 'multiply_overlay': 'Calculates the sum for each grid cell from a group of raster images.',\n",
       " 'multiscale_elevation_percentile': 'Calculates surface roughness over a range of spatial scales.',\n",
       " 'multiscale_roughness': 'Calculates surface roughness over a range of spatial scales.',\n",
       " 'multiscale_roughness_signature': 'Calculates the surface roughness for points over a range of spatial scales.',\n",
       " 'multiscale_std_dev_normals': 'Calculates surface roughness over a range of spatial scales.',\n",
       " 'multiscale_std_dev_normals_signature': 'Calculates the surface roughness for points over a range of spatial scales.',\n",
       " 'multiscale_topographic_position_image': 'Creates a multiscale topographic position image from three DEVmax rasters of differing spatial scale ranges.',\n",
       " 'narrowness_index': 'Calculates the narrowness of raster polygons.',\n",
       " 'natural_neighbour_interpolation': \"Creates a raster grid based on Sibson's natural neighbour method.\",\n",
       " 'nearest_neighbour_gridding': 'Creates a raster grid based on a set of vector points and assigns grid values using the nearest neighbour.',\n",
       " 'negate': 'Changes the sign of values in a raster or the 0-1 values of a Boolean raster.',\n",
       " 'new_raster_from_base': 'Creates a new raster using a base image.',\n",
       " 'normal_vectors': 'Calculates normal vectors for points within a LAS file and stores these data (XYZ vector components) in the RGB field.',\n",
       " 'normalize_lidar': 'Normalizes a LiDAR point cloud.',\n",
       " 'normalized_difference_index': 'Calculate a normalized-difference index (NDI) from two bands of multispectral image data.',\n",
       " 'not': 'Performs a logical NOT operator on two Boolean raster images.',\n",
       " 'not_equal_to': 'Performs a not-equal-to comparison operation on two rasters or a raster and a constant value.',\n",
       " 'num_downslope_neighbours': 'Calculates the number of downslope neighbours to each grid cell in a DEM.',\n",
       " 'num_inflowing_neighbours': 'Computes the number of inflowing neighbours to each cell in an input DEM based on the D8 algorithm.',\n",
       " 'num_upslope_neighbours': 'Calculates the number of upslope neighbours to each grid cell in a DEM.',\n",
       " 'olympic_filter': 'Performs an olympic smoothing filter on an image.',\n",
       " 'opening': 'An opening is a mathematical morphology operation involving a dilation (max filter) of an erosion (min filter) set.',\n",
       " 'or': 'Performs a logical OR operator on two Boolean raster images.',\n",
       " 'paired_sample_t_test': 'Performs a 2-sample K-S test for significant differences on two input rasters.',\n",
       " 'panchromatic_sharpening': 'Increases the spatial resolution of image data by combining multispectral bands with panchromatic data.',\n",
       " 'patch_orientation': 'Calculates the orientation of vector polygons.',\n",
       " 'pennock_landform_class': 'Classifies hillslope zones based on slope, profile curvature, and plan curvature.',\n",
       " 'percent_elev_range': 'Calculates percent of elevation range from a DEM.',\n",
       " 'percent_equal_to': 'Calculates the percentage of a raster stack that have cell values equal to an input on a cell-by-cell basis.',\n",
       " 'percent_greater_than': 'Calculates the percentage of a raster stack that have cell values greater than an input on a cell-by-cell basis.',\n",
       " 'percent_less_than': 'Calculates the percentage of a raster stack that have cell values less than an input on a cell-by-cell basis.',\n",
       " 'percentage_contrast_stretch': 'Performs a percentage linear contrast stretch on input images.',\n",
       " 'percentile_filter': 'Performs a percentile filter on an input image.',\n",
       " 'perimeter_area_ratio': 'Calculates the perimeter-area ratio of vector polygons.',\n",
       " 'pick_from_list': 'Outputs the value from a raster stack specified by a position raster.',\n",
       " 'plan_curvature': 'Calculates a plan (contour) curvature raster from an input DEM.',\n",
       " 'polygon_area': 'Calculates the area of vector polygons.',\n",
       " 'polygon_long_axis': 'Used to map the long axis of polygon features.',\n",
       " 'polygon_perimeter': 'Calculates the perimeter of vector polygons.',\n",
       " 'polygon_short_axis': 'Used to map the short axis of polygon features.',\n",
       " 'polygonize': 'Creates a polygon layer from two or more intersecting line features contained in one or more input vector line files.',\n",
       " 'polygons_to_lines': 'Converts vector polygons to polylines.',\n",
       " 'power': 'Raises the values in grid cells of one rasters, or a constant value, by values in another raster or constant value.',\n",
       " 'prewitt_filter': 'Performs a Prewitt edge-detection filter on an image.',\n",
       " 'principal_component_analysis': 'Performs a principal component analysis (PCA) on a multi-spectral dataset.',\n",
       " 'print_geo_tiff_tags': 'Prints the tags within a GeoTIFF.',\n",
       " 'profile': 'Plots profiles from digital surface models.',\n",
       " 'profile_curvature': 'Calculates a profile curvature raster from an input DEM.',\n",
       " 'qin_flow_accumulation': 'Calculates Qin et al. (2007) flow accumulation.',\n",
       " 'quantiles': 'Transforms raster values into quantiles.',\n",
       " 'quinn_flow_accumulation': 'Calculates Quinn et al. (1995) flow accumulation.',\n",
       " 'radial_basis_function_interpolation': 'Interpolates vector points into a raster surface using a radial basis function scheme.',\n",
       " 'radius_of_gyration': \"Calculates the distance of cells from their polygon's centroid.\",\n",
       " 'raise_walls': 'Raises walls in a DEM along a line or around a polygon, e.g. a watershed.',\n",
       " 'random_field': 'Creates an image containing random values.',\n",
       " 'random_sample': 'Creates an image containing randomly located sample grid cells with unique IDs.',\n",
       " 'range_filter': 'Assigns each cell in the output grid the range of values in a moving window centred on each grid cell in the input raster.',\n",
       " 'raster_area': 'Calculates the area of polygons or classes within a raster image.',\n",
       " 'raster_calculator': 'Performs a complex mathematical operations on one or more input raster images on a cell-to-cell basis.',\n",
       " 'raster_cell_assignment': 'Assign row or column number to cells.',\n",
       " 'raster_histogram': 'Creates a histogram from raster values.',\n",
       " 'raster_perimeter': 'Calculates the perimeters of polygons or classes within a raster image.',\n",
       " 'raster_streams_to_vector': 'Converts a raster stream file into a vector file.',\n",
       " 'raster_summary_stats': 'Measures a rasters min, max, average, standard deviation, num. non-nodata cells, and total.',\n",
       " 'raster_to_vector_lines': 'Converts a raster lines features into a vector of the POLYLINE shapetype',\n",
       " 'raster_to_vector_points': 'Converts a raster dataset to a vector of the POINT shapetype.',\n",
       " 'raster_to_vector_polygons': 'Converts a raster dataset to a vector of the POLYGON shapetype.',\n",
       " 'rasterize_streams': 'Rasterizes vector streams based on Lindsay (2016) method.',\n",
       " 'reciprocal': 'Returns the reciprocal (i.e. 1 / z) of values in a raster.',\n",
       " 'reclass': 'Reclassifies the values in a raster image.',\n",
       " 'reclass_equal_interval': 'Reclassifies the values in a raster image based on equal-ranges.',\n",
       " 'reclass_from_file': 'Reclassifies the values in a raster image using reclass ranges in a text file.',\n",
       " 'reinitialize_attribute_table': \"Reinitializes a vector's attribute table deleting all fields but the feature ID (FID).\",\n",
       " 'related_circumscribing_circle': 'Calculates the related circumscribing circle of vector polygons.',\n",
       " 'relative_aspect': 'Calculates relative aspect (relative to a user-specified direction) from an input DEM.',\n",
       " 'relative_topographic_position': 'Calculates the relative topographic position index from a DEM.',\n",
       " 'remove_off_terrain_objects': 'Removes off-terrain objects from a raster digital elevation model (DEM).',\n",
       " 'remove_polygon_holes': 'Removes holes within the features of a vector polygon file.',\n",
       " 'remove_short_streams': 'Removes short first-order streams from a stream network.',\n",
       " 'remove_spurs': 'Removes the spurs (pruning operation) from a Boolean line image; intended to be used on the output of the LineThinning tool.',\n",
       " 'repair_stream_vector_topology': 'This tool resolves topological errors and inconsistencies associated with digitized vector streams.',\n",
       " 'resample': 'Resamples one or more input images into a destination image.',\n",
       " 'rescale_value_range': 'Performs a min-max contrast stretch on an input greytone image.',\n",
       " 'rgb_to_ihs': 'Converts red, green, and blue (RGB) images into intensity, hue, and saturation (IHS) images.',\n",
       " 'rho8_flow_accumulation': 'Calculates Fairfield and Leymarie (1991) flow accumulation.',\n",
       " 'rho8_pointer': 'Calculates a stochastic Rho8 flow pointer raster from an input DEM.',\n",
       " 'roberts_cross_filter': \"Performs a Robert's cross edge-detection filter on an image.\",\n",
       " 'root_mean_square_error': 'Calculates the RMSE and other accuracy statistics.',\n",
       " 'round': 'Rounds the values in an input raster to the nearest integer value.',\n",
       " 'ruggedness_index': \"Calculates the Riley et al.'s (1999) terrain ruggedness index from an input DEM.\",\n",
       " 'scharr_filter': 'Performs a Scharr edge-detection filter on an image.',\n",
       " 'sediment_transport_index': 'Calculates the sediment transport index.',\n",
       " 'select_tiles_by_polygon': 'Copies LiDAR tiles overlapping with a polygon into an output directory.',\n",
       " 'set_nodata_value': 'Assign the NoData value for an input image.',\n",
       " 'shape_complexity_index': 'Calculates overall polygon shape complexity or irregularity.',\n",
       " 'shape_complexity_index_raster': 'Calculates the complexity of raster polygons or classes.',\n",
       " 'shreve_stream_magnitude': 'Assigns the Shreve stream magnitude to each link in a stream network.',\n",
       " 'sigmoidal_contrast_stretch': 'Performs a sigmoidal contrast stretch on input images.',\n",
       " 'sin': 'Returns the sine (sin) of each values in a raster.',\n",
       " 'single_part_to_multi_part': 'Converts a vector file containing multi-part features into a vector containing only single-part features.',\n",
       " 'sinh': 'Returns the hyperbolic sine (sinh) of each values in a raster.',\n",
       " 'sink': 'Identifies the depressions in a DEM, giving each feature a unique identifier.',\n",
       " 'slope': 'Calculates a slope raster from an input DEM.',\n",
       " 'slope_vs_elevation_plot': 'Creates a slope vs. elevation plot for one or more DEMs.',\n",
       " 'smooth_vectors': 'Smooths a vector coverage of either a POLYLINE or POLYGON base ShapeType.',\n",
       " 'snap_pour_points': 'Moves outlet points used to specify points of interest in a watershedding operation to the cell with the highest flow accumulation in its neighbourhood.',\n",
       " 'sobel_filter': 'Performs a Sobel edge-detection filter on an image.',\n",
       " 'spherical_std_dev_of_normals': 'Calculates the spherical standard deviation of surface normals for a DEM.',\n",
       " 'split_colour_composite': 'Splits an RGB colour composite image into separate multispectral images.',\n",
       " 'split_vector_lines': 'Used to split a vector line coverage into even-lengthed segments.',\n",
       " 'split_with_lines': 'Splits the lines or polygons in one layer using the lines in another layer.',\n",
       " 'square': 'Squares the values in a raster.',\n",
       " 'square_root': 'Returns the square root of the values in a raster.',\n",
       " 'standard_deviation_contrast_stretch': 'Performs a standard-deviation contrast stretch on input images.',\n",
       " 'standard_deviation_filter': 'Assigns each cell in the output grid the standard deviation of values in a moving window centred on each grid cell in the input raster.',\n",
       " 'standard_deviation_of_slope': 'Calculates the standard deviation of slope from an input DEM.',\n",
       " 'stochastic_depression_analysis': 'Performs a stochastic analysis of depressions within a DEM.',\n",
       " 'strahler_order_basins': 'Identifies Strahler-order basins from an input stream network.',\n",
       " 'strahler_stream_order': 'Assigns the Strahler stream order to each link in a stream network.',\n",
       " 'stream_link_class': 'Identifies the exterior/interior links and nodes in a stream network.',\n",
       " 'stream_link_identifier': 'Assigns a unique identifier to each link in a stream network.',\n",
       " 'stream_link_length': 'Estimates the length of each link (or tributary) in a stream network.',\n",
       " 'stream_link_slope': 'Estimates the average slope of each link (or tributary) in a stream network.',\n",
       " 'stream_power_index': 'Calculates the relative stream power index.',\n",
       " 'stream_slope_continuous': 'Estimates the slope of each grid cell in a stream network.',\n",
       " 'subbasins': 'Identifies the catchments, or sub-basin, draining to each link in a stream network.',\n",
       " 'subtract': 'Performs a differencing operation on two rasters or a raster and a constant value.',\n",
       " 'sum_overlay': 'Calculates the sum for each grid cell from a group of raster images.',\n",
       " 'surface_area_ratio': 'Calculates a the surface area ratio of each grid cell in an input DEM.',\n",
       " 'symmetrical_difference': 'Outputs the features that occur in one of the two vector inputs but not both, i.e. no overlapping features.',\n",
       " 'tin_gridding': 'Creates a raster grid based on a triangular irregular network (TIN) fitted to vector points.',\n",
       " 'tan': 'Returns the tangent (tan) of each values in a raster.',\n",
       " 'tangential_curvature': 'Calculates a tangential curvature raster from an input DEM.',\n",
       " 'tanh': 'Returns the hyperbolic tangent (tanh) of each values in a raster.',\n",
       " 'thicken_raster_line': 'Thickens single-cell wide lines within a raster image.',\n",
       " 'time_in_daylight': 'Calculates the proportion of time a location is not within an area of shadow.',\n",
       " 'to_degrees': 'Converts a raster from radians to degrees.',\n",
       " 'to_radians': 'Converts a raster from degrees to radians.',\n",
       " 'tophat_transform': 'Performs either a white or black top-hat transform on an input image.',\n",
       " 'topological_stream_order': 'Assigns each link in a stream network its topological order.',\n",
       " 'total_curvature': 'Calculates a total curvature raster from an input DEM.',\n",
       " 'total_filter': 'Performs a total filter on an input image.',\n",
       " 'trace_downslope_flowpaths': 'Traces downslope flowpaths from one or more target sites (i.e. seed points).',\n",
       " 'travelling_salesman_problem': 'Finds approximate solutions to travelling salesman problems, the goal of which is to identify the shortest route connecting a set of locations.',\n",
       " 'trend_surface': 'Estimates the trend surface of an input raster file.',\n",
       " 'trend_surface_vector_points': 'Estimates a trend surface from vector points.',\n",
       " 'tributary_identifier': 'Assigns a unique identifier to each tributary in a stream network.',\n",
       " 'truncate': 'Truncates the values in a raster to the desired number of decimal places.',\n",
       " 'turning_bands_simulation': 'Creates an image containing random values based on a turning-bands simulation.',\n",
       " 'two_sample_ks_test': 'Performs a 2-sample K-S test for significant differences on two input rasters.',\n",
       " 'union': 'Splits vector layers at their overlaps, creating a layer containing all the portions from both input and overlay layers.',\n",
       " 'unnest_basins': 'Extract whole watersheds for a set of outlet points.',\n",
       " 'unsharp_masking': 'An image sharpening technique that enhances edges.',\n",
       " 'update_nodata_cells': 'Replaces the NoData values in an input raster with the corresponding values contained in a second update layer.',\n",
       " 'upslope_depression_storage': 'Estimates the average upslope depression storage depth.',\n",
       " 'user_defined_weights_filter': 'Performs a user-defined weights filter on an image.',\n",
       " 'vector_hex_binning': 'Hex-bins a set of vector points.',\n",
       " 'vector_lines_to_raster': 'Converts a vector containing polylines into a raster.',\n",
       " 'vector_points_to_raster': 'Converts a vector containing points into a raster.',\n",
       " 'vector_polygons_to_raster': 'Converts a vector containing polygons into a raster.',\n",
       " 'vector_stream_network_analysis': 'This tool performs common stream network analysis operations on an input vector stream file.',\n",
       " 'viewshed': 'Identifies the viewshed for a point or set of points.',\n",
       " 'visibility_index': 'Estimates the relative visibility of sites in a DEM.',\n",
       " 'voronoi_diagram': 'Creates a vector Voronoi diagram for a set of vector points.',\n",
       " 'watershed': 'Identifies the watershed, or drainage basin, draining to a set of target cells.',\n",
       " 'weighted_overlay': 'Performs a weighted sum on multiple input rasters after converting each image to a common scale. The tool performs a multi-criteria evaluation (MCE).',\n",
       " 'weighted_sum': 'Performs a weighted-sum overlay on multiple input raster images.',\n",
       " 'wetness_index': 'Calculates the topographic wetness index, Ln(A / tan(slope)).',\n",
       " 'wilcoxon_signed_rank_test': 'Performs a 2-sample K-S test for significant differences on two input rasters.',\n",
       " 'write_function_memory_insertion': 'Performs a write function memory insertion for single-band multi-date change detection.',\n",
       " 'xor': 'Performs a logical XOR operator on two Boolean raster images.',\n",
       " 'z_scores': 'Standardizes the values in an input raster by converting to z-scores.',\n",
       " 'zlidar_to_las': 'Converts one or more zlidar files into the LAS data format.',\n",
       " 'zonal_statistics': 'Extracts descriptive statistics for a group of patches in a raster.'}"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "wbt.list_tools()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "bizarre-venezuela",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'gaussian_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'hypsometric_analysis': 'Calculates a hypsometric curve for one or more DEMs.',\n",
       " 'maximal_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'mean_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'minimal_curvature': 'Calculates a mean curvature raster from an input DEM.',\n",
       " 'pennock_landform_class': 'Classifies hillslope zones based on slope, profile curvature, and plan curvature.',\n",
       " 'plan_curvature': 'Calculates a plan (contour) curvature raster from an input DEM.',\n",
       " 'profile_curvature': 'Calculates a profile curvature raster from an input DEM.',\n",
       " 'tangential_curvature': 'Calculates a tangential curvature raster from an input DEM.',\n",
       " 'total_curvature': 'Calculates a total curvature raster from an input DEM.'}"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "wbt.list_tools(['curv'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "outdoor-bowling",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Hillshade\n",
      "Description:\n",
      "Calculates a hillshade raster from an input DEM.\n",
      "Toolbox: Geomorphometric Analysis\n",
      "Parameters:\n",
      "\n",
      "Flag               Description\n",
      "-----------------  -----------\n",
      "-i, --dem          Input raster DEM file.\n",
      "-o, --output       Output raster file.\n",
      "--azimuth          Illumination source azimuth in degrees.\n",
      "--altitude         Illumination source altitude in degrees.\n",
      "--zfactor          Optional multiplier for when the vertical and horizontal units are not the same.\n",
      "\n",
      "\n",
      "Example usage:\n",
      ">>./whitebox_tools -r=Hillshade -v --wd=\"/path/to/data/\" -i=DEM.tif -o=output.tif --azimuth=315.0 --altitude=30.0\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(wbt.tool_help(\"hillshade\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "elegant-inspiration",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "urban-certification",
   "metadata": {},
   "source": [
    "## Trabalhando com CopDEM"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "spiritual-banks",
   "metadata": {},
   "source": [
    "O arquivo TIF está disponível em: https://github.com/CarlosGrohmann/whitebox_scripts"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "convenient-coordinator",
   "metadata": {},
   "outputs": [],
   "source": [
    "cop_tif = 'CopDEM_clip_petar.tif'\n",
    "cop = cop_tif[:-4] # nome do arquivo sem extensão '.tif'   "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "id": "spectacular-danger",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"Hillshade\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_shade_315_35_z1' --azimuth=315 --altitude=35 --zfactor='1' -v --compress_rasters=True\n",
      "\n",
      "****************************\n",
      "* Welcome to Hillshade     *\n",
      "* Powered by WhiteboxTools *\n",
      "* www.whiteboxgeo.com      *\n",
      "****************************\n",
      "Reading data...\n",
      "Performing analysis: 0%\n",
      "Performing analysis: 1%\n",
      "Performing analysis: 2%\n",
      "Performing analysis: 3%\n",
      "Performing analysis: 4%\n",
      "Performing analysis: 5%\n",
      "Performing analysis: 6%\n",
      "Performing analysis: 7%\n",
      "Performing analysis: 8%\n",
      "Performing analysis: 9%\n",
      "Performing analysis: 10%\n",
      "Performing analysis: 11%\n",
      "Performing analysis: 12%\n",
      "Performing analysis: 13%\n",
      "Performing analysis: 14%\n",
      "Performing analysis: 15%\n",
      "Performing analysis: 16%\n",
      "Performing analysis: 17%\n",
      "Performing analysis: 18%\n",
      "Performing analysis: 19%\n",
      "Performing analysis: 20%\n",
      "Performing analysis: 21%\n",
      "Performing analysis: 22%\n",
      "Performing analysis: 23%\n",
      "Performing analysis: 24%\n",
      "Performing analysis: 25%\n",
      "Performing analysis: 26%\n",
      "Performing analysis: 27%\n",
      "Performing analysis: 28%\n",
      "Performing analysis: 29%\n",
      "Performing analysis: 30%\n",
      "Performing analysis: 31%\n",
      "Performing analysis: 32%\n",
      "Performing analysis: 33%\n",
      "Performing analysis: 34%\n",
      "Performing analysis: 35%\n",
      "Performing analysis: 36%\n",
      "Performing analysis: 37%\n",
      "Performing analysis: 38%\n",
      "Performing analysis: 39%\n",
      "Performing analysis: 40%\n",
      "Performing analysis: 41%\n",
      "Performing analysis: 42%\n",
      "Performing analysis: 43%\n",
      "Performing analysis: 44%\n",
      "Performing analysis: 45%\n",
      "Performing analysis: 46%\n",
      "Performing analysis: 47%\n",
      "Performing analysis: 48%\n",
      "Performing analysis: 49%\n",
      "Performing analysis: 50%\n",
      "Performing analysis: 51%\n",
      "Performing analysis: 52%\n",
      "Performing analysis: 53%\n",
      "Performing analysis: 54%\n",
      "Performing analysis: 55%\n",
      "Performing analysis: 56%\n",
      "Performing analysis: 57%\n",
      "Performing analysis: 58%\n",
      "Performing analysis: 59%\n",
      "Performing analysis: 60%\n",
      "Performing analysis: 61%\n",
      "Performing analysis: 62%\n",
      "Performing analysis: 63%\n",
      "Performing analysis: 64%\n",
      "Performing analysis: 65%\n",
      "Performing analysis: 66%\n",
      "Performing analysis: 67%\n",
      "Performing analysis: 68%\n",
      "Performing analysis: 69%\n",
      "Performing analysis: 70%\n",
      "Performing analysis: 71%\n",
      "Performing analysis: 72%\n",
      "Performing analysis: 73%\n",
      "Performing analysis: 74%\n",
      "Performing analysis: 75%\n",
      "Performing analysis: 76%\n",
      "Performing analysis: 77%\n",
      "Performing analysis: 78%\n",
      "Performing analysis: 79%\n",
      "Performing analysis: 80%\n",
      "Performing analysis: 81%\n",
      "Performing analysis: 82%\n",
      "Performing analysis: 83%\n",
      "Performing analysis: 84%\n",
      "Performing analysis: 85%\n",
      "Performing analysis: 86%\n",
      "Performing analysis: 87%\n",
      "Performing analysis: 88%\n",
      "Performing analysis: 89%\n",
      "Performing analysis: 90%\n",
      "Performing analysis: 91%\n",
      "Performing analysis: 92%\n",
      "Performing analysis: 93%\n",
      "Performing analysis: 94%\n",
      "Performing analysis: 95%\n",
      "Performing analysis: 96%\n",
      "Performing analysis: 97%\n",
      "Performing analysis: 98%\n",
      "Performing analysis: 99%\n",
      "Performing analysis: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.67s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# relevo sombreado\n",
    "az = 315\n",
    "alt = 35\n",
    "zf = 1\n",
    "shade = f'{cop}_shade_{az}_{alt}_z{zf}'\n",
    "\n",
    "wbt.hillshade(\n",
    "dem = cop_tif,\n",
    "output = shade,\n",
    "azimuth=az,\n",
    "altitude=alt,\n",
    "zfactor=zf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "id": "backed-holiday",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"Slope\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_slope' --units=percent -v --compress_rasters=True\n",
      "\n",
      "****************************\n",
      "* Welcome to Slope         *\n",
      "* Powered by WhiteboxTools *\n",
      "* www.whiteboxgeo.com      *\n",
      "****************************\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.167s\n",
      "./whitebox_tools --run=\"Aspect\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_aspect' -v --compress_rasters=True\n",
      "\n",
      "****************************\n",
      "* Welcome to Aspect        *\n",
      "* Powered by WhiteboxTools *\n",
      "* www.whiteboxgeo.com      *\n",
      "****************************\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.157s\n",
      "./whitebox_tools --run=\"PlanCurvature\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_plan_curvature' -v --compress_rasters=True\n",
      "\n",
      "****************************\n",
      "* Welcome to PlanCurvature *\n",
      "* Powered by WhiteboxTools *\n",
      "* www.whiteboxgeo.com      *\n",
      "****************************\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.166s\n",
      "./whitebox_tools --run=\"ProfileCurvature\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_profile_curvature' -v --compress_rasters=True\n",
      "\n",
      "*******************************\n",
      "* Welcome to ProfileCurvature *\n",
      "* Powered by WhiteboxTools    *\n",
      "* www.whiteboxgeo.com         *\n",
      "*******************************\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.158s\n",
      "./whitebox_tools --run=\"StandardDeviationOfSlope\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --input='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_rough_11x11' --filterx=11 --filtery=11 -v --compress_rasters=True\n",
      "\n",
      "***************************************\n",
      "* Welcome to StandardDeviationOfSlope *\n",
      "* Powered by WhiteboxTools            *\n",
      "* www.whiteboxgeo.com                 *\n",
      "***************************************\n",
      "Reading data...\n",
      "Performing analysis: 0%\n",
      "Performing analysis: 1%\n",
      "Performing analysis: 2%\n",
      "Performing analysis: 3%\n",
      "Performing analysis: 4%\n",
      "Performing analysis: 5%\n",
      "Performing analysis: 6%\n",
      "Performing analysis: 7%\n",
      "Performing analysis: 8%\n",
      "Performing analysis: 9%\n",
      "Performing analysis: 10%\n",
      "Performing analysis: 11%\n",
      "Performing analysis: 12%\n",
      "Performing analysis: 13%\n",
      "Performing analysis: 14%\n",
      "Performing analysis: 15%\n",
      "Performing analysis: 16%\n",
      "Performing analysis: 17%\n",
      "Performing analysis: 18%\n",
      "Performing analysis: 19%\n",
      "Performing analysis: 20%\n",
      "Performing analysis: 21%\n",
      "Performing analysis: 22%\n",
      "Performing analysis: 23%\n",
      "Performing analysis: 24%\n",
      "Performing analysis: 25%\n",
      "Performing analysis: 26%\n",
      "Performing analysis: 27%\n",
      "Performing analysis: 28%\n",
      "Performing analysis: 29%\n",
      "Performing analysis: 30%\n",
      "Performing analysis: 31%\n",
      "Performing analysis: 32%\n",
      "Performing analysis: 33%\n",
      "Performing analysis: 34%\n",
      "Performing analysis: 35%\n",
      "Performing analysis: 36%\n",
      "Performing analysis: 37%\n",
      "Performing analysis: 38%\n",
      "Performing analysis: 39%\n",
      "Performing analysis: 40%\n",
      "Performing analysis: 41%\n",
      "Performing analysis: 42%\n",
      "Performing analysis: 43%\n",
      "Performing analysis: 44%\n",
      "Performing analysis: 45%\n",
      "Performing analysis: 46%\n",
      "Performing analysis: 47%\n",
      "Performing analysis: 48%\n",
      "Performing analysis: 49%\n",
      "Performing analysis: 50%\n",
      "Performing analysis: 51%\n",
      "Performing analysis: 52%\n",
      "Performing analysis: 53%\n",
      "Performing analysis: 54%\n",
      "Performing analysis: 55%\n",
      "Performing analysis: 56%\n",
      "Performing analysis: 57%\n",
      "Performing analysis: 58%\n",
      "Performing analysis: 59%\n",
      "Performing analysis: 60%\n",
      "Performing analysis: 61%\n",
      "Performing analysis: 62%\n",
      "Performing analysis: 63%\n",
      "Performing analysis: 64%\n",
      "Performing analysis: 65%\n",
      "Performing analysis: 66%\n",
      "Performing analysis: 67%\n",
      "Performing analysis: 68%\n",
      "Performing analysis: 69%\n",
      "Performing analysis: 70%\n",
      "Performing analysis: 71%\n",
      "Performing analysis: 72%\n",
      "Performing analysis: 73%\n",
      "Performing analysis: 74%\n",
      "Performing analysis: 75%\n",
      "Performing analysis: 76%\n",
      "Performing analysis: 77%\n",
      "Performing analysis: 78%\n",
      "Performing analysis: 79%\n",
      "Performing analysis: 80%\n",
      "Performing analysis: 81%\n",
      "Performing analysis: 82%\n",
      "Performing analysis: 83%\n",
      "Performing analysis: 84%\n",
      "Performing analysis: 85%\n",
      "Performing analysis: 86%\n",
      "Performing analysis: 87%\n",
      "Performing analysis: 88%\n",
      "Performing analysis: 89%\n",
      "Performing analysis: 90%\n",
      "Performing analysis: 91%\n",
      "Performing analysis: 92%\n",
      "Performing analysis: 93%\n",
      "Performing analysis: 94%\n",
      "Performing analysis: 95%\n",
      "Performing analysis: 96%\n",
      "Performing analysis: 97%\n",
      "Performing analysis: 98%\n",
      "Performing analysis: 99%\n",
      "Performing analysis: 100%\n",
      "Creating integral images: 0%\n",
      "Creating integral images: 1%\n",
      "Creating integral images: 2%\n",
      "Creating integral images: 3%\n",
      "Creating integral images: 4%\n",
      "Creating integral images: 5%\n",
      "Creating integral images: 6%\n",
      "Creating integral images: 7%\n",
      "Creating integral images: 8%\n",
      "Creating integral images: 9%\n",
      "Creating integral images: 10%\n",
      "Creating integral images: 11%\n",
      "Creating integral images: 12%\n",
      "Creating integral images: 13%\n",
      "Creating integral images: 14%\n",
      "Creating integral images: 15%\n",
      "Creating integral images: 16%\n",
      "Creating integral images: 17%\n",
      "Creating integral images: 18%\n",
      "Creating integral images: 19%\n",
      "Creating integral images: 20%\n",
      "Creating integral images: 21%\n",
      "Creating integral images: 22%\n",
      "Creating integral images: 23%\n",
      "Creating integral images: 24%\n",
      "Creating integral images: 25%\n",
      "Creating integral images: 26%\n",
      "Creating integral images: 27%\n",
      "Creating integral images: 28%\n",
      "Creating integral images: 29%\n",
      "Creating integral images: 30%\n",
      "Creating integral images: 31%\n",
      "Creating integral images: 32%\n",
      "Creating integral images: 33%\n",
      "Creating integral images: 34%\n",
      "Creating integral images: 35%\n",
      "Creating integral images: 36%\n",
      "Creating integral images: 37%\n",
      "Creating integral images: 38%\n",
      "Creating integral images: 39%\n",
      "Creating integral images: 40%\n",
      "Creating integral images: 41%\n",
      "Creating integral images: 42%\n",
      "Creating integral images: 43%\n",
      "Creating integral images: 44%\n",
      "Creating integral images: 45%\n",
      "Creating integral images: 46%\n",
      "Creating integral images: 47%\n",
      "Creating integral images: 48%\n",
      "Creating integral images: 49%\n",
      "Creating integral images: 50%\n",
      "Creating integral images: 51%\n",
      "Creating integral images: 52%\n",
      "Creating integral images: 53%\n",
      "Creating integral images: 54%\n",
      "Creating integral images: 55%\n",
      "Creating integral images: 56%\n",
      "Creating integral images: 57%\n",
      "Creating integral images: 58%\n",
      "Creating integral images: 59%\n",
      "Creating integral images: 60%\n",
      "Creating integral images: 61%\n",
      "Creating integral images: 62%\n",
      "Creating integral images: 63%\n",
      "Creating integral images: 64%\n",
      "Creating integral images: 65%\n",
      "Creating integral images: 66%\n",
      "Creating integral images: 67%\n",
      "Creating integral images: 68%\n",
      "Creating integral images: 69%\n",
      "Creating integral images: 70%\n",
      "Creating integral images: 71%\n",
      "Creating integral images: 72%\n",
      "Creating integral images: 73%\n",
      "Creating integral images: 74%\n",
      "Creating integral images: 75%\n",
      "Creating integral images: 76%\n",
      "Creating integral images: 77%\n",
      "Creating integral images: 78%\n",
      "Creating integral images: 79%\n",
      "Creating integral images: 80%\n",
      "Creating integral images: 81%\n",
      "Creating integral images: 82%\n",
      "Creating integral images: 83%\n",
      "Creating integral images: 84%\n",
      "Creating integral images: 85%\n",
      "Creating integral images: 86%\n",
      "Creating integral images: 87%\n",
      "Creating integral images: 88%\n",
      "Creating integral images: 89%\n",
      "Creating integral images: 90%\n",
      "Creating integral images: 91%\n",
      "Creating integral images: 92%\n",
      "Creating integral images: 93%\n",
      "Creating integral images: 94%\n",
      "Creating integral images: 95%\n",
      "Creating integral images: 96%\n",
      "Creating integral images: 97%\n",
      "Creating integral images: 98%\n",
      "Creating integral images: 99%\n",
      "Creating integral images: 100%\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.38s\n",
      "./whitebox_tools --run=\"Geomorphons\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_geomorphons_d25' --search=25 --threshold=0.0 --fdist=0 --skip=0 --forms -v --compress_rasters=True\n",
      "\n",
      "****************************\n",
      "* Welcome to Geomorphons   *\n",
      "* Powered by WhiteboxTools *\n",
      "* www.whiteboxgeo.com      *\n",
      "****************************\n",
      "Reading data...\n",
      "Generating global ternary codes...\n",
      "Computing geomorphons...\n",
      "Progress: 0%\n",
      "Progress: 1%\n",
      "Progress: 2%\n",
      "Progress: 3%\n",
      "Progress: 4%\n",
      "Progress: 5%\n",
      "Progress: 6%\n",
      "Progress: 7%\n",
      "Progress: 8%\n",
      "Progress: 9%\n",
      "Progress: 10%\n",
      "Progress: 11%\n",
      "Progress: 12%\n",
      "Progress: 13%\n",
      "Progress: 14%\n",
      "Progress: 15%\n",
      "Progress: 16%\n",
      "Progress: 17%\n",
      "Progress: 18%\n",
      "Progress: 19%\n",
      "Progress: 20%\n",
      "Progress: 21%\n",
      "Progress: 22%\n",
      "Progress: 23%\n",
      "Progress: 24%\n",
      "Progress: 25%\n",
      "Progress: 26%\n",
      "Progress: 27%\n",
      "Progress: 28%\n",
      "Progress: 29%\n",
      "Progress: 30%\n",
      "Progress: 31%\n",
      "Progress: 32%\n",
      "Progress: 33%\n",
      "Progress: 34%\n",
      "Progress: 35%\n",
      "Progress: 36%\n",
      "Progress: 37%\n",
      "Progress: 38%\n",
      "Progress: 39%\n",
      "Progress: 40%\n",
      "Progress: 41%\n",
      "Progress: 42%\n",
      "Progress: 43%\n",
      "Progress: 44%\n",
      "Progress: 45%\n",
      "Progress: 46%\n",
      "Progress: 47%\n",
      "Progress: 48%\n",
      "Progress: 49%\n",
      "Progress: 50%\n",
      "Progress: 51%\n",
      "Progress: 52%\n",
      "Progress: 53%\n",
      "Progress: 54%\n",
      "Progress: 55%\n",
      "Progress: 56%\n",
      "Progress: 57%\n",
      "Progress: 58%\n",
      "Progress: 59%\n",
      "Progress: 60%\n",
      "Progress: 61%\n",
      "Progress: 62%\n",
      "Progress: 63%\n",
      "Progress: 64%\n",
      "Progress: 65%\n",
      "Progress: 66%\n",
      "Progress: 67%\n",
      "Progress: 68%\n",
      "Progress: 69%\n",
      "Progress: 70%\n",
      "Progress: 71%\n",
      "Progress: 72%\n",
      "Progress: 73%\n",
      "Progress: 74%\n",
      "Progress: 75%\n",
      "Progress: 76%\n",
      "Progress: 77%\n",
      "Progress: 78%\n",
      "Progress: 79%\n",
      "Progress: 80%\n",
      "Progress: 81%\n",
      "Progress: 82%\n",
      "Progress: 83%\n",
      "Progress: 84%\n",
      "Progress: 85%\n",
      "Progress: 86%\n",
      "Progress: 87%\n",
      "Progress: 88%\n",
      "Progress: 89%\n",
      "Progress: 90%\n",
      "Progress: 91%\n",
      "Progress: 92%\n",
      "Progress: 93%\n",
      "Progress: 94%\n",
      "Progress: 95%\n",
      "Progress: 96%\n",
      "Progress: 97%\n",
      "Progress: 98%\n",
      "Progress: 99%\n",
      "Progress: 100%\n",
      "Saving data...\n",
      "Output file written\n",
      "Elapsed Time (excluding I/O): 0.627s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# parâmetros (variáveis?) geomorfométricos\n",
    "\n",
    "# declividade\n",
    "slope = f'{cop}_slope'\n",
    "wbt.slope(dem=cop_tif, output=slope, units=\"percent\")\n",
    "\n",
    "# aspecto\n",
    "aspect = f'{cop}_aspect'\n",
    "wbt.aspect(dem=cop_tif, output=aspect)\n",
    "\n",
    "# curvaturas\n",
    "planc = f'{cop}_plan_curvature'\n",
    "profc = f'{cop}_profile_curvature'\n",
    "meanc = f'{cop}_mean_curvature'\n",
    "wbt.plan_curvature(dem=cop_tif, output=planc)\n",
    "wbt.profile_curvature(dem=cop_tif, output=profc)\n",
    "wbt.mean_curvature(dem=cop_tif, output=meanc)\n",
    "\n",
    "# escala log nas curvaturas\n",
    "planc_log = f'{cop}_plan_curvature_log'\n",
    "profc_log = f'{cop}_profile_curvature_log'\n",
    "meanc_log = f'{cop}_mean_curvature_log'\n",
    "wbt.plan_curvature(dem=cop_tif, output=planc_log,log=True)\n",
    "wbt.profile_curvature(dem=cop_tif, output=profc_log,log=True)\n",
    "wbt.mean_curvature(dem=cop_tif, output=meanc_log,log=True)\n",
    "\n",
    "# rugosidade\n",
    "win = 11\n",
    "rough = f'{cop}_rough_{win}x{win}'\n",
    "wbt.standard_deviation_of_slope(i=cop_tif, output=rough, filterx=win, filtery=win)\n",
    "\n",
    "# geomorphons\n",
    "dist = 35\n",
    "gmorph = f'{cop}_geomorphons_d{dist}'\n",
    "wbt.geomorphons(dem=cop_tif, output=gmorph, search=dist, forms=True)\n",
    "\n",
    "# landforms\n",
    "penn = f'{cop}_pennock'\n",
    "wbt.pennock_landform_class(dem=cop_tif, output=penn, slope=3.0, prof=0.1, plan=0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "novel-timeline",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "id": "compatible-ireland",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"SlopeVsElevationPlot\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --inputs='CopDEM_clip_petar.tif' --output='cop_elev_slope.html' -v --compress_rasters=True\n",
      "\n",
      "***********************************\n",
      "* Welcome to SlopeVsElevationPlot *\n",
      "* Powered by WhiteboxTools        *\n",
      "* www.whiteboxgeo.com             *\n",
      "***********************************\n",
      "Reading data...\n",
      "Slope analysis: 0%\n",
      "Slope analysis: 1%\n",
      "Slope analysis: 2%\n",
      "Slope analysis: 3%\n",
      "Slope analysis: 4%\n",
      "Slope analysis: 5%\n",
      "Slope analysis: 6%\n",
      "Slope analysis: 7%\n",
      "Slope analysis: 8%\n",
      "Slope analysis: 9%\n",
      "Slope analysis: 10%\n",
      "Slope analysis: 11%\n",
      "Slope analysis: 12%\n",
      "Slope analysis: 13%\n",
      "Slope analysis: 14%\n",
      "Slope analysis: 15%\n",
      "Slope analysis: 16%\n",
      "Slope analysis: 17%\n",
      "Slope analysis: 18%\n",
      "Slope analysis: 19%\n",
      "Slope analysis: 20%\n",
      "Slope analysis: 21%\n",
      "Slope analysis: 22%\n",
      "Slope analysis: 23%\n",
      "Slope analysis: 24%\n",
      "Slope analysis: 25%\n",
      "Slope analysis: 26%\n",
      "Slope analysis: 27%\n",
      "Slope analysis: 28%\n",
      "Slope analysis: 29%\n",
      "Slope analysis: 30%\n",
      "Slope analysis: 31%\n",
      "Slope analysis: 32%\n",
      "Slope analysis: 33%\n",
      "Slope analysis: 34%\n",
      "Slope analysis: 35%\n",
      "Slope analysis: 36%\n",
      "Slope analysis: 37%\n",
      "Slope analysis: 38%\n",
      "Slope analysis: 39%\n",
      "Slope analysis: 40%\n",
      "Slope analysis: 41%\n",
      "Slope analysis: 42%\n",
      "Slope analysis: 43%\n",
      "Slope analysis: 44%\n",
      "Slope analysis: 45%\n",
      "Slope analysis: 46%\n",
      "Slope analysis: 47%\n",
      "Slope analysis: 48%\n",
      "Slope analysis: 49%\n",
      "Slope analysis: 50%\n",
      "Slope analysis: 51%\n",
      "Slope analysis: 52%\n",
      "Slope analysis: 53%\n",
      "Slope analysis: 54%\n",
      "Slope analysis: 55%\n",
      "Slope analysis: 56%\n",
      "Slope analysis: 57%\n",
      "Slope analysis: 58%\n",
      "Slope analysis: 59%\n",
      "Slope analysis: 60%\n",
      "Slope analysis: 61%\n",
      "Slope analysis: 62%\n",
      "Slope analysis: 63%\n",
      "Slope analysis: 64%\n",
      "Slope analysis: 65%\n",
      "Slope analysis: 66%\n",
      "Slope analysis: 67%\n",
      "Slope analysis: 68%\n",
      "Slope analysis: 69%\n",
      "Slope analysis: 70%\n",
      "Slope analysis: 71%\n",
      "Slope analysis: 72%\n",
      "Slope analysis: 73%\n",
      "Slope analysis: 74%\n",
      "Slope analysis: 75%\n",
      "Slope analysis: 76%\n",
      "Slope analysis: 77%\n",
      "Slope analysis: 78%\n",
      "Slope analysis: 79%\n",
      "Slope analysis: 80%\n",
      "Slope analysis: 81%\n",
      "Slope analysis: 82%\n",
      "Slope analysis: 83%\n",
      "Slope analysis: 84%\n",
      "Slope analysis: 85%\n",
      "Slope analysis: 86%\n",
      "Slope analysis: 87%\n",
      "Slope analysis: 88%\n",
      "Slope analysis: 89%\n",
      "Slope analysis: 90%\n",
      "Slope analysis: 91%\n",
      "Slope analysis: 92%\n",
      "Slope analysis: 93%\n",
      "Slope analysis: 94%\n",
      "Slope analysis: 95%\n",
      "Slope analysis: 96%\n",
      "Slope analysis: 97%\n",
      "Slope analysis: 98%\n",
      "Slope analysis: 99%\n",
      "Slope analysis: 100%\n",
      "Binning the data: 0%\n",
      "Binning the data: 1%\n",
      "Binning the data: 2%\n",
      "Binning the data: 3%\n",
      "Binning the data: 4%\n",
      "Binning the data: 5%\n",
      "Binning the data: 6%\n",
      "Binning the data: 7%\n",
      "Binning the data: 8%\n",
      "Binning the data: 9%\n",
      "Binning the data: 10%\n",
      "Binning the data: 11%\n",
      "Binning the data: 12%\n",
      "Binning the data: 13%\n",
      "Binning the data: 14%\n",
      "Binning the data: 15%\n",
      "Binning the data: 16%\n",
      "Binning the data: 17%\n",
      "Binning the data: 18%\n",
      "Binning the data: 19%\n",
      "Binning the data: 20%\n",
      "Binning the data: 21%\n",
      "Binning the data: 22%\n",
      "Binning the data: 23%\n",
      "Binning the data: 24%\n",
      "Binning the data: 25%\n",
      "Binning the data: 26%\n",
      "Binning the data: 27%\n",
      "Binning the data: 28%\n",
      "Binning the data: 29%\n",
      "Binning the data: 30%\n",
      "Binning the data: 31%\n",
      "Binning the data: 32%\n",
      "Binning the data: 33%\n",
      "Binning the data: 34%\n",
      "Binning the data: 35%\n",
      "Binning the data: 36%\n",
      "Binning the data: 37%\n",
      "Binning the data: 38%\n",
      "Binning the data: 39%\n",
      "Binning the data: 40%\n",
      "Binning the data: 41%\n",
      "Binning the data: 42%\n",
      "Binning the data: 43%\n",
      "Binning the data: 44%\n",
      "Binning the data: 45%\n",
      "Binning the data: 46%\n",
      "Binning the data: 47%\n",
      "Binning the data: 48%\n",
      "Binning the data: 49%\n",
      "Binning the data: 50%\n",
      "Binning the data: 51%\n",
      "Binning the data: 52%\n",
      "Binning the data: 53%\n",
      "Binning the data: 54%\n",
      "Binning the data: 55%\n",
      "Binning the data: 56%\n",
      "Binning the data: 57%\n",
      "Binning the data: 58%\n",
      "Binning the data: 59%\n",
      "Binning the data: 60%\n",
      "Binning the data: 61%\n",
      "Binning the data: 62%\n",
      "Binning the data: 63%\n",
      "Binning the data: 64%\n",
      "Binning the data: 65%\n",
      "Binning the data: 66%\n",
      "Binning the data: 67%\n",
      "Binning the data: 68%\n",
      "Binning the data: 69%\n",
      "Binning the data: 70%\n",
      "Binning the data: 71%\n",
      "Binning the data: 72%\n",
      "Binning the data: 73%\n",
      "Binning the data: 74%\n",
      "Binning the data: 75%\n",
      "Binning the data: 76%\n",
      "Binning the data: 77%\n",
      "Binning the data: 78%\n",
      "Binning the data: 79%\n",
      "Binning the data: 80%\n",
      "Binning the data: 81%\n",
      "Binning the data: 82%\n",
      "Binning the data: 83%\n",
      "Binning the data: 84%\n",
      "Binning the data: 85%\n",
      "Binning the data: 86%\n",
      "Binning the data: 87%\n",
      "Binning the data: 88%\n",
      "Binning the data: 89%\n",
      "Binning the data: 90%\n",
      "Binning the data: 91%\n",
      "Binning the data: 92%\n",
      "Binning the data: 93%\n",
      "Binning the data: 94%\n",
      "Binning the data: 95%\n",
      "Binning the data: 96%\n",
      "Binning the data: 97%\n",
      "Binning the data: 98%\n",
      "Binning the data: 99%\n",
      "Binning the data: 100%\n",
      "\n",
      "Elapsed Time (excluding I/O): 0.139s\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m/tmp/ipykernel_6568/1758680703.py\u001b[0m in \u001b[0;36m?\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mwbt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mslope_vs_elevation_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minputs\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcop_tif\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moutput\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'cop_elev_slope.html'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mwatershed\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m/mnt/c/Users/PC/Documents/py3/lib/python3.8/site-packages/whitebox/whitebox_tools.py\u001b[0m in \u001b[0;36m?\u001b[0;34m(self, inputs, output, watershed, callback)\u001b[0m\n\u001b[1;32m   4510\u001b[0m         \u001b[0margs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   4511\u001b[0m         \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--inputs='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minputs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   4512\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mwatershed\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--watershed='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mwatershed\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   4513\u001b[0m         \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--output='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 4514\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun_tool\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'slope_vs_elevation_plot'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;31m# returns 1 if error\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m/mnt/c/Users/PC/Documents/py3/lib/python3.8/site-packages/whitebox/whitebox_tools.py\u001b[0m in \u001b[0;36m?\u001b[0;34m(self, tool_name, args, callback)\u001b[0m\n\u001b[1;32m    560\u001b[0m         \u001b[0;32mexcept\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mOSError\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mValueError\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mCalledProcessError\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0merr\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    561\u001b[0m             \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0merr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    562\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    563\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 564\u001b[0;31m             \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchdir\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mwork_dir\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "wbt.slope_vs_elevation_plot(inputs=cop_tif, output='cop_elev_slope.html', watershed=None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "id": "latter-occasion",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"HypsometricAnalysis\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --inputs='CopDEM_clip_petar.tif' --output='cop_hypso.html' -v --compress_rasters=True\n",
      "\n",
      "**********************************\n",
      "* Welcome to HypsometricAnalysis *\n",
      "* Powered by WhiteboxTools       *\n",
      "* www.whiteboxgeo.com            *\n",
      "**********************************\n",
      "Binning the data: 0%\n",
      "Binning the data: 1%\n",
      "Binning the data: 2%\n",
      "Binning the data: 3%\n",
      "Binning the data: 4%\n",
      "Binning the data: 5%\n",
      "Binning the data: 6%\n",
      "Binning the data: 7%\n",
      "Binning the data: 8%\n",
      "Binning the data: 9%\n",
      "Binning the data: 10%\n",
      "Binning the data: 11%\n",
      "Binning the data: 12%\n",
      "Binning the data: 13%\n",
      "Binning the data: 14%\n",
      "Binning the data: 15%\n",
      "Binning the data: 16%\n",
      "Binning the data: 17%\n",
      "Binning the data: 18%\n",
      "Binning the data: 19%\n",
      "Binning the data: 20%\n",
      "Binning the data: 21%\n",
      "Binning the data: 22%\n",
      "Binning the data: 23%\n",
      "Binning the data: 24%\n",
      "Binning the data: 25%\n",
      "Binning the data: 26%\n",
      "Binning the data: 27%\n",
      "Binning the data: 28%\n",
      "Binning the data: 29%\n",
      "Binning the data: 30%\n",
      "Binning the data: 31%\n",
      "Binning the data: 32%\n",
      "Binning the data: 33%\n",
      "Binning the data: 34%\n",
      "Binning the data: 35%\n",
      "Binning the data: 36%\n",
      "Binning the data: 37%\n",
      "Binning the data: 38%\n",
      "Binning the data: 39%\n",
      "Binning the data: 40%\n",
      "Binning the data: 41%\n",
      "Binning the data: 42%\n",
      "Binning the data: 43%\n",
      "Binning the data: 44%\n",
      "Binning the data: 45%\n",
      "Binning the data: 46%\n",
      "Binning the data: 47%\n",
      "Binning the data: 48%\n",
      "Binning the data: 49%\n",
      "Binning the data: 50%\n",
      "Binning the data: 51%\n",
      "Binning the data: 52%\n",
      "Binning the data: 53%\n",
      "Binning the data: 54%\n",
      "Binning the data: 55%\n",
      "Binning the data: 56%\n",
      "Binning the data: 57%\n",
      "Binning the data: 58%\n",
      "Binning the data: 59%\n",
      "Binning the data: 60%\n",
      "Binning the data: 61%\n",
      "Binning the data: 62%\n",
      "Binning the data: 63%\n",
      "Binning the data: 64%\n",
      "Binning the data: 65%\n",
      "Binning the data: 66%\n",
      "Binning the data: 67%\n",
      "Binning the data: 68%\n",
      "Binning the data: 69%\n",
      "Binning the data: 70%\n",
      "Binning the data: 71%\n",
      "Binning the data: 72%\n",
      "Binning the data: 73%\n",
      "Binning the data: 74%\n",
      "Binning the data: 75%\n",
      "Binning the data: 76%\n",
      "Binning the data: 77%\n",
      "Binning the data: 78%\n",
      "Binning the data: 79%\n",
      "Binning the data: 80%\n",
      "Binning the data: 81%\n",
      "Binning the data: 82%\n",
      "Binning the data: 83%\n",
      "Binning the data: 84%\n",
      "Binning the data: 85%\n",
      "Binning the data: 86%\n",
      "Binning the data: 87%\n",
      "Binning the data: 88%\n",
      "Binning the data: 89%\n",
      "Binning the data: 90%\n",
      "Binning the data: 91%\n",
      "Binning the data: 92%\n",
      "Binning the data: 93%\n",
      "Binning the data: 94%\n",
      "Binning the data: 95%\n",
      "Binning the data: 96%\n",
      "Binning the data: 97%\n",
      "Binning the data: 98%\n",
      "Binning the data: 99%\n",
      "Binning the data: 100%\n",
      "\n",
      "Elapsed Time (excluding I/O): 0.119s\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m/tmp/ipykernel_6568/4122187135.py\u001b[0m in \u001b[0;36m?\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mwbt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhypsometric_analysis\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minputs\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcop_tif\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moutput\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'cop_hypso.html'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mwatershed\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m/mnt/c/Users/PC/Documents/py3/lib/python3.8/site-packages/whitebox/whitebox_tools.py\u001b[0m in \u001b[0;36m?\u001b[0;34m(self, inputs, output, watershed, callback)\u001b[0m\n\u001b[1;32m   3624\u001b[0m         \u001b[0margs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   3625\u001b[0m         \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--inputs='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0minputs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   3626\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mwatershed\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--watershed='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mwatershed\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   3627\u001b[0m         \u001b[0margs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"--output='{}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 3628\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun_tool\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'hypsometric_analysis'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;31m# returns 1 if error\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m/mnt/c/Users/PC/Documents/py3/lib/python3.8/site-packages/whitebox/whitebox_tools.py\u001b[0m in \u001b[0;36m?\u001b[0;34m(self, tool_name, args, callback)\u001b[0m\n\u001b[1;32m    560\u001b[0m         \u001b[0;32mexcept\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mOSError\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mValueError\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mCalledProcessError\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0merr\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    561\u001b[0m             \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0merr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    562\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    563\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 564\u001b[0;31m             \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchdir\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mwork_dir\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "wbt.hypsometric_analysis(inputs=cop_tif, output='cop_hypso.html', watershed=None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "romantic-village",
   "metadata": {},
   "outputs": [],
   "source": [
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "institutional-sixth",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"LocalQuadraticRegression\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --dem='CopDEM_clip_petar.tif' --output='CopDEM_clip_petar_LQS_f31.tif' --filter=31 -v --compress_rasters=True\n",
      "\n",
      "***************************************\n",
      "* Welcome to LocalQuadraticRegression *\n",
      "* Powered by WhiteboxTools            *\n",
      "* www.whiteboxgeo.com                 *\n",
      "***************************************\n",
      "Reading data...\n",
      "Performing analysis: 0%\n",
      "Performing analysis: 1%\n",
      "Performing analysis: 2%\n",
      "Performing analysis: 3%\n",
      "Performing analysis: 4%\n",
      "Performing analysis: 5%\n",
      "Performing analysis: 6%\n",
      "Performing analysis: 7%\n",
      "Performing analysis: 8%\n",
      "Performing analysis: 9%\n",
      "Performing analysis: 10%\n",
      "Performing analysis: 11%\n",
      "Performing analysis: 12%\n",
      "Performing analysis: 13%\n",
      "Performing analysis: 14%\n",
      "Performing analysis: 15%\n",
      "Performing analysis: 16%\n",
      "Performing analysis: 17%\n",
      "Performing analysis: 18%\n",
      "Performing analysis: 19%\n",
      "Performing analysis: 20%\n",
      "Performing analysis: 21%\n",
      "Performing analysis: 22%\n",
      "Performing analysis: 23%\n",
      "Performing analysis: 24%\n",
      "Performing analysis: 25%\n",
      "Performing analysis: 26%\n",
      "Performing analysis: 27%\n",
      "Performing analysis: 28%\n",
      "Performing analysis: 29%\n",
      "Performing analysis: 30%\n",
      "Performing analysis: 31%\n",
      "Performing analysis: 32%\n",
      "Performing analysis: 33%\n",
      "Performing analysis: 34%\n",
      "Performing analysis: 35%\n",
      "Performing analysis: 36%\n",
      "Performing analysis: 37%\n",
      "Performing analysis: 38%\n",
      "Performing analysis: 39%\n",
      "Performing analysis: 40%\n",
      "Performing analysis: 41%\n",
      "Performing analysis: 42%\n",
      "Performing analysis: 43%\n",
      "Performing analysis: 44%\n",
      "Performing analysis: 45%\n",
      "Performing analysis: 46%\n",
      "Performing analysis: 47%\n",
      "Performing analysis: 48%\n",
      "Performing analysis: 49%\n",
      "Performing analysis: 50%\n",
      "Performing analysis: 51%\n",
      "Performing analysis: 52%\n",
      "Performing analysis: 53%\n",
      "Performing analysis: 54%\n",
      "Performing analysis: 55%\n",
      "Performing analysis: 56%\n",
      "Performing analysis: 57%\n",
      "Performing analysis: 58%\n",
      "Performing analysis: 59%\n",
      "Performing analysis: 60%\n",
      "Performing analysis: 61%\n",
      "Performing analysis: 62%\n",
      "Performing analysis: 63%\n",
      "Performing analysis: 64%\n",
      "Performing analysis: 65%\n",
      "Performing analysis: 66%\n",
      "Performing analysis: 67%\n",
      "Performing analysis: 68%\n",
      "Performing analysis: 69%\n",
      "Performing analysis: 70%\n",
      "Performing analysis: 71%\n",
      "Performing analysis: 72%\n",
      "Performing analysis: 73%\n",
      "Performing analysis: 74%\n",
      "Performing analysis: 75%\n",
      "Performing analysis: 76%\n",
      "Performing analysis: 77%\n",
      "Performing analysis: 78%\n",
      "Performing analysis: 79%\n",
      "Performing analysis: 80%\n",
      "Performing analysis: 81%\n",
      "Performing analysis: 82%\n",
      "Performing analysis: 83%\n",
      "Performing analysis: 84%\n",
      "Performing analysis: 85%\n",
      "Performing analysis: 86%\n",
      "Performing analysis: 87%\n",
      "Performing analysis: 88%\n",
      "Performing analysis: 89%\n",
      "Performing analysis: 90%\n",
      "Performing analysis: 91%\n",
      "Performing analysis: 92%\n",
      "Performing analysis: 93%\n",
      "Performing analysis: 94%\n",
      "Performing analysis: 95%\n",
      "Performing analysis: 96%\n",
      "Performing analysis: 97%\n",
      "Performing analysis: 98%\n",
      "Performing analysis: 99%\n",
      "Performing analysis: 100%\n",
      "Saving data...\n",
      "Output file 1 written\n",
      "Output file 2 written\n",
      "Output file 3 written\n",
      "Output file 4 written\n",
      "Output file 5 written\n",
      "Output file 6 written\n",
      "Output file 7 written\n",
      "Output file 8 written\n",
      "Elapsed Time (excluding I/O): 5.441s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# LQR - Wood 1996\n",
    "win = 31\n",
    "lwr = f'{cop}_LQS_f{win}.tif'\n",
    "wbt.local_quadratic_regression(dem=cop_tif, output=lqr, filter=win)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "greatest-compromise",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "imposed-award",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "./whitebox_tools --run=\"RasterHistogram\" --wd=\"/mnt/d/Dropbox/USP/extensao/videos_whitebox\" --input='CopDEM_clip_petar_aspect.tif' --output='CopDEM_clip_petar_aspect_hist.html' -v --compress_rasters=True\n",
      "\n",
      "******************************\n",
      "* Welcome to RasterHistogram *\n",
      "* Powered by WhiteboxTools   *\n",
      "* www.whiteboxgeo.com        *\n",
      "******************************\n",
      "Binning the data: 0%\n",
      "Binning the data: 1%\n",
      "Binning the data: 2%\n",
      "Binning the data: 3%\n",
      "Binning the data: 4%\n",
      "Binning the data: 5%\n",
      "Binning the data: 6%\n",
      "Binning the data: 7%\n",
      "Binning the data: 8%\n",
      "Binning the data: 9%\n",
      "Binning the data: 10%\n",
      "Binning the data: 11%\n",
      "Binning the data: 12%\n",
      "Binning the data: 13%\n",
      "Binning the data: 14%\n",
      "Binning the data: 15%\n",
      "Binning the data: 16%\n",
      "Binning the data: 17%\n",
      "Binning the data: 18%\n",
      "Binning the data: 19%\n",
      "Binning the data: 20%\n",
      "Binning the data: 21%\n",
      "Binning the data: 22%\n",
      "Binning the data: 23%\n",
      "Binning the data: 24%\n",
      "Binning the data: 25%\n",
      "Binning the data: 26%\n",
      "Binning the data: 27%\n",
      "Binning the data: 28%\n",
      "Binning the data: 29%\n",
      "Binning the data: 30%\n",
      "Binning the data: 31%\n",
      "Binning the data: 32%\n",
      "Binning the data: 33%\n",
      "Binning the data: 34%\n",
      "Binning the data: 35%\n",
      "Binning the data: 36%\n",
      "Binning the data: 37%\n",
      "Binning the data: 38%\n",
      "Binning the data: 39%\n",
      "Binning the data: 40%\n",
      "Binning the data: 41%\n",
      "Binning the data: 42%\n",
      "Binning the data: 43%\n",
      "Binning the data: 44%\n",
      "Binning the data: 45%\n",
      "Binning the data: 46%\n",
      "Binning the data: 47%\n",
      "Binning the data: 48%\n",
      "Binning the data: 49%\n",
      "Binning the data: 50%\n",
      "Binning the data: 51%\n",
      "Binning the data: 52%\n",
      "Binning the data: 53%\n",
      "Binning the data: 54%\n",
      "Binning the data: 55%\n",
      "Binning the data: 56%\n",
      "Binning the data: 57%\n",
      "Binning the data: 58%\n",
      "Binning the data: 59%\n",
      "Binning the data: 60%\n",
      "Binning the data: 61%\n",
      "Binning the data: 62%\n",
      "Binning the data: 63%\n",
      "Binning the data: 64%\n",
      "Binning the data: 65%\n",
      "Binning the data: 66%\n",
      "Binning the data: 67%\n",
      "Binning the data: 68%\n",
      "Binning the data: 69%\n",
      "Binning the data: 70%\n",
      "Binning the data: 71%\n",
      "Binning the data: 72%\n",
      "Binning the data: 73%\n",
      "Binning the data: 74%\n",
      "Binning the data: 75%\n",
      "Binning the data: 76%\n",
      "Binning the data: 77%\n",
      "Binning the data: 78%\n",
      "Binning the data: 79%\n",
      "Binning the data: 80%\n",
      "Binning the data: 81%\n",
      "Binning the data: 82%\n",
      "Binning the data: 83%\n",
      "Binning the data: 84%\n",
      "Binning the data: 85%\n",
      "Binning the data: 86%\n",
      "Binning the data: 87%\n",
      "Binning the data: 88%\n",
      "Binning the data: 89%\n",
      "Binning the data: 90%\n",
      "Binning the data: 91%\n",
      "Binning the data: 92%\n",
      "Binning the data: 93%\n",
      "Binning the data: 94%\n",
      "Binning the data: 95%\n",
      "Binning the data: 96%\n",
      "Binning the data: 97%\n",
      "Binning the data: 98%\n",
      "Binning the data: 99%\n",
      "Binning the data: 100%\n",
      "\n",
      "Elapsed Time (excluding I/O): 0.248s\n",
      "Complete! Please see /mnt/d/Dropbox/USP/extensao/videos_whitebox/CopDEM_clip_petar_aspect_hist.html for output.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# histograma de raster\n",
    "r_in = aspect\n",
    "wbt.raster_histogram(i=f'{r_in}.tif', output=f'{r_in}_hist.html')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "impaired-reunion",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'CopDEM_clip_petar_aspect'"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "neither-validation",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "logical-monitoring",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "pleased-bouquet",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "least-september",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "lasting-border",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}