空中旋转篮球-CSDN博客

原创 Python遥感图像处理应用篇(二十七)：Python绘制遥感图像各波段热力图（相关系数矩阵）（续）

续-上一篇中使用csv文件计算的相关系数热力图，本篇我们直接使用遥感图像来计算图像波段之间的相关系数。

2022-12-26 21:25:28 2362 7

原创 Arcgis使用教程（十三）ARCGIS地图制图之地图输出参数设置详解

上图中有两个尺寸大小设置，上面Comments中的设置的纸张大小，下面的Page选项中设置的是页面地图页面大小。我们在使用的时候应该这是后面的页面大小，我们使用File主菜单中的Export Map出图后，发现地图输出结果的范围是后面Page设置大小的范围输出的，所以我们选择后面的进行设置。这里地图我们没有缩放，比例尺还是和原来一样，如果我们将页面设置的很大，然后将地图填充页面，那么地图比例尺就会相应的发生变化。我们修改比例尺后，现有的地图页面大小无法显示完整的地图范围，我们需要增大页面尺寸。

2022-12-24 20:52:00 9115 1

原创 Arcgis使用教程（十二）ARCGIS地图制图之比例尺设置参数详解

在Arcmap中加载中国层数据（中国省级行政区图，国界线两个图层），切换到Layout View窗口，设置地图符号要素并添加比例尺等要素，添加比例尺。设置的页面大小是width：33；height：29.7 默认比例尺：15000000。

2022-12-12 00:02:04 47314

原创 Python遥感图像处理应用篇(二十六)：Python+GDAL Sentinel-2数据波段组合

批量读取原始压缩包解压后的数据，然后进行波段组合运算。文件列表如下：之前写过的是基于Landsat8数据，将所有tif文件放置在一个文件夹种实现的波段组合。文件路径很多，我们需要先读取到需要组合波段的完整路径。选取以上9个波段进行波段组合。最后读取到的是全部的波段完整路径，部分显示如下：接着就可以开始按照上一篇文章之中的代码进行波段组合运算了。

2022-11-24 10:29:36 1375 5

原创 Arcgis使用教程（十一）ARCGIS地图制图之经纬网格设置参数详解

本文详细介绍Arcmap中经纬网格参数设置方法……，中国经纬度范围分布很广，我们间隔默认的是10°。在Arcmap中加载中国矢量图层数据（中国省级行政区图，国界线两个图层），切换到Layout View窗口，设置地图符号要素并添加比例尺等要素，默认添加经纬网。(3)Labels：标签，主要设置网格标签字体风格，大小，间距，方向等参数，这个标签打开速度很慢，需要反映较长时间。设置外部主网格显示标签风格，大小，方向以及设置显示度分秒格式等。即不显示内部网格标签。

2022-11-20 21:59:04 28563 9

原创 Python遥感图像处理应用篇(二十二)：Python+GDAL 批量等距离裁剪影像-续

本篇对上一篇博客代码进行少量修改，实现遥感图像按照指定大小等距裁剪。比如我们需要将原始遥感图像按照60km*60km大小标准块裁剪，我们将上一篇中的裁剪函数参数修改一下，将行列号修改为裁剪大小。

2022-11-04 16:51:31 1765 1

原创 ArcGIS应用（二十六）按照属性分割矢量图层要素为新的图层

将一个矢量图层，按照属性分割为多个独立的矢量图层打开ArcToolbox工具，选择“分析工具”->“提取工具”->“按属性分割”；Input Feature和Split Feature不是同一个图层，Split Field 属性中，每个图形是唯一值，我们直接使用工具即可实现图层的分割。比如我们使用的Input数据是某地区的土地利用数据，我们需要按照行政区去裁剪得到各个行政区的土地利用数据，那么我们可以直接使用该功能进行要素分割得到各个行政区的土地利用数据。

2022-11-02 08:21:23 9981

原创 python rasterio库和gdal库匹配与冲突问题安装办法

近期在对栅格数据进行处理的过程中需要用到rasterio库，在使用pycharm安装的过程中出现了各种问题，pycharm无法安装成功。切换了各种网络环境（包括使用代理）都无法解决。后来还是使用直接下载库文件的方法进行安装。），但是出现了新得问题，存在gdal和rasterio版本匹配与兼容问题。即安装的rasterio版本需要匹配对应的gdal版本，两者如果不兼容，则后安装的那个无法安装成功。原始库文件下载地址如下：raterio库下载地址：包。

2022-10-25 10:59:02 4257 2

原创 Python遥感图像处理应用篇(二十五)：Python+GDAL 波段组合

