texts: How do you create a line network with GeoJSON navigation point data

texts:
I have a GeoJSON point file for use in navigation.  I need to create a line network connecting the points to their corresponding 'neighbors' in the data.  I have looked for any sort of plug in or native algorithm but have not found anything that can perform the connection.  Is there a way to do this?  Here is a sample of the data.

 {
"type": "FeatureCollection",
"features": [
    {
        "type": "Feature",
        "geometry": {
            "type": "Point",
            "coordinates": [
                -122.51407721663476,
                37.92935472657569
            ]
        },
        "properties": {
            "id": "n_59d3a0656850d07d450000fc",
            "level": "59d39a1a6850d04750000000",
            "neighbors": [
                "n_59d3a1c66850d07d4500010c",
                "n_59d3a1cd6850d07d4500010d",
                "n_59d3a7f26850d07d4500015e"
            ],
            "weight": 1
        }
    },
    {
        "type": "Feature",
        "geometry": {
            "type": "Point",
            "coordinates": [
                -122.51377017089081,
                37.92918570092726
            ]
        },
        "properties": {
            "id": "n_59d3a1c66850d07d4500010c",
            "level": "59d39a1a6850d04750000000",
            "neighbors": [
                "n_59d3a0656850d07d450000fc",
                "n_59d3a2f66850d07d4500010f"
            ],
            "weight": 1
        }
    },
    {
        "type": "Feature",
        "geometry": {
            "type": "Point",
            "coordinates": [
                -122.51388322052854,
                37.92957582123205
            ]
        },
        "properties": {
            "id": "n_59d3a1cd6850d07d4500010d",
            "level": "59d39a1a6850d04750000000",
            "neighbors": [
                "n_59d3a0656850d07d450000fc"
            ],
            "weight": 1
        }
    },
    {
        "type": "Feature",
        "geometry": {
            "type": "Point",
            "coordinates": [
                -122.51414818592386,
                37.92934461122832
            ]
        },
        "properties": {
            "id": "n_59d3a7f26850d07d4500015e",
            "level": "59d39a1a6850d04750000000",
            "neighbors": [
                "n_59d3a0656850d07d450000fc",
                "n_59d3a0e46850d07d45000106",
                "n_59d3a7f66850d07d4500015f",
                "n_59f354f824ee960562000000"
            ],
            "weight": 1
        }
    }
]

}


texts: Resizing the map item in QGIS print composer without changing scale?

texts: Is it possible to lock the scale in the QGIS 2.0.1 print composer? Every time I resize the map, the scale is adjusted. This makes trying out different map layouts really annoying.


texts: Adding SpatiaLite to Visual Studio project (already has SQLite)

texts:
I am very new to SpatiaLite.

I am working on a large Visual Studio C++ project that already supports SQLite, and I am trying to add SpatiaLite capabilities to it.
After reading the documentation and trying different approaches I am still not sure how to approach this.

The main C++ project loads a separate SQLite Visual Studio project in order to use database functions.
In this SQLite Visual Studio project, we have SQLite headers, DLLs, and LIBs. We use the linker to access them from our main C++ project.

Do you think it is viable to try to add the SpatiaLite extension to the SQLite Visual Studio project?
I am not sure we are approaching it correctly. When trying to simply add header files + spatialite.c, and using OMIT_ICONV, OMIT_GEOS, OMIT_PROJ, OMIT_FREEXL in the preprocessor for building, it has many linker problems.
For example:
1>spatialite.obj : error LNK2019: unresolved external symbol _gaiaDoubleQuotedSql referenced in function _checkSpatialMetaData_ex

What are the minimal dependencies needed, are we missing files?

Is there a purpose of building SpatiaLite in MinGW if our main C++ project is built and run from Visual Studio?

I've tried to build SpatiaLite with MinGW 64 bit instructions, but failed.

On libjpeg step, ran into this error when trying to do the line: cmake -G "MSYS Makefiles" -DCMAKE_INSTALL_PREFIX=/mingw64/local ..

-- The ASM_NASM compiler identification is unknown
-- Didn't find assembler
CMake Error at simd/CMakeLists.txt:41 (enable_language):
  No CMAKE_ASM_NASM_COMPILER could be found.

Tell CMake where to find the compiler by setting either the environment
  variable "ASM_NASM" or the CMake cache entry CMAKE_ASM_NASM_COMPILER to the
  full path to the compiler, or to the compiler name if it is in the PATH.

I tried to skip down to libproj and libgeos, to just build spatialite, and libgeos on "make" step had failures.

error lines:
libtool: compile:  g++ -DHAVE_CONFIG_H -I. -I../../../include -I../../../include                                                                                /geos -I../../../include -DGEOS_INLINE -pedantic -Wall -ansi -Wno-long-long -ffl                                                                                oat-store -03 -MT IndexedPointInAreaLocator.lo -MD -MP -MF .deps/IndexedPointInA                                                                                reaLocator.Tpo -c IndexedPointInAreaLocator.cpp
g++.exe: error: unrecognized command line option '-03'


texts: Calculating area in the attribute field of a QGIS 3.4.3 project with multiple CRS

texts:
The area calculations in my project derived from the attribute field show values way below actual area, because the attribute field assumes the project CRS (EPSG:3857). However, this CRS should remain unchanged as it is the same CRS used for the Bing basemap in the project. The layers in my project are projected using ESPG:4326. This should also remain unchanged as this allows the correct lat/lon to be displayed.

Note - when using the measurement tool or adding geometry attributes the correct area (or close to correct) can be displayed by calculating ellipsoidal area.

My project requires that the attribute fields are pre-determined, so the previously mentioned techniques do not seem viable solutions.

So, given that the project CRS is not suited for calculating area... is there any way to correctly calculate area in the attribute field of a QGIS 3.4.3 project with multiple CRS?


texts: Opening a custom model's geoprocessor window with ArcObjects

texts: I have developed a custom model using ModelBuilder in ArcGIS, and I would like to use a .NET Base Command to open the geoprocessor parameter window for this tool so the user can enter their parameters, same behavior as if I were double-clicking it in ArcToolbox, as opposed to just executing it programmatically without user interaction. Is this possible?


texts: I know that GeoJSON is supported in Leaflet but how about KML? Is there any plugin for reading KML data in Leaflet?


texts: I'm trying to set up a WPS using PyWPS, but when I return GML the QGIS WPS client says process result unsupported MIME type.  Any idea why that is?


texts: I am currently in the process of converting a large quantity of Ordnance Survey Vector Map Local data from .gz to Esri feature class within a file geodatabase. The process is working correctly but there is a limit to the number of .gz files FME will let me load/convert in one go. The error message simply says "there are too many recursive macro subsitutions". My knowledge of FME is limited so I was just wondering if there is a way to get around this as I currently have to keep loading sections of the data into the program.


texts:
I would like to use Google Maps Elevation API to query elevation data for my future free and open source script.

What I would like to know is what is the exact text which I need to include in my license documentation in order to provide the appropriate credit to Google corporation, for using their Google Maps Elevation API?

I checked Google Maps API terms of service page, but I didn't understand which text to include to my license documentation.

Where can I find this info?

I tried contacting Google Maps directly, but never got a reply from them.


texts: Saving outputs of "proccesing" to temp file in PyQGIS

texts:
I'm trying to save the output of a "proccesing" algorithm executed in PyQGIS as a temporal file, and then re-using it with another algorithm. The desired result is to write the output only at the end of the script. I'm working in QGIS 3.2.

    from qgis.core import *
    from PyQt5.QtGui import *
    import processing

path='path/to/output/'

parameters = { '-u' : False,
'GRASS_RASTER_FORMAT_META' : '',
'GRASS_RASTER_FORMAT_OPT' : '',
'GRASS_REGION_CELLSIZE_PARAMETER' : 100,
'GRASS_REGION_PARAMETER' : '565000,600000,6220000,6240000 [EPSG:32721]',
'distance' : 0.5, 'exponent' : 0.5, 'flat' : 0, 'high' : 255,
'output' : path + 'rast.tif', 'seed' : None }

processing.run('grass7:r.random.surface', parameters)

parameters = { 'GRID_IN' : path + 'rast.tif',
'GRID_OUT' : path + 'rast_reclas_v1.sdat',
'LOOKUP' : [0,150,0,151,255,1], 'METHOD' : 2 }

processing.run('saga:reclassifyvaluessimple', parameters)

parameters = { 'COPY_SUBDATASETS' : False,
'DATA_TYPE' : 5,
'INPUT' : path + 'rast_reclas_v1.sdat',
'NODATA' : None, 'OPTIONS' : '',
'OUTPUT' : path + 'rast_v1.tif',
'TARGET_CRS' : 'EPSG:32721' }

processing.run('gdal:translate', parameters)


texts: I am building an app using qgis in Windows 7 with VS 2008.  It was diagnosed that I must have the wrong versions of the libraries because of a potential name mangling problem.  I have 12 unresolved externals and I can find them in gis_core by micreosoft name.  Does anyone have a different library for me to link to?


texts: How to avoid topology errors when digitizing a road network?

texts:
I am working on a road layer project. For the snapping option I'm using the mode to vertex and segment and tolerance 5 pixels. When zooming in (to the max) I do not notice a gap. But when I use the topology checker in QGIS, with the option must not have dangles I notice that there are many dangles (not only on the beginning and end point) of my poly line.
Also I would really like to avoid these type of errors, instead of correcting them in the final phase.
How  can I  avoid these kind of errors?

Currently I'm using QGIS 2.10.1


texts: Exposing function or transformation matrix that QGIS georeferencer develops to georeference raster?

texts:
The problem:
Think of the London Tube Map.  Logically correct, but physically inaccurate.  I want to display such an arist's impression (un-transformed, so recognisable) on a mobile web site and have the phone's GPS accurately place the user on the logical map.  (The tube map is an example - it's not underground :))

I'm currently using CV2's homography functions to achieve this - given some source/destination co-ordinates it emits a transformation matrix which I can use in javascript to take a gps co-ordinate, transform it and display it on the picture.  This works well, but the QGIS georeferencer makes it so much easier to capture co-ordinates and check the result for accuracy, and I think that a gis is going to make a better job of mapping.

So I'd like to use QGIS for this job, but I can't see a way of grabbing the transformation matrix I assume is generated by the georeferencer.
Has anybody got any idea whether this is possible?


texts: Obtain the number of features in a Shapefile using NodeJS

texts: I'm using the shapefile NPM library which works well, but it doesn't provide a count of the number of features in a Shapefile. Is there some other library which can access this information (other than shelling out to ogrinfo), or some other technique?


texts:
I'm struggling to get what seems like a perfectly normal query to execute in MapInfo v10.

In SQLite Manager (an excellent  firefox plugin), this works as expected:

SELECT area_name, street_name, address, disp_lon, disp_lat
FROM mtdArea, pointAddress, streets
WHERE (streets.link_ID = pointAddress.link_ID) AND ((l_area_ID = area_ID) OR (r_area_ID = area_ID))

The equivalent in MapInfo

SELECT RationalisedMtdArea.Area_Name, RationalisedStreets.Street_Name, RationalisedPointAddress.address, RationalisedPointAddress.Disp_Lon, RationalisedPointAddress.Disp_Lat
FROM RationalisedMtdArea, RationalisedPointAddress, RationalisedStreets
WHERE RationalisedStreets.Link_ID = RationalisedPointAddress.Link_ID And ((RationalisedMtdArea.Area_ID = RationalisedStreets.L_Area_ID) Or (RationalisedMtdArea.Area_ID = RationalisedStreets.R_Area_ID))

..gives me a No Join specified between tables RationalisedMtdArea and RationalisedPointAddress.
Invalid join condition in WHERE clause.

Does MapInfo not like logical Ors?
What am I missing?
Thanks.


texts:
I have one table call "LR - Données", with 2 attributes: 'LR nom français' and 'LR nom latin'.
I have a point layer in a geopackage call "Observations" with 2 attributes: 'Nom français' and 'Nom latin'.
The attribute 'Nom français' is a value relation:

The attribute 'Nom latin' is also a value relation:

Until here, everything works fine. But sometimes, the 'LR - Données' table has update. So I set a dependence in my layer 'Observations' with the 'LR- Données' table.
Here is my problem. When I change for example a latin name in my table 'LR - Données', the attribute in layer 'Observations' change, and the old value comes into (). Which is fine in a way, cause it shows me changes. But how can I update all the value in my attribute table to get all the new values?


texts: Is there any sld styles for OpenStreetMap in ArcGis 10.1?

texts:
I want to display OpenStreetMap any possible way.

1) This map must be similar to the official  OpenStreetMap.

