自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

炒茄子的博客

致力于遥感技术的传播,不装逼我浑身难受啊

  • 博客(210)
  • 资源 (1)
  • 收藏
  • 关注

原创 Python:如何对FY3D TSHS的数据集进行重投影并输出为TIFF文件以及批量镶嵌插值?

由于代码中注释较为详细,因此博客中部分操作一笔带过。

2024-04-09 12:28:52 1209

原创 遥感&深度学习:CNN-LSTM模型用于NDVI的预测(Pytorch代码深度剖析)

这是一次完整的关于时空遥感影像预测相关的深度学习项目,后续有时间更新后续部分。pytroch的模型的基本使用包括LSTM和Linear等的框架构建、训练评估和预测;遥感影像如何传入模型中训练,包括数据加载和预输入的要求等,对于大型数据集如何解决内存不足问题(可能会更新);遥感图像的处理例如裁剪掩膜,nc读取相关数据集并glt校正输出为tiff文件等诸多处理。Landuse(年土地利用): 2001 - 2020LST(月均地表温度): 200002 - 202210NDVI。

2024-03-15 22:27:52 2040 2

原创 Python:如何将MCD12Q1\MOD11A2\MOD13A2原始数据集批量输出为TIFF文件(镶嵌/重投影/)?

之前一段时间一直使用ENVI IDL处理遥感数据,但是确实对于一些比较新鲜的东西IDL并没有python那么好的及时性,封装的东西也较为底层需要自己实现的东西相对python还是不那么方便,当然它也训练了我对于数据的处理能力。另外最近一直在探索深度学习,IDL相对小众,在时间不充裕的情况,我将倾向于使用Python进行学习而非IDL。因此,目前相当长一段时间我将同时兼容Python和IDL两门语言,但正如你所见,IDL的趋势正在减弱。

2023-12-16 00:15:58 1472

原创 ENVI IDL:如何MODIS GRID产品进行批量镶嵌、重投影(GLT校正)?

时间类的定义相关概念和基本使用查看:https://blog.csdn.net/m0_63001937/article/details/133975751;初始化初始化父类对象return, 1;1表示成功实例化对象end;若不传入参数则默认该格式化字符串end;更新endif;年积日转化为年月日endend;更新endif;年月日转年积日endend;更新年月日等属性endifendifend;取值;如果用户请求返回某一属性, 那么将其返回。

2023-11-06 12:55:46 885 5

原创 ArcMap:第二届全国大学生GIS技能大赛(广西师范学院)详解-上午题

现有某区域一幅有投影信息的单波段遥感图像(tm00.img)、DEM数据、两幅多光谱遥感图像(tm01.img\tm02.img)及该区域的矢量数据图层文件(bound.shp),其中有一幅多光谱图像无投影信息(tm02.img),现需要对无投影信息的多光谱图像进行几何校正、镶嵌、裁切等处理,并利用矢量数据裁剪镶嵌后的多光谱数据,利用裁剪后的多光谱数据及DEM数据,提取区域植被覆盖信息。朝北缓坡植被:NDVI大于0.25,坡度小于20°,朝北非朝北缓坡植被:NDVI大于0.25,坡度小于20°,非朝北。

2023-10-08 20:59:04 998 2

原创 ENVI IDL:OMI-NO2产品均值计算与可视化(附Python代码)-年、季、月均值计算

月等信息,然后分配到不同的[]中,[]在图示中不是数组而是列表(在IDL中为链表,实际类似于Python的列表), 由于我们事先不知道春季有多少文件,所以定义列表更为合适(只需要add添加元素即可)。其是通过计算每一像元的累积和以及该像元的实际有效天数或者次数进行均值的求取,对于像元的有效次数的求取在我的代码中通过mean函数的/nan参数解决,自动跳过无效值(所以需要事先将所有无效值替换为NAN),其余基本一致。春季为3-5月,夏季为6-8月,秋季为9-11月,冬季为12-2月。

2023-09-19 22:56:41 340

原创 GDAL:Warning 1: All options related to creation ignored in update mode

警告代码在gdal.Warp函数中产生。警告:,很明显,正如它所说,所有和创建相关的options将会在更新模式下忽略。这里的options根据gdal。

