texts: Getting list of regions that lines intersect with using ArcGIS Pro?

texts:
I am using ArcGIS Pro.

I have a bunch of line features with unique ids and a bunch of polygons with unique ids. I want to know how many and which polygons each line interacts with. The only way I know how to do this is to Spatial Join and get the join count. But this wouldnt give me the unique ids of all the polygons the lines intersect with.

I am not a GIS expert, this is something I learned for my job and to me it feels like I want to flip the way database table's x and y axis. I often encounter this problem where I want to have a series of information in a cell. I feel like there must be a way to do this that I have not learned yet.


texts: GeoPDF: How to create a hierachy of layers in which I can select all the sublayers from the parent layer?

texts:
I am trying to create this

Having a look at this example and the PDF specs (optional content configuration dictionary keys), I am quite confused about how to create that vs what I have at the moment

Any suggestion on which dictionary key / setup I should look into?
Note: I manipulate the PDF header programmatically (with PyPDF2) and I am not interested in Acrobat Reader options.


texts: Stacking Landsat 8 bands using semi-classification plugin?

texts:
I am currently using the exercise titled Supervised Land Cover Classification Exercise 1

After I have created the individual band sets using the pre-processing tool, the next step in the exercise (see link above) requires to stack the bands together using the band set icon in SCP. My bands 2-7 were in the band set definition box, (in tandem with the above exercise) but when I click to 'create raster of band set (stack bands)' and run it (in band set icon of SCP plugin), I get a completely pitch black Bstack_raster.tif image. I have tried changing the color style under the properties table but it still remains pitch black. I have tried backtracking to see any mistake but as far as I can see I have followed the correct procedure in the exercise.


texts: Mapbox vector tiles version 2 MVTv2 and Geoserver 2.13.x

texts:
I use WMS on Geoserver to display MVTs in conjunction with mapboxgl.js v0.45.0.
A OSM road layer renders 1:1 on top of the Mapbox base layer, so far so good.

However, I got a warning from the js console:

Vector tile source "geoserver" layer "roads" does not use vector tile spec v2 and therefore may have some rendering errors.

Does anybody know about the protobuf version used in geoserver?
I reviewed the source code of the java tile package used and of the plugin in branch v2.13.x. Seems like MVTv1 is used and triggers the javascript warning. Can someone point me to the right location in the java code? Maybe we can upgrade the MVT package so that MVTv2.1 is used.

Here is the source code of the javascript integration in accordance with a GEOS bug report by another user. That works for me.
see here on GEOS JIRA


texts: Writing ISO 19139 compliant XML metadata from spreadsheet?

texts:
I am looking at a way to migrate the metadata from a spreadsheet into ISO 19139 compliant XML metadata documents for a collection of shapefiles. I can parse the data out of the spreadsheet into strings without issue however I am getting lost on how to take those strings and use them to populate the elements.

This example: http://pymotw.com/2/xml/etree/ElementTree/create.html covers most of what I want to do, however I would need to generate the tree beforehand based on the ISO 19139 schema. PyXB seems to have schema bundles for OpenGIS but I would like some insight on how to connect everything.


texts: Keyboard shortcut to select entire column in attributes table QGIS

texts: How can I select an entire column in the attributes table in QGIS? Ideally I would like to do it using a keyboard shortcut. In a search around the net I cannot find a way to do it.


texts: Add a mosaic to a cached service and update the cache maintaining the existing cache tiles

texts:
Is there any way to add a mosaic to a cached service, and update the cache of this service maintaining the existing cache tiles

I have the following scenario in my work,
I add images (1,2,3) to a mosaic data set(A) then I added  the mosaic (A) to an MXD then I publish and Cache the service using ArcGIS server 10.3.1.
Now I got new Images ( 4,5,6) and added them to mosaic (B) and I need to add them to the same MXD and use publish overwrite the service.

But the question here is there any method to avoid starting the cache form the scratch? Since the mosaic (A) was already cached before, and now I only need the new mosaic (B) to be cached to save the time and space, while both A and B lives in the same service


texts: ArcGIS Geoprocessing Service returning feature class?

texts:
I am trying to send some feature classes from geoprocessing services. I have written Python toolboxes which does some business logic and generates feature classes (both point and polygon for different tools).

I have published these tools as geoprocessing services and want to return these feature classes, but so far I am stuck at points where my tools do all successful processing but are unable to send back the feature classes.

I can generate the JSON output of the result through featureset.JSON and send it across as string but I really don't know how to consume that string to render my map at the JavaScript client side. I tried setting the output parameter type as "feature class", "feature set", "shape file", but the services gives a "invalid return type" error.

Can somebody please help me in this?


texts: I am looking for ideas/libraries that enable snapping when I create a circle with leaflet draw. There is the libaray Leaflet snap, that enables  snapping of draggable markers to polylines and other layers. But I am looking for an option to snap to markers when creating a circle. So basically create a circle around the marker.


texts: Why will QGIS only import the first 20 lines of my delimited text file?

texts:
I have a file of bird species and counts with Lat Long that I want to put into QGIS (2.0). However, when I add the file it only shows the first 20 entries (excluding the field titles). This is the same whether I type the file out by hand from scratch, save it in different formats, export a working attribute table from ArcMap etc... Its even done the same thing with a separate CSV file...

Anyone got any ideas?

Thanks in advance!

Alex

This is a sample of the data:

S_ID,S_NAME,COUNT,DATE,OBSERVER,LOCATION,LATITUDE,LONGITUDE
1,Mallard,6,27-Apr-14,Douglas_Philp,"Ely_Trail,_Pontyclun,_Rhondda,_Cynon,_Taff_County",-3.408753487,51.54959248
2,Grey_Heron,1,27-Apr-14,Douglas_Philp,"Ely_Trail,_Pontyclun,_Rhondda,_Cynon,_Taff_County",-3.408753487,51.54959248
3,White-throated_Dipper,1,27-Apr-14,Douglas_Philp,"Ely_Trail,_Pontyclun,_Rhondda,_Cynon,_Taff_County",-3.408753487,51.54959248
4,Grey_Wagtail,1,27-Apr-14,Douglas_Philp,"Ely_Trail,_Pontyclun,_Rhondda,_Cynon,_Taff_County",-3.408753487,51.54959248
5,Canada_Goose,5,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
6,Ring-necked_Pheasant,2,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
7,Manx_Shearwater,1,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
8,Northern_Gannet,5,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
9,Razorbill,12,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
10,Red-billed_Chough,3,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
11,Willow_Warbler,10,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
12,Northern_Wheatear,12,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
13,Reed_Bunting,2,23-Apr-14,Harry_Hops,"Skomer_Island,_Pembrokeshire_County",-5.30310871,51.73776583
14,Barn_Swallow,2,14-Apr-14,Paul_Owen,"WhiteHill,_Pembrokeshire_County",-4.819516968,51.71073247
15,Common_Chiffchaff,2,14-Apr-14,Paul_Owen,"WhiteHill,_Pembrokeshire_County",-4.819516968,51.71073247
16,European_Goldfinch,4,14-Apr-14,Paul_Owen,"WhiteHill,_Pembrokeshire_County",-4.819516968,51.71073247
17,Common_Merganser,4,12-Apr-14,Howard_Williams,"Bute_Park,_Cardiff,_Cardiff_County",-3.185524062,51.48431536
18,Great_Crested_Grebe,2,12-Apr-14,Howard_Williams,"Cardiff_Bay_Harbour,_Cardiff_County",-3.168233956,51.45798394
19,Barn_Owl,1,07-Apr-14,Toby_Rowland,"Carmel,_Carmarthenshire_County",-4.045078962,51.82685504
20,Carrion_Crow,1,05-Apr-14,Toby_Rowland,"Aberaeron_harbour,_Ceredigion_County",-4.263096052,52.24298592
21,Bar-headed_Goose,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
22,Gadwall,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
23,Eurasian_Wigeon,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
24,Northern_Pintail,2,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
25,Common_Buzzard,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
26,Little_Ringed_Plover,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
27,Spotted_Redshank,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
28,Common_Greenshank,1,05-Apr-14,Simon_Priestnall,"Newport_Wetlands_RSPB_Reserve,_Newport_County",-2.963961861,51.5462345
29,Common_Snipe,1,05-Apr-14,Simon_Priestnall,"Kenfig_NNR,_Bridgend_County",-3.731128585,51.51604494
30,Great_Spotted_Woodpecker,1,05-Apr-14,Simon_Priestnall,"Fendrod_Lake,_Swansea_County",-3.903193858,51.66903139


texts: How to tranform a column field into a list using Field Calculator

texts:
I'm trying to make a simple field calculation with ArcGIS 9.3:

New field = Old field - Old field(first value)

which in Calculate Field with python code should be

Expression:

!Old field!  - first

Code Block:

list= [AngOriz]
first = list[0]

The error is

  ERROR 000539: Runtime error : name
  'AngOriz' is not defined Failed to execute (Calculate Field).

How could I transform a column field into a list?


texts: Getting QGIS / Processing revision number from PyQGIS?

texts:
In QGIS, certain processing algorithms change between revision (in my case the merge function between version 2.14.1 and 2.14.7).

How can I access this information from PyQGIS, in order to make a plugin work on both versions?


texts:
I am working on an rainfall analysis and I have to interpolate some datas with kriging method.

Is co-kriging suitable for the rainfall analysis?


texts: How to append a FeatureCollection to a hosted Feature Layer in ArcGIS Enterprise?

texts:
This question pertains to arcgis.gis in Python, for working with hosted Feature Layers in ArcGIS Enterprise.
I'm using the enrich_data function to calculate a buffer polygon from my input (hosted) point layer, then add population statistics to the resulting polygon - this is working well. I'm not specifying an output_name, so the result is an in-memory FeatureCollection.
Now I wish to append the resulting polygon to an existing polygon layer which is hosted in ArcGIS Enterprise, but I can't figure out how to send the Feature Collection to the existing layer.
I tried the FeatureLayer edit_features command, which has an adds option which expects as input:
adds: FeatureSet/List. The array of features to be added

I can obtain the FeatureSet from FeatureCollection.layer.featureSet but this fails in the edit_features.add command with the below error message, which is unclear to me:

The section self._mapping[key] hints that there is an issue with the field mapping between the input and output polygons, but I can't see any options to remap the field names.
TL;DR - how can I append an in-memory FeatureCollection of polygons to a hosted Feature Layer in ArcGIS Enterprise?


texts: Feature selection in multiple shapefile layers

texts:
I've followed the geotools selection lab tutorial. It works perfectly. Now I would like to add another layer, using another shapefile. So I used JFileDataStoreChooser to select another shapefile and added that layer to the existing map. Both layers are now shown on the map frame.

But when I click on a feature, only one type of feature is selected. For example, I have loaded a nodes layer and a lines layer. When I click on a node, it's selected but when I click on a line it's not selected.

I've seen that SimpleFeatureIterator uses featureSource. What should be the approach I follow in this kind of a scenario?


texts: QGIS Project: where to put CORE_Export definition in Visual studio 2010

texts:
I am trying to develop a QGIS application with the following setup:
OSGE4W qgis Libs
QT 4.8.5
Windows 7
Visual Studio 2010

I have all the references setup in VS, but am not sure where define the following lines:

CORE_EXPORT     __declspec(dllimport);
GUI_EXPORT      __declspec(dllimport);

I cannot compile my code unless these are defined.

