gis | 易学教程

Download and read shapefile function in R

阅读更多关于 Download and read shapefile function in R

问题 I would like to expand on this function. As of now, the function downloads and unzips the shape file from the web. I would like to implement 'rgdal' to read the file into R. library(rgdal) dlshape=function(location) { temp=tempfile() download.file(location, temp) unzip(temp) } I found the following code on SO, but I was unsuccessful in adapting it. It appears that the function still looks at the first file unzipped rather than grep for a file ending with the .shp extension. read.csv.zip <-

ArcGIS实验教程——实验二十一：DEM分析

阅读更多关于 ArcGIS实验教程——实验二十一：DEM分析

一、实验描述表面分析主要通过生成新数据集，如等值线、坡度、坡向、山体阴影等派生数据，获取更多的反应原始数据集中所暗含的空间特征、空间格局等信息。二、实验内容 1、地形因子计算 2、填挖方计算 3、通视分析 4、剖面线计算 5、山体阴影 6、曲率计算 7、等值线三、实验目的 1、掌握表面分析的方法四、实验数据 lanzhouDem 五、实验过程加载数据： dem 1 、地形因子计算（1）坡度计算坡度指过地面人任意一点的切平面与水平地面之间的夹角。 Slope （2）坡向计算坡向指地面上一点的切平面的法线矢量在水平面的投影与该点的正北方向上的夹角。 Aspect 2 、填挖方计算填挖方是一个通过添加或移除表面物质来修改地表高程的过程。原始高程设计高程填挖方 3 、通视分析 4 、剖面线计算 5 、山体阴影山体阴影通过考虑照明源的角度和阴影，根据表面栅格创建晕渲地貌。 6 、曲率计算地面曲率是对地形表面上一点扭曲变化程度的定量化度量因子。曲率计算的输出结果为每个像元的表面曲率，该值通过将像元与八个相邻像元拟合而得。曲率是表面的二阶导数，即坡度的坡度。曲率 7 、等值线等值线就是将表面上相邻的等值点（如高程、温度、降水、大气压力等）连接起来的线。等高线本实验到此结束！专注GIS教育，探索GIS奥秘，分享GIS价值！来源： CSDN 作者：

Loading JSON into a GeoDataFrame

阅读更多关于 Loading JSON into a GeoDataFrame

问题 I'm having difficulty loading the following JSON containing GIS data (https://data.cityofnewyork.us/resource/5rqd-h5ci.json) into a GeoDataFrame. The following code fails when I try to set the geometry. import requests import geopandas as gpd data = requests.get("https://data.cityofnewyork.us/resource/5rqd-h5ci.json") gdf = gpd.GeoDataFrame(data.json()) gdf = gdf.set_geometry('the_geom') gdf.head() 回答1: Setting the geometry fails because the geopandas.GeoDataFrame constructor doesn't appear

How do I plot a “differences” map using geopandas and mapclassify?

阅读更多关于 How do I plot a “differences” map using geopandas and mapclassify?

问题 I have a GeoDataFrame with a "differences" column that stores the delta between two numbers. Sometimes these numbers are positive, negative or zero if there is not difference. I need to make a choropleth map that has some divergent colormap with 0 (or some middle buffer) as the midpoint and non-symetrical colorbar, for example [ -0.02, -0.01, 0., 0.01, 0.02, 0.03 ]. I've experimented with class MidPointNormalize(mp.colors.Normalize): def __init__(self, vmin=None, vmax=None, midpoint=None,

免费GIS平台

阅读更多关于免费GIS平台

http://bgming.cnblogs.com/archive/2005/10/18/257253.aspx 本站主要介绍免费WebGIS平台，它们具有共同特点：免费使用，开放源代码（具体使用规则请参看各个平台的 Lisense文件），希望借此带您进入广阔的开源GIS世界。采用Java开发的几个主要的开放源代码的WEBGIS平台列表（由于时间有限，只选取部分平台）平台当前版本生成日期项目地址授权协议支持格式和接口 ESRI Shape Mapinfo MID/MIF Oracle MySQL PostGIS 其他格式*2 OGC*8 WMS*9 Server OGC WFS*10 Server OGC WCS*11 Server MS .Net Deegree 2.0 2005-2-3 http://deegree.sourceforge.net LGPL √ √ √ √ √ GML2.1.1 √*4 √*5 √*6 Openmap 4.6.2 2005-2-15 http://openmap.bbn.com OpenMap License √ √*1 √*3 √ VPF RPF DTED √ √*7 GeoServer 1.3 2005-6-9 http://geoserver.sourceforge.net GPL 2.0 √ √ √ √ GML2.1.1

修改ArcGIS Server Account / 站点管理员的账号、密码

阅读更多关于修改ArcGIS Server Account / 站点管理员的账号、密码

