[GRASS-user] v.surf.idw Not Working
Micha Silver
tsvibar at gmail.com
Sun Feb 21 05:49:39 PST 2021
Hello:
On 2/19/21 6:46 PM, Michaela Lobato wrote:
> Thanks for all your continued help. I did the steps you recommended
> and was met with a different error this time:
>
> ERROR: Database connection not defined for layer 1
>
> Traceback (most recent call last):
>
> File
> "/Users/michaelalobato/Library/GRASS/7.8/Modules/scripts/r.mblend",
> line 282, in <module>
>
> main()
>
> File
> "/Users/michaelalobato/Library/GRASS/7.8/Modules/scripts/r.mblend",
> line 240, in main
>
> gscript.run_command('v.db.dropcolumn', map=weight_points_edge, layer='1',
>
> File
> "/Applications/GRASS-7.8.app/Contents/Resources/etc/python/grass/script/core.py",
> line 441, in run_command
>
> return handle_errors(returncode, returncode, args, kwargs)
>
> File
> "/Applications/GRASS-7.8.app/Contents/Resources/etc/python/grass/script/core.py",
> line 342, in handle_errors
>
> raise CalledModuleError(module=None, code=code,
>
> grass.exceptions.CalledModuleError: Module run None v.db.dropcolumn
> map=tmp_8758116 layer=1 columns=along ended with error
>
> Process ended with non-zero return code 1. See errors in the (error)
> output.
>
> WARNING: No data base element files found
>
> WARNING: No data base element files found
>
> WARNING: No data base element files found
>
> WARNING: No data base element files found
>
> WARNING: Table <tmp_8758111> linked to vector map <tmp_8758111> does not
>
> exist
>
>
I was finally able to duplicate your error. This happens, it seems, only
when working in a longitude/latitude location. When I used two DEM
rasters in a *projected* location, the interpolation to high-resolution
in the area covered only by the lowres raster worked fine. (Also if the
lowres raster extends in two directions).
However in a Long/lat location, with units of degrees, the addon fails.
I looked thru the source code, and my guess is that there is something
in the conversion from interpolation area to vector lines that is
failing. I'll report to the maintainer.
In any case, would you want to try to transform your DEM rasters to a
projected location and rerun the module?
One caveat: be aware that interpolating a full SRTM tile (approx 110km x
110km) to a resolution of 1 meter will create a huge raster (if my
numbers are right, it would be about 50GB) and would take a very long
time. So I'd suggest to try this out first in a much smaller area.
> On Thu, Feb 18, 2021 at 7:02 AM Micha Silver <tsvibar at gmail.com
> <mailto:tsvibar at gmail.com>> wrote:
>
> No solution yet :-(, only more questions...
>
>
> It looks to me that the plugin is intended to work when the lowres
> raster extends the highres raster only in one direction. In your
> case,
> the 1 degree SRTM tile you are using totally contains the USGS 1
> meter
> region.
>
>
> To test this I would try to clip out a section of the SRTM:
>
>
> # Set the computation region to the highres raster, then extend
> only in
> one direction
>
> g.region -p rast=USGS_one_meter_x39y445_IL_12_County_ChampaignCo_2008
>
> g.region -p w=88:20:00W
>
>
> # Now clip the SRTM to this new region
>
> # But use the coarse resolution
>
> g.region -p res=0:00:01
>
> r.mapcalc "n40_w089_clip = n40_w089_1arc_v3"
>
>
> # Now try r.mblend with this clipped SRTM
>
> r.mblend high=USGS_one_meter_x39y445_IL_12_County_ChampaignCo_2008
> low=n40_w089_clip out=blend
>
>
> Please report back if this helps (or not...)
>
>
>
> On 2/17/2021 2:52 AM, Michaela Lobato wrote:
> > Here is the output for r.info <http://r.info>
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$>
> > for both maps:
> >
> > < r.info <http://r.info>
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$>
> >
> > map=USGS_one_meter_x39y445_IL_12_County_ChampaignCo_2008 at PERMANENT
> >
> >
> +----------------------------------------------------------------------------+
> >
> > | Map:USGS_one_meter_x39y445_IL_12_Date: Tue Feb 16 18:44:21 2021|
> >
> > | Mapset: PERMANENTLogin of Creator: michaelalobato|
> >
> > | Location: champaign|
> >
> > | DataBase: /Users/michaelalobato/Documents/Research/grassdata |
> >
> > | Title:USGS_one_meter_x39y445_IL_12_County_ChampaignCo_2008 |
> >
> > | Timestamp: none|
> >
> >
> |----------------------------------------------------------------------------|
> >
> > ||
> >
> > | Type of Map:raster Number of Categories: 0 |
> >
> > | Data Type:FCELL|
> >
> > | Rows: 9015 |
> >
> > | Columns:11748|
> >
> > | Total Cells:105908220|
> >
> > |Projection: Latitude-Longitude|
> >
> > |N: 40:11:40.173254NS: 40:06:11.010398N Res: 0:00:00.03651 |
> >
> > |E: 88:10:23.5766WW: 88:17:32.519389W Res: 0:00:00.036512|
> >
> > | Range of data:min = 204.0563max = 244.8284 |
> >
> > ||
> >
> > | Data Description:|
> >
> > |generated by r.proj |
> >
> > ||
> >
> > | Comments:|
> >
> > |r.proj --quiet location="temp_import_location_4329"
> mapset="PERMANEN\ |
> >
> > |T" input="USGS_one_meter_x39y445_IL_12_County_ChampaignCo_2008"
> meth\ |
> >
> > |od="nearest" memory=300 resolution=1.0142500028210239e-05 |
> >
> > ||
> >
> >
> +----------------------------------------------------------------------------+
> >
> >
> >
> > < r.info <http://r.info>
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$
> <https://urldefense.com/v3/__http://r.info__;!!DZ3fjg!s23wIXUNs2hY7ZIpBfhbowUyUKf0QNHMqg0IZeaaHL7V9Tw8y2aYLaAANBXlarcJh-k$>
> > map=n40_w089_1arc_v3 at PERMANENT
> >
> >
> +----------------------------------------------------------------------------+
> >
> > | Map:n40_w089_1arc_v3 at PERMANENT Date: Tue Feb 16 18:42:24 2021|
> >
> > | Mapset: PERMANENTLogin of Creator: michaelalobato|
> >
> > | Location: champaign|
> >
> > | DataBase: /Users/michaelalobato/Documents/Research/grassdata |
> >
> > | Title:n40_w089_1arc_v3 |
> >
> > | Timestamp: none|
> >
> >
> |----------------------------------------------------------------------------|
> >
> > ||
> >
> > | Type of Map:raster Number of Categories: 0 |
> >
> > | Data Type:CELL |
> >
> > | Rows: 3601 |
> >
> > | Columns:3601 |
> >
> > | Total Cells:12967201 |
> >
> > |Projection: Latitude-Longitude|
> >
> > |N: 41:00:00.5NS: 39:59:59.5N Res: 0:00:01 |
> >
> > |E: 87:59:59.5WW: 89:00:00.5W Res: 0:00:01 |
> >
> > | Range of data:min = 163max = 305 |
> >
> > ||
> >
> > | Data Description:|
> >
> > |generated by r.in.gdal|
> >
> > ||
> >
> > | Comments:|
> >
> > |r.in.gdal -e
> input="/Users/michaelalobato/Documents/Research/grassda\ |
> >
> > |ta/n40_w089_1arc_v3.tif" output="n40_w089_1arc_v3" memory=300
> offset\ |
> >
> > |=0 num_digits=0 |
> >
> > ||
> >
> >
> +----------------------------------------------------------------------------+
> >
> >
> >
> > Michalea
> >
> --
> Micha Silver
> Remote Sensing Lab, Sde Boker
> Ben Gurion University
> cell: +972-523-665918
>
>
>
> --
> A. Michaela Lobato
> BS Systems Engineering and Design | Autonomous Systems and Robotics |
> December 2019
> MS Systems and Entrepreneurial Engineering | Decision and Control | In
> Progress
> Email: amlobat2 at illinois.edu <mailto:amlobat2 at illinois.edu> | Cell:
> (815) 494-1111
>
--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
More information about the grass-user
mailing list