So I read in a .tif image with rasterio and what I get out when printing dataset.read(1) is what I understood a coordinate system in the "red" band of the image.
The coordinate system is in UTM, the bottom is mostly 0's and to the top I get thousands like so:
array([[1389, 831, 1017, ..., 2221, 2869, 2729],
[1109, 693, 831, ..., 2683, 2499, 2823],
[ 877, 785, 693, ..., 2037, 2683, 2683],
...,
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0]]
I am trying to make sense of it what confuses me more is if I do read(2) which I'm guessing should correspond to the green channel I get:
array([[1249, 831, 925, ..., 2221, 2591, 2453],
[1063, 739, 785, ..., 2591, 2407, 2591],
[ 925, 785, 693, ..., 2129, 2545, 2499],
...,
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0]], dtype=uint16)
Slightly different values, but they should be the same if they correspond to the same location?
This is the images bounding box:
BoundingBox(left=246000.00000000035, bottom=9855098.500000004, right=251000.00000000035, top=9860000.000000004)
It's a tile of Nairobi, Kenya.