wilcoterink / glabtop2-py Goto Github PK

When I ran the model, the following error occurred. I used to run normally.

I want to know what ‘Tau’ represents in the resulting "rgi.csv" file？

Hi Wilco,
I used the GlabTop2 code you posted on GitHub to perform calculations on glacier thickness in the southeastern region of the Tibetan Plateau. But there may be some problems with the results I got (as shown in the picture), you can see several places where the thickness of the glacier is close to 0, such as Galapai Peak and Nanga Barwa Peak on the left. The glacier boundary shapefile I entered was extracted by myself, and I used NASA's 30m dem and set the interval to 50m. Could you please tell me the reason?
I tried to send a mail to you but it didn't work, if you are convenient could you please contact me? My mail address is [email protected].
Thanks and best regards.

Hello Wilco,
Thank you for sharing this interesting model.
A few years ago, it worked fine, but now, with all the new versions of Conda and Python my colleagues and I struggle to make it work.
I tried creating a conda virtual environment with python 3.8 instead of 3.9, but it was not enough. The command conda install -c WilcoTerink GlabTop2-py fails to solve the environment, there is always some missing dependency and, as for now, I'm not very practical of conda.

It would be very useful if you update the documentation to specify what the users are supposed to do to make it work, considering that they will download the latest version of anaconda.

THANK YOU!!
Andrea

Dear Wilco Terink,
You have made a very useful software, thank you very much!
I have a questions. 1) Can glabtop2 assign a value of 0m only for the edges of glacier polygons (QGIS) without doing this with ice divides? Otherwise the depth of the ice is greatly underestimated (photo), this greatly affects the calculation of the total volume of ice.
2) The vertical spatial resolution of digital terrain models is inferior to instrumental measurements in the field, such more accurate data is available for some glaciers. Is it possible, for example, to manually enter these values in the attribute tables of shapefiles so that glabtop2 uses them further, and where they do not exist, them uses DEM?
Sincerely yours,
Semyon


[email protected]

Interesting points. About point 1) You are correct, but from my point of view those glaciers should not be individual glacier polygons but one larger polygon. Why would they need to have a different Glacier ID if they are part of the same ice mass?
Point 2) I have to think about that, but yes, that should not be too difficult. On the other hand, I do not think a single measurement can be representative for the entire grid cell resolution? Unless you do multiple measurements in that grid cell area.

Dear developers,
as for a previous support request, I finally made the GlabTop2 run on my study glacier. Many thanks for your help!!

I would like to tell you one thing that I noticed: in some areas near the glacier boundary, the GlabTop2 estimates a thickness larger than 100 m. My DEM resolution is 20 m. (result did not change much with 5m resolution)

While it may be physically possible, in case of steep slopes, it actually does not correspond to the test site situation: in the areas near the glacier boundary with estimated thickness >100m, the nearby rock slope is not particularly steep.

Am I missing something? e.g. 1) it is actually common for glaciers to have super steep slopes near the boundary, even if nearby slopes are not steep - this is my preferred option. Or 2) am I doing something wrong with the data? (i attach them here). Or 3) is there some parameter in the model that can be tuned to better adjust the interpolation to 0-thickness at the boundary?

Thank you very very much for your insight! I hope to be able to cite GlabTop2 soon :)

have a nice Christmas.

Andrea

test_Rutor.zip

Dear WilcoTerink,
I am getting the above error when I try to run/install GlabTop2-py via Anaconda, How do I deal with it?
I have installed the PC Raster then try to install the GlabTop-2-py. I am new to programming.
My email address is [[email protected]], look forward to your reply.
Yours sincerely,
Shashi

Traceback (most recent call last):
File "D:\ProgramData\Anaconda3\lib\runpy.py", line 193, in _run_module_as_main
"main", mod_spec)
File "D:\ProgramData\Anaconda3\lib\runpy.py", line 85, in _run_code
exec(code, run_globals)
File "D:\ProgramData\Anaconda3\lib\site-packages\GlabTop2\main.py", line 313, in
main(sys.argv[1:])
File "D:\ProgramData\Anaconda3\lib\site-packages\GlabTop2\main.py", line 304, in main
glabtop.run()
File "D:\ProgramData\Anaconda3\lib\site-packages\GlabTop2\main.py", line 253, in run
[glacmask, mCells, innerCells] = self.identifyCells(gID)
File "D:\ProgramData\Anaconda3\lib\site-packages\GlabTop2\main.py", line 102, in identifyCells
if glacMask[r+1,c+1] == self.mv:
IndexError: index 189 is out of bounds for axis 1 with size 189

Dear WilcoTerink,
I am getting above error when I try to run GlabTop2-py via Anaconda, How do I deal with it?
Or, can you provide some input files that can run glatop2-py successful?
My email address is ([email protected]), I will be deeply appreciated that you have done for me and look forward to your reply.
Yours sincerely,
Hongyuan