2) Also I must cache data. So I need the WMTS protocol. This is vital.

I can not use Geoserver for my purpose because I can not find sld styles for points.

So I want use ArcGIS Server. But I don't know about sld styles for ArcGIS.


texts: GeoPandas intersects doesn't find any intersection

texts:
I want to check if points layer that I have intersect with polygons layer that I have, as Boolean column in the points dataframe.
I have two GeoPandas dataframes, the first one is many points and looks like this:
>>>ID  geometry
0  12  POINT (5.0279 7.4547)
1  45  POINT (6.6539 12.139)
...

and the second dataframe is one layer of many different polygons  that looks something like this:
>>>name     code   geometry
0  Desert   12     POLYGON ((5.52013 13.8902, 5.5265 13.892,...)
1  Water    24     POLYGON ((5.53756 13.88472, 5.5291 13.8791,...)
...

I'm trying to check if there is any intersection between the points layer and the regions layer. For that I have determined the crs and then used intersects :
regions=regions.to_crs({'init': 'epsg:4326'})
points=points.set_crs({'init': 'epsg:4326'})
inter=points.geometry.intersects(regions.geometry)

A
The script runs with the following warnings:

FutureWarning: '+init=:' syntax is deprecated.
':' is the preferred initialization method. When
making the change, be mindful of axis order changes:
https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6
return _prepare_from_string(" ".join(pjargs))
/opt/conda/lib/python3.8/site-packages/geopandas/base.py:39:
UserWarning: The indices of the two GeoSeries are different.
warn("The indices of the two GeoSeries are different.")

Then when I check the results the only value is False, like all th epoints do not intersect:
inter.unique().tolist()
>>>[False]

*I saw on QGIS that the there are points that intersect and there are points that do not so there is no way this result is true
*I have checked the dtypes - each one of my geodataframes has one column that is geometry and called geometry.
My end goal: to add new column in the points geodataframe that will tell if it intersects the regions or not.


texts: Format individual line segments of a polygon

texts: Is it possible to apply different formats to different line segments of a polygon in QGIS? I am delineating forest stands within a property, and would like to change the formatting of one of the line segments of a polygon where it meets up against another polygon.


texts:
usually when i load an osm file with lots of features (>10000) the qgis will go not responding,

when the files loaded, layers will show up, but not the features on the map canvas, its just blank

sometimes 1 or more layers will show its content,
but theres always a layers or more wont

i need to restart the notebook so the contents will show, but it will happen again if i try to open another project

problems only occur when add new vector layer with osm files,
loading shp files or old project is okay


texts: Generating random muted/pastel colours using QGIS

texts:
Is there any way of generating random softer colours using QGIS?
I've just spent some while looking, unsuccessfully. There seems to be no way to influence the "seed" colour if they're random unlike other times when the colour/symbol is taken from the semi-default one in the options window. It is common for Landscape Character to be presented with softer colours, usually across the spectrum and there are too many categories to do it manually although this example only shows a small area.


texts: Multi-selection from Combo box in Python AddIn?

texts:
Is it possible to select more than one option from a combo box list in python Add ins for ArcGIS 10.x? Here is part of the code:

def onFocus(self, focused):

    if focused:
        self.mxd = arcpy.mapping.MapDocument('current')
        layer = arcpy.mapping.ListLayers(self.mxd, "layer Name")[0]
        self.items = []
        values = [row[0] for row in arcpy.da.SearchCursor(layer, ("Field Name"))]
        uniqueValues = set(values)

        for value in uniqueValues:
            self.items.append(value)

def onSelChange(self, selection):

        try:
            self.mxd = arcpy.mapping.MapDocument('current')
            layer2 = arcpy.mapping.ListLayers(self.mxd, "Layer Name")[0]
            arcpy.SelectLayerByAttribute_management(layer2, "NEW_SELECTION", "Field Name = '" + selection + "'")

            print "Processing. Please wait..."
            # Process: Calculate Field
            arcpy.CalculateField_management(layer2, "AUTH_STATUS", "\"Authorized\"", "PYTHON", "")
            arcpy.SelectLayerByAttribute_management(layer2, "CLEAR_SELECTION")

            print "Exporting to PDF. Please wait ..."
            arcpy.mapping.ExportToPDF(self.mxd, r"path to pdf output location" + "\\"  + selection + ".pdf")
            arcpy.RefreshActiveView()

        except:

            # By default any other errors will be caught here
            #
            e = sys.exc_info()[1]
            print(e.args[0])


texts: Importing GeoTiff or IMG formats with Ironpython

texts:
McNeel Rhino 5 application is using IronPython 2.7, and this is the application to which I would like to import either a GeoTiff or .img formats.

I googled a bit a found a couple of python modules, which are able to do this, but I couldn't find any information about IronPython ones.

Any recommendations on how could this be done in IronPython?
I found out that I could use some gdal.dll library which contains functionalities for importing both upper formats, but this library is 64 bit only. I need 32 bit one.


texts: Error when installing Geonode: cannot find django-geonode-mapstore-client>=2.1.0

texts:
I am very new to programming. I am trying to install GeoNode on an AWS Windows EC2 image, according to this tutorial (which is super, by the way!)
When I get to the pip install -e geonode stage, the error in my Command Line is this:
ERROR: Could not find a version that satisfies the requirement django-geonode-mapstore-client>=2.1.0 (from geonode) ERROR: No matching distribution found for django-geonode-mapstore-client>=2.1.0
When I look at how to install django-geonode-mapstore-client, the only webpage I can find is this: https://pypi.org/project/django-geonode-mapstore-client/ which works, but it installs django-geonode-mapstore-client 2.0.9. It seems that 2.1.0 has been bumped in a week ago and I'm not sure what to do at this point. When I run pip install -e geonode I get the same error message after installing client 2.0.9 and no other part of the paver setup process works.
I am wondering if anyone else is experiencing the same issue, or has advice on how to proceed? Is there a way to force it to install 2.1.0?


texts: Prevent browser freeze when loading lots of GeoJSON features at once

texts:
In my OpenLayers 3 (v3.15.1) application, I use to load a larger amount of GeoJSON features on my map via an ol.Vector source and a custom loading fuction. This function receives the GeoJSON data from a webserver via Ajax and consequently adds it to the vector layer's source via my_layer.getSource().addFeatures(ol.format.GeoJSON().readFeatures()).

However, for about 30000 features and some custom styling rules applied, this makes the browser freeze for quite a while. While I do not mind that loading the data takes some time, the freeze makes the user feel as if the browser just crashed which is a bad feeling as we all know. Once the data is loaded, however, speed is not a concern since I wrap the vector features in an ImageVector layer which is pretty fast.

I think that the freeze comes from some OpenLayers 3 internals which add the features to the layer synchronously in the main thread.

Is there any possibility how the freeze can be avoided? I saw some other threads on similar topics where people talked about web workers. However, none of them had a solution yet, so maybe there are different ways as well.


texts: In QGIS, how to convert northing and easting from a CSV file into longitude and latitude so it fits with the rest of my data?

texts: I'm adding a CSV file into my map where all the other data comes from shapefiles. I tried messing with the coordinate system, but I can't get the data to show up as datapoints on my map.


texts: Remove non-numeric characters form getLatLng() function

texts:
I'm using Leaflet API. I'm using the Leaflet Draw plugin to select areas on the map. I need to set up a bounding box using leaflet draw to send to GeoServer.

The problem is when I use the following code:

var coords = the_layer.getLatLngs()[0];

I get the coordinates back like this:

LatLng(43.13306, -112.71973),LatLng(45.02695, -112.71973),LatLng(45.02695, -106.3916),LatLng(43.13306, -106.3916)

I need them back in this format (minx,miny,maxx,maxy):

-109.46777343750001,47.96050238891509,-108.58886718750001,48.16608541901253

So I need all of the non-numeric characters (excluding the commas) removed. I also need all spaces removed.

Is this possible?


texts:
I am trying to export a point shapefile to an Excel sheet using the

arcpy.TableToExcel_conversion(FinalClipOutput, ExcelOut, "", "")

The "FinalClipOutput" is the result of a clip that happens earlier in the in process and the output shape is produced.

However when it tried to export that shape to Excel table I receive the following error:

  ERROR 000732: Input Table: Dataset
  C:\Users\AFoddy\Desktop\Working\mycalc does not exist or is not
  supported ERROR 000814: Invalid file type Failed to execute
  (TableToExcel).

If I take the "mycalc" shapefile and use Table to Excel tool separately it works fine.


texts:
So, i had this very weird bug while treating my dem data lately. I used SAGA to create flows channel networks on a dem i acquired a little while ago. Preprocessed it using sink route detection and filling the sinks, then used the "channel network and drainage basins" tool. Everything worked perfectly, i got the data i needed on the area i'm working on.

Then 2 weeks ago, i get some new data covering the same area but with a higher resolution, went from pixels of 25x25 meters to pixels of 5 by 5. The data here was way bigger, but i used the same tools and preprocess. It was long but i got my preprocessed data, and after that everytime i tried using the channel network tool SAGA crashes at the "edge detection" part. Tried using the fill sinks xxl for preprocess, didn't work. Used tiles to divide my data and then applied the process to the samples, worked but what i want is doing this on the whole area in order to have only one vector file with everything i need and the right Strahler orders. Filling in the gaps manually, considering the size of the data and the area afflicted by the discontinuities caused by the tiles, is clearly out of the question.

The original data is over 1 go, i think that might be the issue and what bothers me is that i can't find a fix anywhere, did anyone run into the same issue? (From what i gathered the crash might be related to the ram SAGA is trying to use for the process, attributing to it more ressources than what the computer has.)

EDIT: I did find a way to make the data a little easier to complete and fill by making sure the gaps were as little as possible. For this i used the dem with the 25x25 pixels and extracted the watershed basins from it, then i used thoses to divide my high resolution dem into big basins that i treated on their own afterwards. Stacked the whole vector datas i obtained, added a few lines here and there to make sure that the main lines don't have any gaps and done. It is a lot better than before but obtaining the right Strahler orders and the real basins seems still a little tricky to obtain with thoses. That and using this process actualy erased some of the smaller lines on every edge of thetheorical basins i created.


texts: Vector-Viewing Properties on Multipoint Geometric Correction in Erdas Imagine 2014?

texts:
I am following some old instructions that are not up to date with my latest version of Erdas Imagine. Where can I find the Vector viewing properties within the Multipoint Geometric Viewer?

The instructions state: "The default background color for the vector road coverage is white and the default color for
the roads is black. You can change the default road color by clicking on Vector and then
Viewing Properties."

I have attached 2 screenshots. The first one is a screenshot that came within the instructions which shows the vector option in the viewer, and the second screenshot is what I see in my version of Erdas.


texts: How to create multi-scale variables from landcover data

texts:
I have calculated raster files with values of 0 and 1 for each landcover attribute I am interested in. I would like to make multiple raster files containing different cell sizes. e.g 100m2, 500m2, 1km2, 2km2 5km2 etc. that denote the coverage of each landcover attribute within each cell. How can I do this?

The purpose is for multi-scale distribution modelling using MaxEnt.


texts:
When trying to use the 'Merge' tool in ArcCatalog [EDIT (v. 10.3.1)], I am not seeing all of the merge rules- the options I can select are different than those outlined in the tool's 'About' section. Can anybody tell me why this might be? In particular I would like to use the 'Join' rule to create an output merging individual records with corresponding field values without just appending one table onto the end of another as is currently happening.

In the image below, you can see from the dropdown merge rules which are available in contrast to those that the tool is telling me should be available on the right. The order is different and 'Range' appears to have replaced 'Join'. Do these tools perform the same function?


texts: How to set the tolerance value in pgr_createTopology function of pgrouting?

texts:
I am making a map based on the floor plan of a building.
I made the shapefile and added it to the PostgreSQL database.

I created the topology for the map using the following code:

SELECT pgr_createTopology('ways', 0.00001, 'the_geom', 'gid');

The tolerance provided here is 0.00001. I didn't get many vertices when I checked it up (just 14). How can I set a proper tolerance value? What does it depend on?


texts: Speed up Rendering of Feature Layer Using ArcGIS JavaScript API 4.0

texts:
Are there any tricks to speed up the rendering of a feature layer in the ArcGIS JavaScript API?  Currently one of the layers I am using is:

This takes about 10 seconds to draw but I have used this layer which draws much faster:
Is it just a function of the speed of the ESRI server I am using compared to the "sampleserver6" which the other layer is using?


texts: Limit users' editing privileges by attribute and role

texts:
I have a feature class in an Oracle 12c geodatabase.

The feature class has a type column. The column has a coded domain.

I would like to limit user's ability to edit the features based on the type of feature. I would ideally do this based on their privilege/role in the database.

Example:

Where type = SIDEWALK, only users in the ROLE_EDITOR_GROUP_A can edit
Where type = TRAIL,    only users in the ROLE_EDITOR_GROUP_B can edit
Where type = CYCLING,  only users in the ROLE_EDITOR_GROUP_C can edit

How can I achieve this level of control?

Note:

Due to the dependency that other enterprise systems have on the data, it is our preference to not split out the different types of features into separate feature classes. We feel that single feature class would have these benefits:

Simpler and more robust than unioning together separate feature classes to get all of the records into a single view
Easier to enforce unique IDs
Easier for editors to avoid duplicate or overlapping geometry
Simpler in general to integrate with other critical enterprise systems


texts: Google Earth Engine - qualityMosaic() opposite

texts:
Is there a way to calculate the opposite of what the qualityMosaic() does to an ImageCollection?
That is to minimize a selected column (for example an added layer "ndvi") but keep the rest of the original bands associated with that ndvi (all from the same DATE).
Input= ImageCollection, Output= Single Image.


texts: How to write all attributes when multiple source files have the same name?

texts:
I am using FME to write all the MapInfo TAB files in a directory and it's subdirectories to a File Geodatabase.  I'm using dynamic readers and dynamic writers as I don't know the geometry or schema of the input files beforehand.  So far so good.  Now, suppose that I have a file called "Pipeline.TAB" at the root of the directory and I have another file called "Pipeline.TAB" in a subdirectory.  No problem.  FME can merge the two files together and I get the geometry from both files in my new feature class called "Pipline_line".

My issue is, suppose those two original files have different attributes.  Say, the file in the root folder has a field called "Name" and the second file has a field called "Pipe_Name".  My resulting feature class gets the "Name" field but not the "Pipe_Name" field.  The features that came from the file in the subdirectory have a NULL value in the "Name" field.  Instead, what I would like is to have both fields come through.  I'm happy if the "Name" field is NULL as long as the "Pipe_Name" field is there and has the correct values.

It appears that the schema is determined by the first feature that is being read, whereas, I need it to use the schema from all of the features with the same name.  Is this possible?

EDIT

I have posted this question on FMEPedia as well.  Here is some additional information which may help:

The example I gave was simple in that it only had two files called Pipeline.TAB.  In reality, I'm looking at dozens of folders with thousands of files.  there might be 20 files called "Pipeline" and another twenty called "Valves", etc.  I won't know the schema of all of them at the start.

What I would like the final attribute table of my destination feature class to look like would be something like this:

Name       | Pipe_Name
Pipe1      | NULL
Pipe2      | NULL
Pipe3      | NULL
NULL       | Pipe 4
NULL       | Pipe 5

In the final file, there may actually be many fields.  But I won't know all of them  at the start:

Name      | Pipe_Name  | Diameter | Material | Type
Pipe1     | NULL       | 500      | Steel    | Gas
Pipe2     | NULL       | 1200     | Steel    | Gas
Pipe3     | NULL       | 800      | PVC      | Water
NULL      | Pipe 4     | NULL     |NULL      | Water
NULL      | Pipe 5     | NULL     |NULL      | Water

In response to @BradNesom's answer, I have looked at the SchemaMapper, but I have way too many files to look at for that to be feasible.  I also have looked at the other link you mentioned and I tried adding a reader as a resource but that didn't seem to work either.


texts: What causes Acrobat "An error exists on this page" with PDF maps from ArcGIS 10.3?

texts:
After exporting a map as PDF in ArcGIS 10.3, the following error shows up in Adobe Acrobat:

  "An error exists on this page. Acrobat may not display the page correctly. Please contact the person who created the PDF document to correct the problem."

In the map, there is as well a lack of some features or at least feature parts.

Creating PDFs is a very basic requirement.

Do you know any patch or workaround for this, or do I have to go a version back to 10.2?


texts:
I'm trying to find where multiple lines intersect a polygon for two different geodataframes:

from shapely.geometry import Polygon, LineString
import geopandas as gpd

polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)])
line1 = LineString([(0.5, 0.5), (0.7, 0.7)])
line2 = LineString([(0.9, 0.9), (0.2, 0.6)])

