<br>Hello,<br><br>I am trying to crop and warp th<strong></strong>e enroute chart tif from FAA's <a href="http://aeronav.faa.gov/Enroute/04-05-2012/ENR_H01.zip">http://aeronav.faa.gov/Enroute/04-05-2012/ENR_H01.zip</a> package with the following commands:<br>
<br>1. gdal_translate -of vrt -srcwin 2203 204 21692 7591 ENR_H01.tif ENR_H01.vrt<br><br>2. gdalwarp -wm 1000 -t_srs EPSG:3395 -srcnodata 0 -dstalpha -r bilinear ENR_H01.vrt ENR_H01_tmp.tif<br><br>and when I place the input tif ENR_H01.tif and output tif ENR_H01_tmp.tif into ArcGIS Explorer to compare, I see that black lines on the original map have become transparent on the final map: I can see basemap through them. Also some of the black marks on the map have become dark-blue. So it looks to me that black color is gone or has become transparent and I am not sure how to deal with it since I am new to all of this. <br>
<br>Please advise! Thanks in advance for any tips/hints/directions on how to fix it.<br><br>I tried to run command #1 with -expand rgb but I get the following error message: "Error: band 1 has no color table"<br>
<br>Here's gdalinfo for ENR_H01.tif:<br><br>Driver: GTiff/GeoTIFF<br>Files: ../ENR_H01/ENR_H01.tif<br>Size is 24000, 7998<br>Coordinate System is:<br>PROJCS["US_IFR",<br> GEOGCS["NAD83",<br> DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.2572221010002,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6269"]],<br> PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4269"]],<br> PROJECTION["Lambert_Conformal_Conic_2SP"],<br> PARAMETER["standard_parallel_1",45],<br>
PARAMETER["standard_parallel_2",33],<br> PARAMETER["latitude_of_origin",39],<br> PARAMETER["central_meridian",-95],<br> PARAMETER["false_easting",0],<br> PARAMETER["false_northing",0],<br>
UNIT["meters",1]]<br>GeoTransform =<br> -2843013.39, 91.6788695, -13.0175628<br> 1594997.51, -13.0227943, -91.6740178<br>Metadata:<br> AREA_OR_POINT=Area<br> TIFFTAG_DATETIME=2012:03:15 10:46:33<br> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br>
TIFFTAG_SOFTWARE=GPL Ghostscript 8.71<br> TIFFTAG_XRESOLUTION=400<br> TIFFTAG_YRESOLUTION=400<br>Image Structure Metadata:<br> COMPRESSION=LZW<br> INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left (-2843013.390, 1594997.510) (133d52'54.95"W, 47d48'36.59"N)<br>
Lower Left (-2947127.857, 861788.716) (131d21'53.95"W, 41d23'27.77"N)<br>Upper Right ( -642720.522, 1282450.447) (103d53'35.10"W, 50d14'42.26"N)<br>Lower Right ( -746834.989, 549241.652) (104d17'24.63"W, 43d37'18.75"N)<br>
Center (-1794924.190, 1072119.581) (118d35'41.32"W, 46d33'40.12"N)<br>Band 1 Block=400x400 Type=Byte, ColorInterp=Red<br>Band 2 Block=400x400 Type=Byte, ColorInterp=Green<br>Band 3 Block=400x400 Type=Byte, ColorInterp=Blue<br>
<br><br><br> and gdalinfo for ENR_H01_tmp.tif<br><br><br>Driver: GTiff/GeoTIFF<br>Files: ENR_H01_tmp.tif<br>Size is 22689, 9778<br>Coordinate System is:<br>PROJCS["WGS 84 / World Mercator",<br> GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br> SPHEROID["WGS 84",6378137,298.257223563,<br> AUTHORITY["EPSG","7030"]],<br> AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4326"]],<br> PROJECTION["Mercator_1SP"],<br> PARAMETER["central_meridian",0],<br>
PARAMETER["scale_factor",1],<br> PARAMETER["false_easting",0],<br> PARAMETER["false_northing",0],<br> UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","3395"]]<br>Origin = (-14605257.985280666500330,6436384.405730826780200)<br>Pixel Size = (133.273420682087647,-133.273420682087647)<br>Metadata:<br> AREA_OR_POINT=Area<br>Image Structure Metadata:<br>
INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left (-14605257.985, 6436384.406) (131d12' 4.55"W, 50d 7'55.86"N)<br>Lower Left (-14605257.985, 5133236.898) (131d12' 4.55"W, 42d 0'20.66"N)<br>
Upper Right (-11581417.343, 6436384.406) (104d 2'15.51"W, 50d 7'55.86"N)<br>Lower Right (-11581417.343, 5133236.898) (104d 2'15.51"W, 42d 0'20.66"N)<br>Center (-13093337.664, 5784810.652) (117d37'10.03"W, 46d13'11.99"N)<br>
Band 1 Block=22689x1 Type=Byte, ColorInterp=Red<br> Mask Flags: PER_DATASET ALPHA<br>Band 2 Block=22689x1 Type=Byte, ColorInterp=Green<br> Mask Flags: PER_DATASET ALPHA<br>Band 3 Block=22689x1 Type=Byte, ColorInterp=Blue<br>
Mask Flags: PER_DATASET ALPHA<br>Band 4 Block=22689x1 Type=Byte, ColorInterp=Alpha<br><br><br>Thank you!!!<br>