之前使用arcpy python2.7写了一篇进行遥感数据波段组合计算的博文，使用场景是将所有单波段数据放在一个文件夹中（如下图），文件名称前缀一样，后缀波段编号用来区分不同的波段，数据按照名称默认排序，读取的时候也是按照排列顺序读取。

2022-10-18 22:56:28 2601

原创笔记：错误记录- gdal读取哨兵1雷达数据分贝化数据归一化计算极值为nan，输出相同错误结果

哨兵1雷达数据分贝化数据。对该数据做极值归一化处理的时候，读取的数据最大值和最小值均为nan，归一化的结果都是0值或者无法写入。

2022-10-11 09:44:02 897

原创 ArcGIS应用（二十五）两条直线角度的测算方法

使用ARCGIS软件测算两条之间的夹角是多少。

2022-10-07 17:16:20 4694

原创 ArcGIS应用（二十四）合并多个shapefile文件的几种方法

我们手头有很多个单独的shapefile文件，我们需要将这些文件合并为一个shapefile文件以便于操作和管理。Arcgis软件实现这一目标有一些不同的工具和方法。

2022-09-30 12:09:56 13290

原创 python 经纬度坐标转换为UTM坐标方法与结果验证分析

经纬度坐标是一种地理坐标系，记录的经纬度信息，经度范围是(W)-180~180(E),南半球纬度0~90°S 北半球纬度0~90°N。UTM分带及计算方法，3°分带(高斯克吕格）和6°分带（UTM)带号计算公式如下：Zone = int[(Lon - 1.5)/3] + 1 3°分带Zone = floor(Lon/6)+31 6°分带中央经线计算公式如下：L0 = Zone*3 3°分带根据计算结果L0 = (Zone-1)*3的计算结果和ARCGIS中的分带号一致

2022-09-23 14:25:56 6238 5

原创 ESPA网站中Landsat C2L2数据批量下载网站报错问题500 Internal Server Error

这个问题困扰了我一段时间，期间还向国外网站管理员咨询过（为了咨询，去注册了一个google邮箱，由于网络和google网站注册保护机制等原因，中间问题不断。这也是个技术活，前后花了2-3个小时才弄好），网站管理员回复说我的数据列表有问题，我反复检查了是没有问题的。数据列表txt文件编码问题，我的txt文件默认编码是UTF-16LE,将其另存为utf_8编码保存就好了。

2022-09-22 21:32:01 805

原创 Sentinel-1数据下载与处理

点击上图中的Find Data按钮，搜索Sentinel-1数据点击第二个：检索一下湖北的2020年的Sentinel-1A和Sentinel-1B数据：查看My project，数据量还挺大的：我们在筛选条件里面继续设置格式，选择Geotiff。上图中的Download按钮激活了，可以下载了。数据量太大，我们可以下载自己所需时间的数据，缩小时间范围减少数据。数据的下载先记录到这里，后面再补充数据方法了。3.数据处理3.1软件下载安装桌面如下：正在自动安装一些插件。

2022-09-08 11:45:55 1787

原创错误记录：遥感图像做极值化运算出现：UnboundLocalError: local variable ‘x_nor‘ referenced before assignment

出现以上问题的原因是变量为声明，那么为什么没有声明这个变量，这里是因为if条件语句没有执行，所以x_nor未提前声明。我们这里使用的遥感数据中的最小值大于0值，所以if没有执行。out_ds.GetRasterBand(i).WriteArray(x_nor) # 将每个波段的数据写入内存。

2022-09-07 10:31:34 264

原创分段读取csv文件并可视化处理

使用数据为csv文件，数据有89万条记录，名称字段具体如下：Time (sec)Segment IDDem_heightluccID16126600001612660000161266734190.72271.56288910000161266734190.7330000161266000016126600001612663219390.05000016126600001612660000。

2022-09-04 23:43:27 1080 2

原创 python+matplotlib绘制具有多个子图的图表

图例的位置我们可以采用loc属性设置，该属性包含：loc： supported values are 'best', 'upper right', 'upper left', 'lower left', 'lower right', 'right', 'center left', 'center right', 'lower center', 'upper center', 'center'等位置参数。horizontalalignment设置垂直对齐方式，可选参数：left,right,center。

2022-08-20 22:47:28 3516

原创样本数据异常值处理的三种方法

数据需要服从正态分布，正负3∂的概率是99.7%，那么距离平均值3∂之外的值出现的概率为P(|x-u| 3∂) = 0.003，属于极个别的小概率事件。如果数据不服从正态分布，也可以用远离平均值的多少倍标准差来描述。...

