32得票3回答
错误:isTRUE(gpclibPermitStatus())不为TRUE。

这个问题可能是之前未回答的一个问题的重复。我仍然有这个问题。 我正在尝试使用一个邮政编码shapefile,但出现以下错误:tract <- readOGR(dsn = ".", layer = "cb_2013_us_zcta510_500k") tract<-fortify(...

31得票1回答
将GDAL图层转换成栅格图像

编辑 这里是正确的方法,以及文档:import random from osgeo import gdal, ogr RASTERIZE_COLOR_FIELD = "__color__" def rasterize(pixel_size=25): # ...

31得票13回答
如何在Python或JavaScript中将UTM坐标转换为经纬度

我有一些以UTM格式表示坐标的文件。对于每个坐标,我都有东经、北纬和区域。我需要将其转换为LatLng格式,以便使用Google地图API在地图上显示信息。我找到了一些在线计算器可以完成这项工作,但没有实际的代码或库。http://trac.osgeo.org/proj4js/是Javascr...

31得票7回答
如何在球面上计算点到线段的距离?

我在地球上有一条线段(大圆弧部分),该线段由其端点的坐标定义。显然,两个点定义了两条线段,因此假设我对较短的那条线段感兴趣。 我得到了第三个点,并正在寻找线段和该点之间的(最短)距离。 所有坐标都以经度\纬度(WGS 84)给出。 如何计算距离? 任何合理的编程语言解决方案都可以。

30得票3回答
将geopandas数据框中的Shapely点几何对象转换为纬度/经度列。

我有一个包含shapely点对象列的geopandas数据框。我想从shapely点对象中提取坐标(纬度/经度)以生成纬度和经度列。应该有一种简单的方法可以做到这一点,但我无法想出来。 我知道可以像这样提取单个坐标: lon = df.point_object[0].x lat = df....

30得票3回答
R能用于GIS吗?

我想创建一些GIS图表,想知道R语言能否用于这个任务。以下是一些类似于我想制作的图表示例: 一个美国的温度图(或等高线图),根据各州GDP的大小决定颜色(或高度)。因此,州界将在生成的图中产生不连续性。 一个美国的温度图,其中使用海拔高度进行数据分析。在这种情况下,结果图应该在各州边界之间...

29得票1回答
如何在GIS中选择MongoDB和PostgreSQL之间的区别?

我正在尝试找一些关于这个主题的信息,但大多数文章都是三年以上的。我也还比较新手,不知道该问谁。 就我的使用情况而言,我并不是数据库专家,但我以前用过几次MongoDB。我只在一个小项目中使用过PostgreSQL(没有GIS相关内容)。目前我只有一小部分数据存储在dynamoDB中。虽然它有一...

29得票4回答
使用MongoDB或Cassandra处理空间数据

我正在考虑一个处理大量数据的概念验证,例如大于10G的空间相关数据,需要每秒至少200多次写入和50多次读取。这也是一个不断增长的系统。目前,出于性能原因,我正在考虑将这个大容量数据移动到一个NoSql大表类型的数据库中。 我考虑并仔细研究了MongoDB和Cassandra。就我的阅读而言...

28得票3回答
在Postgis中,使用SRID 4326计算两个点之间的距离(单位为米)

这可能是一个简单的问题,但我并不擅长PostGIS,并且不完全理解所有内容。 基本上,我有一个表(nodes),其中包含一个点列(point)。我已经在这一列上创建了一个索引。create index nodes__points on nodes using gist (point) 该列是...

28得票6回答
如何在“popupopen”事件期间识别Leaflet的标记?

当点击标记时,我需要执行一些代码,找到对应于被点击的标记的id,从后端API检索数据,然后将新检索到的数据添加到将要打开的弹出窗口的content中。 唯一能够监听标记点击事件的方式是:map.on('popupopen', function(e){ // How to retriev...