Bug in point_query: Getting the wrong value

Open yagciali2002 opened this issue 3 years ago • 0 comments

I am getting a wrong value after the point query on tif files below using the shp file below. For validation, I overlapped all of them in QGIS software and got the values manually. I also run "sample raster values tool" on QGIS and I got the same values as the manual one . All files are in the same coordinate system, EPSG 32635 I uploaded all of the files here https://drive.google.com/drive/folders/1rdkYfdjeLTqUZENwZHPNRyvtwOhDuY5X?usp=sharing

Point file: ECOR_Tower_Kirklareli_EPSG32635.shp File 1: LC09_L2SP_181031_20220604_20220606_02_T1_NDVI.tif File 2: LC09_L2SP_181031_20220604_20220606_02_T1_Ts.tif

The expected value of file 1 seen in QGIS software: 2225 (Row: 4127 and Col: 2926) The expected value of file 2 seen in QGIS software: 4212 (Row: 4127 and Col: 2926) I am getting: 2229.435607077935 on File 1 using rasterstats point query I am getting: 4205.417518870619 on File 2 using rasterstats point query

anaconda environment dependencies installed rasterstats.version : '0.17.0' rasterio.version : '1.2.10' shapely.version: '1.8.0' numpy.version : '1.21.5' affine.version: '2.3.0' fiona.version : '1.8.20' simplejson.version: '3.17.6' cligj.version: '0.7.2'

Sep 19 '22 15:09 yagciali2002