Example: Import ERDAS GIS with GDAL and Assign Coordinate System

Manifold has the built-in ability to connect to, link, import and read/write to hundreds of different databases, file databases, file formats, web servers and similar data sources using Manifold's own internal code, with no need to install any additional software.  

 

In addition, Manifold includes a dataport that can launch the GDAL/OGR libraries to connect to any format those libraries support.   

 

Both the Manifold collection of data sources and the GDAL/OGR libraries collection of file formats are large collections so there is considerable overlap between the two.   In many cases there is a native Manifold dataport for a data source and there is also a module available in the GDAL/OGR collections for that format.   However, there are also data sources which Manifold supports that GDAL/OGR does not, and there are formats available within GDAL/OGR that Manifold does not include as a built-in to Manifold.

 

In this example we will use GDAL to get data from an ERDAS .gis file format that encodes a palette image.  Manifold can also read such images with the native ERDAS .gis dataport built into Manifold.   Neither Manifold nor GDAL can import such images with correct coordinate system, since the .gis format fails to provide full coordinate system information.   In both cases we must manually assign the initial coordinate system.

 

 To use GDAL or OGR we must install GDAL in Windows, and then we can use it from Manifold.  If we do not install GDAL, we will not be able to use the capabilities described in this topic. See the the GDAL / OGR topic for installation instructions. The Manifold GDAL dataport allows working with GDAL versions set forth in the Third Party Release Levels topic, and automatically selects the latest version available with automatic adjustments for the call interface.

Import an ERDAS .gis Format File using GDAL

ERDAS is a raster-capable package introduced in 1978 and eventually acquired by Hexagon.  ERDAS .gis format files may be encountered on various government web sites.   In this example we import an ERDAS .gis format file that shows Thematic Mapper imagery acquired by the US Landsat satellite, the file being downloaded from

http://caddolakedata.us/data/general-caddo-data/tm-imagery-1993 with the following information provided as a description:

 

Info: Thematic Mapper (TM) Satellite Imagery - 1993

Scale/Resolution: 30m pixels

Data Type: raster layer, Erdas v7.5 GIS File

Source: U.S. Geological Survey (USGS), National Wetlands Research Center

Description: TM satellite imagery from January 30, 1993. The original TM scene was subsetted into 3 bands (4,5, & 3). The subsetted file was then processed using the RGBCLUS routine in Erdas v7.5. The resulting file is an 8-bit GIS file. Bands 4,5,3 were chosen because of the usefulness of those bands in environmental management.

 

In addition, the website's data page reports:

 

The various raster layers are also in UTM, Zone 15. Each raster layer is only available in NAD27 (CLARKE 1866)

 

To fit into this documentation, illustrations show an artificially small Manifold desktop, with only a few panes, docked to the right side.   In real life we use a much larger Manifold desktop, and all panes would be turned on, with some panes docked to the left and others docked to the right.

 

We launch Manifold.

 

 

 

Choose File - Import.

 

 

In the Import dialog, navigate to the location of the desired file, as downloaded from the source web site.  Choose .Files (GDAL/OGR) (*.*) for the file type.

 

 

Click on the .gis file to load it into the File name box and then press the Open button.

 

 

 

 

The image immediately imports into Manifold, creating a new image and the image's table in the Project pane.  We can double-click open the image to see it in an image window.

 

 

 

 

Choosing the Info pane we see the initial coordinate system displayed in red color using the Pseudo Mercator placeholder default name.   That warns us we must assign an initial coordinate system.

 

 We click the coordinate system picker button, choose Assign Initial Coordinate System and then More... to launch the Coordinate System dialog.

 

 

 

 

From the commentary in the web site, we must assign UTM zone 15.   We click the Standard tab and then we enter universal into the filter box to reduce the incredibly long list of Standard projections to a shorter list.  We scroll down and click on Universal Transverse Mercator Zone 15 (N) to choose that coordinate system.  We know from basic GIS knowledge that UTM is an abbreviation for Universal Transverse Mercator.   The web site did not say the intended Zone 15 is North or South, but we assume from the obviously North American focus of the web site they mean North, so we choose Zone 15 (N) and not Zone 15 (S).   

 

