标签: shapefile

坐标时区

由于标题推断我需要根据一对坐标找到时区(或者可能只是UTC偏移).我一直在寻找不同的解决方案,但有一些Web服务,但我需要能够脱机访问应用程序.由于时区并非完全基于经度,因此看起来并不容易......

我虽然要查询一个ESRI shapefile,但我已经包含了世界上所有国家和它们的时区,但它看起来有点复杂.如果这应该是解决方案,您知道任何提供此功能的.NET库吗?

.net c# timezone esri shapefile

8
推荐指数
2
解决办法
3922
查看次数

将ESRI Shape SHP文件的子集拆分或保存到新文件中?

我正在使用GeoDjango中的形状文件.现在我正在尝试为加载在形状文件中的代码编写测试并将其保存到数据库中.形状文件目前的功能数量为64,118.我想将其减少到极少数,以便测试可以快速加载并确认一切正确.

由于形状文件不是文本格式,是否有一个免费的应用程序或库可以用来提取一些功能并将它们保存到新文件中?

我应该提到我没有许可证也无法访问任何ESRI产品系列.

gis shapefile geodjango

8
推荐指数
1
解决办法
4611
查看次数

列名称的问题在PostgreSQL中包含冒号

我从OSM下载了形状数据.我已经将Shapefile中的数据导入PostgreSQL而没有任何问题,但是当我执行select语句时出现错误.

Select addr:city From location;

Error: syntax error at or near ":"
Run Code Online (Sandbox Code Playgroud)

问题是因为列名包含冒号.任何人都可以帮我解决这个问题吗?我应该在导入过程中拒绝这个shapefile吗?shapefile是否正常?

postgresql postgis geospatial shapefile geoserver

8
推荐指数
1
解决办法
3495
查看次数

如何转换Shapefile的坐标?

我正在尝试将邻居数据存入我的应用程序,并且我遇到了我正在使用的数据的问题,我从这里得到了这些数据.

此文件包含具有旧金山邻域的shapefile.我正在运行Ruby on Rails框架,我目前正在使用GeoRuby来解析shapefile.

代码如下所示:

def self.run_import
  shpfile = '/path/to/realtor_neighborhoods/realtor_neighborhoods'
  ShpFile.open(shpfile) do |shp|
    shp.each do |shape|
      # This gets the first (and only) Polygon from each MultiPolygon
      polygon = shape.geometry.geometries.first 
      puts polygon.inspect
    end
  end
end
Run Code Online (Sandbox Code Playgroud)

代码能够解析文件,但我无法理解解释的坐标.对于有效的纬度和经度,当我期望坐标在-180和180之间时,所有点都具有数百万的值.看一个示例点:

<GeoRuby::SimpleFeatures::Point:0x00000104566a08 @srid=4326, @with_z=false, \
   @with_m=false, @x=6015402.9999795845, @y=2114960.4999904726, @z=0.0, @m=0.0>,
Run Code Online (Sandbox Code Playgroud)

这些坐标值的格式是什么?如何将它们转换为对我有意义的值?(即基于SRID 4326 <=> WGS84空间参考系统的纬度/经度)

先感谢您!

ruby gis postgis ruby-on-rails shapefile

8
推荐指数
1
解决办法
3363
查看次数

在R中将Shapefile转换为光栅?

我有一个从worldwildlife.org下载的shapefile,用于世界地球生态区.该文件可以在这里加载:http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

它是一个标准的形状文件,我想用它做两件事.首先:从我的本地目录中获取shapefile并将其剪辑到北美东部(ext = extent(-95,-50,24,63))

# Read shapefile using package "maptools"
eco_shp <- readShapeLines("F:/01_2013/Ecoregions/Global/wwf_terr_ecos.shp", 
                          proj4string=CRS("+proj=utm +zone=33 +datum=WGS84")) 


# Set the desired extent for the final raster using package "raster" 
ext <- extent(-95, -50, 24, 63)
Run Code Online (Sandbox Code Playgroud)

我确信我必须在"raster"包中使用rasterize功能,但我仍然无法使其正常工作.我很感激有关如何做到这一点的任何建议.

r raster data-conversion shapefile rasterizing

8
推荐指数
1
解决办法
9672
查看次数

shapefile和matplotlib:shapefile坐标的多边形集合

