Image projection saga (continued)
Stephen Davies
sdc at SDC.COM.AU
Wed Jun 20 21:05:57 PDT 2007
Despite feedback from this list, I have still not resolved my problem
where data converted and transformed from a shape file to PostGIS does
not project in the same way as the geotiff on which the shape file was
based.
Previously, I was transforming the shape file coordinates from the
original SRID 28354 to 4283 for storage in the database. When this data
was overlaid on the raster image (SRID=28354), the polygons did not
line up.
The map projection definition is:
PROJECTION
"proj=eqc"
"ellps=GRS80"
"towgs84=0,0,0"
"no_defs"
"lon_0=139.591e"
END
I have now tried creating a polygon (SRID 28354) using the boundary
coordinates of the raster image as determined from gdalinfo and
overlaying that onto the image. (This removes the transform to 4283 and
the subsequent reprojection to eqc.)
Still I see the polygon shift north by about 20m and by a tiny amount
east.
This suggests to me that the method used to project geometries is
somehow different from that used to project rasters as both sources now
have the same SRID and could reasonably be expected to produce the same
results.
To further confuse me, if I change the map projection to either 28354 or
4283, everything lines up (but my scale bar and reference maps break).
I have searched all the documentation I can find but cannot find
anything that describes the relationship between maps, sale bars and
reference maps apart from a statement that reference map extent must be
expressed in map units.
TIA,
Stephen Davies
--
=======================================================================
This email is for the person(s) identified above, and is confidential to
the sender and the person(s). No one else is authorised to use or
disseminate this email or its contents.
Stephen Davies Consulting Voice: 08-8177 1595
Adelaide, South Australia. Fax: 08-8177 0133
Computing & Network solutions. Mobile:0403 0405 83
More information about the MapServer-users
mailing list