poly_gdf = gpd.GeoDataFrame(geometry=[polygon])
line_gdf = gpd.GeoDataFrame(geometry=[line1, line2])

This is what the above geodataframes look like (one has a polygon and the other has two lines). It looks to me as if both lines intersect the polygon:

However, the intersect output is very confusing:

print(line_gdf.intersects(poly_gdf))

    0 True

    1 False

print(line1.intersects(polygon))
print(line2.intersects(polygon))

    True

    True

Why does the geopandas intersect method give a different output to the standard shapely one?

I am using Python 3.5.3 and Geopandas 0.2.1 all on Anaconda.


texts: Invalid Variant Type Issue - Saving from ArcGIS FeatureServer Layer to Shapefile in QGIS

texts:
I'm using QGIS 2.18.10, and I have an ArcGIS server feature layer pulled up that I would like to save locally to my hard drive for editing purposes.

The remote layer with which I am working is the polygon layer available at https://mapservices.nps.gov/arcgis/rest/services/cultural_resources/nrhp_locations/MapServer.

When I right click on the layer in the Layers Panel in QGIS and try to save it as a shapefile, I get an error stating the following:

  Export to vector file failed.  Error: Feature write errors: Invalid
  variant type for field Shape_Area[34]: received QByteArray with
  type......

The Shape_Area field in the source layer is formatted as double. I'm having trouble understanding what's causing this issue.

Can anyone shed some light on this?


texts: How to get the color of a vector layer in pyQGIS?

texts:
I have see this topic (How to change the color of a vector layer in pyQGIS?) but I'm looking for get the color not to change.

How I can resolve my problem ?


texts: Converting ArcGIS Server JSON to GeoJSON?

texts:
I'm creating a web map using Leaflet, and I want to be able to grab feature layers from our ArcServer. I have successfully been able to retrieve a feature class as JSON, but Esri's JSON objects don't follow the GeoJSON standards so they cannot be displayed.

Does anyone know of a conversion script or tool that handles this?

If not, I plan on creating a script to convert ArcServer JSON objects to GeoJSON.


texts: Do older community created widgets for 1.x work with CMV 2.x

texts: I tried adding the "Nearby" widget from the "http://github.com/BrianBunker/cmv-widgets/tree/master//Nearby" user contributions. This appears to be an older widget and doesn't seem to work with CMV 2.0 when JUST editing the viewer.js file for this widget. Is there anything special I am missing? I am new to CMV and only a general hack JS programmer.


texts: GeoKettle - inserting geometry into a table

texts:
I've got some shapefiles, and need to convert the shapes into an SQL Geometry type and insert them into a table in our database. We're using GeoKettle.

I can open the shapefile find and read values in from it into the table (such as a text field in the shapefile into a varchar in the table), but don't know the steps I need to translate the shape, which appears to be coming out of the shapefile as WKT, into a geometry.

Basically, whenever I execute, I get the error that it doesn't seem to be able to insert it. Which is understandable, as I say it seems to be coming out as WKT and trying to insert that into a geometry field. And I know the method to convert it to geometry in SQL (STGeomFromText) - I just essentially need to know how to add that in GeoKettle.


texts: Right-align attribute aliases in QGIS attribute forms

texts:
I have an attribute form created with QGIS 'Drag and drop Designer' that looks like this:

As can be seen, the attribute aliases (the labels on the left) are little longer because of some explanatory text in the round brackets.

It might be a matter of taste and a matter of readability as well, but is there a way to right-align theses labels to make them stand near the related input widet? Can this be achieved with QGIS standard tools or do I have to create an ui file with the QtDesigner?


texts: upgrade from 1.3 to 3.3 javascript api for arcgis with jsviewer

texts: I try to use jsviewer (v.1.0) in my application, but it doesn't work with JSAPI version higher than 1.3. How can I  work with both jsviewer and the 3.3 JSAPI version?


texts:
I'm trying to implement in my plugin such features as "select by freehand" and "select by polygon" for some purposes, and I'm trying to search internet for how they work.