安装ArcGIS Server（10.1之后版本）时会先后创建两个账户： ArcGIS Server Account 和站点管理账户 ArcGIS Server Account：是操作系统账户，Server访问机器上的数据需要通过该账户，因此发布服务时经常需要保证对涉及服务的路径赋予ArcGIS Server Account 读写权限（读写权限的赋予：对应存放数据的文件夹上右键→属性→ 安全赋予arcgis server account（如果安装时没有改名字，默认名字为arcgis）读写权限，如果列表中没有该账户，点编辑→ 添加→ 输入对象名称→ 检查名称）站点管理账户：浏览器端登录Server Manager时的账号密码 ArcGIS Server Account 修改密码先将之前的账户在本地用户和组中删除掉，建议配置时新用户名与原来的相同。第一步：在操作系统:计算机管理→ 本地用户和组→ 用户对ArcGIS Server Account 右键设置新密码第二步：开始→ 所有应用→ ArcGIS→ Configure ArcGIS Server Account→ 设置新密码后面基本一直下一步就行~ 不用填写路径什么的~ 如果没有在计算机管理里先设置新密码，这里操作会显示invalid password..... 注：①因为这里修改密码时无需填入旧密码

ArcGIS 10.2安装之 Server的安装

阅读更多关于 ArcGIS 10.2安装之 Server的安装

第二部分 ArcGIS Server的安装、部署注意：需关闭电脑防火墙，否则arcgis manager可能登陆不上去！ 1 、ArcGIS Server的安装（1）双击打开安装程序，选择ArcGIS for Server，安装接受许可，进入下一步；选择安装目录，进入下一步。给arcgismanager配置密码，此处把账户名修改为了Administrator，密码自行设定，进入下一步，如果下次配置，还需要该文件，则导出配置文件，最好别导出，选择，请勿导出配置文件，进入下一步，点击安装，直到安装完成；安装完毕后，进入授权阶段，如下选择安装目录下的KeyGen.exe文件，打开，选择server标签，Feature处选择arcgisserver，Version处选择102（如果不在列表中，手动输入），表示版本为10.2，点击All生成license文件，将生成的内容，复制，新建一个.txt文本，将其粘贴进去，另存为.ecp文件，在上一步浏览按钮打开后，查找到该ecp文件，即可实现注册。（注：该处为ESRI厂家提供测试使用，如果正式使用，需由厂家提供相应的注册文件）。选择第三项，我已收到来自ESRI的授权文件，并且准备好完成授权过程。点击浏览按钮，选择生成的.ecp文件，点击进入下一步；点击下一步，完成安装；（4）配置ArcGIS Server

GIS网络资源

阅读更多关于 GIS网络资源

目前共发布95个数据集的信息，涉及的总数据量超过25GB (不包括遥感影像数据集)。大部分数据为带密级的数据，不能通过网络直接获取。其中共有49个数据集的数据和24个样例数据可以直接免费下载，数据量达210MB。对于其中带密级的数据，分别给出了数据获取途径、程序以及需要遵照的法规。地理信息系统论坛 http://www.gisforum.net/ 中国GIS咨询网 http://www.gissky.com/index.asp 资源大联盟 http://zydlm.wxhc.com.cn/main.asp FreeGIS组织 http://www.freegis.org 该网站提供FreeGIS平台索引，国际上主要的FreeGIS平台都可以在这里找到 MapTools组织 http://www.maptools.org 该网站提供开放源代码的GIS资源 OGC组织 http://www.opengeospatial.org 国际开放GIS标准的制订者，ESRI,Intergraph,Mapinfo,Oracle,Sun,BEA等都是OGC成员 MySQL数据库 http://www.mysql.com MySQL从4.1开始支持空间数据库，遵循OGC规范，目前正式版本是5.0 PostgreSQL http://www.postgresql.org 中文 http://www

From Voronoi tessellation to Shapely polygons

阅读更多关于 From Voronoi tessellation to Shapely polygons

问题 from a set of points I built the Voronoi tessellation using scipy: from scipy.spatial import Voronoi vor = Voronoi(points) Now I would like to build a Polygon in Shapely from the regions the Voronoi algorithm created. The problem is that the Polygon class requires a list of counter-clockwise vertices. Although I know how to order these vertices, I can't solve the problem because often this is my result: (overlapping polygon). This is the code (ONE RANDOM EXAMPLE): def order_vertices(l): mlat

GIS网址大全

阅读更多关于 GIS网址大全

最近新开通了网站GIS网址大全，收集所有有关GIS资源的权威网站。供大家快速查找对应资源。 IS网址大全,GIS网站导航,GIS学习,GIS教程导航,GIS平台,GIS企业,WebGIS网址同时也需要大家的共同维护，可以在页面上添加网站链接，审核通过后会展示在网站首页。网站首页： GIS网址大全 http://www.gisurl.cn 来源： CSDN 作者： gisdoer 链接： https://blog.csdn.net/gisdoer/article/details/104253242

订阅 gis