学习笔记
GIS小小研究僧
好哈学习,天天向上
欢迎大家关注我的公众号:GISideas
展开
-
地理信息系统教程知识点
10、 投影转换的方式包括正解变换和反解变换,正解是直接建立两个投影间的解析关系式,反解是建立两个投影与地理坐标之间的关系式,先反解出地理坐标,再由地理坐标解出新投影。投影转换的方法包括解析变换(已知新的和原先的投影公示)、数值变换(原先投影公式未知或者难以求得时,利用两个投影间一些离散点的坐标,使用数值逼近的方法建立两个投影间关系式)、解析-数值变换(新的投影公式已知,原先投影公式未知,先反解出地理坐标,再由地理坐标直接代入新的投影公式)。在同一条经线上,维度越低,变形越大,变形最大值位于赤道上;原创 2024-04-28 14:56:32 · 801 阅读 · 0 评论 -
快速获取网页所有图片/获取网页电子资源内的图片
但是今天又发现了一种更快捷方式,是直接使用Ctrl+S,保存到桌面一个文件和一个文件夹,文件夹内有所有的图片,然后直接复制需要的图片就可以了。在之前的文章介绍了使用IDM软件下载所有的图片的方式,这种方式需要获取一个图片的地址并迭代。有时候看一些电子资源/电子教案过程中,想把这些图下载下来,但是不能一个个截图。原创 2024-03-17 21:25:46 · 199 阅读 · 0 评论 -
水土保持与荒漠化防治专业课考试+课本内容整理-西农水保所
只整理了两本书的部分内容,剩下的没有整理完的还包括土壤保持原理那本书的后面几章,水土保持学概论 水土保持工程措施、农业技术措施和林草措施三章。在西北农林科技大学水保所,西农水保所申博考博过程中,整理了课本的内容。感觉还是要有关系,硕士导师托关系找老师或者是有好成果才行。我是得知和本校的竞争一个老师后几乎就想直接放弃了,也就没考上。后面还有可能的真题。原创 2024-03-16 19:23:03 · 543 阅读 · 0 评论 -
Matlab热图制作学习-Zhaoxu Liu / slandarer (2024). special heatmap (https://www.mathworks.com/matlabcentral
Zhaoxu Liu / slandarer (2024). special heatmap (https://www.mathworks.com/matlabcentral/fileexchange/125520-special-heatmap), MATLAB Central File Exchange. 检索来源 2024/1/7.引用也是直接参考的大神提供的现成的文件,略作修改即可。除了修改一些字体字号等,大神提供的SHeatmap.m文件里的内容几乎都不需要修改。还需要使用颜色包工具。原创 2024-01-07 12:43:33 · 486 阅读 · 0 评论 -
arcgis 批量删除Table中的某些Field
当shp或者table文件较少时,可以手动删除每个文件中的某些字段,当文件较多时,就需要使用arcpy或者model进行处理。原创 2023-11-07 16:42:04 · 588 阅读 · 2 评论 -
Arcpy多个栅格数据求平均,使用arcpy.ia.Raster和arcpy.ia.RasterCalculator([RaterA, RasterB], [“x“, “y“], “x+y“)两种方式
这种方式,但是这种方式计算出来的结果加载到arcgis中,Nodata是存在的,但是放大影像后Nodata黑边又不见了。但是两种方式计算的结果值是相同的,只是Nodata的拉伸渲染不一样。有时候需要计算多年的栅格数据平均值,生成一个多年均值栅格,少的时候可以直接使用栅格计算器.进行栅格化后直接相加,而是使用的栅格计算器。在计算平均栅格的时候,一开始不是使用的。但是很多的时候栅格计算器就很麻烦。原创 2023-10-06 12:55:44 · 697 阅读 · 0 评论 -
postgresql,在pgAdmin中修改列名称和列的类型
是否要加引号要看情况,不加引号如果报错,就要加上。无引号:那么所有内容将自动为小写,不区分大小写;带引号:所有内容才能区分大小写。所以这里要带上引号。例如从bigint转为integer。原创 2023-09-26 17:05:32 · 874 阅读 · 0 评论 -
Arcgis pro属性表字段计算生成随机数
【代码】Arcgis pro属性表字段计算生成随机数。原创 2023-09-26 14:21:12 · 313 阅读 · 0 评论 -
SolVES4.1学习2——导入数据运行模型
使用样例数据运行模型很容易,运行自己的数据要根据教程先对数据进行预处理之后根据教程导入数据。首先新建一个solves数据库,之后restore。原创 2023-09-11 14:48:08 · 538 阅读 · 6 评论 -
SolVES4.1学习1——安装与使用教程
SolVES 4版本是QGIS插件,但实际使用过程中发现在最新版的QGIS安装该插件过程中,会报错或异常。因此需安装特定版本的软件。共需安装如下图软件及Java环境等。之后加载上插件,加载过程可能出问题,要卸载插件再重新按一下就没事了。QGIS也要设置Path,添加上Java环境。不知是网络问题还是啥,我的OSM图层没有加载成。根据官方文档安装好后,可以进行相关操作。再点OK,之后可能会卡,运行一会。运行一会后得到另一个界面。运行完后,这个界面会自动关闭。之后根据教程设置参数等。原创 2023-09-04 17:53:10 · 698 阅读 · 2 评论 -
QGIS学习3 - 安装与管理插件
QGIS安装与管理插件主要是使用了菜单栏安装与管理插件这个菜单。原创 2023-09-04 09:19:44 · 1383 阅读 · 0 评论 -
ArcGIS将两个相同范围但不同比例或位置的矢量数据移动到相同位置
有两个市图层,一个是正确经纬度的市行政范围图层,另一个是其他软件导出获取的不正确经纬度信息或缺失信息。如果单纯的依靠移动图层,使不正确的移动到正确位置需要很久。尝试定义投影等也不能解决。使用ArcMap 的空间校正工具条,进行矢量数据的配准。就将两个图层大致校正到一起了,再设置相同的投影就可以了。开始编辑后,激活工具条,使用建立空间连接工具。建立足够连接后,查看误差。原创 2023-09-01 18:43:07 · 1257 阅读 · 0 评论 -
Arcpy 对栅格数据批量归一化
这里直接使用的代数就可以相除,自己还是搞不太懂,什么时候要使用arcpy.ia.raster,arcpy.ia.rastercalculator,还是要多实验一下。这里的InData \MinValue和MaxValue都是栅格数据,但是。,产生的MinValue和MaxValue是临时存储到哪里了呢?有一批栅格数据,最大最小值不一样,但是要进行归一化。原创 2023-07-04 18:32:40 · 466 阅读 · 0 评论 -
使用Python Pandas批量删除EXCEL中的特定列
【代码】使用Python Pandas批量删除EXCEL中的特定列。原创 2023-06-29 09:02:03 · 1491 阅读 · 0 评论 -
Python+ArcPy多个矢量文件批量掩膜多个栅格数据
在使用模型构建器等进行掩膜提取时,有时候会有要使用多个矢量文件去掩膜裁剪多个栅格数据的情况,模型构建器内只能进行一次迭代,无法直接进行多次迭代,在模型里面嵌套模型多次迭代还不如直接使用arcpy方便。下面我有20来个矢量文件,要分别掩膜2000-2021年的栅格文件,基于arcpy。原创 2023-05-02 19:59:51 · 703 阅读 · 5 评论 -
Python3使用sys.argv和os.system 从一个程序调用另一个程序,并将参数传递
由于实验需要,需要从A.py 调用另一个B.py,并将A.py中的参数mean、max、min三个值传递给B。这里参考了其他人的文章。原创 2023-04-20 12:27:27 · 893 阅读 · 0 评论 -
The Copernicus Global Land Service (CGLS)账号注册与数据下载
全球土地服务处系统地生产了一系列关于地表状况和演变的合格生物地球物理产品,这些产品在全球范围内以及在空间分辨率中低的情况下生产,并辅之以长期时间序列的构成。2、如果只需要某些固定范围或时间的数据,直接下载就好了,如果批量数据下载需要在线提交订单,订单多久能呈现complete状态感觉不是很准确,有可能快有可能慢。在订单提交后会获得一个邮件,要从邮件中下载,直接使用网页不一定能获取下载链接,我是没有找到。1、账号注册:注册账号时不能使用国内的网络,注册后的浏览与下载等可以使用国内的的网络。原创 2022-11-16 20:33:31 · 1854 阅读 · 4 评论 -
ArcGIS中的Excel to table等工具无法输出结果
ArcGIS中的Excel to table等和Excel转换相关的程序有时候是无法正常允许,其中很大的原因是没有安装 Microsoft access database engine,再有一部分原因是安装后但是没有启动等,只有安装后的第一次正常使用了,但是关机重启后,engine就自动关闭了。再使用时就需要重新修复或者安装一下才是开启中的,不知道怎么去设置?原创 2022-11-04 09:03:36 · 2082 阅读 · 0 评论 -
使用IDM从网站下载图片——从翻页网站下载图片
有时候一些政府或者事业单位的报告等是翻页的图片,还经常很多页,但是我们需要下载下来看我们可以批量下载一下。原创 2022-10-30 20:31:51 · 2072 阅读 · 0 评论 -
气象站点海拔数据异常情况及处理
首先是将arcgis ele字段乘以了10,因为我前面多除了10,之后减了10000,去掉了异常。从国家气象数据中心获取的部分数据,在使用matlab代码处理过程中,以为海拔高度是厘米,所以都除以了100,后来查说明才发现是0.1m,只需要除以10。比如55294这个站点,海拔高度14800m,明显是错了,在与DEM数据进行对比后发现,是这些数据都大了10000。但是在看arcgis 属性表时发现部分站点如果乘以10,海拔一万多米,明显错误。原创 2022-10-29 11:08:43 · 927 阅读 · 2 评论 -
使用Anusplin对降雨侵蚀等进行气候插值
在以前的论文中,对气候插值使用了kriging和cokriging,但是效果都不是很好,以前也使用过Anuspin,但是因过程太复杂,当时也没有做记录就不想再用,但是无奈其他软件的插值效果实在不好。原创 2022-09-26 17:19:25 · 1812 阅读 · 21 评论 -
ArcGIS栅格数据边缘不匹配,栅格不对齐
很多时候很多分析,并不需要栅格数据分辨率相同,也不需要边缘匹配到一起,但是有的时候边缘不匹配,分析就会报错。在处理时,就遇到了这种情况,看到有人使用clip进行匹配我尝试使用采样的方式进行匹配。原创 2022-09-21 09:17:29 · 3593 阅读 · 0 评论 -
ArcGIS Pro 3.0学习版安装教程
首先复制dll文件到安装路径下,默认路径是:C:\Program Files\ArcGIS\Pro\bin。最后复制liensing文件夹到:C:\Users\ 用户名 \AppData\Local。首先会解压安装包文件 ,之后执行程序。这里解压的安装包文件,最后解压文件可以删掉。包括中文文件包、数据转换拓展模块等等,可以用得到的时候再安装。下载后直接管理员安装即可。原创 2022-09-19 12:42:38 · 6001 阅读 · 6 评论 -
ArcGIS Pro 模型构建器修改输出文件的拓展名
首先是有一部分矢量文件,后缀为.shp。要使用这些矢量文件进行掩膜提取,将后缀修改为.tif。最笨的方法是直接在原.shp后面加上.tif,构成"name".tif,如果需要将.shp删除并修改为.tif,就需要使用parse path。首先导入各省shp的文件夹,之后迭代要素类,将要素类的name输出连接到parse path,作为parse path的输入,之后将name(2)作为掩膜提取的前提条件,设置好栅格和要素类即可。输出结果为各省名称命名的tif文件。原创 2022-09-19 08:57:36 · 1620 阅读 · 5 评论 -
Arcgis处理栅格数据中Nodata值
在前面介绍了对栅格数据Nodata进行补充的一种方法,在实际研究过程中,面临的情况多种多样。下图为LUCC数据,在各省边界上存在较多的数据缺失值,这个缺失值一方面可以使用更高精度的其他数据去填充。在缺失值较少的情况下,也可以使用邻域的数据去填充。以下还有几种方法对数据填充,但都是使用Raster Calculation 工具。原创 2022-09-16 16:05:04 · 3374 阅读 · 2 评论 -
arcgis 栅格数据热点分析
接下来要使用的是:join field工具,注意:直接使用join链接的是临时表,也可以按照官网的建议join后复制要素,进行永久性保存。对每一个栅格进行分析。这个需要将数据先转为nc文件,当然我这个nc文件只有一层。要实现的目的:对栅格数据按照区县进行分区统计,之后对各区县的情况进行热点分析。我这里没有对空间关系的概念化和距离方法进行更改。首先使用分区统计到表,计算各县的平均值。时空模式挖掘是针对每一个栅格进行。,里面有详细的介绍和案例。原创 2022-09-07 08:54:11 · 4535 阅读 · 2 评论 -
Matlab 删除含有特定数值的行 / Matlab定义excel表头
1、删除含有特定值的行,以含有-9999为例。原创 2022-09-06 16:01:33 · 3318 阅读 · 0 评论 -
ArcGIS申请添加天地图服务
之后输入服务与Key,服务只需要复制wmts前面的链接就行。value是申请的key。后来查看别人的介绍才知道,要在wmts前面加上esri才行。如果没有账号需要先注册一个账号,注册账号后,申请Key。这样在线服务和本地的shp文件就能对应起来了。结果加载后和我的矢量数据有很大的偏移。点击创建新应用,并进行相关的设置。之后选择一个自己需要的地图服务。...原创 2022-08-01 14:18:30 · 1766 阅读 · 0 评论 -
【ArcGIS Pro】在不影响其他Map Frame中Map情况下修改Map,断开同一个Map在不同地图框中的关联
在使用ArcGISPro出图过程中,出现一个很无奈的问题如下图所示,对某一个地图框中的地图做出禁用或者修改,另一个地图框中的地图也会做出修改。这个问题是arcmap中不存在的。arcmap中默认的是断开了不同的map中的连接。通过查找官方文档等等,发现这个问题已经存在两三年了,但是官方一直也没有做出修改。...原创 2022-07-20 17:34:20 · 810 阅读 · 3 评论 -
Origin正版申请
使用学校提供的教育邮箱或通过联系软件官方:[email protected],申请免费的origin软件使用同时Origin官方提供了帮助与安装教程:https://shimo.im/docs/e9iqj6EdlXYgoZPB首先访问链接,填写相关信息,使用学校提供的edu教育邮箱https://www.originlab.com/OriginProLearning.aspx官方提供的学习版,许可期限6个月,许可过期后可重新申请许可如果没有教育邮箱,也可以邮件联系一下软件官方。接下来就会收到验证邮原创 2022-07-11 11:02:53 · 2420 阅读 · 0 评论 -
土地利用转移矩阵简易计算方法
土地利用转移矩阵有很多计算方法,如果使用的是二级地类,地类多或数据量大,可以参考这个方法:https://zhuanlan.zhihu.com/p/373096177如果是一级地类,可以直接简单计算。如图所示,一共六种地类打开栅格计算器,将第一年的土地利用10+第二年的当然,这只是一种思路,一种转移的表示方式,乘以100等等都行,甚至转成矢量数据,进行字段运算也能产生想要的结果。打开属性表,11便表示从1这种地类到1这种地类的面积,12表示从地类1到地类2的面积。如果是年份很多,可以将图层转原创 2022-07-11 10:01:12 · 2341 阅读 · 1 评论 -
Matlab 统计每行或每列中,元素出现的个数及比例
使用地理探测器进行探测分析,有X变量是单一值,因为有大量的矩阵,不能手动单一修改,所以要先对矩阵按列,统计每列各个元素的个数,若有一列只有一个或两个值,则要删除该列。tabulate()函数是一个常用的计算矩阵中各个元素出现次数的函数,这里我们只需要稍加变通,设全部的矩阵为dataall,只需要在前面加上一句:data = dataall(:,1),类似于这样的一句,这样table = tabulate(data)计算的就是某一列各元素出现的个数。再通过if函数判断,如果table行数少于3,那么便是不符原创 2022-07-05 20:45:21 · 6025 阅读 · 0 评论 -
R语言使用substr函数、paste函数或str_c函数提取并生成输出文件名
有时候使用R语言处理后的结果需要按照原先文件的名称进行输出保存,可以使用substr函数从原先文件中提取部分需要的信息,使用paste函数或str_c函数拼接输出的路径与输出的文件名、文件格式等。因为自己只是在使用中,学习了这三个函数,没有系统的学习过R语言,所以只能简单的介绍一下。首先,使用这三个函数,需要安装包substr函数可以从一段字符中,截取自己需要的字符。以我自己的使用过程为例:我的原先文件名是一堆城市名称,我的城市名称是:xxxx.xlsx,我要提取的只是城市名称,所以我的代码是:,城市名称原创 2022-07-05 16:37:49 · 2244 阅读 · 0 评论 -
Arcgis pro通过渔网工具生成规则采样点,并对栅格数据进行采样
当有一组栅格数据,需要对栅格数据进行相关性分析等操作时,可以通过采样的方法,将栅格数据的信息提取到EXCEL中,在R或MATLAB中进行数据的相关性分析等。推荐使用Pro,而不是arcmap进行,在采样的能力上,Pro比arcmap快不止几倍。使用创建渔网工具进行采样点的生成,当然也可以使用生成随机点的方式来进行采样,但是毕竟采样点越多,结果与原栅格图层越一致。首先输出要素类需要设置输出的路径与名称渔网的起始坐标点需要设置,这里可以手动设置,也可以使用下拉菜单进行设置推荐使用下拉菜单,会默认填充采样的原创 2022-07-01 20:31:58 · 3828 阅读 · 5 评论 -
基于arcgis 计算RUSLE(修正通用土壤流失方程)的插件GISUS_M介绍
GISUS-M是一个arcgis 插件,通过这个插件可以计算rusle中的五个因子插件的链接:https://www2.ufrb.edu.br/gisus-m/downloadm这个插件在使用过程中,有几个问题:1、LS因子是使用的西北农林张宏鸣老师等人开发的LS_TOOL工具,这个工具对于计算的数据量是有限制的,而且对于计算机的硬件也有一定的要求,如果是大范围高精度的DEM很可能计算不出来结果。2、对于C因子采用的是NDVI进行的计算,但是国内很多研究是采用的植被覆盖度估算C因子。在我个人的使用原创 2022-06-28 13:09:52 · 3334 阅读 · 14 评论 -
arcgis使用栅格计算器con与ISNULL函数,对栅格数据中缺失值/空值进行设置
使用isnull函数,可以将栅格数据中的空值进行判断,使用con条件函数,可以对空值进行设置。isnull函数逐个像元来确定输入栅格中哪些值为 NoData。如果输入值为 NoData,则返回 1,否则返回 0。现在有一个自己计算的栅格数据,因为水面等,在栅格数据中产生了空值。需要在这些空值处补上数据同时又通过其他途径获取了别人计算的数据,通过采样后对数据进行分析,发现数据相关性很高,但其他人计算的数据范围较小,幸运的是正好覆盖了自己的空值区域,故将其他人的计算结果填充到这个空值区域。公式:将K8图层原创 2022-06-28 10:03:29 · 7489 阅读 · 2 评论 -
arcgis pro按照属性字段分割矢量图层的每个要素
在这里主要是使用按属性分割工具(Split by attribute )以行政区划数据为例:设置好输出的路径与用于分割的字段:按照省份名称自动生成了新的shp文件:原创 2022-06-24 17:14:12 · 1329 阅读 · 0 评论 -
HWSD数据处理
数据下载于国家青藏高原科学数据中心HWSD.mdb数据库中的数据可以直接使用excel导入,并进行查看。其中中国区域的数据来源自南京土壤所提供的第二次全国土壤调查数据。核心土壤制度单元唯一验证标识符:MU_GLOBAL-HWSD数据库土壤制图单元标示符,连接了GIS图层。MU_SOURCE1 和 MU_SOURCE2- 源数据库制图单元标识符SEQ-土壤制图单元组成中的土壤单元序列;土壤分类系统利用FAO-7分类系统或 FAO-90分类系统(SU_SYM74 resp. SU_SYM90)或FAO原创 2022-06-23 16:14:51 · 2452 阅读 · 6 评论 -
GIS与空间分析
地理信息系统主要研究在应用计算机技术对地理信息进行获取、存储、管理以及处理和分析过程中所提出的一系列基本问题。具体有:1.地理信息信源理论:地理空间认知,地理信息本体,自然信息源和社会信息源本质等。2.地理空间信息获取与传输:时空信息基准,地理信息传输机理与増衰规律、地理空间信息不确定性等。3.地理信息存储与分析处理:地理信息标准、压缩与检索、变换与处理等。4.地理空间模式信息处理:地理信息分类、空间内插、模式识别等。5.地理信息知识处理:地理知识的表达与获取、空间推理与建模、知识信息处理系统等。原创 2022-06-08 17:03:28 · 5314 阅读 · 2 评论 -
遥感信息分析与应用期末复习题
名词解释:thermal inertia,是一种综合指标,是物质对温度变化热反应的一种度量,即度量物质热惰性(阻止物理温度变化)大小的物理量。高热惯量的物质对温度的变化阻力较大,常用P表示,公式为P=[Kpc]1/2,k是热导率,c是比热容,p是密度。热惯量主要应用于土壤表层水分变化监测、干旱遥感监测等,例如:将土壤水分含量与土壤热惯量建立统计关系,并将遥感卫星的数据与热惯量联系起来,通过遥感数据获得真实热惯量,进一步反演土壤水分含量分布。:Mixed pixel decomposition,混合像元是指原创 2022-06-08 16:54:21 · 2881 阅读 · 0 评论