ChHarding
commented
1 year ago I did run into this just before Christmas. Turns out b/c the Dockerfile uses conda, everything runs in a base env. I had assumed something would activate base but it wasn’t - and still everything runs fine. Except GDAL, the complicated, arcane beast that it is, which needs to go through some sort of init to know where the projection stuff is. (Finding that out was quite a fishing expedition :) For this to happen first the shell script /usr/local/etc/profile.d/conda.sh has to be run, which set up where conda's commands like activate are, which are ofc also shell scripts! Only after does conda activate base work.
Need to activate base conda env so that GDAL projection stuff works
source /usr/local/etc/profile.d/conda.sh # inits conda commands for bash conda activate base
You could run this manually but I’ve also added it to the top of the run_touchterrain.sh script:
[cid:26DA2B7A-C40A-43B0-AB96-F0F924498C16]
Let me know if that solved it.
Cheers
Chris
On Jan 15, 2023, at 19:58, ultrasouls @.**@.>> wrote:
I'm getting an error when I use a .gpx file for args["importedGPX"]
When I get the error, if I go back and comment out that line and use the args["importedGPX"] = None the error is resolved.
Also, my GPX file is working fine as it shows up on the geemap when I'm going to select the area to print. The error only comes up when running the Generating the STL model file code.
That’s b/c gee map can natively read GPX files
Warning 1: PROJ: proj_create_from_database: Open of /usr/local/share/proj failed Warning 1: The definition of projected CRS EPSG:32614 got from GeoTIFF keys is not the same as the one from the EPSG registry, which may cause issues during reprojection operations. Set GTIFF_SRS_SOURCE configuration option to EPSG to use official parameters (overriding the ones from GeoTIFF keys), or to GEOKEYS to use custom values from GeoTIFF keys and drop the EPSG code. ERROR 1: PROJ: proj_create_from_database: Open of /usr/local/share/proj failed
AssertionError Traceback (most recent call last) Cell In[41], line 2 1 # create zipfile under tmp with the STL file of the terrain model ----> 2 totalsize, full_zip_file_name = TouchTerrain.get_zipped_tiles(**args) 3 print("Created zip file", full_zip_file_name, "%.2f" % totalsize, "Mb")
File /TouchTerrain/standalone/touchterrain/common/TouchTerrainEarthEngine.py:902, in get_zipped_tiles(DEM_name, trlat, trlon, bllat, bllon, polygon, polyURL, poly_file, importedDEM, printres, ntilesx, ntilesy, tilewidth, basethick, zscale, fileformat, tile_centered, CPU_cores_to_use, max_cells_for_memory_only, temp_folder, zip_file_name, no_bottom, bottom_image, ignore_leq, lower_leq, unprojected, only, original_query_string, no_normals, projection, use_geo_coords, importedGPX, gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, map_img_filename, smooth_borders, offset_masks_lower, fill_holes, min_elev, **otherargs) 900 if importedGPX != None and importedGPX != []: 901 from touchterrain.common.TouchTerrainGPX import addGPXToModel --> 902 addGPXToModel(pr, npim, dem, importedGPX, 903 gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, 904 trlat, trlon, bllat, bllon) 906 # clip values? 907 if ignore_leq != None:
File /TouchTerrain/standalone/touchterrain/common/TouchTerrainGPX.py:233, in addGPXToModel(pr, npim, dem, importedGPX, gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, trlat, trlon, bllat, bllon) 231 source = osr.SpatialReference() 232 s_res = source.ImportFromEPSG(4326) # This is WGS84, return of 0 => OK --> 233 if s_res != 0: assert False, "addGPXToMode(): source.ImportFromEPSG(4326) returned error" + str(s_res) 235 for gpxFile in importedGPX: 236 pr(f"process gpx file: {gpxFile}")
AssertionError: addGPXToMode(): source.ImportFromEPSG(4326) returned error6
— Reply to this email directly, view it on GitHubhttps://github.com/ChHarding/TouchTerrain_jupyter_docker/issues/8, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AEYDF5PVIYMZMVKOCKEL34LWSSTNBANCNFSM6AAAAAAT4GPQRU. You are receiving this because you are subscribed to this thread.Message ID: @.***>
Chris Harding Associate Professor Department of Geological & Atmospheric Sciences Touchterrain.geol.iastate.eduhttp://Touchterrain.geol.iastate.edu
ultrasouls
I'm getting an error when I use a .gpx file for args["importedGPX"]
When I get the error, if I go back and comment out that line and use the args["importedGPX"] = None
the error is resolved.
Also, my GPX file is working fine as it shows up on the geemap when I'm going to select the area to print. The error only comes up when running the Generating the STL model file code.
Warning 1: PROJ: proj_create_from_database: Open of /usr/local/share/proj failed Warning 1: The definition of projected CRS EPSG:32614 got from GeoTIFF keys is not the same as the one from the EPSG registry, which may cause issues during reprojection operations. Set GTIFF_SRS_SOURCE configuration option to EPSG to use official parameters (overriding the ones from GeoTIFF keys), or to GEOKEYS to use custom values from GeoTIFF keys and drop the EPSG code. ERROR 1: PROJ: proj_create_from_database: Open of /usr/local/share/proj failed
AssertionError Traceback (most recent call last) Cell In[41], line 2 1 # create zipfile under tmp with the STL file of the terrain model ----> 2 totalsize, full_zip_file_name = TouchTerrain.get_zipped_tiles(**args) 3 print("Created zip file", full_zip_file_name, "%.2f" % totalsize, "Mb")
File /TouchTerrain/standalone/touchterrain/common/TouchTerrainEarthEngine.py:902, in get_zipped_tiles(DEM_name, trlat, trlon, bllat, bllon, polygon, polyURL, poly_file, importedDEM, printres, ntilesx, ntilesy, tilewidth, basethick, zscale, fileformat, tile_centered, CPU_cores_to_use, max_cells_for_memory_only, temp_folder, zip_file_name, no_bottom, bottom_image, ignore_leq, lower_leq, unprojected, only, original_query_string, no_normals, projection, use_geo_coords, importedGPX, gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, map_img_filename, smooth_borders, offset_masks_lower, fill_holes, min_elev, **otherargs) 900 if importedGPX != None and importedGPX != []: 901 from touchterrain.common.TouchTerrainGPX import addGPXToModel
--> 902 addGPXToModel(pr, npim, dem, importedGPX, 903 gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, 904 trlat, trlon, bllat, bllon) 906 # clip values? 907 if ignore_leq != None:
File /TouchTerrain/standalone/touchterrain/common/TouchTerrainGPX.py:233, in addGPXToModel(pr, npim, dem, importedGPX, gpxPathHeight, gpxPixelsBetweenPoints, gpxPathThickness, trlat, trlon, bllat, bllon) 231 source = osr.SpatialReference() 232 s_res = source.ImportFromEPSG(4326) # This is WGS84, return of 0 => OK --> 233 if s_res != 0: assert False, "addGPXToMode(): source.ImportFromEPSG(4326) returned error" + str(s_res) 235 for gpxFile in importedGPX: 236 pr(f"process gpx file: {gpxFile}")
AssertionError: addGPXToMode(): source.ImportFromEPSG(4326) returned error6