Me han dado algunos archivos de forma SIN una proyección. Entonces necesito agregar la proyección para mostrarla en el mapa de Google. Encontré la información de proyección que necesito en http://spatialreference.org/ref/epsg/2066 .
Pregunta:
¿Cómo puedo agregar la proyección usando un método de código abierto (no tengo acceso al software de ESRI)?
Gracias
Chris
Respuestas:
ogr2ogr -a_srs EPSG:2066 output.shp input.shp
Lo bueno de ogr es que también puede usar un archivo PRJ existente como entrada:
ogr2ogr -a_srs input.prj output.shp input.shp
La proyección en el archivo .shp está en shapefilename.prj. Básicamente es un archivo de texto, con la información que encontraste. Eche un vistazo a algunos de sus archivos de formas y descubra cómo se escribe la información. O menos cumplido, si tiene otro archivo .shp en la misma proyección, simplemente cópielo y cambie el nombre según corresponda.
Cree un archivo .prj con el mismo nombre base que su archivo de forma (es decir, foo.prj), abra el archivo en un editor de texto y agregue el contenido de, en este caso, http://spatialreference.org/ref/epsg/2066/ prj /
Descargue el archivo .prj en el enlace que publicó y luego cámbiele el nombre para que coincida con su archivo de forma.
Encuentra adjunto un pequeño script que escribe la referencia de proyección desde spatialreference.org en un archivo .prj. Agrega un archivo de proyección a todos los archivos especificados en un directorio. Por ejemplo, todos los archivos de forma en el directorio E: \. Simplemente preocúpese por el código EPSG de la proyección que desea incrustar, la extensión de los archivos a los que desea agregar un archivo de proyección y el directorio donde están. Recorrerá recursivamente todos los subdirectorios, así que úselo con cuidado.
import os
def getWKT_PRJ (epsg_code):
import urllib.request, urllib.parse, urllib.error
# Access projection information
wkt = urllib.request.urlopen("http://spatialreference.org/ref/epsg/{0}/prettywkt/".format(epsg_code))
decoded = (wkt.read().decode('utf-8'))
# Remove spaces between charachters
remove_spaces = decoded.replace(" ","")
# Place all the text on one line
output = remove_spaces.replace("\n","")
return output
def referencer(folder_path, extension):
for path, subdirs, files in os.walk(folder_path):
for name in files:
file_extension = os.path.splitext(name)[-1]
if(extension in file_extension):
file_path = os.path.join(path,name)
file_name = os.path.splitext(file_path)[0]
prj = file_name + ".prj"
projection = open(prj,"w")
projection.write(epsg)
projection.close()
epsg = getWKT_PRJ("25831")
referencer('E:\Testfolder', '.shp')
Intente utilizar Geokettle http://www.spatialytics.org/ , lee archivos de forma y tiene la opción de agregar un SRS en sus transformaciones. De lo contrario, use ogr2ogr gui (solo Windows) para ver la cuadrícula