簡體   English   中英

了解我的shapefile的坐標點(不是緯度/經度)以及為什么空間連接不起作用

[英]Understanding my shapefile's coordinate points (not lat/lon) and why spatial joins own't work

在此處下載了墨西哥城郵政編碼 shapefile。

我正在嘗試進行空間連接,以查看哪些郵政編碼與我感興趣的特定區域相交。

但是當我加載郵政編碼時:

codigos = gpd.read_file('/Users/jaehakyu/Documents/Shapefiles/Mexico Postal Codes/CP_CdMx/CP_09CDMX_v4.shp')
codigos.head()

我注意到這些郵政編碼在多邊形中產生了奇怪的輸出。 它們不包含緯度/經度坐標。 相反,它們看起來像這樣:

POLYGON((2794650.259300001 819469.0967999995,2794601.180299997 819320.3825999983,2794595.038099997 819303.2646000013,2794581.846100003 819266.4985000007,2794552.451700002 819177.5449000001,2794539.257799998 819142.9415000007,2794490.202399999 819013.9316000007,2794461.8829 819006.0326999985,2794333.362800002 818970.185899999,2794287.639799997 818974.6688000001,2794223.908 818980.9173000008,2794210.812299997 818904.4486999999,2794187.823399998 818901.436999999,2794164.019500002 818896.5203​​000009,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 81903 3.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.170 0000018,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.9071 819890.1013999991,2794764.265699998 819868.1314000003,2794762.4494 819852.9858999997,2794736.8922 819766.9114000015,2794727.679099999 819735.4631999992,2794672.272100002 819545.7978999987,2794650.259300001 819469.0967999995))

我仍然嘗試進行空間連接,但是沒有運氣。 我也嘗試過使用CartoDB映射這些多邊形,但是它們沒有顯示出來。 我不確定自己要處理的內容,希望了解一些幫助:

(1)這里的坐標點是什么。 (2)如何使它們達到經/緯度,以便空間連接可以工作

謝謝! 順便說一下,我正在使用geopandas並在Python中進行勻稱處理。

也許不用說,但是要在空間上合並兩個數據集,它們必須處於同一坐標系中。 緯度/經度是某些地理坐標系中的坐標(並且不止一個,因此請確保您知道哪個!)。

Shapefile可以包含各種坐標系中的坐標。 通常,存在一個與shapefile一起分發的.prj文件,該文件描述了坐標系信息。

查看您的數據,我想它是在某些投影坐標系中的東西方向(以米為單位)。 因此,在將數據與地理坐標系中的數據進行空間連接之前,您首先需要取消投影數據(從投影坐標系或平面坐標系到球面坐標系),然后可能需要進行數據移動 (調整不同地理區域之間的差異)坐標系)。

我不了解geopandas或shapep,但是我知道proj4庫具有良好的Python綁定並且可以執行這些操作。 但是,無論哪種方式,您都需要.prj文件中的數據。

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM