The openrspat from ohm-np

A highly motivated and passionate geoinformatician with over three years of working experience in the geospatial field. My research interests include geographic information science (GIS), remote sensing, cloud computing and data analytics.

openrspat's People

Contributors

Watchers

openrspat's Issues

Update the package version

Update the package version with some changes:

well executed download functions
well executed processing routines
technical details of how the package works
necessary vignettes

Progress track

Download functions

create global functions

Create global function zonal_operation for following rasters processing:

Carbon balance
World pop - population count
Terrain Ruggedness Index
Clay content
Accessibility to cities

### rasters only
## routines [crop, mask, rasterize, zonal]
## required arguments [polygon object - sf, raster object - terra, rasterize column, zonal operation]

# create function
library("sf")
library("terra")
library("openRspat")
zonal_operation <- function(pol, rast, opn) {
  
  # transform crs of pol to area_proj()
  pol <- st_transform(pol,
                      "+proj=longlat +datum=WGS84 +no_defs")
  # convert sf to spatvector
  pol_v <- terra::vect(pol)
  print("to vect")
  # crop the raster
  crop <- terra::crop(rast,
                      pol_v)
  print("cropped")
  # mask the raster
  mask <- terra::mask(crop,
                      pol_v)
  print("masked")
  # get unique ID
  r_id <- pol_v$WDPA_PID %>%
    as.numeric()
  print("got id")
  # rasterize the polygon
  r <- terra::rasterize(pol_v,
                        mask,
                        r_id)
  print("rasterized")
  # compute zonal statistics
  z <- terra::zonal(mask, 
                    r,
                    fun=opn,
                    na.rm=T)
  print("zonal statistics")
  # return zonal statistics as data frame
  return(z)
}

pol <- read_sf("../wdpa_kfw_spatial_latinamerica_2021-02-01_supportedPAs_unique.gpkg")
rast <- rast("../../datalake/mapme.protectedareas/input/clay_content/clay_content_30_cm.tif")
opn <- 'mean'

x <- 
  zonal_operation(pol, rast, opn)

Recommend Projects