I love processing with vrt's. you can make lots of interim changes. Evaluate them quickly in QGIS and if you like any of them just translate back to a selfcontained raster format (tif, png etc).
saves lots of time.
U2ros,your uses of vrt's makes perfect sense, to me at least :)
mosaicking and then clipping is what I originally started using vrts for: eliminating interim rasters that I would delete later.
check out this link:
http://www.perrygeo.com/lazy-raster-processing-with-gdal-vrts.html
[link adjusted since the blog has been moved]
hope it has something for you.
I was able to read in the data using gdalwarp after downloading the entire file structure for a .vrt file to my computer. Trying to read the files directly from the WebDAV server must have been too large of a request to complete before timing out/dropping the connection.
I used wget for windows to download the files, and it took about 4 hours for each map (~15000 files, ~5 GB)
wget https://files.isric.org/soilgrids/latest/data/soc/soc_0-5cm_mean/ -r -np -nH --cut-dirs=3 -nc
wget https://files.isric.org/soilgrids/latest/data/soc/soc_0-5cm_mean.vrt -O ./soc/soc_0-5cm_mean.vrt
Then I was able to run the following command without a problem, converting the entire map to a tiff in EPSG 4326 at 7.5" resolution:
gdalUtilities::gdalwarp(t_srs="EPSG:4326", multi=TRUE, wm=200,
co=c("BIGTIFF=YES", "COMPRESS=DEFLATE", "TILED=TRUE"),
tr=c(1/480,1/480), # Desired output resolution
"soc_0-5cm_mean.vrt", # Input VRT
"soc_0-5cm_mean.tif") # Output file
Unsurprisingly, this command also took quite a long time to run (~45 minutes).
Best Answer
You need to change it to the following
Considering a file
debug.vrt
with the above content, it's working using below commandThe changes made are
adresses-33.csv
has to beadresses-33.csv.gz
because url https://adresse.data.gouv.fr/data/ban/adresses/latest/csv/adresses-33.csv returns a 404 whereas https://adresse.data.gouv.fr/data/ban/adresses/latest/csv/adresses-33.csv.gz exists<GeometryField encoding="PointFromColumns" x="lon" y="lat"/>
to make GDAL understand how to interpret columns as geometry