texts: ArcGIS Javascript Legend widget does not show full color ramp

texts:
Has anyone used the legend widget with raster map services? I have created a GP service that returns raster maps, but the legend is incomplete.

The ArcGIS layer properties are shown here, and below how the legend looks in ArcMap TOC (left), REST services directory (middle), and result web map (right).

It looks like by default, there is only a high and low label. I added the Medium label manually before publishing, but the color ramp seems to be recalculated when you run the GP tool through the Javascript API.

Any ideas on how I can get the full color ramp to show up?

Update:
This has been logged as two bugs, see http://support.esri.com for more details:

BUG-000092784 : Running a Geoprocessing service at Rest endpoint and publishing the result as a map service causes the Legend color ramp values to change.

BUG-000092783 : Running a Geoprocessing service using the javascript API and adding the result map and legend to the application causes the medium value of the legend to disappear.


texts: I clip a raster with a vector mask in qgis, the result is correct but its a little bit of, why is this?

texts:
On one side the result is "bleeding" over the vector on the other side its more than one raster pixel from the vector on the inside.

How to correct this?

http://i.imgur.com/d0v6mqt.png

http://i.imgur.com/UKoEepJ.png


texts: Create an animated proportional symbol map

texts:
I have a number of locations (point data)
Each location have an attribute of different values over time (different columns)
I want to make an animated map with proportional symbols showing how these values of the same location change over (i.e. the dots become bigger over time). The below examples are similar to what I want to achieve.

http://bl.ocks.org/lucguillemot/2657e547d5556cc0acbb

https://observablehq.com/@jashkenas/italy-coronavirus-daily-cases-map-covid-19

Unfortunately, I am not good at coding. Can I make an animated proportional symbol map via the existing plugin Time Manager in QGIS?

I can try to reconstruct the data of the same location to different rows instead of the column if that is needed.


texts: Tool for generating Feature Class including field names/type

texts:
I am trying to build a tool that will generate a feature class and then populate the fields with the correct column name and format based off a CSV file.

Has anyone had any luck doing this before?

I'm looking for model or code examples if anyone happens to have anything similar.

Has anyone seen any tutorial documents online that do this type of request?

I am using Arc 10.2. ESRI:Arcpy


texts: gxp_featuremanager nextpage throws Cannot do natural order without a primary key, please add it or specify a manual sort over existing attributes?

texts:
I am using GeoServer and Boundless to build web mapping apps.

I have a feature grid that displays the selected layers' geometry/features.
On load it shows the first 100 features, maxFeatures Number default is 100 - as per the documentation.
So no issue there, the first 100 fetures are loaded and displayed on the map. However, when I click the button to show/retrieve the next 100 features I receive this response:

Status: 200

<ows:ExceptionReport xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:ows="http://www.opengis.net/ows" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.0.0" xsi:schemaLocation="http://www.opengis.net/ows http://localhost:9080/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd">
<ows:Exception exceptionCode="NoApplicableCode">
<ows:ExceptionText>java.lang.RuntimeException: java.io.IOException
java.io.IOExceptionCannot do natural order without a primary key, please add it or specify a manual sort over existing attributes</ows:ExceptionText>
</ows:Exception>
</ows:ExceptionReport>

It does not load any more features, but returns status 200.

The DB table I used does have valid geometry and  a primary key. Can anybody please explain why this  is happening or how I can fix this?


texts: Can the Reporting Tool in ArcMap connect to Auxiliary data, like Excel tables?

texts: Is it possible in ArcMap 10, using the reporting tool, to connect a report to non-spatial data, like from Excel or Access tables?


texts: ArcMap Time Slider-> features disappear when I enable time function!

texts:
When I go to the properties and enable the time function, seen in the screen clip 1, and then click OK or apply, my polygons (and their attributes)disappear and are not visable in the time slider tool, see screen clip 2. I don't understand why?

I can get my data back by disabling the time function. I am using ArcMap 10.3.

screen clip 1

screen clip 2


texts: Automatically moving overlapping features using ArcGIS Desktop?

texts:
I am using ArcGIS 10.4.1. I make a map. I am have two shapefiles: one polyline and one polygon. The polyline is a road and polygon is a building.

The shapefiles overlap. I want to automatically resolve the overlapping features. I want to make a "After" picture.

If polygon is placed at the center of the line, it doesn't matter which direction it moves. You only need to move to either side. I want to moving distance is likely to be 15m-20m away from the road. The building should be located only near the road as shown.

I wonder if ArcGIS can automatically shift buildings by road.

I don't have programming skills. However, the GIS program is well used.I hope to solve this problem with ArcGIS.

How to automatically move overlapping features?

The above figure is one line and the figure below is three lines. The three lines mean the center line of the road and the two outside mean the width of the road


texts: Use ArcGIS for Server geocoding service results with PyQGIS

texts:
An ArcGIS for Server geocoding service can be used via http get (singleline geocoding for 'Jugendamt' [youth office]) http://secure.erlangen.de/arcgiser/rest/services/Comp_fromModel/GeocodeServer/findAddressCandidates?Address=&Postal=&City=&Region=&Country=&SingleLine=Jugendamt&category=&outFields=&maxLocations=&outSR=&searchExtent=&location=&distance=&magicKey=&f=html

and returns some results (below, Address Candidates):

Question: Is it possible to extract the resulting coordinates for further use, e.g. in PyQGIS?


texts: Is there an easy way to get the bounding box for an entire table in PostGIS?


texts:
I'm very beginner in QGIS and I want to erase Features by using PyQGIS. I have written this code. So far, it give me the geometry of each feature polygon layer, but I don't know how to erase it through coding?
layer = iface.activeLayer() #get active layer
for f in layer.getFeatures():
  geom = f.geometry()
  print geom.asPolygon() #get polygon

As picture suggested, a big polygon contains small polygon. I want to erase all small polygons in that same layer so if i move small polygon after erase hallow will be shown in big polygon.


texts: Add parsed GeoJSON points to OpenLayers and make them Clickable

texts:
I have points layer, displayed in OpenLayers map. I want to make them clickable, so when I click on every point I get the property I want (ex: coordinates).

I'm newbie in OpenLayers and don't know where to start from.

Here is my code :

map.on("pointermove", function (evt) {
    var hit = this.forEachFeatureAtPixel(evt.pixel, function(feature, layer) {
        return true;
    });
    if (hit) {
        this.getTargetElement().style.cursor = 'pointer';

        //click event
        map.on('singleclick', function (event) {
        var outpout = event.coordinate;
        console.log(outpout);

});

    } else {
        this.getTargetElement().style.cursor = '';
    }
});

when I click on  a point I get the coordinates of all the points


texts: Does overpassAPI setup supports importing from osm.pbf file? The docs are only mentioning the use of planet.osm file.


texts: How do i convert STDS file to ASCII DEM?

texts:
I've downloaded several 7.5 minute DEM files from USGS and they have file names like:

1661629.DEM.SDTS.TAR

On MAC these decompress into a folder with serval files having the extension '.DDF'. Is there a shareware program that allows me to convert them to a basic ASCII text DEM file under the Mac OS X operating system?


texts:
I already have basemap data. For that Bing Maps or OpenStreetMap are fine.

However, for demonstration purposes I need data for operational layers that I can throw on top of this.

It needs to be free, within Australia, and open for use.

For example:

Flood Zone Data
Historic bush fire data
Historic weather data
Census data
Crime data
Disease data

Are there other sources of Australian spatial data?


texts:
Could someone please help me with my issue? I try to use Mapproxy with OpenStreetMap server of tiles in order to get output WMS for my qgis client. Something like this:
...

sources:
  arcgis_tms:
    type: tile
    url:  http://c.tile.openstreetmap.org/%(tms_path)s.png
grids:
  arcgis_online:
    base: GLOBAL_MERCATOR
    num_levels: 18
    origin: nw

But I see something strange with my tiles(North and West are confused):http://img-fotki.yandex.ru/get/5804/5372191.0/0_76c21_4563963e_orig

I change it to http://c.tile.openstreetmap.org/%s/%s/%s.png and get many mistakes like:

  TypeError: not enough arguments for format string, finally I add
  supported_format: image/png but the problem remained.

Can someone explain how it works with tiles ( %(z)s/%(y)s/%(x)s.png ) ?


texts: How can I get a configured Pop-up in ArcGIS Pro to publish to Portal?

texts:
I am using ArcGIS Pro version 2.5. From the pictures below, I have the pop-up configured and it displays correctly in ArcGIS Pro.

I have read that you can publish layers and the configured pop-ups will follow to Portal which my version is 10.7.1.

https://community.esri.com/thread/246120-pop-up-created-in-pro-does-not-work-when-published-to-portal

https://community.esri.com/thread/212571-how-do-you-display-popups-configured-in-pro-with-expressions-to-display-correctly-in-a-webmap

My workflow is like the first post's for publishing.  Share tab -- Publish web layer -- I am referencing registered data -- Publishing both Feature and Map Image.  I do not have version issues like the second post with my Portal being above 10.6.

After the publish goes though, when I open the web map and load the data in, none of the pop-up configuration transfers.  The custom text, selection of fields, and expressions do not publish.  I have tried this with the image layer, the feature layer, registered data, and hosted data.

Is there something that I am missing?

Thanks for your thoughts!


texts: Layers disappearing in QGIS print composer?

texts:
I am using QGIS 2.18.26.

I have problems with a layer disappearing on map in my composer. First I thought it is a scale based visibility  problem, it is a layer with some rectangles and titles on them on a map, visibility is setup to maximum 9 000 000 and minimum 55 000. Playing around with the scale based visibility made the layer appear sometimes but not in any logical way. Then I noticed that clicking of the scale based visibility setting does not seem to help, the layer is still hidden.

The rectangles does not show up. Never when generating an export / pdf / atlas, sometimes though but not often on the print composer screen before generating.

I have tried many different settings, and cannot see any logics in when the layer is shown and when it is not.

Does this ring any bells, any similar situations which have found a solution?

The layer is a locally stored shape layer with rectangles, 12 of them, covering the whole map (rectangles that show the atlas pages on a minimap on the side of each atlas page/big map). The layer is using the same projection as all other layers in the project. AND the rectangles sometimes appear, only to disappear a bit later. The rectangles / the layer in question can be seen in qqis main view as this: imgur.com/3GrMKzK . But in the print composer, no rectangle layer shown: imgur.com/Z1f05fC . No scaled based visibility check-marked now.


texts: Finding the percent of acres per landuse in ArcMap

texts:
I'm trying to add a field to my attribute table that shows the percent of acres per land use. I tried using this approach enter link description here

but I couldn't seen to get it to work for me. I know that I'm missing something so simple but being that I'm new to using ArcMap, I'm having great difficulty.

Table I have:

Table I want:


texts:
I have installed QGIS Server through OSGEO4W. I followed the instructions by Anita Graser, as well as the instructions on this page but still keep getting a 500 Internal Server Error. The objective was to install a localhost instance to work with QGIS Server in a workshop environment.

