How to convert the convert the FAA NARC.tif file?

0 votes
asked Jun 20, 2016 by ba3user (4,000 points)

I’m trying to convert the FAA NARC.tif (located herehttps://www.faa.gov/air_traffic/flight_info/aeronav/digital_products/ifr/ and under the planning tab) image to be able to use in the BA3 engine but I’m not having any luck converting it to RGBA to run in AltusRaster.  These are the results from the terminal:

Nate:qbert e192619$ gdalinfo /Users/Nate/Downloads/narc_tif/narc.tif

Driver: GTiff/GeoTIFF

Files: /Users/Nate/Downloads/narc_tif/narc.tif

Size is 13688, 9719

Coordinate System is:

PROJCS["NARC",

    GEOGCS["NAD83",

        DATUM["North_American_Datum_1983",

            SPHEROID["GRS 1980",6378137,298.2572219999988,

                AUTHORITY["EPSG","7019"]],

            AUTHORITY["EPSG","6269"]],

        PRIMEM["Greenwich",0],

        UNIT["degree",0.0174532925199433],

        AUTHORITY["EPSG","4269"]],

    PROJECTION["Lambert_Conformal_Conic_2SP"],

    PARAMETER["standard_parallel_1",51],

    PARAMETER["standard_parallel_2",14],

    PARAMETER["latitude_of_origin",32.5],

    PARAMETER["central_meridian",-24],

    PARAMETER["false_easting",0],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

        AUTHORITY["EPSG","9001"]]]

GeoTransform =

  -6093465.39, 598.787949, 0.131915338

  3726966.12, 0.00483412443, -598.773963

Metadata:

  AREA_OR_POINT=Area

  TIFFTAG_DATETIME=2016:01:15 12:18:30

  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)

  TIFFTAG_SOFTWARE=GPL Ghostscript 8.71

  TIFFTAG_XRESOLUTION=350

  TIFFTAG_YRESOLUTION=350

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=PIXEL

Corner Coordinates:

Upper Left  (-6093465.390, 3726966.120) (110d36'22.92"W, 42d27'14.89"N)

Lower Left  (-6092183.305,-2092518.026) ( 75d17'32.97"W,  0d 5'45.05"S)

Upper Right ( 2102744.056, 3727032.289) ( 13d35' 4.81"E, 62d39'59.63"N)

Lower Right ( 2104026.141,-2092451.857) (  4d55'50.37"W, 11d15'50.13"N)

Center      (-1994719.624,  817257.132) ( 48d 5'42.97"W, 38d 4'12.82"N)

Band 1 Block=352x352 Type=Byte, ColorInterp=Red

Band 2 Block=352x352 Type=Byte, ColorInterp=Green

Band 3 Block=352x352 Type=Byte, ColorInterp=Blue

Nate:qbert$ gdal_translate -expand rgba -co COMPRESS=LZW /Users/Nate/Downloads/narc_tif/narc.tif convertedNARC.tif

Input file size is 13688, 9719

Error : band 1 has no color table

1 Answer

0 votes
answered Jun 20, 2016 by ba3user (4,000 points)

These images are 3 band RGB. There are no color tables for them. 

First: Warp them to a non-skewed projection with this command:

gdal_warp -t_srs EPSG:4326 ENR_L01.tif foo.tif

Second: Create a VRT that adds an empty alpha band (so they are indeed RGBA) with this command:

gdal_translate -of vrt -b 1 -b 2 -b 3 -b mask foo.tif foo.tif_RGBA.vrt

Third process that VRT with AltusRaster:

AltusRaster -lk 6F00BBCC-7721-4A19-B964-4D2B975AB5DD -ml 11 -i foo.tif_RGBA.vrt -of foo


You may need to modify the max level, I just chose 11 because that is typical for sectionals.

...