<div dir="ltr">Below is the log of the commands ran through the GRASS GIS in WINDOWS OS.<div>The commands are to find the volume of the dump pile / cavity in the ground.</div><div>I am using everything in the EPSG:4326 CRS.</div><div><br></div><div>Input Files:</div><div><br></div><div>area_polygon.geojson</div><div>```</div><div>{"type": "Feature", "properties": {}, "geometry": {"type": "Polygon", "coordinates": [[[73.31359828017844, 26.222181594890415], [73.31345181478787, 26.222065660491182], [73.31356935634167, 26.221948621838607], [73.31371520633195, 26.22205185876942], [73.31359828017844, 26.222181594890415]]]}}<br></div><div>```</div><div><br></div><div>polygon_points.geojson</div><div>```</div><div>{"type": "FeatureCollection", "features": [{"type": "Feature", "properties": {}, "geometry": {"type": "Point", "coordinates": [73.31359828017844, 26.222181594890415]}}, {"type": "Feature", "properties": {}, "geometry": {"type": "Point", "coordinates": [73.31345181478787, 26.222065660491182]}}, {"type": "Feature", "properties": {}, "geometry": {"type": "Point", "coordinates": [73.31356935634167, 26.221948621838607]}}, {"type": "Feature", "properties": {}, "geometry": {"type": "Point", "coordinates": [73.31371520633195, 26.22205185876942]}}, {"type": "Feature", "properties": {}, "geometry": {"type": "Point", "coordinates": [73.31359828017844, 26.222181594890415]}}]}<br></div><div>```</div><div><br></div><div>DSM file also in EPSG:4326 CRS, containing elevation information.</div><div><br></div><div>```</div><div>(Mon Jul 26 22:21:05 2021)                                                      <br>v.import input=C:\\Users\\Administrator\\Downloads\\area_file.geojson output=polygon_area --overwrite<br>Check if OGR layer <area_file> contains polygons...<br>Creating attribute table for layer <area_file>...<br>Importing 1 features (OGR layer <area_file>)...<br>-----------------------------------------------------<br>Registering primitives...<br>-----------------------------------------------------<br>Cleaning polygons<br>-----------------------------------------------------<br>Breaking polygons...<br>Breaking polygons (pass 1: select break points)...<br>Breaking polygons (pass 2: break at selected points)...<br>-----------------------------------------------------<br>Removing duplicates...<br>-----------------------------------------------------<br>Breaking boundaries...<br>-----------------------------------------------------<br>Removing duplicates...<br>-----------------------------------------------------<br>Cleaning boundaries at nodes...<br>-----------------------------------------------------<br>Merging boundaries...<br>-----------------------------------------------------<br>Removing dangles...<br>-----------------------------------------------------<br>Building areas...<br>-----------------------------------------------------<br>Removing bridges...<br>-----------------------------------------------------<br>Registering primitives...<br>Building areas...<br>Attaching islands...<br>-----------------------------------------------------<br>Finding centroids for OGR layer <area_file>...<br>-----------------------------------------------------<br>Writing centroids...<br>-----------------------------------------------------<br>1 input polygons<br>Total area: 341.906 (1 areas)<br>-----------------------------------------------------<br>Copying features...<br>Building topology for vector map <polygon_area@PERMANENT>...<br>Registering primitives...<br>Building areas...<br>Attaching islands...<br>Attaching centroids...<br>Input <C:\\Users\\Administrator\\Downloads\\area_file.geojson> successfully imported without reprojection<br>(Mon Jul 26 22:21:06 2021) Command finished (1 sec)                             <br>(Mon Jul 26 22:21:13 2021)                                                      <br>v.import input=C:\\Users\\Administrator\\Downloads\\points_file.geojson output=polygon_points --overwrite<br>Check if OGR layer <points_file> contains polygons...<br>Creating attribute table for layer <points_file>...<br>Importing 5 features (OGR layer <points_file>)...<br>-----------------------------------------------------<br>Building topology for vector map <polygon_points@PERMANENT>...<br>Registering primitives...<br>Input <C:\\Users\\Administrator\\Downloads\\points_file.geojson> successfully imported without reprojection<br>(Mon Jul 26 22:21:14 2021) Command finished (0 sec)                             <br>(Mon Jul 26 22:21:34 2021)                                                      <br>v.buffer input=polygon_area -s type=area output=region distance=0.001 minordistance=1 --overwrite<br>Note: In latitude-longitude coordinate system specify distances in degree unit<br>WARNING: Option 'minordistance' is not available with GEOS buffering<br>Buffering areas...<br>Cleaning buffers...<br>Building parts of topology...<br>Building topology for vector map <region@PERMANENT>...<br>Registering primitives...<br>Snapping boundaries...<br>Reading features...<br>Snap vertices Pass 1: select points<br>Snap vertices Pass 2: assign anchor vertices<br>Snap vertices Pass 3: snap to assigned points<br>Breaking polygons...<br>Breaking polygons (pass 1: select break points)...<br>Breaking polygons (pass 2: break at selected points)...<br>Removing duplicates...<br>Breaking boundaries...<br>Removing duplicates...<br>Cleaning boundaries at nodes<br>Building topology for vector map <region@PERMANENT>...<br>Building areas...<br>Removing dangles...<br>Removing bridges...<br>Attaching islands...<br>Building topology for vector map <region@PERMANENT>...<br>Attaching islands...<br>Calculating centroids for all areas...<br>Generating list of boundaries to be deleted...<br>Deleting boundaries...<br>Calculating centroids for areas...<br>Building topology for vector map <region@PERMANENT>...<br>Registering primitives...<br>Building areas...<br>Attaching islands...<br>Attaching centroids...<br>(Mon Jul 26 22:21:34 2021) Command finished (0 sec)                             <br>(Mon Jul 26 22:21:49 2021)                                                      <br>r.external input=C:\\Users\\Administrator\\Downloads\\dsm2.tif output=dsm --overwrite<br>Reading band 1 of 1...<br>Link to raster map <dsm> created.<br>(Mon Jul 26 22:21:50 2021) Command finished (1 sec)                             <br>(Mon Jul 26 22:21:56 2021)                                                      <br>g.region raster=dsm                                                             <br>(Mon Jul 26 22:21:56 2021) Command finished (0 sec)                             <br>(Mon Jul 26 22:22:00 2021)                                                      <br>g.region vector=region                                                          <br>(Mon Jul 26 22:22:00 2021) Command finished (0 sec)                             <br>(Mon Jul 26 22:22:07 2021)                                                      <br>r.mask vector=region                                                            <br>Reading areas...<br>Writing raster map...<br>All subsequent raster operations will be limited to the MASK area. Removing or renaming raster map named 'MASK' will restore raster operations to normal.<br>(Mon Jul 26 22:22:08 2021) Command finished (1 sec)                             <br>(Mon Jul 26 22:22:14 2021)                                                      <br>v.what.rast map=polygon_points raster=dsm column=height                         <br>Column <height> not found in the table <polygon_points>. Creating...<br>Reading features from vector map...<br>Update vector attributes...<br>v.what.rast complete. 5 records updated.<br>(Mon Jul 26 22:22:14 2021) Command finished (0 sec)                                                      <br>(Mon Jul 26 22:32:57 2021)                                                      <br>v.surf.rst input=polygon_points zcolumn=height elevation=dsm_below_pile --overwrite<br>Reading features from vector map ...<br>WARNING: Strip exists with insufficient data<br>Ignoring 1 points (too dense)<br>WARNING: 4 points given for interpolation (after thinning) is less than given NPMIN=300<br>WARNING: There are less than 600 points for interpolation. No segmentation is necessary, to run the program faster set segmax=600 (see manual)<br>Processing all selected output files will require<br>29.39 MB of disk space for temp files.<br>Bitmap mask created<br>Processing segments...<br>v.surf.rst complete.<br>(Mon Jul 26 22:33:00 2021) Command finished (2 sec)                             <br>(Mon Jul 26 22:33:06 2021)                                                      <br>r.mapcalc expression=pile_height_above_dsm=dsm-dsm_below_pile --overwrite       <br>(Mon Jul 26 22:33:07 2021) Command finished (1 sec)                             <br>(Mon Jul 26 22:33:14 2021)                                                      <br>g.region vector=polygon_area                                                    <br>(Mon Jul 26 22:33:15 2021) Command finished (0 sec)                             <br>(Mon Jul 26 22:33:19 2021)                                                      <br>r.volume input=pile_height_above_dsm -f                                         <br>No clump map given, using MASK<br>1:0.73:68050:93627:73.31:26.22:0.00<br>(Mon Jul 26 22:33:19 2021) Command finished (0 sec)                             <br></div><div>```</div><div>The below command in the above output gives some kind of warning, which I am not able to understand.</div><div><br></div><div>`v.surf.rst input=polygon_points zcolumn=height elevation=dsm_below_pile --overwrite`<br clear="all"><div><br></div><div>And thus at the last, the volume output is 0.00.</div><div>Any suggestions or hints, as to what is going wrong?</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Regards,</div><div>Ayush Chauhan</div></div></div></div></div>