# Calculate magnitude and direction of slopes
mag, direction = dem_proc.calc_slopes_directions()
# Calculate UCA/TWI on entire tile
area1 = dem_proc.calc_uca(0)
twi1 = dem_proc.calc_twi()
# Calculate UCA/TWI on chunks
dem_proc.resolve_edges = True # If false, won't do edge resolution
dem_proc.chunk_size_uca = 9 # Chunk size
dem_proc.chunk_overlap_uca = 3 #Overlap between chunks
area2 = dem_proc.calc_uca(0)
twi2 = dem_proc.calc_twi()
# Calculate UCA/TWI on chunk and apply UCA limit on edges
dem_proc.apply_uca_limit_edges = True
area3 = dem_proc.calc_uca(0)
twi3 = dem_proc.calc_twi()
# plot results
import matplotlib.pyplot as plt
plt.matshow(area1); plt.colorbar()
plt.title('UCA Calculated on Full Tile')
plt.matshow(area2 - area1); plt.colorbar()
plt.title('UCA Calculation: (on chunk) - (on full tile)')
plt.matshow(area3 - area2); plt.colorbar()
plt.title('UCA Calculation: (on chunk, limited) - (on chunk)')
plt.show()
plt.show()
# Calculate magnitude and direction of slopes
mag, direction = dem_proc.calc_slopes_directions()
# Calculate UCA/TWI on entire tile
area1 = dem_proc.calc_uca(0)
twi1 = dem_proc.calc_twi()
# Calculate UCA/TWI on chunks
dem_proc.resolve_edges = True # If false, won't do edge resolution
dem_proc.chunk_size_uca = 9 # Chunk size
dem_proc.chunk_overlap_uca = 3 #Overlap between chunks
area2 = dem_proc.calc_uca(0)
twi2 = dem_proc.calc_twi()
# Calculate UCA/TWI on chunk and apply UCA limit on edges
dem_proc.apply_uca_limit_edges = True
area3 = dem_proc.calc_uca(0)
twi3 = dem_proc.calc_twi()
# plot results
import matplotlib.pyplot as plt
plt.matshow(area1)
plt.colorbar()
plt.title('UCA Calculated on Full Tile')
plt.matshow(area2 - area1)
plt.colorbar()
plt.title('UCA Calculation: (on chunk) - (on full tile)')
plt.matshow(area3 - area2)
plt.colorbar()
plt.title('UCA Calculation: (on chunk, limited) - (on chunk)')
plt.show()
