标签: shapefile

我在目录中有一个形状文件列表，我正在尝试使用arcpy.

有任何想法吗？谢谢。

无法删除问题.请参考问题:根据字典值使用Basemap对国家/地区进行阴影状态

我想在每个墨西哥州绘制数据(特定年份的病人数).我正在使用jupyter笔记本.到目前为止,我已经看到了几个选项和教程,但似乎没有一个似乎明确解释如何绘制一个国家的地图.下面我解释一些我见过的选项/教程以及为什么它们没有工作(我只是认为教程不是很直接):

1. Bokeh(http://bokeh.pydata.org/en/latest/docs/gallery/texas.html).在教程中,由于us_counties在bokeh.sampledata中,因此绘制了德州州.但是我没有在抽样数据中找到其他国家.

2. mpl_toolkits.basemap(http://www.geophysique.be/2011/01/27/matplotlib-basemap-tutorial-07-shapefiles-unleached/).虽然我能够导入shapefile,但我无法运行from shapefile import ShapeFile(ImportError:无法导入名称ShapeFile).此外,我无法下载dbflib库.

3. 文森特(为什么Python的文森特地图visuzalization不从数据帧映射数据？在说教程出现没有图像(即使我使用的命令)当我运行从答案的代码vincent.core.initialize_notebook()).

4. Plotly(https://plot.ly/python/choropleth-maps/).本教程绘制了美国从csv表导入信息的地图(没有其他可用国家的信息).如果想要策划另一个国家,是否可以制作表格？

探索了这四个选项我发现教程不是很清楚或不容易理解.我发现很难相信在python中绘制一个国家的地图很困难.我认为必须有一个比过去教程中解释的更简单的方法.

问题是:用python绘制某个国家(任何)地图的最简单(希望是简单的)方法是什么？

我安装了以下软件包:matplotlib,pyshp,mpl_toolkits.basemap,bokeh,pandas,numpy.我还从http://www.gadm.org/下载了墨西哥的地图

提前致谢.

我试图用密歇根州的县数据修改这个例子.简而言之,它是有效的,但它似乎在绘制县的过程中在这里和那里增加了一些额外的形状.我猜测在某些情况下(有岛屿的县),岛屿部分需要列为一个单独的"县",但我不确定另一个案例,例如在较低的韦恩县正确的部分国家.

这是我目前拥有的图片:

这是我到目前为止所做的:

1. 从Bokeh的样本县数据中获取县数据只是为了获得每个州号的州缩写(我的第二个,主数据源只有州号).对于这个例子,我将通过仅过滤状态号26来简化它.
2. 从美国人口普查网站获取各州的州立坐标('500k'文件).
3. 使用以下代码生成密歇根州的"交互式"地图.

注意:要pip install shapefile(真的是pyshp),我想我必须从这里下载.whl文件,然后执行pip install [path to .whl file].

import pandas as pd
import numpy as np
import shapefile
from bokeh.models import HoverTool, ColumnDataSource
from bokeh.palettes import Viridis6
from bokeh.plotting import figure, show, output_notebook
shpfile=r'Path\500K_US_Counties\cb_2015_us_county_500k.shp'
sf = shapefile.Reader(shpfile)
shapes = sf.shapes()

#Here are the rows from the shape file (plus lat/long coordinates)
rows=[]
lenrow=[]
for i,j in zip(sf.shapeRecords(),sf.shapes()):
    rows.append(i.record+[j.points])
    if len(i.record+[j.points])!=10:
           print("Found record with irrular number of …

我对R来说有点新手,所以这对你来说似乎很容易.我有一个加载在R中的shapefile,里面有几个多边形,我只想绘制第一个多边形.我尝试了几个命令,但它没有用.

我的目标是用我的shapefile的第一个元素裁剪一个栅格,再次尝试了几个东西,但它不起作用..

我希望我很清楚,

非常感谢!

裘德

我试图从 shp 文件中绘制俄罗斯，但俄罗斯总是分为两部分。我怎样才能让俄罗斯一气呵成？

我尝试了几个 shp 文件，例如 http://www.naturalearthdata.com/downloads/10m-culture-vectors/ Admin 1 – States, Provinces ；下载无大湖 (14.11 MB) 3.0.0 版

 shp <- readOGR("ne_10m_admin_1_states_provinces_lakes.shp")

子集到俄罗斯

names(shp)
rus <- shp[shp$admin == "Russia" , ]

x11()
plot(rus)

我正在尝试在shapefile中创建一个网格，诸如此类。但是，我无法生成这样的网格。我想知道是否有人对如何实现这一目标有想法。

这是我的代码-

WWWL.Shape<- readOGR("E:/Juan Arango", "WWL_Commerce_OK")
WWWL.Shape
plot(WWWL.Shape)
proj4string(WWWL.Shape)

bb <- bbox(WWWL.Shape)
cs <- c(3.28084, 3.28084)*6000  # cell size 
cc <- bb[, 1] + (cs/2)  # cell offset
cd <- ceiling(diff(t(bb))/cs)  # number of cells per direction
grd <- GridTopology(cellcentre.offset=cc, cellsize=cs, cells.dim=cd)
grd
sp_grd <- SpatialGridDataFrame(grd,
                               data=data.frame(id=1:prod(cd)),
                               proj4string=CRS(proj4string(WWWL.Shape)))
plot(sp_grd)

输出 WWL.Shape

class       : SpatialPolygonsDataFrame 
features    : 1 
extent      : 334367, 334498.7, 4088915, 4089057  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=15 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
variables   : 1 …

在学校，我必须创建一个API，在其中我需要一个shapefile内的几何图形。我使用geotools librairie打开shapefile。这是我的代码：

public class main {

    public static void main(String[] args) throws Exception{
        String filepath = "SHP_test/testSHPindex.shp";
        File shp = new File(filepath);

        Map<String, URL> map = new HashMap<String, URL>();
        map.put( "url", shp.toURI().toURL() );
        DataStore dataStore = DataStoreFinder.getDataStore(map);


        SimpleFeatureSource featureSource = dataStore.getFeatureSource(dataStore.getTypeNames()[0]);        
        SimpleFeatureCollection collection = featureSource.getFeatures();


}

我得到这个错误：

Exception in thread "main" java.lang.IllegalArgumentException: org.geotools.data.DataStoreFactorySpi is not an ImageIO SPI class
    at java.desktop/javax.imageio.spi.ServiceRegistry.checkClassAllowed(ServiceRegistry.java:745)
    at java.desktop/javax.imageio.spi.ServiceRegistry.<init>(ServiceRegistry.java:140)
    at org.geotools.factory.FactoryRegistry.<init>(FactoryRegistry.java:155)
    at org.geotools.factory.FactoryCreator.<init>(FactoryCreator.java:91)
    at org.geotools.data.DataStoreFinder.getServiceRegistry(DataStoreFinder.java:128)
    at org.geotools.data.DataStoreFinder.getAvailableDataStores(DataStoreFinder.java:115)
    at org.geotools.data.DataStoreFinder.getDataStore(DataStoreFinder.java:87)
    at org.ensg.projetperso.llb.maven.api.spatialindexingtest.testshapefile.main.main(main.java:30)

数据存储区不希望被创建，我不知道它是从哪里来的……您能帮我吗？提前致谢 ..

劳尔

有人有过使用 ISOBUS Iso-XML (ISO 11783-10) 的经验吗？

我想将一些 geojson 和/或 shapefile (.shp) 转换为这种格式，但找不到任何信息。

我有三个相互重叠的多边形 shapefile。让我们称他们为：

• file_one.shp（多边形名称为 1）
• file_two.shp（多边形名称为 2）
• file_three.shp（多边形名称为 3）

我想将它们结合起来并保持这样的值。

请问如何在Python中实现结果（如图所示）？

谢谢！

我有一个包含形状文件的文件夹；更正式地说，我拥有的文件是：“.dbf”、“.prj”、“.shp”和“.shx”。现在，我通过以下方式读取形状文件：

shapefile_path = r".\canada.shp"
canada = geopandas.read_file(shapefile_path)

但是其他文件呢？

任何帮助深表感谢！！