Projections with ARCINFO Grids

[Aaron Racicot]

Frank,

Thanks for the quick reply...

It turns out that 'PROJECTION "AUTO" END' does work and my maps are
displaying great.  I just hope the functionality stays working in future
versions!

In the case of my previous example, the new mapfile layer def is:

  LAYER
    DATA "siuslaw/data_cd/region/shd90/"
      METADATA
      END
    NAME "siuslaw_data_cd_region_shd90"
    SIZEUNITS PIXELS
    STATUS ON
    TOLERANCE 0
    TOLERANCEUNITS PIXELS
    TYPE RASTER
    UNITS METERS
    PROJECTION 
      "AUTO" 
    END
  END

And the resulting mapfile after saving it back out via php_mapscript
with the projection filled in automatically by Mapserver:

  LAYER
    DATA "siuslaw/data_cd/region/shd90/"
      METADATA
      END
    NAME "siuslaw_data_cd_region_shd90"
    PROJECTION
      "proj=utm"
      "zone=10"
      "ellps=clrk66"
      "datum=NAD27"
      "units=m"
      "no_defs"
    END
    SIZEUNITS PIXELS
    STATUS ON
    TOLERANCE 0
    TOLERANCEUNITS PIXELS
    TYPE RASTER
    UNITS METERS
  END

Thanks again...

Aaron

-----Original Message-----
From: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU] On
Behalf Of Frank Warmerdam
Sent: Wednesday, May 04, 2005 12:26 PM
To: MAPSERVER-USERS at LISTS.UMN.EDU
Subject: Re: [UMN_MAPSERVER-USERS] Projections with ARCINFO Grids

On 5/4/05, Aaron Racicot <aaronr at ecotrust.org> wrote:
> Hi there,
> 
> I am having issues with displaying some Arc/Info Binary Grid's.  I was
hoping that I would not have to explicitly define the projection
information for layers that I am dynamically adding to my mapfile since
the format they are in has a prj.adf file defining the projection.  I am
displaying my main map in lat-lon, but have a couple of rasters that are
in different projections (such as UTM zone 10N).  If I explicitly define
the projection in the map file things display fine for the layer:

Aaron, 

At one time you could put 'PROJECTION "AUTO" END' in your mapfile
to auto-load the projection of a layer from a GDAL file, but I don't
think
this works anymore.  It certainly isn't the default.  You can file a bug
on this topic if you like.

In general MapServer seems to like you to have to regenerate your
projection definition even when it could figure it out for you.  This is
an
annoying source of errors and confusion.

Best regards,
-- 
---------------------------------------+--------------------------------
------
I set the clouds in motion - turn up   | Frank Warmerdam,
warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent