地理信息系统

我想将不同的CRS（主要是WGS84纬度/经度）中的许多shapefile转换为通用投影（可能是Albers Equal Area Conic），但是一旦我的问题变好，我可能会寻求其他问题的帮助-定义）。 我花了几个月的时间在R中进行空间统计工作，但那是5年前了。对于我的一生，我不记得如何将一个sp对象（例如SpatialPolygonsDataFrame）从一种投影转换为另一种投影。 示例代码： P4S.latlon <- CRS("+proj=longlat +datum=WGS84") hrr.shp <- readShapePoly("HRR_Bdry"), verbose=TRUE, proj4string=P4S.latlon) # Shapefile available at # http://www.dartmouthatlas.org/downloads/geography/hrr_bdry.zip # but you must rename all the filenames to have the same # capitalization for it to work in R 现在，我有了一个SpatialPolygonsDataFrame具有适当的投影信息，但是我想将其转换为所需的投影。我记得有一个不太直观的函数，但是我不记得它是什么。 请注意，我不仅要更改CRS，还希望更改坐标以匹配（“重新投影”，“变换”等）。 编辑 不包括在此shapefile中令人讨厌地放置在墨西哥的AK / HI： library(taRifx.geo) hrr.shp <- …

35 r  coordinate-system 

随着arcpy中数据访问模块的引入（搜索光标速度提高了30倍），我想知道与sql条件匹配的功能计数是否比传统的MakeTableView + GetCount方法更快？

35 arcgis-desktop  arcpy  arcgis-10.1  performance 

我有一个包含2亿个观测值的CSV，格式如下： id,x1,y1,x2,y2,day,color 1,"-105.4652334","39.2586939","-105.4321296","39.2236632","Monday","Black" 2,"-105.3224523","39.1323299","-105.4439944","39.3352235","Tuesday","Green" 3,"-104.4233452","39.0234355","-105.4643990","39.1223435","Wednesday","Blue" 对于每组坐标（x1 / y1和x2 / y2），我想分配它所属的美国人口普查区或人口普查区（我在此处下载了人口普查区TIGER形状文件：ftp : //ftp2.census.gov/ geo / tiger / TIGER2011 / TRACT / tl_2011_08_tract.zip）。因此，我需要为每个观察做两次多边形点操作。比赛必须非常准确。 最快的方法是什么，包括花时间学习软件？我可以访问具有48GB内存的计算机，以防万一这可能是一个相关的限制。 有多个线程建议使用PostGIS或Spatialite（Spatialite看起来更易于使用-但它与PostGIS一样有效吗？）。如果这些是最佳选择，是否必须填充空间索引（RTree）？如果是这样，怎么做（例如使用人口普查道Shapefile）？对于包含示例代码（或示例代码的指针）的任何建议，我将不胜感激。 我的第一次尝试（在找到此站点之前）包括使用ArcGIS对美国人口普查数据块的数据（100,000点）的子样本进行空间连接（仅x1 / y1）。在我终止该过程之前，花了5个多小时。我希望可以在不到40小时的计算时间内对整个数据集实施的解决方案。 很抱歉提出之前提出的问题-我已经阅读了答案，但不知道如何实施建议。我从未使用过SQL，Python，C，并且以前只使用过一次ArcGIS-我是一个完整的初学者。

35 postgis  spatialite  point-in-polygon  large-datasets  optimization 

全球定位系统（GPS）是一项免费服务，由美国政府拥有并运营，并且始终可用。 购买GPS设备时，我们无需支付月费或为GPS支持支付税款。我们只支付设备的价格。 我不知道怎么可能。美国如何赚钱？

35 satellite  gps 

我有一个多边形Shapefile，需要对其进行概括。如何将poligon泛化，并确保没有诸如“条子”之类的拓扑错误。 我已经尝试过： 矢量->几何工具->多边形到线 矢量->几何工具->简单的几何...在这里我已经可以在屏幕上看到，但是我正在产生重叠和间隙，现在我尝试 矢量->几何工具->到线条的线 此处程序显示错误并崩溃

35 qgis  topology  generalization 

我一直在编写自定义代码来处理GPS轨道：手机中的GPX文件和NMEA日志，专用GPS记录器等。我感觉自己在重新发明轮子，并且想知道是否有一些标准的库/工具可以用来工作GPS轨道。我喜欢做的一些事情： 计算在轨道上行驶的距离 计算沿轨道的平均速度和瞬时速度 计算运动时间与静止时间 简化轨道，消除固定路段 将轨道分段为单独的轨道，每趟一次 消除GPS采样误差 通过在DEM中查找点来创建高程剖面 将曲目转换为GPX，KML等以进行可视化 有许多分析GPS轨迹的应用程序。我正在寻找软件库或可编写脚本的工具集，以便可以编写自己的应用程序。我不太在乎输入格式；任何带有时间戳的经纬度点都可以。我主要在Unix上使用Python工作，但几乎所有语言都可以。 我使用的现有工具包括GPSBabel和GDAL / OGR用于转换，以及PostGIS用于将轨迹几何图形作为线串使用。但是，大多数情况下，我最终还是编写了自己的Python代码，该代码可用于经纬度对数组，而且感觉太低了。同样，上述某些任务（例如消除GPS错误）非常微妙。我应该使用一些现有的库吗？