I do not want to use QTCreator as I don't like the IDE.


texts: Adding dropdown menu for QGIS Plugin to the toolbar

texts:
I have set some actions for my plugin and now I want to group all of them in a dropdown menu as one of the image, which is from another plugin. I've been struggling but I can't find a way to do it like that.

I have tried this:
    def initGui(self):
        """Create the menu entries and toolbar icons inside the QGIS GUI."""

        icon_path = ':/plugins/udt_plugin/icon.png'
        self.action1 = QAction(QIcon(icon_path), u"Action 1", self.iface.mainWindow())
        self.action2 = QAction(QIcon(icon_path), u"Action 2", self.iface.mainWindow())
        self.action3 = QAction(QIcon(icon_path), u"Action 3", self.iface.mainWindow())

        self.popupMenu = QMenu(self.iface.mainWindow())
        self.popupMenu.addAction(self.action1)
        self.popupMenu.addAction(self.action2)
        self.popupMenu.addAction(self.action3)

        self.toolButton = QToolButton()

        self.toolButton.setMenu(self.popupMenu)
        self.toolButton.setDefaultAction(self.action1)
        self.toolButton.setPopupMode(QToolButton.InstantPopup)

        self.iface.addToolBarWidget(self.toolButton)

But it creates a more simple dropdown menu like below, which is not exactly what I want. I'm using QGIS 3.10.


texts: ArcMap 10.3.1 Editor Message: "A layer for PROJECT_LINES was not found..."

texts:
When editing a feature class stored in Oracle Spatial, I am receiving 3 warning messages upon stopping the edit session:

There are 2 more after pressing OK - PROJECT_POINTS and PROJECT_POLYGONS.

I wanted to see if it would happen in a File Geodatabase and a Shapefile, so I saved copies of the original data in those formats and tried to start then stop editing. The FGDB throws up the same messages as above while the SHP does not.

Has anyone seen these messages before and can explain them to me or direct me to more information about them. I'd like to resolve the issue entirely but the data seems to be holding onto some reference even after saving as a new file.

Update 1:  I've contacted esri support. If I find out an answer through them, I'll share it here but they seem baffled by this as well.


texts:
I would like to do an unattended update of my QGIS install every night using windows task scheduler.

Is it possible to automatically update all installed packages with the OSGeo4W-Installer without user interaction?

I know that there are options like -k --autoaccept  and -q --quiet-mode and  -A --advanced but do I have to list all packages i want to upgrade or is it possible to get the same result as clicking through the advanced OSGeo4W-Setup where the installer automatically lists all updateble packages?


texts:
3/26/2014 - Since writing this question, a patch has been released.
It can be found here:  http://support.esri.com/en/downloads/patches-servicepacks/view/productid/160/metaid/2003

The bing basemap seems to have been recently removed from the "Add Basemap" dialogue in ArcMap.  I'm using 10.1.  I've also heard reports from a coworker that a bing basemap layer was replaced by the dreaded red exclamation mark in an existing document she was using.
Any idea if bing/ esri  altered their relationship, or is this just a local problem on my end?
Thank you!


texts: Finding xy coordinates of high and low elevation values in DEM

texts:
I have a 32-bit floating point GRID containing elevation data created in Hazus 3.2. The area defines the elevation of a specific boundary for a flood study. I need to find the xy data for the high and low elevations within the DEM. I basically need to do the opposite of the Get Cell Value tool

I can't figure out a good way to do it. I do have the Spatial Analyst extension at my disposal as well.


texts:
When i run a geographical weighted regression in QGIS, this command pops up 'NoneType' object has no attribute 'crs' See log for more details
and no regression analysis will happen. Does anyone know what this means?


texts: No output from GRASS 7 r.reclass in Processing Python script

texts:
I want to run the r.reclass tool in python on a raster layer to reclassify 1 value (and leave the rest the way they are).

Here is my best attempt so far, after trying to reference documentation and other examples. Here is my code:

import processing

input_raster = QgsRasterLayer('/Users/ep9k/Desktop/Key-LogEcovaluator/Rasters/AestheticMax.tif', 'raster')
output_raster = '/Users/ep9k/Desktop/OutputReclass.tif'
rules_file = '/Users/ep9k/Desktop/reclassrules.txt'

parameters = {'Input' : input_raster,
            'rules' : rules_file,
            'Output' : output_raster,
            'GRASS_REGION_PARAMETER' : '1277290, 1314314, -14004, 11845',   #this is the exent of my raster layer
            'GRASS_REGION_CELLSIZE_PARAMETER' : 0
            }

processing.runAndLoadResults('grass7:r.reclass', parameters)

I run the script and I don't get any errors, but nothing happens. No output, nothing is reclassified. I think I am missing something obvious here.

Here is how I define my reclassification rules. As far as I can tell from the grass 7 documentation this is correct:

35 = 99

Really, I am reclassifying one value (35 = 99).

I get the help docs by running the following in the console:

processing.algorithmHelp("grass7:r.reclass")

Two of the input parameters are confusing to me. Specifically GRASS_RASTER_FORMAT_OPT and GRASS_RASTER_FORMAT_META. Are these required?


texts: How to reopen hidden panels in print composer?

texts:
I unfortunately deleted the menu (Qgis 1.8.0) in print composer (settings for labels etc) . Now I am wondering how to display it again? The menu is gone in all my projects and I could not find a solution for this.
Any suggestions?
Thank You!!


texts:
I created TIN from contour lines, then I convert this TIN to Raster.
After that I clipped the raster by Extract by Mask tool,
now the problem is raster does not fill whole boundary.

Below image show the problem, how can I figure out this?


texts: Identifying zip codes that cover multiple townships

texts: I have shapefiles with zip code tabulations in Massachusetts and township areas. Some zip codes fall completely within a single township, while some zip codes cross multiple town borders. I want to identify those multi-township zip codes and shade those zip areas differently than the other zips. What I believe I should be doing is joining by location zip polygons that overlap with town polygons. However, the problem with this is that the polygon borders for matching town/zip borders don't align perfectly in the shapefiles, so that every single zip polygon technically "overlaps" every neighboring township polygon, so I can't use this to identify multi-town zips. Does anyone have any advise on the best way to proceed here?


texts: How to set CRS for a DXF file and prevent any warning on QGIS using pyQGIS?

texts:
I am trying to import dxf file with no warnings or popup from Spatial Reference Selector. However, even after modifying QGIS Spatial Reference Selector setting, I still get a qgisMessagebar warning.
This is what I have done so far.

    dxf_file = 'D:/mydxf.dxf'
    dxf_file_url = dxf_file + '|layername=entities|geometrytype=Point'
    #setting to suppress the popup
    s = QSettings()
    oldValidation = str(s.value( "/Projections/defaultBehaviour", "useGlobal" ))
    s.setValue( "/Projections/defaultBehaviour", "useGlobal" )
    layer = QgsVectorLayer(dxf_file_url, 'layer_name', 'ogr')
    s.setValue( "/Projections/defaultBehaviour", oldValidation )

So what I want to achieve is, either insert projection to the dxf file before it goes into QgsVectorLayer or to hide the QgsMessagebar.


texts: Manually editing y-coordinate in ArcGIS Desktop?

texts:
I am trying to export dwg files in ArcGIS to a shapefile. My problem is that the y-coordinate is 500000 although it should be 3500000 because the dwg file is made to a coordinate system called Finnish Zone 3. So it is missing the zone number 3 which should be in the false easting value in that coordinate system.

How can I manually put the number 3 in the y-coordinate so that I can then easily project the file to EUREF-FIN and use that with other layers that are also in EUREF-FIN?


texts: Intersecting MultiLinestring based geodataframes with geopandas?

texts:
I would like to perform the following operation using geopandas.

The segments are delimited by red points and the blue items are attribute information. My inputs are the first and second line segments and my output is the third line segment.

Initially I thought this would be an intersection operation, but I soon learned that geopandas can only intersect polygons, therefore something like:

intersection = geopandas.overlay(split_lines, original_lines, how='intersection')

returns the following error:

