idl
就是一只白
更新频率取决于干活的频率。
展开
-
ERA5温度、风速风向处理
后面增加个地理校正批处理。原创 2023-09-13 16:40:46 · 1124 阅读 · 0 评论 -
IDL DOY转年月日
基于IDL转DOY为月日原创 2023-06-07 22:01:42 · 121 阅读 · 0 评论 -
IDL 打开ENVI标准文件,并存为tif
打开ENVI标准文件并另存为TIF格式原创 2022-06-09 15:16:40 · 1580 阅读 · 0 评论 -
利用ENVI的ROI统计影像信息
首先在ENVI当中建立ROI,然后保存成XML格式利用下面代码的raster = e.OpenRaster(file)和raster.ExportRois,infilepath+'incidence_angle_sample.csv' , rois, 'CSV'就可以实现按照ROI统计的影像信息的CSV打开生成的CSV后得到RO1 #1 类别 npts(像素个数): 320RO1 #2 类别 npts(像素个数): 320B1、B2、B3为影像的波段1、2和3的像素值..原创 2022-04-17 18:45:22 · 4402 阅读 · 1 评论 -
IDL 对CMIP6数据进行投影+地理信息
在使用这个代码之前,先必须基于IDL读取CMIP6的NC数据,之前那一步生成了经度和纬度及变量三个矩阵tif文件。本代码可以对数据进行重新投影,并且给与地理信息。原创 2022-04-14 11:03:49 · 971 阅读 · 0 评论 -
IDL 读取ICESAT2 HDF5文件
pro read_icesat2_hdf5 infilepath='E:\icesat2\5000001293407\' outfilepath='' cd,infilepath thesefiles = file_search('*',count = numFiles) print,'number of files found:',numFiles FOR fidx=0,numFiles-1 DO BEGIN cd,infilepath+thesefiles[fidx.原创 2021-08-12 15:24:06 · 470 阅读 · 0 评论 -
IDL CMIP6 NC格式数据处理
利用IDL读取CMIP6的NC文件数据,如果想要投影转换,把数据进行地理投影转换可以看另一篇文章原创 2021-06-19 20:27:56 · 3599 阅读 · 4 评论 -
IDL 拆分和合并字符串 STRSPLIT与STRJOIN
s='siconc_SIday_AWI-CM-1-1-MR_historical_r1i1p1f1_gn_18500101-18501231'str=STRSPLIT(s, '_', /EXTRACT)str=原创 2021-06-19 15:45:59 · 660 阅读 · 0 评论 -
IDL 内存释放
最近在写的代码,当中涉及到提取数组中特定位置的数据,然后并保存。大概逻辑就是先从二维数组中提取某一行信息,再根据之前利用where得到的index,对这个行数组提取一些列转载 2021-06-05 16:59:43 · 1454 阅读 · 0 评论 -
PIOMAS Arctic sea ice thickness daily数据下载和IDL处理
网上好像没有PIOMAS数据处理的文章,最近项目需要进行了处理,因此给大家提供以下参考,不过自己摸索一下应该很快就会了PIOMAS数据下载http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model_grid在表格中的Sea ice thickness (Volume per unit Area) , daily mean即为我下载的海冰厚度数据点击右侧的 hiday.H<yyyy&g.原创 2021-03-18 16:08:44 · 1706 阅读 · 9 评论 -
nc、hdf等数据转为带地理信息的geotif
在做工作的时候会碰到很多nc格式的数据,并且读取的数据文件是长时间序列的,具有相同的地理投影与地理范围。通常会用IDL来读取数据,但是读取的数据如何变为带有地理信息的数据,一直都是个问题。以下是我的常用方法和尝试了但个人用起来不方便的方法(很主观,勿喷)一般情况,除了data,也就是你获取的变量数据之外,还会有和data相同矩阵大小的经度和纬度矩阵。这两个是建立地理信息必要的东西。1、利用arcmap打开nc数据,此时如果能够成功获取具有地理信息的数据,再保存成tif文件。在IDL中利用ge原创 2021-02-09 21:55:52 · 5373 阅读 · 3 评论 -
IDL 处理ERA5的NC格式的温度及海平面气压数据
首先是我的数据下载方面,从官网上直接下载没有使用代码,因此所有下载的数据都放到一个NC文件中了(当然2m温度和气压是分开的)我主要下载的区域是北纬30-90,每日只下载四个时间,0 6 12 18,也只下载6个月的,1-3、10-12月nc文件可以尝试使用panoly打开看看,每个月和每个时间是按顺序排好的我利用以下的代码来处理2m温度和海平面气压数据值得注意的点是,在读取数据时,还要读取校正系数和偏差,即下面的a和b变量,并且注意读取时数据格式是int,运算时需要先对数组乘以1.0变成fl原创 2020-12-02 10:13:18 · 2509 阅读 · 0 评论 -
IDL 合并一维二维三维数组
一维数组(行数相同)array3 = [array1,array2]二维数组(列数相同)array3 = [[array1],[array2]]三维数组(行列数都相同)array = [[[array1]], [[array2]]]原创 2019-10-04 12:46:31 · 2367 阅读 · 2 评论 -
IDL将ENVI ROI生成的TXT转为CSV
由于ENVI选取ROI后可以导出一个TXT但是并不是CSV,所以希望生成一个仅有各个特征值和类序号的一个CSV。参看了之前的一个博客,但是我这个可以自动识别特征值有几个构成的,就是ROI输出的波段有几个。希望给大家一点参考吧PRO GETROI ;读取ENVI导出训练区的txt文件 filename = dialog_pickfile(path = 'H:\IceClas...原创 2019-09-24 21:17:55 · 1485 阅读 · 0 评论 -
IDL制legend小整理
今天做了超级多的图,顺手粘一下legend的code,后面再整理 labels = ['< 0','0-50', '50-100', '100-150', '150-200', '200-250', $ '250-300', '300-350', '350-400', '>400'] cgDCBar, NColors=10, Bottom=1, Color='black'...原创 2018-05-17 12:30:51 · 825 阅读 · 0 评论 -
IDL存储和读取类表格的txt
表示不知道怎么存储成csv或者Excel能打开的数据所以写了一个code能存txt格式的表格数据,但是用Excel也能自动打开最后输出的txt和读取的txt格式长这样: 1 151 2 152 3 153 4 154 5 155 6 156 7 157 8 158 9 15...原创 2018-05-06 21:17:47 · 5124 阅读 · 0 评论 -
IDL输出txt换行的问题
直接用printf输出的txt每行的内容有限制,所以百度了一下设置openw函数的width关键字,该关键字决定了打印的宽度(column),默认情况下width=80,也就是打印80个字符后换行。如果你要每行输出六十个数,那么width应该大于等于(数字位数+数字后空格/制表符数)*60....转载 2018-05-07 16:21:50 · 3508 阅读 · 0 评论 -
IDL图形保存图片时文字变黑块的临时解决方法
忙着做毕设没想到出了这么个问题,生成的图被打了码。。。哭泣在网上搜到了这么个办法,转载至连接 测试代码如下,注意其中的xtext_orientation、ytext_orientation、baseline的用法。坐标轴标题使用orientation属性,图形标题使用baseline属性(因为此标题是一个TEXT)。 IDL> p=plot(/test, title...转载 2018-05-20 21:05:29 · 470 阅读 · 0 评论 -
IDL制图配色
之前利用IDL制图老是没有办法把colorbar里面的颜色和实际数值在图里的颜色对应起来,今天好像知道怎么弄了,先记录一下。我分的-4-4区间,分成8级,颜色的读入也是读8个,则那么ncolor就要写8+1B。...原创 2018-09-18 09:55:19 · 818 阅读 · 0 评论 -
IDL 读取nc数据代码
file_ID = NCDF_OPEN(theseFiles[fidx],/NOWRITE) ;open netCDF file for READ only Tag = NCDF_INQUIRE(file_ID) timeid = NCDF_VARID(file_ID,'time') ;initial_time0initial_time0_hours ...原创 2018-10-02 16:41:12 · 3815 阅读 · 0 评论 -
IDL读取写入csv/tif/dat文件
读数据:WRITE_CSV,outfilename,dataread_tiff(outfilename,geotiff=geoinfo);读取csv;也可以这样thesefiles = FILE_SEARCH('ice_multiplier_table.csv') table=intarr(9,7) openr,lun,thesefiles[0],/get_lun...原创 2018-10-04 10:53:46 · 10983 阅读 · 1 评论 -
IDL的mean计算
关于mean的DIMENSION当dimension设为1时,计算行的均值,最后得到一个纵列,当设为2时计算列的均值,最后得到一个行设为3时,可以得到时序(3维)上的均值,结果是一个二位数组。...原创 2018-10-04 13:46:49 · 3420 阅读 · 0 评论 -
利用IDL读取HDF5文件
首先是要了解HDF5文件的结构,利用pannoly文件能快速浏览文件的信息和命名,pannoly里面显示的数组行列和波段和IDL反过来的,IDL是 数组[列][行][波段] pannoly是数组(波段,行,列)然后就是读取file_id = H5F_OPEN('文件地址');打开Data这个群组group_id=H5G_OPEN(file_id,'Data');H5D_OPEN打开...原创 2018-10-17 19:34:50 · 3370 阅读 · 0 评论 -
CGPS2RASTER: ImageMagick cannot be found on this machine.
想利用IDL出图,在把ps文件转成png文件时报错报错的提示:ImageMagick没有安装下载地址:http://www.imagemagick.org/script/download.php但是我第一次安装之后还是显示报错,后来看了http://www.idlcoyote.com/graphics_tips/weboutput.php里面所提到的convert命令我在安装的文件...原创 2018-11-06 09:37:56 · 710 阅读 · 3 评论 -
IDL与geotiff制作相同地理范围的地图
pro geotiff_map_test compile_opt idl2 forward_function TS_De_Anomaly forward_function TS_Anomaly;超快速的geotiff; file = 'G:\sea_ice_study\smos_tb_2018_12_01_geo.tif'; CO=cgGeoMap(file, /Conti...原创 2019-03-04 09:32:04 · 921 阅读 · 0 评论 -
提高IDL最短路径算法运行效率
前两天跑了一个最短路径,在找出需要搜索的点那放了两个循环,因此大约2分钟或者一分半出一张361*361的图,其实搜索的范围就这个图的一半,我当时觉得已经算是在降低运行时间了,但是由于跑的实在太慢,而我毕设的数据量又很大,因此再进行改进。 是想着把那两个循环去掉,其实那两个循环的目的很简单就是找到distance数组里最小的且没有被搜索过的数,因此用了数组运算,将符合条件的像元都存...原创 2018-05-10 08:49:22 · 979 阅读 · 0 评论