This is a backgrounder to things spatial, like point in polygon analysis etcetera. Code links are at the end. ——————-For point in polygon, spatial joins and the like, the normally accepted definitions consider points on the boundary of a polygon are "outside". There are situations where you may want to include them.…

Hi everyone, I’m working on a web-based relationship duration calculator at https://kalkulatorjadian.net/ that calculates how long two people have been together from their relationship start date. The calculator can show the elapsed duration in days, hours, minutes, and seconds. I’m interested in taking this concept a step…

I added an online map to my copy of OSMAnd, using the following address: https://data.geopf.fr/wmts?layer=GEOGRAPHICALGRIDSYSTEMS.COASTALMAPS&style=normal&tilematrixset=PM&Service=WMTS&Request=GetTile&Version=1.0.0&FORMAT=image/png&TileMatrix={0}&TileCol={1}&TileRow={2} This works well (it corresponds to nautical charts of the French coast), the map appears in my map sources and I can use it as a basemap. As far as I understand, this points to a WMTS server. QGIS can, of course, use WMTS. However, I cannot configure my QGIS connections to reach this layer. The obvious things (such as pointing to https://data.geopf.fr/wmts?GetCapacities ) do not seem to work. What would be the correct QGIS settings ?

DCAT-US 3.0 (beta) Now Available in ArcGIS Hub ArcGIS Hub now supports DCAT-US 3.0 (beta)! ArcGIS Hub empowers organizations to share and federate authoritative content through standards-based catalog feeds, like DCAT, without duplicating data. This makes trusted data easier to discover, access, and reuse across open data…

In Experience Builder Developer Edition (1.21), if a widget's files fail to load, the widget shows "Cannot read properties of null (reading 'default')". Closing and reopening the widget shows the same error. The only fix is refreshing the page, which loses any unsaved work like drawings. My hypothesis is that this happens…

I am new to using survey123. I inherited an old observation submission survey that needed a lot of major updates, so I ended up recreating it as a new survey with the features and fields I wanted. I am now ready to roll that new survey out. I have never "unpublished" an existing survey from survey123. Is there a way to…

Building Spatial Thinking with National Geographic MapMaker How do you help students understand not just where places are, but how places relate to one another? For seventh-grade global studies teacher Jess Winkelaar, National Geographic MapMaker has become a way for students to explore those relationships and build…

