Question

我有一些来自gps数据（Lat，Lon）的数据点以及我正在考虑的区域的相对shapefile。所以：

import pandas as pd
import shapefile
sf = shapefile.Reader('roadshape')
df = pd.read_csv('gpsdata.csv')

现在df包含我正在分析的数据是这样的：

 ID          x          y
3447    11.400427   48.816806
3448    11.400759   48.816772
3449    11.401424   48.816684
3450    11.401758   48.816631
3451    11.402090   48.816566
3452    11.402422   48.816490

我想为每个shapefile段分配。我正在尝试执行以下操作。我正在考虑边界框，我想尝试查看数据指向数据栏的边界框。

dfs = pd.DataFrame()
shapes = sf.shapes()
X =list()
Y=list()
for i in range(0,len(shapes)):
      X.append([shapes[i].bbox[0],shapes[i].bbox[2]])
      Y.append([shapes[i].bbox[1],shapes[i].bbox[3]])
dfs['X'] = X
dfs['Y'] = Y

现在我如何查看我的积分是哪个bbox？ dfs就是这样的地方

dfs = 
                 X                             Y
0   [10.9467244189, 10.9704393002]  [48.2671975178, 48.2697440003]
1   [11.5138847999, 11.5143541004]  [48.6497096997, 48.6515363002]
2   [11.4618209998, 11.4620896004]  [48.9305448001, 48.9307776004]
3   [10.6196591004, 10.6207268996]  [48.8635958001, 48.8665684003]
4   [10.652098, 10.6559025999]  [48.8005320998, 48.8042877999]
5   [11.1863882997, 11.1884544004]  [48.3726685999, 48.3738253996]
6   [11.1580075998, 11.1593822] [48.3785226999, 48.3791247996]
7   [11.1077987, 11.1112508996] [48.3829125003, 48.3830440999]
8   [11.0842697004, 11.0886483996]  [48.3840543003, 48.3879626001]
9   [11.0910959001, 11.0926532003]  [48.3903297003, 48.3916850002]
10  [11.4766434001, 11.4822778002]  [49.0389071001, 49.0399456003]
11  [11.7037148998, 11.7073818] [48.6927748996, 48.6961230001]
12  [11.7767894997, 11.7770049998]  [48.6279809001, 48.6279908997]

Answer 1

我会选择fiona，造型和内置csv！

假设我有两个文件：1）一个行shapefile，里面有几行记录。 2）具有三列的csv文件，即id，经度（x），纬度（y）

CSV文件内容（lat，long为UTM） -

id,lat,long

0,207726.012448,2733349.10914
0,197599.591396,2730510.17345
0,203187.5176,2736670.5686
0,207301.877268,2730639.81898
0,200929.610894,2726377.9799
0,204604.214301,2737608.342
0,203780.386032,2734372.2709
0,203077.172106,2731166.44271
0,202477.371994,2728622.46292
0,202249.861606,2734889.33996
0,201794.840831,2732159.21531

现在，下面的代码为line shapefile的每个记录生成bbox，并在打印完状态后检查该bbox文件在该bbox内的哪个点！

import fiona
from shapely import geometry
import csv
csv_file = open(r"C:\data_test.csv",'rb')
reader = csv.reader(csv_file)
pnts = []
for i in reader:
    pnts.append ((i[1],i[2]))
bboxes = []   
fiona_collection = fiona.open(r"C:\Rd.shp")
for i in fiona_collection:
    bboxes.append(geometry.asShape(i['geometry']).bounds)


for i in pnts:
    for j in bboxes:
        shape = geometry.box(*map(float,j))
        pnt = geometry.Point(*map(float,i))
        if shape.contains(pnt):
            print "bbox {0} contains {1} ".format(j,i)

打印 -

bbox (197306.68136754428, 2729718.6493367185, 197823.7504301681, 2732035.118737273) contains ('197599.591396', '2730510.17345') 
bbox (202799.29018572776, 2728807.5534000006, 204604.21430105247, 2737608.342) contains ('203187.5176', '2736670.5686') 
bbox (202799.29018572776, 2728807.5534000006, 204604.21430105247, 2737608.342) contains ('203780.386032', '2734372.2709') 
bbox (202799.29018572776, 2728807.5534000006, 204604.21430105247, 2737608.342) contains ('203077.172106', '2731166.44271') 
bbox (201237.0747999996, 2727558.2885, 203336.7932000002, 2728807.5534000006) contains ('202477.371994', '2728622.46292')

您可以在here 下载测试数据。

Python：为道路shapefile分配点

1 个答案: