GIS stuff

I have incidentally found these two nice things that might interest the GIS community: A huge shapefile with all the municipalities of France (+30,000 polygons), and their names and population. Very accurate, a vertex every few meters. About 7 million edges in one SHP file: http://sandre.eaufrance.fr/ftp/incoming/ONEMA/ADCO_France/ (EPSG:2154, I think) NTv2 grid files to transform from/to the OSGB36 and ETRS89 datums (United Kingdom): http://lists.maptools.org/pipermail/proj/2006-November/002658.html I'm almost sure those grids have been derived by intensively querying Ordnance Survey's conversion service: http://gps.ordnancesurvey.co.uk/etrs89geo_natgrid.asp Have not tested, but I also believe that the accuracy of that grid is not significantly worse than that of the web service itself, so that gives you an error below 10 cm (if I understood it well) for the whole country.

Puesto que gvSIG 1.9 no permite utilizar la transformación Molodensky-Badekas de 10 parámetros , una alternativa consiste en crear una rejilla o grilla de transformación en formato NTv2 utilizando los resultados de la transformación en un gran número de puntos de control repartidos por todo el territorio (en este caso, Venezuela). Hermógenes Suárez y José Nicolás Pérez Duin han publicado en la lista de usuarios de gvSIG una hoja de cálculo con las fórmulas y los parámetros necesarios para realizar la transformación. Para realizar la rejilla de transformación, partimos del rectángulo mínimo que contiene a todo el territorio de Venezuela: Longitud: [-73.5, -59.75] Latitud: [0.5, 12.25] La distancia entre los nodos de la rejilla debe, por un lado, ser lo suficientemente pequeña como para proporcionar una precisión adecuada, pero por otro lado no debemos generar una rejilla con un número desorbitado de nodos. En el primer intento, generaremos una rejilla con una distancia entre nodos de ...

Con la rejilla de transformación que proporciona el IGN, es posible reproyectar de EPSG:23029/30/31 a EPSG:4326 y viceversa. En este vídeo se hace lo siguiente: - Crear una vista en EPSG:23031 y cargar una capa WMS de alta definición del ICC. - Calcar un edificio y almacenarlo como SHP. - Crear una nueva vista en EPSG:4326 y añadir el SHP anterior, reproyectando al vuelo con la rejilla del IGN que viene incorporada en gvSIG. - Exportar la capa para tener un SHP en EPSG:4326 (y se añade a la vista). - Exportar el SHP en EPSG:4326 al formato KML. - Abrir el KML en Google Earth y comprobar que coincide exactamente con el edificio que habíamos calcado al inicio.