-
-
Notifications
You must be signed in to change notification settings - Fork 362
Description
Describe the bug
I have a large map (say, a Sentinel-2 scene) but only want to visualize a small portion.
For that I set the computational region to a map subset. Still InteractiveMap() sometimes crashed (as shown to @wenzeslaus in a video call).
To Reproduce
Steps to reproduce the behavior:
- Go to grass_testcase_interactiveMap.ipynb.gz
- Run the notebook
- See (sometimes) the error in the notebook. Right now I cannot find the python error message any more which I had shown.
However, while running the notebook, in the raleigh_map.add_raster("bigmap")
step r.proj
is run as follows:
ps aux | grep proj
r.proj input=bigmap@user1 output=bigmap_user1 mapset=user1 location=nc_basic_spm_grass7 dbase=/home/mneteler/grassdata resolution=1.8561528186341438
I got the impression that always the full map is fed to r.proj
. This happened with a Sentinel-2 map (cannot easily provide it here as too large).
Expected behavior
InteractiveMap()
should ideally not crash. It may be solved by changing the current usage of region management for the r.proj
call.
System description (please complete the following information):
- Fedora 35
- GRASS GIS version 8.2.dev
g.version -rge
version=8.2.0dev
date=2022
revision=66e80c0e7
build_date=2022-05-11
build_platform=x86_64-pc-linux-gnu
build_off_t_size=8
libgis_revision=b6744f2a1
libgis_date=2022-04-28T21:35:25+00:00
proj=8.2.1
gdal=3.3.3
geos=3.9.2
sqlite=3.36.0
python3 -c "import sys, wx; print(sys.version); print(wx.version())"
3.10.4 (main, Mar 25 2022, 00:00:00) [GCC 11.2.1 20220127 (Red Hat 11.2.1-9)]
4.0.7 gtk3 (phoenix) wxWidgets 3.0.5