python gdal 栅格循环

栅格循环是指在处理栅格数据时，对多个波段数据进行循环操作。在Python GDAL中，可以使用循环语句遍历每个波段，并对其进行相应的处理操作，如写入数据等。下面是一个示例代码，展示了如何使用循环对多个波段数据进行赋值操作：

import numpy as np
from osgeo import gdal

driver = gdal.GetDriverByName("GTiff")
outdata = driver.Create("c:/tmp/new_image.tif", xsize=300, ysize=200, bands=3)

for i in range(3):
    outband = outdata.GetRasterBand(i + 1)
    in_array = np.random.rand(200, 300) * 255
    outband.WriteArray(in_array.astype(int))

outdata.FlushCache()
outdata = None

这段代码首先创建了一个3个波段的栅格数据文件，然后通过循环对每个波段进行赋值操作。循环中，我们使用`np.random.rand(200, 300) * 255`生成一个随机数组，并将其转换为整数类型后写入到对应的波段中。最后，刷新缓存并关闭数据文件。

相关问题

gdal获取栅格最大最小值

使用GDAL库可以方便地获取栅格数据的最大和最小值。下面是一个示例代码：

import gdal

# 打开栅格文件
dataset = gdal.Open('path/to/raster.tif')

# 获取栅格数据的波段数量
band_count = dataset.RasterCount

# 循环遍历每个波段并获取最大最小值
for i in range(1, band_count + 1):
    band = dataset.GetRasterBand(i)
    band_min, band_max, _, _ = band.ComputeStatistics(False)
    print(f'Band {i}: Min={band_min}, Max={band_max}')

# 关闭栅格文件
dataset = None

在这个示例中，我们首先打开要处理的栅格文件。然后，我们使用`RasterCount`属性获取波段数量，并使用循环遍历每个波段。对于每个波段，我们使用`GetRasterBand`方法获取波段对象，并使用`ComputeStatistics`方法计算最小值和最大值。最后，我们输出结果并关闭文件。

需要注意的是，在使用`ComputeStatistics`方法时，我们将参数`approx_ok`设置为`False`，以确保得到准确的结果。

批量裁剪栅格图gdal

要使用GDAL进行批量裁剪栅格图，你可以按照以下步骤进行操作：

1. 导入必要的库和模块，比如`gdal`和`os`。

2. 定义输入栅格数据的路径`input_file`、裁剪后存放的文件夹路径`output_file`以及矢量数据的路径`input_shape`。

3. 获取输入文件夹中所有符合条件的栅格文件，并进行循环处理。

4. 打开栅格数据文件`data`。

5. 使用`gdal.Warp`函数进行裁剪操作，设置输出路径和文件名，以及裁剪的参数，如格式、裁剪边界等。

6. 保存裁剪后的栅格影像。

需要注意的是，使用GDAL进行批量裁剪栅格图的方法可以避免一些arcpy裁剪的缺点，比如不支持中文和对影像重采样导致的偏移。<span class="em">1</span><span class="em">2</span><span class="em">3</span>
#### 引用[.reference_title]
- *1* *2* [gdal矢量裁剪栅格（遥感影像）](https://blog.csdn.net/weixin_44212065/article/details/128006963)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"]
- *3* [批量裁剪栅格 python gdal](https://blog.csdn.net/qq_33339770/article/details/106424980)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"]
[ .reference_list ]