What is peculiar is that in a group of 13 laptops running Win7, 4 loaded fine (with the odd 500 error).

Is there anyone who has managed to overcome this? I copied the dll files into the QGIS Server folder, changed the filename to remove the .exe from the filename qgis_mapserv.fcgi.exe, and then added it back on again. Every change I made, I restarted the apache server, but no change to the end result. Is it perhaps an install issue, which is making this impossible to troubleshoot?


texts: Add 'name' from GeoJSON to Leaflet List Markers control

texts:
Im my web app I use Leaflet List Markers plugin (https://github.com/stefanocudini/leaflet-list-markers). Markers loaded as GeoJSON from GeoServer. I want dispaly name of markers in list control, but unsuccessfully. Please, help me add name of markers in list control.

Example code:


texts: Generate multiple profile charts from multiple lines on multiple DEM

texts:
I'm trying to figure out how to generate automatically (need to be repeated 2k times) a profile cut of 2 DEMs, every 10 meters, with distance from start as a title and customizable colors/scales as an output. In short, Im trying to make something like this :

Honestly, I'm trying to get Civil3D out of my workflow. Would love a ArcGIS solution but doubt it. QGIS maybe? FME?

I know all the process to get the lines at every 10 meters following my centerline. I just don't know where to go from there to get a result as above.


texts: A tool to “display arrows” in a model in geometric network,

texts:
I am trying to find a tool to “display arrows” in a ModelBuilder model in geometric network.

I couldn’t figure out if there is a tool that can be added to the model to display the arrows as we do with the tool available in the “utility network analyst” bar


texts: SOILGRIDS: Timeout error when linking layer as WMS in QGIS

texts:
I have added clay content as a WMS layer in the database manager of QGIS 3.14. I consistently get the following error:
2020-08-29T10:32:40     WARNING    Network request https://dev-maps.isric.org/mapserv?map=/map/clay.map&version=1.3.0&service=WMS&request=GetLegendGraphic&sld_version=1.1.0&layer=clay_0-5cm_mean&format=image/png&STYLE=default&TRANSPARENT=true timed out
2020-08-29T10:32:40     WARNING    Network request https://dev-maps.isric.org/mapserv?map=/map/clay.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-55.97729999999999961,-179.9979999999999905,82.71930000000000405,179.9939999999999998&CRS=EPSG:4326&WIDTH=1204&HEIGHT=465&LAYERS=clay_0-5cm_mean&STYLES=&FORMAT=image/png&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE timed out

My own internet connection fast and stable, I don't think it is on my end. I can download tiles from https://files.isric.org/soilgrids/latest/data/clay/
but then I cannot find how the tiles are coded (I need a part of Haiti).
My students do not know python or R (yet) so I prefer the QGIS method if possible.


texts: Neighborhood analyses with GRASS/QGIS using ASTER GDEM

texts:
I have a problem related with this question (Old question without answer: R.Neighbors with large neighborhoods?).

Basically, I want to estimate roughness of a given area, for which I was using QGIS, but I want to change the neigborhood size, so I'm using GRASS now because I didn't find a way to do that in QGIS. I am working with ASTER GDEM v2.

So I want to explore a circular neighborhood of aprox. 500m around each pixel and, since each pixel of ASTER GDEM v2 has aprox. 30m of width, I believe my neighborhood size should be of 17, right?

I selected "use circular neighborhood" and for method "average". The program runs and creates a new layer but then I cannot open the raster file. The "Map Display" shows a blank screen when I perform "zoom to selected maps" of the result layer. I tried to save the file as .tiff and open it in QGIS but I obtain, again, a blank screen, and the raster has no values (nan-nan).

Can someone help me please? I don't know if this problem is related with the parameters I am using or if it is related with the software!

UPDATE:
After running r.info, I got this:


texts: How to place markers on the outline of a circle in google maps?

texts:
I have the center of a circle and it's radius and I need to place variable amount of markers on the outline/stroke of the circle. The distance between this markers should kept same, so if I have 8 markers to put on the orbit, each of them has 45 degree space (360 / 8) and also it should not be hard-coded as the number of markers is variable.

I wanted to know what is the right way for calculating the Geo-coordinates of each marker?

here you can see an example, the outer marker should move on, however I can't calculate the proper LatLng values.

I finally came up with this formula:

var dg = 90;
var lat = Math.sin(dg * Math.PI / 180) * r + center_geo_lat;
var lng = Math.cos(dg * Math.PI / 180) * r + center_geo_lng;

If you look at snapshots you can see that markers at 0 and 180 degree has placed quite fine, but it seems there is something wrong with other markers. It's gonna works fine for me, because I don't need such a perfect accuracy, however I'm wondering what's cause the problem?


texts: Adding detailed Google Maps baselayer on Leaflet

texts:
I know that Google map baselayer can be used in Leaflet like this:

googleMap = L.tileLayer('http://{s}.google.com/vt/lyrs=s&x={x}&y={y}&z={z}',{ maxZoom: 20, subdomains:['mt0','mt1','mt2','mt3'] });

The Google Maps base layer is showing usually not all the points of interests as you know in Google lot of points of interests are hidden only when you type them in search they will appear. How to have all points of interests even the hidden one all visible in Google Maps baselayer implemented in Leaflet?


texts: How to create a tileset of custom size in geoserver?

texts:
I am trying to create a tileset of size 200X200. But by default it builts 256X256 tiles. The documentation says, the settings is available under gridsets. So I tried creating a new gridset with the tile size of 200X200.

But where do I use this gridset when creating the tiles?

In the tile creation page, the dropdown to select the gridset, does not list the newly created gridset.


texts: Adding Field using C# winform and ArcObjects?

texts:
How to add a field in featureclass by using C# winform,arcobjects?

It means when we click on command winform should open and afetr field entries it should be add in to table.


texts: Postgis - Converting 'geometry_dump' data type to a 'geometry' type

texts:
I am working with a dataset of property boundaries that has had ST_Union() run on it, thus all the polygons appear as a multi polygon in one row of a table. I would like to separate these out so I get a separate row/polygon for each shape. I have tried running ST_Dump() on this and it seems to separate out the polygons in the dataset but I can no longer view this in a GIS (I'm currently using QGIS). The data-type for this column says it's a 'geometry_dump' and I assume this is why I cannot visualize it. When I access my Postgres database from QGIS this new table is listed under the heading of datasets without a geometry.

Does anyone know if there is something I can do to convert a 'geometry_dump' to a 'geometry'?


texts: Creating table using live data from external API using ArcGIS Pro

texts:
I'd like to be able to relate data from an external, web-based, API to an attribute table on a feature layer. I control the external API, so I can format the data in whatever way is necessary (though JSON would be preferred). I'd like the data to be "live", so that if new data is added to the external source it will also be present in ArcGIS Pro.
Is this possible?
If so, how?
If not, are there any workarounds that will get me close?


texts: Null output from field calculator code block

texts:
I have a simple calculation to compare two columns and populate a new column with the lower value but I am getting null output...any ideas?

expression

re_score(!Signal_Strength!,!Signal_Strength_1!,out1)

code block

out1='99'
def re_score(f1,f2,out1):
 if f1<f2:
  out1 = f1
 else:
  out1=f2
 return out1

  WARNING 000405: No records within table Succeeded at Tue Mar 06
  15:13:45 2012 (Elapsed Time: 7.00 seconds)


texts: Is there a way to tell the supported input formats of libmap.dll without the associated .exe?

texts: I have only compiled Mapscript_csharp and supporting dlls used in a .NET application - not the executables. Is there a way to tell what input / output formats are supported by these dlls?


texts: Graduated colouring Voronoi polygons? (QGis)

texts:
I've got a huge database of depth points for the whole Irish sea which I'd like to be a nice smooth surface, or at least have a depth value associated with the entire area (i.e. include the inter-point spaces). I've tried interpolation but can't get it to work well, so I thought to make Voronoi polygons from the points, and then graduated colour them.

Problem: In 2.1 master, a blank column is created in the Voronoi polygons shapefile, and switching from single symbol:  towards graduated results in a complete blank: . This even happens with the csv file, from dots:  towards blank: . Andre Joost got around a similar problem by using semicolon delimited files, but it feels to me like this is a genuine bug in master, especially given the same files don't blank the Voronoi shapefile in DuFour: . Has anyone else experienced this? I'm about to raise it on QGis hub.

Question: can I colour the Voronoi polygons using the depth value? The Voronoi shapefile has depth as a field but it doesn't do anything - possibly it doesn't know how to deal with the fact that each Voronoi polygon has a number of depths equal to its number of edges (usually 4 in my case here). Can I colour them using the depth value from the points file somehow?


texts: How to make background of EMF symbol transparent in ArcGIS?

texts: I have been sent a collection of PNG and SVG images with which I should symbolize my points in ArcGIS 10.0. Unfortunately ArcGIS 10.0 doesn´t support PNG when symbolizing points, so I have converted them to EMF files. The background transparency of the images are now lost and are displayed black. How can I make the background of these emf images transparent in ArcGIS?


texts:
I want to find a way to run a tool's UI using arcpy. I know that I can use this code to execute buffer tool

import arcpy

roads = "c:/base/data.gdb/roads"
output = "c:/base/data.gdb/roads_Buffer"

# Run Buffer using the variables set above and pass the remaining
# parameters in as strings
arcpy.Buffer_analysis(roads, output, "distance", "FULL", "ROUND", "NONE")

This code executes the buffer itself, but I want just display buffer UI just like when we click on buffer. Is there any way to do it?


texts: Identify intersecting but not bordering polygons using SF package in R

texts:
I am trying to identify the set of polygons (census tracts) that intersect another polygon (a police district).  When I use the st_intersects function, I get all the polygons I expect but also a bunch of polygons on the outside the border. (In the example below, the green polygons are tracts, the red outline is the police district and the polygon labeled A is an example of a tract that st_intersect is including but I would like to exclude).

Similarly, if I use st_overlaps I get tracts that are both on the border outside (A) and inside (B) the district.

Finally, if I use st_contains, I don't get any tracts that border the district boundary, but this excludes tracts (like B above) that I would like to include.

Is there a way to select just those tracts that are contained within the police district but exclude those outside the border? Could I, for example, shrink the police district somehow so that the outside tracts don't touch the police district boundary at all?  Or is there better way to achieve what I'm trying to do?


texts:
I have been instructed to create an small offline application in Android that have to:

a) Use our own map in JPEG format as the basemap (Google Maps / OpenStreetMap is not applicable)

b) Show the GPS coordinates from external Bluetooth GPS (Google map services not applicable too!!)

c) Use Spatialite to build network for shortest path calculation

d) Display the route information on top of the JPEG

I am new to the Android application. I am confused with:

After converting GPS data from WGS84 to Local coordinate system to Image coordinate system (x,y pixels), any JavaScript can help to show it on top of the JPEG file by a marker?
Are there any map viewer available for displaying the georeferenced image file?
How to display the route data (Lines) from Spatialite on top of the map?

