
<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Dear mailing list,<div><br></div><div>Apologies I am a new user with GRASS and have searched everywhere online for help with this but to no avail...</div><div><br></div><div><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">I'm trying to run the i.atcorr module on GRASS (version 6.4.3 on a Macbook Pro, Mavericks OSX) with Landsat8 imagery.</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">Parameters didn't seem to be loading from the provided metadata file (MTL_txt) (basically what happened is that when the code was executed it wasn't picking up any of the parameters from the text file), so I went back to the GRASS manuals and tried to input values by hand, here is what I used:</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">8 = geometrical conditions=Landsat ETM+</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">05 04 03.81 98.117 23.131 = month day hh.ddd longitude latitude ("hh.ddd" is in GMT decimal hours)</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">1 = atmospheric mode=tropical</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">0 = aerosols model=none (ASSUMING NO AEROSOLS)</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">0 = visibility [km]</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">0.732 = aerosol optical depth 550 nm</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">-1.408 = mean target elevation above sea level [km]</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">-1000 = sensor on board a satellite</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">0 = THIS IS WHAT I'M WORRIED ABOUT: There is no pre-defined coding for Landsat 8 satellites and I am unsure how to define my own spectral conditions??</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">I checked I was in the correct region (i.e. used g.region) and I've tried using both radiance and reflectance at sensor values (computed using i.landsat.toar, which seems to have worked nicely) but the rasters returned are empty :(</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);"><strong style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">Here is what I'm running (in verbose)</strong></p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">i.atcorr -a --overwrite --verbose iimg=Part1toar_B1@Charlie icnd=/Users/charlottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/.tmp/Charlie.local/819.2 oimg=P1toar_atcorrB1 * <strong style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;"><em style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">*</em>**<em style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">*</em>**<em style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">*</em>**<em style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">*</em>**<em style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">*</em>***</strong> 6s version 4.2b</p><hr style="border: 0px; color: rgb(221, 221, 221); background-color: rgb(221, 221, 221); height: 1px; margin-bottom: 20px; font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; font-size: 14px; line-height: 18px;"><ul style="margin: 0px 0px 1em 30px; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; list-style-position: initial; list-style-image: initial; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);"><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">geometrical conditions identity *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">------------------------------- *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">etm+ observation * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">month: 5 day: 4 *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">solar zenith angle: 23.55 deg solar azimuthal angle: 103.24 deg *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">scattering angle: 156.45 deg azimuthal angle difference: 103.24 deg * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">atmospheric model description *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">----------------------------- *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">atmospheric model identity : *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">aerosols type identity : *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">no aerosols computed * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">spectral condition *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">------------------ *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">constant *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">value of filter function : *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">wl inf=27640754752056197120.000 mic wl sup= 0.000 mic * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">target type *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">----------- *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">homogeneous ground *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">constant reflectance over the spectra 0.000 * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">target elevation description *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">---------------------------- *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">ground pressure [mb] 862.22 *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">ground altitude [km] 1.408 *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">gaseous content at target level: *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">uh2o= 2.056 g/cm2 uo3= 0.244 cm-atm * * *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">atmospheric correction activated *</li><li style="margin: 0px 0px 7px; padding: 0px; border: 0px; vertical-align: baseline; line-height: 16px; word-wrap: break-word;">-------------------------------- * Atmospheric correction... Atmospheric correction complete. (Fri Dec 6 12:16:10 2013) Command finished (22 sec)</li></ul><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;"><strong style="margin: 0px; padding: 0px; border: 0px; vertical-align: baseline;">And here is the output:</strong></p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">| Layer: P1toar_atcorrB1@Charlie Date: Fri Dec 6 12:12:58 2013 | | Mapset: Charlie Login of Creator: charlottewhitha | | Location: NGH | | DataBase: /Users/charlottewhitham/Documents/PhD/My research/my GIS/grassda | | Title: ( P1toar_atcorrB1 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: CELL | | Rows: 7381 | | Columns: 7551 | | Total Cells: 55733931 | | Projection: UTM (zone 47) | | N: 2666715 S: 2445285 Res: 30 | | E: 522315 W: 295785 Res: 30 | | Range of data: min = NULL max = NULL | | | | Data Description: | | generated by i.atcorr | | | | Comments: | | i.atcorr -a iimg="Part1toar_B1@Charlie" iscl=0,255 icnd="/Users/char\ | | lottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/\ | | icnd_lsat1.rtf" oimg="P1toar_atcorrB1" oscl=0,255</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">As you can see the range of data is NULL to NULL... not good.</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">I would be very grateful if somebody could please help me with this - I'm really unsure what to try next and I'm running out of time!</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6);">Thank you very much in advance and please let me know if I need to provide more info,</p><p style="margin: 0px 0px 1em; padding: 0px; border: 0px; font-size: 14px; vertical-align: baseline; clear: both; color: rgb(51, 51, 51); font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif; line-height: 18px; background-color: rgba(248, 248, 248, 0.6); position: static; z-index: auto;">Cheers, Charlie</p></div></body></html>