7 votos

Carga sin éxito

Estoy intente cargar este archivo shapefile con Mapa base con el fin de trazar en Python:

http://geodata.gov.gr/en/dataset/periphereies-elladas/resource/7c80a2c1-93b7-4814-9fc4-245e775acaa6

Pero me sale el siguiente error:

ValueError: shapefile must have lat/lon vertices - it looks like this one has vertices
in map projection coordinates. You can convert the shapefile to geographic
coordinates using the shpproj utility from the shapelib tools

He intentado utilizar shpproj o ogr2ogr como otra respuesta en este sitio sugerido, pero para la vida de mí no puedo averiguar qué parámetros debo pasar a los programas, como tengo cero conocimiento acerca de los SIG en general. Así que si alguien de aquí me podría ayudar con el comando de conversión que sería genial.

Estoy en Ubuntu y el código que estoy usando es bastante simple, supongo:

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
map = Basemap()

map.drawmapboundary(fill_color='aqua')
map.fillcontinents(color='#ddaa66',lake_color='aqua')
map.drawcoastlines()

map.readshapefile('periphereies/periphereies', 'periphereies')

plt.show()

4voto

Nikola Puntos 21

Porque un archivo de proyección (*.prj) se une a lo shapefile (es decir, "GGRS87 / griego Grid" también conocido como EPSG:2100, fuente: EPSG geodésico parámetro registro), fácilmente podemos transformarla a WGS84 (EPSG:4326) utilizar ogr2ogr sin declarar la fuente SRS:

ogr2ogr -t_srs EPSG:4326 periphereies_WGS84.shp periphereies.shp

donde periphereies_WGS84.shp es el shapefile transformada.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X