I am still battleing this conversion. I was just informed, via IRC,
that gdal won't overwrite existing files, so I may have been making a
difference, but not knowing it. In any case, I tossed together a
simple script, that thus far, does nothing...
#!/bin/sh
rm -rf nn.tif n-merc.tif
echo '*** Translate ***'
gdal_translate n.tif nn.tif -of gtiff -CO TFW=YES -mo
TIFFTAG_XRESOLUTION= -mo TIFFTAG_YRESOLUTION= -a_ullr -115d0 56d0
-105d0 54d0 -a_srs "+proj=lcc +lat_0=54 +lon_0=-115 +lat_1=49.33
+lat_2=54.66 +x_0=0 +y_0=0 +datum=NAD27 +ellps=clrk80"
echo '*** Warp ***'
gdalwarp -t_srs "+proj=merc +datum=WGS84" nn.tif n-merc.tif
open n-merc.tif
I have also posted the source file I am working with, albeit scaled
down and converted so its not 140mb.
http://game-sat.com/~brian/north.png
What concerns me is what gdalinfo gives me for the intermediate step,
as the georeferenced coordinates are dead wrong. If I remove the a_srs
to assign it as being an actual lambert conformal projection, it
doesn't even have the georeferenced options, just the first ones.
$ gdalinfo nn.tif
Driver: GTiff/GeoTIFF
Files: nn.tif
Size is 11696, 4208
Coordinate System is:
PROJCS["unnamed",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",49.33],
PARAMETER["standard_parallel_2",54.66],
PARAMETER["latitude_of_origin",54],
PARAMETER["central_meridian",-115],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-115.000000000000000,56.000000000000000)
Pixel Size = (0.000854993160055,-0.000475285171103)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_SOFTWARE=Adobe Photoshop CS3 Macintosh
TIFFTAG_DATETIME=2009:02:21 16:15:49
TIFFTAG_XRESOLUTION=0
TIFFTAG_YRESOLUTION=0
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left (-115.0000000, 56.0000000) (115d 0'6.32"W, 54d 0'1.81"N)
Lower Left (-115.0000000, 54.0000000) (115d 0'6.32"W, 54d 0'1.75"N)
Upper Right (-105.0000000, 56.0000000) (115d 0'5.77"W, 54d 0'1.81"N)
Lower Right (-105.0000000, 54.0000000) (115d 0'5.77"W, 54d 0'1.75"N)
Center (-110.0000000, 55.0000000) (115d 0'6.04"W, 54d 0'1.78"N)
Band 1 Block=11696x1 Type=Byte, ColorInterp=Red
Band 2 Block=11696x1 Type=Byte, ColorInterp=Green
Band 3 Block=11696x1 Type=Byte, ColorInterp=Blue
On Sat, Feb 21, 2009 at 9:55 PM, Brian Murray <idle...@gmail.com> wrote:
I'm trying to get all of it. :)
It appears like a_allr is ok, but then when I apply the definition, it
only grabs a few meters of the real map...
I think I must have to set the scale to something. I just don't know
what yet... I am going to sleep on it. I found some sites that have
some calculations for that now.
On Sat, Feb 21, 2009 at 9:47 PM, Wendy Fay Stevenson <steve...@erau.edu>
wrote:
Um, it looks like it's doing what you're asking it to do.
What corner points are you wanting?
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
