Geopandas의 셰이프 파일 좌표를 변환하십시오.

Question

공간 분석에서 새로운데, 어디서나이 답을 찾을 수 없습니다.

OSN 좌표의 CRS 좌표, 위도 및 경도, 런던의 자치도 모양 파일에있는 게시물 코드 목록이 있는데,이 두 파일을 함께 매핑하고 싶습니다. 이 geopandas

borough = gpd.read_file('London_Borough_Excluding_MHW.shp') 
borough.head() 
borough.head() 

NAME GSS_CODE geometry 
0 Kingston upon Thames E09000021 POLYGON ((516401.6 160201.8, 516407.3 160210.5... 
1 Croydon E09000008 POLYGON ((535009.2 159504.7, 535005.5 159502, ... 
2 Bromley E09000006 POLYGON ((540373.6 157530.4, 540361.2 157551.9... 
3 Hounslow E09000018 POLYGON ((521975.8 178100, 521967.7 178096.8, ... 
4 Ealing E09000009 POLYGON ((510253.5 182881.6, 510249.9 182886, ... 

우리는 다각형의 좌표는 우편 번호의와는 다른 방법을 볼 수 있습니다에서 읽은 모양 파일이있는 동안 이것은 우편 번호

london_post_codes.head() 
Out[81]: 
postcode latitude longitude 
0 WD6 1GS 51.658021 -0.255663 
1 WD17 1LA 51.660366 -0.397525 
2 WC2N 6LE 51.509413 -0.121676 
3 WC2N 6NA 51.508363 -0.124454 
4 WC2N 6ND 51.508216 -0.123829 

의 머리입니다. 나는 내가 그들을 함께 그릴 때 나는

fig, ax = plt.subplots() 
borough.plot(ax = ax) 
borough = gpd.read_file('statistical-gis-boundaries-london/ESRI/London_Borough_Excluding_MHW.shp') 
london_post_codes.plot(kind='scatter',s=10, x='longitude', y='latitude',ax=ax) 

어떤 제안을 얻을?

Answer 1

해결 방법은 간단히 CRS (Coordinate Reference System)를 변경하는 것입니다. 이것들은 epsg라고 불리는 코드에 의해 진행됩니다. 위도/경도 CRS의 코드는 epsg = 4326입니다. 그러므로

borough = gpd.read_file('London_Borough_Excluding_MHW.shp') 
borough = borough.to_crs(epsg=4326) 

그리고 나서 나머지.