Use a random file name for vrt to avoid pb when run in parallel

We used gdal.BuildVRT("/vsimem/var.vrt", ...) but this can lead to
problems if several files are generated at the same time when run in
parallel. We now use a random file name with 6 digits.

forestatrisk/predict/predict_raster.py

@@ -14,6 +14,7 @@
 from glob import glob
 import os
 import sys
+import random
 
 # Third party imports
 import numpy as np
@@ -90,10 +91,12 @@ def predict_raster(
         yRes=-gt[5],
         separate=True,
     )
+    rand_int = random.randint(1, 999999)
+    vrt_file = f"/vsimem/var_{rand_int:06}.vrt"
     cback = gdal.TermProgress if verbose else 0
-    gdal.BuildVRT("/vsimem/var.vrt", raster_list,
+    gdal.BuildVRT(vrt_file, raster_list,
                   options=param, callback=cback)
-    stack = gdal.Open("/vsimem/var.vrt")
+    stack = gdal.Open(vrt_file)
     nband = stack.RasterCount
     proj = stack.GetProjection()
 
@@ -109,7 +112,7 @@ def predict_raster(
     bandND = bandND.astype(np.float32)
 
     # Make blocks
-    blockinfo = makeblock("/vsimem/var.vrt", blk_rows=blk_rows)
+    blockinfo = makeblock(vrt_file, blk_rows=blk_rows)
     nblock = blockinfo[0]
     nblock_x = blockinfo[1]
     x = blockinfo[3]

forestatrisk/predict/predict_raster_binomial_iCAR.py

@@ -14,6 +14,7 @@
 from glob import glob
 import os
 import sys
+import random
 
 # Third party imports
 import numpy as np
@@ -108,10 +109,12 @@ def predict_raster_binomial_iCAR(
         yRes=-gt[5],
         separate=True,
     )
+    rand_int = random.randint(1, 999999)
+    vrt_file = f"/vsimem/var_{rand_int:06}.vrt"
     cback = gdal.TermProgress if verbose else 0
-    gdal.BuildVRT("/vsimem/var.vrt", raster_list,
+    gdal.BuildVRT(vrt_file, raster_list,
                   options=param, callback=cback)
-    stack = gdal.Open("/vsimem/var.vrt")
+    stack = gdal.Open(vrt_file)
     nband = stack.RasterCount
     proj = stack.GetProjection()
 
@@ -127,7 +130,7 @@ def predict_raster_binomial_iCAR(
     bandND = bandND.astype(np.float32)
 
     # Make blocks
-    blockinfo = makeblock("/vsimem/var.vrt", blk_rows=blk_rows)
+    blockinfo = makeblock(vrt_file, blk_rows=blk_rows)
     nblock = blockinfo[0]
     nblock_x = blockinfo[1]
     x = blockinfo[3]