- 博客(108)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 基于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
535
原创 python 已知经纬度获取影像DN值或位置
这两个读取出来的分别就是经度和纬度，但是读取极投影为3413的影像时，读取出来的时投影坐标，因此在程序中多做了一步变换。基本上就是一个经纬度转影像坐标的一个操作。
2023-12-02 01:48:06
235
原创 在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
130
原创 python 调用HEG对MODIS数据进行批处理
只提取了MODIS数据当中的IST一个波段，输出成tif，其他什么都没做。在处理前判断这个文件是否已经被处理出tif。注意出来的tif是没有进行拉伸的，应该乘以0.01才是真正的IST的值。HEG其实可以批处理可以看我另外一篇博文，不需要写代码。但是对于300以上数量的MODIS影像非常容易自动停止，而且越来越慢。还是打算利用python每个调用。处理和直接调用HEG命名是一样的，注意处理前要先在HEG中生成一个参数化文件。
2023-11-15 18:07:38
517
原创 去除地理（遥感）图像的白边（无效值）
之前基于SNAP处理了Sentinel-1数据，发现地理校正之后，影像外围存在大量的nodata，使得文件大小大了一倍。对处理好的影像进行了裁剪，把白边去掉，并保留地理信息。
2023-07-03 11:48:59
364
原创 python pyroSAR对Sentinel-1数据进行预处理
调查了目前能够对Sentinel-1进行处理的python库，发现除了SNAP自带的snappy和另一个完全独立于SNAP的python库外，还有一个基于SNAP本身的XML，但是进行了优化，平均处理每一景2min。目前网上代码很少，因此做一个分享。
2023-07-03 11:18:22
872
2
原创 python rasterio 基于矢量裁剪栅格包括属性筛选
rasterio的矢量裁剪的方法，但是真正运用的时候报错ValueError: No valid geometry objects found for rasterize，但是实际上我的shp和raster是重叠的，发现是没有将矢量数据的投影与栅格的进行统一。另外的一个需求是需要筛选shp的属性，直接shpdata.属性名即可获取属性。
2023-01-21 12:50:17
733
1
原创 Endnote中插入了DOI但是没有显示
除了引用格式的问题，还有就是文章里的引用可能是沿用之前的，并没有更新到本地数据库（意思是之前你使用了别的电脑进行插入，之前无DOI）采用一下方法，注意要把所有的tranveling library的文献都update了才行。...
2022-08-31 22:02:44
6532
原创 AIS数据下载并处理（python）
最近由于项目需求，需要下载AIS数据，其实也就是船舶位置数据。首先通过百度发现目前一些网站下载的AIS数据都是收费的，如船达通等等。船达通是只能查找近一个月的船舶轨迹信息，下载需要支付收费。但是最后我终于找到了一个免费下载AIS数据的网站，并且下面附上了处理说明，可以通过python代码处理。...
2022-08-31 14:48:54
12032
14
原创 利用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
1497
原创 利用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
3941
1
原创 IDL 对CMIP6数据进行投影+地理信息
在使用这个代码之前，先必须基于IDL读取CMIP6的NC数据，之前那一步生成了经度和纬度及变量三个矩阵tif文件。本代码可以对数据进行重新投影，并且给与地理信息。
2022-04-14 11:03:49
794
原创 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
5103
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
5256
13
原创 Arcmap实现航线按日期分段
之前下载了航迹点，然后需要按照日期，对该航迹点进行表示，但是因为我只能使用一种颜色，因此最好是能够对每个日期分辨转为单个的线段，然后每个线段进行分隔，如下图所示因此首先在Arcmap-数据管理系统-要素-点转线，以下的红框选择日期，代表线按照日期生成之后，只需使用以下三类进行表示即可...
2021-09-05 21:00:42
225
原创 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
727
原创 利用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
1303
原创 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
439
原创 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
607
转载 IDL 内存释放
最近在写的代码，当中涉及到提取数组中特定位置的数据，然后并保存。大概逻辑就是先从二维数组中提取某一行信息，再根据之前利用where得到的index，对这个行数组提取一些列
2021-06-05 16:59:43
1390
原创 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
474
原创 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
462
原创 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
1491
2
原创 gdal库安装成功但无法使用
之前能够正常运行，但是突然就不行了。但是利用pip安装显示成功后将源码中的import gdal改为以下from osgeo import gdal就可以运行了
2021-03-24 10:34:18
3740
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
1526
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
1298
原创 python 散点图 自定义x轴标签/经度180到-180连接
最近在做一个散点图，横轴为经度，纵轴为速度，即不是纬度，因此不是那种展现空间关系的。但是数据的经度为70到180，之后-180到-170这样，因此用正常的方法出图会出现，-180和180的散点无法进行连接。于是使用python修改标签的办法。首先将横轴的数据即经度的-180~180变为0-360利用以下函数实现def changto360(data): c = data['lon'] c[c<0]=c+360 data['lon']=c return
2021-03-18 15:46:42
1611
1
原创 MATLAB数组的横向和纵向拼接
横向拼接，因为我这里有元素为元祖 ，用了{}outtemp={date{1} position(i,1) position(i,2) position(i,3) thickness con};纵向拼接result=[result;outtemp];
2021-03-10 10:27:53
2435
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人