I found most of the related Android apps make use of Google Maps or OpenStreetMap, but none of them make use of the image file for positioning and navigation. I am not clearly understanding how can I combine all those of things (JPEG only, Spatialite, GPS data) together, or are there any simple ways to do it?


texts: Importing interchange file (.e00) in ArcGIS Pro?

texts:
I have ArcGIS Pro 1.31 and I found on ESRI website that the "import from E00" tool is no longer available in ArcGIS Pro.

So what's the solution if you want to work with interchange files?


texts: Shapefile for Colombian Cities with District / Borough level outlines

texts:
I am looking for a shapefile for the city of Medellin and Bogota - Where it has the districts/Boroughs outlined

I have been on a both Colombian National gov and Local Gov sites, as well as official websites for the 2 city

But can not find exactly what I am looking for

I want to use it in QGIS, so that I can then place crime data onto it and create a heat map


texts: compute percentage of total area covered by each category of raster

texts:
I have a raster Land use classification file. It has 7 categories: forest, shrubs, snow/glaciers, grassland, water body, barren land, agriculture land. All of them are represented by a number each. So 1 would be forest, 2 would be shurb land and so on.

I need to get the area covered by all of them separately.

  for eg. forest = 40%, shrubs = 10%, agriculture  = 15% and so on or something similar results, in terms of percent of gross value would also work.

I looked into 1 or 2 other questions which recommended to re classify first by assigning 1 for forest and 0 for all others, use zonal stats to find number of cell, then multiply by the cell area. It works but its too lengthy process and besides I have to do this for other rasters as well, not only for this land use file.

Is there a better way to achieve this ?


texts: Using the Make Query Table for Querying a Many-to-Many relationship

texts:
I have two tables that I have connected via an attributed relationship class of the cardinality Many-to-Many.  One is a table of water body (rivers) features and the other is a series of pollutants the water bodies have been assessed for.  All of these tables (river table, attributed relationship class, pollutants table) are located in the same file geodatabase.  The reason for the many-to-many relationship is that one water body can impaired by many pollutants and one pollutant can impair many water bodies.  In ArcGIS 9.3, I was able to use the "Make Query Table" tool to link these tables up and then query based on a single pollutant.  Here is an example of the VBA or SQL code used to do this:

"Rivers_Streams_2014.ComID" = "rc_RiversStreams_2014.ComID_FK" AND
"IR2014_RiversStreams.ID" = "rc_RiversStreams_2014.ID_FK" AND
"IR2014_RiversStreams.Cause_Grouping" = 'Bacteria'

Now in ArcGIS 10.2.1, it keeps giving me an error message stating that

  "An expected field was not found or could not be retrieved properly."

I'm not sure what I'm doing wrong as I've tried tweaking what options I select within the Make Query Table tool about 20 different ways. My guess is that there is a problem with my expression because the "verify" option always tells that it won't work.  (Interestingly enough, it used to do the same thing when I had ArcGIS 9.3 and yet it still worked.)

Basically, what I'm trying to do is have this tool query out only those water bodies that have an assessment for a specified pollutant and be able to display them as a map layer.


texts: ArcGIS Custom Printing Service for a Layout With More Than One Map Frame

texts:
Recently we've been migrating our manual printing workflows to custom printing services at our ArcGIS Server instance. Everything was fun and games until we stumbled upon some Layouts with more than one Map Frame pointing to different Maps. After some googling I find out that the ExportWebMap specification doesn't cover Layouts with multiple Map Frames, forcing us to name the desired one as WEBMAP_MAP_FRAME. We are currently considering writing a more robust Printing Service that will generate its own Web_Map_as_JSON parameters internally resulting in multiple PNGs for each Map Frame and later on binding it all together in a final PDF down the pipeline but... it seems VERY janky!
TL;DR: Is it possible to create a custom printing service for ArcGIS Server to print a Layout with more than one Map Frame pointing to different Map Services?
We are using ArcGIS Pro 2.7.3 and our Enterprise instances are running 10.7.1 and 10.8.1.


texts: Interpolate only no value cells not connected to the egde of the raster

texts:
Is it possible to interpolate geotiff raster no value cells that are not connected to the geotiff border by other no value cells with QGIS or other freely available tools?

The grey areas are gaps. I would like to have the small grey areas surrounded by RGB dots interpolated but leave the grey area in the upper part of the image (the sea) as it is. The image was created with the QGIS Raster -> Analysis -> Fill nodata with max distance=10. If I make the max distance larger, the dots at the coast start to change the sea, too.
QGIS Raster -> Analysis -> Fill nodata does not have such an option.


texts: How to force PostGIS to use a specific GDAL version?

texts:
I have Ubuntu Server (16.04) with installed PostgreSQL (9.5) and PostGIS (2.2).
I have also installed GDAL 2.1 (from ubuntugis/ppa) and I need PostGIS to use it.
But when I check:

SELECT PostGIS_GDAL_Version();

it writes:

GDAL 1.11.3, released 2015/09/16

So, it uses the obsolete version of GDAL that came from somewhere else.

How can I point PostGIS to the proper version of GDAL? Is it a question of environment variables or PostGIS configuration?


texts: Resample fine raster before extracting values to coarser grid?

texts:
I have a raster with resolution 0.05x0.05 deg. I want to extract the values of the raster to a point grid with resolution roughly 14 km (mean latitude is 47 deg), which I have averaged to 0.1575 deg. My concern is that each point represents the whole 14x14 km cell, but extracting values from a finer resolution might not be representative.

Is it advisable to resample the raster to a coarser resolution before extracting the values - even though the data won't fit perfectly anyway due to a different projection - or would it bring even more inaccuracy?

EDIT: my working domain is within 35N 5E (SW corner) to 60N 37E (NE corner). The raster data is in WGS-84 GCS, the point grid has a custom Lambert Conformal Conic projection. I need to assign land cover characteristics to each grid point - I should note that this data is not discrete, it is stored as a multiband raster with percentage values for each band.

EDIT2: Ok, I just noticed that the tool which I am using for extraction (Extract Multi Values to Points) has the option of Bilinear interpolation of values at point locations. Would it possibly be better to use this option instead of interpolating the raster in advance?


texts: Acquire ground coordinates of target while doing aerial survey (airplane or helicopter)

texts:
The Problem:

Each year I board a small airplane and look for diseased individual trees while flying transects over the forest I help manage. The current method for marking a tree's location is manually creating a point in a GIS on a tablet by matching actual ground features with the aerial imagery basemap in tablet. It's not very accurate.

The Question:

Is there a tool that I can point at a ground target (a tree in my case) that will return an accurate pitch and azimuth in a moving airplane?  My plan would be to pair that pitch/azimuth information with the GPS coordinates of the airplane (handheld GPS or embedded in the tool) and lay that offset over a digital elevation model to acquire a ground coordinate for my target.

What I've Already Explored:

TruPulse 360 laser rangefinders: A rep at the company said that their rangefinders would not work from inside a moving airplane.
Drone/UAV-based surveys: A UAV is too slow compared to an airplane, although a good lat/long can be captured by hovering directly over a target and snapping a picture with the camera pointed straight down.


texts: Nominatim reverse geocoding Json format difference

texts:
I am trying to do reverse geocoding using open street maps. I am getting 2 different Json formats in the response.

For, http://nominatim.openstreetmap.org/reverse?json_callback=cb&format=json&lat=12.975360&lon=80.190122&zoom=27&addressdetails=1 The Response contains neigbourhood,suburb etc,.

For,http://nominatim.openstreetmap.org/reverse?json_callback=cb&format=json&lat=40.119198&lon=-75.355764&zoom=27&osm_type=N&addressdetails=1 The Response contains county and other fields.

I do see that the OSM_TYPE is differing in both the response(node/way). Is there any way to restrict the type of response?I tried to add OSM_TYPE in the query string which doesn't seem to work. Is there anyother format that could be received?


texts: ArcGIS for Server - Mosaic Data set layer definition query

texts:
I have a mosaic data set containing four raster images. I have published it on ArcGIS for server as a map service. I am applying Layer definition in export image interface at ArcGIS REST end. I have a text field on which I am applying the definition query. it works successfully in ArcMap but it doesn't work on Arcgis Rest Page of Export Map( showing Definition query is not working). I have tried with correct query and wrong query but it shows the image regardless of definition query.

I am using
ArcGIS API for JavaScript  3.17
ArcGIS Server 10.2
ArcMap 10.4

EDIT

I can not share my arcgis server rest page. however for explaining my problem further i am using esri sample server.

This link shows a basic map export.(without any definition query). Click here to see
The second link shows a map export with definition query applied. Click here to see

See the difference in images.

This is how i want definition query to work, however in my case its giving the same result. Hope now my question is clear.


texts:
I have an image with three bands that are RGB and correspond to percent values from 0 - 100.

I want to create an image where the maximum pixel value amongst all three bands is one single color.

For example, for a specific pixel, if the value in Band 3 (Blue) is highest when compared to the values in Band 1 and 2, then blue is what the color of that pixel will be. And so on for all other pixels.

So far, I understand how to get a single band which simply represents the maximum value for each pixel and I can display that using a singleband gray symbology type.

But, is there a to do is so the maximum value is represented by a color?


texts: How to remove a single marker in Mapbox.js?

texts:
I have a simple Marker added to map like this:
var x = L.marker([50.5, 30.5]);
x.addTo(map);

Now how can I hide/remove x? I searched the API docs for a Marker.remove() method, but couldn't find any.

Even more odd - Google is failing me. All results I can get is how to remove the whole layer with all markers, but not only one.

Am I getting something wrong here about markers? Is there a way to remove only one?


texts: QGIS Plugin Running in QGIS 3.10 but not in QGIS 3.4

texts:
I am new to QGIS, need your help in fixing one problem. I created a Plugin for QGIS. While creating it, I worked on QGIS 3.10 LTR version for which my plugin is working perfectly.

Now when I shared this Plugin folder with my friend who is using QGIS 3.4 LTR version, the Plugin execution is giving error with the message like there was errors executing the algorithm.

I am unable to figure out about exact problem.


texts: Tweaking the table after shapefile import?

texts:
I'm importing shapefiles to Postgresql with the shp2pgsql command line tool.

shp2pgsql -s 3035 -W "latin1" foo.shp | psql

Should I do some table (SQL) tweaks after that? Or is the import fine after that?


texts: Recover M-values from PostGIS `ST_Intersection`?

texts:
I am querying the intersection of a path (LinestringZM) through a region (MultiPolygon) but the ST_Intersection returns only the Lat/Long/Alt coordinates and I need to recover time, stored in the M values.

