Okay, here is what I have done in the past. I haven't done this in a while (e.g., haven't tested it with ArcGIS 9.3, only 9.2), so some steps might have changed. Also, you need ISIS, and since it is free, I would recommend downloading it. You're going to need ISIS for so many other projects anyway that it just makes sense to have it. The BTR values from THMPROC are real and valid. I have tested it by producing my own BTR maps in IDL and get the same values. I hope this helps . . . . IMPORTING THEMIS BAND 9 TEMPERATURE MAPS INTO ARCGIS:1. Order projected, Band 9 ISIS cube from THMPROC (THEMIS Processing Web Interface)
a. Go to http://thmproc.mars.edu/index.php
(you can also get to it from http://themis.asu.edu/thmproc
b. You will need to create a username and password. These are used to log in and check on the status of your order and to retrieve your order when it is finished.
c. Click “New Batch”.
d. Enter the specific THEMIS IR images that you want to process.
e. Enter a job description that will jog your memory about what you ordered.
f. Under “Standard Processing”, make sure to click ALL the boxes (including “Unrectify”). If you don’t click “Unrectify”, your projection will be messed up.
g. Enter the projection you want, and remember that ArcGIS works under a -180:180 longitude system.
h. Cropping is a good (and sometimes necessary) idea, especially if you are ordering polar projected cubes, which end up being huge files. If you don’t crop, your order might not be processed.
i. For “Output”, you generally want just “Brightness Temperature” (ordering more will seriously delay the processing time). Under “Brightness Temperature”, click “32 bit ISIS cube”. Then choose “Band 9” (don’t worry about “Backplanes”).
j. Click “Submit Query”. The rest is straightforward.
k. You will have to check on the progress of your order by logging back in periodically (your order name and status will appear when you log in, and you just click on it to see the product). You won’t receive any kind of notification when the order is processed. When the order is complete, you will be able to preview it to see if it looks okay (this allows you to skip the same check in ISIS qview).
l. Save your cube (.QUB) somewhere that is accessible by ISIS.2. Create a raster file and worldfile from your ISIS cube.
a. To do this, you need the isis2world.pl perl script from the USGS.
• Go to the USGS PIGWAD ftp site http://webgis.wr.usgs.gov/pigwad/tutorials/scripts/perl.htm
• Download isis2world.pl. It will bring up the script in the browser window. Copy-paste it to a text file and make sure to save it with the .pl extension (no .txt extension).
• Obviously, you need perl up and running on your machine, so talk to your sys admin if you don’t know how to do this (too much extra detail for this post).
b. In the ISIS environment, run the perl script on the cube.
• >> perl isis2world.pl –e filename.QUB (the –e tag creates a raw file and is the best option for import into ArcGIS)
c. The perl script will generate 2 files: filename.raw and filename.rww (the world file). filename.raw points to the .QUB for its data.
d. Open the .raw file in a text editor and make sure that the byte order is LSB (not MSB). Change it to LSB if necessary and re-save. If you don’t so this, the USGS tool you need below won’t work.
3. Import the raw file into ArcMap.
a. Create a file in ArcCatalog called themisimagenumber_dataset (obviously you put in the image number and dataset)
b. Copy the .raw, .rww, and .QUB files into this folder (you need the .QUB since the raw file points to .QUB for the data)
c. Add the .raw file to your ArcMap project. It may look really strange at this point, but it should be in the right place.
d. You will need the USGS Image Toolbox v1 for ArcMap 9 to make the raster not look strange anymore.
• Search for this toolbox and download it from the ESRI scripts website (http://arcscripts.esri.com
). Save the script with your other ArcMap toolboxes.
• Add the toolbox to ArcMap by right-clicking on “ArcToolbox” in the toolbox area of ArcMap ? choose “Add Toolbox”.
e. Run the “Covert 32bit NODATA” tool from the USGS Image Toolbox to covert the weird values. Leave the “Expression” box at its default value. Save the “output raster” in the themisimagenumber_dataset folder you created above.
f. Voila! Now it should look like a normal, B/W raster with a scale that corresponds to temperature. You can of course change the color scale. And you can get rid of the black space around it by setting the “nodata” pixels to clear (Under the “Display” tab in the “Layer Properties” box, which you get by double-clicking on the layer).