我正在尝试使用python中的matplotlib在世界地图上绘制填充的多边形国家.

我有一个shapefile,每个国家都有国家边界坐标.现在,我想将这些坐标(对于每个国家/地区)转换为带有matplotlib的多边形.不使用Basemap.不幸的是,零件交叉或重叠.是否有工作,可能使用点到点的距离..或重新排序? 在此输入图像描述

python polygon matplotlib geospatial shapefile

8
推荐指数
1
解决办法
2万
查看次数

如何在Java中运行GDAL(ogr2ogr)将Shapefile转换为GeoJSON

我是编程的初学者,在尝试使用外部库时会感到非常困惑.我有shapefile中的地图,我使用Mapshaper.org网站转换为GeoJSON,然后才能从我的Java应用程序中读取地图.

我希望用户能够直接导入shapefile,因此我发现了具有Java绑定的 GDAL ogr2ogr,并且可以使用单个命令进行此转换,例如:

$ ogr2ogr -f GeoJSON -t_srs crs:84 [name].geojson [name].shp

但我不知道如何从我的程序中使用它.

有人能为我提供一个如何实现这个的例子吗?我该如何将GDAL导入我的项目?我正在使用LibGDX,如果有任何帮助的话.

最后,我需要的是一种直接在我的应用程序上将shapefile转换为GeoJSON的方法,这样用户就可以自己添加新的地图,所以任何建议都是受欢迎的.

非常感谢!

java json shapefile gdal

8
推荐指数
1
解决办法
5546
查看次数

更改Shapefile的投影

我试图更改或指定一个德国-shape文件的投影NA+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0,但不知何故,它不能很好地工作.

可重复的示例:Shapefile和其他文件可以在这里下载:

我尝试的是以下内容:

library(maptools)
library(sp)
library(rgeos)
library(rgdal)
projection.x <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0")
mapG <- readShapePoly("vg2500_lan.shp", verbose=TRUE, proj4string=projection.x)
summary(mapG)
mapG <- spTransform(mapG, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
Run Code Online (Sandbox Code Playgroud)

所以,问题是我无法在地图上绘制我的观察结果.见下文

在此输入图像描述 使用包装中的geocode功能检测到ponits ggmap. 在此输入图像描述 任何想法如何改变shapefile的投影或谷歌坐标的投影将非常感激!

r projection shapefile r-sp

8
推荐指数
1
解决办法
1万
查看次数

Python:如何从文本文件创建点形状文件

我正在编写一个 python 代码来读取多边形形状文件中的并将它们保存在一个点形状文件中。所以首先我制作了一个文本文件并将点的 (x,y) 存储在该 .txt 文件中。然后我试图从文本文件中制作一个点形状文件,但它给出了一个错误。
这是代码(只是最后一部分):

creat point shape-file from text file 
import fileinput
import string
import os
env.overwriteOutput=True
outpath="C:/roadpl"
newfc="newpoint.shp" 
arcpy.CreateFeatureclass_management(outpath, newfc, "Point")
infile="C:/roadpl/roadL5.txt"
cursor=arcpy.da.InsertCursor(newfc, ["SHAPE@"])
array=arcpy.Array()
for line in fileinput.input(infile):
    X, Y=string.split(line, " ")
    array.add(arcpy.Point(X,Y))
cursor.insertRow([arcpy.Point(array)])
fileinput.close()
del cursor
Run Code Online (Sandbox Code Playgroud)

这是错误:

Traceback (most recent call last):
  File "C:\Lab5\P_Code_L5", line 49, in <module>
    point.X, point.Y  = line.split()
  File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\arcobjects\_base.py", line 87, in _set
    return setattr(self._arc_object, attr_name, cval(val))
RuntimeError: …
Run Code Online (Sandbox Code Playgroud)

python polygon arcgis shapefile arcpy

7
推荐指数
1
解决办法
6178
查看次数

如何在Python中使用uber h3将shapefile/geojson转换为六边形?

我想在我的地理地图上创建六边形,并希望保留 shapefile/geojson 指定的数字边界。如何使用 uber 的 h3 python 库来做到这一点?

我对 shapefile 或任何其他地理数据结构不熟悉。我对 python 最满意。

python shapefile geojson h3

7
推荐指数
2
解决办法
1万
查看次数