But all I can find is:

QAction *QgisAppInterface::actionSelectFreehand() {
      return qgis->actionSelectFreehand();
}

I need to know how features are selected under rubberband which is created by this function(s).

It's certainly not  QgsGeometry.intersect, or QgsGeometry.within, or QgisVectorLayer.select()

Is there any secret knowledge how this native built-in functions work?

Doesn't matter if it is cpp or Python documentation.


texts: Can I sell maps produced from ASTER Global DEM V2 data?

texts:
I don't quite understand the term-of-use of ASTER Global DEM V2 data:

  By selecting 'Submit', I agree to the following:

  -I agree to redistribute the ASTER GDEM only to individuals within my organization or project of intended use or in response to disasters in support of the GEO Disaster Theme.

  -When presenting or publishing ASTER GDEM data, I agree to include 'ASTER GDEM is a product of METI and NASA.'

Does it mean I can't sell maps produced using this data?


texts:
Using localhost page generated by osm2po roteamente noticed that the link is made to link even disabling snapping.

Logically thought to create more nodes in the network, but how can I do this without losing the topology of the data?


texts: Copy multiple raster filenames in a single CSV file in R

texts:
Can somebody help me on how to copy multiple filenames in a single CSV file? I have 365 raster files and a polygon dataset with 5 features. What I am trying to do is extract cell values corresponding to polygon features. I am able to do this but not able to copy the corresponding raster filenames in a single CSV. Below is my code:

library(rgdal)
library(raster)

setwd("C:/PHP_Win7/Etios/2007")
polygong <- readOGR("C:/PHP_Win7/Etios/2007", "Ethn")
rasterfiles <- list.files(".", pattern = "tif")

for (j in 1:length(rasterfiles)) {
  curRaster <- raster(paste("C:/PHP_Win7/Etios/2007", rasterfiles[j], sep = "/"))
  csvDatatoWrite <- data.frame(coordinates(polygong),
                               polygong$Name,
                               extract(curRaster, polygong))
  write.table(csvDatatoWrite, file = "output.csv",
              append = TRUE, sep = ",", row.names = F)
}

The output I am getting is

I would like to replace "extract.curRaster..polygong." column heading in output sample with corresponding filenames. Any help on this?


texts: Parser for NOTAM into GIS point objects whith attributes

texts:
I am participating in a gis application development project. There is a task of planning from a drone flight application, taking into account the limitations of NOTAM. Everything would be fine if it were not for the regular NOTAM updates every hour and the need to display them.
Perhaps there is actually a ready NOTAM parser of telegrams for translating them into some GIS format, for example shp or geojson.
An example of the source file from the links:
http://www.caiga.ru/common/AirInter/series_notam/A2005281855_eng.html
http://www.caiga.ru/common/AirInter/series_notam/A2005281855_rus.html


texts:
I'm downloading Sentinel 3 and 5p data using the dhusget.sh script provided here:

https://scihub.copernicus.eu/twiki/do/view/SciHubUserGuide/BatchScripting?redirectedfrom=SciHubUserGuide.8BatchScripting#Query_via_wget

While it works, it is very slow, and downloads data at roughly 1mbps on my 100mbps network. Is this a solvable problem, or is it like this for everyone?


texts: Picking up the coordinates inside this polygon draw

texts:
I already have a geojson being consulted and that is plotting these black dots. And then I inserted the leaflet-geoman plugin to give users the freedom to create polygons on the map. Basically draw something there.

For example: Draw a polygon to get only the information inside that polygon.
When using the tool to draw a polygon on the map I wanted to get which black dot is inside the drawing

in short: I have 4 points within my polygon drawing (which does not come from a geojson, the person draws) and I want to get the information from those points that are within that drawn polygon. Each black dot in this one has an id, and I want to know what points are within the space that the person drew.
My geoman-Leaflet code:
      map.pm.addControls({
        position: 'topright',
        // drawCircle: false,
        drawMarker: false,
        drawCircleMarker: false,
        drawRectangle: false,
        cutPolygon: false,
        editPolygon: false,
      });

My black dots on the map code:
  wellheads.forEach(function (d) {
        d.LatLng = new L.latLng(d.latitude, d.longitude);
      });

  const groupWells = g
        .selectAll('g')
        .data(wellheads)
        .enter()
        .append('g')
        .attr('pointer-events', 'visible');

      const points = groupWells
        .append('circle')
        .attr('class', d => {
          if (wellState.selectedWell.guid === d.wellboreGUID) {
            return 'redW';
          } else {
            return 'blackW';
          }
        })
        .attr('id', d => `id-${d.wellboreGUID}`)
        .attr('r', d => {
          if (wellState.selectedWell.guid === d.wellboreGUID) {
            return 1.7;
          } else {
            return 1.0;
          }
        })


texts: I try to digitize a large number of topographic maps. Digitizes the contour every four meters. I use the program ArcGIS 10.3 and ArcSCAN. The problem is that it has many objects in the map outside of the contour and out of that the biggest problem is that the contours are continuous but is in many pieces. Is there a way to integrate automatically and in which way? It can be solved only by the topology? there is an automatic way to come together (Merge but only the same - meters iso lines)?


texts: Openlayers 4 - Determine if layer is a layer group

texts: Is there a method or property in javascript using the Openlayers (3 or) 4 library which I can use to verify if a given layer is a layer group?


texts: Get waterbody info from base map in Leaflet

