code_01.R

# devtools::install_github("tylermorganwall/rayshader")
library(magrittr)
library(raster)

## Loading required package: sp

## 
## Attaching package: 'raster'

## The following object is masked from 'package:magrittr':
## 
##     extract

library(rayshader)

setwd("/home/michael/Dropbox/BGU/Allan_Just/p_05_shaded_relief_image_mexico_city/")

# Create DEM (from ASTER data)
r1 = raster("ASTGTM2_N19W099_dem.tif")
r2 = raster("ASTGTM2_N19W100_dem.tif")
r = merge(r1, r2)
e = extent(-99.61167649466057, -98.64365263741942, 18.973314712760843, 19.883630141311155)
r = crop(r, e)

# Aggregate
r = aggregate(r, 5)

# Plot DEM
plot(r)

# Calculate RGB for shaded relief
m = as.matrix(r)
s = 
  m %>%
  sphere_shade(remove_edges = FALSE) %>%
  add_water(detect_water(m, remove_edges = FALSE)) %>%
  add_shadow(ray_shade(m, remove_edges = FALSE, lambert = FALSE)) %>%
  add_shadow(lamb_shade(m, remove_edges = FALSE)) %>%
  add_shadow(ambient_shade(m, remove_edges = FALSE))

## 
Calculating Surface Normal [-----------------------------------] ETA:  0s
Calculating Surface Normal [===================================] ETA:  0s
                                                                         

