Python 了解我的shapefile的坐标点(不是lat/lon)以及空间连接不起作用的原因

Python 了解我的shapefile的坐标点(不是lat/lon)以及空间连接不起作用的原因,python,gis,polygon,postgis,shapefile,Python,Gis,Polygon,Postgis,Shapefile,我下载了一个的形状文件 我正在尝试进行空间连接,以查看哪些邮政编码与我感兴趣的特定区域相交 但当我加载邮政编码时: codigos = gpd.read_file('/Users/jaehakyu/Documents/Shapefiles/Mexico Postal Codes/CP_CdMx/CP_09CDMX_v4.shp') codigos.head() 我注意到这些邮政编码的多边形中有奇怪的输出。它们不是由纬度/经度坐标组成的。相反,它们看起来是这样的: 多边形2794650.2593

我下载了一个的形状文件

我正在尝试进行空间连接,以查看哪些邮政编码与我感兴趣的特定区域相交

但当我加载邮政编码时:

codigos = gpd.read_file('/Users/jaehakyu/Documents/Shapefiles/Mexico Postal Codes/CP_CdMx/CP_09CDMX_v4.shp')
codigos.head()
我注意到这些邮政编码的多边形中有奇怪的输出。它们不是由纬度/经度坐标组成的。相反,它们看起来是这样的:

多边形2794650.259300001819469.096799995、2794601.180299977819320.3825999983、2794595.038099977819303.264600013、2794581.84610003 819266.498500007、2794552.451700002 819177.544900001、2794539.257799998819142.941500007、2794490.20239999819013.931600007、279461.8829819006.032699985、27933.362800002.8999、, 2794287.639799997 818974.6688000001, 2794223.908 818980.9173000008, 2794210.812299997 818904.4486999999, 2794187.823399998 818901.436999999, 2794164.019500002 818896.5203000009, 2794128.538500004 818896.0044, 2794122.643299997 818895.9186999984, 2794077.665700004 818894.2122999988, 2794066.248400003 818891.7342999987, 2794046.369499996 818886.5003999993, 2794006.863200001 818880.5144000016, 2793988.908100002 818882.2074999996, 2793965.905199997 818887.1242000014, 2793955.3869 818892.3366999999, 2793903.004799999 818932.752799999, 2793796.985100001 818913.691399999, 2793741.294 818906.6301999986, 2793596.247299999 818889.4554999992, 2793571.957000002 819033.999499999, 2793552.235299997 819034.7226999998, 2793527.926200002 819036.2868999988, 2793444.091899998 818973.02, 2793442.978699997 819118.5722999983, 2793437.997699998 819133.6957999989, 2793430.251699999 819144.4745000005, 2793414.9217 819147.9660999998, 2793334.550800003 819163.0062999986, 2793290.008699998 819170.9169000015, 2793213.590499997 819263.4613000005, 2793110.11 819265.5681000017, 2793117.498800002 819308.4505000003, 2793146.486199997 819475.7197000012, 2793186.218999997 819462.5135999992, 2793249.972900003 819441.3235000001, 2793345.965300001 819414.0322999991, 2793436.257799998 819398.3966999985, 2793517.804700002 819384.2754999995, 2793565.4133 819377.8266999982, 2793601.6972 819374.2901000008, 2793787.7729 819356.7437000014, 2793841.170500003 819343.5784999989, 2793848.679399997 819336.9222000018, 2793862.701300003 819415.7888999991, 2793896.867299996 819558.6755999997, 2793912.430100001 819623.6581000015, 2793934.3768 819672.3953000009, 2793959.285499997 819765.1700000018, 2793964.056900002 819787.6761000007, 2793984.413400002 819883.6871000007, 2793997.486699998 819948.0406999998, 2794004.954099998 819981.5982000008, 2794111.821699999 819973.7314999998, 2794230.936700001 819964.9633999988, 2794347.152999997 819955.9384000003, 2794343.185599998 819899.5188000016,2794458.9071819890.101399991192794764.26569999819868.1314000003、2794762.4494819852.985899997、2794736.8922819766.9114000015、2794727.6790999819735.463199992、2794672.272100002 819545.797899987、2794650.2593000019469.096799995

我试着做一个空间连接,但没有运气。我还尝试使用CartoDB映射这些多边形,但它们没有显示出来。我不确定我在处理什么,希望能得到一些帮助来理解:

1这里的坐标点是什么。 2如何将它们连接到lat/lon以便空间连接可以工作


谢谢大家!!顺便说一句,我正在使用geopandas和Python中的shapely进行这方面的工作。

也许不用说,但是要在空间上连接两个数据集,它们需要在同一坐标系中。纬度/经度是某些地理坐标系中的坐标,并且有多个坐标,因此请确保您知道哪个

形状文件可以包括各种坐标系中的坐标。通常有一个.prj文件,该文件与描述坐标系信息的shapefile一起分发

看看你的数据,我猜它是在东边和北边,以米为单位,在一些投影坐标系中。因此,在与地理坐标系中的数据进行空间连接之前,首先需要将数据从投影坐标系或平面坐标系取消投影到椭球坐标系,然后可能针对不同地理坐标系之间的差异执行基准移动调整


我不知道geopandas或shapely,但我知道proj4库有很好的Python绑定,可以执行这些操作。但无论哪种方式,您都需要.prj文件中的数据。

这就是我看到的:GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_19846378137298.257223563]、PRIMEM[Greenwich,0],UNIT[Degree,0.017453292519943295]]现在正试图阅读更多关于它的信息,但通常对地理空间数据非常陌生,并且会喜欢您可能拥有的任何指针!非常感谢。