GDALRasterizeOptions Failed to process SRS definition

Optical Toolbox sen2like
#1

Hi everyone,
after the fix of the input archive structure the tool finally start working but it gets stucked on another error while running single-tile-mode:

[INFO ] 2021-04-07 13:01:06 - reader - Sentinel2MTL Class
[INFO ] 2021-04-07 13:01:06 - reader - Product: sen2like-master/data/Sentinel2/33TWF/S2A_MSIL1C_20200313T095031_N0209_R079_T33TWF_20200313T102505.SAFE
[INFO ] 2021-04-07 13:01:10 - sentinel2 - SAT_AZ , SAT_ZENITH, SUN_AZ, SUN_ZENITH
[INFO ] 2021-04-07 13:01:10 - sentinel2 - UNIT = DEGREES (scale: x100) :
[INFO ] 2021-04-07 13:01:10 - sentinel2 - /home/ec2-user/wd/14572/S2A_MSIL1C_20200313T095031_N0209_R079_T33TWF_20200313T102505.SAFE/tie_points.tif
[INFO ] 2021-04-07 13:01:10 - sentinel2 - Generating nodata mask from band B01
[INFO ] 2021-04-07 13:01:11 - image_file - Written: /home/ec2-user/wd/14572/S2A_MSIL1C_20200313T095031_N0209_R079_T33TWF_20200313T102505.SAFE/nodata_pixel_mask_B01.tif
[INFO ] 2021-04-07 13:01:11 - sentinel2 - Generating validity mask from cloud mask

Traceback (most recent call last):
File “sen2like-master/sen2like/sen2like/sen2like.py”, line 435, in
sys.exit(main(with_multiprocess_support=True))
File “sen2like-master/sen2like/sen2like/sen2like.py”, line 430, in main
start_process(tile, _products, args, start_date, end_date)
File “sen2like-master/sen2like/sen2like/sen2like.py”, line 253, in start_process
process(_product, args.bands)
File “sen2like-master/sen2like/sen2like/sen2like.py”, line 284, in process
product.mtl.get_valid_pixel_mask(os.path.join(config.get(“wd”), product.name, ‘valid_pixel_mask.tif’))
File “/home/ec2-user/SageMaker/sen2like-master/sen2like/sen2like/core/readers/sentinel2.py”, line 404, in get_valid_pixel_mask
xRes=res, yRes=res)
File “/home/ec2-user/SageMaker/custom-miniconda/miniconda/envs/sen2like_ker/lib/python3.6/site-packages/osgeo/gdal.py”, line 1163, in Rasterize
(opts, callback, callback_data) = RasterizeOptions(**kwargs)
File “/home/ec2-user/SageMaker/custom-miniconda/miniconda/envs/sen2like_ker/lib/python3.6/site-packages/osgeo/gdal.py”, line 1150, in RasterizeOptions
return (GDALRasterizeOptions(new_options), callback, callback_data)
File “/home/ec2-user/SageMaker/custom-miniconda/miniconda/envs/sen2like_ker/lib/python3.6/site-packages/osgeo/gdal.py”, line 3530, in init
this = _gdal.new_GDALRasterizeOptions(*args)
RuntimeError: Failed to process SRS definition: EPSG:32633

I don’t understand the error, the installed version of gdal is the 3.0.1 as required.
gdal_rasterize documentation says “The <srs_def> may be any of the usual GDAL/OGR forms, complete WKT, PROJ.4, EPSG:n” so the parameter should accept “EPSG:32633”

Any idea?
Thank you

#2

Some EPSG:n SRS definitions refer to a specific region, but remote sensing images can cover multiple regions, so in general you need a global SRS.

For EPSG:326330:

Area of use: Between 12°E and 18°E, northern hemisphere between equator and 84°N, onshore and offshore. Austria. Bosnia and Herzegovina. Cameroon. Central African Republic. Chad. Congo. Croatia. Czechia. Democratic Republic of the Congo (Zaire). Gabon. Germany. Hungary. Italy. Libya. Malta. Niger. Nigeria. Norway. Poland. San Marino. Slovakia. Slovenia. Svalbard. Sweden. Vatican City State.

#3

Dear Gnwii, thank you for your answer!

All Landsat8 and Sentinel2 products I have in my input archive are all from a small area in south Italy, so it should be fine right? still I get the error unfortunately.

The thing that confuses me most is that I am running:
python3 sen2like.py single-tile-mode 33TWF --wd ~/wd
So, apart from the organization of the files in my input archive I don’t really know how to handle the error I encountered, what parameters or things should I change?

#4

Solved making sure the virtual environment was properly created, correct miniconda and packages versions