Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I have a basic question about my MapServer .map file, which has two layers, each with a different coordinate system, longlat vs. mercator.

See my .map file below. If I leave it as shown below, with the "version 2" UNITS/EXTENT active, it correctly draws my "MrSid" in the right place on the earth, but doesn't draw my "shape layer" layer. On the other hand, if I comment out the "version 2" UNITS/EXTENT and activate the "version 1" UNITS/EXTENT, it then draws "shape layer" but not "MrSid".

So it's acting like the the UNITS/EXTENT has to match those of the layer for the layer to be drawn. Is that right? If so, how do I mix layers with different UNITS?

(By the way, I also have a secondary question about this .map file. If I don't include any PROJECTION entry at the outer/MAP level, it will draw the "shape layer" layer but not the "MrSid" layer. But if I add a PROJECTION, regardless of whether the projection is "proj=longlat" "ellps=WGS84" "datum=WGS84" or "proj=merc" "lon_0=0" "k=1" "x_0=0" "y_0=0" "a=6378137" "b=6378137" "towgs84=0,0,0,0,0,0,0" "units=m", it will draw either layer -- but not both at the same time, as per the main issue above. So my secondary question is why it seems to be happy using a default projection with my longlat shape file, but not with my mercator MrSid file.)

MAP
  IMAGETYPE   PNG

  # version 1
  # UNITS dd
  # EXTENT      -180 -90 180 90

  # version 2
  UNITS meters
  EXTENT -20037508.34 -20037508.34 20037508.34 20037508.34

  SIZE        550 450
  INTERLACE OFF
  PROJECTION
    # "proj=longlat" "ellps=WGS84" "datum=WGS84"
    "proj=merc" "lon_0=0" "k=1" "x_0=0" "y_0=0" "a=6378137" "b=6378137" "towgs84=0,0,0,0,0,0,0" "units=m"
  END
  SHAPEPATH   "../data"

  OUTPUTFORMAT
    NAME png
    DRIVER "GD/PNG"
    MIMETYPE "image/png"
    IMAGEMODE RGB
    FORMATOPTION  INTERLACE=OFF
    TRANSPARENT OFF
  END

  # Start of LAYER DEFINITIONS-------------------------------
  LAYER
    NAME "shape layer"
    DATA dourokukan
    STATUS DEFAULT
    TYPE LINE

    CLASS
      NAME "“shape class"
      STYLE
        COLOR 227 227 127
      END
    END
  END 
  LAYER
    NAME "MrSid"
    DATA "c:\mrsid\adams\ortho_1-1_1n_s_nd001_2010_1.sid"
    STATUS DEFAULT
    TRANSPARENCY 100
    TYPE RASTER
    PROJECTION
      "proj=utm" "ellps=GRS80" "datum=NAD83" "zone=13" "units=m"
    END
    LABELITEM ""
    TEMPLATE VOID
  END
  # End of LAYER DEFINITIONS -------------------------------
END 
share|improve this question

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.