<div>mapserver rocks on-the-fly!</div>
<div> </div>
<div>now what's the best way to avoid on-the-fly proj because that decreases performance. do i have to manually rectify geotiff to vector in arcmap (think human error)?</div>
<p>thanks a bunch man</p>
<p>jzs</p>
<p><a href="http://mapserver.gis.umn.edu/docs/howto/raster_data/#raster-display-performance-tips">http://mapserver.gis.umn.edu/docs/howto/raster_data/#raster-display-performance-tips</a><br><br></p>
<div><span class="gmail_quote">On 11/16/06, <b class="gmail_sendername">Frank Warmerdam</b> <<a href="mailto:warmerdam@pobox.com">warmerdam@pobox.com</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">John Smith wrote:<br>> guys,<br>> trying to put a geotiff in- works but not right. here's what's happening:
<br>><br>> i download the geotiff from national atlas. when i give no proj info<br>> under LAYER it shows up but not where it should (nwis). so i copy MAP<br>> proj to LAYER proj hoping for some on-the-fly magic- again it shows up
<br>> but nwis.<br>><br>> well i go over to proj_list<br>> <<a href="http://www.remotesensing.org/geotiff/proj_list/">http://www.remotesensing.org/geotiff/proj_list/</a>> and look up the epsg<br>> (EPSG?) code (9820) from metadata but it throws the
<br>> "msProcessProjection(): Projection library error. no options found in<br>> 'init' file" error. then i try all epsg codes given by gdalinfo ;) but<br>> same error. although when i try the MAP epsg code (4326) it shows up but
<br>> nwis.<br>><br>> so what's wrong here? think proj.exe geod.exe cs2cs.exe libgeotiff<br>> listgeo'd help? now i don't think wld is needed because it's included<br>> with the geotiff but i've tried with wld and no luck. i just hope i
<br>> don't have to rectify in arcmap. btw IMAGETYPE=png24 in case that matters.<br><br>John,<br><br>I don't believe there is an EPSG code corresponding to this coordinate<br>system:<br><br>PROJCS["TIFF 6.0, GeoTIFF Revision
1.0",<br> GEOGCS["WGS 84",<br> DATUM["WGS_1984",<br> SPHEROID["WGS 84",6378137,298.2572235630016,<br> AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4326"]],<br> PROJECTION["Lambert_Azimuthal_Equal_Area"],
<br> PARAMETER["latitude_of_center",45],<br> PARAMETER["longitude_of_center",-100],<br> PARAMETER["false_easting",0],<br> PARAMETER["false_northing",0],<br> UNIT["metre",1,
<br> AUTHORITY["EPSG","9001"]]]<br><br>But the corresponding PROJECTION block to put in the LAYER definition<br>is:<br><br>PROJECTION<br> "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +datum=WGS84 +units=m"
<br>END<br><br>At one time you could also put:<br>PROJECTION<br> "AUTO"<br>END<br><br>and MapServer would read it from the geotiff file, but I'm not sure if this<br>still works.<br><br>Best regards,<br>--<br>---------------------------------------+--------------------------------------
<br>I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam">http://pobox.com/~warmerdam
</a><br>and watch the world go round - Rush | President OSGeo, <a href="http://osgeo.org">http://osgeo.org</a></blockquote></div>