35 gps  gpx  tracklog  movement-data 

为什么将数据存储在异常的投影和基准中很有用？ 由于失真等原因，我理解投影作为输出构造的价值：这是合理的。 但是，我不明白为什么例如状态在数据中使用状态平面投影。我们有十进制精度：为什么不只将精确值存储在EPSG：4326中？这是否完全是自动重投影前的剩余时间，还是我缺少一个价值主张？ 我将这个问题的范围限制为矢量数据，只是为了使其更加具体。

35 coordinate-system 

我目前知道以下开源GIS服务器： Geoserver：基于Java，也支持矢量和位图。现在与mapserver相当（性能方面）。 mapserver：c / c ++，曾经是最快的服务器。比矢量更适合栅格吗？ mapguide：不知道这个 mapnik：街上的新孩子？对其了解不多，但看起来很吸引人 mapfish：据我所知，仅检索矢量数据以在openlayers中显示。确实存在Rails实施。 如果可能的话，我想进行某种比较，您选择还是喜欢哪个，为什么？ 我正在寻找一个Rails网站，并且需要某种GIS服务器。我将需要栅格和矢量数据（可单击）。一台服务器可以做到吗？后端数据库将是Postgis。

35 open-source-gis  software-recommendations  osgeo 

我想知道pgRouting工具是否可以替代。 仔细考虑一下，我认为唯一的选择是使用R或RPY进行分析。 基本上，我正在尝试解决基于非PostGIS数据库的路由问题。由于我的大多数数据实际上都以不同的格式存储；但是，转换并不难。我发现很难实现对愚蠢的大脑的pgRouting。因此，也许我可以轻松调整的非数据库解决方案对我来说是理想的选择

35 software-recommendations  open-source-gis  routing  network 

我希望能够将特定区域的openstreetmap数据用作源，并从基础结构改进的地图（主要是道路）中“借鉴”。 我希望它是自托管的（不破坏真实OSM服务器上的任何数据），希望在编辑过程中具有很大的灵活性和细节级别，并且以一种漂亮的方式呈现它会很好。 我对许多GIS软件不熟悉。实际上，以OSM数据作为源并不重要，它也可以是卫星图像。重要的事情-灵活的绘图必须绑定到地图的真实坐标和可视化。 我从哪里开始？ 希望这个问题在这里有意义:)

35 openstreetmap  visualisation  rendering  map-drawing 

我目前正在研究其中一个项目中CartoCSS和TileMill的用法。我已经在http://www.mapbox.com/tilemill/gallery/上看到了令人印象深刻的画廊 画廊看起来令人印象深刻，但是我很难想象它们如何达到某些效果和风格。指南中提供的示例（例如此处）似乎太基础了。 CartoCSS的任何详细样本（而不只是输出）都可以在某处找到吗？

35 tilemill  cartocss 

当使用Leaflet可视化大型数据集（具有10,000个点特征的GeoJSON）时，浏览器崩溃或挂起就不足为奇了。来自同一数据集的1000个要素的子样本可以完美地工作。不幸的是，我无法共享数据集供其他人试用。 有没有人有更好的解决方案来可视化如此大的数据集？（最终目标是将其扩展到200万个功能）我什至愿意考虑离线可视化框架，以防基于浏览器的替代品（例如Polymaps或d3.js等）无法使用。 编辑：忘记提及了，用户需要能够按属性过滤数据集。因此，在N个要素中，可能仅需要动态渲染匹配的n <= N个要素。

35 web-mapping  javascript  visualisation  leaflet  polymaps 

我想从PostGIS获得带有属性的geojson功能。我找到了一个具有功能部件集合的示例，但是我不能使其仅适用于功能部件。 SELECT row_to_json(fc) FROM ( SELECT 'FeatureCollection' As type, array_to_json(array_agg(f)) As features FROM (SELECT 'Feature' As type , ST_AsGeoJSON(lg.geog)::json As geometry , row_to_json(lp) As properties FROM locations As lg INNER JOIN (SELECT loc_id, loc_name FROM locations) As lp ON lg.loc_id = lp.loc_id ) As f ) As fc; 到目前为止，我尝试修改示例的要素收集查询。但输出无效。

35 postgis  postgresql  geojson  sql 

如何NULL在shapefile的属性表中识别哪些特征具有属性？ 当我尝试[field]= 'NULL'它选择与特点value = 0（见图片）

35 qgis  query  attribute-table 

我有一个栅格时间序列，其中每个栅格的像元在特定时间表示某个值。 我想生成一个图表，突出显示每个单元随时间变化的正值或负值趋势的大小。 我有些天真的方法是将简单的线性回归（X =时间和Y =值）拟合到每个像元，并将斜率数组输出到栅格（如下图所示）。只能通过导出重要结果来进行过滤。 在栅格时间序列中，我还能如何表示一段时间内的趋势？ 我对通用技术感兴趣，而不对软件专用说明感兴趣。

34 raster  time  gis-principle