If we have an eye for details, we will note that the JSON description for the chosen coordinate system says it uses the WGS84 datum.   That is wrong for our data, which the originating web site says uses the NAD27 datum.   That is not a problem.   For now we will accept the coordinate system as is, and then in the next step we will change the datum by using the Repair command.

 

Press OK.

 

 

 

 

The new UTM Zone 15 (N) coordinate system appears in the Info pane readout.    However, we know that the UTM Zone 15 (N) system as defined uses the wrong datum so we must change the datum.   

 

 We will change the datum by once again clicking the coordinate system picker button, choosing Repair Initial Coordinate System, and then choosing More... to launch the Coordinate System dialog.

 

 

 

 

In the Coordinate System dialog we click the Custom tab.

 

We press the base coordinate system picker button, and then choose More... to launch the Base Coordinate System dialog.

 

 

 

 

In the Standard tab we want to find something that looks like NAD27, that is, the North American Datum for 1927.   we enter 27 into the filter box to reduce the list to fewer entries.    The originating web site does not say exactly what they mean by "NAD27," but it is probably a safe guess that they mean the North American 1927 datum for the Continental United States, CONUS.   We click on North American 1927 (mean for CONUS) to select that datum and then we press OK.

 

 

 

 

Back in the Coordinate System dialog we press OK.   Done!

 

 

 

 

The image now has the correct initial coordinate system assigned.  It will exactly line up with other, known-good layers such as Bing.   In the illustration above we first created a map using a Bing streets image server as a base layer, and then we dragged and dropped the imported image into the map, zooming it to see how well the image aligns with Bing.   It aligns very well.

 

We float (that is, undock) the map (Shift-click on the Map tab)  so we can resize it to be bigger, to make better use of the space we have available in this documentation for illustrations.

 

 

 

 

We can zoom further into the view to look more closely at alignment between the two layers, noting how roads and rivers align at the edges.

 

 

 

 

In the illustration above we have added a layer above the image that is a data source with Google streets in "transparent" style, where pixels in between the streets and labels are transparent.   That makes it clear that our initial coordinate system assignment is accurate, as all features in the layers align very well with the image.

Problems with GDAL imports

If nothing appears when importing using GDAL, launch the Log Window to see if there are any error messages.   For example, when importing an ESRI .e00 file fails, in the Log Window you might see an error message such as:

 

*** (import) This looks like a compressed E00 file and cannot be processed directly. You may need to uncompress it first using the E00compr library or the e00conv program.

 

Error messages in the Log Window might report other problems, such as the necessary GDAL .dll files not being found.  That usually happens when the PATH environment variable does not include the GDAL installation folders.  See the discussion in the GDAL / OGR topic on setting the PATH environment variable.

Notes

Projections - Take care to find any available commentaries or metadata information that may accompany the .gis files in use.   .gis files do not provide full projection information and will result in imports that do not correctly specify projection information.  In that case, we will have to use whatever documentation we can find on the data set to learn what the correct projection should be so we can specify it using Assign Initial Coordinate System.

 

GDAL is not Manifold - When connecting to GDB using GDAL/OGR we must be aware we are no longer using Manifold code but instead are using GDAL code.  GDAL has earned a good reputation.  It is a tremendous advantage to use GDAL's very broad reach of modules from Manifold to connect to many niche or nearly-extinct formats, and GDAL does a great job at that.  However, GDAL's code in general is not as bulletproof as the Radian technology used in Manifold.  Connections using GDAL fall outside of Manifold's reputation for never crashing.      

See Also

Images

 

Info Pane: Component

 

Assign Initial Coordinate System

 

Repair Initial Coordinate System

 

Coordinate System

 

Base Coordinate System

 

GDAL / OGR

 

Style

 

Palette Images

 

Example: Spectacular Images and Data from Web Servers

 

Example: An Imageserver Tutorial