2022-08-12 23:25:43 7280

转载深度学习网络结构图绘制工具及方法

汇总一些深度学习网络结构图的工具，这些软件都是只能作为模型图的辅助，更常见的模型组合型的图还是要ppt来画了，或许还有visio。

2022-08-07 22:50:20 33832 1

原创 office安装出现了“office对安装源的访问被拒绝30068-4(5)”错误

重新安装visio也报错，无法继续安装，这是可能是因为在安装visio过程中临时多次取消引起的问题。下载一个微软的卸载office 程序，再次对计算机中的office残留进行卸载清理。在安装visio2016和Office2016的时候，反复卸载安装后出现安装问题。.........

2022-08-06 14:15:46 4856 2

原创 Python 绘制线性回归散点图和置信区间线2

fromsklearn.model_selectionimporttrain_test_split库实现数据的分割。使用的数据，采用随机创建方法，不同于上一篇，这里的x，y分别是二维数组。和上一篇类似，绘制线性回归方程和置信区间线。

2022-07-17 23:13:54 4680

原创 Python 绘制线性回归散点图和置信区间线

使用python绘制线性回归散点图和回归线以及置信区间线，使用的数据直接随机写两个一维数组，作为x、y值。

2022-07-17 22:32:27 6833

原创 Python遥感图像处理应用篇(二十四)：Python绘制遥感图像各波段热力图（相关系数矩阵）

给多光谱遥感图像各个波段绘制热力图，首先需要计算波段之间的相关系数矩阵，而计算遥感图像波段相关系数矩阵有不同的方法，常用的我们可以采用遥感图像处理软件计算，比如ENVI软件就可以计算相关系数矩阵，使用工具箱中的Statistics工具即可进行多种统计运算。我们这里不使用遥感图像软件计算，直接使用Python计算。要实现整个目标，我们首先需要提取遥感图像波段中各个像元的信息。这里我们可以将遥感图像读取并保存为csv文件，实现这一步我们可以参照之前的一篇文章，链接如下：Python遥感图像处理应用篇(十四)：G

2022-07-13 10:37:21 2297

原创 ArcGIS应用（二十三）矢量图层要素分割方法（基于属性和空间分割）

对矢量图层数据进行要素分割，包括根据属性进行分割以及根据空间位置进行分割两种情况。这里的分割是指将图层要素处理为多个图层要素。比如根据图层属性要素分割，就是按照属性将原始图层分割为不同属性的几个图层。根据空间位置分割是指按照不同的空间范围进行要素的分割，比如这里建立一个渔网数据对要素进行标准分割。使用数据：点图层矢量数据，包含多个属性。使用到的工具Analysis Tools-》Extract-》split：通常情况下，按照属性分割我们字需要在输入要素和分割要素中输入相同的图层，然后选择所使用的的分割属性字

2022-07-08 23:07:01 6293

原创 Python遥感图像处理应用篇(二十三)：Python+GDAL 批量拼接图像

遥感影像拼接分多种不同的情况，比如比较常用的是两幅影像拼接或者多个图像合并为一个图像。GDAL中实现影像拼接的方式也有多种，比如常用可以通过numpy读取影像数组，并计算各个图像的numpy数组范围，合并多个numpy数组为一个最终的结果，并将其重新输出为一幅新的图像，即可实现图像的拼接。另外比较方便的一种方法便是使用GDAL.warp()函数方法。这个函数的功能很多，可以实现多种不同的目的，比如用于坐标系转换、投影变换、图像合并与镶嵌、地理范围裁剪、更改分辨率、矢量裁剪等多种功能。这些功能可以通过设置Op

2022-07-08 22:10:51 4128 5

原创 Python遥感图像处理应用篇(二十二)：Python+GDAL 批量等距离裁剪影像

给定一幅图像，我们需要将图像等距离分割裁剪，分割行列号相同，比如分割裁剪为N行N列的多幅遥感图像。例如原图如下：裁剪为三行三列之后如下：设置不同波段组合为不同颜色以便区分边界。3*3裁剪图像运行环境：运行环境：运行环境：windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。我们主要是用warp是扭曲或者变形的意思，该函数是对一个数据或者多个数据集进行变形，用法如下：参数包括：关键字参数是：如果选项是作为gdal.War

2022-07-04 09:22:22 2257

原创 ArcGIS应用（二十二）Arcmap加载激光雷达las格式数据

