<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>first lines of test.csv<br></div><div><br></div><div>cat|area|perimeter|compact_square|compact_circle|fd</div><div>1|11876.000000|1398.000000|0.311808|3.618821|1.543932</div><div>2|635805.000000|57572.000000|0.055400|20.367814|1.640512</div><div>3|25898.000000|3322.000000|0.193773|5.823198|1.595825</div><div>4|7767.000000|1702.000000|0.207122|5.447885|1.661054</div><div>5|13605.000000|2200.000000|0.212074|5.320696|1.617158</div><div>6|442366.000000|44230.000000|0.060150|18.759495|1.645730</div><div><br></div><div>--</div><div>Windows 7<br>Dell</div><div>Proc. Intel(R) Core(TM) i7-6700 CPU @ 3.40GHz<br></div><div>RAM: 16 GB</div><div>--</div><div>Could the "ju" be caused my the Mask?</div><div>--</div><div>I used uspo on 3 regions to decide on parameters. 0.24 as treshold and minsize 3500 looked to be the best combination of parameters studying all the results in the output file. <br></div><div>Summary:</div><div><div>regions=regionAB@LUP1,regionD@LUP1,regionRM@LUP1 segmentation_method=region_growing threshold_start=0.18 threshold_stop=0.25 threshold_step=0.02 minsize_start=3000 minsize_stop=6000 minsize_step=500 optimization_function=f number_best=2 memory=8000 processes=4</div><div>Working on region regionAB@LUP1</div><div>Working on region regionD@LUP1</div><div>Working on region regionRM@LUP1</div><div>Best values:</div><div>Region<span style="white-space:pre"> </span>Thresh<span style="white-space:pre"> </span>Minsize<span style="white-space:pre"> </span>Optimization</div><div>regionRM@LUP1<span style="white-space:pre"> </span>0.24<span style="white-space:pre"> </span>3500<span style="white-space:pre"> </span>0.7951456323</div><div>regionRM@LUP1<span style="white-space:pre"> </span>0.24<span style="white-space:pre"> </span>3000<span style="white-space:pre"> </span>0.765849594003</div><div>regionAB@LUP1<span style="white-space:pre"> </span>0.24<span style="white-space:pre"> </span>4000<span style="white-space:pre"> </span>0.630706501774</div><div>regionAB@LUP1<span style="white-space:pre"> </span>0.24<span style="white-space:pre"> </span>3500<span style="white-space:pre"> </span>0.472619083853</div><div>regionD@LUP1<span style="white-space:pre"> </span>0.2<span style="white-space:pre"> </span>3000<span style="white-space:pre"> </span>0.802618475115</div><div>regionD@LUP1<span style="white-space:pre"> </span>0.22<span style="white-space:pre"> </span>3000<span style="white-space:pre"> </span>0.77939730321</div></div><div><br></div><div>I also had other uspo runs but with less useful results. </div><div><br></div><div>The map with segments was created (0.24 treshold minsize 3500). But after the errors with i.segment.stats I reasoned the number of segments (18685) was too much to handle as it failed to create the attribute table. To reduce the number of segments, I used this map as seed to create a new segmentation map with the parameters 0.25 and 5000 and a mask to eliminate the matrix. </div><div>But now I see the number of segments shouldn't have been a problem in this case.</div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Oct 9, 2018 at 9:35 AM Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On 09/10/18 13:48, Jamille Haarloo wrote:<br>
> Hello Moritz,<br>
> <br>
> First a pop-up came about not finding it, so reinstalled it.<br>
> <br>
> g.extension r.object.geometry<br>
> WARNING: Extension <r.object.geometry> already installed. Re-installing...<br>
> Downloading precompiled GRASS Addons <r.object.geometry>...<br>
> Updating addons metadata file...<br>
> Installation of <r.object.geometry> successfully finished<br>
> (Mon Oct 08 16:55:04 2018) Command finished (4 sec)<br>
> (Mon Oct 08 16:55:53 2018)<br>
> r.object.geometry SegmW24IDM4DV4 output=test.csv<br>
> Calculating statistics<br>
> Writing output<br>
> (Mon Oct 08 16:56:14 2018) Command finished (21 sec)<br>
<br>
Could you provide the first lines of content of test.csv ?<br>
<br>
> <br>
> g.region -p<br>
> projection: 1 (UTM)<br>
> zone: 21<br>
> datum: wgs84<br>
> ellipsoid: wgs84<br>
> north: 634826.41313284<br>
> south: 611477.61329757<br>
> west: 775494.11663193<br>
> east: 793630.22784822<br>
> nsres: 0.50000642<br>
> ewres: 0.50000307<br>
> rows: 46697<br>
> cols: 36272<br>
> cells: 1693793584<br>
<br>
<br>
So you have over 1.5 billion cells. That's quite big, but not impossible <br>
to handle. Depends on your machine.<br>
<br>
> <br>
> <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> <<a href="http://r.info" rel="noreferrer" target="_blank">http://r.info</a>> SegmW24IDM4DV4<br>
> +----------------------------------------------------------------------------+<br>
> | Map: SegmW24IDM4DV4 Date: Fri Oct 05 23:46:04 <br>
> 2018 |<br>
> | Mapset: LUP1 Login of Creator: haarlooj <br>
> |<br>
> | Location: Cottica_LUP <br>
> |<br>
> | DataBase: Z:\GRASS_GIS <br>
> |<br>
> | Title: <br>
> |<br>
> | Timestamp: none <br>
> |<br>
> |----------------------------------------------------------------------------|<br>
> | <br>
> |<br>
> | Type of Map: raster Number of Categories: 0 <br>
> |<br>
> | Data Type: CELL <br>
> |<br>
> | Rows: 46697 <br>
> |<br>
> | Columns: 36272 <br>
> |<br>
> | Total Cells: ju <br>
<br>
<br>
This 'ju' is weird, here.<br>
<br>
<br>
> |<br>
> | Projection: UTM (zone 21) <br>
> |<br>
> | N: 634826.41313284 S: 611477.61329757 Res: <br>
> 0.50000642 |<br>
> | E: 793630.22784822 W: 775494.11663193 Res: <br>
> 0.50000307 |<br>
> | Range of data: min = 1 max = 12901 <br>
> |<br>
> | <br>
> |<br>
> | Data Description: <br>
> |<br>
> | generated by i.segment <br>
> |<br>
> | <br>
> |<br>
> | Comments: <br>
> |<br>
> | i.segment --overwrite group="DV4@LUP1,IDM4@LUP1,W2@LUP1,W4@LUP1" <br>
> out\ |<br>
> | put="SegmW24IDM4DV4" threshold=0.25 radius=1.5 <br>
> method="region_growin\ |<br>
> | g" similarity="euclidean" minsize=5000 memory=10000 <br>
> seeds="SegW24IDM\ |<br>
> | 4DV4@LUP1" goodness="SegmW24IDM4DV4_goodn"<br>
<br>
<br>
You only have 12901 segments, so an average size of 32822m2 per segment. <br>
Is this what you aimed for ? As you asked for a minsize of 5000 pixels I <br>
guess this is deliberate.<br>
<br>
You threshold of 0.25 is also very high...<br>
<br>
Moritz<br>
</blockquote></div>