扎不下村村长的博客

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

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>

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

项目需要，在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++实现，在控制文件（文本格式中）按行取数据控制文件格式，类似于这样的，一行一个或者多个数据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上，可以更好的看出效果比如这张速率图（结果不好，见谅）首先在PS中，新建一个像素透明背景的画布，大概大小在8-10个像素，之后使用椭圆工具在画布中心绘制一个6-8个像素大小的黑色圆点击编辑选项，下拉菜单内点击“定义图案”，编辑一个自己的名字从PS打开需要制作的图片，在该图片下方新建两个图层，下面图层填充上背景颜色，上面图层为透明，并点击该图层，按Shift+F5打开填充窗口选择刚才保存

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

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

为了实现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.

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

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

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

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

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

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 ...

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 ...

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);

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_...

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

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

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

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

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

EDM工具箱Imf=emd（X）； X为需要分析的信号的时间序列Emd_visu(X,t,imf); 可视化操作， X为需要分析的信号的时间序列，t为时间的序列，imf为上一步生成的图一是结果 residue是趋势项；图二、图三两个是重构原信号的方式，就是把细节一点一点加上去和把趋势一点一点加上去功率谱、频率图选择需要分析的imf序列[Pxx,F] = pwelch(imf(2...

emd - computes EMD and bivariate/complex EMD with various options计算EMD和双变量/复杂EMD，并提供各种选项emd_online - computes on-line EMD variation. Note that it does not truly apply on-line: the functio...