27、利用NumPy和SciPy进行地图代数运算

tree

于 2025-11-24 12:58:00 发布

阅读量5

点赞数

CC 4.0 BY-SA版权

分类专栏： Python地理处理入门文章标签： NumPy SciPy 地图代数

本文链接：https://blog.youkuaiyun.com/tree/article/details/155434313

Python地理处理入门专栏收录该内容

34 篇文章 ¥499.90

订阅专栏¥69.90

会员秒杀 ¥9.9 重磅福利

超级会员免费看

利用NumPy和SciPy进行地图代数运算

1. 地图代数运算基础代码

以下是一段地图代数运算的基础代码，用于处理栅格数据：

xsize = in_band.XSize
ysize = in_band.YSize
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create(out_fn, xsize, ysize, 1, gdal.GDT_Int32)
out_ds.SetProjection(in_ds.GetProjection())
out_ds.SetGeoTransform(in_ds.GetGeoTransform())
out_band = out_ds.GetRasterBand(1)
out_band.SetNoDataValue(-99)
n = 100
for i in range(0, ysize, n):
    if i + n + 1 < ysize:                                              
        rows = n + 2                                                     
    else:                                                               
        rows = ysize - i                                                
    yoff = max(0, i - 1)