texts: Missing distance field after Spatial Join?

texts:
When I use Spatial Join, the distance attribute is not being added to the resulting shapefile.

I'm trying to use Spatial Join to assign each of 16,151 points to one of three zones (see screenshot), and also assign a distance field so I can include points that fall outside the zones (in this case there is only one, but I will be doing this with many other point shapefiles).

I'm doing a JOIN_ONE_TO_ONE operation, with the CLOSEST match option, a Distance Field Name called DistZone, and a search radius of 2 km.

The tool successfully adds a Zone attribute but no DistZone attribute.

The point shapefile has a lot of attributes (80), but I read that the max is 255, so that should not be the problem. All inputs are in the same coordinate system.

This was working successfully for me yesterday, and now it has stopped working. Any idea why the distance field is no longer being added?

UPDATE: I've answered my own question. Since the shapefile has so many attributes, I didn't realize that the DistZone field was being added all the way on the left, while the Zone field is added all the way on the right.


texts:
I have a large raster stack of about 500 rasterlayers.
What I want to do is to extract specific layers and store them in a different strack.
The layers for which I am interested in, have a step of 8 (e.g. I want layers 1,9,17,25,...).
I have tried to do this with the following loop but I get Error in compareRaster(x, value) : different extent.

What is wrong with my code?

new<-stack()
for (i in seq(1,500,8)){
  new[[i]]<-old[[i]]
}

I managed to solve this by setting new as list() instead of stack().

Can someone explain to me the error?


texts: QGIS does not update data in Geoserver WFS-T


texts:
I’m trying to plot a shapefile of points that correspond for all impact craters on Pluto.
I can't get the points to math with the impact craters from the global mosaic of Pluto's surface. In the best o
QGIS has a CRS named GCS_Pluto_2000, but when I set the global mosaic, the shape of points and the project in QGIS to the same CRS, the points are plotted completely wrong.

I might be missing something but I can’t guess what.

The link for the global mosaic and the CSV with the points:

https://astrogeology.usgs.gov/search/map/Pluto/NewHorizons/Pluto_NewHorizons_Global_Mosaic_300m_Jul2017

https://astrogeology.usgs.gov/search/map/Pluto/Research/Craters/Craters_PlutoCharon_System_Robbins


texts:
Zooming around Northern California on Google Earth reveals that Google have extremely high resolution imagery for not just cities but also rural and wilderness areas.  I would guess it's something around ~10cm, which must be from aerial imagery not satellites, right?  If you'd like see what I mean search for "Big Bar, CA" in Google Maps or Google Earth and fly around a bit.  What's more, it seems to be collected annually (at least for the last two years).

I have no idea how much of the country has this high resolution coverage, I've only just noticed it in my neck of the woods.  And I'm specifically interested in the Northern California imagery.

Does anyone know the source of this data?  Or is there a way to find the source?  Google Earth lacks the familiar copyright notices at the highest resolutions and I'm stumped about how to proceed beyond that. I might be interested in purchasing the data if I could find the source.


texts: calculating differences in species distribution over time

texts:
I have been searching a lot but I didn't find any useful solution yet.
I have a shapefile representing range polygons for different species of primates and rasters representing current and future environmental conditions or predictors. What I would like to do is to assess for each species if there is any difference when subtracting  future - present rasters of each predictor within species range. I have been trying to do this in ArcGIS 10.5 by extracting all predictors using primate shapefile as a mask, using raster calculation and reclassify tools, but it came up with results for primates and not for each primate species.

Any suggestion about how to do this in ArcGIS?

I am a bit familiar with R (more than ArcGIS), so if any script/package is available to do this in R I don't mind to give it a try too.


texts: ArcGIS Pro Segmentation Generates Raster Output with NoData Holes

