Hello,  I've been having a problem using windows, probably a misunderstanding on my part.
My goal is to dump a specific data from a band into a csv file, first I tried to use a window but the results differ from the results that I see on Qgis(the results appear to be shifted 1 position).
So I decided to try to read the whole band, instead to use a window, and the result match with Qgis.
I think that I am doing something wrong when I create/use the window, what I tried to do was:

1. create a window from the bounds of the desired area: 
   bandWin = rasterio.windows.from_bounds(*boundingBox,band.transform)
2. get the transform matrix of that window:
   win_transform = band.window_transform(bandWin)
3. Obtain the window data:
   win_data = band.read(1, window = bandWin)
4. Get the desired window indexes:
   indexes = rasterio.transform.rowcol(win_transform, [Xs], [Ys])
5. Get the values that I want:
   csvData.append(win_data[indexes[0],indexes[1]][0])

The only differences between when I don't use a window is that instead use win_transform I use band.transform and I read the whole band.

All help is appreciated.

I am using rasterio 1.1.5.

--

Hello Sean,

Sorry for abandoning the post.
I installed rasterio via pip, and recently i had to reinstall my OS, and also of course rasterio, and the problem still persists.
I will prepare a demo and a link to a sample data that am I using, soon as I can.

Hello, this is the example:

import rasterio

#data link: https://earthexplorer.usgs.gov/download/external/options/SENTINEL_2A/3022286/INVSVC/

 

pathToImgFolder = ""

pathData = pathToImgFolder + "L1C_T21LXC_A001666_20170701T140052\S2B_MSIL1C_20170701T140049_N0205_R067_T21LXC_20170701T140052.SAFE\GRANULE\L1C_T21LXC_A001666_20170701T140052\IMG_DATA\T21LXC_20170701T140049_B01.jp2" 

boudingBoxCoordinates = (606859.0750363453, 8241169.269917269, 607219.0750363453, 8241529.269917269)

points = [(607014.0750363453,8241374.26991727),(607024.0750363444,8241374.26991727),\

          (607034.0750363453,8241374.26991727),(607044.0750363453,8241374.26991727),\

          (607054.0750363442,8241374.26991727)]

 

with rasterio.open(pathData) as img:

    bandWindow = rasterio.windows.from_bounds(*boudingBoxCoordinates, img.transform)

    winTransform = rasterio.windows.transform(bandWindow,img.transform)

    bandData = img.read(1, window = bandWindow)

    allData = img.read(1)

    for point in points:

        rBand,cBand = rasterio.transform.rowcol(winTransform,point[0],point[1])

        rFull,cFull = img.index(point[0],point[1])

        bandVal = bandData[rBand,cBand]

        fullVal = allData[rFull,cFull]

        print(f"{fullVal} {'=' if bandVal == fullVal else '!='} {bandVal}")

When I execute this example i get this result:

1202 = 1202

1330 != 1202

1330 != 1202

1330 = 1330

1330 = 1330

I asked the same question on stack exchange and now I understand how to obtain the results that once I expected.
But I think the answer for that is a little bit in the dark, should we improve the documentation to have an example that explain this?