- 博客(74)
- 资源 (4)
- 收藏
- 关注
RSS订阅原创 【g2o】 解决 g2o::BlockSolverTraits＜-1, -1＞]: Assertion `_sizePoses ＞ 0 && “allocating with wrong size“‘
Thirdparty/g2o/g2o/core/block_solver_impl.h:70: void g2o::BlockSolver::resize(int*, int, int*, int, int) [with Traits = g2o::BlockSolverTraits]: Assertion `_sizePoses > 0 && “allocating with wrong size”’ failed参考 issue 中的解决方式我在 中修改 进行了修改。...
2022-06-23 22:36:19
518
原创 【git】拉取仓库、本地修改、上传修改
把工程clone到本地查看工程状态查看工程的日志把工程最新的情况pull下来新建想要修改的分支本地修改把修改的东西add为修改的东西注释把修改好的东西push到分支在网页新建merge请求
2022-06-10 17:58:21
237
原创 【视觉SLAM十四讲】学习笔记1: 相机模型和畸变模型
本文是**《视觉SLAM十四讲》（第2版） 第5讲：相机与图像 ** 的阅读笔记，具体内容可查阅书籍和开源代码。相机模型在讲述模型前，首先需要明确三个坐标系：世界坐标系、相机坐标系和像素坐标系。针孔模型空间点P，世界坐标系坐标为[Xw,Yw,Zw]，相机坐标系坐标为[X,Y,Z]Pc = RPw+t = TPw点p的归一化坐标为[X/Z,Y/Z,1]Pnorm=1/Z Pc点P在成像平面上的投影点P’,相机坐标系坐标为[X’,Y’,Z’]X’ = f X/ZY’ = f Y/Z4
2021-12-28 21:14:27
1937
原创 【VMware】虚拟机和主机互传文件
虚拟机和主机互传文件安装完VMware tools后确实能传文件就是容易蓝屏所以，网速好的时候，网盘是个不错的选择
2021-05-24 22:20:10
552
原创 【python】记一次打包经历
1. cython编译(optional)参考链接2. pyinstaller打包参考链接3.遇到的问题'utf-8' codec can't decode byte 0xce in position 138: invalid continuation byte These are the bytes解决方法
2021-05-23 12:05:02
193
原创 【sarscape】使用envi idl 实现哨兵一号批量ComplexDataMultilooking
SARscape是强大的雷达影像处理软件，其中影像导入、影像裁剪等基本常用功能都可以进行批量处理非常方便。但是其中的同轨道影像镶嵌(Slant Range Mosaicing)功能，一次操作 只能完成一个轨道同一时间段下的影像拼接，当需要批量处理时，多次的点选以及时间的把控就很不方便。因此，参考SARscape IDL Scripting中的example以及IDL的Help文档，实现哨兵一号的两景影像批量镶嵌。1.数据准备待多视的影像2.实现过程准备.sml文件在envi软件中，sar
2021-03-24 21:49:34
2056
原创 【sarscape】使用envi idl 实现哨兵一号批量SLC配准到同一主影像
SARscape中的SLC配准(SLC coregistration)功能【注意，并不是coregistration功能】，界面如????。功能界面的样子应该是可以批量配准的，但是一直没有成功，条条大路通罗马，所以决定用idl实现批量配准到同一主影像的功能。1.数据准备本次输入数据，是镶嵌后的VV极化影像存放主影像的文件夹(只有1景)存放待配准影像的文件夹(1景或多景)2.实现过程准备.sml文件在envi软件中，sarscape/preference中设置好参数，然后保存。打
2021-03-24 21:26:42
2508
原创 【sarscape】envi idl 记录
SARscape IDL使用说明查询从软件中进入其中详细介绍了如何查找模块、模块参数及模块调用等方法例子：SARscape IDL 查看模块参数类型保存参数设置SML文件(后续用于IDL中的参数设置)上述保存的sml文件是用户设置的一些参数，可以在处理过程中调用????，使用方法可以参考博文【sarscape】使用envi idl 实现哨兵一号批量镶嵌...
2021-03-24 20:52:09
398
原创 【sarscape】generic error:20000解决办法
1.问题背景：Coherence工具，报错20000本文中提到的数据，是哨兵IW模式的SLC(单视复数影像)。使用/SARscape/Interferometry/Interferometric Tools/Complex data multilooking对复数影像(slc)进行8：1多视后，得到的multilook的结果输入/SARscape/Basic/Feature Extraction/Coherence中，报错generic error:20000.查看报错文档????，文档中提到用于干
2021-03-24 10:08:00
6282
2
原创 【python】scipy包中的BFGS算法
因为需要使用无约束非线性优化问题，所以想用拟牛顿法中的BFGS算法，在此记录一下资料收集的过程。1. 参考链接对于优化算法的介绍，可以参考这篇博文，博主写的比较清楚。常见的几种最优化方法（梯度下降法、牛顿法、拟牛顿法、共轭梯度法等）对于scipy中一些优化算法的介绍，可以参考官方文档。scipy.optimize2.实践使用不同的方法，求解最小值问题，可以使用????scipy.optimize.minimize1. 函数介绍满足初级使用，需要了解以下????几个关键参数，其余细节的参数
2021-03-13 10:51:19
3624
原创 【python】使用gdal存储envi标准格式的complex影像
gdal官方文档自定义函数数据类型：datatype = gdal.GDT_CFloat64文件类型：driver = gdal.GetDriverByName("ENVI")#保存tif文件函数def writeTiffEnvi(data,geotrans,proj,path): if 'int8' in data.dtype.name: datatype = gdal.GDT_Byte elif 'float16' in data.dtype.name:
2021-02-22 13:49:07
1309
原创 【sarscape】使用envi idl 实现哨兵一号批量镶嵌
SARscape是强大的雷达影像处理软件，其中影像导入、影像裁剪等基本常用功能都可以进行批量处理非常方便。但是其中的同轨道影像镶嵌(Slant Range Mosaicing)功能，一次操作 只能完成一个轨道同一时间段下的影像拼接，当需要批量处理时，多次的点选以及时间的把控就很不方便。因此，参考SARscape IDL Scripting中的example以及IDL的Help文档，实现哨兵一号的批量镶嵌。2.实现过程准备.sml文件在sarscape安装目录中找到，SARscape_defaul
2021-01-14 15:20:55
3413
6
原创 【sarscape】报错READF: End of file encountered. Uint: 108 的解决办法
1.问题描述在新机上安装envi5.3和sarscape5.2.1。在安装sarscape的过程中提示，opencl已经有更高版本存在，选择无视。在不做任何操作的情况下，运行envi，报错如下：打开报错的 .txt文件，发现其内容为空。原因依旧是，sarscape没有检测到可以运行的platform。2.解决方法参考 前文【sarscape】无法使用CUDA加速计算的解决方法即可解决。...
2020-12-18 14:27:07
6030
2
原创 【context capture】合并两个block时，高度出现错位的一种解决办法记录
写在前边下文提到的解决办法，仅为一家之言，经测试后暂时解决了目前的问题，从目视效果上看没有明显错误，但是无法保证其精确性。>.<目前并不知道 产生这种错误的原因是什么，可能是我搜索的方向不对，希望有更好解决方法的朋友帮忙纠正和留言交流。>.<问题背景野外拍摄无人机影像时，分块，航高设置相同，但是起飞点不同。因此在空三时，需要分block进行，否则空三无法通过。在merge block时，9个区块中，有2个与其它几个出现了高度错误，如下图。经过查看，发现2个出错的bl.
2020-12-06 17:07:57
1337
1
原创 【sarscape】同一区域，不同轨道数据，时序InSAR中，轨道精炼GCP的选取
1.参考链接SARscape干涉处理中轨道精炼控制点的使用2.用SARscape做时序InSAR提取地表形变时，在轨道精炼的一步，必须输入控制点文件（refinement GCP file），用于修正干涉相位和解缠后的相位，主要包括：– 去除轨道误差（系统的几何误差)– 去除残余地形相位– 减少解缠相位的错误当对同一区域，使用不同观测几何的数据对（升轨数据对和降轨数据对）进行处理时，在轨道精炼时，最好使用同一组地面控制点，这样可以减少两组数据间的相位误差，更方便两组结果的比对。主要实现方法为
2020-11-23 17:57:51
2917
6
原创 【python】sarscape输出的哨兵一号雷达影像干涉图展示
1.参考链接envi sarscape 帮助文档np.arctan() math.tan() math.atan()【python】读取envi输出的复数二进制影像文件2.实现过程读取复数二进制影像文件参考 参考链接3根据sarscape的文档介绍，干涉图相位的计算公式如下图3.使用np.arctan()来计算干涉图复数影像对应的相位，得到相位的弧度值import gdalimport numpy as npdef img2array(path): dataset = gd
2020-11-21 10:56:12
1779
1
原创 【python】读取envi输出的复数二进制影像文件
1.参考链接python复数属性和方法操作实例2.复数二进制影像文件此处读取envi sarscape时序干涉处理过程中生成的复数二进制文件。读取二进制影像文件，使用gdal.Open(path)即可读取，path不需要加入后缀名读取得到的数组 的 dtype为complex64提取复数的实部和虚部进行后续计算import gdaldef img2array(path): dataset = gdal.Open(path) width = dataset.RasterXS
2020-11-21 10:35:40
2115
原创 【sentinel】哨兵二号的归一化积雪指数
1.参考链接ENVI5.4改进的归一化雪盖指数（NDSI）哨兵2号（Sentinel-2）介绍、下载和预处理、批处理积雪指数NDSI2.归一化积雪指数公式：– NDSI = (GREEN - SWIR) / (GREEN + SWIR)– Green代表的是0.5µm到0.6µm范围的绿波段；SWIR1代表的是1.55µm到1.75µm范围的短波红外波段。用途：区分雪和云原理：– 雪和云在可见光波段和远红外波段反射率相似– 雪和云在1.6µm和2.1µm附近的近红外波段反射率差异大，
2020-11-17 18:55:10
4542
4
原创 【精通图像处理经典算法(MATLAB版)】学习笔记【1】
1.图像处理基础图像及图像数字化模拟图像：包括光学图像、照相图像、电视图像等，是连续的。数字图像：将连续的模拟图像经过离散化处理后得到的计算机能够识别的图像，是离散的二维采样数组，由矩阵表示。采样及量化采样：图像信号是二维空间的信号，是一个以平面上的点作为独立变量的函数。二维图像信号的采样遵循采样定理，常用办法是在二维空域内进行均匀采样。量化：模拟图像经过采样后在时间和空间上离散化为像素，但是采样所得像素值仍为连续量。将采样所得的各像素灰度值从模拟量转换为离散量称为图像灰度量化。连
2020-11-15 16:54:54
504
原创 【python】将一个数组中的值替换为另一个数组中的相应值
1.参考链接2.布尔蒙版将数组b中 对应数组a大于50的位置 替换为数组c中的值b[a>50]=c[a>50]
2020-11-10 21:45:10
3611
原创 【python】3. 调用matlab的m文件的一种方法
1.参考链接Python调用matlab脚本终极解决方案MATLAB R2015b使用MATLAB Compiler Runtime生成Python库的方法2.编译.m文件为python库matlab命令行中deploytool→library compiler→python package在exported functions 处选择需要导出的.m文件点击package导出3.python安装activate python36# 包所在路径cd autoFmaskBatch\for
2020-10-15 17:56:35
958
原创 【matlab】3.解决library complier没有编译器的问题
1.问题在matlab命令行输入deploytool，选择 library compiler，发现type处没有可用的编译器选项2.解决重装matlab勾选编译器安装完成后，再次打开，就有了
2020-10-15 09:40:02
1596
原创 【PyQt5】2. 设计界面
1.参考链接2.文件转换拖动左侧小部件，设计界面，比如????保存为.ui文件转换为.py文件切换到.ui所在路径，运行pyuic5 -o gui.py gui.ui，同目录下生成.py文件3. .py文件解析
2020-10-14 12:06:04
229
原创 【matlab】1.使用mex编译.c文件
mexfunction函数为C语言与MATLAB语言的接口函数，调用实例在.c文件中。在MATLAB中调用mex指令编译相关文件，将C语言编译为MEX文件。编译完成后，编写MATLAB函数，调用MEX文件。以MEX文件的形式调用编译完成的C语言函数。在matlab命令行中启动编译，报错提示安装编译器...
2020-10-10 22:40:33
2407
原创 【win10】64位windows安装gcc
1.参考链接Mingw-w64在win10下的安装使用2.安装下载后运行将mingw64\bin所在路径添加到环境变量3.测试gcc -v
2020-10-08 22:39:19
2528
原创 【python】运行报错：Windows MemoryError: Unable to allocate 6.38 GiB for an array with shape
1.参考链接成功解决Windows MemoryError: Unable to allocate 6.38 GiB for an array with shape (38
2020-10-08 22:24:41
16421
原创 【ArcGIS】1.在指定python版本下安装包
1.参考链接python小技巧之：同时安装py2和py3,如何在指定的版本下pip安装包2.描述使用arcgis的IDLE时候，需要tqdm包，但是电脑中有python3和python2两个版本。要在arcgis的python27中安装tqdm，使用参考链接中的方法3。找到python27的python.exe目录运行 python.exe -m pip install tqdm...
2020-10-06 20:18:18
393
原创 【sentinel】3.snap借助gpt命令行实现sentinel-2批量ndvi
1.参考链接Sentinel-2 Toolbox Featuressnap中的s2toolbox提供多种指数的计算。SNAPsnap中可以进行掩膜以及波段计算。在SNAP中利用Sen2Cor进行哨兵二号数据预处理（大气校正，辐射校正）2.安装snap下载链接安装后打开3.安装sen2cor插件tools → plugins...
2020-10-06 14:51:53
619
原创 【sentinel】2.python批量下载哨兵数据报错：HTTP Error 502
1.参考链接【已解决】python执行出错：HTTPError: HTTP Error 502: Bad Gateway2.报错使用python和idm批量下载哨兵offline数据时报错：HTTPError: HTTP Error 502: Bad Gateway3.解决原因是频繁向网站发送请求，过一段时间再次操作即正常。...
2020-10-06 10:21:05
791
原创 【Fmask】2. 启动matlab软件时, 出现函数或变量“ctfroot“无法识别
1.问题掩膜提取软件 fmask-gui版本 由matlab启动时候，报错。2.参考链接解决方法2.实际操作删除该文件后，再次启动软件后即可正常运行。C:\Users\13784\AppData\Local\Temp\13784\mcrCache9.7
2020-10-05 20:44:55
3196
原创 【sentinel】2. sen2cor对L1C进行大气校正
1.参考链接Sen2Cor v2.8官方网站Sentinel-2批量大气校正_基于Sen2Cor2.插件安装windows安装包下载链接解压将包含L2A_Process.bat的目录，加入到系统路径命令行输入L2A_Process.bat --help，检查是否成功3.批处理for /D %s in (I:\S2\0.202003-202006\S2A_MSIL1C*) do L2A_process --resolution=10 %s ...
2020-10-04 16:53:26
1191
2
原创 【google earth engine】11.配置gee 本地环境
为了调用matlab的函数，从云端转本地。环境配置稍有不同，google colab的配置过程见此链接。1.打开anaconda prompt，进入虚拟环境activate python36安装gee的包conda install -c conda-forge earthengine-api检查api是否成功，进入 pythonimport ee授权gee，在终端输入????，跳出网页，复制验证码earthengine authenticate...
2020-10-02 16:06:05
895
SARscape5.2.1
2020-11-15
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人