Guess this Projection

[Nick Dabner]

Used libgeo to make a .tfw file of the raster. Very strange output:

             2.5390625000
             0.0000000000
             0.0000000000
            -2.5390625000
       2478259.8695312501
       4427159.8914062502

Could the results be random. I cannot recognise the format, and proj4
returns a bus error. This raster is a citymap of central Melbourne,
Australia at around 137E -40S WGS84.

ND

On 06/10/2004, at 8:29 PM, Nick Dabner wrote:

> 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
>>
>>
>