Merging different projection

[Sten-Helge Henriksen]

On Fri, 18 Feb 2005 08:50:33 -0500, Ed McNierney <ed at TOPOZONE.COM> wrote:

>Sten-Helge -
>
>Could you post the map file you are using that does not work properly?
>Are you describing a PROJECTION block for each layer, selecting an
>output PROJECTION for the map image, and setting the output map EXTENT
>to match the output projection?
>
>        - Ed
>
>Ed McNierney
>TopoZone.com
>
>-----Original Message-----
>From: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU] On
>Behalf Of Sten-Helge Henriksen
>Sent: Friday, February 18, 2005 6:53 AM
>To: MAPSERVER-USERS at LISTS.UMN.EDU
>Subject: [UMN_MAPSERVER-USERS] Merging different projection
>
>Hi all!
>I have two different types of shape files whice I would like to use in
>the
>same .map.
>The projection from the first corresponding .prj file looks like this;
>
>PROJCS["WGS_1984_UTM_Zone_33N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984"
>,SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0]
>,UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],
>PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],
>PARAMETER["Central_Meridian",15.0],PARAMETER["Scale_Factor",0.9996],
>PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]
>
>while from the second, it looks like this;
>GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,2
>98.257223563]
>],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
>
>I have tried several ways, with differnt types of projection to get
>those
>two filetypes into the same map, without success.
>
>Separately they turn up nicely with respectively extent;
>EXTENT 404929.999959 7516685.438843 424106.999991 7544257.142516
>and
>EXTENT -4.792797 37.909400 51.301297 81.166931
>
>Does anyone have a clue how to solve this?
>
>Thanks in advance
>
>stenhh

Hello there!

Here is a dump of my mapfile with those two layers in qustion. Layer #1
takes data from a dataset based on proj=utm, zone=33, ellps=WGS, datum=WGS84
and units=meters. Together with extent #1, it works perfectly.
Layer #2 with extent #2 does also work perfectly.

However, my objective is to get those two dataset to work together in the
same .map, even though they have different projections, and I can't figure
out which projection to use on layer #2. Even though I have the
corresponding .prj file for the second dataset, it dosn't seem to work. Just
a blank screen with no errors.

MAP
  NAME cmis
  STATUS ON
  EXTENT 404929.999959 7516685.438843 424106.999991 7544257.142516 #extent 1
# Extent -4.792797 37.909400 51.301297 81.166931 #extent 2
  SIZE 800 600
  SYMBOLSET "./etc/symbols.sym"
  UNITS meters
  FONTSET ./etc/fonts.txt
  DEBUG ON
 TEMPLATEPATTERN "cmis"
  WEB
    IMAGEPATH "set in index.html"
    IMAGEURL "set in index.html"
    TEMPLATE test.html84
    LOG "/tmp/mapserv.log"
  END

##### Layer #1

 LAYER
    NAME "Kontur"
    STATUS default
    TYPE line
    CONNECTIONTYPE postgis
    CONNECTION "user=xxxx dbname=xxxx host=xxxxx"
    DATA "the_geom from areal_line"
    FILTER "objtype='Coastal'"
    CLASS
      NAME "Kontur"
      #EXPRESSION "0"
      STYLE
        SYMBOL 0
        COLOR 225 0 0
        OUTLINECOLOR 186 234 245
        SIZE 1
        MINSIZE 1
        MAXSIZE 100
      END
    END

  END

#### Layer #2
  LAYER
    NAME "country_boundaries"
    STATUS default
    TYPE line
    UNITS meters
    CONNECTIONTYPE postgis
    CONNECTION "user=xxxx dbname=xxxx host=xxxx"
    DATA "the_geom from country_boundaries"
    PROJECTION
#      "proj=utm"
#      "zone=33"
#      "ellps=WGS84"
#      "datum=WGS84"
#      "units=m"
#      "no_defs"
       "init=epsg:4326" # I have tried many of these wihtout success
#       "proj=latlong"
    END
    CLASS
      NAME "country_boundaries"
      #EXPRESSION "0"
      STYLE
        SYMBOL 0
        COLOR 225 0 0
        OUTLINECOLOR 186 234 245
        SIZE 1
        MINSIZE 1
        MAXSIZE 100
      END
    END
  END

END #MAP