Global Soil Erodibility - Williams K-Factor - Dataset

Field	Value
Author	Natural Capital Project, Stanford University
Tags	ERODIBILITY ISRIC SOIL SOIL ERODIBILITY
Places	GLOBAL
Last Updated	March 20, 2025, 22:17 (UTC)
Created	March 20, 2025, 22:15 (UTC)
Natcap Last Updated	2025-03-20T22:17:25.400256+00:00

Cite this dataset

Embed this layer into your website

<div id="map" style="width: 100%; height: 400px;"></div>

<link rel="stylesheet" href="https://unpkg.com/leaflet@1.9.4/dist/leaflet.css" integrity="sha256-p4NxAoJBhIIN+hmNHrzRCf9tD/miZyoHS5obTRR9BMY=" crossorigin="" />
<script src="https://unpkg.com/leaflet@1.9.4/dist/leaflet.js" integrity="sha256-20nQCchB9co0qIjJZRGuk2/Z9VM+kNiyxNV1lvTlZBo=" crossorigin=""></script>
<script>
  var map = L.map('map').setView([0, 0], 4);

  L.tileLayer('https://tile.openstreetmap.org/{z}/{x}/{y}.png', {
    attribution: '&copy; <a href="https://www.openstreetmap.org/copyright">OpenStreetMap</a> contributors'
  }).addTo(map);

L.tileLayer('https://titiler-897938321824.us-west1.run.app/cog/tiles/WebMercatorQuad/{z}/{x}/{y}@2x?colormap_name=viridis&bidx=1&url=https%3A%2F%2Fstorage.googleapis.com%2Fnatcap-data-cache%2Fglobal%2Fwilliams-k-factor%2Fwilliams-k-factor.tif&rescale=0.0020859260112047195%2C0.006525821518152952').addTo(map);
</script>

Use this data in Python

from osgeo import gdal

raster = gdal.OpenEx('/vsicurl/https://storage.googleapis.com/natcap-data-cache/global/williams-k-factor/williams-k-factor.tif')
band = raster.GetRasterBand(1)
array = band.ReadAsArray(
    xoff=my_window_x_offset,
    yoff=my_window_y_offset,
    win_xsize=my_window_width_in_pixels,
    win_ysize=my_window_height_in_pixels)