texts:
I am plotting points on a map using Leaflet, and would like to collect any data I can concerning the water bodies that encompass these points.  I have been making use of polygons from HydroSHEDS (https://hydrosheds.org) to identify waterbody outlines, but this is not a complete list of waterbodies (The default base map used in leaflet seems to have polygon info on many more lakes).

As an example, if I make a map like this:

library(leaflet)

Lat<-c(43.847863,43.847199,43.846288,43.852605)
Lon<-c(-74.886135,-74.886111,-74.897731,-74.899169)
LakeID<-c(1,2,3,4)

leaflet() %>%
  addCircleMarkers(lat = Lat, lng = Lon,
                   color = "blue",
                   radius = 2)%>%
  addTiles()

For points that fall in a water body, I'd like to get info on the polygon used for that lake, and the name of the lake as well. I can see that the base map has the larger and smaller waterbodies colored blue, and has the name for the larger one "Twitchell Lake", so I assume that data can be queried in some way.

Is there a way I can access this data from the basemap so I could find out:

Polygon used for individual waterbodies
Name of waterbody


texts: Looking for an efficient way of digitizing building footprints using orthophotos and ArcGIS 10

texts: In the past I've used PCI Geomatica to do something similar using the softwares classification functionality. Are there some tools/tricks in ArcGIS that will partially automate the process. I remember reading somewhere about rasterizing an orthophoto and assigning values to rooftops, which polygons where then created from those values.  I can see there being issues in trying to automate the task as there will surely be some inaccuraces do to color, obstructions, image quality, etc.


texts: Extracting a time interval from NetCDF-data in ArcMap Desktop 10.4

texts:
In ArcMap 10.4 (advanced licence) I have added a netcdf-file with EU-temperature-over-time data (time being every day between 1995 to 2016).

I just need the year 2010 from this data, but I'm still struggling to extract this year. I'm new to python, and the only code I found to extract timeslices was accessing a band to extract the slice. It's from the esri support: http://support.esri.com/en/technical-article/000011318 .
But I can't manage to make a band from the time dimension of my netcdf-file, ArcGIS always crashes somewhere along the way.

I also tried by OPeNDAP-access to get just this one year of data, both downloading it directly from the website as dods, and trying to access the url, directly through the ArcMap tool (make OpENDAP-layer). The download page ( http://opendap.knmi.nl/knmi/thredds/dodsC/e-obs_0.25regular/tn_0.25deg_reg_v15.0.nc.html ) provides the url in a specific slot. But the tool could not connect with the URL, it gives the "error 000241: The specified pathname or file type is invalid for the required input". And I did not find a way to open .dods with ArcMap.

Does anyone know a simple code-snippet I could have a go at editing as a beginner? Or is there a way directly in ArcMap to extract parts of the file? Ore does anyone have a solution to my OPeNDAP-Problems?


texts: Is there any thing in QGIS which can perform operations like ArcScan? I put WMS polygon layer as a base map. I save it as PNG and when I add it as a raster layer, I can vectorize it with ArcScan. I want to know if there is anything in QGIS which will function like this.


texts: Open shapefile in ArcGIS Runtime WPF SDK as Graphics/Feature layer

texts:
I Googled and followed this tutorial and figured out how to import Shape files and Rasters into Runtime.

However, it seems most of the functionality uses a FeatureLayer/GraphicLayer instead an ArcGISLocalDynamicMapServiceLayer that the sample uses.

My question is if there is a way I can import the shapefile as a graphics layer or if I can somehow cast my ArcGISLocalDynamicMapServiceLayer as a Graphics layer because shapefile is essentially the same as a shapefile.

An example is I'm trying to use the FeatureDataGrid() to show the attribute table of my shapefile, but I can't since it requires a graphics layer, not an ArcGISLocalDynamicMapServiceLayer.


texts: I have survey data for contacts listed as eastings and northings - I need to overlay this onto PNG files but the location is quite far north - How do I apply EPSG 1612 transformation to the UTM Zone 32N so they appear as layers in the same place?


texts:
I developped a plugin on QGIS 2.18, and used qwt5 to implement an elevation graph, which works quite well.
Now that a want to migrate my plugin to QGIS3, and knowing that qwt5 is linked to PyQt4, I'm not able to find which module to use in my Python code, or even to know if qwt (6 or else) is still available in QGIS3. In OSGEO4W setup installation, it still seems to be available, but I don't know how to call it in my code.
Answer :
I finally found the python3-qwt packages in OSGeo4W Setup. Cannot say if it's new or if I missed it before.


texts: I have two races in the same town that make use of the same stretches of roads/paths in places (sometimes in the same direction, sometimes in opposing directions).  Given that I know the pace of the fastest and slowest competitor in each race and the start time of each race, is there a QGIS tool that would allow me establish (and show) where the two groups of competitors would be sharing the same stretch of road/path?  That is where a conflict might arise.


texts:
Will the ESRI  ModelBuilder 10.1 iterator select a value from field in a filegeodatabase table?
If so then can this be joined to a feature class and exported or copied?

I am working with property appraiser data and for example  one parcel of land may sell up to 10 times in a year, which is in a table, the parcels are features.  I would like to select, say all the parcels that have sold twice in one year and then join to the a parcel map.

So I would like to use Iterator to select those sales from a table,  and then join to a parcel feature, copy features then move on to next group of Sale, say Sales 3.
I do this somewhat manually know with modelbuilder.


texts: Using elif statements in Python parser of ArcGIS Field Calculator?

texts:
Looking for some help. I have a continuous variable (PCT_B100P) of the percentage of folks in a census tract that are under the 100% federal poverty line. I'm trying to use this variable to assign a categorical variable (POVCAT)to the census tract. If the value is below, between, or above certain limits it will determine what category the record gets. I'm not getting any errors in my code, buy it's only calculating the first category level 1, and nothing else. Here is my code.

def NumValue( PCT_B100P ):

  if (PCT_B100P <5):
    return 1
  elif (PCT_B100P >= 5 and PCT_B100P <10):
    return 2
  elif (PCT_B100P >= 10 and PCT_B100P <20):
    return 3
  elif (PCT_B100P > 20):
    return 4

Expression:

NumValue(!POVCAT!)

any idea as to what I'm doing wrong? I'm quite new to Python. Below is an image of my issue


texts: Loading a web map from ArcGIS Portal on-premise

texts:
In my application, I have a few navigation tabs.  Each tab loads a web map from Portal for ArcGIS on-premise.

When each tab is clicked, it asks for the Portal credentials in order to load the maps.

I want end users to log in once on the main page and be able to access the maps without logging in over and over again.  Are there any ways do so or other alternative ways?


texts: ArcPy reclassify with if statement and cell value

texts:
I have a workspace with NDVI rasters and another workspace with greyscale rasters. The rasters have the same name.

I want to reclassify the NDVI raster while also using the greyscale values in a function. Depending on the greyscale value, another reclass function should be used. So for every pixel, the correct function has to be found.

I need to use an if statement: if the greyscale raster cell value is high, another statement should be used but the code results in an error:

  Invalid use of raster with Boolean operator or function. Check the use of parentheses where applicable.

arcpy.env.workspace = ndviws
ndvirasterlist = arcpy.ListRasters()

for rimage in ndvirasterlist:
    ndvi_raster = arcpy.sa.Raster(rimage)
    grey_raster = arcpy.sa.Raster(os.path.join(greyvalues,rimage))
    if grey_raster <= 28523:
        output_raster = Con(ndvi_raster < -0.0000006 * grey_raster +  0.3448,2,1)
    if grey_raster >28523 and grey_raster < 50000:
        output_raster = Con(ndvi_raster < -0.000009 * grey_raster + 0.437673,2,1)
    if grey_raster >=50000:
        output_raster = Con(ndvi_raster < -0.000007 * grey_raster + 0.5233,2,1)
    output_raster.save(os.path.join(outputws,rimage))

edit: **edited so that there are now 3 if statements. Is this still working with con method?


texts: Convert longitude/latitude to pixels on the map

texts:
I have a map from here. I want to be able to simply convert arbitrary lon/lat pair to pixels on the map (also the ability to do reverse conversion). The maps come with .tfw file and projection info, here it is:

0.02222222222222
0.00000000000000
0.00000000000000
-0.02222222222222
-180.00000000000000
90.00000000000000

and projection info:

Projection: Plate Carree aka Geographic or "LatLong"
Earth ellipsoid: Sphere, radius 6370997 m
Datum: WGS84
Extent: 180 West to 180 East, 90 North to 90 South
Size: 16,200 height samples wide x 8,100 high
Resolution: 2.47 km/pixel

I am completely new to cartography-related stuff, and as far as I figured out I should first do the transformation from WGS84 (lon/lat pair) to Geographic projection (aren't they the same?). It seems to me that they are the same, actually, but sphere radius in projection info above is 6370997 and is different from spatialreference.com page I found for Plate Carree projection.
Anyway, I found DotSpatial.Projections library can do this for me, with the following code:

    var kievCoordinates = new[] { 50.4546600, 30.5238000 };
    var z = new[] { 1.0 };

    var wgs84 = KnownCoordinateSystems.Geographic.World.WGS1984;
    var dest = new ProjectionInfo();
    dest.ParseEsriString(
        "PROJCS[\"WGS 84 / Plate Carree (deprecated)\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]],PROJECTION[\"Equidistant_Cylindrical\"],PARAMETER[\"central_meridian\",0],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1]]");

    Reproject.ReprojectPoints(kievCoordinates, z, wgs84, dest, 0, 1);

And then I must translate resulting coordinates to pixels on the map using world file. I am aware of the following formula:

But it seems that there are degrees not meters in the world file and I don't know what to do with them.
Generally, am I doing right things? Or there is more easy way, given my data?


texts: Python Rasterio not recognized as a supported file format

texts:
I'm trying to run this python code with anaconda.

I have an error that I can not solve.

I try execute this tutorial : https://towardsdatascience.com/satellite-imagery-access-and-analysis-in-python-jupyter-notebooks-387971ece84b

to work with sentinel-2's image.

At the ligne:

b4 = rio.open(R10+'/T31TFJ_20190905T104029_B04_10m.jp2')

I have this error :

RasterioIOError: 'image/S2B_MSIL2A_20190905T104029_N0213_R008_T31TFJ_20190905T135151.SAFE/GRANULE/L2A_T31TFJ_A013047_20190905T104455/IMG_DATA/R10m/T31TFJ_20190905T104029_B04_10m.jp2' not recognized as a supported file format.

I remove gdal with command conda remove gdal and re-install with conda install gdal

Any ideas ?

EDIT: I do every test with GDAL. I Have a error...

from osgeo import gdal

gdal.UseExceptions()
ds = gdal.Open('test.jp2')

This error : `test.jp2' not recognized as a supported file format.
Drivers gdal problems ?

other thing, when I installed Gdal with conda, every conflict found. There is already gdal installed on my computer for Qgis (mapping software). Can this do to problems?

Thanxs.


texts: Where can I find ogr_FileGDB.dll v1.4 compilled for Windows?

texts: Since v1.3 do not have custom CRS support, where can I find the ogr_FileGDB.dll v1.4 compiled for Windows?


texts: Georeference disabled in CartoDB synced dataset?

texts:
If I upload a dataset with synced option enable then I haven't the option to do the geocoding.

Is that normal?

How I can geocode the addresses?

If manually I can't do the geocoding, how I can do programmatically?


texts: Shapefile to KML Conversion Problems in QGIS 2.6.1 on Mac

texts:
I am currently using QGIS 2.6.1 on a Mac with OS 10.9.5.  I am trying to export a shapefile using QGIS to KML format.  I am having success with the export but the data is skewed halfway across the world and is wrong.  Is anyone else having this problem?

Here is my data before the export in shapefile format (WGS84)

Now here is my data after the shapefile to kml conversion

These lines wrap halfway across the world


texts: How can I confirm the existence of pyramids in QGis?

texts:
I have created pyramid layers for my ASTER satellite images (.img) using Gdal Tools "Build Overviews" (Pyramids). The dialog said "Processing completed". I don't have the impression that zooming in or out has been accelerated in any way and under "Properties" / Pyramids there is still the "description", no pyramids are listed.

Does anybody know how I can find out about this?


texts:
I am working with hydraulic network. I have a pipes network. I am trying to check it's geometries.
I have searched many related documents, I know topolygon tool and check geometries, but it only helps to detect error in the case one, the result is duplicate pipes.

Is there any methods to detect error in the case two?


texts: Missing geometry type and SRID in PostGIS

texts:
I created a view with
select a, st_union(geom)
from table
group by a

srid: 25832, geometry type: Multipolygon
The view has no SRID and the geometry type is "GEOMETRY".
Now i want to calculate the area with the geometry in the view. Because of the missing SRID the results are not correct. How can i set the SRID and geometry type for the gemoetries for the view?


texts: How to properly import corine landcover 2012 (vector file) in QGIS?

texts:
I'm fairly new to QGIS (using version 2.18.25), so I may be missing something quite basic. I'm using the Corine landcover (clc12_Version_18_5.gdb) in my project. It's a very large dataset (3GB) so it takes time to render but does eventually in default single colour.

However when I add the legend entries (properties > style > load style > clc_legend) associated with the landcover types the layer fails to render.

Do I need to complete another process for this to work or is something else wrong?


texts: OpenGeo SDK: Add New Server trouble with proxy/?url=

texts:
I'm currently developing an application with the OpenGeo SDK.
I wish to have the "Add New Server" functionality and I've followed the example from the documentation: http://suite.opengeo.org/opengeo-docs/webapps/sdkintro/viewer/addnewserver.html

In my applications file, I've added proxy: "proxy/?url=" and also * @require overrides/override-ext-ajax.js

When I test the application on my server (e.g., suite-sdk debug path/to/myapp) everything works fine.

However, when I deploy the application (I use a Tomcat container) then the "Add New Server" functionalities doesn't work any more. I've always the following error message when I try to add a WMS: Error getting WMS capabilities (Layer source not available.). Please check the URL and try again.

I don't know what I'm doing wrong.


texts: Changing label for record counter in Survey123

texts:
I have a survey where I am collecting nested repeating questions.  I have 3 groups; Stands, Plots, and Trees.  Stands have many Plots.  Plots have many trees.
The display of the records is a little bit confusing.  Instead of showing "1 of 1", I would like the record labels to say "Stand 1 of 1", "Plot 1 of 1" and "Tree 1 of 1".
Is there any way to accomplish this?


texts: Understanding difference between Polygon and Multipolygon for shapefiles in QGIS?

texts:
Creating Scratch layer gives me a choice between Polygon and Multipolygon, so I started wondering how should I make use of them.

Reading a new book "Mastering QGIS" (through this answer by Joseph to What does buffer mean in QGIS Voronoi dialog?), authors' comment was "ESRI shapefiles store geospatial data in multi- type geometry, so using multi- type geometry is a good practice if you plan on converting between formats."

How should I choose one of them (Polygon or Multipolygon) when I create a new polygon layer?
(I had confused multi-polygon with multipart-polygon, so I used to choose Polygon.)

One recent situation I have come to choose Polygon/Multipolygon was Changing parallel line length to split polygon in QGIS? in which I draw a dummy polygon as a background layer. Working mainly on shapefiles it seemed I was supposed to choose Multipolygon, although I needed only one polygon. In this case what I needed was simple polygon geometry without any attribute (Not sure this example make any distinction.)


texts: Atmospheric correction algorithm in Python or R for Landsat 7 ETM+

texts:
I am looking for an algorithm which can be tuned in python or R for atmospheric correction for Landsat 7 ETM.

Searching Google I have found only a few algorithms in ENVI or ERDAS which are not possible to alter the parameters or maybe I am wrong. Can anyone help me to find an algorithm sufficient for Landsat 7 in order to perform atmospheric correction?


texts:
Having used QGIS for quite some time for mapping, I was 'forced' to dig into GRASS as pansharpening is not directly executable in QGIS. After understanding more or less the structure, I did some trials with i.fusion.brovey. I was able to generate the rasters with 15m resolution, but when loading into QGIS I can not achieve/recuperate the nice colour contrast  as is present in the original 30m rasters.

In an earlier post someone confirmed that and proposed to use another function for landsat: 'i.his.rgb' (after using i.landsat.toar, which I don't use). He starts to convert min/max to 1-255 before using i.his.rgb...My original rasters have a range of 0-65535. It did not work out for me...

I quit after some days of struggling to just perform a pansharp, I feel stuck... Why should just applying i.pansharp not work for Landsat as suggested by others? Can someone suggest a more straightforward workflow to execute pansharpening of Landsat 8 images? Looking forward to some advice, thanks in advance...


texts: Surface analysis with R - Shape recognition

texts:
I am trying to analyse a geological surface with R & Python and I would like to know if there are some ways to distinguish, divide and quantify repetitive structures such as advertised by a software company (picture below).

I don't have a precise scope but I would like to understand what shapes are repeating themselves onto the analysed surface.

Are there advised packages that could be useful to accomplish this task?


texts: Calculating area of a lat,lon polygon in some useful units

texts:
How do I tell the default geo spatial reference system? I have this code:
wkt = "POLYGON ((lon1 lat1, lon2 lat2, lon3 lat3))"
poly = ogr.CreateGeometryFromWkt(wkt)
print "Area = %d" % poly.GetArea()

I run it and get back a number.  But what are the units? The documentation for OGRGeometry's get_Area() method says:

the area of the feature in square units of the spatial reference system in use.

So.. how do I tell what spatial reference system is currently in use?
Another way to formulate my question is:
I have a polygon defined as a sequence of lon,lat pairs. How do I get its area in sq kms?
So then I say: probably when people talk about areas of polygons on Earth they don't mean integrating the actual surface integral but area of some projection.  So let me try setting the projection I want.  And I got this code:
wkt = "POLYGON ((lon1 lat1, lon2 lat2, lon3 lat3))"
poly = ogr.CreateGeometryFromWkt(wkt)
sr = osr.SpatialReference()
sr.SetWellKnownGeogCS("NAD83")
sr.SetProjCS("SRS_PT_TRANSVERSE_MERCATOR")
poly.AssignSpatialReference(sr)
poly.GetArea()

but I get the same number as before!
Printing out sr object gives me plausible looking output (it says "metre"!)
PROJCS["SRS_PT_TRANSVERSE_MERCATOR",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["unnamed"],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]

So I tried setting source and target systems and creating a transformation:
poly = ogr.CreateGeometryFromWkt(wkt)
sr = osr.SpatialReference()
sr_proj = osr.SpatialReference()
sr.SetWellKnownGeogCS("EPSG:4326"))
sr_proj.ImportFromEPSG(6991)
print(sr)
print(sr_proj)
poly.AssignSpatialReference(sr)
transform = osr.CoordinateTransformation(sr, sr_proj)
poly.Transform(transform)
print("Area = %f" % poly.GetArea())

gets me square meters!  Yay.
(but the problem now is... how do I pick which EPSG to use?? I don't know where the area is in advance!)


texts: Converting coordinates to pixels with out losing points

texts:
I'm trying to convert coordinates from lat/lon into pixel but I lose points in this process.

The code that I'm using is the following:

from matplotlib import pyplot as plt
import numpy as np
import pandas as pd

df=pd.read_csv('cords.csv')
cords=df.as_matrix(columns=['x','y'])
gt=[7.6445503225, 5.4065168747250134e-06,  0.0,  45.07436634583334,  0.0, -5.406516856707135e-06]

index=np.zeros(cords.shape)
index[:,1]=((cords[:,1] - gt[3]) / gt[5]).round()
index[:,0]=((cords[:,0] - gt[0]) / gt[1]).round()
index=index.astype(int)
index[:,0]=index[:,0]-min(index[:,0])+1
index[:,1]=index[:,1]-min(index[:,1])+1
row=max(index[:,1])
col=max(index[:,0])
image=np.zeros([row+1,col+1])
for i in range(0,len(index)):
    image[index[i,1],index[i,0]]=255

If I plot the cords or the index points I get this:

If I plot the image I get this:

As you can see there are some points that are missing in translating lat/lon in to pixel numbers. Yellow is 255 value and purple is 0 value.
How can this be solved?

Here you find the coordinates that I'm using cords.csv

Here you find the coordinates with the values that need to be set to each pixel. cords_valus.csv

EDIT

In the folowing link you find the raster image from where I have extract the cordinates. It fits to Italy perfectly. raster


texts: Create new custom projections in QGIS 3.12?

texts:
QGIS does not have the coordinate reference system I need, So I decide to define a new custom Geographic coordinate system. To define a CRS, I selected customProjection Custom CRS from the Settings menu. The Custom Coordinate Reference System Definition dialog requires only two parameters to define a user CRS:
1) A descriptive name
2) The cartographic parameters in PROJ.4 format

So I did that as you can see from the picture:

But when I click OK to add the projection within QGIS a pop up message with an error show up and with very different set of parameters in the parameter window. the error message says I have set up a proj string definition equivalent to esri:102218 which is not really the same. I don't really know if this a bug or something went wrong with my parameter setting.


texts:
How can I programmatically determine if a shapefile contains points, polygons, or lines?

ArcEngine 10, VS2010, C#


texts: Losing printer settings when saving copy of mxd from ArcPy?

texts:
Judging from the answer to Changing page layout from portrait to landscape using ArcPy with ArcMap?, one cannot change the page size etc from ArcPy, so I set up my "template" mxds so my script can save a new copy of it, change data sources etc. This methods works perfectly for my A4 maps.

However, my A0 mxds are losing the printer and page settings from the template. My A0_template.mxd is set to print from the plotter on A0 paper using the default printer settings. After the copy of the mxd is produced and I open it, the printer has somehow changed to the normal small printer with the paper size matching the default tray. I am simply using the saveACopy(filename) method which is part of the MapDocument class.

Why would the correct settings for the A4 mxd be saved, but not for the A0?


texts: Creating polyline using mouse clicks with Python Addin tool?

texts:
Can somebody tell me how to create a line geometry by mouse click with python addin?

I am using below code in python addin tool.

class ToolClass8(object):
"""Implementation for Testmsg_addin.tool (Tool)"""
def __init__(self):
    self.enabled = True
    self.cursor = 3
    self.shape = "Line"

def onLine(self, line_geometry):
    NewL = "NewL"
    cur = arcpy.da.InsertCursor(NewL, ["SHAPE@"])
    array = arcpy.Array()
    part = line_geometry.getPart(0)
    for pt in part:
        array.add(pt)

    cur.insertRow(arcpy.Polyline(array))

which throws the error:

cur.insertRow(arcpy.Polyline(array))
TypeError: argument must be sequence of values

Can somebody help me with the code?


texts: While using HEC-RAS 4.0 Beta, how to set PF1 in Steady Data Flow?


texts: Recovering projection information after reduction Google Earth Engine

texts:
Is there a way to reproject the resulting image after a reduction so that we can use the geometry of the image for another operation?

I am using image.reproject() but when calculating the image.geometry() I get the following error:

Here is the link to the script: https://code.earthengine.google.com/1435eedcfff6d9a84e305df2267e2e26

And the script:

var geometry = ee.Geometry.Polygon([[[-93.68530048461537,37.399409695581674],
[-93.69628681274037,36.781581043994294],
[-92.76244892211537,36.78598052408677],
[-92.78442157836537,37.40813695578065],
[-93.68530048461537,37.399409695581674]]])

var landsat = ee.ImageCollection('LANDSAT/LC8_L1T')
  .filterDate('2016-01-01', '2016-12-30')
  .filterBounds(geometry)

//Compute NDVI
var ndvi_images = landsat.map(function(image){
  return image.normalizedDifference(['B5','B4']).rename('ndvi')
})

//Calculate max NDVI
var maxNDVI = ndvi_images.max()

//Reproject
var reprojected = maxNDVI.reproject('EPSG:32615', null, null)

//Remove edges
var clipped_image = reprojected.clip(reprojected.geometry().buffer(-6000))

print(clipped_image)


texts: How to use the ESRI JS API to map SQL Server Spatial Types without ArcGIS Server

texts:
I have heard of people using the ESRI Javascript API to create webmaps with spatial data from SQL Server or other RDBMS's.  From what I can tell, they are not using ArcGIS Server as a middleware.  In particular, I recently heard of a company that uses a SQL Server database with spatial types and Entity Framework to link the database to the app, but then they are using the ESRI JS API on the front end to render the data.

This is something that really interests me as I have worked with SQL Server, Entity Framework, ASP.Net, and the JS API separately and I do not have access to ArcGIS Server.  Does anyone know of a particular architecture or workflow for simply using the JS API as a front end for a non ESRI backend?


texts: Creating web scene from ArcGIS Pro to ArcGIS Online?

texts:
I am trying to create a web scene from ArcGIS Pro to ArcGIS Online. However, I am running into a lot of issues, the most annoying one being that I am running into a very vague warning. Through using the Share tab to share the scene as a web scene and then analysing it, I am able to get the ERRORS to 0 and the WARNINGS to 1, but the warning is a very vague warning that reads Warning 24080: Ensure scene is authored correctly before sharing. The solution points to how to author a web scene, which gives some things to avoid doing. I went through very carefully and reviewed everything step there, and even watered down my scene to something that isn't very useful for my purposes, but just to get something to actually go through, and still no luck. When I try running any web scene to ArcGIS Online from ArcGIS Pro the process will fail. On ArcGIS Online the WFL will show up, but it will be empty. I have no idea what is going wrong with this.

Ultimately I would like to bring a 3D web scene into ArcGIS Online from ArcGIS Pro, but I can't even get a very basic 2D scene to publish. Has anyone found a solution to this problem? Currently I have to create KML features to get anything over to ArcGIS Online, and even then, it is not recognized in the ArcGIS Online scene viewer. The whole intent of the project I am working on is to produce 3D models on ArcGIS Online to be able to share it with a larger audience than keeping it on ArcGIS Pro on the desktop could provide. I have created some web pages using WebGL, ArcGIS v3, ArcGIS v4 Beta, Cesium, and Google Earth APIs with JavaScipt to see how well one of these workarounds would do, but they are a lot of work and they are not as easily integrated with creating story maps from ArcGIS Online in the future.

Recap: Looking for a way to bring my 3D scene in ArcGIS Pro to ArcGIS Online as a web scene. Share As Web Scene in ArcGIS Pro is not working for me, and I can't figure out why.

I am not using multipatch features, extrusions, feature classes with over 2000 features, or 3D symbology. My scene is in global mode as well. As you can see my scene is not very advanced, and even though I want 3D extrusions of the buildings I am not even able to get a scene without 3D Properties to publish.


texts: Python error when batching a Graphical Modeler script [QGIS 3.4.9]

texts:
I have a very simple Model to convert a table in points and to extract some features by their attributes. I can run it without any problem, except if I do it as a batching process.

I received this error:

This error occurs since I upgraded my QGIS version to the last LTR.

EDIT: I've tried to change the lines in the wrappers.py file as proposed by etrimaille and I cannot find the lines in my wrappers.py file


texts:
I need to write a script that updates the "LOCALITY" field in a feature class if it is empty with the suburb it is in.

I have identified the number of properties with select by location and the below script will show that 1340 properties in the first index of the localityList will need to be updated.

PyScripter says that the script runs ok but when I check the attribute table in ArcMap the LOCALITY field has not been updated.

How can I implement the update tool correctly? After that is done all I will need to do is loop through the localityList to move through all of the suburbs.

# Import modules
import arcpy
import os

#Set workspace geodatabase
arcpy.env.workspace = "C:\Users\Harry\Documents\Uni\Semester 1 2016\GIS3407 GIS Programming and Visualisation\Assignment 2\Assignment2.gdb"

# Create a list of suburbs
# Create an empty list
localityList = []
#Create the search cursor
rows = arcpy.SearchCursor("TWB_Suburbs")
#Move to the first row
row = rows.next()
#Append the list with the selected locality name
while row:
    localityList.append(str(row.ADMINAREA))
    #Move to the next row
    row = rows.next()

print(str(len(localityList)) + " <-- # of suburbs")
del row, rows
# Make a layer of all of the properties with no locality value
# Create query
nullQuery = '"LOCALITY" = \' \''
# Create layer
nullLayer = arcpy.MakeFeatureLayer_management("TWB_Property", "TWB_Property_Layer", nullQuery)

# Select by location
indexCounter = 0
localityQuery = '"ADMINAREA" = ' + '\''+ localityList[indexCounter] + '\''
localityLayer = arcpy.MakeFeatureLayer_management("TWB_Suburbs", "TWB_Sububs_Layer", localityQuery)
arcpy.SelectLayerByLocation_management("TWB_Property_Layer", "WITHIN", "TWB_Sububs_Layer")

parcelList = []
rows = arcpy.SearchCursor("TWB_Property_Layer")
row = rows.next()

while row:
    parcelList.append(row.OBJECTID)
    row = rows.next()

print str(len(parcelList))

numberParcel = len(parcelList)

for parcel in range(0, numberParcel):
    rows = arcpy.UpdateCursor("TWB_Property", '"LOCALITY" = ' + "'" + localityList[indexCounter] + "'")
    row = rows.next()

row.setValue("LOCALITY", str(localityList[indexCounter]))
rows.next()

del row, rows

arcpy.Delete_management("TWB_Property_Layer")
arcpy.Delete_management("TWB_Sububs_Layer")


texts: How do I determine appropriate Date query syntax for a layer?

texts:
I'm writing a web mapping query tool to query Esri mapservice layers by date using the REST API.

It has been said that the syntax depends on what kind of database is behind the mapservice layer, as described in the help doc.

For oracle use this:

Datefield = date 'yyyy-mm-dd'

For SQL Server use this:

Datefield = 'yyyy-mm-dd hh:mm:ss'

How do I determine which database is being used by a particular mapservice layer so I can compose the query using the appropriate syntax?


texts:
Is QgsMapLayerRegistry going to be availabe in pyqgis with QGIS3

When I typed it into the QGIS python console (2.99), python seems to recognize it because it auto-completes qgis.core.QgsMapLayerRegistry, but then gives an error that it's not there:

Traceback (most recent call last):
  File "/usr/lib/python3.5/code.py", line 91, in runcode
    exec(code, self.locals)
  File "<input>", line 1, in <module>
AttributeError: module 'qgis.core' has no attribute 'QgsMapLayerRegistry'

I tried this with two versions of QGIS 2.99 (one in 01/07 and on in 06/17). Both versions give the same error. import doesn't work either:

>>>from qgis.core import QgsMapLayerRegistry
Traceback (most recent call last):
  File "/usr/lib/python3.5/code.py", line 91, in runcode
    exec(code, self.locals)
  File "<input>", line 1, in <module>
ImportError: cannot import name 'QgsMapLayerRegistry'

Is QgsMapLayerRegistry removed in QGIS3?


texts: Which books from Esri Press can you recommend?


texts:
I have a classified image that is roughly 400 KB and every time export the image to Google Drive it's roughly 300 KB as a cropped image. Strangely, the first time I exported to drive the image was full size. Thereafter it was cropped. Code below. Any ideas how to export the full size image consistently?

Export.image.toDrive({
  image: classified2011,
  description: 'LandUse2011',
  scale: 30.0
});


texts: In ArcGIS 10.6, how do I build a model to intersect iteratively?

texts:
This is my first time building a model from scratch, so let me know if I've omitted anything in describing this question.

I have a geodatabase that contains shapefiles (polygons) of watersheds. I need to intersect each polygon with a shapefile (polygons) of land use. I can get the following to work, but only for one file at a time.

But I'm not sure I've created a model that will properly loop through all the watersheds in the geodatabase (see below). I used a for loop, but I'm not sure I used it properly. I also need them to be exported as their original name, but in a different geodatabase. Is this how an iterative intersect model would work?


texts: How to do power law test for DMSP/OLS nighttime Satellite data

texts:
I am working with DMSP/OLS satellite's (6 bit) nighttime light data to study the urban expansion of different cities in South Asia. I could not understand the Power law distribution method. I searched over the internet but didn't get how to do that, hence can anyone here help me with that ? The pixel data and the number of a particular pixel occurrences for one city is as follows:

Pixel_VALUE COUNT
3   3
4   14
5   57
6   163
7   182
8   126
9   79
10  80
11  71
12  56
13  67
14  42
15  30
16  30
17  37
18  26
19  23
20  16
21  14
22  16
23  18
24  12
25  20
26  12
27  7
28  10
29  11
30  8
31  21
32  15
33  14
34  12
35  17
36  11
37  9
38  19
39  8
40  12
41  7
42  11
43  16
44  10
45  12
46  8
47  10
48  10
49  7
50  14
51  6
52  12
53  6
54  8
55  13
56  6
57  13
58  10
59  17
60  15
61  23
62  41
63  177


texts: Saving Training Input ROI using the SCP dock within QGIS

texts:
I have just downloaded QGIS and installed the Semi-Automatic Classification Plugin. I am trying to evaluate whether this tool will be of use in creating polygons from high resolution satellite imagery that represent land cover.

Following through the tutorial I am encountering issues when trying to save an ROI.I have successfully created an ROI polygon, and in the classification dock added the MC info and C info when I click on the icon to save the ROI into the training input, the following error message keeps appearing.

2017-05-22T13:29:38 1   Traceback (most recent call last):
              File "C:/Users/sruddick/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\classificationdock\classificationdock.py", line 2240, in saveROItoShapefile cfg.utls.calculateSignature(cfg.shpLay, cfg.rstrNm, [self.ROILastID], cfg.ROIMacroID, cfg.ROIMacroClassInfo, cfg.ROIID, cfg.ROIInfo, 50, 40, "No", "No", UID)
              File "C:/Users/sruddick/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\core\utils.py", line 1661, in calculateSignature bX = cfg.utls.clipRasterByShapefile(tLP, tS, None)
              File "C:/Users/sruddick/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\core\utils.py", line 2779, in clipRasterByShapefile check = cfg.utls.vectorToRaster(cfg.emptyFN, shapefile, cfg.emptyFN, tRxs, raster, conversionType, "GTiff", burnValues)
              File "C:/Users/sruddick/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\core\utils.py", line 5320, in vectorToRaster lPRS = cfg.utls.getEPSGVector(layerPath)
              File "C:/Users/sruddick/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\core\utils.py", line 5222, in getEPSGVector lP.AutoIdentifyEPSG()

AttributeError: 'NoneType' object has no attribute 'AutoIdentifyEPSG'

Has anyone encountered this issue before and can help me resolve this so I can test out the capability of this plugin?


texts: Create GeoJSON Feature from FeatureCollection

texts:
I have a single polygon feature (no holes, no overlaps) which is currently defined as a FeatureCollection in GeoJSON, with type Polygon (not MultiPolygon). I am trying to use it with an online mapping tool, and to do so I need to use a regular Feature in GeoJSON:

I have tried a variety of methods to get it to register and work as a Feature, including: (a) Using the multipart to singlepart tool, which helped it become polygon rather than Multipolygon); (b) using both QGIS and ArcGIS to make the conversion from shapefile to GeoJSON, and (c) editing the GeoJSON file by hand to make it reflect this format:

{
  "type": "Feature",
  "geometry": {
    "type": "Polygon",
    "coordinates": [ [0,0],[0,1],[1,1],[1,0] ]
  },
  "properties": {
    "name": "shape name"
  }
}

For (a) and (b) I still end up with FeatureCollection types, and with (c) nothing maps, despite being recognized as valid by QGIS.

I'd prefer to make this work in QGIS, but am open to suggestions.


texts: I have a map document in ArcMap and I want to find X and Y of my nodes in a special feature layer. I look through my access database and I just find the Shape Index table which has MinGX, MinGY, MaxGX, MaxGY. I will be thankful if you help me how to find the X , Y of my nodes. Generally, from what source is ArcMap locating the nodes on the map?


texts: Designing Spatial Database for temporal Data?

texts:
I am working on a Weather based GIS App.

I have the data from several weather stations, and this data will updated daily (by a web service).

The obstacles I am facing are:

There are currently 40 recording stations, but that might change
Different Stations record different number of parameter, some record 5, some record 7. ect
Some of the parameters are recorded daily (ex: Max Temperature), some are recorded hourly (current temperature) while others are recorded Weekly.
Some facilities at a particular recording station might be decommissioned (for eg: A station that is currently reporting 7 parameters, might report only 5 next year)
Sometimes a parameter might not be reported due to technical issues; Hence I should be able to differentiate between, value=0, null value & value not recorded.

I am working with ArcGIS (desktop & Server) for the spatial part of this web app & using ArcSDE is not a problem , but I think that this kind of database is best left as a pure Database & not bring ArcSDE into it.

Can anybody suggest any books or links that would help me?


texts:
I was googling coordinate related information all the day. Now I am confused with some concepts and terms. I hope some one can help me to make it clear.

Geographic coordinate system

Geographic coordinate is used to describe the location of the earth in a three-dimension which will base on the ellipsoid and a datum. And I always thought that the geographic coordinate are always the latitude and longitude. Am I wrong?
Projected coordinate system

It seems that the projected coordinate system will contains three points:

The projection itself. It means how the earth will be translate to a two-dimension surface. For example the Cylindrical and Pseudocylindrical.
The geographic coordinate system which this pcs working for.
The Mathematics manner to translate the geographic coordinate to projection coordinate.
BTW, is projection coordinate always a point inside the the Cartesian coordinate system? And is projection coordinate only fit for one kind of geographic coordinate system?

Coordinate reference system

Now I have no idea about what a coordinate reference system is. Is the Coordinate reference system has anything to do with the projection?

BTW,is the geodetic latitude is measured on the surface of the referenced ellipsoid?


texts:
How to prevent the map from moving when a user zooms in with scrollwheel?

I would like the center to be preserved in any case.

var mymap = L.map('spotmap', {
   dragging: false
}).setView(spot_LatLng, 18);

var marker = L.marker(spot_LatLng).addTo(mymap);


texts: Group near-points into a single weighted point (dynamically)

texts:
I have a point feature class that is constantly being updated. The points are clustered, so I would like to display near-points as a single, weighted point (and symbolize accordingly). Related information here: Visualize overlapping points

I can do this with the Collect Events geo-processing tool. However, this geo-processing tool produces a static output feature class; it is not dynamic.

Is there a way to dynamically symbolize near-points as a single point in ArcMap?

This is commonly done in web maps, so there must be a way for a platform as powerful as ArcGIS Desktop to do this, right?


texts:
I am in the process of mapping the journey of our data through the different apps used by our service. We use an ArcSDE geodatabase as main data source.
This data is

viewed in ArcMap
edited with forms in ArcMap
published in 3 different webmapping services, one of which has custom webforms.
summary PDFs are also dynamically edited
made available for raw data consultation (data warehouse)
made available for download in 2 or 3 different versions and formats

I'm sure this is quite a frequent case for many of you.
The thing is, it is easy to get lost in the DIFFERENT NAMES used to describe each piece of data (i.e. ATTRIBUTES)
We need to be able to research which attributes are published and where and what are their names in that interface. For instance, what is the attribute called in the web form ? And what is it called in the PDF ?
And we also need to research, for a given final product attribute, of which source attribute(s) it is composed.
Do you know of a tool or of a method to map and research the flow of attributes through different interfaces ?
Ideally, the tool should work both ways.
I've googled a little about this and have come upon articles on Data Lineage tools. But I'm not sure this is the way to go.
I've got a draft of what it could be in Excel but this isn't bi-directional at all.


texts: Displaying Layer Group in QGIS Server with WFS?

texts:
I want to display layers with QGIS server, but the group layer doesn't get displayed, instead I get a bad request response from the server:

The request:

  http://myserver/cgi-bin/91589/qgis_mapserv.fcgi?SERVICE=WFS&VERSION=2.0&REQUEST=GetFeature&typename=94043_CADASTRE&SRSNAME=EPSG:2154&OUTPUTFORMAT=GeoJSON

<ServiceExceptionReport xmlns="http://www.opengis.net/ogc" version="1.2.0">
  <ServiceException code="RequestNotWellFormed">
    TypeName '94043_CADASTRE' unknown
  </ServiceException>
</ServiceExceptionReport>

Server Response

And here is a peek at the layer section in QGIS:

Layers


texts:
I am trying to 'refresh' a cluster layer which is created by L.esri.Cluster.featureLayer and a feature url. I have tried to remove and re-add the layer (with no succees). But the refreshClusters is not a function of a featureLayer, nor of a cluster.featureLayer.

Someone could explain to me?


texts: How to determine temporal land cover change using two rasters? (From-to)

texts:
ArcMap. I am having trouble finding the correct method.
I want to end up with what is essentially a confusion matrix showing what each class from the before image converted to in the after image. I am wanting to create a transition graph like below:

I currently have two rasters with  five qualitative values. Do I need to reclassify the rasters to give them an arbitrary integer value? Would I then need to use the raster calculator? How would I interpret the result?


texts: What is the best open-source web mapping solution for handling 3000+ polygons

texts:
I have 3000+ polygons which represent watersheds (18 main watersheds). At the first stage all the 18 watershed will show and their list will be placed in the left sidebar. If any user mouse over the map only popup will show the watershed name. If any user click on any major watershed on map it will open pop up window with detail. At a larger zoom, all watersheds will show.

I'm asking for your advice which open-source solution is the best for handling that large amount of polygons?

I was thinking about Leaflet or OpenLayers but I don't know can they handle large amounts of polygons.

Also, do you suggest do store them in mysql or?


texts: How to ‘invert’ or ‘reverse’ a vector cutline when gdalwarp'ing a raster?

texts:
I have a raster of North America and I'd like to remove the larger lakes from it, replacing them with the same no-data value as the original raster.

So I downloaded Natural Earth's 10 meter lakes vector data, which has lakes conveniently ranked by scale rank: 0 for large lakes like the Great Lakes, and larger scale ranks for more minor lakes.

The closest I've come with gdalwarp is:

$ gdalwarp -cutline Data/ne_10m_lakes/ne_10m_lakes.shp \
  -csql "SELECT * FROM ne_10m_lakes WHERE scalerank=0" \
  INPUT.tif lakes-only.tif

This actually does the opposite of what I want, it keeps only the Great Lakes portions of my raster, replacing all other non-lake pixels with the no-data value.

My question: how can I invert the cutline, so instead of getting the area inside the lake boundary, I get whatever is outside it? Can I invert the original Natural Earth .shp file, perhaps?

References

In How to invert polygons (or other options for exterior styling)?, GRASS' r.mask was suggested, which has an inverse mode. Is there something similar with OGR/GDAL?

In How to get an ocean shapefile from a very detailed land area shapefile?, the .shp file was converted to a raster, then inverted, and used as a mask. I'd prefer not to generate a raster of non-lakes because my original TIFF raster is very large (high resolution) and I'd prefer not to generate a raster of the same size. As a stopgap though, I was able to do this:

# Create a not-lakes raster of the same spacing and extent as INPUT
$ gdal_rasterize -i -burn 255 -a_nodata 0 -ot Byte \
  -te `gdalinfo INPUT.tif | egrep "Lower Left|Upper Right"|cut -d"(" -f2 | cut -d")" -f1 | sed 's/,//'| tr "\n" " "` \
  -sql "SELECT * FROM ne_10m_lakes WHERE scalerank=0" \
  -tr `gdalinfo INPUT.tif |grep "Pixel Size"|cut -d"(" -f2|sed 's/[,)]/ /g'` \
  Data/ne_10m_lakes/ne_10m_lakes.shp lakes.tif
# Apply the lakes raster to input as some kind of mask
$ gdal_calc.py -A INPUT.tif -B lakes.tif --outfile OUTPUT.tif --calc='A' \
  --NoDataValue=0

Can something similar be done without creating a gigantic non-lake raster?


texts: QGIS Distance Matrix - with line geometry output?

texts: I have a series of Output Area centroids and I want to create a matrix of lines between each point (every point connected to every point), which is effectively what the distance matrix tool does (I think). But I also need the line geometry for visualisation purposes. Can anyone tell me how to do this?


texts: Create a PDF with map imagery and feature data from ArcGIS API for Javascript?

texts:
I would like to allow the user to zoom to a particular zoom, within their preferred extent, and then send this information to a PDF as well as relevant information from visible features. The end result would be a PDF with map imagery, showing the location at the users zoom level, as well as grabbing data from the visible features and filling in a report-like set of data within the PDF as well.

Is something like this possible or maybe even being done already by anyone?


texts:
Let's say that we have a shapefile with a certain projection.

s<-readOGR(dsn=".",layer="Spain")

We also have the airports of Spain as points in another projection.

a<-readOGR(dsn=".",layer="airports")

If we aim to place the points onto the Spain's shapefile we have to arrange the coordinates so they are the same.
Usually, it is done like this:

 a<-spTransform(a,CRS(proj4string(s))

But is it the same with that?

proj4string(a)<-proj4string(s)

If yes, then why isn't the standard way of doing it since it is simpler and have to resort to the spTransform?


texts: How to remove black edges on georeferenced images with GDAL

texts:
After an image has been georeferenced the transformation creates black edges, how could I set these to white instead of black using GDAL?

Would I have to somehow create a colour table and set nodata to white?


texts:
I downloaded data from http://climate.geog.udel.edu/~climate/html_pages/download.html in order to retrieve temperature and precipitation data.

What I actually want is mean temperature and precipitation over a long time-span (e.g. from 1950 to 1990) in gridded format in order to match the grids with other datasets. The problem is that I just do not manage to open the .clim-datafiles provided on the website. The specific datasets I downloaded were air_temp.clim.Z (6.4MB) and precip.clim.Z (7.8MB).

I am using QGIS 3.4. Crayfish (suggested by some other website) does not allow me to open any datasets.
Any solutions?


texts: Topology checker - whole polygon highlighted rather than polygon segment

texts:
When I check for overlaps, topology checker is only identifying errors at the polygon level but not the polygon segment where the errors exist (the whole polygon is highlighted in blue).

I believe this may be happening because the nodes aren't snapped together but the overlap is so small (I checked node coordinates between adjacent polygons and there is no difference in at the 4 decimal places shown). I can fix the error manually but it involves pulling each node apart for the entire polygon as I don't know where the overlap actually exists and then snapping all nodes back but this takes a long time to complete.

Is there any way I can automate or semi-automate this snapping process?

Attached image shows the highlight polygon in blue with the error but not the polygon segment that normally is shown in red:


texts: Calculating accurate distances in OpenLayers 3?

texts:
I'm calculating distance in OpenLayers 3 (3.9.0).

I've got a function included in my project for measuring line/area features, basically lifted from the OpenLayers Examples.  I'm using haversineDistance to calculate accurate lengths/areas, but keep getting an error preventing the numbers from being reported.

Basically, when I start to measure in the browser and console.log the measurements so I can see them, they appear to be reporting accurately; however, a js error in ol.js is preventing the measurements from being returned in the webpage so the viewer can see them.  Error is 'Uncaught TypeError: Cannot read property 'length' of undefined'. Reports as occuring in ol.js on line 69.

My code for the function is:

var formatLengthFt = function(line) {
var length_ft = 0;
var coords = line.getCoordinates();

for (var i = 0; i < coords.length; i++) {
    var c1 = ol.proj.transform(coords[i], 'EPSG:3857', 'EPSG:4326');
    var c2 = ol.proj.transform(coords[i + 1], 'EPSG:3857', 'EPSG:4326');

    length_ft += wgs84Sphere.haversineDistance(c1, c2);
    console.log(length_ft);
}
return length_ft;
}

Before employing haversineDistance, I simply used line.getLength(), but then realized the distances/areas were all inaccurate.  So, tried the solution below with success, just can't report the numbers in the webpage.  I've also tried using turf.js, with some level of success.  At any rate, I just need a method for getting accurate distances/areas and would prefer to get the code above working correctly, but open to other solutions as well.


texts: Adding Esri imagery to QGIS using "services.argisonline.com"

texts:
I'm trying to add this particular imagery in QGIS:

http://services.arcgisonline.com/arcgis/rest/services/World_Imagery/MapServer/3

But when i write the phyton srcipt i get the message that the URL is incorrect.

I'm getting this typical message:
"Cannot open GDAL dataset http://server.arcgisonline.com/ArcGIS/rest/serââvices/World_Imagery/MapServer?f=json&pretty=true:
Failed to create temporary file:/tmp/file.dat"

Anyone know what is the correct URL for High Resolution 30cm Imagery?


texts: How to resolve an Invalid Topology [Incomplete void poly] error from a geoprocessing service?

texts:
I am consuming a geoprocessing service from my Flex application.

My geoprocessing service takes the input geometry, and intersects another layer. It then runs a Summery stats on the intersection output and returns the output.

When I do a query or identify, and send the resulting geometry to the GP service, it runs fine. But when I Keep it in memory and do other things, which take a lot of time, and then send the geometry to the GP service, it errors out, with this error.

The only difference I find in the later geometry, is that the rings are not in the same order as in the geometry I get from the server.

What can I do to preserve the order of the rings?


texts: I need to download a DEM in WGS format for Nepal ( around 86 degrees longitude and 27 degrees latitude), that can be exported to a GIS as raster for analysis. Does anyone know of any good, free data source?


texts:
im currently blueprinting something i wanna develop. Thinking about the stuff that would be needed, i figured out that i would need to have a pgRouting/postGis Database somewhere. The program i wanna make, should calculate data and return the shortest path.

Is there any online service who offer this kind of service (an accessible pgRouting/postGIS Database)?? even for a fee ??
Since i couldn't find one im asking here.

Thanks a lot.

Update from comments:

I want to develop an app that USES a routing service, to display the (elsewhere) routed data on the phone. Like oh i am here, thisService->please tell me how to get from x to y.