<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<html lang="en">
<head><meta http-equiv="Content-Type" content="text/html;">
<title>Scalix message content</title>
<style type="text/css">
<!--
.Normal { text-align: left; text-indent:0pt; margin-top:0pt; margin-bottom:0pt; font-family: Arial; font-size: 12pt; }
.Default_Paragraph_Font {font-family: Arial; color: Black; background-color: white ; font-size: 12pt; font-weight: normal; font-style: normal; font-variant: normal; vertical-align: normal; text-decoration: none;}
div.ltTOCtitle { font-family: Verdana; font-size: 8pt; font-weight: bold; text-align: center;}
div.ltTOCl1 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:32pt; text-indent:-32pt;}
div.ltTOCl2 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:40pt; text-indent:-32pt;}
div.ltTOCl3 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:48pt; text-indent:-32pt;}
div.ltTOCl4 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:56pt; text-indent:-32pt;}
div.ltTOCl5 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:64pt; text-indent:-32pt;}
div.ltTOCl6 { font-family: Verdana; font-size: 8pt; font-weight: bold; margin-left:72pt; text-indent:-32pt;}
-->
</style>
</head>
<body>
<div class="Normal"><span style='font-size: 10pt; '>Hi Markus</span></div>
<br>
<div class="Normal"><span style='font-size: 10pt; '>please check in the proj4
epsg-file the definition of 21781. It needs the +towgs84 parameter:</span></div>
<br>
<div class="Normal"><span style='font-size: 10pt; '>+proj=somerc
+lat_0=46.95240555555556 +lon_0=7.439583333333333 +x_0=600000 +y_0=200000
+ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs
</span></div>
<br>
<div class="Normal"><span style='font-size: 10pt; '>The 200 meters shift sounds
a bit like a forgotten towgs84 parameter but that's just a guess.
</span><br><span style='font-size: 10pt; '>Does the Grass ticket still needs
some clarification?</span></div>
<br>
<div class="Normal"><span style='font-size: 10pt; '>regards</span><br><span style='font-size: 10pt; '>Stefan</span></div>
<br>
<div class="Normal"><span style='font-size: 10pt; '>> -----Ursprüngliche
Nachricht-----</span><br><span style='font-size: 10pt; '>> Von: Markus
Neteler
[mailto:neteler@osgeo.org]</span><br><span style='font-size: 10pt; '>>
Gesendet am: Freitag, 10. Juli 2009
16:18</span><br><span style='font-size: 10pt; '>> An: Ziegler
Stefan</span><br><span style='font-size: 10pt; '>> Cc:
proj</span><br><span style='font-size: 10pt; '>> Betreff: Re: Re: [Proj]
Swiss CORINE: Bessel 1841 Hotine
Oblique</span><br><span style='font-size: 10pt; '>> Mercator Azimuth Natural
Origin?</span><br><span style='font-size: 10pt; '>>
</span><br><span style='font-size: 10pt; '>> Hi
Stefan,</span><br><span style='font-size: 10pt; '>>
</span><br><span style='font-size: 10pt; '>> I should have asked earlier here
:) Thanks a lot.</span><br><span style='font-size: 10pt; '>> However, the map
seems to be shifted for 200-300
meters.</span><br><span style='font-size: 10pt; '>> If there any datum which
I should specify? This reminds me</span><br><span style='font-size: 10pt; '>>
of</span><br><span style='font-size: 10pt; '>>
http://trac.osgeo.org/grass/ticket/35</span><br><span style='font-size: 10pt; '>>
</span><br><span style='font-size: 10pt; '>> Best
regards</span><br><span style='font-size: 10pt; '>> Markus
Neteler</span><br><span style='font-size: 10pt; '>>
</span><br><span style='font-size: 10pt; '>> On Fri, Jul 10, 2009 at 4:07 PM,
Ziegler </span><br><span style='font-size: 10pt; '>>
Stefan<Stefan.Ziegler@bd.so.ch>
wrote:</span><br><span style='font-size: 10pt; '>> > Hi
Markus</span><br><span style='font-size: 10pt; '>> > try
this:</span><br><span style='font-size: 10pt; '>> > ogr2ogr shp_dir/
-s_srs epsg:21781 -t_srs epsg:3035
clcbn85p.shp</span><br><span style='font-size: 10pt; '>> > epsg:21781 is
forcing ogr2ogr to use the real swiss
</span><br><span style='font-size: 10pt; '>> prj-definition.
Forget</span><br><span style='font-size: 10pt; '>> > the strange
definition in the delivered
prj-file.</span><br><span style='font-size: 10pt; '>> >
regards</span><br><span style='font-size: 10pt; '>> >
Stefan</span><br><span style='font-size: 10pt; '>> >>
-----Ursprüngliche
Nachricht-----</span><br><span style='font-size: 10pt; '>> >> Von:
Markus Neteler
[mailto:neteler@osgeo.org]</span><br><span style='font-size: 10pt; '>>
>> Gesendet am: Freitag, 10. Juli 2009
14:59</span><br><span style='font-size: 10pt; '>> >> An: PROJ.4 and
general Projections Discussions</span><br><span style='font-size: 10pt; '>>
>> Betreff: Re: [Proj] Swiss CORINE: Bessel 1841 Hotine
</span><br><span style='font-size: 10pt; '>> Oblique
Mercator</span><br><span style='font-size: 10pt; '>> >> Azimuth Natural
Origin?</span><br><span style='font-size: 10pt; '>>
>></span><br><span style='font-size: 10pt; '>> >> Dear list
members,</span><br><span style='font-size: 10pt; '>>
>></span><br><span style='font-size: 10pt; '>> >> I was earlier
sending this message to the gdal-dev list
without</span><br><span style='font-size: 10pt; '>> >> answer. trying
here now:</span><br><span style='font-size: 10pt; '>>
>></span><br><span style='font-size: 10pt; '>> >> On Mon, Jun 8,
2009 at 5:24 PM, Markus</span><br><span style='font-size: 10pt; '>> >>
Neteler<neteler@osgeo.org>
wrote:</span><br><span style='font-size: 10pt; '>> >> >
hi,</span><br><span style='font-size: 10pt; '>> >>
></span><br><span style='font-size: 10pt; '>> >> > I am trying to
reproject the Swiss CORINE data CLC90 [1]
</span><br><span style='font-size: 10pt; '>> into
LAEA.</span><br><span style='font-size: 10pt; '>> >> > The .prj file
is attached. Unfortunately the </span><br><span style='font-size: 10pt; '>>
transformation fails:</span><br><span style='font-size: 10pt; '>> >>
></span><br><span style='font-size: 10pt; '>> >> > ogr2ogr -t_srs
eu_etrs89_laea.prj
/tmp/clc90_eu_laea.shp</span><br><span style='font-size: 10pt; '>> >>
/geodata/clcbn85p.shp</span><br><span style='font-size: 10pt; '>> >>
> ERROR 6: Failed to initialize PROJ.4 with
`+proj=omerc</span><br><span style='font-size: 10pt; '>> >> >
+lat_0=46.95240555555556 +lonc=7.439583333333333 +alpha=90
+k=1</span><br><span style='font-size: 10pt; '>> >> >
+x_0=-9419820.590700001 +y_0=200000 +ellps=bessel
+units=m</span><br><span style='font-size: 10pt; '>> >> +no_defs
'.</span><br><span style='font-size: 10pt; '>> >> > lat_0 = 0 or 90
or alpha = 90</span><br><span style='font-size: 10pt; '>> >> >
Failed to create coordinate transformation between
the</span><br><span style='font-size: 10pt; '>> >> > following
coordinate systems.  This may be because
they</span><br><span style='font-size: 10pt; '>> >> > are not
transformable, or because projection
services</span><br><span style='font-size: 10pt; '>> >> > (PROJ.4
DLL/.so) could not be loaded.</span><br><span style='font-size: 10pt; '>>
>> > Source:</span><br><span style='font-size: 10pt; '>> >>
> </span><br><span style='font-size: 10pt; '>>
PROJCS["Bessel_1841_Hotine_Oblique_Mercator_Azimuth_Natural_Origin",</span><br><span style='font-size: 10pt; '>>
>> >  
 GEOGCS["GCS_Bessel_1841",</span><br><span style='font-size: 10pt; '>>
>> >      
 DATUM["Not_specified_based_on_Bessel_1841_ellipsoid",</span><br><span style='font-size: 10pt; '>>
>> >          
 SPHEROID["Bessel_1841",6377397.155,299.1528128]],</span><br><span style='font-size: 10pt; '>>
>> >      
 PRIMEM["Greenwich",0.0],</span><br><span style='font-size: 10pt; '>>
>> >      
 UNIT["Degree",0.0174532925199433]],</span><br><span style='font-size: 10pt; '>>
>> >  
 PROJECTION["Hotine_Oblique_Mercator"],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["False_Easting",-9419820.5907],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["False_Northing",200000.0],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["Scale_Factor",1.0],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["Azimuth",90.0],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["Longitude_Of_Center",7.439583333333333],</span><br><span style='font-size: 10pt; '>>
>> >  
 PARAMETER["Latitude_Of_Center",46.95240555555556],</span><br><span style='font-size: 10pt; '>>
>> >  
 UNIT["Meter",1.0]]</span><br><span style='font-size: 10pt; '>> >>
> Target:</span><br><span style='font-size: 10pt; '>> >> >
PROJCS["ETRS89 / ETRS-LAEA",</span><br><span style='font-size: 10pt; '>>
>> >  
 GEOGCS["ETRS89",</span><br><span style='font-size: 10pt; '>> >>
>      
 DATUM["European_Terrestrial_Reference_System_1989",</span><br><span style='font-size: 10pt; '>>
>> >            SPHEROID["GRS
1980",6378137,298.2572221010042,</span><br><span style='font-size: 10pt; '>>
>> >              
 AUTHORITY["EPSG","7019"]],</span><br><span style='font-size: 10pt; '>>
>> >          
 AUTHORITY["EPSG","6258"]],</span><br><span style='font-size: 10pt; '>>
>> >      
 PRIMEM["Greenwich",0],</span><br><span style='font-size: 10pt; '>>
>> >      
 UNIT["degree",0.0174532925199433],</span><br><span style='font-size: 10pt; '>>
>> >      
 AUTHORITY["EPSG","4258"]],</span><br><span style='font-size: 10pt; '>>
>> >  
 UNIT["metre",1,</span><br><span style='font-size: 10pt; '>> >>
>      
 AUTHORITY["EPSG","9001"]],</span><br><span style='font-size: 10pt; '>>
>> >  
 AUTHORITY["EPSG","3035"]]</span><br><span style='font-size: 10pt; '>>
>> ></span><br><span style='font-size: 10pt; '>> >> >
Suggestions welcome to reproject this data
set...</span><br><span style='font-size: 10pt; '>> >>
></span><br><span style='font-size: 10pt; '>> >> >
thanks</span><br><span style='font-size: 10pt; '>> >> >
Markus</span><br><span style='font-size: 10pt; '>> >>
></span><br><span style='font-size: 10pt; '>> >> >
[1]</span><br><span style='font-size: 10pt; '>> >>
http://www.bfs.admin.ch/bfs/portal/de/index/dienstleistungen/g</span><br><span style='font-size: 10pt; '>>
>
eostat/datenbeschreibung/corine_land_cover.html</span><br><span style='font-size: 10pt; '>>
>>    Data: -> GIS-Daten CORINE Land Cover Schweiz 1990 (ZIP
file)</span><br><span style='font-size: 10pt; '>>
>></span><br><span style='font-size: 10pt; '>> > Is there a way to
enable PROJ4 to support this
projection?</span><br><span style='font-size: 10pt; '>> > thank in
advance,</span><br><span style='font-size: 10pt; '>> >
Markus</span><br><span style='font-size: 10pt; '>> >
_______________________________________________</span><br><span style='font-size: 10pt; '>>
> Proj mailing list</span><br><span style='font-size: 10pt; '>> >
Proj@lists.maptools.org</span><br><span style='font-size: 10pt; '>> >
http://lists.maptools.org/mailman/listinfo/proj</span><br><span style='font-size: 10pt; '>>
</span></div>
</body>
</html>