I have a 32-bit GeoTIFF file I'm trying to convert to map tiles using gdal2tiles.py. What the GeoTIFF is supposed to look like when loaded straight into QGIS (no scaling, styling, etc): And here is the output of gdalinfo input.tif: bash-3.2$ gdalinfo input.tif Driver: GTiff/GeoTIFF Files: input.tif input.tif.aux.xml Size is 294, 192 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (-83.400000000000006,35.200000000000003) Pixel Size = (0.016666666666667,-0.016666666666667) Metadata: AREA_OR_POINT=Area Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=BAND Corner Coordinates: Upper Left ( -83.4000000, 35.2000000) ( 83d24' 0.00"W, 35d12' 0.00"N) Lower Left ( -83.4000000, 32.0000000) ( 83d24' 0.00"W, 32d 0' 0.00"N) Upper Right ( -78.5000000, 35.2000000) ( 78d30' 0.00"W, 35d12' 0.00"N) Lower Right ( -78.5000000, 32.0000000) ( 78d30' 0.00"W, 32d 0' 0.00"N) Center ( -80.9500000, 33.6000000) ( 80d57' 0.00"W, 33d36' 0.00"N) Band 1 Block=294x3 Type=Float64, ColorInterp=Gray Min=0.005 Max=0.759 Minimum=0.005, Maximum=0.759, Mean=0.057, StdDev=0.076 NoData Value=-1.69999999999999994e+308 Metadata: STATISTICS_MAXIMUM=0.75856912136078 STATISTICS_MEAN=0.057134358222004 STATISTICS_MINIMUM=0.0052262963727117 STATISTICS_STDDEV=0.075826669687323 What a sample tile looks like when converted using gdal2tiles.py input.tif: OK, so clearly there's something going on with floating points and possibly number signing. So I try to convert it: gdal_translate input.tif output.tif -ot UInt16 ... And then convert to tiles: Result: Only values over 0.5 carry over as black. Everything else is white. I feel like there is something obvious that I'm missing. Does anyone have any insights into what's going on here? Things I've done: -Try all the different convert options: Byte: Monochrome conversion Int16: All black UInt16: Monochrome conversion UInt32: Monochrome conversion Int32: All black Float32: TV snow noise Float64: TV snow noise CInt16: All black CInt32: All black CFloat32: TV snow noise CFloat64: TV snow noise

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo = wms.getfeatureinfo( layers=['0'], styles=['default'], srs='EPSG:4326', bbox=(-68.68876,-23.09997,-67.82525,-21.12957), # manually set feature coordinates size=(500,500), format='image/jpeg', query_layers=['0'], info_format='application/geojson', xy=(250,250) ) out = open(wms['0'].title.replace(" ", "_")+".json", "wb") out.write(fInfo.read()) out.close() So my problem is exactly that: I need to know beforehand where are all the features and specify a bounding box inside them. Is there a way to list all features or objects from a WMS and download their metadata?

Is it a direct projection on 20 faces or on 32 faces? Can I use PostGIS (with ISEA) to calculate and obtain cell-identifiers? For hexagons (like Uber H3) or for diamonds (like DGGRID's ISEA4D). The online literature say nothing: * the Wikipedia short page about ISEA. No clues about algorithm. Citation into the article: Begin with a Snyder Equal Area projection to a regular icosahedron (...) inscribed in a sphere (...) Then project the hexagon back onto the sphere using the Inverse Snyder Icosahedral equal area projection. * the PROJ page for ISEA, used by PostGIS and many other, would be the best source, but it doesn't say much, on the contrary, it creates some confusion by saying "... map projections on polyhedral globes for the dodecahedron and truncated icosahedron". And not shows how to use each polyhedre. * This paper explains how to transform regular icosahedron faces into projection planes and, them, into a Truncated icosahedron of 32 faces... But I am not a good reader of mathematical formalism, nor of Chinese. --- NOTES Clues by illustrations: * ... * Step n+1: Lambert azimuthal equal-area projection over each 20 faces of the icosahedron. * Step n+2: "project back onto the sphere using the Inverse Snyder" each 5 faces to truncate, generating 12 new pentagonal faces. So, the algorithm has a reprojection step? * ...

I am going through the gdal Syntax documentation and am trying to break down the components of a command sent to the CLI. I understand that there is the gdal program name, commonly followed by a top-level command, then often followed by a subcommand. Synopsis for each command name (such as for hillshade) break the rest of command syntax into positional arguments, common options, options, and advanced options, etc. The Syntax documentation page mentions arguments (referenced with either long or short names), some of which are flags if they can be called without a value. But there is no mention of options. So is it safe to say then that everything other than the command names are arguments? In other words, that options are types of arguments that are optional?

A comprehensive look at every location service in ArcGIS Location Platform, from basemaps to elevation, with capabilities and real use cases.

Learn how to use imagery embeddings in ArcGIS Pro to identify features in imagery and transfer the results to classify overlapping lidar point clouds without training a custom AI mod

Good afternoon all, I'm hoping to get some people's opinions on single ExB application with lots of pages/views vs creating multiple ExB applications in AGOL. I am building an application for data editing and visualisation for my organisation. There are several stages to this, each with different web maps, widgets and…

I previously mentioned this issue in a comment on another discussion , but thought I had better raise it in this more appropriate forum instead. We are seeing failing requests to hosted feature services where a large number of values are passed in the where clause. The error is 'SQL query is nested too deeply. Rewrite the…

ArcGIS Hub boosts data discoverability with standards-based catalog feeds for open data sharing and federation.

I have published a free, MIT-licensed Python toolbox for ArcGIS Pro. One coordinate in, a full river-basin analysis out. Fourteen tools — one runs everything, the other thirteen do a single piece each, so they work in ModelBuilder. It delineates the basin, downloads elevation, land cover, soil, rainfall, surface water and…

Anyone experience an issue with pop ups on a feature layer getting mixed up and showing info on the wrong features? It happens when the layer is set to automatically refresh. In my example i have a wildfire layer set to refresh after 5 minutes. After that refresh occurs i can click a wildfire from TX and it shows info for…

This idea is a minor change to nothing more than the messaging generated by the tool as it runs. It would be useful and informative that when the tool runs it states in the messaging what was the initial extent and what is the newly calculated extent. This of cause could be the same set of coordinates if there was no…

I have a GeoTIFF file which when imported using rasterio get flipped on the y-axis, and now I am trying to rotate this GeoTIFF file using NumPy rotate, it gives an error (Axes=(0, 1) out of range for array of ndim=0) rf = rasterio.open(r'C:\Users\NESAC\tp.geotif') fig, ax= plt.subplots(1,1) show(rf, ax = ax, title = 'temperature') rot = np.rot90(rf, k=3)

I have a shapefile that contains a column of orientation. I use that column to assign the orientation of the marker. My problem is that when I transfer the file into a leaflet webmap, I get an error regarding the orientation: "line 217, in getSymbolAsStyle if symbol.dataDefinedAngle().useExpression(): AttributeError: 'QgsProperty' object has no attribute 'useExpression'" Is there a workaround I could do?

With PyQGIS how can I convert absolute screen coordinates (say, 123,456) into map coordinates of the iface.mapCanvas? It would be no problem if I had screen coordinates relative to the iface.mapCanvas.viewport(), since iface.mapCanvas.viewport().geometry() gives me sth like QRect(2,2,2198, 1090) and iface.mapCanvas().extent() gives me a QgsRectangle with xmin,xmax,ymin,ymax, so (leaving projection issues out of the scope) the conversion would be simple. But since the absolute position of the mapCanvas may vary with toolbar positions or layer windows width, the question reduces to 'how to derive the absolute screen position of the map canvas?'.

I would like to move the cursor to a given x,y coordinates of the Map Canvas. Is there an easy way to do this using PyQGIS?

I have a question regarding the aggregation of a binary raster to a dataset with neighbourhood census tract polygons to research how much of the population in neighbourhoods complies with the 3-30-300 rule which means a person has to ive in a neighbourhood with 3 visible trees, 30% tree canopy cover / water surface, live 300 metres from a high quality green space. for each part of the rule i will create a 10 x 10 raster with binary scores 0 being non-compliant and 1 being compliant. Then I want to calculate the amount of people that live within those 10 x 10 cells and sum them up for the neighbourhood polygons. What is the best way to accomplish this. I was thinking either aggregating the mean of the binary raster tothe neighbourhood polygons directly using zonal statistics, this will however result in less accurate estimations of how many people really live in accordance with the rule and then would do (total population of neighbourhood * mean compliance) Or i could aggregate the binary raster to a 100 x 100 census tract raster database and then summing all the 100 x 100 cells up within the neighbourhood polygon. Is there any other way I could calculate how many people within the neighbourhood polygon comply with the rules? Keep in mind that the R script runs for 43 cities in the netherlands, so it would be nice to keep a fast computation time. thanks in advance

I am curious what the 'fixed' status means on this bug that is affecting us. BUG-000186104 for ArcGIS Velocity for ArcGIS Enterprise I am super confused because it says it was fixed in Velocity for Enterprise 12.1.1 but I don't get what that means. I don't see any patches available on the Server or a new version of…

How can I rename a field in PyQGIS? I tried: for field in layer.pendingFields(): if field.name() == 'old': layer.startEditing() oldname = field.name() field.setName('new') newname = field.name() print(u'Old name: {}, New name: {}'.format(oldname,newname)) layer.commitChanges() This code prints 'Old name: old, New name: new', but in the attribute table the field name is still 'old'.

Learn about ArcGIS Earth desktop 3.0 and mobile 2.8, introducing Local Scenes, Building Scene Layers, Gaussian splats, and Projects.

Good afternoon folks! We have an old ArcGIS Enterprise 10.9.1 deployment that we're intending to upgrade to 11.5 soon. The ArcGIS Servers are federated with the Portal and anonymous access is turned off. When investigating user access while planning the upgrade, we noticed that our ArcGIS Servers allow us to access…

Here's how you can give contours a better sense of uppy-downiness.

When using the Varying Distances option in the Divide Tool for lines, add the capability to insert a row in between existing rows instead of just adding at the end. Similar to the insert before in the Traverse tool. It would also be nice to be able to re-order the rows as well.

The 2027 Esri Partner Conference is our largest annual event for our partners and distributors. A reminder... this year's event is going to be held March 7–8, 2027. This is a great opportunity to put your organization in front of the right audiences, build new connections, and strengthen existing relationships. At EPC,…

I've tried to debug this without success, so am wondering if anyone else has run into this issue. Here's what I'm trying to do: * Use the OOTB UPDM model - with a few additions, but no changes to asset type or CP fields * Use the Cathodic_Protection_Traceability_Distribution_Pipe attribute rule to update the cptraceability…

Hi all! My name is Haley, and I work as a GIS Specialist at an engineering design firm in Tennessee, primarily working with water utilities and local governments. I do a lot of utility mapping, data validation and quality assurance, and asset management! I previously worked as a water engineer in my firm, but upon finding…

Hi everyone, I’m using the Utility Network Trace widget in ArcGIS Experience Builder on ArcGIS Enterprise 12.1. When I run a downstream trace, the results highlight on the map, but the corresponding records do not appear as selected in the Table widget. What I’ve checked: * Select records is checked in the Trace widget’s…

I could able to remove the recent modified vertex of a feature using removePoint method of modify interaction, as well as using deleteCondition event of modify interaction. However, that is not what I am trying to achieve. The use case for this this problem is as below. I have modify interaction added on map. Now when user tries to modify a feature over map, I am highlighting the vertices, like this so that user can directly go to these vertices and then do right click which will result in removal of vertex. For this, I had to implement the contextmenu to handle the right click over map. But, now the challenge is to remove the received coordinate from the modifying geometry. I tried to implement it as follows, with simple .splice() method on coordinate array of target geometry if(this.modifyInteraction.features_.getLength() && this.modifyInteraction.features_.getArray()[0]){ //Feature being modified let modifyft = this.modifyInteraction.features_.getArray()[0]; let coord = modifyft.getGeometry().getCoordinates()[0]; let foundIndex; for(var k = 0; k < coord.length; k++){ if(coord[k][0] == rightClickedVertex.getGeometry().getCoordinates()[0] && coord[k][1] == rightClickedVertex.getGeometry().getCoordinates()[1]){ foundIndex = k; break; } } if(foundIndex && foundIndex > -1){ modifyft.getGeometry().getCoordinates()[0].splice(foundIndex, 1); //nothing changes } } is there any other way to remove the vertex of geometry being edited?

I get this error when I run this: vector_layer = "mask.shp" raster_layer = "raster.tif" target_layer = "mask.tif" # open the raster layer and get its relevant properties raster_ds = gdal.Open(raster_layer, gdal.GA_ReadOnly) xSize = raster_ds.RasterXSize ySize = raster_ds.RasterYSize geotransform = raster_ds.GetGeoTransform() projection = raster_ds.GetProjection() # create the target layer (1 band) driver = gdal.GetDriverByName('GTiff') target_ds = driver.Create(target_layer, xSize, ySize, bands=1, eType=gdal.GDT_Byte) target_ds.SetGeoTransform(geotransform) target_ds.SetProjection(projection) # rasterize the vector layer into the target one ds = gdal.RasterizeLayer(target_ds, [1], ogr.Open(vector_layer).GetLayer(), burn_values=[1]) target_ds = None It happens when I call target_ds = None after RasterizeLayer. If I don't run RasterizeLayer it does not fail. So what I am trying to do is create a mask of the polygons in mask.shp using the projection and geotransform information from raster.tif. Note I am using gdal version 1.10.1.

Learn how ArcGIS Urban metrics can model parking requirements, assess development potential and support zoning policy decisions.

Image analysis in ArcGIS turns imagery into answers with raster math, SAR, hyperspectral analysis, and GeoAI across Pro, Enterprise, and Online.

Learn about an updated user experience and new capabilities in the Crime Analysis solution for ArcGIS Pro 3.7.

I have used the Minimalist Plugin Skeleton to create some very basic plugins which now need updating to be ready for versions of QGIS using Qt6. I have made some changes so I now use the version independent qgis.PyQt.QtWidgets rather than PyQt5 classes, and the plugin runs successfully when re-loaded using the Plugin Reloader plugin. My issue is that QGIS thinks: This plugin is incompatible with this version of QGIS and I am not sure why. The plugin metadata only specifies the minimum version required: [general] name=Run a processing model description=Does some magic using a processing model. version=0.1 qgisMinimumVersion=3.0 author= me email=me@home.com icon=logo.png I have tried the minimalist plugin without any of my code and that isn't compatible, even after I changed to import the widgets from independent qgis.PyQt.QtWidgets.

I would like to understand whether the Enterprise Utility Network Associate (EUNA) 2026 exam includes the Telecom domain as part of its scope. Based on the exam syllabus and study materials that I found in the Esri Community, the topics appear to focus on Utility Network concepts and workflows but do not explicitly mention…

After publishing new services to ArcGIS Server, the corresponding items do not appear in the Portal's Content or search results. The services themselves are available and return correctly through the ArcGIS Server Services Directory.

I have a polygon and a line feature class inside of a file geodatabase. Both the feature classes, and the map reference Washington State Plane North NAD83. After doing a spatial join moving attributes from the polygons to the lines I notice the line shape_length attribute had taken the values of the shape_length from the…

AEC (Architecture, Engineering and Construction) Documents

I have been trying to optimize a series of snow routes and I keep getting a gap where the route should continue. I've told it no u-turns at the order location and sometimes it closes the gap, but it still makes gaps in other places when I re-run it. The plow has to go over all the street segments that are a snow route. How…

It would be useful to have an optional confirmation prompt when stopping an active line capture in QuickCapture. Our biosecurity team uses QuickCapture to record GPS tracks while inspecting properties, alongside buttons for rapidly capturing pest locations. Currently, a single tap on the active tracking button stops the…