- 博客(177)
- 资源 (1)
- 收藏
- 关注
转载 构建栅格属性表
有时候会遇到将32位栅格数据提取属性表的操作,但是一般此类数据都是浮点型,是无法计算得到属性表的。因此我们可以利用数据管理工具下的:复制栅格工具,在最下面选择16位即可,看自己数据情况选择signed或者unsigned类型。导出后的数据可以计算属性表,两种方式,一种是利用工具进行;另外一种是右键属性后,找到颜色条显示,默认是拉伸,选择唯一值显示,即可计算属性表,关闭后即可看到属性表。*****需要注意的是,你直接转出时,会将你原先的浮点型数据四舍五入处理,因此,如果你需要保留2..
2021-03-16 18:54:15 4125
原创 envi拉伸(归一化)功能&ArcGIS栅格计算器&工具箱的模糊隶属度
(1)ENVI可以通过拉伸实现归一化操作。如图设置最小值和最大值就可以拉伸至该区间。(2)ArcGIS栅格计算器&工具箱的模糊隶属度也可以实现归一化处理。
2021-03-13 22:34:19 6290 2
原创 shp文件导入mike zero文件需注意,要将M Z值去除
shp文件导入mikezero文件需注意,要将M Z值去除:Google earth中绘制的kml或kmz文件转的shp文件会存在m z值,这表示该shp文件是三维球形文件而arcgis中的shp文件是二维文件,没有高度信息Mike 11支持arcgis的二维shp文件导入因此,将来自Googleearth的shp文件导入Mike时,需要通过以下方式去除M Z值,即可。...
2021-03-11 10:51:39 2334
转载 ArcGIS python计算长时间序列多个栅格数据的平均值
import arcpyfrom arcpy.sa import *arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput=1 #customarcpy.env.workspace="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\1Length\\"#customoutpath="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\2mea.
2021-03-08 21:16:17 1789 2
原创 arcgis上下标问题
在使用ArcGIS软件出图时,经常遇到图例中计量单位上下标的问题,然而ArcMap中没有像Word方便的进行上下标设定,使用一个小小的函数便可轻松搞定,如下上标:a下标:a在ArcMap中首先进行如下设定:然后插入Legen,就会显示出上标了,...
2021-03-06 23:34:33 674
原创 ArcGIS栅格数据图层空间相关性分析方法
ArcGIS栅格数据图层空间相关性分析方法矢量的空间相关分析可以通过Spatial Join以后,导出属性表到EXCEL中,进行相关分析Correlation Analysis. 而两个栅格数据图层直接能不能直接在ArcMap中进行相关分析呢?答案是肯定的。在Toolbox工具箱里面找到Spatial Analyst Tools->Multivariate,选择band collection statistics工具,选中进行相关分析的两个栅格图层,并且勾选需...
2021-03-05 14:59:07 20969 5
原创 ArcGIS批量栅格影像归一化操作
代码实现arcgis栅格数据归一化:import arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.CheckOutExtension("Spatial") #这个是工作空间目录,待处理的所有影像都在这个文件型数据库中env.workspace = "E:/poet/clip/clip.gdb" #设置读取的栅格类型,可以参考ArcGIS的官方文档,下面的是读取目录下所有类型的栅格文件rasterList ...
2021-03-04 22:50:34 3204 2
转载 如何确定Z检验的值(查正态分布表时要注意中间的数字都是面积,最左边一列和最上面一行都是Z值)
一、当给定了检验的显著性水平a=0.05时,如果检验时要检验是否相等,就是双侧检验,允许左右各有误差,即a/2=0.025。此时要查尾部面积是0.025时的Z值。但是我们参考书中说明表中间的数字是指从最左面一直到右侧某一点的面积,而Z值是指从中间均值所在的位置往右计算的长度。所以当Z=0时,中间的面积=0.50就是这个道理。现在我们要的是从右边尾部面积查Z值。当右边尾部面积是0.025时,左边的面积应是1-0.025=0.975。所以我们查表时要在表中间找到0.975。从这一行水平往左得到1.9,往上对得到
2021-01-19 15:42:53 16263
原创 Mann–Kendall test-Sen’s slope-Quantile regression (quantreg.m)
By “Geographically heterogeneous temporal trends of extreme precipitation in Wisconsin, USA during 1950–2006”(1950-2006年美国威斯康辛州极端降水的时空变化趋势)论文中关于mk趋势检验等的代码:
2021-01-15 17:30:38 499 1
转载 关于显著性检验,你想要的都在这儿了!!(基础篇)
这里有必要提一下anova1函数中的参数displayopt 的作用。在大规模的anova1调用中(例如把anova1放在for循环中反复调用),需要把displayopt设置为'off',否则anova1每调用一次就会绘制两幅图,这样会迅速的耗费计算机的内存,容易造成程序崩溃。除了上文中介绍的第一种调用anova1的方式,还有一种方式用于均衡的方差分析。所谓均衡就是要求不同的组别内的统计数据个数必须相同。在上例中出现的各个组的统计个数分别为{8,6,6}就属于非均衡。在均衡状态...
2021-01-15 15:36:25 2546
原创 基于matlab的长时间栅格序列逐像元多元回归
在长时间序列的栅格数据分析过程中,有时候需要建立多元回归分析来分析多种要素对单个要素的影响。本文以降水和温度为x变量,以NDVI为y变量,逐像元建立NDVI与降水和温度的关系。具体代码如下% @author yinlichang3064@163.com[aa,R]=geotiffread('D:\qixiang\年尺度数据\插值的结果\平均温度\tem2000.tif');%先导入纬度数据info=geotiffinfo('D:\qixiang\年尺度数据\插值的结果\平均温度\tem2000.tif
2021-01-06 16:56:15 6330 19
原创 ArcGIS皮尔逊相关系数
http://blog.sina.com.cn/s/blog_4ba16d110101d8wt.html#cmt_54509B26-7F000001-E29D94AE-8F2-8A0
2021-01-04 18:39:16 5043
原创 Excel趋势曲线中R平方值的含义
R平方值是趋势线拟合程度的指标,它的数值大小可以反映趋势线的估计值与对应的实际数据之间的拟合程度,拟合程度越高,趋势线的可靠性就越高。R平方值是取值范围在0~1之间的数值,当趋势线的 R 平方值等于 1 或接近 1 时,其可靠性最高,反之则可靠性较低。R平方值也称为决定系数。by:http://blog.sina.com.cn/s/blog_671c54fe01012lzs.html...
2021-01-04 12:14:46 21714
原创 arcgis fishnet渔网功能学习
听说要来跟我请教fishnet功能,吓我一跳,arcgis不就只有个创建渔网功能吗??纳尼?!赶紧捡捡知识点,避免一问三不知=.=。版本12.4后工具位于data management/sampling/create fishnet,12.4之前应该是在data management/要素类/创建渔网。放弃了。需要的时候我直接看帮助手册吧。...
2021-01-04 12:14:23 3413 1
原创 基于matlab 的长时间栅格数据的Sen+MK显著性检验趋势分析
在前一篇文章中讲述了用sen法进行长时间的趋势分析,但并未对结果进行显著性检验,通常Sen与MK检验是结合在一起的,因此本文主要讲述如何进行MK检验。具体代码如下% @author yinlichang3064@163.comclear[a,R]=geotiffread('D:\GIS\vegetation\output\yearmax\1982.tif'); %先导入投影信息info=geotiffinfo('D:\GIS\vegetation\output\yearmax\1982.tif
2021-01-03 11:38:34 3158 9
转载 学习笔记之 基于matlab的GIMMS3G-NDVI数据处理-从nc文件解析到月和年最大合成tif
GIMMS的NDVI数据常常用来进行长时间序列的植被情况分析,最新的NDVI3g数据集,时间跨度为1982-2015年,半月合成,空间分辨率大致为8km(0.08333度),下载地址https://nex.nasa.gov/nex由于下载下来的是nc文件,部分用户可能不熟悉该数据格式,无法很好的使用该数据集,本文提供一套基于matlab的处理方式。% @author yinlichang3064@163.comlon=ncread('F:\Global\NDVI3g\ndvi3g_geo_v1_20
2021-01-03 11:35:33 2599 6
转载 学习笔记之 基于Matlab的NDVI最大合成
采用最大合成法将月数据ndvi合成年数据NDVI,本文以spot的月ndvi为例进行,利用matlab平台,代码如下所示:%author yinlichang3064@163.com [a,R]=geotiffread('F:\项目\data\spot_月ndvi\ndvi200001.tif');%先投影信息info=geotiffinfo('F:\项目\data\spot_月ndvi\ndvi200001.tif');for i=2000:2015datasum=zeros(size(a,
2021-01-03 11:24:06 2200 2
转载 多波段栅格数据在arcgis中的显示
还有一种情况是这样,我们在拿到数据方提供的卫星影像时,是7个分开存储的单波段的影像,在ArcMap中全部加载进来也还是灰色的。如何将多个单独的灰色影像以真彩色的方式渲染出来呢?首先要注意的一点是,与专业的遥感软件不同,ArcGIS中不提供直接在一个栅格上叠加其他栅格来进行渲染,也就是说不能简单在Symbology中设置一下来实现RGB的渲染,只能通过事先将不同波段的影像数据组合在一起,然后在进行加载。另外,由于不能在指定的波段上选择Load Data,所以也不能通过导入数据库的方式将不同波段的.
2020-12-29 20:32:04 5714
原创 平均气温
在气象学上,平均温度是一个专有的名词。它分年平均气温月平均气温和日平均气温三种。我们可以按定义日平均气温和月平均气温的方法来定义年平均气温。既然日平均气温指一天中最高气温和最低气温的几何平均值。月平均气温是指一个月中最高温度与最低温度的几何平均值。那么年平均气温就是指一年最高温度与最低温度的几何平均值。...
2020-12-28 17:02:28 860
原创 ArcGIS 栅格图层相关性分析
Spatial Analyst Tools——Multivariate(多元分析)——Band Collection Statistics(波段集统计)。添加图层,勾选Compute covariance and correlation matrices以输出相关第分析结果,结果保存成txt。使用的是皮尔逊相关系数(Pearson Correlation Coefficient)。by:http://blog.sina.com.cn/s/blog_6215b1010102xvt7.html.
2020-12-28 16:58:24 9550
原创 你为什么要选R?
关于学术绘图。具有工科背景的人员常使用matlab,具有计算机背景的人员常使用Python,具有统计学科背景的人员常使用R。????软件绘图便捷,但灵活性略显不足,只能去适应其操作规则,而不能够创造新的图表。matlab、Python和R,则需要编程才能实现图表的绘制。灵活性大,自足设计图表。三款选啥不要紧。所谓“心中有剑,落叶飞花,皆是兵器”。by摘抄...
2020-12-25 22:07:35 160
原创 深入学习投影转换那点事儿(1)
地理坐标系=大地坐标系地理坐标系:geographic coordinate systems投影坐标系:projected coordinate systems空间参考:空间参考大于坐标系,包括坐标系、空间域、精度和容差等。数据框坐标系和数据本身坐标系。动态坐标转换(动态投影) :将不同坐标系的数据添加到地图中,ArcMap会自动进行坐标转换。坐标转换(生成新的数据):①raster-project;②create custom geographic transformation....
2020-12-17 15:57:24 141
原创 R语言实现NetCDF气象数据的逐月解析提取
NetCDF格式数据在气象领域应用较多,处理此类数据的思路或软件有:ArcGIS/QGIS/R语言。以全球降水量数据集“cru_ts4.04.2001.2010.pre.dat.nc”为例。ArcGIS:利用“多维工具”进行导入,在“属性-NetCDF”一栏选取对应月份,逐一导出。QGIS:查看方便,导出tiff格式后数据容易丢失。R语言:方便/快捷/数据无丢失。代码如下:#打开nc文件,查看nc文件数据结构library(ncdf4)ncdata <- nc_op.
2020-12-16 12:25:31 4870 12
原创 QGIS导出NC月份数据
1.下载/安装QGIS,加载NC数据。2.属性-符号化,可以看到“渲染类型”选择。3.选择好渲染类型,即可进行影像逐月导出。
2020-12-15 20:10:11 1442
原创 世界数据图谱分析平台 www.knoema.com
「世界数据图谱分析平台是一个收集整理世界多方面数据而后进行统计分析并制作成简单易懂的图谱报告的平台,该平台世界各国的关键统计数据可进行分类浏览,搜索相关数据并可视化,可将其导出为Excel,PowerPoint和PDF等。特点:1、全世界190多个国家详细统计的关键数据;2、 世界各国的区域统计数据,如美国,英国,加拿大等- 经济,人口,医疗,犯罪,能源,农业,教育和贸易指标;3、国家排名- 交互式的世界地图- 国家和地区间的并行对比- 相关信息的可视化。Knoema is a knowledge pla
2020-12-14 15:27:05 1164 3
原创 五个公开分布在美国的统计缩小的气候数据集08(待完善)
五个公开分布在美国的统计缩小的气候数据集:ClimateNA,NASA NEX-DCP30,MACAv2-METDATA,MACAv2LIVNEH,WorldClim。ClimateNA (https://adaptwest.databasin.org/pages/adaptwest-climatena/),美国国家航空航天局NEX-DCP30 (https://cds.nccs.nasa.gov/nex/),MACAv2-METDATA (https://climat.
2020-12-08 22:09:48 181
原创 什么是R&D?
https://wk.baidu.com/view/a71f320852ea551810a68702?pcf=26、 对应届生来说,考研不能解决你的问题,对90%的人来说,考研大概率就是得到一个心安理得继续浪费几年时间的理由而已。这一点在2008年已经被你们的前辈验证过了。7、 对那10%能走科研路线的考研者不适用以上结论,有机会尽量向科研路线上靠拢。现在低技术含量的消费互联网已经进入尾声,下个风口会分散在各行各业的企业服务里,从头部创投今年新投的项目可以窥见一二。R&D可能是你突破简单重复劳动
2020-12-06 22:57:24 321
原创 解读一段读取NC文件的R代码
library(raster)petbr <- brick("TerraClimate_pet_2003.nc", varname="pet")writeRaster(petbr, filename =petbr@data@names, bylayer=T, format="GTiff")以上:1)<-位“=”,赋值;2)若提示存在“巴拉巴拉”,但没有这种溢出号。是因为R中索引用/不能用\,用\就会提示这玩意。...
2020-12-05 21:55:42 2273 1
原创 XML文件是干什么用的?
答:xml最大的用处就是存储数据,因此可以说,xml就是纯文本型数据库!可以跨平台、跨系统的使用。一个好的Xml文件可以作为各种软件、Web页面和程序的数据库。与html的区别在于:html的用处是怎么样显示数据,而xml的用处就是存储数据并供软件和页面使用。用记事本和浏览器(最好是IE5.5或以上版本)都可以打开它。更多见帖子:https://segmentfault.com/a/1190000013252686...
2020-12-05 11:29:47 6835
原创 originpro实现纵坐标断点
https://jingyan.baidu.com/article/fdffd1f843ae63f3e88ca150.html
2020-11-16 23:32:25 2792
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人