GEOTIFF generation from BYU-MERS "SIR" image format

The BYU-MERS "sir" image format was developed by the Brigham Young University (BYU) Microwave Earth Remote Sensing (MERS) research group to store images of the earth along with the information required to easily earth-locate the image pixels. The binary SIR file includes a header with all the information to geolocate each pixel in the image. The origin of the .SIR file image is in the lower left corner. The earth location of a pixel is identified with its lower-left corner. Most sir files are in either polar stereographic or Lambert equal-area projections, though some EASE grid images are regularly produced.

Two experimental procedures have been developed to generate GEOTIFF-formatted files from .SIR files. One is based on compiling a custom program. The second, developed by Doug Nebert, is described briefly below .

The two packages are used to convert the GIF to GeoTIFF (or other formats) are the GDAL transformation library and its utility programs and the PROJ.4 projection library used by GDAL to perform the projections. These programs were run on Windows platform in command-line form though they can be used in the Unix/Linux environments.PROJ.4: http://www.remotesensing.org/proj/
Contains the public domain libraries for the conversion of hundreds of known projections and the ability to pass parameters to create custom projections. Download and install this such that the projection libraries are installed as public on the system (or in windows, be sure the proj4.dll is in the WINDOWS directory).

GDAL: http://trac.osgeo.org/gdal/wiki/DownloadingGdalBinaries
This page allows for the download of the GDAL software for various platforms. Download and install as per instructions to allow the execution of GDAL commands as part of the PATH variable.

Once these two programs are installed (preferably in this order) one can execute the following command in the command shell to convert images to other formats, or in so doing, insert coordinate reference details into the header of capable files (e.g. GeoTIFF). For projections not well-known to PROJ.4 the parameters are passed to assign the output properties.

First, create a full scale GIF image from the SIR file. Use a one-to-one pixel version. Then, run gdal_translate to create a .tif file that has the geotiff header. You will have to input the the corner coordinates (in meters) of the SIR image. This can be done by using a SIR viewer tool and looking at the corner positions.

Consider conversion of a QuickSAT slice Arctic image (qush-a-ArcYY-ddd-ddd.sir). Examinging the .SIR file header, .SIR image is in polar stereographic form with reference latitude 70 deg and vertical longitude -45 deg. The image corner coordinates are 3400000 3400000 3400000 -3400000. (While all q*s*-*-Arc* have the same header, images of other areas will have different corners and projections). A full scale .gif image is generated, then gdal_translate is run.

gdal_translate -a_srs "+datum=wgs84 +proj=stere +lat_ts=70.0 +lat_0=90 +lon_0=-45" -a_ullr -3400000 3400000 3400000 -3400000 -of GTiff -co compress=LZW INPUT.gif OUTPUT.tif where
-a_srs (assigned Spatial Reference System, followed by proj4 params)
+datum (wgs84)
+proj (projection stereographic method)
+lat_ts (latitude of true scale 1:1)
+lat_0 (the latitude of origin)
+lon_0 (the longitude of the origin)
-a ullr (the projected coordinates of the upper left and lower right)
-of (output format)
-co (conversion options relevant to a format)

You can then script the batch conversion of images quickly and make practically any format you want. It will perform coordinate transformations to other projections also, like Universal Polar Stereo N (-t_SRS EPSG:32661).