<div dir="ltr">Hello all,<div><div>I am using i.atcorr module of GRASS 7.0.4 release. I compiled the release, as I wanted to add the Sentinel-2 sensor. So far, I was successful in compiling the release. However, when I run atcorr, I get an empty raster file. I have included the reports from <a href="http://r.info">r.info</a>, i.atcorr and the parameters below:</div><div><br></div></div><div><b><a href="http://r.info">r.info</a> output of the image</b></div><div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div><div>(Wed Mar 16 16:38:39 2016)                                                      </div><div><a href="http://r.info">r.info</a> map=S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32TPS_B02@PERMANENT</div><div> +----------------------------------------------------------------------------+</div><div> | Map:      S2A_OPER_MSI_L1C_TL_EPA__2016  Date: Wed Mar 16 16:29:33 2016    |</div><div> | Mapset:   PERMANENT                      Login of Creator: abheeman        |</div><div> | Location: munich                                                           |</div><div> | DataBase: /media/abheeman/data/grassdata                                   |</div><div> | Title:     ( S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32TPS_B02 ) |</div><div> | Timestamp: none                                                            |</div><div> |----------------------------------------------------------------------------|</div><div> |                                                                            |</div><div> |   Type of Map:  raster               Number of Categories: 0               |</div><div> |   Data Type:    CELL                                                       |</div><div> |   Rows:         10980                                                      |</div><div> |   Columns:      10980                                                      |</div><div> |   Total Cells:  120560400                                                  |</div><div> |        Projection: UTM (zone 32)                                           |</div><div> |            N:    5200020    S:    5090220   Res:    10                     |</div><div> |            E:     709800    W:     600000   Res:    10                     |</div><div> |   Range of data:    min = 0  max = <b>23562</b>                                   |</div><div> |                                                                            |</div><div> |   Data Description:                                                        |</div><div> |    generated by r.in.gdal                                                  |</div><div> |                                                                            |</div><div> |   Comments:                                                                |</div><div> |    r.in.gdal -e input="/media/abheeman/data/grassdata/S2A_OPER_PRD_MSIL\   |</div><div> |    1C_PDMC_20160315T115050_R065_V20150826T102655_20150826T102655.SAFE/G\   |</div><div> |    RANULE/S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32TPS_N02.0\   |</div><div> |    1/IMG_DATA/S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32TPS_B\   |</div><div> |    02.jp2" output="S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32\   |</div><div> |    TPS_B02" memory=300 offset=0                                            |</div><div> |                                                                            |</div><div> +----------------------------------------------------------------------------+</div><div>(Wed Mar 16 16:38:39 2016) Command finished (0 sec)   </div></div></blockquote></div><div><div><b><br></b></div><div><b>i.atcorr verbose return</b></div></div><div><br></div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div>(Wed Mar 16 16:39:48 2016)                                                      </div><div>i.atcorr -i -r --overwrite --verbose input=S2A_OPER_MSI_L1C_TL_EPA__20160315T065417_A000920_T32TPS_B02@PERMANENT range=0,10000 parameters=/media/abheeman/data/grassdata/conf output=munich_test_atcorr</div><div>* ****************************** 6s version 4.2b</div><div>****************************** *</div><div>*                        geometrical conditions identity</div><div>*</div><div>*                        -------------------------------</div><div>*</div><div>*                        Sentinel 2  observation</div><div>*</div><div>*</div><div>*</div><div>*    month: 8 day: 26</div><div>*</div><div>*    solar zenith angle:   39.16 deg  solar azimuthal angle: 161.21 deg   </div><div>*</div><div>*    view zenith angle:     0.00 deg  view azimuthal angle:0.00 deg   </div><div>*</div><div>*    scattering angle:    140.84 deg  azimuthal angle difference: 161.21 deg   </div><div>*</div><div>*</div><div>*</div><div>*                        atmospheric model description</div><div>*</div><div>*                        -----------------------------</div><div>*</div><div>*            atmospheric model identity :</div><div>*</div><div>*                tropical            (uh2o=4.12g/cm2,uo3=.247cm-atm)            </div><div>*</div><div>*</div><div>*</div><div>*            aerosols type identity :</div><div>*</div><div>*                    Continental aerosols model</div><div>*</div><div>*</div><div>*</div><div>*            optical condition identity :</div><div>*</div><div>*                 visibility :   13.00 km  opt. thick. 550nm:   0.3516       </div><div>*</div><div>*</div><div>*</div><div>*                        spectral condition</div><div>*</div><div>*                        ------------------</div><div>*</div><div>*            sentinel_2 band1</div><div>*</div><div>*                value of filter function :</div><div>*</div><div>*                 wl inf=    0.300 mic   wl sup=    2.600mic                  </div><div>*</div><div>*</div><div>*</div><div>*                        target type</div><div>*</div><div>*                        -----------</div><div>*</div><div>*            homogeneous ground</div><div>*</div><div>*              constant reflectance over the spectra 0.000                 </div><div>*</div><div>*</div><div>*</div><div>*                        target elevation description</div><div>*</div><div>*                        ----------------------------</div><div>*</div><div>*            ground pressure  [mb]        954.77</div><div>*</div><div>*            ground altitude  [km]        0.520</div><div>*</div><div>*                 gaseous content at target level:</div><div>*</div><div>*                 uh2o=    3.227 g/cm2        uo3=    0.246 cm-atm             </div><div>*</div><div>*</div><div>*</div><div>*                         atmospheric correction activated</div><div>*</div><div>*                         --------------------------------</div><div>*</div><div>Atmospheric correction...</div><div>Atmospheric correction complete.</div><div>(Wed Mar 16 16:39:59 2016) Command finished (10 sec) <b>          </b></div></blockquote><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div><div><br></div></div></blockquote><div><div><b>Parameters used</b></div></div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div><div>20                            - geometrical conditions=Sentinel-2</div></div><div><div>08 26 10.48 11.353693 48.200905 - month day hh.ddd longitude latitude ("hh.ddd" is in decimal hours GMT)</div></div><div><div>1                            - atmospheric mode=tropical</div></div><div><div>1                            - aerosols model=continental</div></div><div><div>13                           - visibility [km] (aerosol model concentration)</div></div><div><div>-0.520                       - mean target elevation above sea level [km] (here 600m asl)</div></div><div><div>-786                        - sensor height (here, sensor on board a satellite)</div></div><div><div>129                           - 1nd band of Sentinel-2A</div></div></blockquote><div><br></div><div><b>Note</b></div><div>Before I invest more time in finding process to derive correct AOT and chose appropriate atmospheric model, I wanted to give it a try with fake parameters. As far as I understand, with fake parameter there should be some short of false atmospheric correction. Therefore, there are false values for atmospheric mode, aerosol model and visibility just for test purpose.</div><div><i>From <a href="http://r.info">r.info</a>, maximum value seems to be 23562. But the value in sentinel bands can range only between 0-10000 currently, which brings to the assumption that it probably could be noise. Will it affect the result?</i></div><div><br></div><div><b>Some general questions:</b></div><div><ol><li><i>For 6s parameter, it take geographic coordinate (I guess EPSG:4326) as input. So, I was wondering if we also have transform the coordinate system of the image to the same?</i></li><li><i>In my system, it took 10 seconds to process one image. But, I am not sure if the processing time will remain same when the process outputs non-empty file. My image size is 10980x10980 pixels. What is general processing time it generally takes, for may be single Landsat-8 image ?</i></li><li><i>Regarding the raster with empty output, does anybody have experienced such behavior recently? The discussions I searched mostly date back to the time when Landsat-8 was not supported.</i></li></ol></div><div><br></div><div>If it helps, I am running grass in the computer with Ubuntu 14.04 system, with 8 cores,16 gb ram and Intel Core i7-4790.</div><div><br></div><div>Regards,</div><div>Abhishek</div><br></div>