texts:
I am trying to segment an image using ArcGIS Pro and I am sometimes getting NoData holes in my output raster. This has happened via Segment Mean Shift geoprocessing tool (and coded), and via Classification Tools > Segment. The raster I'm using to segment is an output from Extract by mask.
I've tried a work around of segmenting the original raster instead of a subset (i.e. before Extract by Mask) and then use Extract by Mask (I don't want to classify the whole image is it's much larger than my AOI). However, it doesn't recognize the output from Extract by Mask as a segmented image.
Is there a way to avoid NoData holes (i.e. ensure that the output segmented raster covers the original raster completely?

Subset Image for Segmentation

Segmented Raster with NoData Holes (white color inside pink borders)
Please note question also posted on ESRI Community Site: https://community.esri.com/t5/imagery-and-remote-sensing/inconsistent-segment-mean-shift-output-with-nodata-holes/td-p/1027816


texts:
I was debating with one of my collegues on OpenLayers v/s Leaflet. I made a point that OpenLayers is much better API if we wish to build a project, where you need direct connectivity to the Geoserver and PostGIS.

Then I found Open Data Kit, which looks pretty new but has the features of connectivity with the Geoserver and PostGIS.

So my project details are as follows,

Use the map interface to fetch Feature Info
Create a customized tool that takes the lat/lon from user as to where he/she clicks on the map and then fetches the Climate Data from the raster (which is handled by a py script on the server)
Allows user to upload excel, which is sent to the py script, which returns a GeoJSON, which creates Vector Features on the map
Allow user to create vector polygons, which will fetch the Features it intersects from the WFS Layer.
Fetches Layer from the PostGIS Datastore on GeoServer and displays the layers on the map.

So now I am confused on which is better and why using OpenLayers over Leaflet makes more sense or not?


texts: Selecting markers within a GeoJSON polygon (Leaflet)

texts:
I have a GeoJSON FeatureCollection object made up of several polygons and I would like to select the markers which lie within each polygon. The information contained within the markers is to be used to show average stats for each polygon.

I have seen some suggesting using http://turfjs.org/docs/#within or http://turfjs.org/docs/#inside.

Markers are simply in L.Marker form. I intend to use the onEachFeature function to have a popup for each polygon. Has anyone done this before and can recommend the best method?


texts:
I want to create sankey diagrams (see: https://en.wikipedia.org/wiki/Sankey_diagram) from date like:
zipcodepoly1 zipcodepoly2 value
12345        45789        250000,00
....

I have polygons with attributes for the zipcodes and I have the values.
It seems that it should be possible with geoplot.sankey? But I have to define a line layer.
It should also be possible with geometry generator in qgis. Any Hints?
How to create Sankey diagrams on a map in QGIS or with Python?
Because of the inspiration from https://anitagraser.com/2019/05/04/flow-maps-in-qgis-no-plugins-needed/ I create a flowmap that is nearley a sankey. I am not realy satisfied by the scale but as first effort it is OK.


texts: Changing color between polygon gaps (QGIS 2.18.3)

texts:
Trying to figure out how to solve this. I'd like to change all of the white area some other color. Looking at the space spaces closely shows an optical illusion going on.

Anyone know how to change the white areas any other color?


texts: Python code execution error in QGIS WIEN

texts:
I have tried to make several different actions in QGIS but keep getting similar error messages when I try to select the output filename. The error message reads as follows:

An error has occurred while executing Python code:

TypeError: QgsEncodingFileDialog(QWidget parent=None, QString caption=QString(), QString directory=QString(), QString filter=QString(), QString encoding=QString()): argument 3 has unexpected type 'QPyNullVariant'
Traceback (most recent call last):
  File "/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools/doVectorGrid.py", line 333, in outFile
    (self.shapefileName, self.encoding) = ftools_utils.saveDialog(self)
  File "/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools/ftools_utils.py", line 327, in saveDialog
    fileDialog = QgsEncodingFileDialog(parent, QCoreApplication.translate("fTools", "Save output shapefile"), dirName, filtering, encode)
TypeError: QgsEncodingFileDialog(QWidget parent=None, QString caption=QString(), QString directory=QString(), QString filter=QString(), QString encoding=QString()): argument 3 has unexpected type 'QPyNullVariant'

Python version: 2.7.10 (default, Oct 23 2015, 19:19:21) [GCC 4.2.1 Compatible Apple LLVM 7.0.0 (clang-700.0.59.5)]
QGIS version: 2.14.3-Essen Essen, exported

Python Path:
/Users/nickjames/.qgis2/python/plugins/postgis_geoprocessing
/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins/processing
/Users/nickjames/.qgis2/python/plugins/pktools
/Applications/QGIS.app/Contents/MacOS/../Resources/python
/Users/nickjames/.qgis2/python
/Users/nickjames/.qgis2/python/plugins
/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins
/Library/Frameworks/SQLite3.framework/Versions/C/Python/2.7
/Library/Frameworks/GEOS.framework/Versions/3/Python/2.7/site-packages
/Library/Python/2.7/site-packages/numpy-override
/Library/Python/2.7/site-packages/matplotlib-override
/Library/Frameworks/GDAL.framework/Versions/1.11/Python/2.7/site-packages
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python27.zip
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-darwin
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac/lib-scriptpackages
/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-tk
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-old
/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-dynload
/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/PyObjC
/Library/Python/2.7/site-packages
/Users/nickjames/.qgis2//python
/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools
/Users/nickjames/Desktop/qgis

This particular error was in response to trying to create a vector grid and when clicking the browse button for the 'output shapefile' button this error came up. This may seem like a repeat question but on other threads no working solution was given for my situation.

I am running QGIS ESSEN (2.14) on a macbook pro which I just re-installed in an attempt to fix the problem but to no avail


texts: The Logic Marker Moving with the Offset Function in QGIS 3.6.2

texts:
In QGIS the location of the marker can be changed, without changing the coordinates in the Attribute Table. This is done by changing the X or Y values in the Offset function, located in the Layer Styling panel.

When a positive value is adding to the X bar the Marker moves to the left, and vice versa. Just like a X-axis standard system.

However, adding a negative value to the Y bar will move the Marker up. Which is opposite to the a Y-axis.

What is the logic that the offset function is working by?


texts: Batch processing a folder of pairs of grids in QGIS-python

texts:
I am creating a script to run some processing tools on pairs of rasters - in this case VelocityX and VelocityY raster grids.  The files are called .Vx and .Vy and are uniquely identified by names in this format:

StandardBaseName_0001.Vx
StandardBaseName_0001.Vy
StandardBaseName_0002.Vx
StandardBaseName_0002.Vy
etc etc

This question has some code to loop over a folder of rasters - can this be adapted for identifying each 'pair' of rasters and feeding their names into variables that I can give to the processing algorithms?

def findRasters (path, filter):
    for root, dirs, files in os.walk(path):
        for file in fnmatch.filter(files, filter):
            yield file

for raster in findRasters(INPUT_FOLDER, '*.tif'):

I was thinking that matching the name (without extension) of the .Vx raster to the name of the .Vy raster might work but I'm not sure how to do it.


texts: Pivot table output name and ModelBuilder failing?

texts:
I have created a model that utilises the Pivot Table tool. When I run the model, the tool creates a DBF file in a file GDB. However the naming convention of the file is incorrect. It just has the name "dbf", rather than the output table name stated in the parameters. It subsequently stops the model from continuing giving me an error:

  ERROR 000229: Cannot open
  C:\AIS\Model_testing_area\AIS_2013_Density.gdb\PivotTable.dbf Failed
  to execute (Add Join). I have tried to trim the name, change the
  workspaces (in_memory, %currentworkspace%) however to no avail.


texts: Attribute Query for First/Last record and every nth inbetween

texts: I have a feature class with over 5000 points.  There are too many to label so I would like to be able to grab the first, the last and every 30th or 50th in between.


texts:
I have many polygons that present wind directional movement.

I use Arc-Map 10.3.1.

How can I add arrows next to each polygon, to show the wind direction in any polygon?


texts: How do I determine which polygon feature with an SDO_geometry doesn't partially overlap a polyline feature with an SDO_geometry?

texts:
I have 2 datasets in the Oracle spatial format, built from Smallworld data;

One with polygon features
One with polyline features.

I'm looking to determine which polygon features do not contain (part of) a polyline feature, but the Smallworld application the organization uses has just been redeveloped and doesn't really allow spatial queries yet.
Both features have an SDO_geometry attribute, but I'm a bit stuck on how exactly to approach the comparison. Both are pretty massive datasets (>1.000.000 records each), so I'm thinking FME should be able to be faster than an Oracle query, but I might just be overthinking it. How should I go about this comparison?


texts:
Is it possible to disable log messages in QGIS Dufour?

When I work in QGIS 2.0, log messages increase to more than 2000 and it makes QGIS very slow. Then I save project and open it again, then QGIS work fine and everything is fast. But during work in QGIS, logged messages iNcrease and it works slowly.

Here I have enclosed the screen shot of the log messages from QGIS


texts:
I have the following problems with Surfer:

1, I have a grid file and a bln but i'm not sure that they overlap each other exactly. I set the coordinate system with grd and bln! (WGS84) but there is no changes happend...

2, I tryed to blank the grid file around bln but i got this message:

Does anyone can help me how to fix the exact position to overlap each other and blank the grid around bln?


texts: I am in the process of developing predicting map in MAXENT. I have current data of distribution of my trees species in tropical forest. I need to use Maxent software and Arcgis 10.3 to find the potential area where this tree species can be found. My preoccupation is to know which environmental data and where can I find it to realize my predicting map, and how can I use Maxent software to do it.


texts:
I previously had 10.5 and opened up an ArcGIS 10.2 versions MXD.  This week, I upgraded to 10.6 and tried to open the MXD from that same 10.2.  That MXD would not open. All it would do is "not responding".  I copied, renamed it and even re-downloaded it from the email it was sent from.  Nothing.  Corrupted?

I opened up older mxd's of mine when I had 10.1.  Another Mapper used different versions 10.5 and 10.2, could not open that one either.  But, could open others.

Are there any solutions?


texts:
I would like to change some default QGIS options.
For example, the default font used by QGIS composer on my system is ubuntu. I would like to change that.

I would also like to get no border when creating an object.
Is there a config file I can create or modify in my .qgis user folder?

A user option file (like in scite editor) for example, would be valuable if it does not exists yet.


texts: QGIS 2.18.17: python warning missing init

texts:
Every time I open QGIS 2.18.17, I get the following python warning in the protocol window:

2018-03-02T12:30:52 1   warning:C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py:607: ImportWarning: Not importing directory 'C:\PROGRA~1\QGIS2~1.18\apps\Python27\lib\site-packages\win32\test': missing __init__.py
          mod = _builtin_import(name, globals, locals, fromlist, level)

        traceback:  File "<string>", line 1, in <module>
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 298, in loadPlugin
            __import__(packageName)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 607, in _import
            mod = _builtin_import(name, globals, locals, fromlist, level)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python/plugins\processing\__init__.py", line 29, in <module>
            from processing.tools.general import *              # NOQA
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 607, in _import
            mod = _builtin_import(name, globals, locals, fromlist, level)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python/plugins\processing\tools\general.py", line 31, in <module>
            from processing.core.Processing import Processing
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 607, in _import
            mod = _builtin_import(name, globals, locals, fromlist, level)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python/plugins\processing\core\Processing.py", line 62, in <module>
            from processing.algs.saga.SagaAlgorithmProvider import SagaAlgorithmProvider
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 607, in _import
            mod = _builtin_import(name, globals, locals, fromlist, level)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python/plugins\processing\algs\saga\SagaAlgorithmProvider.py", line 33, in <module>
            from SagaAlgorithm import SagaAlgorithm
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python\qgis\utils.py", line 607, in _import
            mod = _builtin_import(name, globals, locals, fromlist, level)
          File "C:/PROGRA~1/QGIS2~1.18/apps/qgis-ltr/./python/plugins\processing\algs\saga\SagaAlgorithm.py", line 20, in <module>
            standard_library.install_aliases()
          File "C:\PROGRA~1\QGIS2~1.18\apps\Python27\lib\site-packages\future\standard_library\__init__.py", line 483, in install_aliases
            import test

How can I fix this warning? I think it may be the reason why I am not able to execute the v.net.connect algorithm.


texts: Republic of Ireland .csv to mark out points on site. Which map/co-ordinate system to use?

texts:
I need to be able to create points in Arcmap 10.1 (or 10.4) in a geo-referenced map layer which I can then export as a .csv file and then use in my gps (GeoExplorer 6000) to mark out those points on site.

I have a process that works for this in England/Wales/Scotland but I am not sure how to do this in the Republic of Ireland or Northern Ireland as I don't know which co-ordinate system to use or even which maps to download/buy.

Any advice?


texts: How to render iD editor-style streetmaps?

texts:
Is there a simple way to render OSM data (let's start with just streets and names for a single metro area) client side as nodes and edges like the new iD editor?

I'm comfortable with d3 but am having trouble seeing exactly how they're pulling it off.


texts:
How do I export to PDF or SVG file? I tried SimpleSVG and it consistently gives me an error message.

Then I tried to print to PDF, but it comes up blank and I can't figure out what I'm doing wrong.

I'm on a Mac running 10.11.16 "El Capitan" and using QGIS 2.18.2 "Las Palmas"

End goal is to get something I can edit in Adobe Illustrator (I have Illustrator CS5 and IllustratorCC and IllustratorCS2)


texts:
SQL statements in Python continue to plague me. I'm trying to get the following expression written in a stand alone script:

This is what I've got so far but I don't think it's even close to the right syntax. I continue to get a syntax error at each iteration of this statement I try. Is there a guide to placing operands, field names, values, double and single quotes?

QF5 = "Week"
QF6 = "Wk"
WC4 = '" QF5" + ' <> "'" + QF6 + "'" AND '" QF6" + ' <> "''''"


texts: Geographic Transformations using ArcPy vs. Coordinate Systems used by ArcGis online basemaps

texts:
I have three shapefiles with undefined projection, but I know which one was used (Lisboa_Hayford_Gauss_IPCC). I need to first Define Projection. It works. I then need to transform it to WGS_1984, so I've used:

import sys, arcpy
#globals set env.workspace and other used paths
import globals as g

def main():
    # set local variables
    inDataset = ["rede_Coimbra_node.shp", "rede_Coimbra_link.shp", "rede_Coimbra_zone.shp"]

    prjFile = g.coveragePath + "Lisboa Hayford Gauss IPCC.prj"
    prjFileDatum73 = g.coveragePath + "Datum 73 Hayford Gauss IPCC.prj"
    prjFileWGS = g.coveragePath + "WGS 1984.prj"
    # Geographic transformation -
    transformation1 = "Datum_Lisboa_Hayford_To_Datum_73_1"
    transformation2 = "Datum_73_To_WGS_1984_4"

    sr1 = arcpy.SpatialReference(prjFile)
    print sr1.name
    #"spatialRef.name = Lisboa_Hayford_Gauss_IPCC"
    #sr2= arcpy.SpatialReference(prjFileDatum73)
    #print sr2.name
    #print (arcpy.ListTransformations(sr1, sr2))
    #spatialRef.name = "Datum_73_Hayford_Gauss_IPCC"
    sr3 = arcpy.SpatialReference(prjFileWGS)
    print (arcpy.ListTransformations(sr1, sr3))
    #spatialRef.name = "GCS_WGS_1984"

    try:
        print "..."
        for fc in inDataset:
            fcFullPath = g.coveragePath+fc
            #arcpy.DefineProjection_management(fcFullPath, prjFile)
            outfc= fc.strip(".shp")
            outfc = arcpy.env.workspace+"/"+outfc+"_wgs1984"
            print outfc

            arcpy.Project_management (fcFullPath, outfc,prjFileWGS, u'Datum_Lisboa_Hayford_To_Datum_73_1 + Datum_73_To_WGS_1984_4',prjFile)
            #arcpy.Project_management (fcFullPath, outfc,prjFileDatum73, transformation1,prjFile)
            #arcpy.Project_management (outfc, outfc+"_wgs84",prjFileWGS, transformation2,prjFileDatum73)

    except Exception, e:
        # If an error occurred, print line number and error message
        tb = sys.exc_info()[2]
        print "-->Line %i" % tb.tb_lineno
        print e.message

if __name__ == '__main__':
    main()

I also have a layer from OSM, which I also projected to WGS_1984 GCS and it fits perfectly to the World Imagery provided by ArcGIS Online or OSM basemap.

However, my projected three shapefiles (now converted to feature classes) don't seem to fit so well. They have a deviation from 1-2m to the raster basemap. Did I do something wrong or is it because the online basemaps are in other CRSs and I can't change that?


texts:
I want to save some geometry to a shapefile with the python27 module fiona and manually specify the coordinate reference system (CRS).

Here's what I tried:

import fiona
import fiona.crs

outpath = "shape.shp"

schema = {'geometry': 'Point', 'properties': {'place': 'str'}}
point = {u"type": u"Point", u"coordinates": [0.0, 51.4769]}
properties = {'place': 'Greenwich'}
wgs84 = fiona.crs.from_epsg(4326)

with fiona.open(outpath, 'w', crs=wgs84, driver='ESRI Shapefile',
                schema=schema) as layer:
    layer.write({'geometry': point, 'properties': properties})

Unfortunately the shape.prj file that's produced is empty, hence the CRS is not specified in the output file.

A little diagnostics: After running the script, the layer has the CRS specified:

In [10]: layer.crs
Out[10]: {'init': 'epsg:4326', 'no_defs': True}

Any idea what's going wrong there?

I can work around this issue by downloading the .PRJ file from spatialreference.org, placing it in the same folder and renaming it shape.prj. I guess I'll figure out how to do that programmatically if there's no direct solution (should be easy).


texts: Why NASA world view (EOSDIS) have black strips?

texts:
There is the famous NASA World View site based on EOSDIS.

I relatively new to this and I wonder why & what are those black strips on the existing coverage? I also notice that the position of the strips is not the same on all days.

(This is different from the black part which is due to delay in processing time)


texts: Create line layer around selected polygons

texts:
I want to select numerous polygons out of a polygon layer and create an outer line surrounding only selected features. All polygons are part of the same layer.

For example, I want to have an outer line as a new layer around the selected features in photo.


texts: How to create an error map to support an average kernel density map?

texts:
I've created an average kernel density map by running KDEs on points stacked within the same spatial extent.  For example, say we have three point shapefiles representing seedlings in three different forest gaps of the same shape and size.  I ran a KDE for each point shapefile.  The output from the KDE were then stacked based on spatial extent in order to calculate the average in Arc's raster calculator, for example: Float(("KDE1"+"KDE2"+"KDE3")/3).  Here is the final product:

Now I'm interested in creating a map depicting the error associated with the averaged KDEs.  I hope to use the error map to visually depict how much error is associated with the hotspots (e.g. is the SW hotspot due entirely to the points in one gap?).  How should I go about creating a map of the error associated with the averaged KDEs?  Would MSE be the most appropriate measure of error in this case?


texts: Adding Oblique imagery basemap to ArcMap?

texts:
Is it possible to include oblique imagery in ArcMap?

The default basemaps only offer orthographic view basemaps.

I need to add an oblique view and show my other layer shapefiles in this oblique view.


texts: How to get FGDB support in GDAL 1.9 without compiling?

texts:
I was hoping that I could use the 'master' (nightly-build) version of QGIS to get fgdb support without compiliing (I'm not a developer and frankly have wasted days of productive time trying to get things to compile with no luck).

Since it seems that gdal 1.9 support fgdb, I thought I could just download the developer version. However, the latest version of QGIS (191a229), even with gdal 1.9 does not seem to allow me to add my gdb file.

Is the gdal version included with the master build just not compiled with fgdb support? Does anyone know if there's a way to do this without compiling a new version?


texts:
For an iPAD application (based on WhirlyGlobe) we've two modes: offline and online.

It's a meteorological application so each day around 1600 layers (worldwide) per day are created on the servers. We've a backend service (based on Spark) which is generating mbtiles.

The iPAD loads these layers using a tile web service (online) or just loading the MBTiles (offline).

For the offline we need around 3GB per day and this is quite difficult to maintain for an iOS App.

We're using MBTiles because we use Mapnik for styling the layers. The simbology of the layers are quite specials so we need a really good/flexible render engine.

We've tried to use Mapnik Vector Tiles (MVT) but surprisingly the size of mbtiles of images are quite similar to mbtiles of MVT.

For offline, the max_zoom of the layers is 4.

Should I consider other format of data? Does anybody have an idea about how the size could be reduce?


texts: I have been trying to impot netcdf files following this guide online http://www.ggiuliani.ch/download/netcdf_qgis_GG.pdf into qgis 2.10 with no success. ALthough the file is imported, the dimensions are wrong and I am not sure how to assign lat/lon information to the selected var layer (eg temperature). I have tried plotting the grid in python using matplotlib and everything is fine with the file.


texts: How to use list broken Data sources for layers saved in Windows not within an MXD

texts:
I have seen a few question related to listing broken data sources  of layers within an MXD  MXD list broken links  See Link  .

I am interested in layers that are not within an MXD. That are sitting within a windows file directory

It seems the syntax is the same , maybe I am passing a file string that is not in the correct format my code is as follows

folderPath = r"C:\PROJECT\02_USER_FOLDERS\JOHN_SMITH\LAYERS"

def Get_Layer_files(Folderpath):

    import arcpy,os

    ws = arcpy.env.workspace = folderPath
    arcpy.env.overwriteOutput = True
    Layer_Container =[]

    Layer_List = arcpy.ListFiles("*.lyr")

    for layer in Layer_List:
        Full_Name = folderPath + "\\" + layer
        Layer_Container.append(Full_Name)
    return Layer_Container

import arcpy

my_list = Get_Layer_files(folderPath)

for file in my_list:

    lyrfile = arcpy.mapping.Layer(file)
    lyrlist = arcpy.mapping.ListLayers(lyrfile)
    for lyr in lyrlist:

        brokenlist = arcpy.mapping.ListBrokenDataSources(file)
        for lyr in brokenlist:

            print lyr

At the moment I think I am passing my layers file that are saved here

folderPath = r"C:\PROJECT\02_USER_FOLDERS\JOHN_SMITH\LAYERS"

That are each listed as

C:\PROJECT\02_USER_FOLDERS\JOHN_SMITH\LAYERS\a_layer.lyr

C:\PROJECT\02_USER_FOLDERS\JOHN_SMITH\LAYERS\b_layer.lyr

Is this what should be passed into the

brokenlist = arcpy.mapping.ListBrokenDataSources('C:\PROJECT\02_USER_FOLDERS\JOHN_SMITH\LAYERS\a_layer.lyr')

I am I missing something simple ?


texts: I'm using QGIS and I would like to measure the distance between the centroid of King County in the US and the centroids of all other US counties. I have imported the shapefile and I have executed the "centroid" tool. However, I am not able to set "King County" as a starting point. When I do "Distance Matrix", there is no way to set King County as the starting point. I would like to have "King County" as the reference point and then measure the distance from it to all other counties.


texts: Arcpy Toolbox sending and error message in specific condition

texts:
I'm trying to understand how the ordering works in the Python Toolbox.

What I need is, everytime the user updates the parameter, the tool will check if a feature exists there, and display an error immediately if it doesn't.

But the setErrorMessage on my parameter inside updateParameters is not actually putting the parameter in error state, so I'm not sure what to do here, the documentation even shows that this is the way to do it.

Tool code :

class CustomTool(object):

expected_class = 'fclass.shp'

def __init__(self):
    """Define the tool (tool name is the name of the class)."""
    self.label = "Custom Tool"
    self.description = '''
    Automatization script
    '''
    self.canRunInBackground = False

def getParameterInfo(self):
    """Define parameter definitions"""
    params = [
        #Only parameter is a multivalued workspace
        arcpy.Parameter(
            'datasets',
            'Geodatabase(s) containing the data',
            'Input',
            'DEWorkspace',
            'Required',
            True,
            None,
            None,
            True
        )
    ]
    return params

def isLicensed(self):
    """Set whether tool is licensed to execute."""
    return True
    def updateParameters(self, parameters):
        """Modify the values and properties of parameters before internal
        validation is performed.  This method is called whenever a parameter
        has been changed."""
        gdbs = parameters[0].valueAsText.split(';')
        for gdb in gdbs:
            arcpy.env.workspace = gdb
            #If the expected class does not exists
            if not arcpy.Exists(diretiz_class):
                #This will not work :
                return parameters[0].setErrorMessage('{0} is not present in the {1} workspace'.format(expected_class,gdb))
        parameters[0].clearMessage()
        return

    def updateMessages(self, parameters):
        """Modify the messages created by internal validation for each tool
        parameter.  This method is called after internal validation."""
        return

    def execute(self, parameters, messages):
        """The source code of the tool."""
        #WIP
        return


texts: Filter negative SO2 image collection in GEE

texts:
I'm collecting Sentinel-5p tropospheric SO2 column data. I want to remove negative SO2 band values from my collection. How do I achieve this?
Here is my code:
// define a 'collection' of images
var collection = ee.ImageCollection('COPERNICUS/S5P/NRTI/L3_SO2')
.filter(ee.Filter.calendarRange(12,12,'month'))

.select('SO2_column_number_density')

.filterBounds(table)
.filterDate('2018-05-01', '2021-05-31')

.filterMetadata("SO2_column_number_density", 'less_than', 0);

// print(collection)
// var col=collection.values.filter(ee.Filter.notNull(['SO2_column_number_density']))

print(collection)

var no2Image = collection.mean(); // mean across each period
print('NO2 image', no2Image); // debug info in console
var band_viz = {
min: 0,
max: 0.0002,
opacity: 0.45,
palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};
// The region of interest - a planar rectangle around London/S.England.
var rect = table

Map.setOptions('hybrid'); // set map to hybrid backdrop
Map.centerObject(rect, 7); // zoom in to the defined area
Map.addLayer(no2Image, band_viz, 'Sentinel 5P NRT N02 - Nov 2019');

// Create a task that you can launch from the Tasks tab.
// commented out here - but saves off geoTIFF files
Export.image.toDrive({
image: no2Image.clip(rect),
description: 'SO2_red_sea_december',
folder: 'GEE',
fileNamePrefix: 'SO2_change',
region: rect,
fileFormat: 'GEOTIFF',
scale: 1000
});

//


texts: Encoding user-defined datum and projection in GeoTIFF for a non-Earth image?

texts:
I'm writing code to export maps (photos, elevation models, other measurement data) of fictional celestial bodies from a videogame into GeoTIFF with correct georeferencing, so that they could be imported and processed in GIS software. To that end, I figured out I need to specify a user-defined datum, in which I provide the axes of an ellipsoid that correspond with a planet/moon size.
I want to export two types of maps that are generated within the game - a plate carré of the entire body, and an orthographic map "zoomed in" on some area of the body. With the former, I can deal in a pinch by setting GTModelTypeGeoKey to ModelTypeGeographic, encoding the datum and encoding tiepoints in degrees (lat/long). But for orthographic image, I assume I should be using ModelTypeProjected.
When I try to use GTModelTypeGeoKey = ModelTypeProjected, and then set the datum (Geog* parameters) and user-defined projection, I find that gdalinfo only shows the Projected CS definition (with no custom datum in sight), and QGIS doesn't know how to line up such an image with a ModelTypeGeographic-encoded plate carré.
My questions thus are: Is it possible to encode both a custom datum and a custom projection in GeoTIFF simultaneously? If so, how? If not, what is the recommended approach of georeferencing projected images of non-Earth bodies, for which no standard datum can be used?


texts: Automating Redshift data dump to a PostgreSQL DB every 12 hours

texts: I have a redshift database that changes multiple times a day. I want to use this "live" data to make web maps and I have a Postgres DB where all my web-mapping data is stored. I want to continuously dump data from the ever-changing Redshift DB to my maps via the Postgres DB. Ideally, I'd want it to update every 12 hours at most since I think redshift refreshes every 6hrs. What is the simplest way to do this? I have experience with SQL so if that's the route it would be helpful. No extensive experience with complex server-side systems.


texts: Displaying TinyOWS layer in OpenLayers 4.2

texts:
I went through all the other posts related to this and it got me pretty far, in terms of no errors. But my xml response doesn't contain any geometry points from my source. I'm using the demo from TinyOWS to keep this generic for everyone and since the packaged demo is for openlayers 2.
What I have as my Layer:

new ol.layer.Vector({
              source: new ol.source.Vector({
                format: new ol.format.WFS({schemaLocation: "/ms4w/apps/tinyows/schema/"}),
                url: function(extent) {
                   return 'http://localhost/cgi-bin/tinyows.exe?service=wfs&request=DescribeFeatureType&version=1.1.0&typename=tows:frida&srsname=EPSG:3857'

                }
              })
      })

And the response:

    <xs:schema xmlns:tows="http://localhost/tinyows/" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns="http://www.w3.org/2001/XMLSchema" xmlns:gml="http://www.opengis.net/gml" targetNamespace="http://localhost/tinyows/" elementFormDefault="qualified" version="1.1">
<xs:import namespace="http://www.opengis.net/gml" schemaLocation="http://schemas.opengis.net/gml/3.1.1/base/gml.xsd"/>
<xs:element name="frida" type="tows:fridaType" substitutionGroup="gml:_Feature"/>
<xs:complexType name="fridaType">
<xs:complexContent>
<xs:extension base="gml:AbstractFeatureType">
<xs:sequence>
<xs:element name="gfshapeid" type="short" nillable="true" minOccurs="0" maxOccurs="1"/>
<xs:element name="gftypid" type="short" nillable="true" minOccurs="0" maxOccurs="1"/>
<xs:element name="gfname" nillable="true" minOccurs="0" maxOccurs="1">
<xs:simpleType>
<xs:restriction base="string">
<xs:maxLength value="50"/>
</xs:restriction>
</xs:simpleType>
</xs:element>
<xs:element name="geom" type="gml:MultiSurfacePropertyType" nillable="true" minOccurs="0" maxOccurs="1"/>
</xs:sequence>
</xs:extension>
</xs:complexContent>
</xs:complexType>
</xs:schema>

So given that it pulls the layer from the config file and the correct DB columns, I know I'm close. But I don't see what's missing.


texts:
I try to add labels to polygons in ArcMap 10 based on a field in the corresponding attribute table.

However, the labels always appear inside the polygon, and I can't find an option that lets me put the labels beside the polygons. How do I achieve this? Thanks.


texts: How to set a specific ArcSDE version for tabular reporting?

texts:
I have created a multiversioned view in Oracle Database 11g Enterprise Edition Release 11.2.0.3.0 - Production. I can access a specific version using a SQL script in SQL Developer. However, using other tools, I can only access DEFAULT.

The following script works fine:

call sde.version_util.set_current_version('MY_VERSION');
select * from mv_my_view;

Encapsulating this logic inside an Oracle View (something similar to Jay Cummins' answer to this question would be a great solution, but I cannot see a way to put a call inside a view).

In Crystal Reports, I can add SDE.VERSION_UTIL.SET_CURRENT_VERSION as a table, but cannot use it to set the version. If I specify a known good version, I get a Query Engine Error: 'Cannot obtain error message from server.'

In contrast, if I specify a non-existent version (with a typo, for example), I get an 'HY000:[Oracle][ODBC][ORA]ORA-20126: Version x not found... ORA-06512'

How can I set the current version to something other than DEFAULT for tabular reporting purposes - either in an Oracle View, or using Crystal Reports, or using some other tool?


texts: 6 degree to 3 degree Transverse Mercator transformation

texts:
I am looking to plot coordinates that were given to me but I am having some issues projections.

I currently have my data projected to NAD83 MTM 1, but the coordinates that were given to me are based in a 51 Degree West longitude of the 6 Degree Transverse Mercator Projection, NAD 83.

Using ArcGIS, what projection would I use to do the transformation.


texts: Editing shapefile in ArcMap programmatically?

texts:
Is this something that can be done?

For instance if I created an overlay on a map and saved it as a shapefile, I want to go and edit this (pre-existing file.) overlay... Either add points/remove points change color/fillstyle/etc.

Ideally I want to avoid creating a new shapefile and deleting the old one.

Thoughts?

Edit 1:  OK, now that someone has said it can be done.  The obvious question is, how do you edit the shapefile (in ArcObjects pro-grammatically)?

Edit 2: Looking at comments it appears I need to do something like this, but I'm not sure.

I currently have a basic shapefile that I want to edit in ArcObjects in C#.

I've created a shapefile by creating a FeatureClass , Workspace and then obviously passing in my shape.  In this case it's just a simple Polyline with a few points on it.  What I'd like to do is have the ability to update this shapefile.

I'm assuming I need to re-use that FeatureClass, Workspace and pass in a new shape.  However, what I'm doing doesn't seem to work.  This is where I'm at on the edit - any advice would be great.

IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)inWorkspace;
workspaceEdit.StartEditing(true);
workspaceEdit.StartEditOperation();

ComReleaser comReleaser = new ComReleaser();

IFeatureCursor featureCursor = inFeatureClass.Update(null , true);

IFeature feature = null;
while((feature = featureCursor.NextFeature()) != null)
{
   int id = feature.Fields.FindField("shape");
   IFeatureBuffer featureBuffer = inFeatureClass.CreateFeatureBuffer();
   featureCursor.DeleteFeature(); //Delete the old
   featureBuffer.Shape = (IGeometry)inShape;
   featureCursor.InsertFeature(featureBuffer);
}

workspaceEdit.StopEdidtOperation();
workspaceEdit.StopEditing(true);

EDIT 3:
Final Solution

featureCursor = inFeatureClass.Update(null ,true);
//DELETE old shape
featureCursor = inFeatureClass.Insert(true);
//Insert new
workspaceEdit.StopEditOperation();
workspaceEdit.StopEditing(true);


texts:
I am creating a website for a non-profit organization and would like to plot points on a map similar to this: http://www.arcgis.com/home/item.html?id=931d892ac7a843d7ba29d085e0433465

I am currently plotting points on google maps through its API, but it isn't a great solution and I have heard that Google has rights to that point. (Correct me if I'm wrong) Anyway, I was wondering if there was a good way to use the USGS Topo layer on the USGS website to plot points dynamically through javascript with an API or some other means... Is there anyway I can set this up through mapnik, similar to this (http://www.toposm.com/us/index.html)?

On another note I have to be "licensed" in order to access ArcGIS javascript API, don't I? Or is this incorrect?


texts: Shp's will NOT line up  (coordsys) - shps made using GPSeismic sotftware (client made shp's)

texts:
OK, I asked this previously but solved it on my own last time.  Same project, but differant shp files.

I recieved from a client shp files that were made in GPSeismic 2011. Upon recieving the shp files from client, he sent the 5 files (dbf, shp, shx, txt and sp1) as you can see in the screen shot.

The POREC did not have the prj (coordinates) so I put in the correct coordinates that are being used. I've tried the coords that we are currently using, "undefined" coords, even did the same with the surface shp file prj, switched that to 'undefined'.  Nothing seems to work.  The postplot keeps falling into the same place, just to the East even with differant coords.

What is going on with this, I make the coordinate changes in ArcCatalog as I have ArcView 9.3. I have two projects as this is the 2d one, the first one I had no problem with the coords as the client sent the shp files made with the GPSeismic software program.  My Proj Manager could not get this to line up on her computer either.

Waht could be the problem AND SOLUTION?  Could the SP1 file be the problem?

sp1 file below (partial)


texts:
I just recently experimented with CartoDB. I have an existing Leaflet Map with lots of Plugins and dependencies. Just adding the CartoDB.js does not work (I think there is a version conflict, but I am not sure - I think it does not work with L.Control).

Therefore, I was wondering if it's possible to just use the CartoDB Core with leaflet. I only want to show layers I created with CartoDB. Popups for the layer would be super, but not mandatory.

Is this possible somehow?


texts: Is it possible to create a popup window in Openlayers2.12 without setting any coordinates on it?

texts:
I use Openlayers 2.12 and create a popup window like this:

popup = new OpenLayers.Popup.FramedCloud("chicken",
              new OpenLayers.LonLat(latlotArr[0],latlotArr[1]).transform( fromProjection, toProjection),
              new OpenLayers.Size(300, 300),htmlImg,null, true);

I get the coordinates from some images. But some of the images don't have any coords. In this case I just want to display the popup without a specific location.

How can I do that? I tried this:

popup = new OpenLayers.Popup.FramedCloud("chicken",
                     new OpenLayers.LonLat(null,null),
                     new OpenLayers.Size(300, 300),
                     htmlImg,
                     null, true);

But ofcource it didn't work. I get an error: invalid value for variable lat.


texts: Iteratively Filling a Dictionary and Updating a Feature Class

texts:
import arcpy

fiberCable = r'C:\Users\broadhurstn\Documents\ArcGIS\Default.gdb\FIBERCABLE'
fibercableToConduit = r'C:\Users\broadhurstn\Documents\ArcGIS\Default.gdb\FibercableToConduit'
conduit = r'C:\Users\broadhurstn\Documents\ArcGIS\Default.gdb\CONDUIT'

# Populate a dictionary that maps segment ID to the placed/planned date
segment_dates = {}
with arcpy.da.SearchCursor(fiberCable, ['segmentid', 'bv_cable_placed_planned']) as cursor:
    for row in cursor:
        segment_dates[row[0]] = row[1]

# Populate a dictionary that maps conduit ID to a two-item list containing segment ID and placed/planned date
ftcs = {}
with arcpy.da.SearchCursor(fibercableToConduit, ['conduitipid', 'segmentid']) as cursor:
    for row in cursor:
        # Use .get() so that it won't fail if a segment ID isn't in the dictionary
        ftcs[row[0]] = [row[1], segment_dates.get(row[1])]

# Populate conduit 'bv_placed_est' field
with arcpy.da.UpdateCursor(conduit, ['ipid', 'bv_placed_est_1']) as cursor:
    for row in cursor:
        for k, v in ftcs.items():
            if row[0] == k:
                row[1] = str(v[1])
                cursor.updateRow(row)

The above code and feature class table produce a dictionary that fills as follows:
{'daac02d0-bf62-4968-b1cd-0fe36c391b2f': ['55a16ddd-2d1d-443b-9935-ecafd2bc7d3b']

How do I add a second value to the dictionary using the below feature class table?

Notice how the value in the dictionary matches "segment ID" in this feature class. I want to have the the SearchCursor look at the current dictionary, see if there is a value that matches the "SegmentID" field values and append the field "Cable Placed Planned" to it's corresponding dictionary item as a second value.

Expected output for the first item:

{'daac02d0-bf62-4968-b1cd-0fe36c391b2f': ['55a16ddd-2d1d-443b-9935-ecafd2bc7d3b', '3/28/2020']

What I've created so far, which doesn't work correctly:

with arcpy.da.SearchCursor(FIBERCABLE, ['segmentid', 'bv_cable_placed_planned']) as cursor:
    for row in cursor:
        if FTCDICT[str(row[0])] == str(row[0]):
            FTCDICT[str(row[1])] = str(row[1])


texts:
I have read the questions about creating heatmaps but I cannot find the Heatmap Plugin in any of the repositories or using the fetch python plugins, it simply is not there.
Could someone please tell me how to install it.


texts: QGIS - Lakes - Change Polygon Z values, using Single Field Value

texts:
I want to change all the elevation of the polygons of a shapefile using the field value ELEV from each polygon. So each vertex of each polygons as all the same Z (elevation). Anyone knows how to do that?

Please

Steeve


texts: Get map from OpenStreetMap and use it as a background with basemap

texts: I want to take a bounding box, use the OpenStreetMap API and get a tile, then use it as a background in basemap. This should be easy, but I haven't been able to figure out the API.


texts: Parse multiple field in field calculator using python

texts:
I'm trying to fill a blank field depending on the contents of three other fields using python in the field calculator.

Sample code is something like:

def getClass(Text1,Text2,Text3):
  if Text1 is None:
    return Text2
  elif Text2 is None:
    return Text3
  else:
    return Text1

I can retrieve the value of the first field but never the 2nd two. I know this isn't the right way but I can't figure out the nesting levels... Any ideas?

UPDATE: Fixed by stripping out the white space. Credit to this post Using Calculate Field tool to calculate on null fields?

def getClass(Text1,Text2,Text3):
  if Text1.strip():
    x = Text1
  elif Text2.strip():
    x = Text2
  elif Text3.strip():
    x = Text3
  else:
    x = "NA";
  return x


texts: Interpreting result of Geographically Weighted Regression (GWR)?

texts:
I want to use the Geographically Weighted Regression (GWR) to model local relationships between my dependent variable and a set of independent variables.

When running GWR in ArcGIS, the coefficients with the parameter estimates can be mapped, which is also recommended. But I am not sure, how to interpret these values correctly.

The following graphic shows a sample (!) how a result could look like. Assuming the dependent variable is auto theft and the visualized independent variable are parking lots. I want to understand how to interpret the parameter estimates.

I do not understand the connection between the parameter estimate values and what they mean in regard to the correlation, meaning if for higher values it can be said that "many parking lots lead to many auto thefts" and for lower values "few parking lot lead to many/few auto thefts".

How can this result be interpreted correctly?

The aim is to find out the local regions, where parking lots correlate with auto thefts.


texts:
I'm having trouble getting labels to appear in a web map, both Leaflet and OpenLayers, using qgis2web.  I can't get either label buffers or backgrounds to show up and even a Bold style only works in Leaflet (not OL).  And to top it off, the placement, around centroid, only works at higher zoom levels.

Are there ways around this if the qgis2web plugin is deficient?  Maybe modifying the output code (CSS or javascript) after the export?


texts:
I am using QGIS 2.18.4 Las Palmas under Ubuntu 16.04. Xenial. At startup, but particularly, when using Orfeo, I get error messages tagged as ERROR 4. My console output. At startup:

ERROR 4: /initrd.img: Keine Berechtigung. Most of all, Orfeo obviously can't save its outputs.

When using Orfeo (Supervised LibSVM classification):

ERROR 4: /tmp/processing0808ccc1678d4ad68ac22abbd3d1a967/4fa15f1ea4d24fb3b868f39e6e85e908/out.tif: Datei oder Verzeichnis nicht gefunden

ERROR 4: /home/andreas/Schreibtisch/output.tif: Datei oder Verzeichnis nicht gefunden

I am fully aware that this may be trivial, however, I haven't been able to solve it anyway.

I could chown on this particular file. However 1. running sudo qgis returns: ERROR 4: `/initrd.img' not recognized as a supported file format. And even if changing the owner solved the problem, it wouldn't help me out regarding the outputs of the toolboxes. For instance, in my third error message, a new output.tif has to be created. For such situations, since files are new, chown woulnd't help. Am I right? Do you have any suggestions?


texts: postgres command to save the line coordinates in openlayers

texts:
I am working on project that we create one map ant editing toolbar. Then by using that editing tool bar we can draw point, line, polygon on the map. Next we are converting the line, point, polygon objects into coordinates using geometry.toString() function. After that i want to save that coordinates in postgres. I am not able to come out from this dig. If you know the route please let me know. I will follow the steps.

Regards,
Kiran


texts:
In QGIS 3.4 I'm trying to add a custom projection that was used for topographic maps of Suriname in the 1950s. I have the parameters from an article about historical CRS used in the country. Using the proj documentation I defined the custom CRS. I even tested it and the returned coordinates confirm that the projection does what I want. However, when I click OK to add the custom CRS, I get the popup: "The proj4 definition of 'Roussilhe Suriname' is not valid."

Here is the proj string:

 +proj=rouss +lat_0=4.11667 +lon_0=-55.68333 +k_0=0.9996 +x_0=300000 +y_0=775000 +ellps=bessel +units=km +no_defs

How can I test it, but then get the message that the definition is invalid? In any case, I now very little about proj4 and can't find the error. Below is a snip of the article containing the parameters.


texts: Create a single line from a resultset of line segments

texts:
How can I create a single line using the line segments geometry found in each row in resultset?

I followed this tutorial: http://anitagraser.com/2013/07/06/pgrouting-2-0-for-windows-quick-guide/

I found roads in a route from point A to point B.
This roads are rows in a resultset and now I need to create a single line from point A to point B representing the entire route.


texts: How to choose a language for web map development?

texts:
I am building a web map application using:

GeoServer
OpenLayers
Tomcat

I would like to build a web application, where the user can log in and have the map arranged according to their profile. As I am already using tomcat, I would like to stick with java.

However, I have found that python is more suitable to web programming, so is jython an alternative?

This thread provided some light, but is not complete.

Are there any examples?
I have thought of using JSF and implementing wms service from geoserver. Any ideas and suggestions to get started. I am ready to learn new technologies.

Features:

  User Should be able to log in to an application.
  User Should be able to only view maps that their role/user permits.
  The application should only make use of the technology listed above, if possible.


texts:
With QGIS3 came python3.6 (which are really nice!),
I'm just curious of the plan ahead, python3.7 is schedule to be released the 15 of June, is QGIS going to follow the python releases? Or will 3.6 be the python version for the next couple of years until "QGIS4"?

Where can I read more on the plan ahead for QGIS (like what's aiming to be implemented in 3.2)?


texts: Converting raster image to feature in GEE

texts:
I am new to GEE, and aim to use reduceToVectors() to convert pixel-value = 1 in a binary image into a featureCollection, calculate the areas of all polygons and finally select polygons above a certain size.
What are the best way to do this?


texts: Where can I find the 1-band color raster plugin for QGIS?

texts:
Possible Duplicate:
  How to install “1-band color raster color table”?

I'm searching for the 1-band color raster plugin. What is the essential repository? Or is there another tool for creating raster maps?


texts: GeoTIFF image draping in Leapfrog Geothermal

texts: I am trying to drape a GeoTIFF image onto my topography. I used manual georeferencing since my image and topography have different CRS. But my image keeps tilting. The referenced points are correct but my image won't lay flat, and I have no idea why or how to fix this.


texts: How do I create a symbol on qgis that has a number assigned to it? I am trying to make a square symbol with a number inside. everytime I place the marker I would like it to ask me what the number is everytime so I can place specific number symbols in different areas on my map. Anybody know how to do this?


texts: Is it possible to sort the features by an attribute programmatically?

texts:
I have written a PyQGIS script that loops through an attribute table and finds duplicated values in a target attribute field.

My technique requires that the attribute table is first sorted by the target field.

Is that possible with PyQGIS? I've only been able to have the table sorted by the row number.


texts: Spatial Join a List of Feature Class Files from Seperate Directories

texts:
I am trying to join two lists of feature class files from seperate folders into a new folder. I've been successful spatially joining a single file but cannot figure out how to list files from folder 1, files form folder 2 - and then join them together. He is what I have tried so far:

# Import Packages
import arcpy,os,sys,string
from arcpy import env
from arcpy.sa import *

#Call the Spatial license
arcpy.CheckOutExtension("Spatial")

#Define the outfolder
outFolder = r"C:\OutFolder"

#Set working environment and list files from folder one
env.workspace = r"C:\Folder1"

#List Files from Folder 1
featureclasses = arcpy.ListFeatureClasses()
for inFeature1 in featureclasses:
    print inFeature1

#Set working environment and list files from folder two
env.workspace = r"C:\Folder2"

#List Files from Folder 2
featureclasses = arcpy.ListFeatureClasses()
for inFeature2 in featureclasses:
    print inFeature2

    #Create file names for tne joined file
    outFeature = outFolder + "/" + inFeature2
    arcpy.SpatialJoin_analysis(inFeature1, inFeature2, outFeature, "JOIN_ONE_TO_ONE", "KEEP_ALL", "", "INTERSECT", "", "")

print "done"

Here is the error I am getting:

  Target Features: Dataset 20180201.shp does not exist or is not
  supported Failed to execute (SpatialJoin).

I think the error has something to so with the filenames possible being the same between the two folders. At the very least, I know the ordering of my lists and the application of the function are incorrect.

Any suggestions?

New error after code update:

  Target Features: Dataset \wrrs1\Student Homes\Group
  Projects\TRMM_GPM_Showdown\DiffFiles\Joined Final\20180201.shp does
  not exist or is not supported Failed to execute (SpatialJoin).


texts: Problems creating TIN from British Columbia DEM data

texts:
Is anyone familiar with creating a TIN using ArcGIS9.3 and British Columbia DEM data (.GEN files)?
I can't figure out how to get the breaklines to work during the creation of the TIN.  Using Workstation/AML it works perfectly.
I can provide more details if anyone's familiar with the problem I'm encountering.  See the images below.

Thanks.


texts: GRASS 7.6 v.net.iso not providing output in QGIS 2.18

texts:
There are simply two layers in EPSG:31469 both are located on the Server side, see image below:

points 'HST_Bestand_GK5'
polylines 'vwESKN_1_25000_Teil'

On my Windows 10 with QGIS 2.18 and GRASS 7.6, I am trying to execute the v.net.iso algorithm. However, there are no outputs, logs, errors etc. What could be a problem?

GRASS 7.6 setting looks as following

I tried what was written in this thread, but my Log console looks okay, see below. Any ideas?


texts: Error memory on Viewshed Analysis plugin QGIS

texts:
My question is exactly the same as TTome! (Viewshed analysis plugin on QGIS keeps resulting in errors)

"I've been trying to test the Viewshed analysis plugin on QGIS 2.10, but it keeps giving back error messages. I've tried to reduce the radius area and have also reprojected both my DEM and my point shapefile to UTM, in order to avoid a possible memory error mentioned in another question regarding this plugin (that the using of geographic coordinates might cause a memory error).
Still, I keep getting Python error messages when I try to run it. Does anyone know whats going on? Is this plugin "broken"?"

Error message:


texts:
i have a layer of points, and a vector layer of county boundaries for the UK. I would like to be able to find the nearest county boundary to any given point and the distance using QGIS. any ideas?

I've installed processing, and trying to use those tools -- but no luck so far. I tried converting the polygon to nodes, but the nodes are spaced too widely for this to be accurate.


texts: How to generate colorful images using GDAL?

texts:
I am using gdal_translate -of JPEG -scale -co worldfile=yes ndvi.tiff output.jpg command to convert tiff images to jpg and resize it. But it generates only black and white images.

How can I generate colorful images when converting from tiff to jpg/jpeg?


texts: Reshaping polygons while preserving area

texts: I'm looking to manually skew the shapes of areal polygons, but need to ensure their areas are not altered.  Is anyone aware of any GIS tools or programming workflows that exist for manually reshaping spatial polygons while preserving areas?  I wonder if there is an approach that utilises anything like the force-directed algorithms used in graph drawing.  Grateful for any advice.


texts: Catch error of empty response of Planet API

texts:
I am automating the download of some small polygons of PlanetScope4 images.

It works well, but some times, at a random polygon, it breaks down with this error:

Traceback (most recent call last):
  File "planet_api/search_asset.py", line 105, in search_images_in_interval
    json_res = search_result.json()
  File "lib/python3.6/site-packages/requests/models.py", line 892, in json
    return complexjson.loads(self.text, **kwargs)
  File "lib/python3.6/json/__init__.py", line 354, in loads
    return _default_decoder.decode(s)
  File "lib/python3.6/json/decoder.py", line 339, in decode
    obj, end = self.raw_decode(s, idx=_w(s, 0).end())
  File "lib/python3.6/json/decoder.py", line 357, in raw_decode
    raise JSONDecodeError("Expecting value", s, err.value) from None
json.decoder.JSONDecodeError: Expecting value: line 2 column 1 (char 1)

It looks like an empty response, not a response of empty features as I find normally.

I tried to catch it like this:

while True:
    try:
        search_result = requests.post(
            PLANET_API_SEARCH_URL,
            auth=(PLANET_API_KEY, ''),
            json=search_request)
        json_res = search_result.json()
        break
    except json.decoder.JSONDecodeError:
        logger.critical("API failed")
        time.sleep(1)
logger.debug(json_res)

But the error still happened.

So the questions are: How is the proper way of catching this error? Any idea why this error happens and how to prevent it so I don't need to catch it?

EDIT:
I just saw that I'm using a not recommended method anymore in this post Planet API: Download a Subarea vs. Clip Images to AOI.

Will my problem be solved if I change to the Subarea method?

I know there are 3 questions in the same post, but I think all of them point to solve the same specific problem.


texts: How to get the tiles on the tileload event and store it locally ?

texts:
I am using openlayers3 to make a small mapping demo , now what i want to do is when the user visits the map demo the first time , ofcourse the tiles are loaded from over the internet , so the 1st time around i would really like to store the map tiles locally , so that when the user visits the map the secound time aronnd the tiles are loaded from the localstorage of the users device rather than from over the internet . So far i have made THIS demo .

Below is how i am initilizing the map:

var map = new ol.Map({
    layers: [
        new ol.layer.Tile({
            source: new ol.source.OSM()
        })
    ],
    target: 'map',
    controls: ol.control.defaults({
        attributionOptions: /** @type {olx.control.AttributionOptions} */ ({
            collapsible: false
        })
    }),
    view: view
});

Below is how i get notified, everytime a tile is loaded,

var source = map.getLayers().item(0).getSource();

source.on('tileloadend', function (evt) {
    console.log('loaded ...');
    console.log(source.getUrls());
    //console.log(tileUrlFunction(evt.tile.getTileCoord(), 1, ol.proj.get('EPSG:3857')));
});

If you check the url i have posted , you will see  console.log('loaded ...'); in the console.  But now what function/method do i use to get the tile , so that i can store it locally ? I.E. store the tile locally ?


texts: How to set initial basemap in web appbuilter developer edition?

texts:
I wan to set a basemap in my basemap gallery. This is my web appbuilder developer edition configuration file.
{
  "theme": {
    "name": "LaunchpadTheme",
    "styles": ["default", "style2", "style3"],
    "version": "2.4"
  },
  "portalUrl": "https://arcgis.swissmediacast.ch/portal",
  "appId": "",
  "authorizedCrossOriginDomains": ["geodata.state.nj.us", "ahocevar.com", "wms.geo.admin.ch"],
  "title": "Wallis",
  "subtitle": "",
  "keepAppState": true,
  "logo": "",
  "geometryService": "https://utility.arcgisonline.com/arcgis/rest/services/Geometry/GeometryServer",
  "links": [],
  "map": {
    "3D": false,
    "2D": true,
    "position": {
      "left": 0,
      "top": 0,
      "right": 0,
      "bottom": 0
    },
    "itemId": "3f2202ad86704ca595280484910dfb2d",
    "mapOptions": {
      "extent": {
        "xmin": 2491135.701601292,
        "ymin": 1060085.625618667,
        "xmax": 2755886.37406769,
        "ymax": 1269336.157114772,
        "spatialReference": {
          "wkid": 2056
        }
      }
    },
    "id": "map",
    "portalUrl": "https://arcgis.swissmediacast.ch/portal/"
  }
}

And I have a basemap config json.
{
  "basemapGallery": {
    "showArcGISBasemaps": false,
    "mode": 0,
    "basemaps": [
      {
        "title": "my map",
        "thumbnailUrl": "images/img.png",
        "layers": [
            {
                "url": "...",
                "id": "m1",
                "type": "WebTiledLayer",
             }
         ]
      }
     ]
   }
}

So I want to start application using m1 basemap.
How can I do it?


texts: Maximum GPS-based location Sampling rate of smartphones

texts:
I am working on Android application that needs to take GPS locations at the rate of 10Hz to 1Hz. I tried it by setting 100ms for time interval for requestLocationUpdates(_, 100, _, _) function. But GPS reading is changed after 10 consecutive reading which means the update is available every  second. Does it mean I can't get location update at higher frequency than 1Hz? The smartphone I am working is Samsung Galaxy A5(2017), Model SM-A520F.
Similar questions are asked at the following links but I couldn't get exact answer.

getting the highest GPS update rate from the GPS Hardware in my
Android
Nexus 5 gps update rate
What is the maximum sampling rate for GPS on Android
GPS Refresh Rate


texts:
Trying to use the following code with an ArcGlobe Control

Type t = Type.GetTypeFromProgID("esriFramework.AppRef");
System.Object obj = Activator.CreateInstance(t);

Howevewr, t is coming back as null.  Any ideas as to why?  I do have the ESRI.ArcGIS.Framework reference in my workspace.  I was hoping to access this so I could then use some ArcGlobe Commands via ICommandBar which I wanted to access via

IApplication myApp = obj as IApplication;
ICommandBars myCommandBars = myApplciation.Document.CommandBars;
myCommandItem = commandBars.Find("ID" , false);
myCommandItem.Execute();


texts: Importing NAS Data into SQLite database?

texts:
I'm trying to import German NAS-Data into a SQLite- Spatialite-database. For this I use the following command:

ogr2ogr -f SQLITE --config OGR_SQLITE_CACHE 8192 -gt 65536 -skipfailures -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=500000 +y_0=0 +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs" -t_srs EPSG:31468 data.sqlite data.xml

This takes quite some time, but finally a SQLite database will be produced. But if I try to open it with QGis I'm told that QGis faild to check metadata. I'm not sure what is missing in the above command.

For dataprotection reasons I cannot provide real data - Testdata can be downloaded from this (official) link: https://ldbv.bayern.de/file/zip/2407/Testdaten_Lurchingen_ALKIS.zip

Can anybody help?


texts: Filtering across image collections based on a minimum pixel value

texts:
I am using the MODIS Terra Continuous Vegetation Fields Global Yearly 250m dataset to derive PTC values for forests. For this, I need to filter pixels that have PTC values that are at least 10%. While I know how to do it on individual images, I am not sure as to how I can apply the same condition for a pixel throughout the image collection (10 years).[for all 10 years, the pixel must be at least 10%] I've tried the filter.Metadata() but I don't think the dataset has a metadata property for me to use. Any ideas as to how I can go about implementing the condition?
var PTC = ee.ImageCollection('MODIS/006/MOD44B')
  .filterDate('2001-01-01','2011-01-01')
  .filterBounds(MYR) //you could use a sample region
  .select('Percent_Tree_Cover')
  .map(function(image) { return image.clip(MYR)});

print(PTC);

var PTC_forests = PTC
  .filterMetadata('system:bands', 'not_less_than', 10); // not too sure which property I should use

Will mapping a function work instead? I believe this might only filter out pixels that fulfil the criteria within the image, rather than applying the same criteria on a pixel throughout the 10 years.
var PTC_forests = PTC.map(function (image){
    var forest = image.gte(10)
     return image.mask(forest)
});


texts:
I'm sorry, I know this is a pretty simple question but am having trouble getting this to work.  Basically, I have 2 linear featureclasses and they both have the same fields and spatial reference.  FC1 has 675 records and FC2 only has 505 records.  So how can I check what Prop_ID(s) are missing from FC2?  I have some code that will copy the missing record from FC1 to FC2 but must first determine the missing Prop_ID.

many thanks


texts: How much of territory do buffers cover using QGIS?

texts:
Is there a way to obtain how much % of the country covers that buffers using QGIS?
This is a country (Ecuador)

And these are buffers


texts: Where is the default QGIS Point layer style .qml file installed in Win 32?

texts: Where is the default QGIS Point layer style .qml file installed in a Windows 32 bit installation? I have done a search on *.qml and can only find four python plugins styles.


texts: How to find layer names within vector tiles, without TileJSON or the .mbtiles file?

texts:
I'm trying to use an endpoint in Mapbox vector tile format. I can access the .pbf's ok, but there's no documentation or TileJSON endpoint.

How can I find the names of layers that exist within the tiles? That is, the value that I would set to source-layer in a vector source in Mapbox-GL-JS.

Theoretically this must be doable, but I don't know any tool that lets me achieve this.

(In this particular case, it's a WMTS endpoint on GeoServer, but I'd be more interested in solutions that work generally.)


texts:
I'm having a lot of trouble with GDAL. Besides just the at times lacking documentation, there seems to be little support in Python. Anyways this is just a batch reproject script that goes thusly:

source_file = gdal.Open(filepath.encode("ascii"))
source_wkt = source_file.GetProjectionRef()
source_srs = osr.SpatialReference()
source_srs.ImportFromWkt(source_wkt)

reproj_file = gdal.AutoCreateWarpedVRT(source_file, source_wkt, dest_wkt)
gdal.ReprojectImage(source_file, reproj_file, source_wkt, dest_wkt)
reproj_attributes = reproj_file.GetGeoTransform()

driver = gdal.GetDriverByName("GTiff")
dest_file = driver.CreateCopy(outputpath.encode("ascii"), reproj_file, 0)

The dest_srs and des_wkt part are not defined in this chunk of code but it's somewhere outside the loop (since it's only needs to be defined once). It seems to work once, I can get one nice looking tif out of it, then give me an 'ERROR 6 WriteBlock() not supported' and python crashes. They are all GeoTIFFs created in the same way with the same basic data (just different times).

Also due to the nature of the reprojection (from GCS to PCS), the AutoCreateWarpedVRT tends to create a lot of blank space, but gives it a value of 0, which is an issue since that could be a real data value. Is there someway to set the nodata value to -99 instead?


texts: Altas with rules on different layer than the cover layer

texts: I need to generate an atlas based on a polygon shapefile, but I'd like to see just what is inside the polygon itself. Maybe with rules that involves other shape's attributes. Any suggestion?


texts: Add text in the body of the a polygon in googlemap

texts: I've created some polygons on a map and wish to have the labels displayed with them .How can i add a name or label to each polygon in the google map?


texts: Unable to see the featureLayer but it is present on the map

texts:
I have created a featureLayer from featureCollection and I am able to add it to the map and when I execute, map.graphicsLayerIds, it also returns the id of the featureLayer, ["resultLayer"]. But I am not able to see it on the map. I have also set the renderer but still it makes no difference. This is the application.

It takes some time to execute (around a minute). This is the code I am using to create the featureLayer,

function addResults(results) {
        //console.log(results);
        var features = results.value.features;

        var layerDefinition = {
            "geometryType": "esriGeometryPolygon",
            "fields": [{
                "name": "FID",
                "alias": "ID",
                "type": "esriFieldTypeOID"
            }, {
                "name": "STATE_NAME",
                "alias": "State Name: ",
                "type": "esriFieldTypeString"
            }, {
                "name": "NAME",
                "alias": "County Name: ",
                "type": "esriFieldTypeString"
            }, {
                "name": "value",
                "type": "esriFieldTypeDouble",
                "alias": "Absolute Value"
            }, {
                "name": "RESULT2",
                "type": "esriFieldTypeDouble",
                "alias": "Percentage Change"
            }]
}
        var featureCollection = {
          "layerDefinition": layerDefinition,
          "featureSet": {
            "features": features
          }
        };
        var states_Symbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, //EVENTS
                    new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID,  //EVENTS
                    new dojo.Color([255,255,255,0.35]), 1),  //EVENTS
                    new dojo.Color([125,125,125,0.35]));

        var popupTemplate = new esri.dijit.PopupTemplate({
                title: "{NAME}",
                fieldInfos: [{
                    fieldName: "STATE_NAME",
                    visible: true,
                    label: "State Name: "
                },
                {
                    fieldName: "value",
                    visible: true,
                    label: "Absolute Change: $"
                },
                {
                    fieldName: "RESULT2",
                    visible: true,
                    label: "Percentage Change: "
                }
                ]
            });

        resultLyr = new esri.layers.FeatureLayer(featureCollection, {
          id: 'resultLayer',
          infoTemplate: popupTemplate
        });
        resultLyr.setRenderer(new esri.renderer.SimpleRenderer(states_Symbol));
        map.addLayers([resultLyr]);

        dojo.connect(resultLyr,"onClick",function(evt){
           map.infoWindow.setFeatures([evt.graphic]);
        });

    }

Any idea what I am missing?


texts:
I need help in identifying the feature attributes displaying a pdf file when click in that feature. I see many examples but I don't have success.
I try this code:

command = "firefox";
imagerelpath = "images_test/test_image.jpg";
layer = qgis.utils.iface.activeLayer();
import os.path;
layerpath = layer.source() if layer.providerType() == 'ogr'
  else (qgis.core.QgsDataSourceURI(layer.source()).database()
  if layer.providerType() == 'spatialite' else None);
path = os.path.dirname(str(layerpath));
image = os.path.join(path,imagerelpath);
import subprocess;
subprocess.Popen( [command, image ] );

but it gives me always an error:

WindowsError: [Error 2] The system cannot find the file specified

Basically, I have a shapefile and I want to click in a feature and open a pdf file. I want to develop that in a python plugin.


texts: Polygon fill outside of stroke in print composer

texts: For some reason the polygon fill is outside of the stroke in QGIS print composer. My stroke width is set to 0.18 which I would like to keep constant. Is there anyway to fix so that the stroke goes on top of the polygon fill?


texts:
I have two raster maps (source: http://ghsl.jrc.ec.europa.eu/ghs_pop.php),
one is with data from 1975, one is with data from 2015.

I want to substract the value of each individual pixel from the map of 2015 with each individual pixel from the map of 1975. That will get me the difference in population for each pixel.

How do I do this?


texts: What are differences between ArcGIS Online, ArcGIS for Server and ArcGIS Publisher?

texts: Can someone help me tell the different between ArcGIS online, ArcGIS server and ArcGIS publisher?


texts: how to remove the color formed on the edge of satellite image?

texts: my AOI requires two scene of satellite image but when I try to mosaic the images using ERDAS imagine, the color of the edge is seen forming black line after mosaic! so please how can I remove the image edge color?


texts:
I am running a script in ArcMap that will update a text box. Some of the text that I am updating it with wraps to two lines, while others are on a single line. What I want to happen is that regardless of the new text being a single or double line, is that the text is left-center adjusted in the box at all times. Currently, it looks like it's left-upper adjusted.

The text box was created from the drawing tool bar as a rectangle text box. I've tried playing around with the Size and Position tab of the box properties window by changing the anchor point, but it doesn't change it when the new text is a single line versus double.


texts: Test if a spatial query contains another

texts:
I have two tables : country and regions.
Country contains all the Europeans countries and regions contains some subdivisions of those countries.
I created a small Union of some countries and regions.

SELECT ST_UNION(geom) AS poly
FROM (
    SELECT reg."id", reg."geom"
    FROM "public"."regions" AS reg
    WHERE reg."id" in (1, 10, 12)
    UNION
    SELECT c."id", c."geom"
    FROM "public"."country" AS c
    WHERE c."id" IN (14, 15, 13)) AS geom_union

And now I want to test if a point (for instance (3, 46) roughly the middle of France) is in the result polygon of my union.

So I tried :

SELECT st_contains(
    SELECT ST_UNION(geom) AS poly
    FROM (
        SELECT reg."id", reg."geom"
        FROM "public"."regions" AS reg
        WHERE reg."id" in (1, 10, 12)
        UNION
        SELECT c."id", c."geom"
        FROM "public"."country" AS c
        WHERE c."id" IN (14, 15, 13)) AS geom_union,
    SELECT ST_GeomFromText('POINT(3 46)') AS point
)

But for some reason I don't understand, i get this error :

syntax error at or near "SELECT"
LINE 2:  SELECT ST_UNION(geom) AS poly
         ^


texts: How to perform statistical test for two sets of points?

texts:
(I have asked this question originally on Cross Validated; however, no good answer and someone suggested me to ask the question here).

We know that we can perform statistical test (e.g., t-test) for two sets of measurements (e.g., some distance measurement) stored in vectors A and B easily.

However, currently I want to perform a similar statistical test for two sets of points (for example, with x-y coordinates stored in A and B). May I ask how to do this in Matlab?

Here is a picture showing what I want to do:

'A' stores the coordinates of the red stars, and 'B' stores the coordinates of the blue circles. These points will be some manual annotations from two human operators. The idea is trying to use statistical test to test whether the annotations from these operators are significantly different or not.

N.B. The points are not 1 to 1 mapped. For 1 point in set A, it can have several corresponding point (nearest ones) in set B.


texts: Adding a new CRS to QGIS 3.10 with custom coordinate reference system

texts:
I am using a Shapefile of which the projection is not recognized by QGIS. The .prj file of the shapefile has the following information:
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
I want to add this custom coordinate reference system, but I can't get it to work. My parameters are validated, but the "name" keeps causing an error (see added picture)
Does anyone know how to solve this?


texts: Responsive design OpenLayers - Bootstrap or Foundation

texts:
I'm looking to use ZURB Foundation or Twitter Bootstrap with my map app, but before I pick a platform I was trying to find out if people had done this already. It seems lots of people have had issues making the OL maps responsive. Full screen maps seem to work but my maps are in a container with a header and side nav.

Has anyone successfully made their maps responsive?

We have been using JavaScript on window resize but its horribly slow, so I've removed that code and don't really want to use JavaScript if I can do it with responsive design, or do I need to wait until version 3 of OpenLayers?

Suggestions on RD framework welcomed too.


texts:
I need to use some of the new functionality in the Dev build and just can not install it on our windows 7, 64 bit machines. It works fine on my personal windows 8 laptop but I can't get it to work at work. I have tried unistalling / reinstalling many times and Lisboa is fine but Dev won't work.

I tried the fix in Why is QGIS not launching?
but no luck.

Problem signature:
  Problem Event Name:   APPCRASH
  Application Name: qgis.exe
  Application Version:  0.0.0.0
  Application Timestamp:    4fdce845
  Fault Module Name:    QtCore4.dll
  Fault Module Version: 4.7.1.0
  Fault Module Timestamp:   4d530715
  Exception Code:   c0000005
  Exception Offset: 000ee27b
  OS Version:   6.1.7601.2.1.0.256.4
  Locale ID:    3081
  Additional Information 1: 0a9e
  Additional Information 2: 0a9e372d3b4ad19135b953a78882e789
  Additional Information 3: 0a9e
  Additional Information 4: 0a9e372d3b4ad19135b953a78882e789

This is on multiple machines with Win 7 -64 bit.

I have been trying to install qgis-dev but I keep getting qgis_core.dll not found issues and I have looked online but can't find anything that works. You get the Dev screen (image) for a millisecond and then no output. The paths seem to load fine.

OSGEO4W home is C:\OSGeo4W

C:\>path C:\OSGeo4W\apps\grass\grass-6.4.3RC3\bin;C:\OSGeo4W\apps\msys\bin;C:\OS
Geo4W\bin;C:\Program Files (x86)\Common Files\Intergraph\Grid Analysis SDK\1.0\P
rogram;C:\Windows\system32;C:\Windows;C:\Windows\System32\Wbem;C:\Windows\System
32\WindowsPowerShell\v1.0\;c:\lastools;C:\apps\FME\;C:\Program Files (x86)\Quick
Time\QTSystem\;C:\OSGeo4W\apps\msys\bin;C:\OSGeo4W\apps\Python27\Scripts;C:\OSGe
o4W\apps\qgis-dev\bin;C:\OSGeo4W\apps\grass\grass-6.4.3RC3\lib

C:\>set QGIS_PREFIX_PATH=C:/OSGeo4W/apps/qgis-dev

C:\>rem Set VSI cache to be used as buffer, see #6448

C:\>set VSI_CACHE=TRUE

C:\>set VSI_CACHE_SIZE=1000000

C:\>start "QGIS" "C:\OSGeo4W"\bin\qgis-dev-bin.exe

C:\>pause
Press any key to continue . . .


texts:
We have a dataset of approximately 30 million new points per month (latitude, longitude, and a measurement associated with that point).
So far, we have applied DBSCAN to generate a monthly map such that the plotted points are:

A "centroid" for each group generated by DBSCAN (calculated by us from the points of each group).
A centroid color equal to the mean of the measurement associated with all the points used to calculate that centroid.

The problem is that from now instead of making a different map each month with the new data, we want to make a map that includes the new data.
Therefore, I have one question:
Would it make sense to apply DBSCAN using both the new data from one month plus the previously calculated centroids, as if these centroids were a raw dataset? (would be an incremental approximation)


texts: How to console.log the imagerySet of a layer?

texts:
I found this simple openlayers example online , the map initialization looks like below:

var styles = [
  'Road',
  'Aerial',
  'AerialWithLabels',
  'collinsBart',
  'ordnanceSurvey'
];
var layers = [];
var i, ii;
for (i = 0, ii = styles.length; i < ii; ++i) {
  layers.push(new ol.layer.Tile({
    visible: false,
    preload: Infinity,
    source: new ol.source.BingMaps({
      key: 'Ak-dzM4wZjSqTlzveKz5u0d4IQ4bRzVI309GxmkgSVr1ewS6iPSrOvOKhA-CJlm3',
      imagerySet: styles[i]
      // use maxZoom 19 to see stretched tiles instead of the BingMaps
      // "no photos at this zoom level" tiles
      // maxZoom: 19
    })
  }));
}
var map = new ol.Map({
  layers: layers,
  // Improve user experience by loading tiles while dragging/zooming. Will make
  // zooming choppy on mobile or slow devices.
  loadTilesWhileInteracting: true,
  target: 'map',
  view: new ol.View({
    center: [-6655.5402445057125, 6709968.258934638],
    extent: [-16655.5402445057125, 6609968.258934638,0,6809968],
    maxZoom: 20,
    minZoom:13,
    zoom: 13
  })
});

var select = document.getElementById('layer-select');
function onChange() {
  var style = select.value;
  for (var i = 0, ii = layers.length; i < ii; ++i) {

    console.log( layers[i].getSource().item(1) +  " " + styles[i] );

    layers[i].setVisible(styles[i] === style);
  }
}

select.addEventListener('change', onChange);
onChange();

Now what i really was trying to do was , i wanted to console.log the imagerySet of every layer , inside the below peice of code:

function onChange() {
  var style = select.value;
  for (var i = 0, ii = layers.length; i < ii; ++i) {

    console.log( layers[i].getSource().item(1) +  " " + styles[i] );

    layers[i].setVisible(styles[i] === style);
  }
}

I tried the below few method:

layers[i].getSource().item(1);

and

layers[i].getSource().getImagerySet();

but it does't give me the desired result , can anybody corrent me and tell me whats the correct way to imagerySet inside the loop ?

FIDDLE HERE


texts: Deleting a common word in attribute table

texts:
I have a dense attribute table of all the districts in the UK. I want to label the districts but each row in my 'NAME' field has 'District' after it. I wish to take this off, only leaving the name of the district e.g, Wycombe, South Bucks, Chiltern etc...

Is it possible to create a new field by writing an expression to eliminate all values that match 'district'?


texts:
For a dev environment, I created a fake Certificate Authority using OpenSSL.
I first generated a key and certificate PEM then used the files to export a PFX. I then created a server key and request file. I signed the request file with the CA -adding a config file with subject alternative names- that I had created and exported a PFX.
For the next steps I added the CA to my trusted root of CAs on windows, added the server PFX in IIS to my site, then bound port 443 to the certificate.
Despite this, I still get the error that the server could not prove its identity because of its name.
Am I missing something? Do I need to do use the same server certificate I created to replace the one by default in server admin?
Note that I am operating on a VM in Azure. I added localhost and the server hostname on Azure to the list of SAN.


texts: I work with ArcGis Explorer and the coordinate system default is WGS 1984. Is it possible to change the coordinate system?


texts:
When I try to load qgis within an external Python instance, I get the following error:

>>> from PyQt4.QtCore import *
>>> from PyQt4.QtGui import *
>>> import qgis
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "E:\QGIS\apps\qgis\python\qgis\__init__.py", line 27, in <module>
    from qgis.core import QgsFeature, QgsGeometry
  File "E:\QGIS\apps\qgis\python\qgis\core\__init__.py", line 30, in <module>
    from qgis._core import *
ImportError: DLL load failed: The specified procedure could not be found.

QGIS version: 2.16.1.

Python in QGIS:

>>> print(sys.version)
2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]

Python installation:

>>> print(sys.version)
2.7.12 (v2.7.12:d33e0cf91556, Jun 27 2016, 15:24:40) [MSC v.1500 64 bit (AMD64)]

Paths in Python installation:

>>> print(sys.path)
['', 'E:\\QGIS\\apps\\qgis\\python', ...

Additionally the following paths are set in the environment of Windows:

E:\QGIS\apps\qgis\python;E:\QGIS\apps\qgis\bin;

Any idea how to fix this? Or how shall I setup my system to be able to access QGIS "interactively" from external Python IDEs ?


texts:
I use data driven pages to produce several maps. In the end product, all maps should show point data on a Country of europe.

So far the country is shifting around each time I change the data shown in the map. However, I want the country to stay in the middle of the map when I export the maps with a scale of 1:750'000. I tried all Options in the Data Driven Pages Set Up (best fit, Center and maintain current scale and data driven scale).

What settings do I have to use so that the Country will stay in the Center of the map? (I work with ArcGIS 10.2.2)


texts: WKT MultiPoint string not working in GeoServer?

texts:
I have a problem with a MULTIPOINT WKT string (SLD file for GeoServer) that is not displaying at all:

<WellKnownName>wkt://MULTIPOINT((-1.5 1), (0 1), (1.5 1), (-1 -1), (1 -1))</WellKnownName>

I have also tried without nested parenthesis but same issue:

<WellKnownName>wkt://MULTIPOINT(-1.5 1, 0 1, 1.5 1, -1 -1, 1 -1)</WellKnownName>

The features are not visible at all, whereas they are with a simple WellKnownName parameter such as circle.

Here is the symbol I need to produce:


texts:
I have QGIS 2.4.0. Until just today, the "Delete ring" feature worked just fine. Now it no longer works at all. I rebooted and even reinstalled QGIS 2.4, and it didn't fix the problem.

But I also still have QGIS 2.2.0. "Delete ring" works fine on that.

Why did "Delete ring" on QGIS 2.4 stop working?


texts: Can't put a basemap in its Precise location

texts:
I am a total beginner with ArcGIS.

( My software info: ArcGIS Desktop Version 10.1 ) (Windows 8.1 Pro).

I want to use a raster as a Basemap for my city, I know the coordinates of my city. my question is how can I Set this raster in its Precise location? I tried to set up the exact coordinates of my city but when I try to project a point shapefile to it, the points shows up somewhere ELSE and not on the raster.(even though they have the same coordinates of the raster which is my city's coordinate).


texts:
I'm trying to use OSM maps in a simple html page, just to learn the bases and evaluate if it will allow us to represent specific data on a map.
Our needs are to have a web map of a specific highway and a vector representation of the two directions, bifurcations, exits and so on. The top would be to get a representation of the different lanes to mark work in progress along the road, but I don't understand if I can get to a so deep level.

As a simulation, I downloaded the highway osm file map from OSM map site and edited it, using osmosis, removing all the data I don't need.

I copied the example at https://wiki.openstreetmap.org/wiki/OpenLayers_osm_file_example in a local file and updated the lat and lon coordinates and the osm file name. What I obtained was the normal map of open street map without any of the content of my osm file. Is there something more I need to do to get it works?

The example above uses Openlayer 2, so I thought probably it is too obsolete to properly work, then I took a look at Openlayer 5 vector example. I simply copied and pasted the code on a local html file, without changing a word, and what I obtained was a perfectly blank page and the browser console was full of complains due to a not correct syntax. So or the example on line was really full of errors or I'm missing something...

I really don't understand where I'm wrong in both cases.

Can anyone give me a suggestion?


texts: Different SRID shown in QGIS and PostGIS for the same data

texts:
I have data stored in PostGIS. When I view its properties in QGIS, it says:

SRS: EPSG:4326 - WGS 84

However, when I want retrieve the SRID value directly using the Find_SRID function, it returns 0. And by looking at the data, SRID is really 0 in PostGIS.

Query:

select * from geometry_columns;

Returns:

f_table_name | f_geometry_column | coord_dimension | srid |  type

-----------------+------------------+--------------+-------------------

vector       | wkb_geometry      |               2 |    0 | POINT

vector       | wkb_geometry      |               2 |    0 | POINT

buff_out     | wkb_geometry      |               2 |    0 | POLYGON

buff_out     | wkb_geometry      |               2 |    0 | POLYGON

Question: Why is SRID for the same data different in QGIS and in PostGIS? How do I get the correct SRID value from PostGIS?

I am using PostGIS 9.5.


texts: How to use Select by Attributes & LIKE operator to find a %

texts:
I have a large data set of points which has the field 'Name'.

Some of the points have an attribute such as 'Slope >40%' or 'Grade 10-30%' for example.

I would like to use Select By Attributes and the LIKE operator to select all points which have a % within the attribute.

I've tried different combinations such as:

WHERE NAME LIKE '%[%]%'

But no luck on trying several combinations...


texts:
I am a new user in PostGIS.
I need to create a table with the ID_DIR of table A, ID_SQ of table B and have a new column called "Distance" that would be the distance between the ID_DIR and the ID_SQ in a zone of influence of 500 meters of each ID_SQ.

CREATE TABLE RESULT AS (
   SELECT a.ID_DIR,
          ST_AsTEXT(ST_Shortesline(a.geom,b.geom) AS geom,
          a.ID_SQ,
          ST_Length(
            St_Shortesline(a.geom,b.geom)::geography/1000) AS Distance
   FROM TABLE_A AS a,
        TABLE_B AS b
   WHERE ST_Dwithin (a.geom,b.geom,500)
   ORDER BY a.ID_DIR,Distance ASC
);

I have two problems:
1) This query does not work
2) I have seen that the processing time in these examples is very high. Table A has 4 million data and it is duplicated depending on the times it crosses with the ID_SQ buffer
Help!


texts: Clipping single tile layer with georeferenced polygon

texts:
I have very similar question to this one. The difference is I would like to clip only one layer that is not a base layer. I like the solution purposed in mentioned question, but unfortunately it clips all layers.

Does any one have any clue how to achieve this?


texts: QGIS Join attribute by most similar lines

texts:
Using QGIS 3's (new?) Join attributes by nearest feature on two multiline layers, I get unexpected results... Namely, if a line crosses another, that is favored over one that is almost exactly the same, but shifted slightly:

While this seems like unexpected results to me, I may be doing it wrong... is there a "Join attributes by most similar" tool that will get what I'm looking for and match the two parallel lines over the two perpendicular ones?


texts:
I'm working with QGIS 2.16.2 and I'm trying to use the MMQGIS hub lines plugin from the python console.
For some reason the output of the process only connects a few point in the hubs shapefile with a few point in the spokes shapefile.
Both shapefiles share point locations. 'centroids_distr1' contains a 'from_district' column with unique values from 1 to 12 and 'centroids_distr2' contains a 'to_district' column with repeated values from 1 to 12.

I'm new using python in QGIS so I don't know if it's a bug or it is that I am doing something the wrong way.
This is the code lines I'm using:

import processing
processing.runalg('qgis:hublines', 'C:/BUS/centroids_distr1.shp', 'from_district', 'C:/Documents_AGP/TUS/centroids_distr2.shp', 'to_district', 'C:/BUS/HubLines_TEST.shp')

Is there anything wrong?


texts: Returning a point by a specific pixel value or range using Google Earth Engine

texts: I have an image classified with pixel values equals to one and zero. I want to return the lat and long from pixels values equals to one, would be possible, using GEE?


texts:
Using code I found at Go through items in a composer legend (QGIS 2.6) I can access the  the layers in my legend fine. However, I have layers that are styled by category and I can't figure out how to change the category labels. I just want to replace the underscores in the category labels with a space.

After looking through the API I think I may be able to use the rendererV2 class. The code below is giving strange results. the first set of print commands return 'foo', but the second set return the original labels and these are also printed on the output pdf.

Is there something I'm missing with the renderer? or, is there a better way to achieve this task?

I'm using QGIS 2.18.7 on windows 8

    for legendLyr in iface.mapCanvas().layers():
        if legendLyr.name() != "os1250_line" and legendLyr.name() != "os1250_text":
            renderer = legendLyr.rendererV2()
            if renderer.type() == "categorizedSymbol":
                myRenderer = renderer.clone()
                for cat in myRenderer.categories():
                    cat.setLabel("foo")
                    print(cat.label())
                legendLyr.setRendererV2(myRenderer)
                for cat in legendLyr.rendererV2().categories():
                    print (cat.label())
                legendLyr.triggerRepaint()


texts: Install QGIS 3 and QGIS 2.0 Ubuntu 18.04

texts: I have successfully installed  QGIS-3 on my ubuntu 18.04. I want to install QGIS-2 as well along with QGSI-3. How can i install QGIS-2.18 as well. So that QGIS 2 and QGIS 3 both would be running on my Ubuntu 18.04.