KilllerQueen-CSDN博客

遥感影像叠加，栅格图像叠加，不同尺寸和不同波段数背景高分四号卫星搭载两个不同的传感器，获取的PMS数据包含5个波段，IRS只有1个中红外波段，他们的尺寸和波段数量是不同的，IRS数据包含在PMS数据内部。叠加效果如下图：思路考虑到它们的左上角并不是对齐的，所以想要实现波段叠加，不能直接将它们的数据矩阵叠加，否则必然错位。而是要找到IRS左上角在PMS中对应的位置（x，y），然后再写到同PMS一样大的array中。代码实现def layerstaking(PMS,IRS,merge):

2022-05-30 14:23:54 972 3

原创 python+GDAL给图像设置地理参考（GeoTransform和Projection的设置）

使用GDAL给遥感影像设置地理参考和投影，分别设置两个参数：GeoTransform和Projection。对于没有附带地理信息的影像，需要从其它途径获取参数，如xml等元数据。GDAL 中，GeoTransform是一个六个元素的元组。六个参数分别为：（左上角x坐标，水平分辨率，旋转参数，左上角y坐标，旋转参数，-垂直分辨率）应该注意，最后一个参数是负值，旋转参数一般是0。而对于没有投影的图像（如WGS 84只有地理坐标系），它的分辨率不是以米位单位，而是经纬度，所以不能简单地指定分辨率这两个

2022-05-18 15:58:14 5992 1

原创解决GDAL写投影问题——proj.db报错

处理高分影像时，因为原数据没有存储投影信息，所以需要自己定义。使用GDAL（osr）对遥感影像定义投影。使用以下代码from osgeo import gdal, osr’‘’其他读取影像数据的程序‘’‘# 获取地理坐标系统信息，用于选取需要的地理坐标系统srs = osr.SpatialReference()srs.ImportFromEPSG(4326) # 定义输出的坐标系为"WGS 84"out_IRS.SetProjection(srs.ExportToWkt()) # 给新建

2022-05-18 15:38:24 7022 10

原创使用imgaug进行语义分割数据集增强

使用imgaug批量进行图像增强，实现图片和标签mask同步转换，基于三波段png格式遥感影像数据集

2022-04-09 12:04:24 3903 12

原创远程服务器搭建MySQL数据库

背景在客户提供的虚拟机搭建数据库，完成建表和数据入库。系统：目标机CentOS 7，本地机Win10软件：Xshell 7、navicat操作过程由于是第一次远程在客户机器上建库，且目标机没有安装图形化界面，也不能连接互联网，所以首先需要客户开通维护端口22和数据库端口3306。经过沟通开通了可以连接客户内网的VPN，可以通过VPN访问22端口；3306端口外网可以访问到，使用navicat可以实现可视化建表和数据入库。通过Xshell远程连接首先连接VPN访问内网，然后使用Xshell

2022-03-28 11:59:02 2559

原创 json上传数据库的一些mysql命令

背景：Win10系统，通过Xshell连接远程数据库，上传数据。建表和简单数据录入直接通过navicat完成。问题：使用Navicat粘贴比较大的json数据到表字段中时，出现长度受限的情况，无法将完整的json数据放入数据库表字段。查了mysql给的指南，JSON的长度应该是和LONGTEXT一样的，所以理论上能存下完整的数据，但navicat无法完整的粘贴进来，也不能完整显示已经存好的数据。所以只能通过mysql命令行来插入和查看了。JSON Storage RequirementsIn

2022-03-21 16:19:00 1129

原创配置多GPU训练的笔记-xview2挑战项目复现1

最近尝试复现一个xview2挑战（xiew2）的冠军项目，这个挑战的目标是基于遥感影像对自然灾害前后建筑物的损坏程度进行定位和评分，官方提供了非常庞大的数据集，原项目代码在Github能找到（xview2_1st_place_solution）。遇到的坑主要有以下几个：1.APEX安装项目用到了NVIDIA的APEX库，需要先从Github把这个项目克隆下来（https://github.com/NVIDIA/apex），再进到apex项目里安装。官方给的Linux系统安装方式是：git clone

2022-03-04 11:22:27 1870 8

原创解决在Word安装Mendeley失败，提示“错误429：activex部件不能创建对象“问题

找到Mendeley的安装路径，注册MendeleyWordPlugin.exe.如图所示操作即可。

2021-12-08 23:25:42 1671 3

转载转载：遥感图像空间索引和提取像元值

遥感影像有两组坐标系。一组是图片的行列坐标，另一组是空间坐标。遥感影像中的每一个像元都对应地球上的一个地点。使用rasterio可以方便的对两组坐标系进行操作。假设我们需要获取某一地理坐标的NDVI数值，如经纬度(-119.770163586, 36.741997032)，但是我们数据是UTM坐标系的，所以我们首先要将经纬度坐标投影为UTM，我们可以使用pyproj完成投影。作者：雪涵链接：https://zhuanlan.zhihu.com/p/412865295来源：知乎著作权归作者所有。商业转

2021-11-24 14:51:28 941

原创上传本地文件到远程仓库-Gitlab

