对上面的代码进行修改,获取卫星影像边界范围
时间: 2024-03-17 11:44:48 浏览: 49
要获取卫星影像的边界范围,可以使用GDAL库的`GetGeoTransform()`方法获取影像的地理变换参数,并根据这些参数计算影像的范围。以下是修改后的示例代码:
```python
from osgeo import gdal
# 读取卫星影像
src_ds = gdal.Open('image.tif')
# 获取影像的地理变换参数
transform = src_ds.GetGeoTransform()
x_origin = transform[0]
y_origin = transform[3]
pixel_width = transform[1]
pixel_height = transform[5]
cols = src_ds.RasterXSize
rows = src_ds.RasterYSize
# 计算影像的范围
x_min = x_origin
y_max = y_origin
x_max = x_origin + cols * pixel_width
y_min = y_origin + rows * pixel_height
# 输出影像的范围
print((x_min, y_min, x_max, y_max))
```
这段代码将会读取名为`image.tif`的卫星影像,获取影像的地理变换参数,根据这些参数计算影像的范围,并输出影像的范围。
阅读全文