2024-05-09 21:58:30 586

原创 简记:机器学习中关于训练集和测试集的标准化问题

错误处理:在未进行数据集划分之前,对所有特征项(X/因子项)进行标准化例如MinMaxScaler\Normalizer等方法,然后再进行数据集划分;或者在数据集划分之后分别对训练集和测试集进行独立的标准化。

2024-04-29 16:39:30 472

原创 Python:如何使用GDAL实现多波段遥感影像的批量拼接?

设置拼接好的遥感影像的缺失值(这里就还是照旧0.0吧,如果说你觉得0.0可能是有意义的,那么你可以设为例如-9999等)本打算自己写一个拼接函数,但是感觉太折腾了,查了一下gdal的实现方法,发现有封装的非常好的一个方法(:原始数据集(待拼接影像)的缺失值数值。(例如你的遥感影像中0.0为缺失值,那么就传入0.0):列表形式,传入所有需要进行拼接的tiff文件的路径;例如此处需要拼接的tiff文件有两个。一行代码就可以解决,需要注意的。)可以解决,直接就拿来用了。:拼接好的影像的输出路径。

2024-03-17 10:52:47 604

原创 Python:如何统计特定返回周期下的GEV分布值和实际观测值的超越概率?

例如,某个站点的在某天发生了降水事件,日降水值是300mm,那么我们通过GEV分布得知它的非超越概率是95%,那么超越概率就是5%,也就是说,这次降水事件发生的概率是5%,换言之,在所有降水事件,几乎有95%的降水事件的强度都小于此次降水的强度。累积分布函数基于得到的GEV分布特征,对传入的单列数据的每一值进行计算,给出对于当前数据值的非超越概率(*100 百分化),再用1减去意为超越概率,也就是说,当前值发生的概率有多大。在对单列数据进行GEV分布拟合,得到该分布的分布特征,然后使用。

2024-03-15 17:39:59 927

原创 R语言:如何基于地球外辐射(Ra)和相对日照(n/N)计算太阳辐射Rs?

这些常量会在后续的计算中使用到。pi,G

2024-03-15 15:37:21 1259

原创 Python: 如何绘制核密度散点图和箱线图?

这是数据样式:Note:数据中存在无效值NA(包括后续的DEM),需要注意。

2024-03-12 18:20:15 639 1

原创 SNAP:如何批量预处理Sentinel2 L2A数据集并输出为TIFF文件?

我目前就是希望下载哨兵2号数据,然后在SNAP中进行批量提取真彩色波段并输出为TIFF文件。

2024-03-03 22:59:26 986 1

原创 IDL:如何对地表温度数据进行分组统计并生成温度分布直方图?

获取一张地表温度数据,对其进行分组统计,并制作间隔为3K的温度分布频数和频率直方图。

2024-02-28 23:53:35 511

原创 Wget: 如何批量下载风云卫星和NASA数据?

首先进入风云卫星服务官网: https://satellite.nsmc.org.cn/portalsite/default.aspx假定我们下载火点数据,如下先找到下载界面:接着依据数据所对应卫星\传感器\产品名称进行检索:选定需要的数据加入到购物车中:接着提交订单等待所有数据回调完成就可以下载:找到下载连接并另存为txt文件:下面为txt文件的链接:意为唤醒wget程序起来工作了;-i意为对文本文件中的各个下载链接进行迭代逐个下载;-P。

2024-02-26 23:27:36 729

原创 论文解析:NDVI Forecasting Model Based on the Combination of Time Series Decomposition and CNN – LSTM

论文地址:https://link.springer.com/10.1007/s11269-022-03419-3这篇论文是想要基于时间序列分解和CNN-LSTM两种方法的混合进行NDVI的预测,需要注意的是此处的CNN和LSTM你可以拆分看,他们并不是一个完整的概念,只是二者的简单叠加即CNN的输出作为LSTM的输入。我就不讲这些算法到底好亦或者差,仅仅从实现原理上进行解析。仅个人解读。

2023-12-30 22:41:18 1009

原创 Python:(Sentinel-1)如何解析SNAP输出的HDF5文件并输出为GeoTIFF?

