Problem with Oblique and Hotine oblique projections in proj?
Clever, Max
Maxc at SPICERGROUP.COM
Fri Dec 17 08:16:14 PST 2004
Hello,
I have been trying to project Michigan state plane coordinate shapefiles (EPSG:2253 or EPSG:102690 both are LCC) onto a Mapfile defined in Michigan Georef (EPSG:102123) which is a hotine oblique mercator. I think Proj is miscalculating the reprojection because the state plane data is showing up in the area of the false northing and easting of the georef projection. I have checked the parameters of the projection in the EPSG file and they are correct according to http://www.michigan.gov/documents/DNR_Map_Proj_and_MI_Georef_Info_20889_7.pdf. Also, I know that these parameters work when applied correctly, because I've written my own small VB program to convert coordinates from state plane to georef and back. To check that proj was actually calculating the projections, I checked it by assigning the State Plane projection to the mapfile then projecting a shapefile in UTM zone 17 onto it. That was flawless, but once I switch the Mapfile back to Georef proj miscalculates. To see where everything has been going I have set my extents to gargantuon numbers to encompass any possible placement of the projections and has helped me find where everything was going.
So now, if it really comes down to getting dirty I think I might need to open up PROJ4's guts and check out what its doing wrong. Hopefully though, I am the one that is screwing things up, so if anybody knows about how proj handles oblique mercators, I would be happy to hear any comments.
Also here is the portions of the mapfile that are pertinent to this message:
MAP
NAME zooba.map
STATUS ON
SIZE 400 400
EXTENT 129827.580 41110.205 823355.654 889381.112
#EXTENT 13067766 660421 13155653 722085
UNITS meters
SHAPEPATH "c:\webgisdata\zooba\v1c\"
SYMBOLSET "c:\ms4w\Apache\htdocs\zooba.org\zoobastart\gis\symbols.sym"
IMAGETYPE png24
FONTSET "c:\ms4w\Apache\htdocs\zooba.org\zoobastart\gis\fonts.list"
#GEOREF PROJECTION BASED USING GRS80 ELLIPSOID IN METERS
PROJECTION
"init=EPSG:102123"
# "proj=omerc"
# "lat_0=45.30916666666666"
# "lonc=-86"
# "alpha=337.255555555556"
# "k_0=0.9996"
# "x_0=2546731.496"
# "y_0=-4354009.816"
# "ellps=GRS80"
# "datum=NAD83"
# "units=m"
# " no_defs"
# "init=EPSG:2253"
END #PROJECTION
QUERYMAP
STATUS ON
STYLE SELECTED
COLOR 255 255 0
END
LAYER
NAME easements_4_27_04.shp
TYPE Polygon
TRANSPARENCY 50
# maxscale 1500000
STATUS on
PROJECTION
"init=epsg:102123"
END #PROJECTION
DATA easements_4_27_04
HEADER 'easements_4_27_04.shp_header.html'
TEMPLATE easements_4_27_04.shp_query.html
CLASS
Name 'Easements'
COLOR 100 100 231
OUTLINECOLOR 0 0 0
END # CLASS
END # LAYER
LAYER
NAME geoeasements_4_27_04.shp
TYPE Polygon
TRANSPARENCY 50
maxscale 1500000
STATUS on
PROJECTION
"init=EPSG:102690"
#"+proj=lcc"
#"+lat_1=43.66666666666666"
#"+lat_2=42.1"
#"+lat_0=41.5"
#"+lon_0=-84.36666666666666"
#"+x_0=3999999.999984"
#"+y_0=0 +ellps=GRS80"
#"+datum=NAD83"
#"+to_meter=0.3048"
#"+no_defs"
#"no_defs"
END #PROJECTION
DATA Geoeasements_4_27_04
HEADER 'easements_4_27_04.shp_header.html'
TEMPLATE easements_4_27_04.shp_query.html
CLASS
Name 'Geoeasements'
COLOR 255 138 231
OUTLINECOLOR 0 0 0
END # CLASS
END # LAYER
More information about the MapServer-users
mailing list