- 博客(110)
- 资源 (1)
- 收藏
- 关注
原创 基于MOD02/MYD02获得亮度温度再转冰温
radiance_scales,和radiance_offset这两项参数代表波段的增益和偏移量,是辐射定标的系数。比如要计算31波段的辐射亮度,读取到scales为0.00084002,offsets为1577.33972168,带入MODIS辐射定标的通用公式:Radiance=scales*(DN-offsets),即可以得到该波段的辐射亮度。用HEG处理MOD02/MYD02,提取里面的EV_1KM_Emissive波段,band为11和12(其实就是band 31和32)。2. 辐射亮度转亮温。
2024-01-15 05:10:58
827
原创 python 已知经纬度获取影像DN值或位置
这两个读取出来的分别就是经度和纬度,但是读取极投影为3413的影像时,读取出来的时投影坐标,因此在程序中多做了一步变换。基本上就是一个经纬度转影像坐标的一个操作。
2023-12-02 01:48:06
325
原创 在import pyproj时候报错FileNotFoundError: [Errno 2] No such file or directory
在import pyproj时候报错FileNotFoundError: [Errno 2] No such file or directory: 'C:\\Users\\Asus\\AppData\\Roaming\\Python\\Python38\\site-packages\\pyproj.libs\\.load-order-pyproj-3.4.1'
2023-11-21 18:31:44
171
原创 python 调用HEG对MODIS数据进行批处理
只提取了MODIS数据当中的IST一个波段,输出成tif,其他什么都没做。在处理前判断这个文件是否已经被处理出tif。注意出来的tif是没有进行拉伸的,应该乘以0.01才是真正的IST的值。HEG其实可以批处理可以看我另外一篇博文,不需要写代码。但是对于300以上数量的MODIS影像非常容易自动停止,而且越来越慢。还是打算利用python每个调用。处理和直接调用HEG命名是一样的,注意处理前要先在HEG中生成一个参数化文件。
2023-11-15 18:07:38
624
原创 去除地理(遥感)图像的白边(无效值)
之前基于SNAP处理了Sentinel-1数据,发现地理校正之后,影像外围存在大量的nodata,使得文件大小大了一倍。对处理好的影像进行了裁剪,把白边去掉,并保留地理信息。
2023-07-03 11:48:59
518
原创 python pyroSAR对Sentinel-1数据进行预处理
调查了目前能够对Sentinel-1进行处理的python库,发现除了SNAP自带的snappy和另一个完全独立于SNAP的python库外,还有一个基于SNAP本身的XML,但是进行了优化,平均处理每一景2min。目前网上代码很少,因此做一个分享。
2023-07-03 11:18:22
1196
2
原创 python rasterio 基于矢量裁剪栅格包括属性筛选
rasterio的矢量裁剪的方法,但是真正运用的时候报错ValueError: No valid geometry objects found for rasterize,但是实际上我的shp和raster是重叠的,发现是没有将矢量数据的投影与栅格的进行统一。另外的一个需求是需要筛选shp的属性,直接shpdata.属性名即可获取属性。
2023-01-21 12:50:17
865
1
原创 Endnote中插入了DOI但是没有显示
除了引用格式的问题,还有就是文章里的引用可能是沿用之前的,并没有更新到本地数据库(意思是之前你使用了别的电脑进行插入,之前无DOI)采用一下方法,注意要把所有的tranveling library的文献都update了才行。...
2022-08-31 22:02:44
7461
原创 AIS数据下载并处理(python)
最近由于项目需求,需要下载AIS数据,其实也就是船舶位置数据。首先通过百度发现目前一些网站下载的AIS数据都是收费的,如船达通等等。船达通是只能查找近一个月的船舶轨迹信息,下载需要支付收费。但是最后我终于找到了一个免费下载AIS数据的网站,并且下面附上了处理说明,可以通过python代码处理。...
2022-08-31 14:48:54
14594
18
原创 利用ENVI的ROI建立mask图像
首先需要打开影像,然后在上面建立ROI利用Basic Tool-Region of Intertrest-Create Class Image from ROI可以将基于ROI的mask建立,具有和影像相同的大小和投影信息。生成的mask,maks值为1,非mask值为0
2022-04-20 16:31:37
1704
原创 利用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
4442
1
原创 IDL 对CMIP6数据进行投影+地理信息
在使用这个代码之前,先必须基于IDL读取CMIP6的NC数据,之前那一步生成了经度和纬度及变量三个矩阵tif文件。本代码可以对数据进行重新投影,并且给与地理信息。
2022-04-14 11:03:49
994
原创 MATLAB读取及存储geotiff
读取geotiff,对该数组进行操作之后进行存储,并且将地理信息一并存储进去[I, R] = geotiffread('H:\20170409\test\HHjiajihua\20170409_clip2_tif.tif');info = geotiffinfo('H:\20170409\test\HHjiajihua\20170409_clip2_tif.tif');geotiffwrite('H:\20170409\test\HHjiajihua\20170409_clip2_tif_geo
2022-04-07 19:51:03
5413
2
原创 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
5933
15
原创 Arcmap实现航线按日期分段
之前下载了航迹点,然后需要按照日期,对该航迹点进行表示,但是因为我只能使用一种颜色,因此最好是能够对每个日期分辨转为单个的线段,然后每个线段进行分隔,如下图所示因此首先在Arcmap-数据管理系统-要素-点转线,以下的红框选择日期,代表线按照日期生成之后,只需使用以下三类进行表示即可...
2021-09-05 21:00:42
253
原创 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
781
原创 利用python调用wget批量下载ftp文件(有筛选)
我的目的是对下面的这个ftp://sidads.colorado.edu/DATASETS/NOAA/G02202_V4/north/daily每年的文件夹进行遍历,下载1-2月的数据因此需要筛选出1-2月的文件利用了-A 这个参数,后面用-A "*198901*","*198902*" 这样就可以筛选出1989年的1-2月文件了,而且注意用逗号进行分割-r这个参数是递归下载-c这个参数是断点续传我需要对每年进行遍历,所以采用python调用cmd进行批量下载因为我的wget并.
2021-08-27 21:50:14
1407
原创 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
474
原创 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
675
转载 IDL 内存释放
最近在写的代码,当中涉及到提取数组中特定位置的数据,然后并保存。大概逻辑就是先从二维数组中提取某一行信息,再根据之前利用where得到的index,对这个行数组提取一些列
2021-06-05 16:59:43
1472
原创 NSIDC Data Download Script修改为批量下载
最近从NSIDC下载数据,想批量下载每一年的1-2月数据因此对他官方的代码修改如下只需要修改main方法即可其实主要是在每个循环中对filename_filter进行修改并且每次循环都要将url_list=None,否则后面就不会重新搜索文件def main(argv=None): global short_name, version, time_start, time_end, bounding_box, \ polygon, filename_filter,
2021-05-07 10:19:11
542
原创 python热力图
其实之前网上有很多代码主要是记录一下,首先是对日期的处理,把csv里面的日期,转换成python可以识别的日期,之后又把它转换成月份/日期我的数据是第一列为日期,其他列的列名都为年份sb.heatmap之后报错,后来发现需要将index设为日期以此记录import matplotlib.pyplot as pltimport numpy as npimport pandas as pdimport seaborn as sbfrom matplotlib import colo
2021-05-03 21:03:39
493
原创 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
1619
2
原创 gdal库安装成功但无法使用
之前能够正常运行,但是突然就不行了。但是利用pip安装显示成功后将源码中的import gdal改为以下from osgeo import gdal就可以运行了
2021-03-24 10:34:18
4168
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
1722
9
原创 python 散点图 多个Y轴
#该程序为了绘制散点图,即横坐标为时间,纵坐标有4个import matplotlib.pyplot as pltimport numpy as npimport pandas as pddata=pd.read_csv(r"H:\mission\speed\3_AIS_position_value.csv") #extent_LS.csvdata['datetime'] = pd.to_datetime(data['datetime'])c = data['piomas']c[c==-1.
2021-03-18 15:51:06
1366
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人