<html>
<p>I have found out that the DIgital orthos are Virginia State Plane Coordinate system. Through the metadata on the dvd, I found that although the discs read VASPCS North, The metadata actually refers to the VASPCS south. I metadata also refers to the units in US survey foot. whereas the epsg file refers to metre. I have tried using the information from the meta to create a separate projection, but it is not working correctly. Can you see anything wrong in my translation.</p>
<p>Spatial_Reference_Information:<br/>
Horizontal_Coordinate_System_Definition:<br/>
Planar:<br/>
Grid_Coordinate_System:<br/>
Grid_Coordinate_System_Name: State Plane Coordinate System 1983<br/>
State_Plane_Coordinate_System:<br/>
SPCS_Zone_Identifier: 4502<br/>
Lambert_Conformal_Conic:<br/>
Standard_Parallel: 36.766667<br/>
Standard_Parallel: 37.966667<br/>
Longitude_of_Central_Meridian: -78.500000<br/>
Latitude_of_Projection_Origin: 36.333333<br/>
False_Easting: 3500000.000000<br/>
False_Northing: 1000000.000000<br/>
Planar_Coordinate_Information:<br/>
Planar_Coordinate_Encoding_Method: row and column<br/>
Coordinate_Representation:<br/>
Abscissa_Resolution: 1.000000<br/>
Ordinate_Resolution: 1.000000<br/>
Planar_Distance_Units: US Survey Foot<br/>
Geodetic_Model:<br/>
Horizontal_Datum_Name: North American Datum of 1983 (HARN)<br/>
Ellipsoid_Name: Geodetic Reference System 80<br/>
Semi-major_Axis: 6378137.000000<br/>
Denominator_of_Flattening_Ratio: 298.257222<br/>
</p>
<p>I set up in the layer record:(here's the testing map file)</p>
<p>MAP<br/>
NAME "City of Harrisonburg, Virginia"<br/>
SIZE 400 400<br/>
STATUS ON<br/>
# SHAPEFILE EXTEXTS<br/>
EXTENT 1864266.6510595274 259839.01912492423 1907053.0983945231 299828.4636574<br/>
4545<br/>
#RASTER EXTENTS<br/>
# EXTENT 11360000.0 6820000.0 11390000.0 6870000.0<br/>
UNITS FEET<br/>
SHAPEPATH "/opt/mapserver/www/htdocs/harrisonburg/data"<br/>
PROJECTION<br/>
"init=epsg:102746"<br/>
END</p>
<p> WEB<br/>
IMAGEPATH "/opt/mapserver/www/htdocs/tmp/"<br/>
IMAGEURL "/tmp/"<br/>
LOG "/opt/mapserver/www/htdocs/tmp/ms.log"<br/>
END #end web</p>
<p> # 3876.sid<br/>
LAYER<br/>
NAME 'layer13'<br/>
DATA '/opt/mapserver/www/htdocs/harrisonburg/data/digitalorthos/3885.sid'<br/>
TYPE raster<br/>
#MAXSCALE 25000<br/>
PROJECTION<br/>
"proj=latlong"<br/>
"ellps=GRS80"<br/>
"datum=NAD83"<br/>
"lat_0=36.333333334"<br/>
"lat_1=36.766666667"<br/>
"lat_2=37.966666667"<br/>
"lon_0=-78.5"<br/>
"units=ft"<br/>
END<br/>
STATUS ON<br/>
END #end layer</p>
<p> # RAILROADS<br/>
LAYER<br/>
NAME 'Railroads'<br/>
DATA '/opt/mapserver/www/htdocs/harrisonburg/data/railroad'<br/>
CLASS<br/>
STYLE<br/>
COLOR 255 102 200<br/>
SYMBOL 0<br/>
END #end style<br/>
END #end class<br/>
TRANSPARENCY 80<br/>
TYPE LINE<br/>
STATUS ON<br/>
END #end layer</p>
<p>END #end map</p>
<p><br/>
Paul A. Malabad<br/>
Director of Information Technology<br/>
City of Harrisonburg<br/>
(540) 432-7706<br/>
<br/>
------ Original Message ------<br/>
From: Frank Warmerdam <warmerdam@pobox.com><br/>
Date: Saturday, July 8th, 2006 10:24 AM EDT<br/>
To: Paul A Malabad <malabapa@CI.HARRISONBURG.VA.US><br/>
Subject: Re: [UMN_MAPSERVER-USERS] raster vs shp files extents don't match<br/>
<br/>
Paul A Malabad wrote:<br/>
> Corner Coordinates:<br/>
> Upper Left (11370000.000, 6860000.000)<br/>
> Lower Left (11370000.000, 6850000.000)<br/>
> Upper Right (11380000.000, 6860000.000)<br/>
> Lower Right (11380000.000, 6850000.000)<br/>
...<br/>
> Extent: (1880340.635992, 264114.932237) - (1906873.587035, 299715.644746)<br/>
...<br/>
> I have tried using the different HARN and NAD27 epsg projections in the <br/>
> layer, but I think the problem is more basic (at least I hope it is.) <br/>
> since the extents are so skewed. Any help or a pointer to the right <br/>
> manual would be greatly appreciated.<br/>
<br/>
Paul,<br/>
<br/>
The problem here is that without any coordinate system information for<br/>
your MrSIDs you are left playing the "guess the coordinate system" game.<br/>
And I would suggest that the rest of us out here in internet land have even<br/>
less to go on than you.<br/>
<br/>
If you can't get information from the original MrSID producers, then you<br/>
will basically have to mess around with a variety of guesses. The coordinates<br/>
for your raster are enormous. If they are in meters then 11 million meters<br/>
is roughly 1/4 of the way around the world from the projection origin.<br/>
<br/>
The ratio of x to y coordinates is quite different between the vector and<br/>
raster data, so if they are near each other in the real world, it is unlikely<br/>
this is just a units issue.<br/>
<br/>
I hate playing guess the coordinate system.<br/>
<br/>
Best regards,<br/>
-- <br/>
---------------------------------------+--------------------------------------<br/>
I set the clouds in motion - turn up | Frank Warmerdam, warmerdam@pobox.com<br/>
light and sound - activate the windows | http://pobox.com/~warmerdam<br/>
and watch the world go round - Rush | President OSGF, http://osgeo.org<br/>
<br/>
</p>
</html>