Rasterio Source shape is inconsistent with given indexes 1

Question

I need to save out a 3 band geotiff to file. I am currently using rasterio and when I go to write out the 3 band imagery I get the error Source shape (1, 3445, 4703, 4) is inconsistent with given indexes 1.

My final goal is to be able to perform some analysis on an image and write this out to file.

I have already tried reshape_as_raster and reshape_as_image. I have tried a few other combinations as well as .transpose(arr, (0,1,2))

https://rasterio.readthedocs.io/en/stable/topics/image_processing.html#imageorder

with rio.open(r"C:\Users\name\Documents\project\name.tif") as src:
naip_data = src.read()
naip_meta = src.profile

image = reshape_as_raster(naip_data)

with rio.open('C:\\Users\\name\\Documents\\UAV_test_save\\filename.tif',     'w',**naip_meta) as dst:
        dst.write(image, 3)

I am expecting to have a geotiff saved in the file. instead I get:

ValueError rasterio_io.pyx in rasterio._io.DatasetWriterBase.write()

ValueError: Source shape (1, 3445, 4, 4703) is inconsistent with given indexes 1

I cant find anyone with this issue, Any help out there?

freegnome
– freegnome

2019-10-08 22:27:26 +00:00
Commented Oct 8, 2019 at 22:27 — freegnome
– freegnome, Commented Oct 8, 2019 at 22:27

sehan2 · Accepted Answer · 2020-03-03 17:05:37Z

10

I have had the same issue. I solved it by removing the 'indexes' argument.

so instead of:

dst.write(image,indexes=3)

use this:

dst.write(image)

answered Mar 3, 2020 at 17:05

sehan2

1,8551 gold badge14 silver badges27 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

mikey Over a year ago

Thanks. I was able to solve the issue with the opposite approach (I added , indexes = 1 to dst.write)

Marius Van Nieuwenhuyse Over a year ago

To clarify a bit, if you have a 2d array, you should include the indexes and if you have a 3d one, you should not

befpy · Accepted Answer · 2020-02-20 20:45:46Z

2

It looks like you have a 4 band raster that you are trying to write.

rio is reading (1, 3445, 4703, 4) as 1 band, 3445 rows, 4703 columns, and 4 of another dimension...I'm not a rio expert.

If you want to keep your 4 bands:

naip_data2 = np.moveaxis(naip_data.squeeze(),-1,0) # move axis with 4 entries to beginning and remove extra dimension

print(naip_data2.shape) # check new shape, should be (4,3445,4703)
print(naip_meta) # check that naip_data2.shape == (naip_meta['count'],naip_meta['height'],naip_meta['width'])

# if they are equal, naip_meta doesn't need to be updated, otherwise update it
with rio.open(fname,'w',**naip_meta) as dst:
    dst.write(naip_data2)

answered Feb 20, 2020 at 20:45

befpy

413 bronze badges

Comments

Ek00f · Accepted Answer · 2021-12-11 17:41:18Z

2

Source metadata contains field count and it is equal to 1. You need to update that value to the desired number of output bands.

source_meta.update({
    "count": 3
})

edited Dec 11, 2021 at 17:41

Ek00f

32 bronze badges

answered Aug 28, 2021 at 7:38

Казыбек Аскарбек

312 bronze badges

Comments

CalebJ · Accepted Answer · 2022-08-11 10:45:05Z

0

Good source

with rio.open(r"C:\Users\name\Documents\project\name.tif") as src:
naip_data = src.read()
naip_meta = src.profile

image = naip_data.reshape(3445,4703,4)

with rio.open('C:\\Users\\name\\Documents\\UAV_test_save\\filename.tif',     
'w',**naip_meta) as dst:
    dst.write(image, [0,1,2],[2,1,0])

answered Aug 11, 2022 at 10:45

CalebJ

111 bronze badge

Collectives™ on Stack Overflow

Rasterio Source shape is inconsistent with given indexes 1

4 Answers 4

2 Comments

Comments

Comments

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

4 Answers 4

2 Comments

Comments

Comments

Comments

Your Answer

Sign up or log in

Post as a guest

Related