Как добавить проекцию (т.е. epsg = 102001) в DF Geopandas, который находится в файле esri.extra, а НЕ в файле proj_def.dat?

Question

У меня есть кадр данных Geopandas с CRS epsg: 4326, который я хотел бы преобразовать в CRC espg: 102001, который имеет определение proj.4 следующим образом:

proj4_102001 = '+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs' 

В моем /В каталоге usr / share / proj проекция находится в файле esri.extra, но не в файле proj_def.dat, как показано на снимке экрана ниже:

Снимок экрана 1) результатыкоманды Grep для поиска 102001 и 2) ls каталога / usr / share / proj

Согласно документации Geopandas , правильные способы определения проекцииследующим образом:

world = world.to_crs({'init': 'epsg:3395'})

Я пробовал "world = world.to_crs ({'init': 'esri.extra: 102001'})" и ряд других вариантов, но ни один из них не работал.Я пытался найти в Интернете решения, включая документацию Proj.4 , но не нашел ни одного.Пожалуйста, дайте мне знать, если вы знаете, как это сделать.Я использую Ubuntu (Ubuntu-18.04), Python (3.7.1), Geopandas (0.4.0) и Proj.4 (5.2.0).

Заранее спасибо.

Том

Answer 1

Вы должны предоставить полную строку проекции в функции to_crs.

world = world.to_crs('+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs')