遥感概论
遥感是探测和监测一个地区的物理特征的过程，方法是测量该地区在一定距离(通常是从卫星或飞机)反射和发射的辐射。特殊摄像机收集远程遥感图像，这有助于研究人员“感知”地球上的事物。
由带电粒子振动产生的电磁能量，以波的形式在大气和真空空间中传播。这些波有不同的波长和频率，波长越短，频率就越高。有些电波、微波和红外线的波长较长。而紫外线、x射线和伽马射线的波长要短得多。可见光位于长到短波辐射范围的中间。这一小部分能量是人类肉眼所能探测到的。需要仪器来探测所有其他形式的电磁能。在不同卫星的帮助下，我们利用全范围的光谱来探索和理解发生在地球和其他行星上的过程。
一些例子是:
- 卫星和飞机上的相机拍摄了地球表面大面积的图像，使我们看到的东西比我们站在地面上看到的多得多。
- 船上的声纳系统可以在不需要潜入海底的情况下生成海底图像。
- 卫星上的照相机可以用来拍摄海洋温度变化的图像。
地球遥感图像的一些具体用途包括:
- 大型森林火灾可以从太空中绘制地图，使得护林员可以看到比地面大得多的区域。
- 追踪云层来帮助预测天气或观察喷发的火山，并帮助观察沙尘暴。
- 在数年或数十年的时间里跟踪城市的发展和农田或森林的变化。
- 发现并绘制海底崎岖的地形(例如巨大的山脉、深谷和海底的“磁条”)。
在本文中，我们将使用Sentinel-2卫星收集的数据。以下是简要信息:
Sentinel-2任务由两颗卫星组成，用于支持植被、土地覆盖和环境监测。2015年6月23日，欧洲航天局(ESA)发射了Sentinel-2A卫星，在太阳同步轨道上运行，周期为10天。第二颗相同的卫星Sentinel-2B于2017年3月7日发射，并可在EarthExplorer上获取数据。它们每隔五天就会覆盖地球所有的陆地表面、大型岛屿以及内陆和沿海水域。
Sentinel-2多光谱仪器(MSI)沿着290公里的轨道带获得了从可见和近红外(VNIR)到短波红外(SWIR)波长的13个光谱波段。MSI测量了在13个光谱波段内通过大气的反射辐射。空间分辨率取决于特定的光谱波段:
10米波段:蓝色(490 nm)，绿色(560 nm)，红色(665 nm)，近红外(842 nm)。在20米处有6个波段:4个窄波段用于植被表征(705 nm、740 nm、783 nm和865 nm)和2个较大的SWIR波段(1610 nm和2190 nm)用于雪/冰/云探测或植被水分胁迫评估。60米波段:主要用于云层遮挡和大气校正(气溶胶为443 nm，水蒸气为945 nm，卷云探测为1375 nm)。
土地覆盖分类的深度学习
深度学习。神经网络。反向传播。在过去的一两年里，我经常听到这些流行词，最近它也确实抓住了我的好奇心。深度学习是当今一个活跃的研究领域，如果你一直关注计算机科学领域，我相信你至少一次遇到过其中的一些术语。
深度学习可能是一个令人生畏的概念，但如今它变得越来越重要。谷歌已经通过谷歌Brain项目和最近收购伦敦的深度学习初创公司DeepMind在这一领域取得了巨大进展。此外，深度学习方法几乎在每一个指标上都击败了传统的机器学习方法。
深度神经网络(Deep Neural network, DNNs)和卷积神经网络(Convolutional Neural network, CNNs)是目前广泛应用于土地覆盖分类的方法。
深度神经网络(DNN)
一般来说，神经网络是一种用来模拟人类大脑活动的技术——具体来说，就是模式识别和输入通过模拟神经连接的各个层次。
许多专家将深度神经网络定义为具有输入层、输出层和中间至少有一层隐藏层的网络。每一层在一个被称为“特性层次结构”的过程中执行特定类型的排序和排序。“这些复杂的神经网络的关键用途之一是处理未标记或非结构化数据。“深度学习”一词也被用来描述这些深度神经网络，因为深度学习代表了一种特定形式的机器学习，在机器学习中，技术使用人工智能方面寻求以超越简单输入/输出协议的方式对信息进行分类和排序。
本文对实现高光谱图像的土地覆盖分类进行了验证。
卷积神经网络(CNN)
根据《CS231n:用于视觉识别的卷积神经网络》，卷积神经网络与普通神经网络非常相似:它们由具有可学习权值和偏差的神经元组成。每个神经元接收一些输入，执行点积并有选择地进行非线性操作（激活）。整个网络仍然表示一个单一的可微分数函数:从一端的原始图像像素到另一端的分类分数。它们在最后一层(完全连接)上仍然有一个损失函数(例如SVM/Softmax)，我们为学习规则神经网络开发的所有技巧/技巧仍然适用。
改变什么呢?ConvNet架构明确假设输入是图像，这允许我们将某些属性编码到架构中。这样可以使转发功能更有效地实现，并极大地减少网络中的参数数量。
池化或下采样层负责减少激活映射的空间大小。通常，在其他层(即卷积层和非线性层)的多个阶段之后使用它们，通过网络逐步减少计算需求，并最大限度地减少过拟合的可能性。
卷积网络中完全连接的层实际上是一个多层感知器(通常是两层或三层MLP)，其目标是将先前不同层的组合激活量映射到一个类概率分布中。
用于卫星影像土地覆盖分类的架构有AlexNet、ResNet、Inception、U-Net等。
Sundarbans国家公园卫星图像
孙德尔本斯（Sundarbans）是孟加拉湾恒河、雅鲁藏布江和梅克纳河汇流形成的三角洲中最大的红树林地区之一。孙德尔本斯森林横跨印度和孟加拉国，面积约为1万平方公里，其中40%位于印度境内，是许多珍稀和全球濒危野生物种的家园。在本文中，我们将使用2020年1月27日Sentinel-2卫星获取的孙德本斯卫星数据的一部分。
让我们开始编码。
读取数据
让我们使用rasterio读取12个波段，并使用numpy.stack()方法将它们堆叠成一个n维数组。叠加后得到的数据形状为(12,954,298)。使用loadmat方法从scipy.io包中读取卫星图像的地面真值。真值包含6类，包括水，植物，树木，土地等。
fromglobimportglobimportnumpyasnpfromscipy.ioimportloadmatimportrasterioasrioS_sentinel_bands=glob("/content/drive/MyDrive/Satellite_data/sundarbans_data/*B?*.tiff") S_sentinel_bands.sort() l= [] foriinS_sentinel_bands: withrio.open(i, 'r') asf: l.append(f.read(1)) #Dataarr_st=np.stack(l) #GroundTruthy_data=loadmat('Sundarbands_gt.mat')['gt']
数据可视化
这些数据具有多个波段，其中包含从可见光到红外光的数据。因此，很难对人类的数据进行可视化。通过创建RGB复合图像，可以更轻松地有效理解数据。要绘制RGB复合图像，您将绘制红色，绿色和蓝色波段，分别是波段4、3和2。由于Python使用从零开始的索引系统，因此您需要从每个索引中减去1的值。因此，红色带的索引为3，绿色为2，蓝色为1。
如果像素亮度值偏向零值，则我们创建的合成图像有时可能会很暗。可以通过使用参数Stretch = True拉伸图像中的像素亮度值以将这些值扩展到电位值的整个0-255范围以增加图像的视觉对比度来解决此类问题。另外，str_clip参数允许您指定要裁剪的数据尾部的数量。数字越大，数据将被拉伸或变亮的越多。
我们来看一下绘制RGB复合图像以及所应用的拉伸的代码。
ep.plot_rgb( arr_st, rgb=(3, 2, 1), stretch=True, str_clip=0.02, figsize=(12, 16), #title="RGB Composite Image with Stretch Applied", ) plt.show()
让我们使用eathpy.plot包中的plot_bands方法可视化。
ep.plot_bands(y_data, cmap=ListedColormap(['darkgreen', 'green', 'black', '#CA6F1E', 'navy', 'forestgreen'])) plt.show()
下图显示了Sundarbans卫星数据的合成图像和地面真实情况。
正如我们所讨论的，数据包含12个波段。让我们使用EarhPy套件可视化每个波段。plot_bands（）方法获取带和图的堆栈以及自定义标题，这可以通过使用title =参数 将每个图像的唯一标题作为标题列表传递来完成。
ep.plot_bands(arr_st, cmap='gist_earth', figsize= (20, 12), cols=6, cbar=False) plt.show()
