RasterioIOError: not recognized as a supported file format.

[swhisnant5201]

Hi All,

I am currently working on a project with the Tobler package of PYSAL. When it comes to Dasymetric Interpolation, I need to use a raster to create a mask. That said, I am encountering an error when I attempt to read my raster.

My code is below:

import os
import geopandas as gpd
import matplotlib.pyplot as plt
import tobler
import rasterio
from zipfile import ZipFile

os.chdir("C:\\gispy\\Tobler_Tutorial")

def unzip_data(zip_file): 
    my_zip = ZipFile(zip_file, "r")
    my_zip.extractall()
    del(my_zip)
#reads and extracts files to working directory

def areal_inter(source_file, target_file, variable, title=""):
   result = tobler.area_weighted.area_interpolate(source_df=source_file, target_df=target_file, extensive_variables = [variable]) 
   fig, ax = plt.subplots(1,2, figsize=(14,7))

   result.plot(variable, scheme='quantiles',  ax=ax[0], legend = True, alpha = 0.5, k = 5, edgecolor = 'black' )
   source_file.plot(variable, scheme='quantiles',  ax=ax[1], legend = True, alpha = 0.5, k = 5, edgecolor = "black")

   ax[0].set_title('interpolated')
   ax[1].set_title('original')

   for ax in ax:
       ax.axis('off')
   fig.suptitle(title)
#calls tobler's area weighted areal interpolation method, plots output data
   
def dasy_inter(source_file, target_rast, target_file, variable, title=""):
   results = tobler.dasymetric.masked_area_interpolate(raster=target_rast, source_df=source_file, target_df=target_file, intensive_variables=variable, pixel_values = [17]) 
   fig, ax = plt.subplots(1,2, figsize=(14,7))

   results.plot(variable, scheme='quantiles',   ax=ax[0])
   source_file.plot(variable, scheme='quantiles',  ax=ax[1])

   ax[0].set_title('interpolated')
   ax[1].set_title('original')
   for ax in ax:
       ax.axis('off')
   fig.suptitle(title) 
#calls tobler's dasymetric areal interpolation, plots output data
    
    
nhood = gpd.read_file("Neighborhoods_Philadelphia")
pop_stat = gpd.read_file("Vital_Population_CT")
my_rast = rasterio.open("C:\\gispy\\Tobler_Tutorial\\land_cover_2020_30m_tif\\NA_NALCMS_landcover_2020_30m\\data\\NA_NALCMS_landcover_2020_30m.tif")

print(pop_stat.columns)

pop_stat = pop_stat[pop_stat["COUNT_ALL_"] != 0]
#remove census tracts where no one lives, since they can't have a percentage

crs = 2272
#set the crs for Pennsylvania South USft

pop_stat = pop_stat.to_crs(crs)
#reproject the necessary data to the correct CRS for Philadelphia
  
areal_inter(pop_stat, nhood, "PERCENT_AS", "Asian Percentage")
areal_inter(pop_stat, nhood, "PERCENT_BL", "African American Percentage")
dasy_inter(pop_stat, my_rast ,nhood,"PERCENT_AS", "Asian Percentage")

and the error traceback is:

Traceback (most recent call last):

File C:\miniconda3\envs\tutorials\Lib\encodings\cp1252.py:23 in decode
return codecs.charmap_decode(input,self.errors,decoding_table)[0]

UnicodeDecodeError: 'charmap' codec can't decode byte 0x8d in position 229: character maps to

Exception ignored in: 'rasterio._filepath.filepath_read'
Traceback (most recent call last):
File "C:\miniconda3\envs\tutorials\Lib\encodings\cp1252.py", line 23, in decode
return codecs.charmap_decode(input,self.errors,decoding_table)[0]
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
UnicodeDecodeError: 'charmap' codec can't decode byte 0x8d in position 229: character maps to
Traceback (most recent call last):

File rasterio_base.pyx:310 in rasterio._base.DatasetBase.init

File rasterio_base.pyx:221 in rasterio._base.open_dataset

File rasterio_err.pyx:221 in rasterio._err.exc_wrap_pointer

CPLE_OpenFailedError: '/vsipythonfilelike/dad878e5-6eea-4d45-9fad-1048e9f7f89f/dad878e5-6eea-4d45-9fad-1048e9f7f89f' not recognized as a supported file format.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):

File C:\miniconda3\envs\tutorials\Lib\site-packages\spyder_kernels\py3compat.py:356 in compat_exec
exec(code, globals, locals)

File c:\gispy\tobler_tutorial\tobler tutorial.py:78
dasy_inter(pop_stat, my_rast ,nhood,"PERCENT_AS", "Asian Percentage")

File c:\gispy\tobler_tutorial\tobler tutorial.py:44 in dasy_inter
results = tobler.dasymetric.masked_area_interpolate(raster=target_rast, source_df=source_file, target_df=target_file, intensive_variables=variable, pixel_values = [17])

File C:\miniconda3\envs\tutorials\Lib\site-packages\tobler\dasymetric\masked_area_interpolate.py:69 in masked_area_interpolate
raster_mask = extract_raster_features(

File C:\miniconda3\envs\tutorials\Lib\site-packages\tobler\dasymetric\raster_tools.py:106 in extract_raster_features
with rio.open(raster_path) as src:

File C:\miniconda3\envs\tutorials\Lib\site-packages\rasterio\env.py:451 in wrapper
return f(*args, **kwds)

File C:\miniconda3\envs\tutorials\Lib\site-packages\rasterio_init_.py:242 in open
return FilePath(fp).open(driver=driver, sharing=sharing, **kwargs)

File C:\miniconda3\envs\tutorials\Lib\site-packages\rasterio\env.py:398 in wrapper
return f(*args, **kwds)

File C:\miniconda3\envs\tutorials\Lib\site-packages\rasterio\io.py:236 in open
return DatasetReader(mempath, driver=driver, sharing=sharing, **kwargs)

File rasterio_base.pyx:312 in rasterio._base.DatasetBase.init

RasterioIOError: '/vsipythonfilelike/dad878e5-6eea-4d45-9fad-1048e9f7f89f/dad878e5-6eea-4d45-9fad-1048e9f7f89f' not recognized as a supported file format.

Any guidance on this is greatly appreciated.