首次上传本地文件到远程仓库，由于远程仓库不是自己建的，而将本地内容push上去必须让本地与远程版本匹配，在此基础上进行修改。而我是自己git init新建的本地仓库，所以肯定是不行的。所以对于非自己建立的远程仓库，首先要把它克隆到本地。git clone 远程仓库的SSH然后整个远程仓库的文件被克隆到本地了，再把自己的文件放到这个本地仓库中想要的位置，然后更新远程仓库的版本即可。git add . 这里“.”表示上传所有文件，会检测所有的更改。也可以把“.”换成文件名。然后将add的文件co

2021-11-10 14:31:56 617

原创 pandas表合并——按指定列合并两个表里同样名称的行

在处理观测站气象数据时，遇到了需要将两个表的数据整理到一起的问题，两个表的数据结构如下：现在需要将两个表按第一列‘Sta_ID’拼接，将第二个表中‘Sta_ID’值与第一个表相同的行拼接到同一行，效果如下：需要的库基本有pandas，xlrd，xlwt这几个，首先使用pandas.read_excel(）将数据分别读入pandas。A = pd.read_excel('E:\\MODIS\\MCD19A2\\AOD_Max.xls', sheet_name='AOD_Max') # 读取 exc

2021-11-03 14:31:19 13572

原创 python 导入gdal失败：DLL load failed 找不到指定的模块。

最近在conda环境中安装了pyqt和labelme等一些库之后，打算运行一个geotif转jpg的程序，但一开始运行就报错，显示imp.py中import importlib.machinery这句出了问题，最后一行return _load(spec) ImportError: DLL load failed: 找不到指定的模块。虽然import不提示语法错误，但运行就报错，推测是gdal库导入出了问题，由于之前安装的gdal是可用的，现在不能用了，捣鼓了半天，最后直接重装解决问题，与版本没关系。c

2021-10-08 14:19:40 8659 7

原创使用Python读写netCDF4（.nc）格式的VIIRS数据

**使用Python读写netCDF4（.nc）格式的VIIRS数据**使用的主要工具：python 3.7gdal 3.2.1netcdf4 1.5.5.1首先介绍netCDF4格式的数据存储方式，它是一种多维数据格式，类似数组组织形式的数据集合。在python中可以直接用data = netCDF4.Dataset(data_path)读取数据集，然后用data.dimensions查看数据集的维度和详细信息。本次用的地表反射率数据，文件名为SurfRefl_v1r1_j01_s20

2021-06-24 15:43:12 1819 1

原创 Python脚本批量读取哨兵2号（Sentinel2）影像并另存为Geotiff格式

由于Sentinel影像的存储格式特殊，通常需要特殊的软件（如SNAP）来处理，有局限性。于是使用GDAL对其进行读取，并重新写为.tif格式，方便进一步处理。影像数据用的是大气校正后的L2A级数据，当然L1C数据也是也行的。打开SAFE后缀的文件夹，可以看到里面的内容，jp2格式的图片数据是存储在GRANULE文件夹中，这里不需要知道它的存储规则，因为GDAL可以直接读取最下面的.xml文件（注意GDAL的版本，2.4以下的版本可能不支持）。文件的打开方式与其他格式的遥感影像相同，用gdal.Ope

2021-03-10 16:38:32 7245 36

原创使用IDM实现批量下载Sentinel-2影像（Python）

使用IDM实现批量下载Sentinel-2影像（Python）需要的工具：IDM下载器：http://www.internetdownloadmanager.com/sentinelsat包：命令行pip安装即可 pip install sentinelsat主要流程是，从数据网站选择数据并产生链接文件，使用python解析资源地址，并调用IDM创建下载队列。首先要获取数据源地址，需要在欧航局官网手动选择自己需要的数据，然后加入购物车，进入购物车后点击右下角下载按钮，可以下载到一个.meta4文件

2021-03-09 17:06:44 1558 3

原创遥感图像匹配(高分4匹配Landsat8)-Python+Opencv/SIFT+GDAL

高分系列遥感影像的一个问题是难以进行精准的几何校正，通常是在专业软件（如ENVI）中处理，需要借助标准影像手动选择地面控制点来进行校正，操作比较复杂。本文的思路是借助Landsat8影像来进行关键点匹配，这样就能自动得到控制点。难点有：①两种不同卫星的影像分辨率和细节方面差异较大，难以进行准确的匹配。为解决此问题，将分辨率较高的Landsat影像进行了均值滤波。②遥感影像的存储格式不能直接用opencv处理，需要转换成8位整型灰度图，这里只用了一个全色波段。# -*- coding: utf-8 -*

2021-03-01 14:03:03 3280 1

KilllerQueen的博客

原创 python 3维矩阵转为2维

原创遥感影像叠加，任意波段数tif影像组合

原创 python 科学计数法的小数转换

原创 python中datetime不同时区转换

原创设置python脚本开机自启--anaconda虚拟环境、windows系统

原创 Python调用webapi，使用Token认证

原创 Sqlsugar中AndIF()条件查询C#

原创实现不同尺寸和不同波段数的图像叠加(layerstaking)——按地理位置对准GDAL,Numpy