
地信GIS
忠言睿长
从事3S技术与水文气象研究
展开
-
ArcScene三维显示DEM
ArcScene三维显示DEM1.打开ArcScene软件2.加载DEM数据3.原创 2020-11-13 10:56:16 · 6016 阅读 · 3 评论 -
win10安装arcgis10.2问题
1.开始安装时,提示必须安装.net framework 3.5以上 提示框内容:ArcGIS 10.2 for Desktop requires Microsoft .NET Framework 3.5 sp1 or equivalent. Please see the installation Guide for more information. To download and install Microsoft .Net Framework 3.5 sp1, visit the Microsof原创 2020-09-24 11:44:56 · 4786 阅读 · 0 评论 -
arcpython设置处理边界extent
原创 2016-08-02 16:13:23 · 1914 阅读 · 0 评论 -
Arcgis+Python实现对栅格归一化处理
图像归一化就不多说了,就是(数值-min)/(max-min),把结果都划归到0-1范围,便于不同变量之间的比较,取消了不同数量差别。第一个方法,需要对栅格数据预先知道取值范围。第二种方法,比较好点,直接读取属性原创 2016-01-15 10:25:53 · 12465 阅读 · 12 评论 -
MATLAB计算宇宙外太阳辐射Ra并存tif格式
原创 2015-10-12 22:37:36 · 1907 阅读 · 1 评论 -
ArcGIS中制DEM渲染图
1、DEM计算山体阴影,注意在Hillshade对话框中,zfactor设置为5,即夸张效果。2、在DEM图层的Property-Symbology中选择下面所示两者中的一个,因为这两个色彩系比较适合以DEM制图效果,同时Stretch的n设置为5,也是夸张效果,注意这里有个Use hillshade effect效果不明显,可以不选择3、将DEM图层放置到上面,将Hillsh转载 2014-04-14 16:55:13 · 5245 阅读 · 0 评论 -
ArcGIS栅格数据SetNull IsNull Con函数
当两个栅格进行叠加,有时会有一部分没有数据,即用identify点击该区域,Value为NoData,而不是像其他非空区域一样有值。此时注意nodata区域要赋予0值,因为nodata任何数=nodata,因此当需要对栅格空值数据赋值时,要采用条件查询函数将NoData的地方进行赋值。(1)方法是ArcTools->Spatial AnalystTools->Map Algebra->Si转载 2014-12-14 20:59:12 · 52041 阅读 · 2 评论 -
Arcgis重分类无法赋新值为小数
问题如题目所示,无论待分类栅格数据时int还是float类型重分类都不能赋值小数,这是arcgis本人问题。解决方法是对你要分类的所有数值乘上一个常数,如100,1000,保证所有值都为整数了,分类成功后再使用栅格计算器除以你的乘数。如我要分类为0.01, 0.5,1, 2四类,那么我在分类器raclassifiy中写成1,50,100,200,分类完成后再除以100.0(加上小数保证得到的是原创 2015-07-13 11:09:35 · 10224 阅读 · 4 评论 -
Arcgis实现为栅格中nodata或者数据质量差的点赋邻域值Nibble工具
工作中遇到一个栅格图中有一些区域或者几个别的点为nodata值,或者已知数据中某个点的值质量较差想用邻近值替换掉。这时可以使用nibble工具填充1、首先使用栅格计算器中的con和isnull函数组合,为原始图像A中nodata赋一个值(否则nodata不参与运算),如9999或其他,得到图像B。2、在nibble工具中Input raster中输入B图像作为待处理的图像,inpu原创 2015-07-13 11:26:08 · 10190 阅读 · 1 评论 -
在ArcGIS中利用“行政单元面积权重法”实现人口数据格网化
1.在ArcGIS中,基于行政区划图建立一张Fishnet(注:格网大小应根据数据尺度确定,本例为1KM*1KM),具体方法见前一篇文章。2.将该Fishnet(为Polyline类型)转换为Polygon类型(Data Management Tools-Features-Feature to polygon)。3.将转换后的Fishnet按照行政区划裁剪(注:可利用Sel转载 2015-06-04 15:28:19 · 11173 阅读 · 2 评论 -
Python常用常见知识汇总
1、获取一个文件的子串如: str1="myfirst.tif",想获取不带扩展名的文件名最简单方式为str2=str1[0:-4]2、获取指定文件夹下的所有指定类型文件arcpy.env.workspace=inputDirfileList=arcpy.ListFiles("*.tif")即可获取所有tif文件的链表3、获取指定文件路径下的所有文件名称d原创 2015-06-16 23:30:07 · 739 阅读 · 0 评论 -
MATLAB读取nc文件写成TIF文件arcgis统计错误 -不能正常显示统计信息,空值处理
问题描述:MATLAB读取nc文件很方便,然而arcgis或者ENVI等地理软件往往是GeoTif。本文在MATLAB读取nc文件并写成tif后发现能在arcgis中正常显示,黑色区域显示的为NoData,感觉一切正常。但是统计信息中的mean和stdeve显示是-1.int, 如下图所示:有人说,既然arcgis显示没问题,那就行了呗。不然,如果统计信息有误,那么你在进行区域统计的时候就原创 2016-07-15 10:26:50 · 6365 阅读 · 4 评论 -
MATLAB海量数据读取access数据到excel并保存为shapefile文件
(1) 通过access数据库通过联立查询并生成新表由于access只限制5万条记录,所以要通过MATLAB读取所有数据并保存到xlsx中去。SELECT * INTO AA IN 'E:\aa\aa.accdb'FROM (SELECT MetroSta742Info(modify).*, MetroStaDataPET_2013.*FROM MetroStaDat原创 2016-07-31 20:10:58 · 4390 阅读 · 0 评论 -
arggis大批量导入坐标点报错
经常在野外或者获取的数据经纬度坐标,需要导入到arcgis中,以点的形式显示,那么就需要通过arcgis的"add X Y"功能,通过制定经度, 纬度等对应的变量,把带坐标的文件导入。那么需要注意以下事项:1.arcgis中导入带坐标的文件不可以是 .xlsx格式,可以使.xls格式(即1999-2003格式)或者.csv格式。所以,需要先把 .xlsx格式转换为 .xls格式。但是问题是.x...原创 2019-04-14 23:21:00 · 1847 阅读 · 0 评论 -
arcPython细节汇总-setll函数和地图计算Raster Calculator函数
1、 setnull函数 arcpy.CheckOutExtension("Spatial") outSetNull = SetNull(inputRaster,inputRaster,"VALUE < 0") outSetNull.save(outputRaster)注意:一定要检查许可否则会报错;另外,其中的whereclaus...原创 2015-05-30 06:56:51 · 3858 阅读 · 0 评论 -
地理学与水利工程学科体系示意图
原创 2017-11-05 16:50:33 · 797 阅读 · 0 评论 -
Arcgis中添加XY坐标的数据
通过arcgis添加具有XY坐标的excel文件,要求具有经纬度坐标数据不再多说。(1)但有时候会报错,比如加载 .xlsx时说未加载项目吧啦吧啦,这时候需要转换呈.csv或者.xls格式。(2)要求加载的.csv文件的名字不要太长,不要"—”存在;另外,文件夹的路径不要太长。(3)字段的名称要为字母或者字母+数字,不可以为存数字。原创 2017-06-24 20:33:04 · 28248 阅读 · 1 评论 -
AHP方法计算权重流程示意图与代码实现
关于AHP的方法不再多说原创 2014-07-27 11:07:24 · 16510 阅读 · 0 评论 -
Matlab读取带有地理信息的图像
MATLAB 遥感影像的读写操作,带空间坐标信息!filepath=''E:\Sciencetific Work\北洛河_状头.tif';'[Data, R] = geotiffread(filepath); %读取带有空间信息的地理数据,返回数据,和地理栅格空间参考信息spatialref.GeoRasterReferencei转载 2016-10-18 17:21:43 · 20473 阅读 · 1 评论 -
雪水当量数据(EASE-Grid Data)处理方法
雪水当量数据(EASE-GridData)处理流程1、 数据描述该数据描述的是冰雪覆盖融化后相当于的水量值。是由欧空局(ESA)数据用户要素支持的全球雪盖项目支持,覆盖了北半球雪要素:雪范围(SE)和雪水当量(SWE),时尺度分别为17年和34年。其中SWE是被动微波和地面气象观测融合而成,时间覆盖为1979-2013日尺度。空间分辨率是25km,全球覆盖北纬35-85。原创 2016-09-16 20:41:10 · 11099 阅读 · 7 评论 -
python实现栅格计算器功能
对于Python不能直接调用arcgis中的栅格计算器工具。但是可以通过使用Raster函数实现。如下面代码实现了每年12个月做一个累加,然后输出一个年值的功能:for year in range(1960,1965): print year temp=0 #必须提前定义要不报错 for month in range(1,13): p原创 2015-05-21 22:20:39 · 6171 阅读 · 1 评论 -
python字符串函数
生成字符串变量str='python String function' 字符串长度获取:len(str)例:print '%s length=%d' % (str,len(str))连接字符串sStr1 = 'strcat'sStr2 = 'append'sStr1 += sStr2print sStr1复制字符串#strcpy(sStr1,sS转载 2015-01-21 17:55:15 · 864 阅读 · 1 评论 -
tiff与tfw
TIFF World File (TFW)格式说明摘要:TFW文件包含相关的TIFF文件的空间参数(spatial reference)数据。本文详细描述了TFW文件中定义的参数。 概述 两个文件构成TFW格式。 ASCII 头文件 图象文件 图象文件为TIFF格式(未压缩)。 详细转载 2014-11-19 20:51:04 · 5103 阅读 · 0 评论 -
遥感数据下载地址
1、全球陆表特征参量产品生成与应用研究-GLASSGeneration and Applications of Global LAnd Surface Satellite (GLASS) Products介绍:原创 2014-10-15 17:38:42 · 2615 阅读 · 1 评论 -
oracle 11g 用户密码过期问题
1、Oracle 11g 之前默认的用户时是没有密码过期的限制的,在Oracle 11g 中默认的profile启用了密码过期时间是180天。如下查看有效期:select * from dba_profiles where profile='DEFAULT' and resource_name='PASSWORD_LIFE_TIME';如果想设置密码不过期,可用管理员登陆,然后执行: AL转载 2014-01-16 09:51:00 · 775 阅读 · 0 评论 -
ArcGIS二次开发读写access,gdb数据库
privatevoid StartOper() { IWorkspaceFactoryWSF1 =null; IFeatureWorkspacepFeatureWorkspace =null; IWorkspaceFactoryWSF2 =null; IFeat原创 2014-01-17 10:32:04 · 2926 阅读 · 0 评论 -
ArcGIS AO二次开发环境搭建及入门
每次做arcgis的开发时都为开发环境的搭建花费好长时间,可以说是出现多种错误,现把环境搭建到初步读取要素的步骤总结如下,对于会出现的错误已经尽量标出。希望对初学者有帮助。 首先确保安装了arcgis desktop(版本越新,权限越高越好。如arcgis10,arcinfo权限);安装了ArcObject SDK for the Microsoft .NET Framework开原创 2014-01-17 14:13:47 · 3850 阅读 · 0 评论 -
ARCGIS常见拓扑错误的修改方法
一、先看看在ArcInfo里面的拓扑的实现: 其中在ArcGIS 8以后,Geodatabase成为一种全新的空间数据模型。它采用一种开放的结构将空间数据(包括:矢量、栅格、影像、三维地形等)及其相关的属性数据统一存放在工业标准的数据库管理系统DBMS中。空间要素类(Feature Class),如:河流、国界、宗地、电杆等,对应了DBMS中的表,而具体的一个要素(Feature)则是表转载 2014-01-10 15:16:18 · 18527 阅读 · 0 评论 -
WGS1984 -UTM投影问题
关于WGS1984投影坐标系的相关知识:打开WGS1984投影文件,仔细看看,我们可以发现里面中有三种不同的投影文件:如下:WGS1984 BLM Zone 14N(ftvs).prjWGS 1984 Complex UTM Zone 20N.prj (该处由20N——30N)WGS 1984 UTM Zone 9s.prj(该处由9s——60s)此处的S代表南半球,同样北半球有同转载 2014-03-21 09:05:31 · 23212 阅读 · 0 评论 -
2014年中科院地理所博士考试真题-地理信息系统
1、名词解释矢量数据模型;DEM;行程编码;TIN;地址匹配 2、 简单题地理空间及其存在形式;空间拓扑分析;最多路径分析;地表坡度计算方法;地理信息服务质量 3、论述题矢量与栅格数据的优缺点;互联网地图服务发展现状与瓶颈和突破点分析;遥感、GIS融合发展趋势分析原创 2014-07-30 12:13:44 · 1577 阅读 · 0 评论 -
根据DEM生成连续曲面栅格(提取DEM的经纬度)
1、原始的栅格高程图2、通过栅格转换点图层,提取栅格的中心点的高程值3、为点图层添加经纬度X,Y坐标属性(同理获取Y的属性坐标值。此外也可以通过datamanagement tools->featureàadd XYcoordinates方法添加经纬度坐标信息)4、点图层x,y,z值分别转为栅格图层 5、借助栅格计算器计算关于经纬度和高程的函原创 2014-05-20 15:36:28 · 9376 阅读 · 8 评论 -
arcgis的python入门
入门首先安装了arcgis的原创 2014-10-10 20:59:04 · 2261 阅读 · 0 评论 -
ArcGIS+Python读取flt文件并进行区域统计
1、http://sourceforge.net/projects/dbfpy/files/dbfpy/ 由于统计后生成的为dbf文件,孤下载dbfpy-2.3.0.win32.exe,专门读取写dbf文件的包原创 2014-10-21 11:18:43 · 3610 阅读 · 0 评论 -
python获取栅格点和面值
1、 获取指定位置的点值:OutputFile = open(statisticResultTXT, 'w') cellvalue=arcpy.GetCellValue_management(inputfilepath+filenname+".tif",staionXY, "1")OutputFile.write(stationID+"_"+filenname+""+str(cel原创 2015-04-08 08:55:14 · 4640 阅读 · 0 评论 -
arcgis脚本Python批量裁切与统计实现
首先通过批量裁切为arcgis的extractbymask工具,代码如下: 其次,利用arcgis的Python脚本的GetRasterProperties_management方法,代码如下:原创 2015-01-17 11:33:43 · 2655 阅读 · 0 评论 -
arcgis图例横放效果实现
默认出图时图例是竖直的,有时候是想达到图例横放的效果。如下图所示:方法步骤描述如下:1、 右击图层属性,设置符号化symbology为拉伸(stretch)样式。为了显示多个过渡色彩,然后点击标签(labeling),产生5个色彩即可。2、 插入图例,可以看到默认是竖向的。通过打开属性,点击样式style,然后选择第一个即为横向显示。3、原创 2015-01-03 23:41:23 · 43041 阅读 · 13 评论 -
ArcGIS操作之创建要素
1、属性添加自动提醒:设置的目的就是让创建要素时,提醒输入对应的属性信息,以防遗忘输入属性值。位置:编辑工具条—选项,然后编辑选项下的属性选项卡,可对单个图层或多个图层进行设置。2、创建标签create label:设置标签的目的是为了快速定位到想要的显示区域,方便操作。位置:标签标题—创建,即为当前视图创建了一个标签,创建成功后可在“标签“下看到。管理标签的选项可以重命名、删除、上下移动已有原创 2015-01-04 21:16:40 · 2930 阅读 · 0 评论 -
arcgis缓冲区获取的外边界方法
1、通过找到所有节点的链表,删除内环的所有节点坐标,仅保留外部节点即可2、通过多边形转为多义线,然后select外面的多义线,然后把该polyline转为polygon即可。原创 2015-02-07 00:25:28 · 6033 阅读 · 0 评论 -
学习arcgis的脚本Python入门技巧
学习arcgis的脚本Python的快捷方式是构建自己的modelbuilder,然后导出为Python脚本,通过查看其脚本就知道如何调用arcgis的工具箱了。然而arcgis10.0后的栅格计算器不再使用9.3的用法,如下在10.0中会报错。# Process: Raster Calculatorarcpy.gp.RasterCalculator_sa("\"%yubidu%\"* 5原创 2015-01-17 16:03:13 · 4604 阅读 · 0 评论 -
栅格数据挑选子区域
对于一个大的栅格数据,有时候想提取出其子区域,例如四川省栅格和成都市的栅格,如果想获取除去成都市的四川栅格图层。那么,首先想到的是获取到除去成都的四川范围的矢量图层,然后使用extract by mask方法裁切出来。但是如果只有两个大小栅格图层,想获取mask图层估计需要一些额外操作(如栅格转矢量,然后再合并等才能得到两个边界,还要做clip才能得到mask图层,此处不再赘述)。原创 2014-12-14 21:04:49 · 1533 阅读 · 0 评论