使用Icesat-2星载激光雷达ATL03数据，数据格式为las格式。Argis工具箱中的LAS Dataset工具模块，位于Data Management Tools里面，这个工具包含一下工具：Add Files to LAS DatasetCreate LAS Dataset LAS Dataset StatisticsLAS Point Statistics as RasterRemove Files From LAS Dataset使用Create LAS Dataset 工具先转件一个LAS 数据

2022-07-03 22:50:14 4340

原创 Python遥感图像处理应用篇(二十一)：Python+GDAL 批量计算遥感图像NDVI指数

之前也写过NDVI批量计算的实现方法，采用的是Acrpy python2.7实现的，这里我们采用GDAL实现NDVI指数的批量计算。使用数据，多波段遥感图像数据。Landsat05 Colection2 Level2数据。NDVI=(NIR-RED)/(NIR+RED) 就算之前需要搞清楚多光谱数据中的波段对应关系。影像数据地表反射率修正系数：0.0000275+（-0.2）运行环境：运行环境：windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-w

2022-07-03 20:34:48 2139

原创 ArcGIS应用（二十一）Arcmap删除图层指定要素的方法

针对Arcmap删除矢量图层指定条件要素的问题，对于一个GISer来说，首先想到的应该就是使用Arcmap中编辑工具进行处理。即启动编辑工具之后，采用属性选择工具，选择要删除的要素，然后直接Delete之后保存数据即可，这是最直接和便捷的方法。当然还有很多其他的方法，比如采用AE编程实现要素删除，可以定制自己需要的各种功能，也可以是使用Python库-Gdal编程实现要素删除，比如这篇：python批量删除shapefile文件中某字段的指定属性值_空中旋转篮球的博客-CSDN博客https://soder

2022-07-02 22:46:37 24205 3

原创 Python GDAL+numpy遥感图像处理过程中背景像元处理方法

我们在处理遥感图像过程中，经常要涉及到背景像元的处理，遥感影像通常采用规则矩形数组存储，而有效影像区域往往是不规则数据，这就涉及到背景像元值的处理问题。大多软件对于遥感图像背景值处理是记为0值，也有的是-9999、-19999、nodata、nan等。有些计算，背景值不影响计算结果，这里的计算结果是指影像区域的计算结果，背景值往往还是会产生变化的，但是有些变化是非0值计算为0值了，所以我们打开影像的时候显示也没多大变化。有些计算背景值会对计算产生错误结果，后面将详细介绍。打个比方，以影像背景值为-19999

2022-06-29 13:35:33 1712

原创 Python遥感图像处理应用篇(二十)：Python+GDAL 批量提取多波段图像为单波段图像

1.实现目标批量将多波段图像提取为各个单波段图像运行环境：windows10 pycharmpython3.7.7GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。2.实现代码"""此段代码将遥感图像背景值去除处理"""import numpy as npfrom osgeo import gdal, gdalconstimport os#将遥感影像归一化处理写成函数def GetEnvolopePoint(inputpath,outp...

2022-06-26 22:31:48 2233 4

原创 Python+gdal裁剪遥感图像出现“Warning 1: TIFFReadDirectory:Sum of Photometric type-related color channels ……”

运行环境：windows10、pycharm、python3.7、gdal等python使用GDAL 裁剪遥感图像的时候出现以下警告：“Warning 1: TIFFReadDirectory:Sum of Photometric type-related color channels and ExtraSamples doesn't match SamplesPerPixel. Defining non-color channels as ExtraSamples.”使用的代码段：Warp函数有很多...

2022-06-24 09:22:56 7793 3

原创 Python遥感图像处理应用篇(十九)：GDAL +numpy批量对遥感图像外围背景值进行处理

最近下载了一些遥感影像数据，这些数据都包含大量的外围背景数据，如下图所示：外围背景值都为0值。本文描述的是采用python批量处理外围背景，只保留最小外围背景区域。如下图：基本思路如下：首先获取遥感影像有数据区域四个角的坐标，这里的坐标我们用行列号表示，即左上角（x1,y1),右上角（x2,y2），左下角（x3,y3），右下角（x4,y4）。我们需要通过遍历遥感数据的方法获取到四个角点的行列号。然后，根据四角坐标计算四至范围，即left，top，right，bottom值。之后，分别获取左上角坐标和右下角坐

2022-06-23 17:51:29 1288 2

原创 ArcGIS应用（二十）Arcgis 栅格图像符号系统提示“This dataset does not have valid histogram required for classificati…”

