李百柏-CSDN博客

原创 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 1373

原创用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 253

原创生成冰川面积统计结果

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 363

原创波段比值法统计冰川

pro Part_8endifendifendforend。

2024-06-06 20:36:41 288

原创改善质量后的云层统计，第二部分

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 213

原创统计质量改善后的遥感影像的结果——第一部分，云层检测二值化

pro Part_6endforendforend。

2024-06-06 20:33:42 376

原创接上一个，依据云量统计结果，完成图像拼接改善图像质量

由于上边的图像给出的是云检测结果的文件名，我们需要改成裁切后结果的文件名。开始冒泡排序法，对第一组数据，按照图像质量进行排序。开始冒泡排序法，对第二组数据，按照图像质量进行排序。开始冒泡排序法，按照图像质量进行排序。合并两个排序之后的数组得到运算结果。读取文件名、含云量、图像完整度。，删除前三个个数组之中的空元素。，删除前三个个数组之中的空元素。依据前方计数的有效元素数量。依据前方计数的有效元素数量。将第一张影像作为底图，后续。波段比值细化，选择更佳像素。次，按照波段比值公式。

2024-06-06 20:23:41 1503

原创（接上一个程序）统计云层检测的二值化结果

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 209

原创使用波段运算阈值的方法，结合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 299

原创用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 291

原创 IDL&ENVI示例程序——批量将单波段的TIFF合成一个多波段TIFF

下载Landsat 8&9的数据的过程中，很多人往往选择一次性下载整个压缩包，包含这景影像的所有波段和文档。个文件夹中挑选出其中的3、4、、6波段，合成n张3波段遥感影像。个文件夹包含了Landsat 8下载的数据压缩包。但在某次研究中，我仅仅需要其中的三个波段即可。D:\Lansat_Data_原始数据文件夹，求出该路径下对应格式的影像文件路径。类型的文件，请检查路径是否正确。数据，各个波段单独成文件。若没找到对应文件，直接返回。求出波段合成顺序数组的大小。如果没找到文件，则返回。中的顺序构造波段数组。

2024-06-06 11:46:10 875

aoaoooo的博客