最近在了解sentinel-1的预处理过程,但是由于影响太大了,常规的GeoTIFF无法输出预处理结果,BigTIFF输出时似乎也遇到了一些问题(好在后面解决了,所以正好做一下HDF5文件输出的TIFF文件与BigTIFF文件的对比),对于输出的HDF5文件则完全没有问题。但是问题在于HDF5文件的结构尚不了解,因此对于其中的地理信息如何提取很关键(当然你可以使用ArcGIS或者ENVI打开其中的VV和VH波段,但是都无法自动读取到其中的地理信息或者坐标系信息)。

2023-12-18 23:01:17 1299 2

原创 Python:如何解决MODIS GRID(正弦投影/GCTP_SNSOID)的重投影问题?

GIS茄子;若博客出现纰漏或有更多问题交流欢迎关注GIS茄子,或者邮箱联系(推荐-见主页).

2023-12-15 12:26:42 1447 3

原创 ENVI IDL:如何提高程序运行效率?

最近突然(并不是)想深入了解IDL底层的一些原理,因此记录一下感受和体会。

2023-12-12 13:01:36 770

原创 遥感&深度学习:如何读取TIFF文件和切片成256*256?

最近打算认真从头开始学习深度学习和遥感结合的相关内容,主要通过Python进行处理。此前用深度学习进行遥感相关的学习一直都是用框架,但是考虑很多因素我后面打算换用Pytorch进行学习。好在tensorflow我只是浅尝辄止,不然还真有点不舍。

2023-12-08 15:34:14 881

原创 ArcEngine:如何进行缩放图层、属性信息显示、状态栏显示?

如果以后不是工作相关,我或许以后再难了解ArcEngine相关的GIS桌面软件二次开发的内容?

2023-11-20 22:21:47 510 1

原创 ENVI IDL:如何生成FY4A快照

FY4A全圆盘(DISK,全球)多光谱影像,panoply软件打开数据层次结构如下:我们生成快照主要使用到其中的进行快照显示,注意我并没有进行辐射定标。

2023-11-13 18:19:21 226

原创 ENVI IDL:如何编写多IF-ELSE结构?

最近在整理代码框架结构,对于之前的一些逻辑框架进行重新梳理,我一度以为在IDL中并没有设计多IF-ELSE结构,只能单IF结构或者IF-ELSE结构,我尝试从帮助中寻找相关多IF-ELSE结构,但似乎并没有结果,暂作罢。

2023-11-13 10:53:10 261

原创 ENVI IDL:如何将txt文本文件转化为GeoTIFF文件?

此处的文本文件形式如下:里面包含了众多点位信息(不是站点数据),我们需要依据上述点的经纬度信息放到对应位置的像素点位置,放置完后如下:可以发现,还存在部分缺失值,我们还需要进行缺失值的填补。

2023-11-12 17:19:13 403

原创 ArcGIS:如何迭代Shp文件所有要素并分别导出为Shp文件?

尝试用IDL实现,奈何又涉及新的类IDLffShape,觉得实在没有必要学习的必要,毕竟不是搞开发,只是做做数据处理,没必要拿IDL不擅长的且底层的东西自己造轮子。这里想到使用Python去解决,gdal太久没用有些生疏了;Arcpy虽然也可以,但是似乎也太麻烦了,且不常用。就用ArcGIS模型构建器吧,无脑操作就好了,不用动脑子的快乐。

2023-11-12 00:30:17 1241 5

原创 ENVI IDL:如何基于气象站点数据进行反距离权重插值?

仅仅练习,大可使用ArcGIS或者已经封装好的python模块进行插值,此处仅仅从底层理解如何从公式和代码理解反距离权重插值的过程,从而更深刻的理解IDL的使用和插值的理解。

2023-11-11 21:23:25 738

原创 ENVI IDL:如何解析XML文件(以Landsat9-MTL.xml文件为例)

我们原本是打算对Landsat9文件进行辐射定标,但是辐射定标的参数在MTL文件中,从文件中查看参数直接复制到IDL中固然可行,但是当我们对Landsat9文件进行批量辐射定标时,这种方法就将失效了。因此我们需要自动从MTL文件中读取相关参数,这里的相关参数实际上只包含两个参数(对于一个波段),一个是比例系数,一个是偏置量。对于Landsat9,给出三种MTL形式:这里我们只讨论txt文本文件和XML文件的解析和提取。