Arcgis打开单波段的遥感图像进行符号色设置，选择Classified选项的时候出现以下问题：打开工具Calculates statistics工具，位置在Data Management Tools->Raster->Raster Properties->Calculates statistics.打开图像计算一下：参数都默认，直接点击确定。忽略值可以选0.计算完之后，重新打开数据，设置符号：可以正常选择Claasified选项了。这样就可以自由设计符号系统了。...

2022-06-22 14:54:12 1475

原创 Python遥感图像处理应用篇(十八)：GDAL +numpy对遥感图像归一化处理

所用输入数据为遥感影像数据，可以是单波段也可以是多波段数据。运行环境：windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。完整代码如下：采用批处理方式，输入和输出为文件路径。只需要修改输入输出文件路径即可计算所需结果。3.数据运行测试与结果比如采用Landsat8中的NDVI数据进行批量计算. 图像对比：图像数值范围：图像背景值变化：原图是nodata计算之后背景值出现变化，变为0值。像元值查看对比：..

2022-06-21 11:01:31 2288 4

原创 Python遥感图像处理应用篇(十七)：GDAL 将归一化处理csv数据转化为多波段遥感影像

本文所使用数据和上一篇数据相同，上一篇链接：Python遥感图像处理应用篇(十六)：GDAL 将归一化处理csv数据转化为遥感影像_空中旋转篮球的博客-CSDN博客上一篇文章将csv文件转化为单波段图像，这一篇将csv文件合成为一个多波段图像。注意一个前提是需要知道原图像的基本信息，如数据行列号，投影信息，数据类型信息等。实现写入多光谱数据的时候，核心问题在于使用循环写入各个通道数据。这里直接写入的是二维数组数据，看别人的文章的时候，有的使用的三维数组，如格式（channel,row,col）,我在使用的时

2022-06-19 18:28:06 1029

原创 Python遥感图像处理应用篇(十六)：GDAL 将归一化处理csv数据转化为遥感影像

将上一篇文章中得到的计算结果作为转换数据。链接如下：Python遥感图像处理应用篇(十五)：GDAL 读取多光谱数据做归一化处理_空中旋转篮球的博客-CSDN博客基本思路，先读取csv数据，现读取每一列数据为一维数组，然后根据原始影像数据的大小（row，col）将一维数组塑形，即使用reshape函数变换为（row，col）形状对的二维数组。使用gdal将二维数组转换为tif图像。完整代码如下：修改一下以上代码，可以将数据转化为和原始数据相同波段数的多波段遥感影像。可以使用for循环一个通道一个通道的

2022-06-16 15:49:30 1549 1

原创 Python遥感图像处理应用篇(十五)：GDAL 读取多光谱数据做归一化处理

上一篇文中计算得到的csv文件：Python遥感图像处理应用篇(十四)：GDAL 读取多光谱数据为二维数组并存入csv文件_空中旋转篮球的博客-CSDN博客我们在做归一化处理的时候，遥感软件往往是针对单波段图像的，或者不方便批处理图像。我们可以使用python自己实现计算处理。首先将数据读取为numpy数组，然后根据数组运算，计算得到归一化计算结果；这里采用（X-Xmin）/（Xmax-Xmin）计算公式归一化处理。计算完整代码如下： 2.3 批量计算代码将以上代码修改为批量计算代码，可以一次性

2022-06-14 11:33:14 1881 1

原创 Python遥感图像处理应用篇(十四)：GDAL 读取多光谱数据为二维数组并存入csv文件

之前写过一遍文章是采用Arcgis读取多光谱数据为csv文件的，相对来说有一些麻烦，本文有个更简单一点的方法，使用python直接读取。本文使用Landsat05数据，采用GEE下载的多光谱数据，该数据下载下来包含19个波段。这些波段信息可以查看数据该数据详细说明。1-7是其对应的多光谱波段，第六波段是地表温度数据（ST),其他的是(SR) 地表反射率产品数据。.....................

2022-06-13 12:07:06 2363 10

ICESat2_ATL03_ATBD_r004.pdf

ICESAT-2 ATLAS 产品说明,比如介绍各级数据级处理算法和理论基础。该产品由ICESat-2任务生产。它包括参数的描述和这些ICESat-2参数的推导。在ICESat-2科学研究者领导的处理系统(SIPS)上实现的待定bd处理的算法理论基础文件。 sip支持在ICESat- 2空间飞船上的ATLAS(先进地形激光高度计系统)仪器，包括ATLAS科学算法软件(ASAS)和调度和数据管理系统(SDMS)。科学算法软件会生成0级到0级 4级标准数据产品以及相关产品质量评估和元数据信息等。

2021-12-23