Calculating Surface Normal [-----------------------------------] ETA:  0s
Calculating Surface Normal [===================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [===------------------------------------------------] ETA:  3s
Raytracing [====-----------------------------------------------] ETA:  3s
Raytracing [====-----------------------------------------------] ETA:  4s
Raytracing [=====----------------------------------------------] ETA:  4s
Raytracing [=====----------------------------------------------] ETA:  5s
Raytracing [======---------------------------------------------] ETA:  5s
Raytracing [=======--------------------------------------------] ETA:  5s
Raytracing [=======--------------------------------------------] ETA:  6s
Raytracing [========-------------------------------------------] ETA:  6s
Raytracing [=========------------------------------------------] ETA:  6s
Raytracing [==========-----------------------------------------] ETA:  6s
Raytracing [===========----------------------------------------] ETA:  6s
Raytracing [============---------------------------------------] ETA:  6s
Raytracing [=============--------------------------------------] ETA:  6s
Raytracing [==============-------------------------------------] ETA:  6s
Raytracing [===============------------------------------------] ETA:  6s
Raytracing [================-----------------------------------] ETA:  6s
Raytracing [=================----------------------------------] ETA:  6s
Raytracing [==================---------------------------------] ETA:  6s
Raytracing [===================--------------------------------] ETA:  6s
Raytracing [===================--------------------------------] ETA:  5s
Raytracing [====================-------------------------------] ETA:  5s
Raytracing [=====================------------------------------] ETA:  5s
Raytracing [======================-----------------------------] ETA:  5s
Raytracing [=======================----------------------------] ETA:  5s
Raytracing [========================---------------------------] ETA:  5s
Raytracing [=========================--------------------------] ETA:  5s
Raytracing [=========================--------------------------] ETA:  4s
Raytracing [==========================-------------------------] ETA:  4s
Raytracing [===========================------------------------] ETA:  4s
Raytracing [============================-----------------------] ETA:  4s
Raytracing [=============================----------------------] ETA:  4s
Raytracing [==============================---------------------] ETA:  4s
Raytracing [===============================--------------------] ETA:  4s
Raytracing [===============================--------------------] ETA:  3s
Raytracing [================================-------------------] ETA:  3s
Raytracing [=================================------------------] ETA:  3s
Raytracing [==================================-----------------] ETA:  3s
Raytracing [===================================----------------] ETA:  3s
Raytracing [====================================---------------] ETA:  3s
Raytracing [=====================================--------------] ETA:  3s
Raytracing [=====================================--------------] ETA:  2s
Raytracing [======================================-------------] ETA:  2s
Raytracing [=======================================------------] ETA:  2s
Raytracing [========================================-----------] ETA:  2s
Raytracing [=========================================----------] ETA:  2s
Raytracing [==========================================---------] ETA:  2s
Raytracing [===========================================--------] ETA:  2s
Raytracing [===========================================--------] ETA:  1s
Raytracing [============================================-------] ETA:  1s
Raytracing [=============================================------] ETA:  1s
Raytracing [==============================================-----] ETA:  1s
Raytracing [===============================================----] ETA:  1s
Raytracing [================================================---] ETA:  1s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  0s
Raytracing [=================----------------------------------] ETA:  0s
Raytracing [==================---------------------------------] ETA:  0s
Raytracing [===================--------------------------------] ETA:  0s
Raytracing [====================-------------------------------] ETA:  0s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [===============------------------------------------] ETA:  0s
Raytracing [================-----------------------------------] ETA:  0s
Raytracing [=================----------------------------------] ETA:  0s
Raytracing [==================---------------------------------] ETA:  0s
Raytracing [===================--------------------------------] ETA:  0s
Raytracing [====================-------------------------------] ETA:  0s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  0s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  0s
Raytracing [=================----------------------------------] ETA:  0s
Raytracing [==================---------------------------------] ETA:  0s
Raytracing [===================--------------------------------] ETA:  0s
Raytracing [====================-------------------------------] ETA:  0s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [===========----------------------------------------] ETA:  1s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [======================-----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  0s
Raytracing [====================-------------------------------] ETA:  0s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [======================-----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [======================-----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  0s
Raytracing [======================-----------------------------] ETA:  0s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [======================-----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  1s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s
                                                                         

Raytracing [---------------------------------------------------] ETA:  0s
Raytracing [============---------------------------------------] ETA:  1s
Raytracing [=============--------------------------------------] ETA:  1s
Raytracing [==============-------------------------------------] ETA:  1s
Raytracing [===============------------------------------------] ETA:  1s
Raytracing [================-----------------------------------] ETA:  1s
Raytracing [=================----------------------------------] ETA:  1s
Raytracing [==================---------------------------------] ETA:  1s
Raytracing [===================--------------------------------] ETA:  1s
Raytracing [====================-------------------------------] ETA:  1s
Raytracing [=====================------------------------------] ETA:  1s
Raytracing [======================-----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  1s
Raytracing [=======================----------------------------] ETA:  0s
Raytracing [========================---------------------------] ETA:  0s
Raytracing [=========================--------------------------] ETA:  0s
Raytracing [==========================-------------------------] ETA:  0s
Raytracing [===========================------------------------] ETA:  0s
Raytracing [============================-----------------------] ETA:  0s
Raytracing [=============================----------------------] ETA:  0s
Raytracing [==============================---------------------] ETA:  0s
Raytracing [===============================--------------------] ETA:  0s
Raytracing [================================-------------------] ETA:  0s
Raytracing [=================================------------------] ETA:  0s
Raytracing [==================================-----------------] ETA:  0s
Raytracing [===================================----------------] ETA:  0s
Raytracing [====================================---------------] ETA:  0s
Raytracing [=====================================--------------] ETA:  0s
Raytracing [======================================-------------] ETA:  0s
Raytracing [=======================================------------] ETA:  0s
Raytracing [========================================-----------] ETA:  0s
Raytracing [=========================================----------] ETA:  0s
Raytracing [==========================================---------] ETA:  0s
Raytracing [===========================================--------] ETA:  0s
Raytracing [============================================-------] ETA:  0s
Raytracing [=============================================------] ETA:  0s
Raytracing [==============================================-----] ETA:  0s
Raytracing [===============================================----] ETA:  0s
Raytracing [================================================---] ETA:  0s
Raytracing [=================================================--] ETA:  0s
Raytracing [==================================================-] ETA:  0s
Raytracing [===================================================] ETA:  0s

s = aperm(s, c(2, 1, 3)) # Transpose 'x' and 'y'

# To RGB raster
u = brick(s)
extent(u) = extent(r)
proj4string(u) = proj4string(r)

# Plot Shaded relief
# png("mexico_city_shaded_relief.png", width = 7, height = 7, units = "in", res = 400)
plotRGB(u, scale = 1, maxpixels = ncell(u))

# dev.off()

code_01.R

michael

Tue Jul 10 13:56:13 2018