- 博客(19)
- 收藏
- 关注
RSS订阅原创 对已转化成工程坐标后的.COR文件进行进一步统计、剔除错误值、采样以及重采样，并输出用于探地雷达数据水平归一化、高程校正的文件
（4）每10个RTK定位进行一次采样，进行校正；每5米进行一次采样，进行校正；输出实际轨迹方向、布置测线方向的归一化文件，共4个。Result1，布置测线方向上，每10秒采一次样；Result2，将Result1重采样为每5米一个样，再放大10倍；Result3，在实际轨迹上，每10秒采一个样；Result4，将Result3重采样为每5米一个样，再放大10倍；（2）采用机器分类、人眼识别的方法剔除RTK浮动解；（3）输出高程校正文件correct.MAR；（1）统计定位信息；
2025-07-29 19:11:03
179
原创 经纬度+椭球高，以厘米级精度转化工程坐标系
当工程坐标系定义为能明确时，7参数方法不可用。我们仅仅需要得知中央经线、纬度原点，以及3个控制点的经纬度静态测量结果及其相应的工程坐标，就可以用4参数方法，以经纬度+椭球高，求解厘米级精度的工程坐标系，满足GPR定位的要求。（3) 利用至少3个控制点，求解4参数，并利用4参数，将经纬度的横轴墨卡托投影结果转换到目标工程坐标系中；.COR文件记录的是某一道RTK定位的经度、纬度、椭球高。在进行轨迹分析前，必须将其转换成工程坐标系。（4）利用控制点进行多项式拟合，求解工程坐标正高。
2025-07-29 18:50:02
310
原创 对COR文件进行重排序、修改时间轴、重定义道数
在将两份COR文件拼接成一份文件的过程中，应该先将C->A的定位数据进行反转，然后再将C->B的定位数据添加到A->C测线的末端。（1）由于在拼接.RD3文件的过程中，雷达的道数被改变，因此在拼接COR文件时，应该重新定义定位数据所对应的雷达道数；处理雷达.RD3文件时，我们应将C->A测线反转成A->C，再将C->B添加到A->C测线的末端。（C->A.COR、C-->B.COR）B——年-日-月（也有可能是年-月-日，依据具体时间而定），UTC时间戳。C——时：分：秒，UTC时间戳。
2025-07-29 18:10:30
341
原创 虚拟串口+Python3，打通iRTK2到Mala地质雷达的数据链
（2）兼容转译层的解决思路为——（1）Python读取电脑串口的“GNGGA”格式数据；（2）“GNGGA”转译成“GPGGA”（3）使用软件“com0com”创建虚拟串口对A与B，需要注意的是，在Windows11系统下，经检查，问题为：iRTK2的5针串口输出“GNGGA”格式的定位信息，这是GNSS的数据格式；（1）iRTK2到Mala地质雷达的数据链为——iRTK2的5针串口输出（GNGAA）->（5针串口转7针串口）（GNGAA）->进入电脑串口（GNGAA）->软件读取（仅仅支持GPGAA）；
2025-05-25 20:33:33
428
原创 ISMIP-HOM冰川模型，以ISSM提供的程序为例
输入模型表面（DEM）、冰川厚度分布、冰川基底（DEM-冰川厚度）。需要注意的是，在模型表面的任意一处，冰川厚度都不能定义为0，这会导致模型不连续，所以我的建议是，冰川厚度为0的区域，都取值为0.0000001。（5）准备好数据DEM、冰川厚度、冰川表面速度长（请重采样到一致的行列数）。下面这个是依据冰川基底剪应力，结合冰川表面地形和DEM反演冰川厚度的代码，运行前请准备好SLOPE.tif、2013_V_Final.tif、DEM.tif——下面这是主程序的参数文件(ismipF.par)的代码——
2024-06-09 21:05:47
1797
原创 用shpfile切割envi 的Dat
shp_name="C:\Users\l_j_w\OneDrive\毕业设计\2013年\喀喇昆仑边界\喀喇昆仑边界.shp"y = (file_shp.data_range[2]-file_shp.data_range[0])/30;envi_file="D:\Lansat_Data_成品数据\2013_class"PIXEL_SIZE=[30.0D,30.0D], $;
2024-06-06 23:07:37
400
原创 生成冰川面积统计结果
pro Part_9 compile_opt idl2 e=envi(/headless) file_mkdir,"D:\Lansat_Data_综合报告" glacier_all = float(0.0) pixel_all_69 = float(0.0) folder = "D:\Lansat_Data_Ratio\" envi_arr = file_search(folder,'*.dat',count=num) arr_out_put = strarr(4,70)
2024-06-06 20:37:16
427
1原创 改善质量后的云层统计，第二部分
pro Part_7 compile_opt idl2 e=envi(/headless) cloud_all = float(0.0) pixel_all_69 = float(0.0) folder = "D:\Lansat_Data_Cloud_Final" envi_arr = file_search(folder,'*.dat',count=num) arr_out_put = strarr(4,70) for i=0,num-1,1 do begin en
2024-06-06 20:35:47
267
原创 接上一个，依据云量统计结果，完成图像拼接改善图像质量
由于上边的图像给出的是云检测结果的文件名，我们需要改成裁切后结果的文件名。开始冒泡排序法，对第一组数据，按照图像质量进行排序。开始冒泡排序法，对第二组数据，按照图像质量进行排序。开始冒泡排序法，按照图像质量进行排序。合并两个排序之后的数组得到运算结果。读取文件名、含云量、图像完整度。，删除前三个个数组之中的空元素。，删除前三个个数组之中的空元素。依据前方计数的有效元素数量。依据前方计数的有效元素数量。将第一张影像作为底图，后续。波段比值细化，选择更佳像素。次，按照波段比值公式。
2024-06-06 20:23:41
1585
原创 （接上一个程序）统计云层检测的二值化结果
proPart_4idl2e=envifork=1911DOBEGINkk=strtrimstring(k),1file_mkdir+kk+"\"endforfori=1911dobegin;if(ine86thencontinueii=strtrimstring(i),1folder=+ii+"\"filearr =(folder,'*.dat'long0long0forj=
2024-06-06 20:21:58
262
原创 使用波段运算阈值的方法，结合ENVI 的Band Math工具，批量，识别高原冰川地区的云层，并生成二值化图片
proPart_3idl2e=envifork=1911DOBEGINkk=strtrimstring(k),1file_mkdir+kk+"\"endforfori=1911dobegin;;if(ine86thencontinueii=strtrimstring(i),1folder=分区切割\"+ii+"\"filearr =(folder,'*.dat'forj=0。
2024-06-06 11:55:02
733
原创 用M个shapefile，切割N景遥感影像（两次循环）
proPart_2fork=1911DOBEGINkk=strtrimstring(k),1file_mkdir分区切割\"+kk+"\"endfore=envi波段合成'*.dat'fori=0,num-11dobegin;;OpenRasterforj=1911dobeginjj=strtrimstring(j),1shp_name=毕业设计\ROI_shp\"+jj+"\"
2024-06-06 11:50:43
360
原创 IDL&ENVI示例程序——批量将单波段的TIFF合成一个多波段TIFF
下载Landsat 8&9的数据的过程中，很多人往往选择一次性下载整个压缩包，包含这景影像的所有波段和文档。个文件夹中挑选出其中的3、4、、6波段，合成n张3波段遥感影像。个文件夹包含了Landsat 8下载的数据压缩包。但在某次研究中，我仅仅需要其中的三个波段即可。D:\Lansat_Data_原始数据文件夹，求出该路径下对应格式的影像文件路径。类型的文件，请检查路径是否正确。数据，各个波段单独成文件。若没找到对应文件，直接返回。求出波段合成顺序数组的大小。如果没找到文件，则返回。中的顺序构造波段数组。
2024-06-06 11:46:10
1840
原创 基于IDL语言和ENVI API的遥感数据批处理快速入门
pro”是procedure，意为“程序”，相当于C语言的“main”函数。其定义方法如下：proPart_1;这里面是一个空的主函数，分号是注释符号;IDL代码变量名不区分字母大小写;主函数名定义为”Part_1”的过程必须将代码文档命名为end以下是定义一个名为LayerStack、以parameter作为输入参数的函数：FunctionLayerStackreturn;返回参数ENDarray=[457。
2024-06-05 23:41:06
629
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人