pysal

1熱度

1回答

對於以統一的分配，我們必須使用PySal計算莫蘭的我。我發現了一個examplar代碼，並從我去，但我收到錯誤消息「‘文件’對象有沒有屬性‘by_col’」任何想法將不勝感激。數據集是紐約的文本文件，「T0_MINC」屬性是包含不同地區家庭收入中位數的列。該腳本被附加以下 import os import pysal import numpy as np os.chdir = ("C:\

1熱度

1回答

Python - 如何使用PySAL計算交互式空間自相關（Moran I）？

我有一個積分表my_table我在PostgreSQL數據庫中有geometry列和其他屬性。我有一些my_table的樣本數據如下（屬性my_table）。 id val1 1 72.54513286 2 73.67371014 3 74.204424 4 73.76017279 5 77.7912762 6 77.78789496 7 65.51822878 8 65.5182

5熱度

2回答

使用shapefile多邊形查找一階鄰居

我正在尋找一種有效的方法來查找給定多邊形的一階鄰居。我的數據格式爲shapefile。我的第一個想法是計算多邊形的質心的x和y座標，以找到鄰居的質心。 import pysal from pysal.common import * import pysal.weights import numpy as np from scipy import sparse,float32 impor

2熱度

2回答

Tokenising整數

我有一個包含的形式座標的文本文件： [-1.38795678, 54.90352965] [-3.2115, 55.95530556] [0.00315428, 51.50285246] 我希望能夠通過每個迭代協調，以檢查它是在多邊形（UK縣shape文件），但我不知道如何tokenise的號碼，這樣我可以有沿行代碼... for line in coordinates: for

2熱度

1回答

無法在pysal中設置距離單位

對於Pysal來說是新手，並且只是掌握加載，閱讀等操作。我有一些測試點數據，在Decimal Degrees（WGS84）中，並且想要使用get_point_dist函數。運行後，我得到正確的響應，但希望以米/ KM等的距離而不是DD。任何想法如何做到這一點，或者是答案，「投影輸入數據」？

6熱度

1回答

優化scipy最近鄰居搜索

我正在嘗試查找距離1 KM範圍內的所有最近鄰居。這裏是我的腳本來構建樹和搜索最近點， from pysal.cg.kdtree import KDTree def construct_tree(s): data_geopoints = [tuple(x) for x in s[['longitude','latitude']].to_records(index=False)]

2熱度

1回答

如何在Pysal中使用Shapefile

我希望在給定的多邊形中相交的點的結果，但是出現錯誤。我的代碼是： from pysal.cg.standalone import get_polygon_point_intersect poly=pysal.open('Busroute_buffer.shp') point=pysal.open('pmpml_24.shp') i=get_polygon_point_intersect(