<div dir="ltr"><div>Hi Robi, <br></div><div><br></div><div>My bad, region was wrong... anyway, trying to replicate from the beginning I see that i.sentinel.import with -r is importing empty maps. However, the jp2 files have data. I can visualize them in QGIS. I'll report in a different thread. Once I can figure it out and reach this point again, I'll report back to you. Thanks for the clarification!</div><div><br></div><div>Cheers<br></div>Vero<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mié., 29 ene. 2020 a las 9:08, roberta fagandini (<<a href="mailto:robifagandini@hotmail.it">robifagandini@hotmail.it</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Vero,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
i.sentinel.mask doesn't retrieve the maximum from the xml file but computes the value from the image using r.univar. The error could be due to the computational region..has the region been set correctly?</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Can you send the <span style="color:rgb(50,49,48);font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255);display:inline">MTD_TL.xml file? I can download the image and make some tests</span><span style="color:rgb(50,49,48);font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255);display:inline">.</span></div>
<div>
<div id="gmail-m_7941873122284512215appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Roberta</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_7941873122284512215divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>Da:</b> grass-user <<a href="mailto:grass-user-bounces@lists.osgeo.org" target="_blank">grass-user-bounces@lists.osgeo.org</a>> per conto di Veronica Andreo <<a href="mailto:veroandreo@gmail.com" target="_blank">veroandreo@gmail.com</a>><br>
<b>Inviato:</b> martedì 28 gennaio 2020 19:57<br>
<b>A:</b> grass-user <<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>><br>
<b>Oggetto:</b> [GRASS-user] error in i.sentinel.mask: nan values in xml file</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>Hello</div>
<div><br>
</div>
<div>I try to use i.sentinel.mask but I get the following error:</div>
<div><br>
</div>
<div>i.sentinel.mask blue=$blue green=$green red=$red swir11=$swir11 nir=$nir swir12=$swir12 nir8a=$nir8a cloud_mask=cloud shadow_mask=shadow mtd_file=`find $out_dir -name 'MTD_TL.xml'`</div>
<div><br>
ADVERTENCIA: All subsequent operations will be limited to the current<br>
computational region<br>
ADVERTENCIA: Any rescale factor has been applied<br>
--- Start computing maximum values of bands ---<br>
T21JYM_20200126T134201_B02_10m<br>
T21JYM_20200126T134201_B03_10m<br>
T21JYM_20200126T134201_B04_10m<br>
T21JYM_20200126T134201_B08_10m<br>
T21JYM_20200126T134201_B8A_20m<br>
T21JYM_20200126T134201_B11_20m<br>
T21JYM_20200126T134201_B12_20m<br>
--- Computed maximum value: dict_values([nan, nan, nan, nan, nan, nan,<br>
nan]) ---<br>
--- Statistics have been computed! ---<br>
--- Start clouds detection procedure ---<br>
--- Computing cloud mask... ---<br>
Mapa no válido <nan><br>
Parse error<br>
ERROR: error al parsear<br>
ERROR: An error occurred while running r.mapcalc with expression:<br>
cloud_def82115 = if((((T21JYM_20200126T134201_B02_10m > (0.08*nan))<br>
&& (T21JYM_20200126T134201_B03_10m > (0.08*nan)) &&<br>
(T21JYM_20200126T134201_B04_10m > (0.08*nan))) == 1) &&<br>
(((T21JYM_20200126T134201_B04_10m < ((0.08*nan)*1.5)) &&<br>
(T21JYM_20200126T134201_B04_10m ><br>
T21JYM_20200126T134201_B12_20m*1.3)) == 0) &&<br>
(((T21JYM_20200126T134201_B11_20m < (0.1*nan)) &&<br>
(T21JYM_20200126T134201_B12_20m < (0.1*nan))) == 0) &&<br>
((if(T21JYM_20200126T134201_B8A_20m ==<br>
max(T21JYM_20200126T134201_B8A_20m, 2 *<br>
T21JYM_20200126T134201_B02_10m, 2 * T21JYM_20200126T134201_B03_10m,<br>
2 * T21JYM_20200126T134201_B04_10m))) == 0) &&<br>
((T21JYM_20200126T134201_B02_10m > 0.2) == 1), 0, null())<br>
</div>
<div><br>
</div>
<div>I checked the xml file and there are a lot of nan's yes, but there are also real values. Why doesn't the max function in this part:
<a href="https://github.com/OSGeo/grass-addons/blob/master/grass7/imagery/i.sentinel/i.sentinel.mask/i.sentinel.mask.py#L321" target="_blank">
https://github.com/OSGeo/grass-addons/blob/master/grass7/imagery/i.sentinel/i.sentinel.mask/i.sentinel.mask.py#L321</a> recognize them? Isn't there a way to avoid nan while computing max?<br>
</div>
<div><br>
</div>
<div>Thanks much in advance</div>
<div>Vero</div>
<div><br>
</div>
<div>ps: I can provide the xml if needed. <br>
</div>
</div>
</div>
</div>
</div>
</blockquote></div>