-
Notifications
You must be signed in to change notification settings - Fork 551
Description
Expected behavior and actual behavior.
A seg fault occurs when I open tif files with fsspec and set boundless=True
Steps to reproduce the problem.
This causes a segmentation fault
import rasterio
import fsspec
with fsspec.open('https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/13/S/DV/2020/4/S2B_13SDV_20200428_0_L2A/TCI.tif') as fobj:
with rasterio.open(fobj) as src:
window = rasterio.windows.Window(100,100,100,100)
data = src.read(1, window=window, boundless=True)
Setting boundless=False works fine. Also opening the tif url directly with rasterio works fine with boundless=True
import rasterio
with rasterio.open('https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/13/S/DV/2020/4/S2B_13SDV_20200428_0_L2A/TCI.tif') as src:
window = rasterio.windows.Window(100,100,100,100)
data = src.read(1, window=window, boundless=True)
Environment Information
rasterio info:
rasterio: 1.3.7
GDAL: 3.6.4
PROJ: 9.0.1
GEOS: 3.11.1
PROJ DATA: /home/shawn/miniconda3/envs/geo2/lib/python3.9/site-packages/rasterio/proj_data
GDAL DATA: /home/shawn/miniconda3/envs/geo2/lib/python3.9/site-packages/rasterio/gdal_data
System:
python: 3.9.16 (main, Mar 8 2023, 14:00:05) [GCC 11.2.0]
executable: /home/shawn/miniconda3/envs/geo2/bin/python
machine: Linux-5.15.0-73-generic-x86_64-with-glibc2.35
Python deps:
affine: 2.4.0
attrs: 23.1.0
certifi: 2023.05.07
click: 8.1.3
cligj: 0.7.2
cython: None
numpy: 1.24.3
snuggs: 1.4.7
click-plugins: None
setuptools: 67.8.0
Installation Method
In a fresh conda environment on linux using pip 23.01