扎不下村村长的博客

原创 利用ISCE2处理时序COSMO-SkyMed数据

COSMO-SkyMed卫星是意大利航天局和意大利国防部共同研发的COSMO-SkyMed高分辨率雷达卫星星座的第二颗卫星。COSMO-SkyMed卫星是X波段合成孔径雷达（SAR）卫星，波长约3.11cm。

原创 Linux离线安装Armadillo以及各种依赖软件

Armadillo是一个线性代数C++库，其中不光包含了各种线性代数的基本运算，还包含了一些常用的线性运算程序（如PCA等）以及稀疏矩阵方法等。

原创 Linux安装MiaplPy

MiaplPy是迈阿密大学（University of Miami）InSARLab团队开发的一款基于相位连接（Phase Linking）的InSAR处理软件包。适用于全分辨率（即多视比为1:1）的非线性（Non linear）相位反演，采用的输入数据是ISCE2生成的共配准（Co-registered）之后的SLC数据。

原创 c++程序在Windows和Liunx中运行不同代码

#if defined (WIN32)||defined (_WIN32)|| defined(WIN64)||defined(_WIN64)&& !defined(__CYGWIN__)//Windows代码cout<<"this is Windows!"<<endl;#else//Linux代码cout<<"this is Linux!"<<endl;#endif

原创 c++获取当前程序运行的目录

string running_path = _pgmptr;

原创 c++获取文件路径

在c++中用string获取输入或者输出目录的不包含文件名的路径out_path = out_path.substr(0, out_path.find_last_of("\\") + 1);

原创 c++在Windows中避免使用/，导致找不到文件

