C#本程序是实现的单层栅格数据的演示，若改为多层（波段，band）栅格数据，数据的读取（存储变量设置）数据的输出该如何改进。

要实现多层栅格数据的读取和存储，可以使用GDAL库来处理栅格数据。GDAL是一个开源的栅格数据处理库，可以读取和写入各种常见的栅格数据格式。

在C#中可以使用GDAL的.NET封装库：GDAL.NET来进行栅格数据的读写。下面是一个简单的示例代码，演示如何读取一个多波段（多层）的栅格数据：

using OSGeo.GDAL;

// 打开栅格数据
Dataset ds = Gdal.Open("path/to/raster.tif", Access.GA_ReadOnly);

// 获取栅格数据的宽度和高度
int width = ds.RasterXSize;
int height = ds.RasterYSize;

// 获取栅格数据的波段数（层数）
int bandCount = ds.RasterCount;

// 依次读取每个波段
for (int i = 1; i <= bandCount; i++)
{
    // 获取当前波段
    Band band = ds.GetRasterBand(i);

    // 定义一个数组来存储当前波段的像素值
    float[] buffer = new float[width * height];

    // 读取当前波段的像素值
    band.ReadRaster(0, 0, width, height, buffer, width, height, 0, 0);

    // 输出当前波段的像素值
    Console.WriteLine("Band {0}:", i);
    for (int j = 0; j < width * height; j++)
    {
        Console.Write("{0} ", buffer[j]);
    }
    Console.WriteLine();
}

// 关闭栅格数据
ds.Dispose();

在上面的代码中，我们首先打开了一个栅格数据文件，并获取了宽度、高度、波段数等信息。然后，我们使用一个循环来依次读取每个波段的像素值，并输出到控制台上。

如果你需要将多层栅格数据存储到变量中，可以使用一个多维数组来存储每个波段的像素值。例如，可以定义一个三维数组来存储一个三层的栅格数据：

float[,,] data = new float[width, height, 3];

for (int i = 0; i < 3; i++)
{
    Band band = ds.GetRasterBand(i + 1);

    float[] buffer = new float[width * height];
    band.ReadRaster(0, 0, width, height, buffer, width, height, 0, 0);

    for (int j = 0; j < width * height; j++)
    {
        data[j % width, j / width, i] = buffer[j];
    }
}

在上面的代码中，我们使用一个三维数组来存储三个波段的像素值。我们首先定义一个三维数组，然后使用一个循环来依次读取每个波段的像素值，并将其存储到对应的数组位置中。注意，在存储像素值时，我们使用了一个简单的算法来计算像素在数组中的位置：第一维是列数，第二维是行数，第三维是波段数。