Questions tagged «shapefile»

矢量数据存储格式，用于存储地理特征的位置，形状和属性。当ESRI开发它时，通常称为“ ESRI shapefile”。

1
如何重新投影shapefile？
我想重新投影一个QGIS项目，但是没有这样的选择。当寻找答案时ogr2ogr，通过重新投影每个矢量层（在此示例中为shapefile），我看到了该任务的建议。但是我没有得到预期的结果： C:\temp>ogrinfo -al -so misc.shp INFO: Open of `misc.shp' using driver `ESRI Shapefile' successful. Layer name: misc Geometry: Line String Feature Count: 10 Extent: (21.267388, 42.015857) - (21.270225, 42.017470) Layer SRS WKT: PROJCS["MGI_Balkans_zone_7_deprecated", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",21], PARAMETER["scale_factor",0.9999], PARAMETER["false_easting",7500000], PARAMETER["false_northing",0], UNIT["Meter",1]] id: Integer (10.0) C:\temp>ogr2ogr -t_srs …
5
如何从非常详细的陆地区域shapefile获取海洋shapefile？
我正在寻找“反转” shapefile的可能性。我想要做的是将具有大陆和岛屿的shapefile转换为具有水多边形的shapefile。来自自然地球的现有海洋shapefile太不准确。 作为OSM的积极贡献者，我认为使用OSM的现有shapefile对我来说很好。从GIS的角度来看，这些方法更为准确，但存在很多错误。但是，mapserver可以处理它，所以我认为这些应该没问题。这些shapefile，其中一些您可以在http://metro.teczno.com/#coastline 上找到，它们 非常大。 我创建了一个shapefile，其中包含一个带世界范围（-180到180，-90到90）的多边形。然后，我尝试从我的世界范围shapefile中裁剪出陆地区域在另一个shapefile中的区域。我尝试使用QGIS，矢量工具->差函数。在100％CPU使用率的12小时后，我停止了它。它始终显示进度为100％。RAM不是问题。 我正在寻找另一种可能性，可以从“世界多边形”（最好是Linux命令行）中裁剪出陆地区域。
3
使用QGIS在给定多边形上的Voronoi图
我正在尝试在特定的多边形上创建Voronoi图，但是在处理边缘时会出现一些问题。这是我的方法： 为Voronoi和图的多边形准备点 创建Vornoi 使用多边形切割边缘 使用Union组合切割的Voronoi和多边形（不确定这是否是最好的方法），这是我的边缘问题。我想：1）从原始Voronoi地图中删除框架（标记为蓝色）2）将边缘扩展到区域多边形的边缘（标记为红色），以使边缘侧的服务区域更加逼真。 我怎样才能做到这一点？
1
在QGIS中对shapefile运行SQL查询吗？
有什么方法可以将shapefile导入QGIS，然后对它运行空间SQL查询？ 我显然ogr2ogr -sql可以这样做，但是每次都要将输出导入QGIS来查看它，这很痛苦。 如果我可以将shapefile导入QGIS，然后对它运行SQL并立即查看结果，那将是很好的。
9 qgis  shapefile 
2
ggmap：从shapefile绘制多边形
使用ggmap，我想在包含一些位置点的地图上包括来自shapefile的市政边界（多边形）。除了绘制多边形外，此脚本会执行所有操作： library(rgdal) library(ggmap) # Get shapefile with Drammen municipality borders tmpzip<-tempfile() tmpdir<-tempfile() dir.create(tmpdir) download.file("http://www.kartverket.no/Documents/Kart/N50-N5000%20Kartdata/33_N5000_shape.zip",tmpzip) unzip(tmpzip, exdir=tmpdir) kommune <- readOGR(dsn=tmpdir, layer="NO_AdminOmrader_pol") kommune<-kommune[kommune$NAVN=="Drammen",] kommune<-spTransform(kommune, CRS("+init=epsg:4326")) # Get location point data subscr<-data.frame(lon=c(10.1237,10.2161,10.2993),lat=c(59.7567,59.7527,59.6863), pop=c(58,12,150)) coordinates(subscr)<-~lon+lat proj4string(subscr)<-CRS("+init=epsg:4326") lon <- c(10.0937,10.3293) lat <- c(59.7916,59.6563) map <- get_map(location = c(lon[1], lat[2], lon[2], lat[1]), maptype = "roadmap", source …
9 shapefile  r  ggmap 
1
将R中的readShapePoly坐标转换为经纬度坐标
我正在尝试为丹麦市政当局查找shapefile的质心，然后找到它们之间的行驶时间。我使用R的readShapePoly从功能maptools与组合gCentroid，从功能rgeos，和一切正常。但是，我得到诸如 SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA 我显然不能在Google中使用这来获取旅行时间。我正在寻找一种将这些数字转换为经度-纬度的方法，但不知道如何。 当我readOGR从rgdal库中读取数据时，我得到了相同的坐标，但是它告诉我以下有关我假设是投影的信息（但是坐标是相同的） Slot "proj4string": CRS arguments: +proj=utm +zone=32 +ellps=intl +units=m +no_defs 可重现的示例：我将示例的数据放在这里：https : //github.com/sebastianbarfort/shapefiles 这应该重现该问题： library(maptools) library(rgdal) library(rgeos) map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp") centroid = gCentroid(map) centroid
2
使用OGR / GDAL确定shapefile和栅格在Python中是否重叠？[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我正在使用OGR / GDAL在python中构建脚本。 我有一组shapefile和一组GeoTiff栅格文件。 如果脚本不与栅格区域相交，我希望脚本忽略它们。 shapefile不是矩形，因此我不能简单地比较layer.GetExtent（）返回的xmin / xmax，ymin / ymax值。我需要代表其整体形状的实际多边形，然后通过某种方式确定该多边形是否与栅格正方形相交。 我当时想我可以以某种方式将shapefile中的所有多边形合并为一个要素，然后读取该要素上的几何图形，然后将该信息与栅格范围进行比较。但是，我不确定具体如何执行此操作。 如何从shapefile中提取边界多边形信息？ 如何确定该多边形是否与给定的正方形相交？
5
GIS分析以查找重复的几何
我有一个很大的shapefile，其中包含我工作所在的城镇的所有建筑物和房屋（约90,000个要素）。建筑物/房屋的数据由镇的测量工程师保存，并且由于不当做法以及不同测量师对数据的访问，许多建筑物/房屋已被保存两次并在地图上显示为重复项。 其中一些是完全重复的（它们一个出现在另一个之上），而另一些则是在两个对象之间留有间隔的（例如，一个对象在另一个对象内部，请参见所附的屏幕截图）。 我想清理这些数据，以便在城镇中只有正确的建筑物/房屋，所以我的问题是： 我是否可以运行任何GIS分析或SQL表达式来查找所有重复的功能（确切的功能和位于其他功能内的功能）？我同时拥有ArcGIS和QGIS，因此欢迎您提出所有建议。
3
如何从命令行查找shapefile的属性？
我正在关注Mike Bostock的教程“让我们绘制地图”。他介绍了如何ogr2ogr根据特征的属性来过滤shapefile 。我可以在QGIS中打开shapefile，然后查看属性表以找出那些功能。有没有办法从命令行中找到该信息？
3
将shapefile转换为SVG并保留内部数据？
我对GIS非常陌生。我找到了有关如何使用Python工具制作合奏的教程。我想遵循它，但是使用我自己的数据和地图。但是，本教程使用了美国地图SVG文件，我需要使用CA邮政编码地图。我在美国人口普查的shapefile中找到的唯一地图。我尝试使用kartograph.py将其转换为SVG，但是这样做时，它会丢失文件中的邮政编码。如何将文件转换为SVG，并保持每个路径的ID为适当的邮政编码。 我正在尝试遵循“ 如何使用免费工具制作美国县专题地图”教程。只有我没有SVG格式的地图，这应该可以使您了解我要做什么。 我看过对javascript实用程序的引用。我从未使用过Javascript，但想想我可以使用简单的东西。
6
如何以编程方式检查形状的数量=表记录的数量？
我有大约1000个shapefile已损坏（请参阅附件的错误消息）。shapefile是从eCognition Developer 8生成的。有一个脚本工具似乎可以在将shapefile识别为损坏后对其进行修复。 编辑： 我想创建一个快速脚本来遍历我所有的shapefile，并检查形状的数量是否与表记录匹配。我可以使用以下方法对表记录进行计数： # Name: fcCount.py # Purpose: calculate the number of features in a featureclass # Import system modules import arcpy from arcpy import env env.workspace = "C:/data" Sample = "MyShp.shp" result_dbf = int(arcpy.GetCount_management(Sample).getOutput(0)) print result_dbf 我最终想创建某种逻辑检查，例如： if result_dbf = result_shp: pass else: print "There is a …
2
如何在没有ArcGIS的情况下将GDB数据转换为shapefile？
我们从GDB文件夹中的第3方获得了一些数据（应该是点的简单集合，或者至少是我们需要提取的全部数据）（我不知道如何知道哪个版本），但是我们没有没有ESRI Arc *产品。我们主要不是GIS公司，我们的工具基本上仅支持shapefile。我们如何将GDB转换为shapefile？ 注意：我尝试按照说明在QGIS中获得GDB支持（在卸载QGIS并通过osgeo4w-setup重新安装后），但是它不起作用。QGIS仍说“ Xyz.GDB不是有效或公认的数据源”。
4
USGS数据有哪些替代来源？
我一直在使用USGS国家地图查看器（http://viewer.nationalmap.gov/viewer/）来获取QGIS的数据，但是下载系统使我发疯。我输入我的电子邮件地址，并[有时]通知我，我的数据将在最多2天内可用（为什么这么长？我很不幸。 还有其他地方可以获取数据吗？这个“点击并希望您的数据能在2天之内到达”的系统使我麻烦不断，工作陷入困境。我不明白为什么我可以按需从USGS下载多个GB的海拔数据，但是水力数据的大小却需要等待一段时间。 编辑：我专门在寻找华盛顿州基蒂塔斯山谷的水利，运输和土地覆盖。感谢您的出色回答。
9 shapefile  data 
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.