<HTML dir=ltr><HEAD><TITLE>Re: [Proj] Swiss CORINE: Bessel 1841 Hotine Oblique MercatorAzimuth Natural Origin?</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.6000.16850" name=GENERATOR></HEAD>
<BODY>
<DIV id=idOWAReplyText61119 dir=ltr>
<DIV dir=ltr><FONT face="Times New Roman" color=#000000 size=2>I believe the current datums are: "CH1903+" and "CHTRS95."</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>See my August 2001 column on Switzerland at: <A href="http://www.ASPRS.org/resources/grids">www.ASPRS.org/resources/grids</A> .</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>C. Mugnier</FONT></DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>From:</B> proj-bounces@lists.maptools.org on behalf of Markus Neteler<BR><B>Sent:</B> Fri 10-Jul-09 09:18<BR><B>To:</B> Ziegler Stefan<BR><B>Cc:</B> proj<BR><B>Subject:</B> Re: [Proj] Swiss CORINE: Bessel 1841 Hotine Oblique MercatorAzimuth Natural Origin?<BR></FONT><BR></DIV>
<DIV>
<P><FONT size=2>Hi Stefan,<BR><BR>I should have asked earlier here :) Thanks a lot.<BR>However, the map seems to be shifted for 200-300 meters.<BR>If there any datum which I should specify? This reminds me<BR>of<BR><A href="http://trac.osgeo.org/grass/ticket/35">http://trac.osgeo.org/grass/ticket/35</A><BR><BR>Best regards<BR>Markus Neteler<BR><BR>On Fri, Jul 10, 2009 at 4:07 PM, Ziegler Stefan<Stefan.Ziegler@bd.so.ch> wrote:<BR>> Hi Markus<BR>> try this:<BR>> ogr2ogr shp_dir/ -s_srs epsg:21781 -t_srs epsg:3035 clcbn85p.shp<BR>> epsg:21781 is forcing ogr2ogr to use the real swiss prj-definition. Forget<BR>> the strange definition in the delivered prj-file.<BR>> regards<BR>> Stefan<BR>>> -----Ursprüngliche Nachricht-----<BR>>> Von: Markus Neteler [<A href="mailto:neteler@osgeo.org">mailto:neteler@osgeo.org</A>]<BR>>> Gesendet am: Freitag, 10. Juli 2009 14:59<BR>>> An: PROJ.4 and general Projections Discussions<BR>>> Betreff: Re: [Proj] Swiss CORINE: Bessel 1841 Hotine Oblique Mercator<BR>>> Azimuth Natural Origin?<BR>>><BR>>> Dear list members,<BR>>><BR>>> I was earlier sending this message to the gdal-dev list without<BR>>> answer. trying here now:<BR>>><BR>>> On Mon, Jun 8, 2009 at 5:24 PM, Markus<BR>>> Neteler<neteler@osgeo.org> wrote:<BR>>> > hi,<BR>>> ><BR>>> > I am trying to reproject the Swiss CORINE data CLC90 [1] into LAEA.<BR>>> > The .prj file is attached. Unfortunately the transformation fails:<BR>>> ><BR>>> > ogr2ogr -t_srs eu_etrs89_laea.prj /tmp/clc90_eu_laea.shp<BR>>> /geodata/clcbn85p.shp<BR>>> > ERROR 6: Failed to initialize PROJ.4 with `+proj=omerc<BR>>> > +lat_0=46.95240555555556 +lonc=7.439583333333333 +alpha=90 +k=1<BR>>> > +x_0=-9419820.590700001 +y_0=200000 +ellps=bessel +units=m<BR>>> +no_defs '.<BR>>> > lat_0 = 0 or 90 or alpha = 90<BR>>> > Failed to create coordinate transformation between the<BR>>> > following coordinate systems. This may be because they<BR>>> > are not transformable, or because projection services<BR>>> > (PROJ.4 DLL/.so) could not be loaded.<BR>>> > Source:<BR>>> > PROJCS["Bessel_1841_Hotine_Oblique_Mercator_Azimuth_Natural_Origin",<BR>>> > GEOGCS["GCS_Bessel_1841",<BR>>> > DATUM["Not_specified_based_on_Bessel_1841_ellipsoid",<BR>>> > SPHEROID["Bessel_1841",6377397.155,299.1528128]],<BR>>> > PRIMEM["Greenwich",0.0],<BR>>> > UNIT["Degree",0.0174532925199433]],<BR>>> > PROJECTION["Hotine_Oblique_Mercator"],<BR>>> > PARAMETER["False_Easting",-9419820.5907],<BR>>> > PARAMETER["False_Northing",200000.0],<BR>>> > PARAMETER["Scale_Factor",1.0],<BR>>> > PARAMETER["Azimuth",90.0],<BR>>> > PARAMETER["Longitude_Of_Center",7.439583333333333],<BR>>> > PARAMETER["Latitude_Of_Center",46.95240555555556],<BR>>> > UNIT["Meter",1.0]]<BR>>> > Target:<BR>>> > PROJCS["ETRS89 / ETRS-LAEA",<BR>>> > GEOGCS["ETRS89",<BR>>> > DATUM["European_Terrestrial_Reference_System_1989",<BR>>> > SPHEROID["GRS 1980",6378137,298.2572221010042,<BR>>> > AUTHORITY["EPSG","7019"]],<BR>>> > AUTHORITY["EPSG","6258"]],<BR>>> > PRIMEM["Greenwich",0],<BR>>> > UNIT["degree",0.0174532925199433],<BR>>> > AUTHORITY["EPSG","4258"]],<BR>>> > UNIT["metre",1,<BR>>> > AUTHORITY["EPSG","9001"]],<BR>>> > AUTHORITY["EPSG","3035"]]<BR>>> ><BR>>> > Suggestions welcome to reproject this data set...<BR>>> ><BR>>> > thanks<BR>>> > Markus<BR>>> ><BR>>> > [1]<BR>>> <A href="http://www.bfs.admin.ch/bfs/portal/de/index/dienstleistungen/g">http://www.bfs.admin.ch/bfs/portal/de/index/dienstleistungen/g</A><BR>> eostat/datenbeschreibung/corine_land_cover.html<BR>>> Data: -> GIS-Daten CORINE Land Cover Schweiz 1990 (ZIP file)<BR>>><BR>> Is there a way to enable PROJ4 to support this projection?<BR>> thank in advance,<BR>> Markus<BR>> _______________________________________________<BR>> Proj mailing list<BR>> Proj@lists.maptools.org<BR>> <A href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</A><BR>_______________________________________________<BR>Proj mailing list<BR>Proj@lists.maptools.org<BR><A href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</A><BR></FONT></P></DIV></BODY></HTML>