wilcoterink / glabtop2-py Goto Github PK
View Code? Open in Web Editor NEWGlabTop2 (Glacier bed Topography) model
License: GNU General Public License v3.0
GlabTop2 (Glacier bed Topography) model
License: GNU General Public License v3.0
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
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
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.
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
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.
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
my issue is that, i started processing and waited for 6-7 hrs but still got no output. Also, i was not able to figure out whether the processing is going on or not or how much processing is completed or left. Please help me in this matter and provide some instructions so that i can complete the processing of glacier with few time left.
please reply me , i have provided my email id below.
Name :- Akash Basotra
Email:- [email protected]
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.