gdal栅格数据处理注意事项-备忘录

• gdal读取栅格数据，其中行列号和返回的数组索引是相反的，所以赋值的时候一定是行列互换，否则读取的时候就会越界

  array_band = raster_band.ReadAsArray(0, 0, raster_ds.RasterXSize, raster_ds.RasterYSize).astype(np.float64)
  
  for i in range(x):
      for j in range(y):
          array_band[y-1, x-1] = 100.0

• gdal中SetNoDataValue()函数说明
  
  例如：band.SetNoDataValue(-32767)，表示在tiff文件中，所有像素值为-32767栅格均默认为该栅格处的像素属于无数据的含义