第9章 读写栅格数据

本章内容：
（1）理解栅格数据的基本概念
（2）GDAL介绍
（3）读写栅格
（4）数据重采样

9.1 栅格数据介绍

栅格数据集： 由波段组成，每个波段都是一个二维数组，波段的集合组成三维数组。栅格数据集可以描述具有连续变化的数据，如高程、坡度、温度等。当然，栅格数据集不必包含连续的数据，如土地覆盖分类——由栅格数据制作而成。
栅格数据集的空间参考系统： 栅格数据集通过一组坐标、像素的大小和数据集的旋转量，来确定图像其余部分的坐标（也称为仿射变换，是栅格数据集几何校正的一种）。该组坐标一般为图像左上角，称为原点。

栅格数据集的内存： 栅格数据集的所占内存大小取决于像素大小和影像大小。因此，在有限的CPU中，往往会采用概视图层（也可称为金字塔涂层或递减分辨率）的方法来缓解受限的CPU内存。
概视图层： 通过递减分辨率的方法获得与原始影像面积相同，但重采样后的像素尺寸更大。当在相同内存的情况下，1）全局查看整图时，将粗糙的显示全图，此时像素尺寸很大；2）精细查看时，只显示部分区域，但像素尺寸是原始大小，可以看到更精细的细节。
栅格数据压缩： 有损压缩，会丢失部分信息，如.jpeg。无损压缩，不丢失分和信息，课重建为原始未压缩数据，如.png和.GeoTIFF。
栅格重采样： 栅格重采样到大小不同的格网单元或重投影到另一个空间参考系统时，没有像素一对一映射，因此需要计算新的像素值。

9.2 GDAL介绍

例子：利用GDAL将独立的栅格波段合成一张图像

###------代码清单9.1 将独立的栅格波段合成一张图像------

import os
from osgeo import gdal

os.chdir(r'/home/daidai/bigdata/Python-for-GIS/osgeopy-data/Landsat/Washington/')
band1_fn = &