2023-11-11 19:58:36 837

原创 ENVI IDL:如何监测代码运行时间(计时器函数实现)?

毫无疑问,timer_keep函数应该存储全局变量,否则无法记录相邻两次调用的间隔时长。查找IDL的全局变量设置方法,似乎并没有全局变量的设置方法,约莫半晌。关键字,称为公共块,不同IDL进程或者函数均可以访问该块内的变量,公共块在进程结束时不会被销毁。需要使用到两次systime,分别为循环开始和结尾各一次,是否可以只调用一次函数解决呢?如果想要解决第一次运行为0的问题,那么解决办法就是循环开始和结束都使用一次。缺点就是,第一次调用该函数,输出是0s,这似乎并不好解决。暂且如此,bye~.

2023-11-11 09:59:03 123

原创 ENVI IDL: 如何进行对MODIS GRID数据进行重投影?

其包含三个pro过程,其中为主程序,img_warp为常规pro过程,read_h4为自定义函数。基于角点信息进行投影坐标数据集的生成并转化为经纬度数据集,img_warp基于经纬度数据集进行重投影(GLT校正)

2023-10-22 20:45:19 443

原创 ENVI IDL:如何基于面向对象思想进行编程?

最近打算使用markdown语法进行博客的编写,所以风格和格式方面会有区别,见谅。

2023-10-22 20:32:01 266

原创 ENVI IDL:对于GEOTIFF结构体的说明

其中最关键的只有两个标签Tag,一个是,一个是。至于我没用过不了解,但是应该是关于仿射变换相关的,用于将像素坐标与地理/投影坐标进行转换的矩阵。

2023-10-20 17:47:40 262

原创 ArcGIS: 第二届全国大学生GIS技能大赛(广西师范学院)详解-下午题

某地区拟对X流域综合治理规划设计,需进行流域分区,并进行流域地形特征分析。现有该区域的地形数据,请按照要求进行数据处理,并制作专题图;同时梳理上述问题及其解决过程,规划出科学合理的GIS应用系统,进行系统功能及界面设计。

2023-10-10 00:45:18 1989 1

原创 ArcGIS Engine:鹰眼图的拓展功能-点击和矩形+坐标状态栏

说是拓展,不过是忘记了实验还有附加实验.这里补上.ArcGIS Engine:视图菜单的创建和鹰眼图的实现_炒茄子的博客-CSDN博客这里加上三个功能:为程序增加鹰眼(缩略图)显示功能,并扩展鹰眼功能1)增加鹰眼图控制功能:在鹰眼图中点击鼠标,跳转主图到鼠标点击位置;在鹰眼图中绘制矩形,跳转主图到该矩形范围;2)增加地图控件操作响应鼠标在地图上移动事件,将当前地理坐标显示在状态栏;

2023-10-05 23:43:31 1240

原创 ArcGIS Engine:视图菜单的创建和鹰眼图的实现

因为我们的代码一定意义上是从头开始写的,所以每一个控件都是互不关联的,例如你的ToolBar控件中的工具,正常人都知道这个工具应该是需要在MapControl控件中进行使用(或多或i少),例如TOC控件中展示的图层信息应该是在MapControl控件中显示的那些文件,再例如你使用ToolBar控件中的工具进行数据的加载那么加载好的数据是不是应该显示在MapControl中呢,所以这又是联系。方法加载一个新的MXD(Map Exchange Document,ArcMap的地图文档格式)时,当前的。

2023-10-05 20:06:43 2058

原创 ArcGIS Engine:报错无法嵌入互操作类型“ESRI.ArcGIS.Geometry.EnvelopeClass”。请改用适用的接口。

此错误是由于尝试直接实例化ArcGIS COM组件的某些互操作类引起的。这在.NET Framework 4及更高版本中是不被推荐的。为了解决此问题,你需要确保在工程的引用中将ArcGIS的互操作类型设置为“不嵌入”。

2023-10-04 10:12:04 967

原创 ArcGIS Engine:实现Shp/Mxd数据的加载、图层的简单查询

