data process
就是一只白
更新频率取决于干活的频率。
展开
-
gdal转投影
【代码】gdal转投影。原创 2023-12-12 17:21:57 · 145 阅读 · 0 评论 -
python 已知经纬度获取影像DN值或位置
这两个读取出来的分别就是经度和纬度,但是读取极投影为3413的影像时,读取出来的时投影坐标,因此在程序中多做了一步变换。基本上就是一个经纬度转影像坐标的一个操作。原创 2023-12-02 01:48:06 · 345 阅读 · 0 评论 -
python 调用HEG对MODIS数据进行批处理
只提取了MODIS数据当中的IST一个波段,输出成tif,其他什么都没做。在处理前判断这个文件是否已经被处理出tif。注意出来的tif是没有进行拉伸的,应该乘以0.01才是真正的IST的值。HEG其实可以批处理可以看我另外一篇博文,不需要写代码。但是对于300以上数量的MODIS影像非常容易自动停止,而且越来越慢。还是打算利用python每个调用。处理和直接调用HEG命名是一样的,注意处理前要先在HEG中生成一个参数化文件。原创 2023-11-15 18:07:38 · 636 阅读 · 0 评论 -
ERA5温度、风速风向处理
后面增加个地理校正批处理。原创 2023-09-13 16:40:46 · 1243 阅读 · 0 评论 -
去除地理(遥感)图像的白边(无效值)
之前基于SNAP处理了Sentinel-1数据,发现地理校正之后,影像外围存在大量的nodata,使得文件大小大了一倍。对处理好的影像进行了裁剪,把白边去掉,并保留地理信息。原创 2023-07-03 11:48:59 · 539 阅读 · 0 评论 -
python pyroSAR对Sentinel-1数据进行预处理
调查了目前能够对Sentinel-1进行处理的python库,发现除了SNAP自带的snappy和另一个完全独立于SNAP的python库外,还有一个基于SNAP本身的XML,但是进行了优化,平均处理每一景2min。目前网上代码很少,因此做一个分享。原创 2023-07-03 11:18:22 · 1266 阅读 · 2 评论 -
建立矢量缓冲区+合并矢量要素
程序包含把单个shapefile所有的要素都合并在一起,基于shapefile建立缓冲区两个功能原创 2023-03-03 09:42:02 · 455 阅读 · 1 评论 -
python rasterio 基于矢量裁剪栅格包括属性筛选
rasterio的矢量裁剪的方法,但是真正运用的时候报错ValueError: No valid geometry objects found for rasterize,但是实际上我的shp和raster是重叠的,发现是没有将矢量数据的投影与栅格的进行统一。另外的一个需求是需要筛选shp的属性,直接shpdata.属性名即可获取属性。原创 2023-01-21 12:50:17 · 886 阅读 · 1 评论 -
IDL 对CMIP6数据进行投影+地理信息
在使用这个代码之前,先必须基于IDL读取CMIP6的NC数据,之前那一步生成了经度和纬度及变量三个矩阵tif文件。本代码可以对数据进行重新投影,并且给与地理信息。原创 2022-04-14 11:03:49 · 1039 阅读 · 0 评论 -
GEE Sentinel-1影像数据下载及处理
单独下载一张影像,影像信息如代码里的//导出HH HV 和angle//下载指定文件名称的S1数据,导出为北极极投影,40m,双极化数据,波段为HH HV 入射角var image = ee.Image.load('COPERNICUS/S1_GRD/S1A_EW_GRDM_1SDH_20210429T152824_20210429T152924_037667_0471AC_7DBB');// image = image.select("VH");//转成float32位image = im原创 2022-02-19 21:29:09 · 6020 阅读 · 15 评论 -
dataframe按照时间均匀采样
目标是希望一个dataframe按照其中的一列,datatime,按照6个小时为间隔进行均匀采样,即按照0、6、12、18点,进行时间段内的求取均值data=data.resample('6H', on='datetime').mean()data=data.reset_index(drop=False)重新将index去掉,因为后面需要访问data['datatime']报错了...原创 2021-09-05 20:49:46 · 792 阅读 · 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 · 484 阅读 · 0 评论 -
IDL CMIP6 NC格式数据处理
利用IDL读取CMIP6的NC文件数据,如果想要投影转换,把数据进行地理投影转换可以看另一篇文章原创 2021-06-19 20:27:56 · 3668 阅读 · 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 · 692 阅读 · 0 评论 -
IDL 内存释放
最近在写的代码,当中涉及到提取数组中特定位置的数据,然后并保存。大概逻辑就是先从二维数组中提取某一行信息,再根据之前利用where得到的index,对这个行数组提取一些列转载 2021-06-05 16:59:43 · 1490 阅读 · 0 评论 -
PIOMAS月海冰厚度数据批量处理+建立地理信息(IDL的IGM建立地理信息)
首先我的程序分为两个部分,先把PIOMAS的数据提取出来pro read_piomas_monthlynv = envi(/HEADLESS) infilepath = 'H:\mission\PIOMAS\monthly\' outfilename='H:\mission\PIOMAS\monthly\hdr\' CD,infilepath thesefiles = FILE_SEARCH('heff*',count = numFiles) ;解压 FOR fidx=...原创 2021-03-25 09:58:20 · 1643 阅读 · 2 评论 -
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 · 1752 阅读 · 9 评论 -
AIS数据下载并处理(python)
最近由于项目需求,需要下载AIS数据,其实也就是船舶位置数据。首先通过百度发现目前一些网站下载的AIS数据都是收费的,如船达通等等。船达通是只能查找近一个月的船舶轨迹信息,下载需要支付收费。但是最后我终于找到了一个免费下载AIS数据的网站,并且下面附上了处理说明,可以通过python代码处理。...原创 2022-08-31 14:48:54 · 15257 阅读 · 18 评论 -
nc、hdf等数据转为带地理信息的geotif
在做工作的时候会碰到很多nc格式的数据,并且读取的数据文件是长时间序列的,具有相同的地理投影与地理范围。通常会用IDL来读取数据,但是读取的数据如何变为带有地理信息的数据,一直都是个问题。以下是我的常用方法和尝试了但个人用起来不方便的方法(很主观,勿喷)一般情况,除了data,也就是你获取的变量数据之外,还会有和data相同矩阵大小的经度和纬度矩阵。这两个是建立地理信息必要的东西。1、利用arcmap打开nc数据,此时如果能够成功获取具有地理信息的数据,再保存成tif文件。在IDL中利用ge原创 2021-02-09 21:55:52 · 5500 阅读 · 3 评论 -
MATLAB将多个NETCDF数据合并/MATLAB写NETCDF数据
帮同学写的一个程序,目的是:把很多个NC数据合并在一起,并且每个.nc数据的变量名称是相同,我平常读.nc数据比较多,较少涉及写入。我先使用了IDL但是由于原数据的属性信息也很多,因此最后放弃。而MATLAB有基于现有的数据创建一个模板,对数据进行生成。代码如下:1.先对需要的.NC数据进行提取2.提取所需的数据属性等3.对数据进行赋值ticclose all;clear...原创 2019-12-06 12:28:35 · 2832 阅读 · 0 评论 -
Sentinel 1文件名命名规则
Sentinel 1文件名命名规则如下图:官网的详细说明:https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1-sar/naming-conventions以下为例子:S1A_EW_GRDM_1SDH_20180112T082556_20180112T082700_020119_0224E3_7F9BS1A...原创 2019-03-28 14:21:06 · 10133 阅读 · 11 评论 -
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 · 3853 阅读 · 0 评论 -
关于处理长时间序列遥感数据的小贴士——不定时更新
0.同期数据替换的时候,一般会先从文档找可信度比较高的那一期数据。之后利用另一期数据填补选择的那期数据的缺失值,此时一定要对比两个数据的差= =然后再考虑一下1.空间滤波有时候会扩大之前数据异常值的的范围,特别是均值滤波,要小心。2.拿到数据的时候,检查投影和缺失值异常值的情况3.S-G滤波本就已经避开了异常值,因此滤波只用对检查出的异常值进行取代就好,不用对所有时间序列进行滤波处理。...原创 2018-05-14 10:28:53 · 2518 阅读 · 0 评论