[Mapserver-users] Once again: no options found in 'init' file

[Bjoern Platzen]

Hello Andreas,

now I had the time to try if your advice works for me...

... it does!

Thanks a lot,

bye,

Bjoern.


Am Mittwoch, 6. August 2003 18:03 schrieb Maul, Andreas-Alexander:
> Hi Bjoern,
> I do in my application a projection from lat/lon-WGS84 coordinates to
> Gauss-Krueger (9 deg meridian) coordinates because the coordinates in
> my shapefiles are related to WGS84 whereas the raster topography
> bought from Bundesamt für Kartographie und Geodaesie (BKG Germany)
> refers to Gauss-Krueger coordinates stripe 3 (9 deg east). I specify
> the DHDN parameters without use of the epsg file. You can try to do
> your projection with similar parameters. My projection section of the
> map file looks like:
>
> PROJECTION
> "proj=tmerc"
> "ellps=bessel"
> "towgs84=dx,dy,dz,rx,ry,rz,m" # see notes below
> # datum shift DHDN Zone 3,
> # note: if datum shift is applied here, then it must be specified in
> the # layer section too, eg. "datum=WGS84"
> #
> "lat_0=0n"     # latitude of origin
> "lon_0=9e"     # meridian for stripe number 3
> "x_0=3500000." # Stripe number 3 for 9 deg E, Offset: 500000
> "y_0=0."       # offset y
> "k=1."         # scale factor
> "units=m"
> END
>
> In the layer section I have:
> LAYER
>   NAME "Border"
>   DATA "..."
>   TYPE LINE
>   STATUS ON
>   PROJECTION
>    "proj=latlong"
>    "datum=WGS84"
>   END
>   ...
> END
>
> Note 1: If you need to perform a datum shift you have to take care
> about the signs of the translation and rotation parameters of the
> ellipsoid.
>
> There were some discussions about the signs of translation and
> rotation parameters and the scale factor used by the Europeans and
> Americans (i.e. as defined in the proj.4 package) in the user list of
> the proj archive. See e.g.
> http://www.remotesensing.org/lists/proj_archive/msg00298.html
>
> Note 2: If you need the Potsdam Datum parameters, you can find
> approximate numbers in http://ottmarlabonde.de/L1/datum.html The
> official datum parameters for the Potsdam Datum from BKG are
> confidential.
>
> You can try the following numbers:
> dx = 586 m, dy = 87 m, dz = 409 m, rx = 0.52 ", ry = 0.15 ", rz =
> -2.82 " m = 1.000009 (9ppm)
>
> therefore you can insert the following statement in the projection
> section: "towgs84=586, 87, 409, 0.52, 0.15, -2.82, 1.000009"
> and look whether the accuracy is sufficient for you.
>
> Regards,
> Andreas
>
>
> -----Ursprüngliche Nachricht-----
> Von: Bjoern Platzen [mailto:bplatzen at sosnetz.de]
> Gesendet: Dienstag, 5. August 2003 18:30
> An: mapserver-users at lists.gis.umn.edu
> Betreff: [Mapserver-users] Once again: no options found in 'init'
> file
>
> Hello List,
>
> by searching through the archives, I saw that this seems to be a
> common problem. But I did not find any solution that helped me...
>
> I'm trying to project a lat/lon-WGS84-Layer on Gauss-Kruger-Layers.
> In my mapfile it looks like this:
>
> # Start of map file
> NAME mgla2
> STATUS ON
> SIZE 6400 6400
> SYMBOLSET symbols/symbols.sym
> EXTENT 2523000 5680000 2539000 5663000
> UNITS METERS
> SHAPEPATH "gis_data"
> FONTSET /usr/local/httpd/htdocs/parkinfoneuss/font/font.list
> IMAGECOLOR 255 250 238
> IMAGETYPE png
>
> #
> # Start of web interface definition
> #
> WEB
>   TEMPLATE rahmen.html
>   MINSCALE 3500
>   MAXSCALE 200000
>   IMAGEPATH "/usr/local/httpd/htdocs/parkinfo/tmp/"
>   IMAGEURL "tmp/"
>   LOG "/usr/local/httpd/htdocs/parkinfo/mg.log"
>   error "http://moers.office.sos/parkinfo/error.htm"
> END
>
> PROJECTION
>   "init=epsg:31466"
> END
>
> LAYER
>   NAME bebauung_de012w
>   TYPE Polygon
>   STATUS on
>   DATA bebauung_de012w
>       CLASS
>         COLOR 255 165 165
>         OUTLINECOLOR -1 -1 -1
>      END  # CLASS
> END  # LAYER
>
> LAYER
>   NAME industrie_de012w
>   TYPE Polygon
>   STATUS on
>   DATA industrie_de012w
>       CLASS
>         COLOR 165 165 165
>         OUTLINECOLOR -1 -1 -1
>      END  # CLASS
> END  # LAYER
>
> #...some more layers...
>
> LAYER
>   NAME 'oe'#this is a lat/long-Layer
>   TYPE Polygon
>   STATUS on
>   DATA oeffentl_einr_de012w
>   CLASS
>     COLOR 205 44 65
>   END#CLASS
>   PROJECTION
>     "init=epsg:4326"
>   END#projection
> END#Layer
>
> #...some more layers...
> END
> #End of mapfile
>
> Now I get the Error
> msProcessProjection(): Projection library error. no options found in
> 'init' file
>
> I tried as well
> platzen at codefabrik:~/parkinfo/mg/karte> proj -v +init=epsg:31466
> #Transverse Mercator
> #       Cyl, Sph&Ell
> # +init=epsg:31466 +proj=tmerc +lat_0=0 +lon_0=6 +k=1.000000
> +x_0=2500000
> # +y_0=0 +ellps=bessel +units=m +no_defs
> -> epsg:31466 seems to work
>
> and
>
> platzen at codefabrik:~/parkinfo/mg/karte> proj -v +init=epsg:4326
> Rel. 4.4.7, 31 March 2003
> <proj>:
> +proj=latlong unsuitable for use with proj program.
> program abnormally terminated
> -> epsg:4326 does _not_ seem to work
>
> Then I checked my epsg-file unter /usr/local/share/proj/ and i found
> # WGS 84
> <4326> +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs  no_defs <>
>
> So I think, the error comes because of epsg:4326 but I have no idea,
> why it comes and how to workaround...
>
> Any help would be _very_ appreciated...
>
> thanks,
>
> Bjoern.

-- 
small office solutions
info at sosnetz.de  -  http://www.sosnetz.de