Hello, I have a script that has ran without failure for the last 6 months but very recently started failing. The script includes lines of code that pass in a hosted feature layer from our ArcGIS Online instance into various geoprocessing tools. Here is an example:well_points = "https://services3.arc...

I am trying to work with some coordinates which are in EPSG:3346 format. Now my understanding is that if I want to do something with these coordinates in Google Maps API, the format i need is WGS84 (Is this correct?) Is there library in PHP that can make these conversions?

Business Analyst Enterprise's November 2025 update delivers new features and enhancements to various software components.

I am working in QGIS and have been wanting to do a simple Least Cost Path Analysis using the SAGA tool. I have created a Slope raster for it, as well as an Accumulated Cost Surface. I have two points which I have imported from Google Earth and transformed into two separate point layers. I have tried numerous times to run the plug in, but I keep getting the same error: ERROR: The start-point layer contains no legal point. I tried to re import the points, change the CRS, researching on the internet and watching tutorials, but nothing seem to work or to give a reason for this issue. Is there a way to fix this?

ArcGIS Enterprise 12.0 for Windows and Linux is here, bringing an array of impactful enhancements designed to improve the way you secure, scale, and manage your self-hosted GIS. Here are a few of the highlights from this release: The ability to create custom sign-in experiences for your organizatio...

Hi everyone, I developed a Python script to automate the generation of maps in ArcGIS Pro based on various parameters. In my layout, I want to be able to modify the LEGEND_Item, but all my attempts have failed so far... In fact, in my legend, I only want one layer to appear, while the others should ...

As of Oct 2025, ArcGIS Hub allows site admins to set choose how a layer is symbololized in Hub at the layer level. This means if an organization prefers to use their own symbology across multiple sites, potentially with hundreds of layers, they will need to visit each layer in each Hub content libra...

I'm trying to load a shapefile to R using the following code: state_shape

Join us for ArcGIS StoryMaps Live on December 3! Topic: Remaking classic Esri Story Map using ArcGIS StoryMaps

Hello Community Members! Today marks the 26th annual GIS Day celebration! This year's theme is "Geo-Generalist Era: Where Spatial Meets Everything" – celebrating how GIS professionals are becoming geo-generalists who apply spatial thinking across many disciplines, from public health to marketing, cl...

I am classifying a PlanetScope Imagery into 3 classes(water, non-water and mixed) using NDWI. Natural breaks (Jenks)​​ worked the best for me when I tried different data classification methods in ArcGIS Pro. Now, I need to automate this process using python. I used 'jenkspy' and it took forever to classify even a single image. When I only use sample size of 100k pixels to find the class intervals, it is faster but the accuracy is really bad. Are there other better ways to do it? Suggestions please.

Hi Community! I have problems with adding, editing, modifying and generally doing anything with layers in Map Viewer. I am unable to edit anything. I am the administrator of an online organization, and this issue affects all organization members, too. The layers and maps that I want to adjust belong...

I'd like to enable a draw function in an Openlayers 3 application for drawing a rectangle - up to a certain size (in the defined projected coordinate system, in square meters). For this ol.interaction.draw instance the draw.createBox geometryFunction would be used and the 'Circle' type, as per http://openlayers.org/en/latest/examples/draw-shapes.html?q=draw the geometryFunction is where to "modify the geometry during drawing" as per a prior thread: Undoing last point when drawing linestring in OpenLayers 3? ) but for drawing a rectangle the predefined createBox geometryFunction is used and not a custom one that might perhaps continually check the size of the rectangle. Or is this not the place to put in a check of the current size of the rectangle being drawn? Is what I wish realistic?

I want to automatically read and write settings of my QGIS-project with PyQGIS. Especially I would like to do this for the QGIS Server Service Properties (Top of QGIS-Server Tab). I found this Link but I can't solve my problem with it... I also found this thread, but I want more variables and I can't find the names of them. Someone have a full list?