raise TypeError("overlay only takes GeoDataFrames with (multi)polygon "
TypeError: overlay only takes GeoDataFrames with (multi)polygon

This to me looks like a standard geoprocessing operation and I really hope I won't have to code this up from scratch. Are there any simplified ways to come up with the following result without having to code a custom function?

EDIT
Unfortunately I cannot get this to work for more complicated geometries such as

ORIGINAL_LINES

                                            geometry property
0  (LINESTRING (0 0, 6.656423206909781 4.43757029...        a
1  (LINESTRING (6.656423206909781 4.4375702913320...        b
2  (LINESTRING (8.070636769282876 5.8517838537051...        c
3  (LINESTRING (6.656423206909781 4.4375702913320...        d
4  (LINESTRING (10.98655022583197 1.9375702913320...        e
5  (LINESTRING (13.68293236472948 0.6224568509648...        a
6  (LINESTRING (17.54663566988575 -0.412819329445...        a

SPLIT_LINES

                                            geometry  susc
0  LINESTRING (0 0, 4.160264504318614 2.773481432...     1
1  LINESTRING (4.160264504318614 2.77348143208253...     2
2  LINESTRING (6.656423206909781 4.43757029133205...     3
3  LINESTRING (9.950815132334437 8.18950268263064...     4
4  LINESTRING (13.08444573742037 12.0857007308397...     5
5  LINESTRING (6.656423206909781 4.43757029133205...     4
6  LINESTRING (10.98655022583197 1.93757029133205...     3
7  LINESTRING (15.61478401730761 0.10481876075978...     2

The output seem to be a 1D line... which is incorrect for this application.


texts: How to calculate mean for multiple raster files and multiple variables using R

texts:
I have a folder with multiple raster files (.asc) which contains the biomass of different fish species (e.g., Anglerfish, Bathydemersal etc) for different years (e.g., 00012, 00024, 00036 etc)
I want to calculate the mean per species using R
So far i have done the following:
files_path <- "~/new/asc/" #path where my files are
all_files <- list.files(files_path,
                            full.names = TRUE,
                            pattern = ".asc$") #take all the ascii files in
files_stack <- stack(all_files) #stack them together

Here comes the problem: mean_ <- calc(files_stack, mean) code makes one mean for the entire files_stack.
I want to have mean per each species. Any tips on how to do it?


texts: GDAL_translate: converting ESRI GRID to Geotiff in batch

texts:
The Challenge:
To batch convert 1000 rasters in ESRI Grid binary to Geotiff using either using a batch file or Python 2.7 without ArcGIS/Arcpy module.

The Tools:
Win7Pro, GDAL 1.9.2 (via OSGeo4W), python 2.6/2.7, 1000 binary ESRI Grid rasters (see image)

File structure:

Single file conversion: (piece of cake)

gdal_translate GRID1 -of GTiff GRID1.tif

Batch Route:

for %i in (*.img) do gdal_translate -of GTiff %i %~SourceFilename.tif

The Problem:
ESRI Grids don't have a file extension and are instead treated as file directories by the operating system. I need to find an alternative way to iterate. My only option it seems is to iterate by directory as (to my knowledge) there is no method (to my knowledge) to read the GRIDs into a variable like using the ArcPy method: rasters = arcpy.listRasters("*", "GRID")

Does anyone have any suggestions about how to modify the (*.img) portion of the code so that the loop will read each GRID directory as a file object so I can iterate using the for loop posted above? Note: I'm also aware that I will need to add the parameter for recursion to the batch file and possibly ignore all directories with the name "info".

I know I could just use arcpy and python to accomplish this task, but I wasn't satisfied with the results on google of how to batch convert ESRI GRIDs with only GDAL and maybe python. I appreciate any and all advice!


texts: Grouping of polygons or overlaying layers in same shapefile? Problems calculating the correct polygon areas

texts:
I am having trouble calculating polygon areas of my shapefile. When opening the shapefile in QGIS it seems like there are multiple polygons overlapping within the same shapefile or that polygon areas are grouped together somehow. Then when calculating area with $area, this gives too large areas. What I want to calculate is the area of all the black smaller polygons seperately and not the red somehow "grouped" polygons.

For clarification, CRS of the shapefile is: EPSG:32630 - WGS 84 / UTM zone 30N - Projected (unit is meters).
I have already tried to fix invalid geometry following this guide: https://www.qgistutorials.com/en/docs/3/handling_invalid_geometries.html . But it did not change the way the polygon areas are calculated.


texts: Converting river mile points into a polyline that closely matches the curvature of the stream

texts:
I have a dataset of river miles for all the rivers in Oregon. I am particularly interested in the estuaries on the west coast of Oregon, as I would like to compare sinuosity of the portions of the estuaries that are within tidal reach, and the portions that are entirely fluvial. In the future I may also look into curvature and asymmetry in the meanders as well.
The dataset consists of points at every mile along a river/estuary. Converting these points to a polyline yields a rough trace of the river, but the distances between the vertices means that there are a significant number of straight segments. This is problematic because the vertices are also points of artificially high curvature when they are around meanders, meaning that any value of curvature for a meander would not be accurate.
Is there an approach to creating a polyline from these points, either through interpolation or otherwise, that follows the meanders of the river/estuary accurately enough and has enough points to calculate curvature and asymmetry?
I'm fairly new to coding, but want to build my skillset and understand how to approach problems like this from a computing perspective.


texts: Export list of coordinates in a points layer to spreadsheet?

texts:
This is similar to my previous question, but I would like to ask it for ArcGIS, instead of QGIS. If I am supposed to somehow merge the questions, please let me know.

With an existing point feature class as a layer in an ArcMap 10 document, how do I obtain a list of the points' x y coordinates (ideally in a csv file)? The attribute table for this layer does not have that data.

What I would like is to export those points to a spreadsheet, with the existing attribute table, but also with columns representing each point's corresponding x, y coordinates. Thanks.


texts: How to listen for the completion of a dragpan event in OpenLayers v6.4.3?

texts:
Is there a way to listen for the completion of a drag/pan event in OpenLayers v6.4.3? The 'moveend' event listener fires whenever the map move ends but is not specific to the map being moved by a drag/pan. If this event listener is used it fires when other actions cause the map to move, e.g., changing the center point of the map programmatically. The 'pointerdrag' event fires with every movement of the pointer drag not just when the drag/pan ends.
What we need is to specifically listen for the end of a drag/pan by the user. Is there an OpenLayers event listener for this purpose?


texts:
I am trying to execute the pythonscript code for the automation of a mapservice. While execution tempdraft.sddraft, and updatedDraft.sddraft files are automatically generated but .SD file has not generated. And I am getting the error :

  Runtime error  Traceback (most recent call last):   File "",
  line 1, in    File "c:\program files
  (x86)\arcgis\desktop10.2\arcpy\arcpy\server.py", line 1293, in
  UploadServiceDefinition     raise e ExecuteError: Failed to execute.
  Parameters are not valid. ERROR 000732: Service Definition: Dataset
  C:\msurya\teritory.sd does not exist or is not supported Failed to
  execute (UploadServiceDefinition).

Can somebody help how to solve this problem for automation of publishing a service?

import arcpy, os, sys
... import xml.dom.minidom as DOM
...
... arcpy.env.overwriteOutput = True
... serviceName = "teritory"
... tempPath = r"C:\msurya"
...
>>> path2MXD = r"C:\msurya\teritory.mxd"
... userName = " "
... passWord = " "
...
>>> SDdraft = os.path.join(tempPath, "tempdraft.sddraft")
... newSDdraft = os.path.join(tempPath, "updatedDraft.sddraft")
... SD = os.path.join(tempPath, serviceName + ".sd")
...
>>> arcpy.SignInToPortal_server(userName, passWord, "https://sce2.maps.arcgis.com/")
<Result 'true'>
>>> mxd = arcpy.mapping.MapDocument(path2MXD)
... arcpy.mapping.CreateMapSDDraft(mxd, SDdraft, serviceName, "MY_HOSTED_SERVICES")
...
>>> doc = DOM.parse(SDdraft)
>>> tagsType = doc.getElementsByTagName('Type')
... for tagType in tagsType:
...     if tagType.parentNode.tagName == 'SVCManifest':
...         if tagType.hasChildNodes():
...             tagType.firstChild.data = "esriServiceDefinitionType_Replacement"
...             tagsState = doc.getElementsByTagName('State')
... for tagState in tagsState:
...     if tagState.parentNode.tagName == 'SVCManifest':
...         if tagState.hasChildNodes():
...             tagState.firstChild.data = "esriSDState_Published"
...
... # Change service type from map service to feature service
... typeNames = doc.getElementsByTagName('TypeName')
... for typeName in typeNames:
...     if typeName.firstChild.data == "MapServer":
...         typeName.firstChild.data = "FeatureServer"
...         configProps = doc.getElementsByTagName('ConfigurationProperties')[0]
... propArray = configProps.firstChild
... propSets = propArray.childNodes
... for propSet in propSets:
...     keyValues = propSet.childNodes
...     for keyValue in keyValues:
...         if keyValue.tagName == 'Key':
...             if keyValue.firstChild.data == "isCached":
...                 keyValue.nextSibling.firstChild.data = "false"
...                 configProps = doc.getElementsByTagName('Info')[0]
... propArray = configProps.firstChild
... propSets = propArray.childNodes
... for propSet in propSets:
...     keyValues = propSet.childNodes
...     for keyValue in keyValues:
...         if keyValue.tagName == 'Key':
...             if keyValue.firstChild.data == "WebCapabilities":
...                 keyValue.nextSibling.firstChild.data = "Query,Create,Update,Delete,Uploads,Editing"
...
>>> f = open(newSDdraft, 'w')
... doc.writexml( f )
... f.close()
... analysis = arcpy.mapping.AnalyzeForSD(newSDdraft)
...
>>> if analysis['errors'] == {}:
...     # Stage the service
...     arcpy.StageService_server(newSDdraft, SD)
...
>>> arcpy.UploadServiceDefinition_server(SD,"My Hosted Services",serviceName,"","","","","OVERRIDE_DEFINITION","SHARE_ONLINE","PUBLIC","SHARE_ORGANIZATION","")

  Runtime error  Traceback (most recent call last):   File "",
  line 1, in    File "c:\program files
  (x86)\arcgis\desktop10.2\arcpy\arcpy\server.py", line 1293, in
  UploadServiceDefinition     raise e ExecuteError: Failed to execute.
  Parameters are not valid. ERROR 000732: Service Definition: Dataset
  C:\msurya\teritory.sd does not exist or is not supported Failed to
  execute (UploadServiceDefinition).


texts: QGIS 64-bit ODBC set up link to 32-bit Access on win7

texts:
I am trying to set up a an ODBC link to Access 2007 32-bit but not having any success.  My system is a win7 pro 64-bit and I have QGIS 64-bit installed.  Here are the steps I have already tried.

Set up an ODBC link using c:\windows\syswow64\odbcad32.exe.  However, when I try to use that, the error message states that there is an architecture mismatch, which I assume means between 64- and 32-bit architecture.
Tried to set up the ODBC link from the standard location in Control Panel, but the error is that the Microsoft Access driver isn't installed.
Downloaded the Microsoft Access 64-bit Runtime, but when installing the error comes up that the 32-bit Access is installed and the 64-bit Runtime cannot be installed.

I am not sure where to go from here.

I know that I can change databases and go with a spatial system, but this is a regular download from a government site, which would mean a lot of extra work.  Any assistance would be much appreciated.


texts: Elevation profile that can follow the change in view of Atlas page in QGIS

texts: I think I have read every thread on QGIS profiles in the stack and I cannot find a way to create a profile that can be produced in context on each page of an atlas print.  I can create an image for each segment and call the image as per the appropriate atlas page but none of them allow any control of the y-axis label so they all restart at 0.  I want the profile on each page to be a continuation of the one that proceeded it.  I have a multi page pipeline alignment and I want the appropriate profile for the shown segment to appear on each page.  I hope this makes sense.


texts: Get all polylines from a shapefile (containing polylines)

texts:
I have a *.shp file which contains only polyline features. How can I iterate over all the point features. Given ESRI Polyline example only uses already defined points to make polylines. But I want to iterate over all the polylines found on a shapefile.

How can I do this using ArcPy?


texts: GeoTiff to 16 bit tiff, png or bmp image for heightmap

texts:
I have a GeoTiff that acquired through processing Mars Hirise stereo images with ISIS3 and ASP. I want to convert it to a 16bit greyscale tiff, png or bmp file to use as a heightmap in a game engine, unreal engine 4. I tried converting it with gdal_translalte, first to 16 bit .raw but turns out photoshop can't open it in 16 bit only in 8bit, I get a message saying that the size is to big for the file in 16 bit, something that seems to have to do with the file size in byte, the two channels and the image size I try to import. Then I tried converting it to tiff and png, also with gdal_translate but both images are just black when I open them. I tried opening my GeoTiff with BelnderGIS but it is just a flat plane. I opened the Geotiff in QGIS and then I see that it does have information and is not just flat and black. I made a smaller cropped version of the DEM and opened it in MicroDEM which works but when I try to open the original file it says that it is too big.

Does anybody know how I could do this, with gdal, qgis or something else? I'm very new to using both gdal and and qgis, I just followed examples from others but there might be a way to do it so it works for me in the end.


texts: Differences between capabilities of QGIS and ArcGIS in 3D modeling

texts:
I am learning to use QGIS and ArcGIS for Desktop.
What are the differences between the capabilities of QGIS and ArcGIS in term of processing & representing 3D model?


texts: Cut Tool or Process to mimic cutting in a "cookie cutter" like fashion

texts:
Say I have a sidewalk polygon, and on that segment of sidewalk I have 5 tree grates that I want to cut out of the sidewalk polygon. Is there a tool or an easy process that I could cut out one tree grate in the sidewalk polygon once and be able to cut out 4 more with same exact dimensions as the first cut? As if I was using a cookie cutter to cut out consistent shapes? Besides just doing it manually, which leads to inconsistent shapes and takes a bit of time, I haven't been able to find anything that will recreate another cut that was done to a polygon.

I am using ArcGIS 10.2.2 for Desktop.


texts:
I've been trying to upload my satellite imagery which is 1.5 GB to my PostgreSQL database,  using the code below, the error for each command is stated below;
raster2pgsql -s 4326 -t auto -F -I -C -Y '.\Ceva I Ra 2020.tif' public.Ceva2020 | psql -d MB_database -U postgres

raster2pgsql -s 4326 -t auto -F -I -l 2,4,8,16,32,64,128,256 -C -Y '.\Ceva I Ra 2020.tif' public.Ceva2020 | psql -d MB_database -U postgres

Program 'psql.exe' failed to run: Insufficient memory to continue the execution of the program. At line:1 char:109
+ ... Ceva I Ra 2020.tif' public.Ceva2020 | psql -d MB_database -U postgres
+                                           ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~.

At line:1 char:1
+ raster2pgsql -s 4326 -t auto -F -I -l 2,4,8,16,32,64,128,256 -C -Y '. ...
+ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    + CategoryInfo          : ResourceUnavailable: (:) [], ApplicationFailedException
    + FullyQualifiedErrorId : NativeCommandFailed


texts: Finding Stream Order from DEM using ArcGIS Desktop?

texts: I have DEM Data for the study area, I want to know how to find out the stream order using DEM. I saw Stream order tool in Hydrology toolset in ArcGIS. in Stream order tool what is input stream raster and how to create it, which method of stream ordering is mostly used


texts: PyQGIS script processing: 'unicode' object has no attribute 'getFeatures'

texts:
I'm beginner.
I try to use a code that works to create a script processing.
The code creates folder and subfolder using the value os a column.
This is then processing script

import os
##Crea folder and subfolder=name
##Vector=group
##layer=vector
##idx=field layer
##path=folder

#layer = iface.activeLayer()
#idx = layer.fieldNameIndex('NAME')
#path = "D:\\GIS\\_Temp\\"

for feat in layer.getFeatures():
    attrs = feat.attributes()
    if not os.path.exists(path + attrs[idx]):
        os.makedirs(path + attrs[idx])

for feat in layer.getFeatures():
    attrs = feat.attributes()
    foto=os.path.join(attrs[idx], attrs[idx]+"_foto")
#    foto= attrs[idx]+'\\'+attrs[idx]+'_foto'
    if not os.path.exists(path + foto):
        os.makedirs(path + foto)
    else:
        continue

for feat in layer.getFeatures():
    attrs = feat.attributes()
    doc= os.path.join(attrs[idx], attrs[idx]+"_documenti")
    if not os.path.exists(path + doc):
        os.makedirs(path + doc)
    else:
        continue

for feat in layer.getFeatures():
    attrs = feat.attributes()
    vid= os.path.join(attrs[idx], attrs[idx]+"_video")
    if not os.path.exists(path + vid):
        os.makedirs(path + vid)
    else:
        continue

but It returns this error: 'unicode' object has no attribute 'getFeatures'  .
and I don't understand what the problem is.


texts: Split a linestring by distance every x meters using PostGIS

texts: I'm looking to split a linestring every x meters. I don't know which way to start, and I'm going in circles! I can find a point at 20 % from the start with ST_LineInterpolatePoint, but not at 20 m. The idea is to determine the set of points at 20, 40, 60 m and use ST_LineSubstring then. Does that seem correct to you?


texts: gdalwarp won't run with "cubic" resampling method

texts:
I am having an issue with running the gdalwarp command from command prompt.
Basically I have a raster WGS84 file, which I would like to reproject to spatial reference which is using the Azimuthal Equidistant projection, but with "cubic" resampling method.
This is the syntax that I am using:

"C:/gdalwin32-1.4.1/bin/gdalwarp.exe" -s_srs EPSG:4326 -t_srs "+proj=aeqd +datum=WGS84 +lat_0=40.81266 +lon_0=14.414252" -r cubic -of GTiff "C:/vesuvius_wgs84.tif" "C:/vesuvius_cubic_aeqd.tif"

(40.81266, 14.414252 is the center of the raster)

For some reason when I press "Enter", the vesuvius_cubic_aeqd.tif file is not created, but instead I get the following text printed on my command prompt:

gdalwarp [--help-general] [--formats]
    [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]
    [-order n | -tps | -rpc | -geoloc] [-et err_threshold]
    [-refine_gcps tolerance [minimum_gcps]]
    [-te xmin ymin xmax ymax] [-te_srs srs_def]
    [-tr xres yres] [-tap] [-ts width height]
    [-ovr level|AUTO|AUTO-n|NONE] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]
    [-srcnodata "value [value...]"] [-dstnodata "value [value...]"] -dstalpha
    [-r resampling_method] [-wm memory_in_mb] [-multi] [-q]
    [-cutline datasource] [-cl layer] [-cwhere expression]
    [-csql statement] [-cblend dist_in_pixels] [-crop_to_cutline]
    [-of format] [-co "NAME=VALUE"]* [-overwrite]
    [-nomd] [-cvmd meta_conflict_value] [-setci] [-oo NAME=VALUE]*
    [-doo NAME=VALUE]*
    srcfile* dstfile

What is the problem?

Interestingly, if I remove the "-r cubic" part like so:

"C:/gdalwin32-1.4.1/bin/gdalwarp.exe" -s_srs EPSG:4326 -t_srs "+proj=aeqd +datum=WGS84 +lat_0=40.81266 +lon_0=14.414252" -of GTiff "C:/vesuvius_wgs84.tif" "C:/vesuvius_cubic_aeqd.tif"

everything works correctly, and vesuvius_cubic_aeqd.tif file is created.
To make it more strange, I tried replacing the "-r cubic" with "-r bilinear" or "-r average", but again neither of these worked. I also tried, removing the "-of GTiff", and that did not work either.

What am I doing wrong in the upper gdalwarp command syntax?


texts: How to calculate size of a pixel for a scanned map?

texts: A land use map has a scale of 1:150000. The map is scanned at a resolution of 150 DPI. What will be the surface area in m square represented by each pixel in the resulting easter map?


texts: Transforming raster to vector in ArcGIS for Desktop?

texts: I want to convert from raster data to vector dataset using ArcGIS. Raster data consists of pixels and center of z values. Vector data sets consist of points, lines and polygons and x,y values. When I convert between data sets, how is the programme converting and describing from z value to x,y values? Shortly, how can the programme describe from raster to feature class?


texts: How to autofill output column in arcmap batch process menu?

texts: Referring to the screenshot, I'd like to know how to fill the output column with unique values that either have the same name as the input files or otherwise have some unique value in each cell. If I fill in the first output cell, and then right-click and use fill, it fills the entire column with the same file name as used in the first cell, therefore I only get one output file. If I don't fill in the output column at all, I get unique file names for each cell, but they are saved to the default gdb which gives me difficulty. What I'd like to do is to save unique file names to a directory that I choose. Is there a way to do this?


texts:
I have a census tiger line file for counties available from the Census website.

I am trying to display this file on a ggmap, but fortify seems to jumble the regions. I can't find anyone else who has run into this error, but I'm certainly not doing anything out of the ordinary...

library(maptools)
library(ggmap)
library(rgeos)

# read shapefile
WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0")
Counties <- readShapeSpatial("tl_2013_us_county.shp", proj4string = WGS84)

# limit to Ohio
Counties <- Counties[ Counties$STATEFP == 39, ]

# fortify
Counties.fort <- fortify(Counties, region = "GEOID")

# get base map and plot
OhioMap <- get_map("Ohio", source = 'google', maptype = 'terrain', zoom=7)
ggmap(OhioMap) + geom_polygon(aes(x=long, y=lat), fill = 'grey', alpha=0.4,
                              color = 'black', data = Counties.fort)


texts: Check if property exists on features in layer

texts:
I was wondering if there is a good way to check to see if a property exists on a particular feature layer in GeoServer? My use case is that I have a bunch of features I want to display on a map. Some of the feature layers will have a date property associated with them. If they have a date property I want to be able to display a control to filter the features shown by date. I noticed that WFS has an operation called GetPropertyValue which will return an exception if the property doesn't exist. However, if it does it (as expected) returns the value for that property for all features which can be very slow. My general algorithm is as follows:

var wfsLayer = new OpenLayers.Layer.Vector('MyVectorLayer', {
    strategies: [new Openlayers.Strategy.BBOX()],
    protocol: new OpenLayers.Protocol.WFS({
        url: 'http://mygeoserver/geoserver/wfs',
        featureType: 'MyVectorLayer',
    })
});

// if layer has date property (this is the part I need help with)
    var dateFilter = new OpenLayers.Filter.Comparison({
        type: OpenLayers.Filter.Comparison.BETWEEN,
        property: 'date',
        lowerBoundary: startDateVal,
        upperBoundary: endDateVal
    });
    wfsLayer.filter = dateFilter;
// end if


texts: How to change the colormap in QGIS 2.01?

texts:
I decided to say goodbye to ArcGIS and switch to QGIS which is nice, but it has given me a hard time to handle various problems.
Few months ago, I made a map using QGIS 1.8 with a relief background having a colormap from green to brown. Later on, I upgraded to QGIS to 2.01. Now, when I reopen the same map project with the intent to change the colormap of the map because that one was very pale. But, it seems that either I forgot how a colormap can be changed or my QGIS has lost that capability or there is a problem with my tif file which has three bands.

This is what I get when I open the Colormap tab in the Layer Properties window:

Any idea of how I can change the colormap?


texts:
I've tried to install Nominatim with a planet.pnf file to a server which is hosted at Dacentec(36gb ram, 16 threads, 1x2TB 7200 RPM hard).
The installation was going well, it was already indexing data at Rank 26, when suddenly the server was remotely shut down. I"ve received the following message from Dacentec

  Your server with IP x.x.x.x has been flooding the network with traffic in excess of 833315000 bits/sec and 109198 packets/sec. This server has been shutdown.

What could have been the problem? I tought that indexing the data shouldn't even require network access, am I wrong in this matter?

As the server was rented Monday, and only the Nominatim installation ran on it, I am totally buffled by this issue.

Since then I've learned that the traffic was outgoing traffic.


texts: Cannot produce x,y coordinates in ArcPy?

texts:
ArcMap version 10.3.1
Python version 2.7.12

I'm trying to produce a points table with one row. The row should include x,y coordinates for the point and optionally a name.

I have tried this two different ways and failed, each yielding separate results.

Attempt 1:

(Creates a shapefile, but it is empty and unprojected)

temp is a variable holding '[37.00, -79.00]'

pointslist = []
pointslist.append(temp)

pt = arcpy.Point()
ptGeoms = []
for p in pointslist:
    pt.X = p[0]
    pt.Y = p[1]
    ptGeoms.append(arcpy.PointGeometry(pt))

arcpy.CopyFeatures_management(ptGeoms, r"C:\Temp\test.shp")

Attempt 2:

(Always results in the name field being populated, but the x,y remaining null)

fc = "p:/435final/standin.gdb/userx"
cursor = arcpy.da.InsertCursor(fc, ["SHAPE@XY"])
xy = (39.55, -79.00)
cursor.insertRow(["USER", xy])
del cursor


texts: Plant area estimation and phenological map creation

texts:
I am thinking to estimate the area of sugarcane plants. For that, firstly, I have to extract only the sugarcane area from the image. I shall use NDVI/EVI/Chlorophyll Index to delineate the same but tye problem for me lies in the delineation of each growth stage from a single image, instead of time-series.
Then, if I try to find all the growth stages of sugarcane in one image, at a time and try to estimate the area for each life cycle, how to execute the same?  Signature separability and building spectral library is a way to classify but I don't know the different spectral profile for each growth stage of sugarcane.
Germination, Tillering, Maturation and grand growth- all have to be classified from a single imagery (Sentinel-2) and visualise it on a map.


texts: Creating one buffer in ArcMap for different distances on same polygon?

texts: I am using ArcMap 10.2 and I have a polygon feature that has four sides. On the south side, I would like a buffer of 500 feet and west side, I would like a buffer of 300 ft. I do not need a buffer of the remaining sides (I realize I can just delete those vertices later).


texts: Registering SQL Spatial View in ArcSDE and updating Row_ID?

texts:
I am having issues registering an inner join SQL view with a spatial column in ArcSDE. Any help would be much appreciated.

Here is the background info:

ArcSDE 10, SQL Server Enterprise 2008 R2, SQL 'Geometry' Spatial Type, Non-Versioned database, DBO schema, My Windows login is a sysadmin user

Here is the goal:

The view I am trying to register is an inner join of a polygon feature class to a table. The left of the join is the table, which in the database has a ONE_TO_MANY relationship with the feature class, where the feature class is the primary key as ONE and the table is the foreign key as MANY. I need to display in ArcMap the polygon geometries of all the feature class records as joined to the table. Since this is a ONE_TO_MANY relationship, it will result in coincident/duplicate polygon geometries, but that's OK.

Initially I am adding the view to SDE using sdetable -o create_view:

sdetable -o create_view -T Int_Veg_TEST -t "Interaction,Vegetation" -c
"Interaction.ObjectID AS INT_OID, Interaction.INT_UUID, Interaction.EVENT_UUID,
Vegetation.ObjectID AS VEG_OID, Vegetation.VEG_UUID, Vegetation.Shape" -w
"Interaction.EVENT_UUID = Vegetation.VEG_UUID" -i -D

This finishes OK, but the resulting view only displays the count and geometries from the right of the join (feature class) and I need it to display the feature class geometry and count as joined to the left (table). I would love to set the Row_ID when creating the view but that option is not available in -o create_view. When using the sdetable -o describe option, it shows that ArcSDE has pulled the feature class ObjectID field as the Row_ID:

Table Int_Veg_TEST:
Column name    Attribute    type    Null?  Length,DPs   RowID Column?
-------------------------------------------------------------------------------
INT_OID                  SE_INT32  NOT NULL   10
INT_UUID                 SE_UUID   NULL       38
EVENT_UUID               SE_UUID   NULL       38
VEG_OID                  SE_INT32  NOT NULL   10
VEG_UUID                 SE_UUID   NULL       38
SHAPE                    SE_SHAPE  NULL        0
OBJECTID                 SE_INT32  NOT NULL   10          SDE Set

From here I am having trouble using command line sdetable or sdelayer to update the Row_ID of an existing spatial view. Anyone know if this is possible using command line?

I can update the Row_ID manually in the sde_table_registry rowid_column, and it works. All polygon geometries for the join are displaying in ArcMap.

2 problems here though:

A: SDE still thinks it is maintaining the Row_ID. And I can't find where in the sde registry to manually update this to USER.

Table Int_Veg_TEST:
Column name   Attribute   type Null?   Length,DPs   RowID Column?

INT_OID       SE_INT32     NOT NULL        10          SDE Set
INT_UUID      SE_UUID      NULL            38
EVENT_UUID    SE_UUID      NULL            38
VEG_UUID      SE_UUID      NULL            38
VEG_OID       SE_INT32     NOT NULL        10
SHAPE         SE_SHAPE     NULL            0

B. I'd rather not update the SQL tables directly.

I've also tried creating a new view using SQL Server Management Studio and registering it in SDE using sdelayer -o register, but this also does not work:

SQL:

SET QUOTED_IDENTIFIER ON
GO

CREATE VIEW [Int_Veg_TEST]

WITH SCHEMABINDING

AS

SELECT dbo.INTERACTION.OBJECTID AS INT_OID, dbo.INTERACTION.INT_UUID,
dbo.INTERACTION.EVENT_UUID, dbo.VEGETATION.OBJECTID AS VEG_OID,
dbo.VEGETATION.VEG_UUID, dbo.VEGETATION.SHAPE

FROM dbo.INTERACTION
INNER JOIN dbo.VEGETATION ON dbo.INTERACTION.EVENT_UUID = dbo.VEGETATION.VEG_UUID

GO

CREATE UNIQUE CLUSTERED INDEX Ix_INT_OID_UNIQ_CLUSTERED
ON Int_Veg_TEST (INT_OID)

GO

command line:

sdelayer -o register -l Int_Veg_TEST,Shape -e a -t Geometry -C INT_OID,USER -i -D

Error: Underlying DBMS error (-51).
Error: Cannot Create Layer.
Microsoft SQL Server Native Client 10.0: Cannot alter 'INT_VEG_TEST' because it is not
a table.

More information is available.


texts: Can't install QGIS on Linux / "E: broken packages"

texts:
I am trying to install QGIS the first time on Linux Mint 20.1 Cinnamon by following the official QGIS-Installation instructions for Debian/Ubuntu.
(Please note, that I am Linux beginner and wouldn't be asking here if I had found any proper solution in the internet myself.)
I got stuck here:
    apt install qgis qgis-plugin-grass
Reading package lists... Done
Building dependency tree
Reading state information... Done
Some packages could not be installed. This may mean that you have
requested an impossible situation or if you are using the unstable
distribution that some required packages have not yet been created
or been moved out of Incoming.
The following information may help to resolve the situation:

The following packages have unmet dependencies:
 qgis : Depends: libgdal20 (>= 1.8.0) but it is not installable
        Depends: libqgis-3d3.16.3 but it is not going to be installed
        Depends: libqgis-analysis3.16.3 but it is not going to be installed
        Depends: libqgis-app3.16.3 but it is not going to be installed
        Depends: libqgis-core3.16.3 but it is not going to be installed
        Depends: libqgis-gui3.16.3 but it is not going to be installed
        Depends: python3-qgis (= 1:3.16.3+28bionic) but it is not going to be installed
        Depends: qgis-providers (= 1:3.16.3+28bionic) but it is not going to be installed
 qgis-plugin-grass : Depends: qgis-provider-grass (= 1:3.16.3+28bionic) but it is not going to be installed
                     Depends: grass740 but it is not installable
                     Depends: libgdal20 (>= 1.8.0) but it is not installable
                     Depends: libqgis-app3.16.3 but it is not going to be installed
                     Depends: libqgis-core3.16.3 but it is not going to be installed
                     Depends: libqgis-gui3.16.3 but it is not going to be installed
                     Depends: libqgisgrass7-3.16.3 but it is not going to be installed
**E: Unable to correct problems, you have held broken packages.**


texts:
I have ArcView 10.1 and am, without luck, searching for/how to "dock the TOC" on the left hand side of the screen.

When I grab the TOC, I see the buttons light up with arrows for up/down, left/right etc.

How do I dock the TOC?


texts: Cesium - Camera.computeViewRectangle to get current view bounds

texts:
I'm trying to get custom west-south-east-north bounds in Cesium's HelloWorld.html example, using Camera.computeViewRectangle, as suggested in this answer and its elaboration.

Perhaps I misunderstand the purpose of the function, but which input parameter should I provide, if I'm trying to obtain the rectangle bounds of the view currently displayed?

As illustration here is the adapted HelloWorld.html example embedded in a custom web application:

cesiumContainer


texts:
I need to create a model. I have a group of regions (A, B, C)
each one of those has multiple rasters and some other shapefiles. For every region is needed to calculate the area of every region content in the attribute table of each shape, then to export every table to a .xls file. Also, I have to display certain symbology, add the legend and export it as a png image.
So, I have region "A, B, C"; each of those has raster "1, 2, 3" and shapefiles   "I, II, III".In the end, there has to be a .png image with its own legend and a .xls table of every 1,2,3, and I, II, III.

I started with shapefiles because i found it easier, but i dont know how to get the .png files.

Ignore the names within the figures.


texts: nx_spatial does not read all shapefile features

texts:
I have a polyline shapefile for a road network. The roads are split at every intersection (node), e.g. if two roads intersect at one point, then they are split into four links(edges). So using the attribute table of the shapefile, I know this shapefile has 600 edges. But when the shapefile is read in nx_spatial, total edges found is only 497. I would like to know why it is missing other 103 road links (edges). Thanks

import nx_spatial as ns
roads = ns.read_shp(path)   #path = path to polyline shapefile
roads.number_of_edges()     #output = 497

Hope these info helps: Windows7, Python 2.7 Anaconda 64bit, Eclipse. The problem is related to nx_spatial. But I could not create such tag for low reputation.


texts: Changing deselected highlight color in Attribute Editor window

texts:
Does anyone know if you can change the highlighted color or deselected color in the Attribute Editor window? My eyes are about to bleed from the strain. Here's what I'm talking about:

This light gray color after I click off of it is nearly invisible to my eyes. I'm running ArcMap 10.2.1


texts: What XY coordinate system goes with EGM96 when projecting elevation raster in ArcGIS Pro?

texts:
I am using ArcGIS Pro 2.3.1.

I have an elevation raster in WGS84 and a national Z coordinate system. I need to project this DEM raster to the EGM96 Z coordinate system.

Which XY coordinate system goes with that?


texts: How to disable 'Use JAI ImageRead (deferred loading)' on an Image Mosaic layer?

texts:
I'm running GeoServer from the binary files on Ubuntu 19.10. I'm trying to optimise Image Mosaic layers for speedy display via WMS. To do this I would like to enable Multithreaded granule loading (disable JAI ImageRead to use it), however this requires Use JAI ImageRead to be set to off. The tick box beside this option on the layer settings page is greyed out and cannot be unchecked. How can I set this option to off?

Locks   0
Free locks
Connections     0
Memory Usage    343 MB / 672 MB
Free memory
JVM Version     Private Build: 1.8.0_252 (OpenJDK 64-Bit Server VM)
Java Rendering Engine   sun.java2d.pisces.PiscesRenderingEngine
Available Fonts     GeoServer can access 314 different fonts. Full list of available fonts
Native JAI  true
Native JAI ImageIO  true
JAI Maximum Memory  336 MB
JAI Memory Usage    5 MB
Free memory
JAI Memory Threshold    75%
Number of JAI Tile Threads  7
JAI Tile Thread Priority    5
ThreadPoolExecutor Core Pool Size   5
ThreadPoolExecutor Max Pool Size    10
ThreadPoolExecutor Keep Alive Time (ms)     30000
Update Sequence     588```


texts: Using gpkg layer as input in rasterize with PyQGIS?

texts:
I try to rasterize a gpkg layer but it does not work efficiently. It always return empty raster data. However, if I run it by using processing toolbox it works. To see the code, I look the processing history and it gives me that:

processing.run("gdal:rasterize", {'INPUT':'D:/6/6.gpkg
layername=Nehir_Hidrojeolojik','FIELD':'fid','BURN':0,'UNITS':1,'WIDTH':524.3125958528195,'HEIGHT':904.5392875018064,'EXTENT':'590591.444106,601077.696023,4725233.720396,4743324.506146','NODATA':0,'OPTIONS':'','DATA_TYPE':5,'INIT':0,'INVERT':False,'OUTPUT':'D:/deneme shp/2.tif'})

However, as you see, after 'INPUT':'D:/6/6.gpkg it passes the new line and begin with layername=Nehir_Hidrojeolojik'. As you guess, it does not work in python console when I use it directly.

So I try to add | between 'INPUT':'D:/6/6.gpkgand 'layername=Nehir_Hidrojeolojik' something like that:

processing.run("gdal:rasterize", {'INPUT':'D:/6/6.gpkg|layername=Nehir_Hidrojeolojik','FIELD':'fid','BURN':0,'UNITS':1,'WIDTH':524.3125958528195,'HEIGHT':904.5392875018064,'EXTENT':'590591.444106,601077.696023,4725233.720396,4743324.506146','NODATA':0,'OPTIONS':'','DATA_TYPE':5,'INIT':0,'INVERT':False,'OUTPUT':'D:/deneme shp/2.tif'})

But it does not work again. Is there anybody who have an idea about that?


texts: How to rotate polygon to find best fit in Mapnik?

texts:
I have a quite simple problem but somehow I just fail to solve it.

There is a polygon and a range xsize and ysize. The x and y size determines the mapnik generated image.

The polygon should now be made to fit best, by deciding if:

it is good already or
a x and y swap and later rotation by 90 degress would make it fit best.

Boundary box and Envelope is given, so no problems there.

The problem is realy just the rotation criterion which I'm lost in.

I hope you understand the simple problem and can show me an idea how to solve it.

Thanks!


texts: Sending parameters to ArcGIS script through web interface

texts:
I have written a simple script in ArcGIS environment and I want to run the script through a PHP program.

At present, I am trying this within localhost. In my script, I have defined two parameters. My script is mentioned below.

import arcpy
from arcpy import env
import os
import sys
print("Hellow")
E=sys.argv[1]
E1=sys.argv[2]
arcpy.AddMessage(E+"It is running")
print(E1)

My expectation is to send 2 parameters to this arc script through a  web interface.
In my PHP programme, I have called the following command to execute my script.

$output = shell_exec("C:\wamp\www\ssh\php_arc\1\test.py '". $_GET['EmployeeID']."' '".$meek1."'");
echo $output;

This is working. But I am just getting the "Hellow"^^ as a result. I  am not getting the other results (data given as parameters).

Can anyone help me to make this successful?


texts: Runtime error GDALWarp command in C++ programme

texts:
I am trying to use GDALWarp in VS2019: A snapshot of the code is
#include <gdal_utils.h>
#include <gda_priv.h>
#include <cpl_string.h>

int main()  { GDALDataSetH* Dataset=GDALOpen(); char ** argv= NULL;
argv = CSLAddString(argv,"-r"); argv = CSLAddString(argv,
"cubicspline"); //I have added more options ...

GDALWarp("output.tif",NULL,1,Dataset,opt,NULL);

}

I am getting an error on the GDALWarp command, the error states that "exception thrown at (gdal301.dll) in GDALWarp". The programme builds correctly.
Can you help me fix this?


texts: Adding point at the beginning of call out line

texts:
I am trying to add a point at the start of my callout lines in QGIS but I don't manage.
I create an additional geometry generator to my callout line and use this code:
make_point(closest_point($geometry, make_point("X_Core", "Y_Core")))

What am I doing wrong?


texts: arcpy.mapping exportPDF with ONLY one buffer ring turned on for associated iterated feature

texts:
I am exporting X number of maps iterating through a point feature class. Each feature class has a buffer drawn around them. I know how to export a PDF for each feature class at a fixed scale (if I lock extent the map just pans from one feature to another and works; however, if I don't lock the extent and write zoomtofeature code then the map never refreshes and I get the same extent for every iterated feature). I can get by that for now although it is frustrating.

How can I export a PDF for each feature but with only the associated buffer ring turned on and all others turned off? (many buffer rings overlap each other and I've been asked to produce just one ring per feature regardless of overlapping buffers from nearby features).


texts:
I have access to some open data for which I want to obtain some coordinates to put onto a map. Unfortunately, the data is in the Shapefile format which is rather difficult to work with. I am trying to convert it to something easier, preferably geoJSON or some kind of XML-based document.

A requirement is to do this within the client's application in the event that the data changes, so unfortunately I am unable to use some online tool nor can I just open the data in QGIS and export the data in another format.

I attempted to write my own Shapefile parser however I ran into difficulty with IndexOutOfBoundsException while reading some of the records despite having shapes of type 1 (Point). Furthermore, there is a projection file that I have to deal with and I haven't been able to find much on how to map the coordinates according to a PRJ.

Is there any open-source library I can include in my application that can take a Shapefile collection (SHP, SHX, DBF, PRJ, SHP.XML) and convert it to geoJSON? I'm particularly looking into Java for Android/BlackBerry development, possibly Windows Phone in the far future.


texts: Search feature property in a leaflet layerGroup not producing a result

texts:
Background
I am trying to build a web map (using leaflet) that shows the location and extent of orienteering maps owned by a number of clubs/organisations in my state. Each club/organisation has a layer that is populated with their maps. To make finding these maps easier, the ability to search by name is required.
The Problem
To be able to search all layers at once I have grouped them into a layerGroup (as suggested in this answer). This does not produce a result when I attempt to search for a map (the developer console throws the error TypeError: a.getLatLng is not a function in relation to leaflet-search.js). Searching with only one layer does produce a satisfactory result (and doesn't throw an error).
Here is my current code
The webmap with search through multiple layers
The webmap with search through only one layer
How can this error be solved? Is there a better method to search through multiple layers?


texts:
I am looking for an accurate Map for visualizing positional data collected in Kenya. Right now I am using EXOMAPS (https://tiles.maps.eox.at/wms), in particular, a Sentinel2 cloudless view, but I am getting as you can see a poor resolution.

The points have an average distance of 1 meter.

I was checking at the WMS Global Mosaic (https://gcmd.nasa.gov/KeywordSearch/Metadata.do?Portal=NASA&KeywordPath=Projects%7CM+-+O%7COGC%2FWMS&EntryId=NASA_JPL_WMS_GMR&MetadataView=Full&MetadataType=0&lbnode=mdlb1), but I cannot find the WMS server APIs.

I am open to suggestions. I would have like to have at least a resolution of 2 or 3 meters (as the one of Google Maps).


texts: How to specify the desired WMS version in Quantum GIS?

texts:
How can i set a WMS version of layer which i want to display. I am asking because some WMS addresses are still served via older versions like 1.1.1.

In short, I am looking for the equivalent of this one from ArcGIS:


texts: Measure the angle of each corner of a polygon

texts: How can I automatically measure the angle for each corner of a polygon (in a vector GIS) using VBA and ArcObjets?


texts: Finding pseudo nodes in free GIS software?

texts:
Software gvSIG OA Digital Edition 2010 have tools topology for finding pseudo nodes in linear geometry. I set cluster tolerance 0.00002 and maximum number of errors -10000 for 20000 link count linear geometry. But unsuccessful result.

Are there any solutions that find pseudo nodes in free GIS software?

I need to layer pseudo nodes (one solution to this problem - to use tools topology of ArcInfo, but priority for me is to use free software). Linear geometry created several users in QGIS 1.8.0 in PostGIS (v. 2.0.1) database.

Add new image: 12 linear features with three pseudo nodes in A (line 4/5), B (line 6/7), C (line 9/10). Pseudo nodes should be points instead - two linear features with intersection in one point (node) should be one linear feature (line 4/5 - line 4, ...).

Is it possible to make a request in PostGIS, which will result in a layer of pseudo nodes?

Add new image of examples pseudo nodes: if I receive for linear layer point layer pseudo nodes (blue rects) I corrected following errors in linear layer: A - add missing geometry, B - snapped line in intersection, C - remove pseudo node.


texts:
I want to load mvt tiles into QGIS through a URL service

  https://abc.com:8080/buildings/{z}/{x}/{y}.mvt

But, the menu in the Vector Tile Reader plugin allows TileJSON URL.

Is there any way I can load mvt tiles in QGIS environment?


texts: Python folium package for 'satellite' map

texts:
I am trying to have a zoomable satellite map for my Transportation network analysis. Other free maps all OK but could not get satellite to work. I am using:

token = 'pk.eyJ1IjoiYXNhZGkxOTc4IiwiYSI6ImNqanphbHVhMDA4ODgzd25yYzY1MW1nc2oifQ.dDDcWq2-tUC-XXXXXXX'
tiletxt = 'https://api.mapbox.com/v4/mapbox.satellite/3/2/3@2x.png?access_token=' + str(token)

my_map = folium.Map(location = [lat,lon], tiles = tiletxt, attr='WHAT SHOULD I PUT HERE?', zoom_start = 7)

But give a grid static map which cannot be zoomed etc.

Can I get some instructions how I can get a nice satellite map?


texts: How Do I search for MODIS data using pyModis?

texts:
I want to download the latest Modis vegetation Index Data(MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006) , whenever it is available for my Study Area.

I know the Tile that my data falls in, but I can't figure out how to search, and get a list of available tiles in a Given Time period.

Is there a sample script which shows how you can search for Modis data using pyModis?


texts: How to fetch multiple images using Sentinel-Hub's EO Browser Custom script?

texts:
I'd like to use the floods Sentinel-hub Playground custom script in EO Browser so that I may download the result as geoTiff. https://github.com/sentinel-hub/custom-scripts/tree/master/sentinel-1/flood_mapping#
The Playground script compares two S1-GRD VV images at different points in time and returns the difference.
Bringing it into EO Browser's Custom script area gives the error:
An error has occurred while fetching images:
Dataset with id: 0 not found.
There is a checkbox named "[X] Use additional datasets (advanced)" where you can select "S1GRD" as a source, which I'd expect would allow me to fetch additional images but I still get the error.
Is it possible to use multitemporal data in EO Browser's custom scripts?


texts: Getting extent of a point cloud using laspy

texts:
As part of a geo-indexing script, I'd like to add in support for numerous .laz files we have sitting around on our SAN.  The only intent of the handler is to grab the filename/path, the SRS, and the extent of the point cloud file.  I do not need to do any processing or analysis on the dots.

So I'm able to do

from laspy.file import File

my_file = File("path/to/file.laz", mode='r')

and it loads the files without error.  I can call the File.header property, but have no idea how to "get into" it to extract values or calculate the min/max X/Y.


texts: Apply changes to style-file to all layers using that style

texts: I'm using QGIS. I have one qml-file which I apply to many layers of different projects (so all gas-pipes have the same style). As for now, I can make a change to that qml, save it and than I have to go through all projects which contain a gas layer relying on that qml and manually load that updated qml. Is there a way to apply changes to that qml and all gas-pipe layers in all different projects get this updated style automatically?


texts: Kriging in QGIS 3.4 Problem with Variogram Model

texts: In QGIS 3.4.13 - Madeira there seems to be a bug in the Ordinary Kriging Modul (Processing Toolbox > SAGA). When I enter the modell "1600000*x+19202*x^2" QGIS will add a second "^" leading to wrong results. Am I doing something wrong? Is this bug known?


texts: ST_DWithin: Why does query not use the Spatial Index?

texts:
I want my query to output all street segments, that are within a distance of 20 meters of those segments which have the attribute clazz of 11 (motorways).

I have this query from a user on this site:

SELECT DISTINCT s.id, s.clazz
FROM berlin_2po_4pgr s
INNER JOIN berlin_2po_4pgr m
ON ST_DWithin(s.geom_way::geography, m.geom_way::geography, 20)
WHERE m.clazz = '11' AND s.clazz != '11';

I have a index:

CREATE INDEX gist_index_berlin_2po_4pgr ON berlin_2po_4pgr USING gist(geom_way);

The query takes a very long time (more than two minutes). Probably the reason is that the index isn't used:

"HashAggregate  (cost=8076781.04..8076781.13 rows=9 width=8)"
"  ->  Nested Loop  (cost=0.00..8076781.00 rows=9 width=8)"
"        Join Filter: (((s.geom_way)::geography && _st_expand((m.geom_way)::geography, 20::double precision)) AND ((m.geom_way)::geography && _st_expand((s.geom_way)::geography, 20::double precision)) AND _st_dwithin((s.geom_way)::geography, (m.geom_way)::g (...)"
"        ->  Seq Scan on berlin_2po_4pgr s  (cost=0.00..2267.43 rows=49871 width=101)"
"              Filter: (clazz <> 11)"
"        ->  Materialize  (cost=0.00..2270.24 rows=563 width=93)"
"              ->  Seq Scan on berlin_2po_4pgr m  (cost=0.00..2267.43 rows=563 width=93)"
"                    Filter: (clazz = 11)"

Why? The SRID is 4326, so unit is degrees. But the cast to geography changes the units to meters, right?

The output is correct, so unit can't be the problem.


texts: If-then-else Statements and string errors?

texts:
I'm defining a string field based on contents of another string field.
Simply put, if the contents of field A contain this value, then fill field B with this text.

It throws a string related error message which I don't understand how it fits into the error.

What am I not seeing here?


texts: How to get mapPoint z coordinate by Lat and Lon

texts:
Please see the map on:
https://codepen.io/stzdzyhs/pen/WamYav

the whole code is also inserted below.

When click on the scene(A IntegratedMeshLayer), it can show the Lon, Lat, and Z values in the SceneView onclick event.

however, how to query the Z value by inputting (Lon,Lat), e.g:
whether it has a function like:
double getZ(lon, lat);

the code:

<meta charset="utf-8">
<meta name="viewport" content="initial-scale=1,maximum-scale=1,user-scalable=no">

<title>IntegratedMeshLayer - 4.6</title>
<link rel="stylesheet" href="https://js.arcgis.com/4.6/esri/css/main.css">


texts: QGIS 2.14 - Can't read message while editing

texts:
While creating a new feature (polygon), I get a blue message for some vertices but the pop-up message is just too fast for me to read or click on any pause button.

How can I set the time a message is displayed? Is there a log for all messages?


texts: Any plug-in that allows layer editing in the MirrorMap window?

texts:
QGIS 'Dockable MirrorMap' I am using a plugin. This plugin is able to add one or more extra panels in the main window with additional map canvases.

The 'Dockable MirrorMap' plug-in simply has the ability to synchronize its location in the expanded window, is there a plug-in that allows layer editing in the expanded window?


texts: Which Opensource products to use for Maprendering/Hosting?

texts:
I'm trying to setup a server for hosting maps for a gis application. The path I have taken so far is using an UBUNTU 10.04LTS installation. I'm using PostgreSQL 8.4 and PostGIS 1.5. I already made a few attempts at using Mapnik and right now I'm struggling with TileCache to cache my maps.

I recently asked myself the question - why?
Has anyone got an overview with different server-tools to setup such a server. Why mapnik and not mapserver or why tilecache and not tirex etc. Is there a proven strategy or concept to follow?

My only information is in this PDF

Right now, I just need to be able to host maps for my old prototype app in js/openlayers. I'd like be able to show that I am doing all I can now to help decide in which direction the appication I'm working on should go.

Thanks for your help!
Martin


texts: Reordering layers in GeoPackage and changing the default naming in QGIS?

texts:
I have successfully saved several layers in one GeoPackage in QGIS 3.8. This is how it looks like in DB manager:

I found out that the ordering of the layers is somehow strange when loading the GPKG file in QGIS:

Of course, I can manually change the layer ordering by dragging them in the Layers panel up or down but it would be nice if the default order could be correct by default. Same with the layer names in QGIS - is there any way to exclude the filename from the layer name?

That means get the layers loaded this way by default?

Because the support for multiple layers in one GeoPackage file is the much less useful if I have to manually adjust the names and ordering each time.

I tried to import the layers from my GPKG to another one but the resulting order was the same... Is there any way to solve this?


texts: How to Make QGIS on Windows use Spatialite 3.0 instead of 2.4?

texts:
I spent all afternoon trying to figure out why my complex triggers are working fine in Spatialite-GUI but QGIS refuses to write to the layer.  Finally remembered I ran into a similar issue a few months back.  Problem is, there's no workaround this time.

I'm using triggers to force a point to snap to the nearest line on another layer, using the midpoint if no geometry is supplied, or the nearest point if it is.
Works great in 3.0, but sadly, almost all of the geo functions I'm using (Interpolate Line, Shortest Distance, etc) are not on the list of functions in 2.4.

So... I know there's been discussion about removing the built-in Spatialite from QGIS.  That seems to imply there'd be a way to point QGIS at Spatialite installed somewhere else.  Does anyone know how to make that happen?  Running Windows 7, QGIS 1.7.4.  Sorry if it sounds stupid, but I've never compiled anything so I have pretty much no idea what I'm talking about.


texts:
When I try to bring fieldwork back to my desktop, QField seems to be asking for an import folder. Exploring the problem a little, it seems like folks copy/paste their field work data into an existing folder on their home computer, rather than using the import folder.

Could someone please explain the steps/workflow for getting fieldwork updates off my mobile device and back into the QGIS project?


texts: Displaying information on a map using an InfoWindow

texts:
I am wanting to display information in an InfoWindow when the user clicks on a point on the map. I have implemented the InfoWindow in XAML, but I do not understand the code behind it. Any information will be appreciated. Here is my (what I feel is relevant) XAML:

 <esri:InfoWindow x:Name="MyInfoWindow"
                             Padding="2"
                             CornerRadius="20"
                             Background="{StaticResource PanelGradient}"
                             Map="{Binding MyMap}"
                             ContentTemplate="{StaticResource MyFeatureLayerInfoWindowTemplate}"
                         MouseLeftButtonUp="MyInfoWindow_MouseLeftButtonUp" />

If you need anymore information, let me know. I have the information being displayed correctly at the moment in a MapTip. Here is the XAML for that:

<esri:FeatureLayer.MapTip>
                        <Border CornerRadius="10" BorderBrush="Black" BorderThickness="2" Margin="0,0,15,15" Background="{StaticResource PanelGradient}">
                            <StackPanel Margin="7">
                                <TextBlock Text="{Binding [WVAL_ID], StringFormat='ID: \{0\}'}" FontWeight="Bold" Foreground="White"  />
                                <TextBlock Text="{Binding [SIZE], StringFormat='Size: \{0\}'}" FontWeight="Bold" Foreground="White" />
                                <TextBlock Text="{Binding [WVAL_USE], StringFormat='Use: \{0\}'}" FontWeight="Bold" Foreground="White" />
                                <TextBlock Text="{Binding [WVAL_TYPE], StringFormat='Type: \{0\}'}" FontWeight="Bold" Foreground="White"  />
                                <TextBlock Text="{Binding [OPENS], StringFormat='Opens: \{0\}'}" FontWeight="Bold" Foreground="White"  />
                                <TextBlock Text="{Binding [TURNS], StringFormat='Turns: \{0\}'}" FontWeight="Bold" Foreground="White"  />
                                <StackPanel Orientation="Horizontal" >

                                </StackPanel>
                            </StackPanel>
                        </Border>
                    </esri:FeatureLayer.MapTip>


texts: Why is there no state with code "03" (among other missing numbers)?

texts:
The US census documentation for state codes is missing several numbers. For example, there is no state with code "03", "07" and others.

Here is the start of the list...

STATE|STUSAB|STATE_NAME|STATENS
01|AL|Alabama|01779775
02|AK|Alaska|01785533
04|AZ|Arizona|01779777
05|AR|Arkansas|00068085
06|CA|California|01779778
08|CO|Colorado|01779779
09|CT|Connecticut|01779780
10|DE|Delaware|01779781
...

What exactly happened here?
It's not as though there are constantly new states joining and old states leaving. This list should be pretty stable and predictable.


texts: Comparing two maps with two different libraries

texts:
I'm currently making a demonstration to compare the two OpenLayers libraries (2 and 3). For doing this, I have made two div elements and put the different maps on each other. The goal is to make a map which starts in the first div and continues in the second div, fluently. The result can be seen on the following fiddle: http://jsfiddle.net/GFarkas/xxqu97nc/.

To synchronize the two maps, I have to bind events to each one and pass the corresponding center coordinates to the other map.

My basic concept is, that if we want to get the maps aligned, they have to share lower left coordinates (because one of the divs is smaller than the other). Because they are the same height, the latitude won't change, it can be passed from the event object. However, because of the difference in their width, the longitude have to be calculated. As they have to share the same lower left coordinates, the center can be calculated by averaging the longitude at the sides of the corresponding bounding box:

view.setCenter([(map.getLonLatFromPixel(new OpenLayers.Pixel([0,0])).lon +
map.getLonLatFromPixel(new OpenLayers.Pixel([map2.getSize(),0])).lon)/2,
evt.object.center.lat]);

to change the OpenLayers 3 div, and

map.setCenter([(map2.getCoordinateFromPixel([0, 0])[0] +
map2.getCoordinateFromPixel([640, 0])[0])/2,
evt.map.getView().getCenter()[1]], evt.map.getView().getZoom());

to change the OpenLayers 2 one.

This method works as intended until I get near to the date line (where the positive and the negative longitudes meet as we're in Pseudo Mercator). When I get between the -19000000m and 20000000m longitudes, this simple equation to calculate the center coordinates fails and my view goes near to the origin.

The question is, does anybody have a better and safer method to calculate coordinates to synchronize the two maps? I have tried to create an OpenLayers.Bounds to find out the center coordinates, but no chance, it still fails.

Oh, and the divs must go on each other, because I'm swiping between them, just JSFiddle refused to maintain that function in my code.


texts: Coordinates of an area on a raster in QGIS

texts:
I am new to GIS and I have to do the following:

I have a GeoTIFF raster file of a map with some features like clouds and water bodies. It is a satellite image. My images are taken from a snapshot of https://worldview.earthdata.nasa.gov/ which gives me a GeoTIFF file.
I have to encircle a particular feature like a particular water body or some area of the cloud coverage.
The area I encircle would give me the coordinates in (lon/lat) of the area so that I can use the coordinates in programming contour plots.
It is better if I can get the coordinates as a CSV file.

What I have tried so far:

I added the raster TIFF image.
Used the New Shapefile Layer Icon to create a disastrous polygon to encircle a particular area.
A shapefile (vector layer) gets added on top of the raster map.
I tried exporting the geometry coordinates to a CSV file, but was unsuccessful.
I tried the MMQGIS plugin and export the shp coordinates to a CSV file, but it was of no avail, because I got only the four boundary coordinates.

Is there any better way in getting around this stuff? Especially in drawing a freeform vector and getting the coordinates of the encircled area.


texts:
I have read al lot of the pgrouting, but i still not quite sure if it can do what i want.

I have 20 stop points and i have 1 starting point and 1 finish point.
I want to calculate the shortest path between starting point, all 20 stop point and finish point.

Its something like delivering packects, so from starting point to point x then point y and last finish point,

Maybe iets possible with a dijkstra algorimte ?


texts: Changing values in attribute table and split it into different columns?

texts:
How to change the values  in the table and divide it into different columns?
I have an example

I want to get a result


texts: How to Send a Feature to a Different Layer

texts: I have a polygon that I made using the buffer tool on a layer. The problem is, I can't figure out how to send this polygon to the layer I would like it to be on. The yellow polyon is the one I'm trying to move, and I'm trying to move it to the red layer. How can I do this using ArcGIS 10.1 for desktop?


texts: How to Join two feature classes without common attribute

texts:
I have two shapefiles of essentially the same objects however, one shapefile has terrible spatial data but great attribute data and the other is the reverse. What I'm looking to do is essentially a join but there are no common fields that I can use as the joining attribute.

Without that common field, I don't believe this is possible, but I thought I'd ask in case someone has done this in the past.


texts: For multiple rows of values on singular route feature, how does ArcMap symbolize by quantities?

texts:
So here is my situation:

I have a road segment that is one feature. Along this road I have data that is captured via GPS pings (which generally occur every 5 seconds), and are delineated by 'START_DATE' and 'END_DATE'; so there are multiple pings that make up the information per 'capture' or row.

The question I have is with the fact that I have multiple rows of values that encapsulate the entire segment and therefore occupy the exact spatial area of the line.

One of the columns in the table is 'QUANTITY_SOLID', which contains the total amount of solid materials that were disseminated during that 'START_DATE / 'END_DATE' instance.

I then attempt to symbolize by 'Quantities > Graduated Colors'

I understand that because it is only one feature class, and all of the information contained in the rows along that line segment all cover the same singular feature class, that there will not be any differentiation between the colors when symbolized.

However, this leads to my actual question: How does ArcMap choose what color to make the line if there are multiple values over the same spatial area? Does it take the sum of all the 'QUANTITY_SOLID' and average them? Or something else?


texts:
There is such a task: you need to display one of the raster layer and  one or several vector objects (use layer CQL of filter).

I use geoserver to store the geodata and leaflet js on the client to display geodata.

One "beautiful" solutions is to use WMS with several layers and filters for vector layer.

For it i must define bbox. How can i do it?

I have input data:

size of the picture
bounding box of geoobjects

I need result: bbox

For examples


texts: How to get the maximum value of a vector layer without filters?

texts:
I'm trying to get a maximum value of a field in a vector layer without using the filter that is active in QGIS. So far I've tried:
layer = iface.activeLayer()
idx = layer.fields().indexFromName('fid')
print layer.maximumValue(idx)

This sends the following SQL query to my SQL Server database (I have filtered the layer by fid in QGIS):
select max([fid]) from [mySchema].[myLayer] where ("fid" < 1000)

Is it possible to send
select max([fid]) from [mySchema].[myLayer]

instead? I'd still like to have the filter on, so that not all features are shown on the map.


texts: Approach to Geographic Location Optimization in PostGIS

texts:
My goal is to take a geocoded list of targeted consumers, a shapefile of "regions", and create an optimized number of "district" shapes (a geographic child to or subregion of larger regions). I'd like to ensure districts fit within their parent regions and create the smallest number of districts while maximizing the density of consumers inside of the districts.
Let's call the consumer shapefile consumer_points, the region shapefile office_regions, and the generated district shapes office_districts. The scope is within a single US state.
So far, what I've got:

Assign consumer to Region and append region to consumer_points using ST_WITHIN()
A zip_code shapefile to use as a building block for districts (i.e. the most granular assembly shape for a district).
Assigned consumers to their corresponding zip shape using ST_WITHIN()

consumer_points now has: person_id, latitude, longitude, zip_assignment, zip_geom, region_assignment, region_geom
zip_codes now has: zip_code, zip_geom, region_assignment, region_geom
I'm struggling to figure out the next step or potential solutions.


texts: Showing differents label's position for the same vector object in two different view maps in the QGIS Print Composer?

texts:
I´m trying to do a "detail map" (one window with a bigger scale of one part of the map) over my "main map" in the print composer. But in this detail map I need to change the position of the labels in some specific objects of my vector's layer, but the problem is that the "Lock layers for map item" option does not lock this.

For example in these images:

1) The position of the labels in red is this for the "main map" (scale 1:1000)

2) The position of the labels in the "detail map" should be more close to the objects, like this:

3) But if I change the position like in the image above, in the "main map" the labels will be shown like that:

The labels are over the objects and even some of them are hidden

There is any way to do that?


texts:
I have a query that returns the closest distance between a set of points and multiline:

SELECT
   DISTINCT ON (h.field_1) h.field_1, b.id, ST_Distance(h.wkb_geometry, b.wkb_geometry)  as dist
FROM points As h, lines As b

However, given that the projection is EPSG 4326, the distance returned is in degrees which doesn't make much sense.  How could one rewrite the query and cast to ::geography so that the distances returned are in meters ?


texts:
I use the ModelBuilder to run my model
and I have question

Sub M()
' Create the Geoprocessor object
Dim gp As Object
Set gp = CreateObject("esriGeoprocessing.GPDispatch.1")

' Check out any necessary licenses
gp.CheckOutExtension "spatial"

' Load required toolboxes...
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx"
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx"
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx"

' Local variables...
Dim Output_Features As String
Dim Output_point_features As String
Dim Layer_Name_or_Table_View As String
Dim Kriging_poin3 As String
Dim Output_variance_of_prediction_raster As String
Dim dbtexts_mdb As String
Dim ws As String
Dim dbtexts_mdb__2_ As String
Output_Features = "C:\temp2\RasterT_Kriging1.shp"
Output_point_features = "C:\temp2\RasterT_Kriging1.shp"
Layer_Name_or_Table_View = "ws_Layer"
Kriging_poin3 = "C:\temp2\kriging_poin3"
Output_variance_of_prediction_raster = ""
dbtexts_mdb = "C:\dbtexts.mdb"
ws = "C:\dbtexts.mdb\ws"
dbtexts_mdb__2_ = "C:\dbtexts.mdb"

' Process: Make XY Event Layer...
gp.MakeXYEventLayer_management ws, "x", "y", Layer_Name_or_Table_View, ""

' Process: Kriging...
gp.Kriging_sa Layer_Name_or_Table_View, "data", Kriging_poin3, "Spherical 249.101224",     "250", "VARIABLE 12", Output_variance_of_prediction_raster

' Process: Raster to Point...
gp.RasterToPoint_conversion Kriging_poin3, Output_point_features, "Value"

' Process: Add XY Coordinates...
gp.AddXY_management Output_point_features

' Process: Table to Geodatabase (multiple)...
gp.TableToGeodatabase_conversion "C:\temp2\RasterT_Kriging1.shp", dbtexts_mdb__2_

End Sub

in the sub I have parameters but I want it to run with the data I have

can I just like VB give it a do loop and make it run itself?


texts:
I have a feeling I know the answer to this but want to make sure before I start redoing all the work.

I've recently started using QGIS to help with mapping of an estate.  Last week a created a new vector polygon layer (and associated shapefile) to map the woodland.  By Friday afternoon I'd completed the mapping and entry of all the attributes for each polygon.

I opened up QGIS this morning to show a colleague the progress and the layer is at a point I reached on Thursday afternoon.  All my updates (the bulk of the work) made on Friday aren't there.  The confusing bit is that when I check the underlying shapefiles in Windows Explorer the .dbf .shp .sxh files are all showing a last modified date of 1638 on Friday, when I last saved before the weekend.

I've tried creating a new vector layer based on these but with no success.

Update - I've now recreated all the polygons and associated attributes.  But have been thinking more about what could have caused the original issue.  If I had only been saving layer edits during the Thursday and Friday and then not saved the project before going home on Friday, could this have potentially caused the issue? I can't be certain this is what I did but trying to find an explanation so I can avoid it in the future.


texts: How to generate dynamic legend using ArcPy?

texts: I used this  example to produce printing service. As the result I get map with too many elements in the legend. Is it possible ti show classes that are visible in the current map extent using python script?


texts: Python GDAL: Rasterising point data, raster and points are not geolocated

texts:
I am trying to rasterise some point data using Python GDAL module. But when I plot the data there is an offset between the two (original points and raster). I am sure it has to do with the geotransform values.

Here is my code and data can be downloaded from here.

import pandas as pd
import gdal
import numpy as np
import matplotlib.pyplot as plt

pts = pd.read_csv("SAMPLE_POINTS.csv")

# here I grid the data and count how many points are in each cell
pts['XX'] = pts.X // 10000
pts['YY'] = pts.Y // 10000

COUNT = np.empty((len(pts.XX.unique()), len(pts.YY.unique())))

for i, XX in enumerate(np.sort(pts.XX.unique())):
    for j, YY in enumerate(np.sort(pts.YY.unique())):
        c = pts.ix[(pts.XX > XX) & (pts.XX <= XX + 1) &
                   (pts.YY > YY) & (pts.YY <= YY + 1)]
        COUNT[i, j] = len(c)

# Setting gt is the bit I am sure is wrong
# This is taken from SetGeoTransform() I would use to export raster
gt = [pts.X.min(), 10000, 0, pts.Y.max(), 0, -10000]

# plot raster
cmap = plt.imshow(np.ma.masked_equal(np.rot90(COUNT), 0), interpolation="none",
                  extent=(gt[0], gt[0] + (gt[1] * COUNT.shape[0]),
                          gt[3] - (gt[1] * COUNT.shape[1]), gt[3]))
plt.colorbar(cmap)

# plot points
plt.scatter(pts.X, pts.Y, s=1, c='k', edgecolor="none")

# close up of misalignment
plt.xlim(450000, 500000)
plt.ylim(5850000, 5900000)
plt.show()

I don't think it is plotting the data in Python that is wrong because when I export the raster and display in a GIS there is the same offset.

Here are a few of images illustrating the offset. The points and the rasters are orientated correctly but the raster is offset south by maybe 2 pixels.


texts:
I want to merge the selected features. I found this post but I can not do it.
 Can anybody help me?

layer = iface.activeLayer()
obj1 = 24
obj2= 25
expr = QgsExpression( " \"ID_PARCEL\" = '{}' OR \"ID_PARCEL\"='{}' ".format(obj1,obj2))
it = layer.getFeatures( QgsFeatureRequest( expr ) )
ids = [i.id() for i in it]
if ids == []:
    print 'ERROR - WAS NOT FOUND'
else:
    layer.setSelectedFeatures( ids )
    selectFeatures = layer.selectedFeatures()

#MERGE SELECTED FEATURES


texts: Java library to break up a big polygon into many parts

texts: I need to break up large polygons like these into smaller manageable pieces that can be stored in my DB. I tried searching for tools to do this and mostly come across using software such as ArcGIS or QGIS. Is there no support for this in JTS or Geotools? I need to do this in a java program. As this is a complex task, I would much prefer to use a high quality tested and tried library rather than develop my own code to do this. My requirement is that if I have a JTS geometry that is more than 1MB when serialized into WKB, I need a method which can split it into multiple geometries each no bigger than 1MB. Other than that I have no constraint.


texts:
I have an ArcPad 10.2 project that uses several Shapefile layers.

How can I filter the features in one of those layers?

According to my searches, I can only filter if it is an AXF layer.

I'm happy to resort to scripting in vbs if required.

The layer I have has features that each have a status field, and it is this I wish to be able to quickly filter on.  I have to keep it as a shapefile really, as they are getting auto generated out of another system.


texts: Adaptive raster to points sampling based on change of slope in QGIS?

texts:
I have a high resolution elevation map and need to convert it to XYZ points. Simple conversion generates a lot of points in flat / constant slope areas which are redundant for 3D terrain generation. The goal is to reduce number of points by sampling only the points when there is a significant "break" in terrain surface. This would be equivalent to mesh simplification in other software.

I was looking into filtering points by Terrain Ruggedness Index, but from what I understand this index indicates the uneveness of the terrain, not change in gradient.