[Mapserver-users] proj.4 stuff

Bob Basques bob.basques at ci.stpaul.mn.us
Thu Apr 22 15:14:14 PDT 2004


Frank Warmerdam wrote:

> Bob Basques wrote:
>
>> Ok, I can't get it to work:
>>
>> HELP!!
>>
>>
>
> This corresponds to the PROJ.4 definition:
>
>   +proj=lcc +lat_1=44.88333333333333 +lat_2=45.13333333333333 
> +lat_0=44.79111111111111 +lon_0=-93.38333333333334 +x_0=152400.3048 
> +y_0=30480.06096 +a=6378418.941 +b=6357033.31 
> +to_meter=0.3048006096012192 +no_defs
>
> I got this by putting your definition in a text file, trimming some of 
> the
> HTML junk out and then doing the following.  The testepsg program in 
> gdal/ogr
> (needs to be forced to build ... not a default utility) converts 
> between different
> coordinate system formats.  The ESRI:: prefix tells the program the. 
> prj file has
> ESRI WKT instead of OGR WKT.
>
>
Ok, that's cool, so that's in the GDAL stuff, interesting.  Learn 
something everyday.  Although I didn't see any reference to it anywhere 
in my searching.

So I took the above and made this Map File, is the correct way to 
implement (it didn't work, I get a blank page still):

MAP
  NAME 'test'
  SIZE 800 650
  STATUS ON
  EXTENT 332000 -18000 678000 328000
  UNITS FEET

  PROJECTION
    "init=epsg:26915"
  END

  LAYER
    NAME 'tlg01'
    DATA './data/tlg01/tlg_roads_l'
    STATUS DEFAULT
    TYPE LINE

    PROJECTION
      "proj=lcc"
      "lat_1=44.88333333333333"
      "lat_2=45.13333333333333"
      "lat_0=44.79111111111111"
      "lon_0=-93.38333333333334"
      "x_0=152400.3048"
      "y_0=30480.06096"
      "a=6378418.941"
      "b=6357033.31"
      "to_meter=0.3048006096012192"
      "no_defs"
    END
  END

END ## end Map

Thanbks for the help on this

bobb





More information about the MapServer-users mailing list