在Windows中经常会遇到输入输出时，存在/导致系统无法找到文件，所以在打开或者输出的时候，需要把目录中所有的/都换成\string& replace_all_distinct(string& str, const string& old_value, const string& new_value){ string::size_type pos = 0; while ((pos = str.find(old_value, pos)) != string::npos

原创 c++中使用GDAL获取图像上某个点的值

获取一系列点的值点的信息存在txt中，使用read_point 类读取txt分为三列 点名 经度 纬度n1 3.1 49.991n2 3.2 49.992n3 3.3 49.993n4 3.4 49.994n5 3.6 49.996n2 3.5 49.992n2 3.33 49.9923#include "read_point.h"#include<fstream>#include<iostream>#include<sstream>void

原创 c++中使用GDAL将两张影像相减

void image_process::do_diff(int master_date, int slave_date, string img_path, string out_path){ img_path.append("\\"); out_path.append("\\"); string filepath1 = img_path + to_string(master_date) + ".tif"; string filepath2 = img_path + to_string(slave

原创 使用dos或cmd删除文件

在Windows中自带的cmd中使用del /F /S /Q .\\file需要在c++中调用的，可以使用string mvfile1 = "del /F /S /Q .\\" + filename;cout << mvfile1 << endl;std::system(mvfile1.c_str());

原创 c++借助已有软件，实现压缩文件的解压

使用c++打包解压软件，实现压缩文件的解压常用WinRAR Bandizip等软件 注意版权这类软件除了平时使用的图形界面形式，通常可以使用命令行CMD调用，WinRAR自带说明书、Bandizip的说明Bandizip推荐使用Bandizip，使用方便、通用性强、还可以避免不必要的问题string winrar_path = "WinRAR.exe";unzip = "\""+winrar_path + "\"" + " x C:\\gacos_temp\\" + txt_file.get_n

原创 c++读取XML文件

在c++中使用tinyXml类库，来读取或写XML文件在头文件中引入#include"tinyxml2.h"void read_xml_file::read_xml_inf(string file_path){ XMLDocument doc; doc.LoadFile(file_path.c_str()); XMLElement* root=doc.RootElement(); XMLElement* element = root->FirstChildEle

原创 c++ 简单方法来判断文件是否存在

inline bool exists_file (const std::string& name) { struct stat buffer; return (stat ((char *)name.c_str(), &buffer) == 0); }返回1表示存在，0表示不存在需要引入头文件#include <sys/stat.h>

原创 c++分离文件名和路径

string outname_base = filename; //这里使用string或者char数组都可以 //string用的方便，但是很多时候需要转换 char *p = (char*)outname_base.c_str(); string outname = basename(p); //最后斜杠之后的文件名 string outdir = dirname(p); //文件名之前的所有路径...

原创 c++ 实现简单的ftp上传和下载

项目需要，在windows下使用c++利用ftp将文件上传/下载到远程服务器不需要写ftp客户端，能够借助windows中ftp.exe即可原理：在cmd中执行，需要使用bat批处理，命名为b.bat在批处理中写入命令ftp.exe -n -s:a.txt而a.txt中就是需要执行的命令这里只写了几个基本的，登录，cd，binary等输入的参数host 主机地址，默认端口21username 登录ftp的用户名pass 登录ftp的密码file 上传或下载的文件

原创 c++ 从txt文本中按行取数据

项目需求，用c++实现，在控制文件（文本格式中）按行取数据控制文件格式，类似于这样的，一行一个或者多个数据Is/20190704-20190716/result/20190704-20190716.aIs/20190716-20190728/result/20190716-20190728.aIs/20190728-20190815/result/20190728-20190815.aIs/20190809-20190815/result/20190809-20190815.aIs/201908

原创 影像结果制作点状图，并绘制到Google Earth

为了更好地展示影像计算的结果，可以制作点状图，将整幅影像的结果打散为一个一个的点，再展到Google Earth上，可以更好的看出效果比如这张速率图（结果不好，见谅）首先在PS中，新建一个像素透明背景的画布，大概大小在8-10个像素，之后使用椭圆工具在画布中心绘制一个6-8个像素大小的黑色圆点击编辑选项，下拉菜单内点击“定义图案”，编辑一个自己的名字从PS打开需要制作的图片，在该图片下方新建两个图层，下面图层填充上背景颜色，上面图层为透明，并点击该图层，按Shift+F5打开填充窗口选择刚才保存

原创 如何理解同震、震后、震间的含义

这个三个词都是用来表示地震发生前后地面的形变的时间段同震：表示震前的一段时间到震后的一段时间的形变震后：表示地震发生后的形变震间：表示两次地震（较大地震）发生之间的地面变化...

原创 将InSAR处理结果绘制为带有等高线的网格图

为了将InSAR处理的结果更加形象的显示，可以利用Matlab中自带的一些函数和方法对InSAR处理的结果图像进行绘制，绘制成为带有等高线的网格图a=imread('E:\vel_los_utm.tiff');%%读取tiff影像数据，这里只能读取到值，其他的坐标等信息需要后面人为添加b=a(6001:6500,6001:6500);%%选择绘制的区域，这里需要注意使用Matlab进行绘图时，数据量不能太大，我的电脑只能绘制1000*1000的格网，所以选择小一些的范围b=smoothdata(b,'

原创 Matlab保存动态图像为gif

为了实现matlab中能够将动态的三维图像保存执行以下命令pic_num = 1;for i=-1:- 1:-89 %%% for循环 view(i,18); %%%% 移动视角 pause(0.09); %%%% 暂停时间 F=getframe(gcf); I=frame2im(F); [I,map]=rgb2ind(I,256); if pic_num == 1 imwrite(I,map,'test1.

原创 Matlab移动色带位置

caxis %%显示最大最小值通过最大最小值来移动colorbar位置caxis([-10,10]); %%设置色带范围

原创 Matlab将底色改为白色

该开始时在MATLAB命令行输入并执行set(0,‘defaultfigurecolor’,‘w’)之后无需设置，均为白色或者在每次绘图之前执行figure(‘color’,[1 1 1]);

原创 第一个bash脚本，实现批量画图

为了实现在Linux服务器上能够批量的处理数据写了一个bash脚本实现读取文件目录里的经纬度坐标和名词，然后利用软件功能进行绘图，生成图片首先构建文件mk_local，第一列经度，第二列纬度，第三列名词112.4463842385 35.3966815865 固隆煤业112.4609471513 35.5214877910 晋圣松峪煤业112.5722615384 35.5067995420 润东煤业112.5722615384 35.5067995420 凤红煤业112.572261538

原创 Origin导入数据时设置对应的时间

时间格式按照yyyymmdd排列的数据格式，导入到Origin中显示是数值&文本在属性中选择日期，在显示中选择自定义，在自定义显示中选择yyyymmdd

原创 vscode中的${workspaceFolder}等变量

${workspaceFolder} :表示当前workspace文件夹路径，也即/home/Coding/Test${workspaceRootFolderName}:表示workspace的文件夹名，也即Test${file}:文件自身的绝对路径，也即/home/Coding/Test/.vscode/tasks.json${relativeFile}:文件在workspace中的路径，...

原创 软件推荐--Files2Folder

将文件放入新建文件夹内选中要打包的文件，右键选择Files2Folder，输入文件夹名，即可打包到文件夹里官网：http://www.dcmembers.com/skwire/download/files-2-folder/经测试，windows10 不可用...

原创 你的“开始”菜单出现了问题，我们将尝试在下一次登录时修复它？

今天遇到了一个问题，当我点击开始菜单，居然出现了关键错误你的“开始”菜单出现了问题，我们将尝试在下一次登录时修复它只能选择注销登录但是注销登录进去还是同样的问题网上能够找到的方法主要是把windows服务中的，User Manager修改为自动，但是我的本来就是自动所以问题不出现在User Manager服务上这时我还可以打开CMD尝试用chkdsk来修复C盘chkdsk c: ...

转载 Git使用教程

网上大部分的git教程只是从简单的第一次操作为例子介绍的，而实际使用时会出现很多问题。熟练掌握Git命令，可以提高开发的效率。安装GitWindows直接在官网上去下载。下载完成后，随便在某个文件下右键如果有Git Bash Here就安装成功。安装后，还要在命令行输入：$git config --global user.name"你的名字"$git config --global us...

原创 用于处理CRS发布的GRACE RL06版本的数据的Matlab程序

function [cs,cs_sigma,int_year,int_month,time] = gmt_readgsm_RL06(dir_in,file_name,lmax)% Read the gravity filed files%% INPUT:% dir_in full path% file_name filename % lmax ...

原创 GRACE RL06版本的数据预处理

function GRACE_RL06_preprocessing(controlfile_path)% Read the Control Filefid=fopen(controlfile_path,'r');num_file = fscanf(fid,'%d',1);radius_filter = fscanf(fid,'%d',1);destrip_method ...

原创 画出GLDAS数据的matlab程序

function plotgldas(gldas,time,num_file)gldas(find(isnan(gldas)==1)) = 0;for k=1:num_fileave_gldas(k,1)=mean(mean(gldas(k,:,:)));endtime=time-min(time);plot(time,ave_gldas);

原创 处理GRACE数据而进行的水文模型GLDAS的处理matlab程序代码

function [gldas,time,num_file]=readgldas(controlfile_path)% Read the Control Filefid=fopen(controlfile_path,'r');num_file= fscanf(fid,'%d',1);lat= fscanf(fid,'%d',1);lon= fscanf(fid,'%d',1);dir_...

原创 GPS掩星技术

原创 sentinel1精密轨道下载

https://qc.sentinel1.eo.esa.int/POD Precise Orbit Ephemerides [AUX_POEORB]

原创 InSAR数据的获取模式

原创 用Global Mapper拼接多个DEM

导入下载的DEM 数据选择文件导出 Elevation Grid Format其他选项默认

原创 GRAMAT控制文件

GSM文件数高斯平滑的距离（单位km）从选项中指定去条带方法（NONE/SWENSON/CHAMBERS2007/CHAMBERS2012/CHENP3M4/CHENP4M6/DUAN）指定GIA效果是否将被移除（GIA_notRemoved/GIA_Removed_Geru）输入格式（ ICGEM/GRACE）输出格式、最大阶数/次数和输出文件名（分辨率），部分选项：SH_MAT 60...

原创 Matlab-mean函数

M = mean(A)返回沿数组中不同维的元素的平均值。如果A是一个向量，mean(A)返回A中元素的平均值。如果A是一个矩阵，mean(A)将其中的各列视为向量，把矩阵中的每列看成一个向量，返回一个包含每一列所有元素的平均值的行向量。如果A是一个多元数组，mean(A)将数组中第一个非单一维的值看成一个向量，返回每个向量的平均值。...

原创 Matlab-cell数组

cell数组一般被叫做元胞数组，它的每个单元可以储存不同的数据类型，可以是数值，字符或矩阵或元胞数组等，类似于学过的c语言里的结构体1.cell数组的创建可以直接通过{}类似于矩阵的直接赋值：a={‘winter’,123,‘coming’,‘哈哈’};也可以通过cell函数预分配内存，再赋值：a=cell(1,4);a={‘winter’,123,‘coming’,‘哈哈’};2.cel...

原创 GRACE处理流程preprocessing_core

打开控制文件读取文件数读取高斯平滑半径读取去条带方法检查所选方法(NONE,SWENSON,CHAMBERS2007,CHAMBERS2012,CHENP3M6,CHENP4M6,DUAN)读取GIA选项检查GIA选项（GIA_notRemoved，GIA_Removed_Geru）读取文件类型检查文件类型（GRACE,ICGEM）读取输出文件类型 SH_MAT + 最...