Uli, Unless your neighbors are very localized (in an undistorted projection), I would recommend using some geodesic methods (in ArcMap or OGR).
For ArcMap, you can beat this free add-on (Tools for Graphics and Shapes) which does geodesic lengths, areas, and azimuths no matter the input projection: http://www.jennessent.com/arcgis/shapes_graphics.htm (very good documentation) For OGR, I would recommend running the lat/lon pairs through or Haversine or Vincenty's equations. examples: Haversine: http://stackoverflow.com/questions/11703407/how-to-calculate-distance-between-points-on-python-side-of-my-appplication-in-wa Vincenty's: https://code.google.com/p/pyroms/source/browse/trunk/pyroms/greatcircle.py?r=58 Good luck, Trent On Wed, Mar 6, 2013 at 5:33 PM, Uli Strötz <[email protected]> wrote: > Thanks for your answers. I tried number two and the files are projected in > Albers (meters) (see attachment). The results are still the same as being > unprojected. Is there a way to tell GDAL to only use the Projected > Coordinate System instead of the CRS? > > > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