SELECT path.name AS path_name,
  ST_Length(path.geom::geography) AS distance,
  ST_Length(geography(ST_Intersection(path.geom,region.geom))) AS region_distance,
  ST_M(ST_EndPoint(path.geom)) - ST_M(ST_StartPoint(path.geom)) AS time,
  ST_M(ST_EndPoint(ST_Intersection(path.geom,region.geom))) - ST_M(ST_StartPoint(ST_Intersection(path.geom,region.geom))) AS region_time
  region.name AS region_name
FROM path INNER JOIN region ON ST_Intersects(path.geom, region.geom)
ORDER BY vehicle_name, region_distance DESC;

How can I get the exact elapsed time within a region? (the above region_time is null because M coordinate is dropped. I was thinking of using LATERAL select to get the equivalent matching segment of the original path like so:
ST_SharedPaths(path.geom,rel.intersection_path)
...
LATERAL (SELECT ST_Intersection(path.geom,region.geom) AS intersection_path) AS rel

But that gives me ERROR:  GEOSSharedPaths: TopologyException: found non-noded intersection between LINESTRING (...

I'm using PostGIS 2.4 and PostgreSQL 10.1 beta, in case it matters.


texts: hiding some "feature"s in GeoJSON file on map

texts:
I have geojson file including boundaries of most countries in world.

When I clicked on a "feature" (aka country, or polygon) I want to hide it while keep showing the others. because I am planing to showing more detailed geojson (states, cities) for that clicked country.

I am planning to use Leaflet JS, but I am open to any other libraries as well.

===================

This is the first screen, apply the world countries boundaries.

This is after clicking somewhere on map, I added to US boundaries geoJson to main map variable.

As you can see, I put US geojson upon the world geojson (you can realize that there is still some green area(which comes from world geojson) on the top left part of USA).

What i want is to hide only clicked country, USA in this case, and show more detailed geojson belong to that country.

PS: I dont use tiles, I need to clarify countries by geojson


texts: Clipping street layer with district layer in QGIS?

texts:
What I am trying to do is clip a streets layer within a district using Qgis.

They are in different CRS's and I've tried saving them into the same CRS.

I right click both of the layers and in the properties the show the same CRS.

I then go to "vector" in the menu and then to "geoprocessing tools">> clip.

I specify the layers and the result is another layer that does not show the intended results.

I've been trying to save them in the same CRS and remove the OTF projection.

Is there something that I am missing?


texts: Using the Live tracking with internal GPS?

texts:
We have some Panasonic ToughPads (FZ-G1) with a GPS a module built in.

I know where it is and can locate in it in COM 3 and it connects, registers a signal, shows satellite locations, but refuses to give me a any details. It also doesn't place a cursor where I am located on the map.

How does QGIS decide what Baud rate to use?

Is there a way I can manually set it?

or can anyone recommend a better plug-in?


texts: How do I resolve line outlines jutting into other lines in Mapnik/Cascadenik

texts:
We are designing a basemap, and I am trying to resolve a problem with outlines that comes up in roads. The idea is to have lighter colors roads, but they don't get pushed out far enough with our light background. I've tried various combinations of outline-cap and outline-join, but still run into tiny outlines sticking into the joints of lines/roads.

is there something I am missing to resolve this, or should i just choose a color that stands out above the land?


texts: ArcGIS AddIn with WPF DockableWindow causes ArcMap to not even start up anymore

texts:
I have a simple ArcGIS AddIn with a toolbar and a button on it, which opens a DockableWindow holding a UserControl with some UI Elements on it (a textbox and a listbox). The whole thing was made with the ArcGIS VisualStudio Integration Wizard with the "Windows Forms" Designer type and working as expected.

Now I wanted to port this to WPF Designer type because of improvements in data binding and so on. So I fired up ArcGIS VS Wizard again and created a new DockableWindow with WPF. In Config.esriaddinx I changed the config to point to the new WPF Control instead of the Windows Forms Control, just to test if it works (no more functionality implemented yet).

But when I start ArcMap now, it immediately crashes and Visual Studio Output Windows tells me

The program '[6296] ArcMap.exe: Managed (v2.0.50727)' has exited with code -1 (0xffffffff).

Nothing more there to pin down the reason for this. Any ideas what I am doing wrong?


texts: QGIS Python Composer Attribute Table has no content

texts:
I add to the composer an attribute table and in this table there is no content.

My Code:

    mapRenderer = self.iface.mapCanvas().mapRenderer()
    c = QgsComposition(mapRenderer)
    c.setPlotStyle(QgsComposition.Print)
    x, y = 0, 0
    w, h = c.paperWidth(), c.paperHeight()
    composerMap = QgsComposerMap(c, x,y,w,h)
    c.addItem(composerMap)

    composerLabel = QgsComposerLabel(c)
    composerLabel.setText("Hello world")
    composerLabel.adjustSizeToText()
    c.addItem(composerLabel)

    tableGes = QgsComposerAttributeTable(c)
    tableGes.setVectorLayer(vlayer)
    tableGes.setComposerMap(composerMap)
    tableGes.setDisplayOnlyVisibleFeatures(True)
    c.addItem(tableGes)
    ## Here i try this Methode also:
    ##c.addComposerTable(tableGes)
    ## Same result!

    printer = QPrinter()
    printer.setOutputFormat(QPrinter.PdfFormat)
    printer.setOutputFileName("E:\***\out.pdf")
    printer.setPaperSize(QSizeF(c.paperWidth(), c.paperHeight()),
    QPrinter.Millimeter)
    printer.setFullPage(True)
    printer.setColorMode(QPrinter.Color)
    printer.setResolution(150)

    pdfPainter = QPainter(printer)
    c.renderPage(pdfPainter, 0)
    pdfPainter.end()

Additional I post the out.pdf as a image.

Can someone try the code in their own plugin?

If there is the same result i can make a bug report.

If not, I need to search for the bug in my plugin.


texts: Installing Pythonwin for use with ArcPy?

texts:
I have ArcGIS 10.0
And installed pythonwin 2.6.5

When wrote a script in Python Shell - it works just fine
But in pythonwin - i get import arcpy error as below

>>> import arcpy
Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
ImportError: No module named arcpy


texts:
I'd like to know how do I export a map from Gis to a dwg file (to use it in autocad for architectural needs).
I've tried to watch already a couple of tutorials and also registered to ARCGIS free trial but still can't get it.
Somehow, after I've registered to ARCGIS, I have only a user for the site but don't have any software which I can access from my desktop.
Does someone can help me understand how to I do it?
How can I save a map from Gis and to what format?
And later on, how do I transform this format to a dwg?
Finally, is it necessary to install same plugin in autocad for this purpose?


texts: Why are my points offset when using the Go to XY tool?

texts:
I am creating points in ArcGIS 10.1 to add to a GPS so I can navigate to them in the field.

I add my points by editing an empty shapefile created in ArcCatalog. I then Add Field: Lat and Long as "Float" points. I then use the "Calculate Geometry" to add X and Y coordinates in decimal degrees. The points layer and map layer are all WGS1984.

To check the coordinates I used the "Go to XY" toolbar and the points created by this tool are always offset from the shapefile points by 20-40 m.

Does anyone know the reason for this and how it can be fixed?


texts: Make selected polygon different color without covering everything?

texts:
Working in ArcGIS Desktop 9.3

I'm selecting polygons by attributes and trying to have the selection show up as a light yellow. I set the selection color to the light yellow and it's the 'lowest/last' layer but it covers all of the other layers. Is there a way to change it?


texts: Using .NET Generics to store ArcObjects?

texts:
What are the implications of using Generic types in .NET to store ArcObjects, such as ILayer, IField, etc?

The compiler throws me a warning when using these kind of values

public class Foo
{
    private List<ILayer> fooLayers;

    public List<ILayer> FooLayers
    {
        get { ... }
        set { ... }
    }
}

Here is the warning:

  Warning   15  Type library exporter warning processing
  'ArcMemorialCore.Topography.IMemorialDocument.set_ProfessionalsEnvolved(value),
  ArcMemorialCore'. Warning: Type library exporter encountered a generic
  type instance in a signature. Generic code may not be exported to
  COM.  ArcMemorialCore

My concern is not only the warning's existence, but of good design practices, performance, etc.


texts: Querying a Styled Layer Descriptors file for OS MasterMap area fill data

texts:
I am attempting to write a software to import OS MasterMaps into our product and am having difficulty getting the area styles to import correctly.

How do I get the fill color for a combination of
osgb:descriptiveGroup + osgb:descriptiveTerm + osgb:make?

When I look at the public domain style files Ordnance Survey provides (e.g. topographicarea-outdoor.sld) I do not find these three descriptions.

How would I determine what fill color to use for a combination of these three values? The .sld files make no reference to osgb:descriptiveGroup, osgb:descriptiveTerm, or osgb:make. All I see is a  tag with a description of the Feature Type and its  color.

For Example, the OS MasterMap (a GML format) has an entry with this information

 <osgb:descriptiveGroup>Tidal Water</osgb:descriptiveGroup>
 <osgb:descriptiveTerm>Foreshore</osgb:descriptiveTerm>
 <osgb:make>Natural</osgb:make>

My quandary is how to find out what color to fill the area by querying a Styled Layer Descriptors (SLD) file.


texts: OpenLayers 3 / Geoserver - add popup w/ attributes to Image layer?

texts:
I have a web app using OpenLayers 3, with several layers added as image layers, like so:

 mylayer= new ol.layer.Image({
     source: new ol.source.ImageWMS({
         url: 'http://server/geoserver/opengeo/wms',
         params: { 'LAYERS': 'layername' },
         serverType: 'geoserver'
     })
 });

I would like to have an info window popup on the click of a feature , I am following the example here : http://openlayers.org/en/v3.11.1/examples/icon.html

I have Bootstrap and JQuery in my app (these are both needed for this example) and am not getting any errors, but no results either..

..this example is for a vector layer, and so far I haven't been able to get this to work for me...Is there a way I can do this with an image layer? or do I need to add the layer as a vector layer?


texts: Is it possible to calculate statistics across multiple fields in arcpy?

texts:
The summary statistics command can be used in arcpy to calculate statistics - this tool runs against a single field, eg:

  SUM—Adds the total value for the specified field.

Is there an out-of-the-box method to calculate statistics across multiple fields?

eg I'd like to calculate the MIN, MAX, MEAN and STDEV of Field1, Field2 and Field3, combined:

OID  Field1  Field2 Field3 Field4
 1     23      24     25     26
 2     12      13     14     15
 3     17      18     19     20

One approach is to read through the table with a cursor, and calculate these values manually in a script. Is there a better/faster way?

Edit: I need this to run on only the selected records, and these will change frequently.


texts: Get the spatial reference of an ArcGIS layer with C#

texts:
I'm using the following C# code to add my layer dynamically on an ArcGIS silverlight application:

ArcGISDynamicMapServiceLayer mylayer = new ArcGISDynamicMapServiceLayer()
{
   ID = layerName,
   Url = mapServiceUrl,
 };
 mylayer.Initialize();
 mylayer.Initialized += mylayer _Initialized;
 mylayer.InitializationFailed += mylayer _Failed;
 myMap.Layers.Add(mylayer );

But when I try to get the spatial reference of the layer using:

SpatialReference sp = mylayer.SpatialReference;
int testww = sp.WKID;

I get an error "Object reference not set to an instance of an object." because mylayer.SpatialReference is Null.

Where am I going wrong?


texts: Unable to delete .shp and .dbf files from directory when on canvas PyQGIS

texts:
I am creating a standalone PyQGIS application and is met with a problem.
When a shapefile is shown on my canvas, I am unable to delete the corresponding .shp and .dbf files in my directory.
This is the error I get:
PermissionError: [WinError 32] The process cannot access the file because it is being used by another process: 'Data\\tmpy187h2zs\\0-Boundary\\0-Boundary.dbf'

This does not happen when I try to close raster files.
I am using PyQt for the GUI. And there's a workaround to this problem by having 2 buttons. Remove from canvas and Delete.
After pressing the Remove from canvas button, I am able to delete the layer. However, when I combine both of these buttons together, an error occurs.
Here is the function for remove layer and delete:
Remove layer:
def removeAllLayers(Main):
    layers = Main.layers
    registry = QgsProject.instance()
    if len(layers) != 0:
        layers = []
        registry.removeAllMapLayers()
    Main.canvas.refresh()

Delete:
def deleteFiles(Main):
    for files in os.listdir(Main.tempPath + '/' + button):
        os.remove(Main.tempPath + '/' + button + "/" + files)
    print("Files Deleted")

When combined:
def removeAndDeleteFile(Main):
    removeAllLayers(Main)
    deleteFiles(Main)

and these functions are just embedded into QPushButton:
pushButtonRemoveFromCanvas.clicked.connect(lambda: removeAllLayers(self))
pushButtonDeleteFiles.clicked.connect(lambda: deleteFiles(self))
pushButtonRemoveAndDelete.clicked.connect(lambda: removeAndDeleteFile(self))

This works:

Click pushButtonRemoveFromCanvas
Click pushButtonDeletFiles

This does not work:

Click pushButtonRemoveAndDelete
Returns PermissionError as mentioned above

The main problem stems from that the .shp and .dbf files are not deleted. I have browsed other threads such as can't erase dbf and shp pyQGIS and Deleting raster file after doing some processing in QGIS? but they do not directly answer the question and is a workaround specific to their question.
I have been on this problem for weeks and is unable to solve it, wonder if anyone can help me.


texts: Feature class and field names contains blank records in text file

texts:
I want to print all feature class names along with field names from a GDB, which are having empty/blank records in text file for user. Fields in feature classes are not same. So I want to skip some particular fields which are not requited to check for empty records. Below is my code and I don't know how to skip some fileds
import arcpy
from arcpy import env

env.workspace =  r"C:\Users\Test\Basemap.gdb\ROAD"

fclist = arcpy.ListFeatureClasses()

# want to find empty records for below list of feature classes
fields = ["S_CODE","CATEGORY"]

f=open(r'C:\Users\test\List_GDB\Empty_rec.txt','w')
for fc in fclist:

    bad_list = []

    for field in fields:
        #print (fields)
        with arcpy.da.SearchCursor(fc, field) as cursor:
            for row in cursor:
                if row[0] in ["", None, " "]:
                    bad_list.append(field)
                    break

    f.writelines(fc+"contains empty records in"+str(bad_list)+"\n")

f.close()

Above is GDB structure. Please suggest  how to skip particular fields (Ex: CODE, ADDRESS etc) in all feature classes while checking for empty records.


texts: Getting selected column/field for graduated color via expression or other routine in QGIS?

texts:
I have choropleth map with graduated colors. The column - the field name - that I choose to define the classes and the colors is important for my work on the maps. Is there a method to make it visible or refer to it at any other part of QGIS. It would be especially useful in the print composer, where I can check on other map objects and labels when I use several print composers. As far as I explored, there are only expressions to get values from the attribute table of certain layers or to get other layer properties than the style.

Any ideas?


texts: Hosted tiled mapping (Portal) VS. Cached map service (Server)

texts:
I'm trying to create a detailed building basemap for the state of Washington in the Portal for ArcGIS on premises. In order to make the basemap moving faster, and considering it's large size and static, I simply chose to use tiled mapping service or cached map service. However, there are some concepts I am not very clear, especially about the difference between those two services.
Both services have a function to meet my requirement from Esri's documentation:

Hosted tiled service:

  Tile layers published to and hosted on your portal support fast visualization of large datasets using a collection of predrawn map images or tiles. Hosted tiles provide geographic context for your operational layers. For example, you might include tiles of streets in your neighborhood to provide a visual reference for the street signs in your feature layer. (http://server.arcgis.com/en/portal/latest/use/publish-tiles.htm)

Cached map service:

  The best way to create fast map services is to cache them. When you cache a map service, the server draws the map at a set of scale levels that you define and saves the images. When the server receives a request for a map, it's much quicker to return one of these cached images than to draw the map again. (http://server.arcgis.com/en/server/latest/get-started/windows/tutorial-creating-a-cached-map-service.htm)

Here are five (5) questions:

What's the difference between those two methods? It seems both of
them can manage the caching properties in the Service Editor; while
the hosted tiled mapping is published thru the Portal, while the
cached map service thru the server. Some people say they are the
same (https://geonet.esri.com/thread/79128).
Moreover, publishing thru the Portal only has two options for the
Capabilities: Tiled Mapping and Feature Access; while there are more
capabilities choices by publishing thru the Server. Does that mean
publishing thru the Portal is less functional than publishing thru
the Server?
And why there are both "Tield Mapping" choice and "Caching" for
publishing thru the Portal; while only one choice, Caching, is
availabe for publishing thru Server? Are they the same?
Does the federation affect these two methods, publishing thru
Portal or Server?

Service Editor when publishing a hosted tiled mapping (Portal):

Service Editor when publishing to a Server:

After publishing the tiled/cached service, there is an option to View Cache Status in "My Hosted Services"; while there is no such option in "GIS Servers" in the Catalog. How to understand this "inconsistence"?

My Hosted Services (Portal for ArcGIS):

GIS Servers/Hosted folder:

Here are answers I got from Esri support, and you can reopen this post if it will be helpful for other GISers.

One is hosted, the other is published to server. Essentially the same, depending on how you want to consume them.
No. Server has more capabilities because Portal does not support certain capabilities for hosted services (i.e. KML, WMS, WFS, etc).
Portal is preconfigured with basemaps (like AGOL).
Usually use feature services to work with Portal (add FS to web map)
Tiled mapping would be to include a custom basemap
Server can have a map service enabled without caching (dynamic) – uses DataStore.
Portal will host it on Server.
It can. All existing map services will appear on Portal and you can use those map services dynamically.
Cached map is only being hosted on Server. Server will not manage items published to Portal. Only serves as a reference for the service, not as a manager. Portal overrides most Server functionality.


texts:
I need set the new created fGDB as current for save all FC:

import arcpy
import os
out_folder_path = arcpy.GetParameterAsText (0)
out_name = arcpy.GetParameterAsText (1)
arcpy.CreateFileGDB_management (out_folder_path, out_name)
arcpy.AddMessage ("fGDB created")
# Set fGDB created as scratch workspace
arcpy.env.scratchGDB = out_folder_path + "/" + out_name + ".gdb"
arcpy.AddMessage ("fGDB set current (scratchGDB)")

But I receive the next error:

Traceback (most recent call last):
  File "H:\Python_ArcGIS_MasterSIG\Proyecto_Final\Codigo\Crear_fGDB.py", line 30, in <module>
    arcpy.env.scratchGDB = out_folder_path + "/" + out_name + ".gdb"
  File "c:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\geoprocessing\_base.py", line 541, in set_
    self[env] = val
  File "c:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\geoprocessing\_base.py", line 601, in __setitem__
    ret_ = setattr(self._gp, item, value)
AttributeError: Object: Environment <scratchGDB> cannot be set

Failed to execute (CrearfGDB)

What is the mistake?


texts: Controlling line width at different points in QGIS

texts:
I have vector data for rivers including their width at various points. From what I have seen so far, QGIS 3.8 can set the line width only per linestring, no per segment.

Is it possible (maybe with a Plugin) to control the width of my rivers as they become larger (after inflows), or do I have to split them up into several segments (which messes with feature selection, which is why I'd like to keep each river as one linestring).


texts: Boolean Intersection Filtering in Carto.js

texts:
Context
I have a neighborhood dataset indicating boolean values by building lots for various indicators (i.e. within a special district, 10-minute walking distance to a subway station, below the average median income, etc.). The goal of this interactive map is to create a check-box for each of these indicators, and once checked, will display all TRUE values for that indicator. When more than one indicator is checked, it will find the boolean intersection value and display building lots where TRUE for all checked indicators.
Steps to Reproduce

Create input type in HTML <input type=”checkbox” name=”[column name]” id=”[column name]” value=”T” unchecked>
Create filter.Category(‘[feature name]’, {in: })
Add filters source.addFilters([feature1, feature2, feature3, ...])
Apply filters for each feature
Apply registerListeners()

Current Result
Console.log prints the selection of features only selected from the check-box and pulls in TRUE values for these features as an array. When 2 features are selected:
(2) [Array(2), Array(2)]
> 0: [“featureName1”, “T”]
> 1: [“featureName2”, “T”]

However, the map keeps showing the intersection of all features (whether selected or not selected) in the database.
Expected Result
When one check-box is selected, the map will display all TRUE values for that feature. When multiple check-boxes are selected, the map will only display values that are TRUE for all checked features -- intersection of TRUE for all features.

Console.log seems to be separating out each feature that is selected but the map does not show the corresponding results.
Is the issue with the way addFilters function is formatted?

Code
HTML:
<input type="checkbox" name="featureName1" id="featureName1" value="T" unchecked>
<input type="checkbox" name="featureName2" id="featureName2" value="T" unchecked>

Script:
const source = new carto.source.SQL('SELECT * FROM dataset');

function getSelectedTypes () {
      const envCheckboxes = [...document.querySelectorAll('#controls input[type=checkbox]')];

      const values = [];

      envCheckboxes.forEach(input => input.checked ?
                                 values.push([input.name, input.value]):null);

      return values;
    }

const filter1 = new carto.filter.Category('featureName1', {in: getSelectedTypes()});
const filter2 = new carto.filter.Category('featureName2', {in: getSelectedTypes()});

source.addFilters([filter1, filter2]);

function applyFilter1 () {
    filter1.setFilters({in: getSelectedTypes()})};
function applyFilter2 () {
    filter2.setFilters({in: getSelectedTypes()})};

function registerListers() {
    const envCheckboxes = [...document.querySelectorAll('#controls input[type=checkbox]')];

    const value1 = [ ]
    envCheckboxes.forEach(input => {
    if (input.name === ‘featureName1’) {
    value1.push(input)
    }});
    const value2 = [ ]
    envCheckboxes.forEach(input => {
    if (input.name === ‘featureName2’) {
    value1.push(input)
    }});

    value1.forEach(
    input => input.addEventListener(‘click’, () => applyFilter1()));
    value2.forEach(
    input => input.addEventListener(‘click’, () => applyFilter2()));
}

const layer = new carto.layer.Layer(source,style)

client.addLayer(layer);
client.getLeafletLayer().addTo(map);
registerListeners();


texts: Drawing points in layer given its coordinates using PyQGIS

texts:
I want to draw points in a certain polygon layer. Specifically, I want to display the centroid over each respective feature. All this using PyQGIS.
So far, none of the answers I found helped me. I am using QGIS 3.10.
I am a newbie here. Should I create another layer? Should I edit the current layer? How do I do this?


texts: Delineating small scale / subwatersheds using ArcGIS Desktop

texts:
I have a number of data points showing the placement of stream gauges. What I'd like to do is to delineate the subwatershed for each point, so that when rainfall events occur, I can see the drainage area for water that will travel through that stream gauge.

I have my clipped DEM, to my area of study (in red), as well as the first stream gauge point I'd like to use for delineation (green):

However, I've tried a handful of guides for this delineation, and each time the output is not correct.

I believe that my main issue is that when the flow accumulation is performed, it is only picking up on the largest nearby river, and not some of the smaller tributaries. The gauge/pour point I'm interested in is very far from the nearest accumulation point:

Some of the steps I've tried:
Clip DEM, Fill DEM, Flow Direction, Flow Accumulation, Pour Point (using the green gauge point as the pour point), and Watershed.

I've also followed this guide before:
https://imap.maryland.gov/Documents/Training/LiDAR_Mod7_DelineatingWatersheds.pdf

Any suggestions on what I'm doing wrong, or maybe someone has a better guide/ set of steps to get this done?


texts: Determining if mouse click falls within polygon using ArcObjects?

texts:
I have a polygon drawn on top of a layer using the graphicsContainer of the current activeView, then I want to be able to show a custom context menu if the user right clicks INSIDE the polygon. So I have to figure out a way to determine if the mouse click falls within the polygon. I have tried to use the IRelationalOperator2 as suggested in this question Does anyone know a function to determine if a Point exists within a polygon?, but I can't seem to get it working properly. Here is what I have so far:

System.Drawing.Point screenMousePoint = System.Windows.Forms.Form.MousePosition;
IPoint screenPoint = new Point();
screenPoint.X = screenMousePoint.X;
screenPoint.Y = screenMousePoint.Y;

IPoint mapPoint = this.controller.GetMapCoordinatesFromScreenCoordinates(screenPoint, activeView);

// this.currentDefinedArea is a polygon
IRelationalOperator2 relation = this.currentDefinedArea as IRelationalOperator2;

if (relation.Contains(mapPoint))
{
    // Show context menu
}

The "GetMapCoordinatesFromScreenCoordinates" method looks like this, and it's one of the snippets that is included in arcobjects:

public IPoint GetMapCoordinatesFromScreenCoordinates(IPoint screenPoint, IActiveView activeView)
    {
        if (screenPoint == null || screenPoint.IsEmpty || activeView == null)
        {
            return null;
        }

        IScreenDisplay screenDisplay = activeView.ScreenDisplay;
        IDisplayTransformation displayTransformation = screenDisplay.DisplayTransformation;

        return displayTransformation.ToMapPoint((System.Int32)screenPoint.X, (System.Int32)screenPoint.Y);
    }

I found out what was causing the above code to fail. Using:

System.Drawing.Point screenMousePoint = System.Windows.Forms.Form.MousePosition;

...gives a different set of X, Y coordinates than using the MouseEventArgs provided by the OnMouseDown method. Don't know why I wasn't using them in the first place...here's how my working code looks:

protected override void OnMouseDown(ESRI.ArcGIS.Desktop.AddIns.Tool.MouseEventArgs arg)
{
    if (arg.Button == System.Windows.Forms.MouseButtons.Right)
    {
        if (this.controller.CurrentDefinedArea != null)
        {
            IPoint screenPoint = new Point();
            screenPoint.X = arg.X;
            screenPoint.Y = arg.Y;

            IPoint mapPoint = this.controller.GetMapCoordinatesFromScreenCoordinates(screenPoint, activeView);

            IRelationalOperator2 relation = this.controller.CurrentDefinedArea as IRelationalOperator2;

            if (relation.Contains(mapPoint))
            {
                this.controller.CreateContextMenu(ArcMap.Application);
            }
        }
    }
}


texts: Updating sources with arcpy in ArcGIS Pro

texts:
I'm attempting to update a script that I used in ArcMap for ArcGIS Pro. Getting an AttributeError: 'list' object has no attribute 'ListLayers'
Any ideas? I based this on the code snippet here: https://pro.arcgis.com/en/pro-app/arcpy/mapping/updatingandfixingdatasources.htm
import arcpy

aprx = arcpy.mp.ArcGISProject('CURRENT')

lyr = aprx.listMaps('GeologicHazardMap').listLayers('qfaults')[0]
cp = lyr.connectionProperties
cp['connection_info']['database'] = 'N:\\ReferenceMaterials\\ArcGIS Python\\Current Catalog\\Test Data'
cp['dataset'] = 'FaultsReprojected.shp'
lyr.updateConnectionProperties(lyr.connectionProperties, cp)


texts:
I'm fiddling around trying to create a thumbnail of a fairly large (originally 27MB) raster, chunked into 64x64 blocks. I've got the following and it works, quite well, but I'm wondering if there is a way to make it more efficient, or if it is just totally mental and shouldn't be attempted.

It suffers from various flaws - that it is not really so reusable, because it is quite specific, for example.

SELECT ST_AsPNG(ST_Reclass(ST_Union(rastnxn), 1, '[0-10]:1-65535', '16BUI', 0))
FROM (
        SELECT filename, ST_Resample(rast, 8, 8) as rastnxn
        FROM rasters
        WHERE filename = 'someraster.tif'
) as subset
GROUP BY filename;

What do you think? Are there other (better) ways to quickly get the original raster back out in a 'thumbnail' style?


texts: arcpy.FieldMappings renames input fields internally

texts:
I'm trying to make use of Field mappings to merge many datasets into a new one... as per the manual.
   ...
   print (type(idealFieldList))   # class 'list'

   # let's prune the fluff

   fieldmappings=arcpy.FieldMappings()
   fcSources = arcpy.ListFeatureClasses("Source*", "Polygon")
   for fc in fcSources:
       print (fc + " is being absorbed..", end='')
       fieldmappings.addTable(fc)       # but something else is happening here....

   fieldMaps=fieldmappings.fieldMappings
   print ("Done. \nAll fieldmaps have been acquired")

   # let's drop the unwanted fields
   for eachfld in fieldMaps:
       outf=eachfld.outputField
       if outf.name in idealFieldList:
            #do nothing
            print (outf.name + " ------------> is ok")
       else:
            print (outf.name + " must go")

OUTPUT:
ContactName_1 must go
ContactName_12 must go
ContactName_12_13 must go
ContactName_12_13_14 must go
ContactName_12_13_14_15 must go
ContactName_12_13_14_15_16 must go

The problem is this line fieldmappings.addTable(fc) # but something else is happening here....
I have 6 source feature classes, all of them have a field called "ContactName". In my fieldMappings, after they have been added, I find the field name has been altered as shown above in the output.
I thought the whole point of FieldMappings was that I would have a single output fieldname "ContactName" and that single output would have 6 inputs that are it's origin.
How do I do it differently?


texts: Cannot install ArcGIS for Desktop 10.2 following install guide on Windows 8

texts:
I'm trying to install ArcGIS for DEsktop 10.2, I get to the following screen:

But any of the instructions highlighted in blue do not open, the mouse hovers over them changing their colour to darker blue as if you can click and open, but nothing opens. I can't go through with the install. I thought there might be a problem with the download so I downloaded the file 'Arc GIS for desktop 10.2, from esri online again, using my authorisation codes. The file is 2GB. After downloading the program for install, the only .exe file is written ESRI.exe. I click this and receive the screen as shown in the screenshot above in my post. Both times I have not benn able to continue with the install as not setup links open. There has never been another version of ArcGIS for desktop on this computer.

I have CORE i7, 64 bit, Windows 8, 16GB RAM.

I did simulataneously download the ArcGIS desktop tutorial and the DBMS files, I didn't run any of those. I read the .net 35 has to be installed first, when I ran that, no options opened in any windows, other than to say can this program modify your computer? to which I answered yes. Here is a screen shot of all the downloads I did, I just want to get the ArcGIS for desktop going. I'm completely inexperienced:


texts: Google Earth: Circular polygons that follow earth curvature

texts:
I'd like to draw filled circles in Google Earth which follow the curvature of the earth, i.e. for 1000ft altitude, not only the circle's edge but also the center is required to be 1000ft above reference. At the moment, the circle's edge is at the required altitude, but - not unexpectedly - the earth's surface "pokes out" in the center if the circle gets too wide.

I think I can approximate what I need by creating polygons which are made of concentric circles every so-many distance units, so that all in all, they create a shell-like or hat-like shape which fits above the local earth surface at the required altitude.

However, I can't seem to get the kml geometry for my polygon to work for this. I gather that I don't need an inner boundary, as there will be no holes. Creating outer boundaries with linear rings in various combinations only accepts the last entered linear ring, though, which is the outermost; the inner ones are ignored.

Can someone help? Is it clear what I am trying to achieve?


texts: Creating many lines from point feature using ArcGIS Desktop?

texts:
I would like to create many line from a point feature, or a table (I created the feature with 'Display XY Data'). I was using 'points to line', but it works only for a single line. In the point feature I have a field where is described for each point to which line it correspont. But if I put this lineID in the 'Line_Field', it says me:

  ERROR 000972: Table must have Object IDs in order to be sorted.

but the table has the ObjectID, and if I let the 'Line_Field' empty it works but than I have only a line.


texts: How to calculate average grade over roads from DXF input?

texts:
I receive a dxf file containing percentage grades of a road. It comes from lidar data and it's very accurate. The problem is percentage is included into points and the line is just a line. What I would like to achieve is to calculate longer stretches of a road and get the average grade.

I don't know where to start.


texts:
I am a bit confused on the releases naming of QGIS.

On the website you read everywhere QGIS 2.18 = new LTR (that I suppose stands for "Long Term Release"), QGIS 2.14 = previous LTR.
Neverthless the *-ltr linux repositories point to QGIS 2.14, and so do all the ubuntugis repositories.

So, is it safe to install on ubuntu 16.04 the 2.18 version?

Will it be supported for some time ?

Which repository should I use (I would like to install also GRASS)?

If I use the *-ltr repository, will I soon (when?~) get the 2.18 with a simple apt-get update ?


texts: How to activate GRASS GIS 7 plugin in QGIS 2.8?

texts: Apparently, my linux ubuntu setup of QGIS doesn't contain grass gis, although when I initially installed it last week I added grass gis to the command line. Anyway, I typed in once again and installed grass gis 7.0, however while GRASS starts, the plugin QGIS for it doesn't show up. Any ideas?


texts: How to create river polygon from border's polylines?

texts:
I got hydrography data where rivers are represented by polylines of their borders :

I would like to symbolize these rivers with blue filled polygons, so I want to convert into a polygon layer.
Does anyone have an idea of how to do it?

In the feature class attributes, there's only the name of the river in common between both borders. Polylines are not closed as the image above, so a simple Feature To Polygon doesn't work.


texts: proper use of python subprocess for GDAL

texts:
I need to process some sentinel-2 files which are in jp2 format using python. I first want to clip the jp2 files and save the result as geotiff with the same CRS. But some of the modules I use (e.g., rasterio for clip) don't like jp2 format so I try to convert the jp2 to geotiff. I am able to do so using QGIS translate tool which is essentially GDAL translate. However, when I try to use the translate tool in python:

import geojson, gdal, subprocess

tif_file = r"C:\Users\Administrator\L2A_B04_10m.tif"
jp2_file = r"C:\Users\Administrator\L2A_T29RNP_20170422T110651_B04_10m.jp2"

args = ['gdal_translate', '-of', 'Gtiff', jp2_file , tif_file]
subprocess.Popen(args)

I get a file without CRS. I tried to embed the CRS in the args but I can't figure out a way to do so, I don't get any error but it seems that the code is ignoring my tries, like:

args = ['gdal_translate', '-a_srs','ESPG:4326', '-of', 'Gtiff', jp2_file , tif_file]
subprocess.Popen(args)

results: nothing happens, no error, and no file is created

also:

args = ['gdal_translate', '-a_srs ESPG:4326', '-of', 'Gtiff', jp2_file , tif_file]
subprocess.Popen(args)

also:

args = ['gdal_translate', '-a_srs', '+proj=longlat +ellps=WGS84', '-of', 'Gtiff', jp2_file , tif_file]
subprocess.Popen(args)

and:

args = ['gdal_translate', '-a_srs +proj=longlat +ellps=WGS84', '-of', 'Gtiff', jp2_file , tif_file]
subprocess.Popen(args)

What am I missing? or, is there another easier way to clip jp2 format using python?

I'm using python 2.7 on anaconda on win10, gdal 2.2.3


texts: which framework is good for closed intranet

texts:
My requirement is simple. By latitude/longitude, draw icon on the map, draw a line between two icons.I may need put a customer picture over the map.

It is not difficult by google map javascript or bing map javascript.

The problem is that the system is running in a closed intranet environment. No internet access.

Which framework shall I choose? Does it provide map? The system only focus on a specific area, such like a national park, it does not need the map for the whole world.


texts: Getting Mapbox GL addlayer 'in' filter to work when match target is comma separated list?


texts: Finding correlation between two rasters in ArcGIS for Desktop?

texts:
I am working with a NDVI image and a an interpolated yield map in ArcGIS 10.3.1. They are both projected in the same coordinate system and the cell sizes are equal. I interpolated the yield harvester points using IDW interpolation method. The final product I want to achieve is the correlation (R-squared). I am new to GIS and have been stuck on this for a while trying to use various spatial analyst tools in ArcToolbox but have been unsuccessful. I have clipped the interpolated yield map to another polygon that I have created, however it does not overlay the NDVI completely.

Is it possible to extract the overlaying pixels values and import this into R software where I can perform further stats?


texts: Importing long text from Excel using ArcGIS Desktop always shows up as BLOB?

texts:
ArcGIS Desktop 10.6.1

I am trying to load long text (>255 characters) into my GIS data. The destination field (let's call it "Field Alpha") is a string field with a length of 2000 characters, in a feature class in an enterprise GDB.

My source data is from an Excel file. Let's call the source long text field "Field Bravo".

I intend to load the long text data by joining a table (Excel document or GDB table) to the feature class, and then using field calculator to populate Field Alpha from Field Bravo.

However, no matter how I go about it, when I get the data into ArcGIS it converts to BLOB. If I simply add the Excel doc to an MXD, Field Bravo shows up as BLOB. If I try to import the Excel document into a new geodatabase table (in a separate, file Geodatabase), it will import as BLOB, even if I stipulate in the import that the field is "Text" and designate a length longer than 255.

Is there any way to import text longer than 255 characters into ArcGIS 10.6.1?


texts: What are the min and max values of Map.addLayer on Google Earth Engine?

texts:
I'm starting to play with Google Earth Engine. Following https://developers.google.com/earth-engine/image_visualization I see its possible to set min/max values.

May I ask the meaning of this and when should I use them?


texts: Searching semantically for symbols in ArcGIS Desktop?

texts:
Sorry for the noob question, but is there any resource (on the web) where I could more or less search semantically for symbols to use in an ArcMap layout?

Using the Symbol Selector, I have to scroll through all the symbols from different Style References. I assume there might be one which is fitting my purpose, so is there a more intelligent way of finding it?

(Current case: have to add some measurement towers to a map.)


texts: Defining projection, reprojecting, and wrong area values

texts:
My data: http://water.usgs.gov/GIS/dsdl/aquifrp025.tar.gz (USGS circa 2003 data on aquifers in the contiguous U.S.) Program: ArcGIS 10.1

The shapefile has no coordinate system or spatial reference. The metadata specifies NAD 1983 and in decimal degrees. My goal is to render this in such a way as to obtain 'normal' area/length/etc. values in meters for the purposes of calculating overlays with other spatial data. I cannot project the data because it has no coordinate system (error: no input coordinate system specified; paraphrasing)...batch/define doesn't work correctly either (the former same error, the latter (technically for use with rasters?) defines the projection, but reprojecting using the project tool doesn't do anything substantive (i.e. areas for huge swaths of the U.S. "46.1" etc. etc.).


texts: Set mandatory output format and folder in graphic modeler

texts:
I'm working with the graphic modeler in QGIS3.2 Bonn. I have created some models and they work quite well. In order to make them more easy to use for others:

Is it possible to set the output format of the result of a model mandatory to ods or shp format?
And is it possible to set an output folder mandatory?


texts: POSTGIS: lwgeom_clip_to_ordinate_range clipping routine returned NULL

texts:
I'm trying to do linear referencing with st_locatebetween.

In the past everything went well.  But from now on, i get following issue:

  lwgeom_clip_to_ordinate_range clipping routine returned NULL

Searching on the internet gives me a few results pointing to: results.
I don't even know what it means.

I'm using following code to generate the linear reference:

SELECT
a.geometry,
a.ws_oidn,
a.type,
a.morf,
a.wegcat,
a.rstrnm
   FROM ( SELECT st_locatebetween(ws.geom, wv.vanpos, (wv.totpos + 10)) AS geometry,
        ws.ws_oidn,
        wv.type,
        ws.morf,
        ws.wegcat,
        ws.rstrnm
       FROM wegsegment ws
         LEFT JOIN attwegverharding wv ON ws.ws_oidn = wv.ws_oidn) a


texts: Scaling SVG symbols which are not label background on zoom level in QGIS?

texts:
I need to scale the SVG symbols on zoom level.

@Underdarks answer in the posted question did not help me thus I can't set the SVG symbols as label background for I have labels which should not be overlaying the SVG.

How to scale SVG symbol based on zoom level in QGIS?


texts: How to "clean" and vectorize a raster map?

texts:
I'd like to make use of a geological map as shown in the example below, which is only a pretty small and downsized part of the map, so editing/redrawing doesn't seem to make sense in the end.

My idea would be to have it first converted to vector with only bedrock and tectonic faults left, thereafter having the possibility to optionally use a "cleaned" raster file.

Any hints on a "how to do?"


texts: Determining distances between many individual points using ArcGIS for Desktop?

texts:
I am trying to find an easy way to calculate linear distances between an array of points.  I have an array of points that border a rectangular area that represent a dewatering well array.  As of now, I simply export the XY Coordinates and drop them into an excel sheet to calculate this but was hoping an all-in-one use might be possible within ARCGIS 10.  I need each point to have a distance to every other point in the array.

Example.  5 points, P1 to P5.  Need distance from P1 to P2, P1 to P3, P1 to P4, etc. etc.

Is there a simple solution that I am overlooking?  NEAR function seems to only let me get one distance and I'm hoping I don't have to create individual shapefiles for each well point.


texts: Will storing geospatial data in a database improve performance of R

texts:
When working with sp objects in R (large SpatialPolygonsDataFrame to be exact) I run into performance issues (mostly memory related - e.g. save_tmap() from tmap package fails).

Will the performance improve if I store the polygons in a database, e.g. Postgres, instead of regular R objects? I have the rgdal and rgeos packages, but no experience with PostGIS.


texts: Reducing images to list and array size mismatch

texts:
I want to convert a google earth image to a numpy array and use Google Earth Engine's Reducer.toList operation. But the size of array I get is different than the size of image as I see in image metadata. Below is the code. I have a sentinel 2 image that I clip it to an area and then look at clipped image dimensions, which will be 11x12. However, reducing it will give me an array of 11x11 elements (the resolution of sentinel band B2 is 10 meters and I mentioned in in reducer parameters). What is the problem with different results?

area = ee.Geometry.Polygon(
    [[[105.532, 19.059], [105.533, 19.059], [105.533, 19.060], [105.532, 19.060], [105.532, 19.059]]])

image = ee.Image("COPERNICUS/S2/20160209T034234_20160209T090731_T48QWG").select('B2')
image_clip = image.clip(area)
info = image_clip.getInfo()
dims = info['bands'][0].get('dimensions')
print(dims)
# Will give us [11, 12] -> 132 pixels

reduced = image.reduceRegion(
    reducer=ee.Reducer.toList(),
    geometry=area,
    maxPixels=1e8,
    scale=10);
data = np.array(reduced.get('B2').getInfo())
print(data.size)
# Will give us a list of 11x11 = 121 pixels


texts:
How to highlight area in KML googlemaps on mouse-over event?

I am trying the same thing as this example

I saw the example given with polylines on the same link above.
But is the same thing possible with KML. I have KML for various regions, But now what i need is to display them in the same way as shown in the wiki-link(wikimapia) or may be in this way - http://livelovely.com

According to what i have found out, the mouse-over event of a KML does not work as of now. So we have to find a way around to this, as in the example of polylines  mouse-over event of tables is used as well if we go through the source code of "lovely.com" they have used mouseover event on div tag, but i don't understand how.

Does any one know how this should be done? Any suggestions, highly appreciated.


texts: How to test if a point shapefile is in a raster in ArcMap using Python?

texts:
I have a shapefile where the features are points. It should be about 500 points. All of this points are on top 5 different rasters. Some points are in one raster, but some others are not.

I need to do batch clipping, but before that, I need to test if a given point XY is in a raster A, otherwise I get an error.

Any ideas how to implement this?

I saw something similar for polygons, but I couldn't find something for points.


texts: Unexpected Polygon Zoom Extents with a Rotated Dataframe

texts:
I have created a series of roated Grid Index Features.  I wrote an ArcPy script to generate the grids and adds a rotation value field.  My plan is to run another ArcPy script that zooms to each grid and exports a PDF.  I am using ArcPy instead of Data Driven Pages because of some additional requirements ArcPy fulfills.

The problem I am experiencing is that I'm getting unexpected zoom results when the dataframe is rotated.  In the example below, I rotate the dataframe 47 degrees to "straighten" the grid.  When I use Zoom to Layer or Zoom to Selected Feature it doesn't provide a tight zoom around the polygon.  If the dataframe is not rotated, the zoom works as expected.  Rotating the dataframe causes the view to zoom out too far.

Interestingly, when using Data Driven Pages and providing the rotation field, the zoom behaves as I expect.  Below are some examples of this phenomenon. How can  I get the  results of third image (the Data Driven method) using ArcPy?

EDIT 1:

The ArcPy code is simply using the Grid Index Features Geoprocessor to build the grids.  Where my code comes into play is when I rotate that geometry.  I write the results of the following function to the geometry of each grid with an Update Cursor.  The code I'm using for calculating the rotation is this:

def RotateAxis(AnchorX,AnchorY,inputx,inputy, rAngle):
    x = inputx - AnchorX
    y = inputy - AnchorY
    resultx = ((x * math.cos(rAngle)) + (y * math.sin(rAngle))) + AnchorX
    resulty = ((x * math.sin(rAngle)) - (y * math.cos(rAngle))) + AnchorY
    return (resultx,resulty)

However, even when I manually rotate the grid in ArcMap with the rotate tool, it provides the same results as the code.

Furthermore, it seems like the actual extent doesn't wrap around these polygons.  When you use the edit tool and click on the geometry you can see the outline of the extent.  I've included another screenshot of this. (I did thicken the outline in the screenshot so it would show up better, its really faint in ArcMap.)  I suspect this extent may have something to do with it, but again, I get the same issue even with manually rotated geometry.

I have not found a way to create geometry that is rotated at the time of its creation.  Grid Index Features only creates unrotated geometry and thus my work around in using Arcpy to rotate after the fact.

Edit 2 :

Just to add further clarification, all of the grids are different sizes and rotations, so defining a set scale and moving the dataframe extent will not work.  The aspect ratio of them however is consistent.

Zoom to Layer - No Dataframe Rotation:

Zoom to Layer - 47 Degree Dataframe Rotation:

Using Data Driven Pages using the Rotation Field:

Rotated Dataframe (317 degrees) with Polygon Extent Visible:


texts: QGIS 2.2 Openlayers Plugin shift between shapefile and base layer for large format

texts: I am attempting to create a 24x36 map of New York City, and overlay it with a shapefile from City Planning. I'm using the Openlayers plugin, and I've tried all of the different types of base street map layers. I have been able to successfully convert the shapefile from a NAD 83 projection (EPSG 2263) to a WGS 84 projection (EPSG 3857), but I'm still suffering from some type of shift whenever I draw the map out on the print composer on a large scale piece of paper. The larger I try to print, the more dramatic the shift the misalignment.


texts: update spatial data from excel to ms-sql server 2008

texts:
I have 1000's of data with 50's of column. It is point data. The data
is in MS-SQL server.Now each day I need to load 100's of data similar
in structure. How can I achieve it ? I tried to append table but it
was lengthy task. I need to change column type since mapinfo does not
properly recognize every field.Also when I change the column it cannot
convert some type.
Is it possible or has anybody done this can of stuff! Is there any solution to this pain?
I was thinking of doing direct update from excel but before I get
started I was wondering advice from experts!!


texts: The My Maps solution from Google is great, especially for quick adding layers and styling. Is it possible to export the map to my own website  as a full sized map and without the My Maps logo? There´s the embed solution but that´s not what I´m looking for. In Fusion Tables this is easy, when complete html code is generated. I imagine that i have to integrate my API key as well.


texts: Large queries in Overpass Turbo ran out of memory

texts:
I would like to query a large bounding box in Openstreetmap using the overpass API from within python. A solution similar to this one would be ideal.

My extent, however is very large, and I would like to determine if there is any native "overpass API" method to loop through "sub-bounding boxes" and extract data for each smaller bounding box.

I can think of ways to do this in python but I would rather take care of this at the Overpass API level. Surely they have devised a way to avoid errors like these:

  runtime error: Query ran out of memory in "query" at line 1. It would
  need at least 531 MB of RAM to continue.

Are there any Overpass-specific capabilities which will successfully download data for large extents? Which commands/queries are of use here?


texts: Equation or approximation of short distance path on Mercator

texts:
I use osm map in my project and I want to draw shortest line between two point like this:

How do they draw this path? Give me please a clue )


texts: carto.com when 2 points have the same location

texts: Carto.com: We are mapping homicides and when two deaths occur at the same location, how do I prevent 2 points with the same lat, lon from overwriting each other?  I want to be able to view each point and the associated data.


texts: Selecting parts of a raster within certain parameters in ArcGIS Pro

texts: I am new to ArcGIS Pro and trying to learn how to do analysis with rasters. I have two raster datasets: one with an elevation profile of a US state, and another showing land use. I am trying to create a new raster that shows an intersection between land use of a specific category and elevation <1000. In ArcGIS Pro, what tool should I use to create the new raster?


texts: How do you create a column computing area in SF with a SpatialPolygonDataFrame?

texts:
I am attempting to add a column that provides the area of my HUC_12 watershed polygons to a watershed data set, but the only code online that demonstrates how to do this is using the raster library, which is incompatible with the sf package and tidyverse. If you have an idea of how to create this is sf, please let me know.
Here is the code I was using for reference:
library(raster)
shape = shapefile("~/Downloads/all_borders.shp")
shape$Area_sqm = area(shape)
raster::area(shape)

Thank you for your assistance.


texts: Webgis with leaflet or ... with another simple solution

texts:
I need to create a static webgis (no connection with db!) to display my interferometric points on a map using Leaflet or similar tools. My points are in geojson but:

How can load these points on a map?
How can I wear my points with a graduated classification and colors?
How can I open their attribute tables with a popup?

I'm not a programmer but I'm sure that could be a tutorial or a simple way to start it with your help.


texts: ClassNotFound Exception: org.postgresql.Driver Using Geotool: gt-imagemosiac-jdbc-x.x.jar

texts:
I am a newbie to GIS, coming from a strict MIS background. Anyway, I have been having a little challenge for the past couple of days trying to import tiles and pyramids into my database. Using a tutorial I found here. Whenever I run this command on windows (8.1);

java -jar <your_geoserver_install_dir>/webapps/geoserver/WEB-INF/lib/gt-imagemosaic-jdbc-{version}.jar import  -config <your geoserver data dir>/coverages/osm.postgis.xml -spatialTNPrefix tileosm -tileTNPrefix tileosm -dir tiles -ext png

I get this exception

java.lang.ClassNotFoundException: org.postgresql.Drive
        at java.net.URLClassLoader$1.run(Unknown Source)
        at java.net.URLClassLoader$1.run(Unknown Source)
        at java.security.AccessController.doPrivileged(Native Method)
        at java.net.URLClassLoader.findClass(Unknown Source)
        at java.lang.ClassLoader.loadClass(Unknown Source)
        at sun.misc.Launcher$AppClassLoader.loadClass(Unknown Source)
        at java.lang.ClassLoader.loadClass(Unknown Source)
        at java.lang.Class.forName0(Native Method)
        at java.lang.Class.forName(Unknown Source)
        at org.geotools.gce.imagemosaic.jdbc.Import.start(Import.java:366)
        at org.geotools.gce.imagemosaic.jdbc.Toolbox.main(Toolbox.java:46)

I have already copied the postgressql driver jar file to my java runtime, lib/ext directory. which is found here; C:\java\jre7\lib\ext. I have also tried to write a java program to test for the driver and it works fine, and by the way I also have Netbeans(8.0.2) installed and the driver is working fine in there too. But I don't know why I still get that exception.


texts: Pivot sf geometry column such that each row contains a single polygon

texts:
I am trying to calculate the average distance to all habitat patches within a landscape from a site point (sf object). To do this, I have polygonized my habitat raster (dissolving to get patches) and am hoping to use st_distance() to calculate the distance from my points to each habitat patch (with 4 possible habitat types).
When I polygonize my raster and convert to an sf object, the resulting data frame looks like this:
habitats.sf <- rasterToPolygons(habitat.raster, dissolve = TRUE) %>%
  st_as_sf()
view(habitats.sf)

Each layer has a number of geometries corresponding to each patch within the layer.
When I calculate the distance between my points and these habitat polygons, the result looks like this:
st_distance(sites.sf, habitats.sf, by_element = FALSE)

with each matrix row corresponding to a site, and each column corresponding to a habitat type.
Is the st_distance() calculation calculating the distance to the closest polygon of each class, or the mean distance of all of the polygons of a given class? If the former, how can I pivot my habitats.sf data frame such that each row only contains a single polygon geometry rather than a list so that I can calculate the distance to ALL patch geometries within each habitat type?


texts: I will often create many temporary layer files with MakeFeatureLayer_management. I was wondering if there was a means to list the temporary layer files I've created. I've tried setting the workspace to "in_memory" and applying ListDatasets () and ListFiles (), but since these layers don't seem to actually be stored in "in_memory", this doesn't produce my desired results. Is there another location these files are stored that I can access? it would be nice to just be able to list all my created layers and delete them without having to keep track of each, to free up locks and what have you.


texts:
I have a little lack, because when my colleagues use QGIS, sometimes they make a mistake with merge function, but singlepart/multipart isn't solution.

A simple example:
Somebody accidentally merged two polygons

The split tool cannot be use in this situation, because these polygons not neighbouring. I see all vertices in vertex editor dialog, but dont separate polygons.

Opposite if I use ArcGIS sketch properties, which shows separate the polygons.

ArcGIS use part(x) separate lines, which help to the identification of the polygons, hereby separates objects.

I want to similar solution in QGIS, that capable this task but didn't find suitable tool/plugin


texts: Seeking collaborative map editor tool for non-technical users with OpenStreetMap basemap

texts:
I'm looking for a tool to allow non-technical users to create map polygons based off printed maps, and export them into a geographiclal format.

It should be based on OSM base tiles (to avoid licensing issues) and ideally offer open satellite imagery too.

The basic requirements are: easy to sign up for, easy to use, and allow data export in an open format.

I have found UMap, which basically does what I need, but I would ideally like something with satellite tiles (a lot of my data is in the countryside, where OSM is less detailed, and it would be useful to check underlying geographical features).

I don't mind asking my users to install a piece of software, so I could use desktop software if that's easier.