双击进入的点击事件中(实际上就是点击了该工具会发生什么事情,这一部分事情放在一个函数中让你自由发挥)由于时间问题,实在没有精力一一说明代码的详情,可以多看视频(视频和本博客代码稍有出入)。这样用户才可以选择想要的文件,但是这个不需要我们自己从头去弄,调用人家写好的就好了。传入四个参数(也可以就传入一个字符串,其他都不管等等也可行,看重载方法有没有写)与前面类似,这里不再详细描述,对于其中的。时间精力有限,其他功能不再一一演示。请看后续代码中该函数的编写.

2023-10-04 00:22:19 2047

原创 ENVI IDL:MODIS SWATH产品的点位-像元提取(另附Python代码)

如果是WGS84等常规一些的坐标系(常规就在于你仅仅依据左上角的经纬度和经纬度分辨率、行列号就可以简单得知每一个像元的经纬度,反过来就是你知道一个像元的经纬度那么就知道这个像元的行列号,再类推一下,如果你知道点位的经纬度,那么你可以近似知道该点位距离哪个经纬度最靠近然后一直使用那个像元的经纬度)。如果距离指定点位(实际解决问题时此处应该就是气象站点)的最近像元(而且是有效的)的距离是10°呢,我个人认为这个最近像元的栅格值与该点位进行对应是没有意义的,所以我们还需要进行范围的限制。

2023-09-26 22:32:59 341

原创 ArcGIS Engine:C#基础语法的了解

学了一些语言了,确实发现,语言只是工具,只有你作为初恋的语言值得你花大量时间去深究。而对于其他作为应用或者说只是作为桥梁的语言去学习其他内容的课程,没有必要过分深究语言的基础语法,基本上了解了一些数据类型、简单逻辑运算、算术运算以及简单的判断和循环就够了。至于其他的,要不遇到了再说,要不报错了再深究,不然永远患得患失或者沉迷于语言的深究中而不能将语言真正应用于我们的专业中。

2023-09-18 22:04:52 412

原创 ENVI_IDL: 基础语法详解

string(a[3, 4], format='%.2f'):格式化输出,或许已经发现了,前面的format实际上是fortran语言的格式化输出样式,但是我不太了解该语言所以用不顺手,我更倾向于C语言的格式化的输出(Yeah,IDL支持C语言和Fortran两种的语言的格式化输出)。:这是定义的一个函数(详见下文),其中的arr是位置参数,prefix,format都是关键字参数;: 这是自定义一个二维数组,注意,数组的定义均是使用中括号[],高维亦是如此;:定义整型变量,s可不写;

2023-09-12 23:54:52 595

原创 ArcGIS Engine10.2 Setup 报错

由于 ArcGIS Engine 10.2 需要 Microsoft .NET Framework 3.5 SP1 来正常工作,但它没有在你的计算机上找到该框架。为了解决这个问题,你需要安装 Microsoft .NET Framework 3.5 SP1。通过微软官方下载Microsoft .NET Framework 3.5 SP1并安装;Ctrl + R输入control打开控制面板;通过控制面板下载 .Net3.5框架。时间有限,具体操作略;

2023-09-05 21:26:28 1313

(永不付费)ArcGIS实操数据

主要是关于如何进行离散点数据插值分析(IDW)、栅格数据的重分类、栅格计算器的简单使用、缓冲区分析、掩膜? 这是部分数据的学习要求: 1)巩固学生掌握矢量数据插值分析、栅格数据重分类、叠加分 析的基本原理; 2)熟悉 ArcGis 中离散点数据插值分析的基本方法; 3)熟悉 ArcGis 中栅格数据重分类、栅格计算器的基本操作; 4)熟悉 ArcGis 中栅格数据分区统计的基本方法; 5)了解 ArcGis 中缓冲区分析、按掩膜提取的基本方法。

2022-12-17

专题地图制作相关文件_rar

(注:以下内容均为公开资源,可在国家地理空间数据网站查找并下载) (比例1:400万) 主要包括 省级行政区 省会城市 全国县级统计数据 线状省界 地级城市驻地 国界线 经纬网 中国地州界 中国湖泊 中国县界 主要河流 主要铁路

2022-05-22

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除