Hi All I have published the feature services from utility network . These services are branch versioned. I have changed the version to branch version from default in ArcGIS Pro. I am able to update the existing data via services when doing it from ArcGIS Pro. However, while doing the add operation i...

In time for GIS Day 2025 we are releasing a new set of Guided Lessons for National Geographic MapMaker. Focusing on topics from the World Cultures and Geography curriculum and aimed at Grades 6-8, your students can explore a wide range of topics using rich interactive content. The new resources cons...

A while ago I wrote the following Java 8 code on Win 10 to create an index file for a .laz file: Process p = null; ArrayList lines = new ArrayList(); String lasbinpath = "somepath\\LAStools\\bin"; String command = "lasindex64 -i someotherpath\\myfile.laz"; try { ProcessBuilder builder = new ProcessBuilder("cmd.exe", "/c", command); Map env = builder.environment(); env.put("path", lasbinpath); builder.redirectErrorStream(true); p = builder.start(); BufferedReader br = new BufferedReader(new InputStreamReader(p.getInputStream())); String line; while(true) { line = br.readLine(); if (line==null) { break; } lines.add(line); } br.close(); } catch(IOException e) { e.printStackTrace(); } finally { if(p!=null) { p.destroy(); } } Recently I got a new Win 11 PC, installed the latest LAStools version (241017), the same Eclipse version (2020-06) and Java 8 on it and copied the old project. Yesterday I ran it once and it was working fine. Then I installed Miniconda but now it always returns an error: The "DOSKEY" command is either misspelled or could not be found. done with 'someotherpath\myfile.laz'. took 0.006 sec. It does create a .lax file (can't check if it's done correctly atm) but I've never seen this error before and definitely want to fix it. What I've already tried/checked: * Cleaned the project * Restarted Eclipse and my PC * The PATH (environmental variable) looks fine imo, only the "condabin" folder was added to it * The "doskey.exe" file exists in the "system32" folder * env.get("path") is null before I set it and it's got the correct path afterwards * When I copy this path into file explorer, it does lead to the "bin" folder, so the path is definitely correct. * The lasinfo64 command opens the lasinfo window but doesn't display anything (not sure what it should look like, as I've never used it before) and the same error is printed * Java is working fine otherwise and e.g. java -version prints the java version without any errors * I tested the same code in Eclipse 2025-09 with OpenSDK 24 and it's printing the same error * I added the path to PATH (system, not user) but it's still printing the same error The weird thing is, if I call it directly from cmd like this: cd someotherpath set path=C:\Sonstiges\RestBuero\LfU\LAZ\LAStools\bin;%path% lasindex64 -i myfile.laz ...then it works - no error. With the path added to PATH, I can omit the second line and it's still working (of course). What am I missing? Did Win 11 change how Eclipse can access the environmental variables? Is there anything I have to do/add to make it use LAStools correctly? Is this version of LAStools possibly bugged in a way?

Introduction This question is indirectly related to a previous post of mine (Calculating accumulated walking pace from starting location outwards using R?), but asks guidance about a different aspect, that's why I am creating a new post. I have put togheter some doce (below) to provide a reproducible example. Aim Calculating cumulative cost-surface expressing walking time from a given location outwards. Background I have followed the code provided by this source (https://mran.microsoft.com/snapshot/2014-11-17/web/packages/gdistance/vignettes/gdistance.pdf) to calculate an accumulated cost surface from a location outwards, using the 'gdistance' package. The cost surface has to represent walking time according to the Tobler's hiking function (https://en.wikipedia.org/wiki/Tobler%27s_hiking_function). ## EXAMPLE ## library(gdistance) # create the data to work with dtm

I have a standalone-script that should create a new QgsVectorLayer for an Oracle-connection: QgsApplication.setPrefixPath(r"C:\OSGeo4W\apps\qgis-ltr", True) # matches QGIS_PREFIX_PATH-env-var qgis_app = QgsApplication([], False) qgis_app.initQgis() layer = QgsVectorLayer('dbname=\'mydatabase\' host=\'mymachine\' port=1521 user=\'myuser\' password=\'mypassword\' srid=25832 type=MultiCurve table="myschema"."mytable" (GEOMETRIE)', 'oraLayer', 'oracle') Once I check this layer vor validitiy, I get False: layer.isValid() The exact same code works directly in the python-console in QGIS, though. Furthermore I can easily open a geopackage like this: layer = QgsVectorLayer('mygeopackage.gpkg|layername=myTable', 'gpkLayer', 'ogr') layer2.isValid() I also set the env-var QGIS_PREFIX_PATH explicitly as mentioned in this post: os.environ['QGIS_PREFIX_PATH'] = r'C:\OSGeo4W\apps\qgis-ltr'

I am exporting a vector layer (single polygon) to an ESRI Shapefile in QGIS, because Shapefile is the only supported GIS file format accepted in a data viewer I need to use. However, when exporting to ShapeFile in QGIS (version 3.40.4-Bratislava), 'only' the .cpg, .dbf, .prj, .qmd and .dwg files are created, but not the .shx file. The data viewer I am using expects .shp, .prj, .dbf AND .shx. I have already used the 'Repair Shapefile' tool (under Vector general) to generate this missing .shx file after exporting, but without success. Does anyone know how I can generate this missing .shx file, either in QGIS or with another tool?

I am trying to construct an XML file to give as input to the "Statistics and Histogram" raster function. I found that the correct formatting for the statistics was :        ...     ...     ...

If multiple people populate different attributes of the same shapefile on separate computers, is there a way to merge the attribute table without creating copies of the polygons or fields? Each file will have the same polygons and fields but everyone will be contributing to populating the fields for only some of the polygons. We want to be able to combine all of the files at the end to create one shapefile with all of the populated attributes. In testing the basic merge functions on QGIS, I found that merging created copies of the fields rather than combining them.

I would like to know how I can obtain the projection generated by a raster formed by an MDT and a raster created by me, consisting of four buildings of different heights. The problem is that it does not generate the shadow of my buildings.

Creating a line feature in ArcPro with snapping on:  My question is:  when I create a line in ArcPro, with snapping on, I have the option when digitising to snap to any other part of the existing line network (excepting the piece I am digitising) – Pic A below left – you can see the snap “test edge”...

I have a table in my Oracle database with multiple different geometry types (ie, points and lines). Now, I want to show the data in QGIS. I can create a new vector layer, however, its geometry type is arbitrarily defined as either point or line. Although all features are displayed in the attribute table, they aren't displayed in the map when their geometry type is different from the previously mentioned type chosen by QGIS. I also tried a filter expression on the table that checks for the SDO_GTYPE within the layers' filter expression: However, that prompts an error, as GET_GTYPE isn't defined (it's defined on Oracle's SDO_GEOMETRY type, but that seems inaccessible in the filter expression). Ultimately, I want to create a Python script that does this: writer_options = QgsVectorFileWriter.SaveVectorOptions() writer_options.driverName = 'GPKG' result = QgsVectorFileWriter.writeAsVectorFormatV3(layerFromOracle, outFile, QgsProject.instance().transformContext(), options=writer_options) EDIT: From https://gis.stackexchange.com/a/197800/20297 I got, that I can create a layer from any arbitrary SQL-expression using this: uri = QgsDataSourceURI() uri.setConnection('host', 'port', 'database', 'username', 'password') sql = 'SELECT * FROM myschema.mytable t WHERE t.geometrie.get_gtype() = 2' uri.setDataSource('', '(' + sql + ')', 'geometrie', '', 'id') vlayer = QgsVectorLayer(uri.uri(), 'LayerName','oracle') which I'd like to save then to my output-file. However when I check if that layer is valid, it returns False. Furthermore, when adding it to my map, I get a warning that the datasource for the layer is broken. However when opening that layers' properties, it seems, QGIS is able to query the correct data, but is unable to get the appropriate geometry-type and thus to display it properly: