Commit 76e9e751 authored by Poggio, Laura's avatar Poggio, Laura

Update markdown/webdav_from_R.md, markdown/webdav_from_Python.md,...

Update markdown/webdav_from_R.md, markdown/webdav_from_Python.md, markdown/webdav_from_bash.md files
parent 2eb2f05a
<a href="https://www.isric.org" rel="isric.org"> <img src="https://www.isric.org/themes/custom/basic/logo.svg" height="130"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<a href="https://soilgrids.org" rel="soilgrids.org"> <img src="https://www.isric.org/sites/default/files/styles/gallery_big_image_900x700/public/SoilGrids_banner_web.png" height="130">
# [WebDAV](https://files.isric.org/soilgrids/data/recent/): direct access to the maps in VRT format.
### Python3
#### Import packages, select boundary box and cell size
Initially you need to import the following package and select the bounding box of the area you are interested in:
```python
from osgeo import gdal,ogr,osr
bb=-337500.000,1242500.000,152500.000,527500.000 # Example bounding box (homolosine)
igh='+proj=igh +lat_0=0 +lon_0=0 +datum=WGS84 +units=m +no_defs' # proj string for Homolosine projection
res=250
```
#### To a geotiff in Homolosine
This GDAL command will create a local geotiff in the Homolosine projection
``` python
kwargs = {'format': 'GTiff', 'projWin': bb, 'projWinSRS': igh, 'xRes': res, 'yRes': res, 'creationOptions': ["TILED=YES", "COMPRESS=DEFLATE", "PREDICTOR=2", "BIGTIFF=YES"]}
ds = gdal.Translate('./crop_roi_igh_py.tif',
'/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs/ocs_0-30cm_mean.vrt',
**kwargs)
del ds
```
#### To a geotiff in a different projection
The following commands describe a workflow to obtain a VRT or a GeoTiff for an area of interest in a projection of your choice. In this example we will use EPSG=4326
##### To local VRT in homolosine (directly from the webdav connection)
The first step is to obtain a VRT for the area of interest in the Homolosine projection. We suggest to use VRT for the intermediate steps to save space and computation times.
```python
kwargs = {'format': 'GTiff', 'projWin': bb, 'projWinSRS': igh, 'xRes': res, 'yRes': res}
ds = gdal.Translate('./crop_roi_igh_py.vrt',
'/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs/ocs_0-30cm_mean.vrt',
**kwargs)
del ds
```
##### To a VRT in, for example, LatLong
The following command will generate a VRT in the projection of your choice:
```python
ds = gdal.Warp('./crop_roi_ll_py.vrt',
'./crop_roi_igh_py.vrt',
dstSRS='EPSG:4326')
del ds
```
##### To a final Geotiff
The following command will generate a Geotiff in the projection of your choice for the area of interest defined above
```python
kwargs = {'format': 'GTiff', 'creationOptions': ["TILED=YES", "COMPRESS=DEFLATE", "PREDICTOR=2", "BIGTIFF=YES"] }
ds = gdal.Translate('./crop_roi_ll_py.tif',
'./crop_roi_ll_py.vrt',
**kwargs)
del ds
```
<a href="https://www.isric.org" rel="isric.org"> <img src="https://www.isric.org/themes/custom/basic/logo.svg" height="130"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<a href="https://soilgrids.org" rel="soilgrids.org"> <img src="https://www.isric.org/sites/default/files/styles/gallery_big_image_900x700/public/SoilGrids_banner_web.png" height="130">
# [WebDAV](https://files.isric.org/soilgrids/data/recent/): direct access to the maps in VRT format.
### R
#### Load libraries, select boundary box and cell size
Initially you need to load the following libraries and select the bounding box of the area you are interested in:
```R
library(rgdal)
library(gdalUtils)
bb=c(-337500.000,1242500.000,152500.000,527500.000) # Example bounding box (homolosine)
igh='+proj=igh +lat_0=0 +lon_0=0 +datum=WGS84 +units=m +no_defs' # proj string for Homolosine projection
sgurl=
```
#### To a geotiff in Homolosine
This GDAL command will create a local geotiff in the Homolosine projection
```
gdal_translate('/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs/ocs_0-30cm_mean.vrt',
"./crop_roi_igh_r.tif",
tr=c(250,250),
projwin=bb,
projwin_srs =igh,
verbose=TRUE)
```
#### To a geotiff in a different projection
The following commands describe a workflow to obtain a VRT or a GeoTiff for an area of interest in a projection of your choice. In this example we will use EPSG=4326
##### To local VRT in homolosine (directly from the webdav connection)
The first step is to obtain a VRT for the area of interest in the Homolosine projection. We suggest to use VRT for the intermediate steps to save space and computation times.
```
gdal_translate('/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs/ocs_0-30cm_mean.vrt',
"./crop_roi_igh_r.vrt",
of="VRT",tr=c(250,250),
projwin=bb,
projwin_srs =igh,
verbose=TRUE)
```
##### To a VRT in, for example, LatLong
The following command will generate a VRT in the projection of your choice:
```
gdalwarp("./crop_roi_igh_r.vrt",
"./crop_roi_ll_r.vrt",
s_src=igh,
t_srs="EPSG:4326",
of="VRT")
```
##### To a final Geotiff
The following command will generate a Geotiff in the projection of your choice for the area of interest defined above
```
gdal_translate("./crop_roi_ll_r.vrt",
"./crop_roi_ll_r.tif",
co=c("TILED=YES","COMPRESS=DEFLATE","PREDICTOR=2","BIGTIFF=YES"))
```
#### Read in R
Finally you can read any of the generated VRTs or GeoTiffs in R for further analysis
```
r=readGDAL("./crop_roi_igh_r.vrt") # Or any other of the files produce above
```
<a href="https://www.isric.org" rel="isric.org"> <img src="https://www.isric.org/themes/custom/basic/logo.svg" height="130"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<a href="https://soilgrids.org" rel="soilgrids.org"> <img src="https://www.isric.org/sites/default/files/styles/gallery_big_image_900x700/public/SoilGrids_banner_web.png" height="130">
# [WebDAV](https://files.isric.org/soilgrids/data/recent/): direct access to the maps in VRT format.
### Linux Bash
#### Select boundary box and cell size
Initially you need to import the following package and select the bounding box of the area you are interested in:
```Bash
BOUNDS=-337500.000 1242500.000 152500.000 527500.000 # Example bounding box (homolosine)
CELL_SIZE=250 250
```
#### To a geotiff in Homolosine
This GDAL command will create a local geotiff in the Homolosine projection
``` bash
gdal_translate -projwin $BOUNDS -tr $CELL_SIZE \
-co "TILED=YES" -co "COMPRESS=DEFLATE" -co "PREDICTOR=2" -co "BIGTIFF=YES" \
"/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs_0-30cm_mean.vrt" \
"ocs_0-5cm_mean.tif"
```
#### To a geotiff in a different projection
The following commands describe a workflow to obtain a VRT or a GeoTiff for an area of interest in a projection of your choice. In this example we will use EPSG=4326
##### To local VRT in homolosine (directly from the webdav connection)
The first step is to obtain a VRT for the area of interest in the Homolosine projection. We suggest to use VRT for the intermediate steps to save space and computation times.
``` bash
gdal_translate -of VRT -projwin $BOUNDS -tr $CELL_SIZE \
-co "TILED=YES" -co "COMPRESS=DEFLATE" -co "PREDICTOR=2" -co "BIGTIFF=YES" \
"/vsicurl?max_retry=3&retry_delay=1&list_dir=no&url=https://files.isric.org/soilgrids/latest/data/ocs_0-30cm_mean.vrt" \
"ocs_0-5cm_mean.vrt"
```
##### To a VRT in, for example, LatLong
The following command will generate a VRT in the projection of your choice:
```bash
gdalwarp -overwrite -t_srs EPSG:4326 -of VRT "ocs_0-5cm_mean.vrt" "ocs_0-5cm_mean_4326.vrt"
```
##### To a final Geotiff
The following command will generate a Geotiff in the projection of your choice for the area of interest defined above
```bash
gdal_translate ocs_0-5cm_mean_4326.vrt ocs_0-5cm_mean_4326.tif \
-co "TILED=YES" -co "COMPRESS=DEFLATE" -co "PREDICTOR=2" -co "BIGTIFF=YES"
```
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment