Непонятно, как работать с «urn: ogc: def: crs: OGC: 1.3: CRS84» в геопандах - PullRequest
Новая
       26

Непонятно, как работать с «urn: ogc: def: crs: OGC: 1.3: CRS84» в геопандах

1 голос
/ 26 октября 2019

У меня есть этот файл геоджона 

{
"type": "FeatureCollection",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:OGC:1.3:CRS84" } },
"features": [
{ "type": "Feature", "properties": { "visit_date": "2013-03-27Z", "name": "Mayi-Tatu", "n_workers": 150.0, "mineral": "Gold" }, "geometry": { "type": "Point", "coordinates": [ 29.66033, 1.01089 ] } },
{ "type": "Feature", "properties": { "visit_date": "2013-03-27Z", "name": "Mabanga", "n_workers": 115.0, "mineral": "Gold" }, "geometry": { "type": "Point", "coordinates": [ 29.65862, 1.00308 ] } },
{ "type": "Feature", "properties": { "visit_date": "2013-03-27Z", "name": "Molende", "n_workers": 130.0, "mineral": "Gold" }, "geometry": { "type": "Point", "coordinates": [ 29.65629, 0.98563 ] } },
...
{ "type": "Feature", "properties": { "visit_date": "2017-08-31Z", "name": "Kambasha", "n_workers": 37.0, "mineral": "Cassiterite" }, "geometry": { "type": "Point", "coordinates": [ 29.05973167, -2.25938167 ] } }
]
}

Я читаю этот файл со следующим кодом: 

filename = "ipis_cod_mines.geojson"
df_congomines_crs84_geo = gpd.read_file(filename)

Но когда я проверяю свойство crs в df_congomines_crs84_geo, 

df_congomines_crs84_geo.crs

Я получил "{'init': 'epsg: 4326'}", я не понимаю, почему я не получаю правильные crs. (первый вопрос)

После того, как я прочитал другой набор данных для той же области (оба данных принадлежат Конго) 

df_countries_4326_geo = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))

Этот набор данных имеет crs, равный {'init': 'epsg:4326}. Когда я строю оба набора данных (без изменения crs), 

ax = congo_df.plot(alpha=0.5, color='brown', figsize=(11,4))
df_congomines_crs84_geo.plot(ax=ax, column='mineral')
plt.show()

я получил следующее изображение: Результат изображения

Почему оба изображения не перекрываются, если они принадлежатв ту же область ??? Как я могу это исправить ??? Связана ли эта проблема с зоной UTM ??? (второй вопрос)

...