Guess this Projection

[Nick Dabner]

Thanks for your help Peter. I misquoted the zone - it is actually 53,
and I adjusted your parameters to suit. No luck unfortunately. The tiff
is sourced straight of the Melways (eg Ausways / Sydways) CD.

As I said, I can provide the sample file to a worthy challenger ;-)

Nick

On 06/10/2004, at 3:15 PM, Peter Bayley wrote:

> Nick,
>
> The AMG66 Zone 55 Projection needs to be placed in the LAYER object in
> the
> map file:
>
> Taking the proj4text vales from the spatial_ref_sys table in postgis,
> you
> can construct the
> projection of the layer as:
>
> LAYER
>     PROJECTION
>         "proj=utm"
>         "zone=55"
>         "ellps=aust_SA"
>         "south"
>         "units=m"
>     END
>      CONNECTIONTYPE OGR
>      CONNECTION "ADELBRIS.TIF"
>    TYPE RASTER
>    STATUS DEFAULT
>    UNITS METERS
> END
>
> The Projection in the LAYER should be the projection of its dataset.
> By the
> way, 20253 is Zone 53 not 55
> The projection you want everything delivered in should be in the MAP
> section
>
> Regards
>
> Peter Bayley
> Sydney
> Australia
>
> ----- Original Message -----
> From: "Nick Dabner" <dabner at OPTUSNET.COM.AU>
> To: <MAPSERVER-USERS at LISTS.UMN.EDU>
> Sent: Wednesday, October 06, 2004 12:55 PM
> Subject: [UMN_MAPSERVER-USERS] Guess this Projection
>
>
>> Hello all,
>>
>> Have one nasty little raster here that is apparently in AMG66 Zone 55
>> (epsg:20253), according to Mapinfo. I cannot get it to display in
>> mapserver unless I warp it, but doing so results in an image that is
>> incorrectly registered and skewed (also out by 1500km or so). I am
>> trying to output in WGS84 / GDA94 in line with my vectors. Gdalinfo
>> output is shown below. I can provide a copy by email if anyone would
>> care to take a closer look. This has been causing me problems for
>> months now. Mapfile is below. Any help greatly appreciated, Nick
>>
>> Driver: GTiff/GeoTIFF
>> Size is 4809, 4075
>> Coordinate System is:
>> PROJCS["unnamed",
>>      GEOGCS["unnamed",
>>          DATUM["unknown",
>>              SPHEROID["unretrievable - using
>> WGS84",6378137,298.257223563]],
>>          PRIMEM["Greenwich",0],
>>          UNIT[,0.0174532925199433],
>>          AUTHORITY["EPSG","0"]],
>>      UNIT["unknown",1],
>>      AUTHORITY["EPSG","0"]]
>> Origin = (1857716.251568,5544261.870844)
>> Pixel Size = (317.52333125,-317.52333125)
>> Metadata:
>>    TIFFTAG_IMAGEDESCRIPTION=
>>    TIFFTAG_DATETIME=2004:01:05 22:39:05
>>    TIFFTAG_XRESOLUTION=300
>>    TIFFTAG_YRESOLUTION=300
>>    TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>> Corner Coordinates:
>> Upper Left  ( 1857716.252, 5544261.871)
>> (106439300d44'382252089395.85"E,317662805d42'1142461300788.73"N)
>> Lower Left  ( 1857716.252, 4250354.296)
>> (106439300d44'382252089395.85"E,243527362d35'876173328429.93"N)
>> Upper Right ( 3384685.952, 5544261.871) (193928220d
>> 0'700079669251.66"E,317662805d42'1142461300788.73"N)
>> Lower Right ( 3384685.952, 4250354.296) (193928220d
>> 0'700079669251.66"E,243527362d35'876173328429.93"N)
>> Center      ( 2621201.102, 4897308.083)
>> (150183760d22'541165879323.76"E,280595084d 9'1009317314579.33"N)
>> Band 1 Block=4809x1 Type=Byte, ColorInterp=Palette
>>
>>
>> The mapfile looks roughly like this (originally based on the itasca
>> demo):
>>
>> #
>> # Start of map file
>> #
>>
>> NAME DEMO
>> STATUS ON
>> SIZE 600 600
>> EXTENT 144.5 -37.5 145.5 -38
>> UNITS KILOMETERS
>> SHAPEPATH "data"
>> IMAGECOLOR 255 255 255
>>
>>      PROJECTION
>>        "init=epsg:4326"
>>      END
>>
>> WEB
>>    HEADER demo_header.html
>>    TEMPLATE demo.html
>>    FOOTER demo_footer.html
>>    IMAGEPATH "set in demo_init.html"
>>    IMAGEURL "set in demo_init.html"
>> END
>>
>> QUERYMAP
>>    SIZE 200 200
>>    STATUS ON
>>    STYLE HILITE
>>    COLOR 255 0 0
>> END
>>
>> SYMBOL
>>    NAME 'circle'
>>    TYPE ELLIPSE
>>    POINTS 1 1 END
>>    FILLED TRUE
>> END
>>
>> LAYER
>>      PROJECTION
>>        "init=epsg:4202" # Also Zone 53 "20253"
>>      END
>>      CONNECTIONTYPE OGR
>>      CONNECTION "ADELBRIS.TIF"
>>    TYPE RASTER
>>    STATUS DEFAULT
>>    UNITS METERS
>> END
>>
>> END # Map File
>
>