Python function to determine x and y extent of a raster

Question

I have been using gdal_grid to convert xyz files to a raster

I am wondering is there anyway to automatically detect the extents to use from the xyz file txe and tye using python so that I can use a more general code every time I grid?

This way I will be stating my input xyz file and interpolation method but the raster extents will be pre-calculated for best results using python??

e.g.

gdal_grid  -a nearest:radius1=1.0:radius2=1.0:angle=0.0:min_points=0:nodata=0.0 -ot Float32 -of GTiff -txe [pre-calculated] -tye [pre-calculated] -tr 1 1 "input" output

Answer 1

You can use Rasterio to read the xyz file and the bounds method to get the bounding box of the dataset. For example (sample data):

import rasterio

xyz = '/temp/small.xyz'

dataset = rasterio.open(xyz)
dataset.bounds

Which would yield:

Out[1]: BoundingBox(left=172762.5, bottom=210787.5, right=172937.5, top=210637.5)

Then you can pull individual coords from the bounding box object:

dataset.bounds[3]

Which would yield:

Out[1]: 210637.5