本文介绍基于R语言中的raster包，读取单张或批量读取多张栅格图像，并对栅格图像数据加以基本处理的方法。 1包的安装与导入   首先，我们需要配置好对应的R语言包；前面也提到，我们这里选择基于raster包来实现栅格图像数据的读取与处理工作。首先，如果有需要的话，我们可以先到raster包在R语言的官方网站中，查阅raster包的基本情况，比如其作者信息、当前的版本、所依赖的其他包等等；如下图所示。   当然，这些内容看不看都不影响我们接下来的操作。接下来，我们开始安装raster包；这里我是在RStudio中进行代码的撰写的。   首先，我们输入如下的代码，从而开始raster包的下...

  本文介绍基于ENVI软件，对不含有任何地理参考信息的栅格遥感影像添加地理坐标系或投影坐标系等地理参考信息的方法。   我们先来看一下本文需要实现的需求。现有以下两景遥感影像，其位于不同的空间位置；但由于二者均不含任何地理参考信息，导致其在ENVI软件中打开后会自动重叠在一起；如下图所示。   那么我们就以其中一景遥感影像为例，对其添加地理参考信息。   明确了具体需求，接下来就可以开始操作。首先，我们在ENVI软件中打开对应的两景遥感影像；其次，在需要添加地理参考信息的图像名称处右键，选择“ViewMetadata”。   弹出如下所示的元数据浏览窗口。   这里我们需要注意：如果大...

  本文介绍基于ENVI软件，利用“PixelBasedMosaicking”工具实现栅格遥感影像镶嵌拼接的方法。   首先需要说明的是，本文需要镶嵌的遥感影像并不含地理参考信息，因此仅可以使用ENVI中的“PixelBasedMosaicking”工具（该工具可以对含有或不含有地理参考信息的图像进行镶嵌），而不是更为先进的“SeamlessMosaic”工具（该工具仅可对含有地理参考信息的图像进行镶嵌）；针对后者，大家可以查看我们后续的博客。   我们先来看一下本文需要实现的需求。现有以下两景遥感影像，在ArcMap软件中打开，其中一景如下图所示。   另一景则如下图所示，可以看到两景遥感...

  本文介绍Landsat7遥感影像数据中B61、B62两个热红外波段的区别，以及研究应用时二者选择的依据。   Landsat7遥感影像数据具有2个热红外波段，分别是Band61与Band62这两个波段；有些地方也将这两个波段写作Band6L与Band6H、Band6a与Band6b等等。因此，这两个波段的具体区别，以及在实际应用时究竟该选择哪一个波段，成为一个需要考虑的问题。   首先，关于这两个波段的区别，其数据官网已经给出了具体的解释，如下图所示。   简单来说，这两个热红外波段就是由同一个遥感传感器（也就是增强型专题制图仪，ETM+）获取得到的两个具有不同增益（Gain）的热红外波...

  本文介绍在C语言中，矩阵库Armadillo的mat、vec格式数据与计算机视觉库OpenCV的Mat格式数据相互转换的方法。   在C语言的矩阵库Armadillo与计算机视觉库OpenCV中，都有矩阵格式的数据类型；而这两个库在运行能力方面各有千秋，因此实际应用过程中，难免会遇到需要将二者的矩阵格式数据类型加以相互转换的情况。本文就对其相互转换的具体方法加以介绍。   首先，二者相互转换需要用到的代码如下。 include<iostream> include<armadillo> include<opencv2/opencv.hpp> usingn...

  本文介绍在VisualStudio2022中配置、编译C计算机视觉库OpenCV的方法。 1OpenCV库配置   首先，我们进行OpenCV库的下载与安装。作为一个开源的库，我们直接在其官方下载网站中进行下载即可；如下图所示，我们首先选择需要下载的操作系统。   随后，即可在弹出的新界面中自动开始OpenCV库的下载。   下载完毕后，可以得到OpenCV库的.exe格式文件。   我们双击这一.exe格式文件，即可开始OpenCV库的文件提取过程，也就相当于是安装过程；其中，我们首先需要选择OpenCV库安装的路径。   随后，即可开始OpenCV库的文件提取。   Open...

  本文介绍基于Python中ArcPy模块，实现基于栅格图像批量裁剪栅格图像，同时对齐各个栅格图像的空间范围，统一其各自行数与列数的方法。   首先明确一下我们的需求。现有某一地区的多张栅格遥感影像，其虽然都大致对应着同样的地物范围，但不同栅格影像之间的空间范围、行数与列数、像元的位置等都不完全一致；例如，某一景栅格影像会比其他栅格影像多出一行，而另一景栅格影像可能又会比其他栅格影像少一列等等。我们希望可以以其中某一景栅格影像为标准，将全部的栅格影像的具体范围、行数、列数等加以统一。   本文所用到的具体代码如下。 -coding:utf-8-- """ CreatedonThuDec292...

  本文介绍R语言及其集成开发环境RStudio的下载、安装方法。   R语言是一个属于GNU操作系统的开源软件，在数据统计与分析、可视化等方面具有优秀的表现；而RStudio则是R语言的集成开发环境（IDE），可以帮助我们更好地编辑、调试R语言的代码。这二者的关系有点类似于Python与Spyder的关系——我们可以只下载R语言，用其自带的原生编辑窗口来完成代码的撰写与运行等工作；而如果想提高代码的撰写、调试效率，就可以通过RStudio来完成。 1R语言的下载与安装   首先，我们进行R语言的下载。我们可以在TheComprehensiveRArchiveNetwork（CRAN），也就是...

  本文介绍在Win10电脑中，安装Anaconda环境与Python语言的方法。   在这里需要注意，本文介绍的方法是在电脑自身原本不含有Python的情况下进行的；如果大家电脑中原本就下载、安装过Python语言，需要首先将原本安装Python时的设置的环境变量删除。   首先，先进行Anaconda环境的安装。   我们可以在Anaconda的官方网站进行下载，但是官网下载时可能网速较慢。   因此，我们可以在清华大学开源软件镜像站进行下载。   下载后，双击打开下载好的.exe文件，即可开始安装。   其中，在“SelectInstallationType”一栏中，建议大家选择“...

  本文介绍基于Python中gdal模块，对大量多时相栅格图像，批量绘制像元时间序列折线图的方法。   首先，明确一下本文需要实现的需求：现有三个文件夹，其中第一个文件夹存放了某一研究区域原始的多时相栅格遥感影像数据（每一景遥感影像对应一个时相，文件夹中有多景遥感影像），每一景遥感影像都是.tif格式；第二个文件夹与第三个文件夹则分别存放了前述第一个文件夹中原始遥感影像基于2种不同滤波方法处理后的遥感影像（同样是每一景遥感影像对应一个时相，文件夹中有多景遥感影像），每一景遥感影像同样也都是.tif格式。我们希望分别针对这三个文件夹中的多张遥感影像数据，随机绘制部分像元对应的时间序列曲线图（每...

  本文介绍基于C语言GDAL库，为CreateCopy()函数创建的栅格图像添加更多波段的方法。   在C语言的GDAL库中，我们可以基于CreateCopy()函数与Create()函数创建新的栅格图像文件。其中，CreateCopy()函数需要基于一个已有的栅格图像文件作为模板，将模板文件的各项属性信息（例如空间参考信息、像元个数、像元大小、波段数量等），自动作为新创建的栅格图像文件的属性信息；而Create()函数则是仅仅新建立一个栅格图像，需要我们自行定义新栅格图像的各类属性信息。   因此，一般我们选择CreateCopy()函数来创建栅格图像文件较为方便，因为其不需要我们手动为所...

  本文介绍VisualStudio2022软件Community（社区版）的下载、安装、运行与使用方法。   首先需要提一句，本文介绍的是VisualStudio2022软件的下载；而其它版本的VisualStudio软件（例如VisualStudio2019软件）下载方法与此大同小异，都可以参考本文所述的步骤。   首先，我们在软件官网中，进行VisualStudio2022软件的下载。VisualStudio2022软件一共有Community（社区版）、Professional（专业版）与Enterprise（企业版）三个版本。其中，社区版是免费的，主要面向学生、开源软件开发者等个人；...

  本文介绍查看VisualStudio软件_MSC_VER值的方法。   _MSC_VER是微软公司推出的C/C编译器——MSVC编译器的一个内置宏，其值表示当前VisualStudio软件中MSVC编译器的具体版本。不同的VisualStudio软件版本对应着不同的MSVC编译器版本——无论是不同发布年份的版本（例如VisualStudio2017与VisualStudio2019），还是同一发布年份的不同版本（例如VisualStudio2017version15.8与VisualStudio2017version15.9），其MSVC编译器版本都有差异，因此其_MSC_VER值各不相同...

  本文介绍在VisualStudio软件中配置C环境下线性代数运算库Armadillo的方法。   首先，我们需要在Armadillo库官网下载其源代码，直接点击下图所示红色框内部分即可。   点击上图所示位置后，将弹出一个新的下载界面；Armadillo库的源代码将随后自动下载。   接下来，我们在VisualStudio软件中新建一个项目。   其中，模板选择“空项目”。   项目的名称与存储位置大家可以自行设定，但存储路径建议选择在某个盘符下的第一个子文件夹中（即路径不要设置的太深即可）。   随后，将本文开头我们下载好的Armadillo库的源代码解压到上述新创建项目的存储...

  本文介绍在VisualStudio软件中调用C各种配置、编译完毕的第三方库的方法。   在撰写C代码时，如果需要用到他人撰写的第三方库（例如地理数据处理库GDAL、矩阵运算库Armadillo等），并不能像Python等语言那样，安装好库后直接在不同代码文件中使用；而是需要每一次新建一个代码文件（在VisualStudio软件中就是新建一个项目）后，还需要额外配置一下才可以使用。   本文就以之前的文章VisualStudio部署C矩阵库Armadillo的方法中介绍的矩阵运算库Armadillo为例，介绍安装完某一个第三方库后，如何在VisualStudio软件新的项目中调用这个库。  ...

  本文介绍在VisualStudio软件中配置、编译C环境下GDAL库、SQLite环境与PROJ库的详细方法。   GDAL库是一个非常方便的地理数据处理库，但其在C环境下的配置与编译流程较为复杂；尤其是最新的GDAL3及以上版本，其在C环境中的配置更是首先需要满足许多其他的环境配置条件（包括SQLite环境、CMake软件与PROJ库），使得其配置变得更复杂一些。本文就详细介绍在VisualStudio软件中，配置GDAL库所需的环境，并最终配置、编译GDAL库的方法。 目录 1写在前面 2具体配置 2.1SQLite配置 2.2CMake配置 2.3PROJ配置 2.3....

  本文介绍基于C语言，遍历文件夹中的全部文件，并从中获取指定类型的文件的方法。   首先，我们来明确一下本文所需实现的需求。现在有一个文件夹，其中包含了很多文件，如下图所示；我们如果想获取其中所有类型为.bmp格式的文件的名称，如果文件数量比较多的话，手动筛选就会很麻烦。而借助C代码就可以简单地实现这一需求。如果需要借助Python代码来实现同样的需求，可以参考文章PythonArcPy批量掩膜、重采样大量遥感影像，基于其中提到的arcpy.ListRasters()函数来实现。   首先需要说明的是，本文代码只能实现对某一文件夹下的文件进行遍历并筛选；如果是当前文件夹下的子文件夹中的文件...

  本文介绍基于C语言GDAL库，批量读取大量栅格遥感影像文件，并生成各像元数值的时间序列数组的方法。   首先，我们来明确一下本文所需实现的需求。现在有一个文件夹，其中包含了很多不同格式的文件，如下图所示。   其中，我们首先需要遍历这一文件夹，遴选出其中所有类型为.bmp格式的栅格遥感影像文件（一共有6个），并分别读取文件（已知这些遥感影像的行数、列数都是一致的）；随后，将不同遥感影像的同一个位置的像素的数值进行分别读取，并存储在一个数组中。例如，最终我们生成的第一个数组，其中共有6个元素，分别就是上图所示文件夹中6景遥感影像各自(0,0)位置的像元数值；生成的第二个数组，其中也是6个元...

  本文介绍基于C语言GDAL库，批量创建大量栅格遥感影像文件，并将数据批量写入其中的方法。   首先，我们来明确一下本文所需实现的需求。已知我们对大量遥感影像进行了批量读取与数据处理操作——具体过程可以参考文章CGDAL提取多时相遥感影像中像素随时间变化的数值数组；而随后，就需要对我们处理后的栅格数据再进行输出，即建立新的大量的栅格遥感影像，并将我们处理后的像元数据依次输入进去。   明确了具体需求，接下来就可以开始代码的实践；本文所用到的具体代码如下。这里需要注意，在这里就仅将与本文需求有关的代码放了上来，其他无关的代码就省略了（所以以下代码只是程序主函数中的一部分）；大家在实践过程中，依...

  本文介绍在VisualStudio软件中配置、编译C环境下matplotlibcpp库的详细方法。   matplotlibcpp库是一个C环境下的绘图工具，其通过调用Python接口，实现在C代码中通过matplotlib库的命令绘制各类图像。由于其需要调用Python接口，因此在配置matplotlibcpp库时有些较为麻烦的操作。本文就将matplotlibcpp库的具体配置方法进行详细介绍。 1Git配置   Git是一个分布式开源版本控制系统，在后期我们需要基于其完成vcpkg包管理器的下载与安装，因此需要首先完成Git的配置；具体方法大家可以参考下载、安装Git并拷贝GitHu...