texts: Render raster from vector tiles in WGS84 (EPSG:4326) using Tileserver GL

texts:
I have vector tiles that are in the WGS84 CRS (EPSG:4326). I can confirm that they are in WGS84 because on a Web Mercator basemap (pg_tileserv preview) they display in Svalbard (should be in Austria):

When I load them in Tileserver GL using the following style.json, they do show up:
{
  "version": 8,
  "name": "Demo",
  "metadata": {
    "maputnik:license": "https://github.com/maputnik/osm-liberty/blob/gh-pages/LICENSE.md",
    "maputnik:renderer": "mbgljs",
    "openmaptiles:version": "3.x"
  },
  "sources": {
    "public-parcels": {
      "type": "vector",
      "tiles": [
        "localhost:8888/public.declarations_4326/{z}/{x}/{y}.pbf"
      ],
      "minZoom": 0,
      "maxZoom": 14
    }
  },
  "sprite": "https://maputnik.github.io/osm-liberty/sprites/osm-liberty",
  "glyphs": "https://api.maptiler.com/fonts/{fontstack}/{range}.pbf?key={key}",
  "layers": [
    {
      "id": "public-parcels",
      "type": "fill",
      "source": "public-parcels",
      "source-layer": "public.declarations_4326",
      "paint": {"fill-color": "rgba(41, 99, 39, 1)"}
    }
  ],
  "id": "demo"
}

When I try to render the tiles as rasters (requesting the ?raster endpoint), nothing shows up. I can confirm that the same configuration works for raster rendering web mercator tiles.

Has anyone successfully rendered WGS84 raster tiles from vector tiles in Tileserver GL?


texts: How do I create a leader line for a label in my bar graph in ArcMap?

texts: I am displaying a bar graph of various tornado watches and warnings along with the actual tornado touchdowns themselves. Since the tornado touchdowns are such a small number compared to the large quantity of watches/warnings I would like to label the number of tornado touchdowns with their specific number. The problem is that when I display the label (called a mark within the properties window) it covers my points on the graph. I would like to create some type of "leader line" to move the labels away from my points. I have enclosed a picture of something kind of like how i want it - I drew the sloppy little black lines is if they were leader lines and somewhere near the end is where I would like to move the labels so they don't cover my points on the graph. If there are any suggestions I would appreciate it. Thanks a bunch.


texts:
Does anyone out there know how to perform radiometric corrections on a single image taken with a multispectral sensor?

I know there is proprietary software (such as PiX4D) that can do this; however, they can only process images in series (which I don't have and can't obtain).


texts: Qgis: Fill a field depending from 2 others fields

texts:
I want to automatically fill a field based on the information from two others.

I have a column/field 1(integer) with values from 1 to 10...

...and, for each of those ten values in column 1, on the column 2 might appear 100 "different"/possible strings ( made up by associations like in the figure below)

I have to replace those 100, with 100 others,and i want the result to be written on a 3rd  new field/collumn

(for example, if i get 1 in collumn1 i want to have this replacement :"1 and 2" -> "YPC5 plus YPC6"...". And 100+ like this one. Each cell/possibility will have a match, but also each value that made up the association/string will have one and only one possible match. 1 will always be replaced by YPC5.s ( note: observe column 2, bellow)


texts: Using pyshp to convert .csv file to .shp?

texts:
I am trying to understand how I can use the csv module in python to open a csv file in the same folder as the python script, and then create a shapefile using the shapefile module pyshp.

The csv file looks like this, but can have a couple of thousand rows of records:

id_nr;date;target;start_lat;start_lon
1;2012-05-21;navpoint 25x;55.123654;13.456954
1;2012-05-23;navpoint 11f;55.143654;12.456954


texts: Creating Image Mosaic database table in GeoServer?

texts:
When I create an image mosaic store in geoserver, the table is automatically created in the PostGIS database.

If it is possible to alter the table name, how will GeoServer know what table to look to?


texts: Unable to create Dimension Feature Class

texts:
I am trying to create a Dimension feature class within a new geodatabase.

I right-click on the geodatabase and select the New, then Feature Class.

The New Feature Class window comes up and I can put in the Name, but when I use the drop down to select Dimension Features, I only get Polygon, Line, Point, Multipoint, Multipatch, and Annotation Features.

I am using Arcmap 10.2 standard license. As I understand, I should be able to do this.

Is this a bug? Or, am I doing something wrong?

How do I make a new Dimension Feature?


texts:
I am trying to display raster image using MapServer but it is not displaying anything, page is blank

My map file code is:

MAP
NAME "ny_pop"
SHAPEPATH "C:\ms4w\apps\webgis\mpgaps"
EXTENT 8240841.450745 2400319.900386 9218007.463612 3107183.406947
SIZE  450 500
STATUS OFF
DEBUG ON
IMAGETYPE "png"

PROJECTION
    "init=epsg:900913"
END

OUTPUTFORMAT
  NAME "png"
  DRIVER AGG/PNG
  MIMETYPE "image/png"
  IMAGEMODE RGB
  EXTENSION "png"
  FORMATOPTION "GAMMA=0.75"
END

WEB
    METADATA
        "wms_title" "WMS Pop Tif"
        "wms_onlineresource"  "http://localhost:8085/cgi-bin/mapserv.exe?mode=map&map=C:\ms4w\apps\webgis\mpgaps\raster.map"
        "wms_enable_request"  "GetMap GetFeatureInfo GetCapabilities"
    END
END

###########
LAYER
    NAME "ny_pop"
    DATA "094_056_d_26mar2012_516_12-mar-08.tif"
    STATUS OFF
    TYPE RASTER
    METADATA
        "wms_title" "WMS NY Raster"
        "wms_srs" "EPSG:900913"
    END
END

###########

END

Data path is C:\ms4w\apps\webgis\mpgaps


texts: Converting DateTime field to another time zone?

texts:
I'm using QGIS 2.16.3 and attempting to convert a number of datetime fields that are appearing in the incorrect default GMT +0 timezone. I would like to add, for example, 10 hours to accurately reflect the time zone of data collection.

So far I've attempted using the field calculator:

"Date_Field" +   ('10:00:00')

which just adds a string to the end (i.e. '2017-4-06 14:00:001')

I've also tried:

( "Date_Field") + to_time('10:00:00')
( "Date_Field") + to_datetime('10:00:00')
to_datetime( "Date_Field") + to_datetime('10:00:00')

which gives me an error "Cannot convert '2017-04-06 14:00:00' to double", although I am able to convert the datetime field alone.

This seems like a simple, common operation for anyone working with different time zones. My only workaround is to create a CSV and edit the fields in a spreadsheet.


texts: Write referenced GeoTiff with 4 corners on matlab

texts:
I have some matlab raster matrix to export to geotif.
I also have the location of each corner in Decimal degrees (WGS84).
I know how to export and write the data into a geotiff but i'm having a problem to use the 4 corners for geolocation.
The raster is a strip map and it's not oriented to the north (it's slanted).

I can georefrence it afterwards on arcmap but i'm trying to avoid it and add it to the matlab script.


texts: Output raster from IDW interpolation tool is from 3.40282e+038 to -3.40282e+038 in ArcPy

texts:
I am not familiar with Python in ArcGIS, so I used ModelBuilder first to set the environment and put the tool in. The model worked right and the output raster seemed good.

Then I exported the model to a Python script. The raster output is generated, but the value of raster is from from 3.40282e+038 to -3.40282e+038. Why did this happen and how can I fix it?

My script is like this:

# Import arcpy module
import arcpy

# Check out any necessary licenses
arcpy.CheckOutExtension("3D")

# Set Geoprocessing environments
arcpy.env.scratchWorkspace = "C:\\Documents\\ArcGIS\\biomass.gdb"
arcpy.env.outputCoordinateSystem = "PROJCS['WGS_1984_Plate_Carree',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Plate_Carree'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',0.0],UNIT['Meter',1.0]]"

arcpy.env.extent = "60.8785969898536 -11.0048500049116 148.898540012983 53.560734"
arcpy.env.cellSize = "MAXOF"
arcpy.env.workspace = "C:\\Desktop\\shape"
arcpy.env.mask = "newasia.shp"
# Local variables:
xyz_shp = "xyz.shp"
b1_ald22222 = "C:\\Documents\\ArcGIS\\mass.gdb\\b1_ald22222"
# Process: IDW
arcpy.Idw_3d(xyz_shp, "ALDX", b1_ald22222, "27000", "2", "VARIABLE 8", "")


texts: What format to use for highly compressed vector data?

texts:
Does anyone have information about possible formats for highly compressed vector data?

Some basic requirements.

Must be possible to read using the larger commercially available platforms. (ArcGIS, Luciad, MapInfo etc.)
Must have tools that allows us to create files/format from other data formats (i.e. Shape)
Speed when reading the format is important as the data is among other things used as a live map in mobile platforms. Compression speed is less important (although nice)
Must be possible to use a map product in  a wide scale range, a 50k map shall be possible to read at a 1 000 k zoom level
Ability to do selective updates on area, feature-ID or other search criteria would be nice, but not a must.
Memory footprint must be small, it needs to work with Windows XP and/or Windows 7 allocating less than 1GB RAM so some kind of smart cache management must be implemented

To give an idea of the required compression and as input for memory footprint requirements, I need to store approximately 1TB of shape data to a 100GB disk.

Does anyone have suggestions for good products?


texts: Is there an equivalent argument to index_col in geopandas.read_file() function?

texts:
I am looking for the equivalent argument to pandas.read_csv(index_col='column') in the geopandas.read_file() function. My goal is to read a shapefile and specify the index column. However, the following does not set any index column and returns the same GeoDataframe as calling the function without the argument:

import geopandas as gpd

gdf = gpd.read_file('test.shp', index_col='column')

The documentation of the read_file function does not list index_col as an argument and states that one can specify keyword arguments that are going to be passed to the fiona.open function (but not to any pandas.read_* function).

A workaround is to use the set_index() method after reading the file:

import geopandas as gpd

gdf = gpd.read_file('test.shp')
gdf = gdf.set_index('column')

# or

gdf = gpd.read_file('test.shp').set_index('column')

However, it feels like specifying the index should be possible directly from the read_file function.


texts: How to merge raster with FME with different spatial resulution?

texts:
I have two sets of raster (tiff, aerial photos) which I want to merge to one raster file with FME.
The two sets have different spatial resolutions. Set A has a better resolution as B. The two set do overlap.
How can I merge these raster with FME? Requirements:

Use Resolution (size of pixel) from set A
Where the raster do overlap: use pixels from set A
Optional requirement: Position of pixels from set A should not change (nor resampling of A. only resampling of B)

Example with 4 raster files from set A (blue) and 4 raster files from set B (red).

I tried to use FMEs RasterMosaiker and set the order in Navigator but I can not meet requirement 1 and 2 simultaneously:

first A then B: I loss requirement 2
first B then A: I loss
requirement 1


texts: Seeking books for learning ArcGIS Desktop?

texts:
I was wondering what book(s) is/are available for learning ArcGIS Desktop functionality.

I would like to use one in conjunction with books on the theory of GIS since even though I can understand the concepts I still have trouble with how to use the more advanced ArcGIS Desktop features.

I understand there are quite a few online tutorials but I do my best learning if I have a book in front of me.


texts: Checking the opening of an attribute table with pyqgis

texts: I know how to show the attribute table of a specific layer (with showAttributeTable), but how to check if the attribute table of a specific layer is already opened ? I don't want to show it again if it's already visible.


texts:
Are there any available examples that show how python can be integrated into a .NET ArcGIS Addin or Extension? I have a series of python scripts that I would like to call from form events (e.g. button clicks) in .NET. I would also like to be able to set variables in the script from a .NET form.

Should I bite the bullet and re-write these scripts in .NET, or is there a simple way to wrap the python code?


texts: Synchronizing KML file in Google Earth with ArcMap?

texts:
I have a KML file in Google Earth that I keep updated.

Usually every update I save this KML file then I import it to ArcMap.

Is there any way to sync this KML file with ArcMap?

I saw this Link but it didn't work with ArcGIS 10.3 for Desktop (i didn't find the file ESRIRegAsm.exe and I couldn't add it as plugin).


texts: Snap a point to a line and get the line ID (Oracle Spatial or SDE.ST_GEOMETRY)

texts:
I have an Oracle 18c/10.7.1 geodatabase that has:

A sidewalk polyline feature class (existing features)
And a sidewalk inspection point feature class (the features are created via ArcGIS Collector/Portal)

As soon as the sidewalk inspection points are created, I want the database to automatically snap the points to the closest line (and insert the line's ID into the point FC).

It would be ideal if this could happen instantly -- instead than doing it with post-processing.
I'm not interested in cheating by snapping to the line vertices (sidewalk deficiencies/inspections don't necessarily occur at vertices). I need to be able to snap to the closest position along the line (even if the closest position is mid-segment).
The search tolerance would be 5 metres.

Is there a way to snap a point to a line (and get the line's ID) via logic in the database? (Oracle Spatial or SDE.ST_GEOMETRY)


texts: Editing attribute cell data/text in "edit" mode using QGIS

texts:
I am using QGIS3.4.
I turn my layer on "edit" mode.  I want to type text into a blank cell or paste text.  I don't know why it will not allow me to enter anything while the layer is in "edit" mode.  I have been able to do this in the past.
I am fairly new to QGIS so some of the following answers might not make sense due to my lack of knowledge but here we go. I think I may already know what is wrong but I am answering the questions fully to assist other users like me in the future and also I do not fully understand GIS as a whole at this point.
My layers are from multiple origins.  Some layers are shapefiles (via zip files downloads) from my local government agency website.  The layers that I am trying to edit are shapefiles created internally from a previous a software called UDIG.
I am not sure I understand the difference between a file or a database.  I imported shapefiles into a new QGIS project.  Then I saved the QGIS project to my local drive.  From there, I downloaded the QGIS Cloud plugin to create some online maps.
At first, the file was saved locally on my drive.  After using the QGIS Cloud plugin, I believe the file/database is stored basically on QGIS Cloud.  I have lost some data and had problems not knowing how QGIS cloud is handling my data.
I think the file type is all Shapefiles.  I really couldn't explain the differences in each of a Shapefile, a Geopackage, a PostGIS database, a WMS connection etc.  I need to read more on this.  I assume Shapefiles and Geopackage files are basically the same but store my vector data a little different.  I don't know really anything about a WMS connection yet but I assume it is a way to connect remotely to a database of information without storing the information on my computer.  Basically, I open the QGIS file and it opens the WMS connection.
I have included a screenshot of what I am trying to do.

I have included a screenshot of my layer properties / Information tab to show the source of the layer.


texts:
I have a 1/3 arc-second DEM downloaded from USGS that I am trying to use to create a slope analysis layer.

The DEM is in ESPG 4269 to start and I've re-projected it into NAD83 UTM Zone 11N (the DEM is for eastern Oregon) so the vertical and horizontal units are in meters.

The output looks correct except all of the lines through it (see screenshot).  Maybe I went wrong in how I projected it?


texts: Determining Equal-Area Slope of Longest Water Course in Catchment?

texts:
I'm still very new to python and need some direction in developing a ArcPy\Python script to determine the Equal-Area Slope of the Longest Water Course within a catchment. The Longest Water Course is a 3D polyline that represents the the longest river centre line in a catchment. The two jpegs explain the formula that I've trying to regenerate using ArcPy and Python:

I've figured out that I need to append the X-coordinate; Y-coordinate; and Z-value from the geometry into a python list. The I then need to step through the python list to determine the following: Area = [(0 + 1)] X 500 + ....

The 0 is the difference in height based on the first vertex Z value, which allways will be 0, the 1 is the difference in height of the second vertex Z value from the first vertex Z value, in this case 1m, the 500 is the distance between the first and second vertex.

I need assistance\guidance in how to step through the python list and calculate the following and then stick the results back into a python list or dictionary. I can then step through the python list or dictionary to determine the equal-area = [(0 + 1)/2]x 500 + [(1 + 2)/2]x 300 ....


texts:
I got the request to transform a PNG in RD_new to WGS84. Is this possible?

I was thinking of changing the accompanied world file its coordinates. I do not think it is possible to load in a png in a transformation function right?

I normally would use spTransform from RStudio for it but I am expecting it wont work.


texts: Some layers won't load when running application out of QGIS. "Item has already been added to this scene"?

texts:
I'm developing a basic application out of QGIS and am having trouble getting my code to generate all the layers I have set. Currently it only generates two of the four layers: "self.basemap_layer", which is a TIFF raster layer, and "self.pointlayer", which is a point layer generated by a text-delimited CSV. It does not generate the two vector shapefiles left over.

Instead, I get two errors:

QEventLoop: Cannot be used without QApplication
QGraphicsScene::addItem: item has already been added to this scene

I can't figure out why it does connect to two layers, but does not connect to the other two...

My code is below:

def loadMap(self):

    cur_dir = os.path.dirname(os.path.realpath(__file__))

    #DOES LOAD
    filename = os.path.join(cur_dir, "data","NE1_50M_SR_W.tif")
    self.basemap_layer = QgsRasterLayer(filename, "basemap")
    QgsMapLayerRegistry.instance().addMapLayer(
        self.basemap_layer)

    #DOESN'T LOAD
    filename = os.path.join(cur_dir, "data", "ne_10m_roads.shp")
    self.roads_layer = QgsVectorLayer(filename, "roads", "ogr")
    QgsMapLayerRegistry.instance().addMapLayer(self.roads_layer)

    #DOES LOAD
    filename = "/home/me/Desktop/PyQGIS/appfiles/data/thiefsymbolcsv.csv?delimiter=%s&crs=epsg:3857&xField=%s&yField=%s" % (
    ",", "X", "Y")
    self.pointlayer = QgsVectorLayer(filename, "point", "delimitedtext")
    QgsMapLayerRegistry.instance().addMapLayer(self.pointlayer)

    #DOESN'T LOAD
    filename = os.path.join(cur_dir, "data", "line_points.shp")
    self.loblayer = QgsVectorLayer(filename, "lobs", "ogr")
    QgsMapLayerRegistry.instance().addMapLayer(self.loblayer)

    self.showVisibleMapLayers()
    self.mapCanvas.setExtent(QgsRectangle(-127.7, 24.4, -79.3, 49.1))

It's not giving me errors related to my imported PyQT and PyQGIS libraries, and the application window works just as inteded, so I don't think I've set my environment variables incorrectly. But I did read elsewhere that these error might have to do with an incorrect wrapper script, so I'm not entirely sure.

Why are these errors occurring, and how do I get rid of them? Seeing as both the problem layers are .SHP files, and the ones that work are .TIF and .CSV, could their file type somehow be the issue?

My full code is here for reference

UPDATE: After closer examination, I've found that all layers do indeed load, just in the wrong order. My point feature is at the top, followed by the basemap, then the line and road layers on the bottom. Trying to figure out how to change the drawing order. I also moved my "watcher" portion, which was below my csvLobRefresh method, to the bottom of my __init__ method in my MapExplorer class, which got rid of the "QEventLoop" error. the "QGraphicsScene" error remains, however.


texts: Apply graduated symbology renderer works in Python console but not in standalone script

texts:
I use the following PyQGIS script to load a vector layer from a shapefile and apply graduated symbology to it (tried to make it as minimalistic as possible). This works well in the QGIS 2.4 Chugiak Python console - the specified vector layer is loaded and graduated symbology is applied according to the arguments given to the function (vectorLayer, targetField, classes, mode).

However, when using this in a standalone script, I always get

layer loaded succesfully
ERROR 1: Argument of column Summary Function 'MIN' should be a column.
applying graduated symbology
ERROR 1: Argument of column Summary Function 'MAX' should be a column.

and no graduated rendering is applied. Any suggestions why this happens or how I can use QgsGraduatedSymbolRendererV2 in a standalone script (preferably by using one of the existing rules instead of creating all classes on my own)?

Other functions like a more or less selfmade appliance of categorized rendering work well, so this does not seem to be a general issue with my Python/IDE/QGIS environment. I already made sure that my columns are numeric (and otherwise it would not work in the Python console).

Here is the code I use. When trying this out, make sure to adapt the parameters at the top to match your file. When trying it in the Python console, comment/delete the part at the top concerning set up qgis app and the last line app.exitQgis(), otherwise QGIS will probably crash.

# imports
from qgis.core import *
from qgis.utils import *
from PyQt4.QtCore import *
from PyQt4.QtGui import *

import os, glob

###########
# parameters
###########
pathToLayer = '/path/to/your/favorite/vectorLayer.shp'
targetField = 'someNumericColumnInTheVectorLayer'

pathToQGIS = '/path/to/your/qgis/folder/'

###########
# set up qgis app
###########
QgsApplication.setPrefixPath(pathToQGIS, True)
QgsApplication.initQgis()
app = QgsApplication([], True)

###########
# load a layer
###########
def loadVectorLayer(pathToLayer):
    vectorLayer = QgsVectorLayer(pathToLayer, 'ThisIsAVectorLayer', 'ogr')

    # validity check. if layer is valid, add to canvas
    if not vectorLayer.isValid():
        print "Layer failed to load!"
    else:
        QgsMapLayerRegistry.instance().addMapLayer(vectorLayer, True)
        print 'layer loaded succesfully'
        return vectorLayer

###################################
# function to apply graduated symbology to a layer
###################################
def applyGraduatedSymbology(vectorLayer, targetField, classes, mode):

    print 'applying graduated symbology'

    # set up renderer
    symbol =  QgsLineSymbolV2()

    colorRamp = QgsVectorGradientColorRampV2.create({'color1' : '0,255,0,255',
                                                     'color2' : '255,0,0,255',
                                                     'stops' : '0.5;255,255,0,255'})

    renderer = QgsGraduatedSymbolRendererV2.createRenderer(vectorLayer, targetField, classes, mode, symbol, colorRamp)

    vectorLayer.setRendererV2(renderer)
    return

###################################
# main
###################################
def main():

    vectorLayer = loadVectorLayer(pathToLayer)

    applyGraduatedSymbology(vectorLayer, targetField, 3, QgsGraduatedSymbolRendererV2.Jenks)

    print 'finished'

    return

main()

app.exitQgis()


texts: Calculating XY in Model Builder using Calculate field

texts:
I'm trying to calculate the lat/long of a feature class in GCS NAD83.  The picture shows what goes below.

I can get it to work for distances and areas, but not a centroid.  What looks incorrect?


texts: Raster Subtraction with 2 different sized rasters

texts: I am a beginner to ArcGIS and am using it for my dissertation. I have two raster files of different areas and subtracted one from the other to find the difference between the two, however, it has given me a result which falls outside the boundaries of the smaller raster, why?


texts: Importing waypoints from OziExplorer wpt files

texts:
I do not succeed in importing waypoints from Ozi files (extension .wpt). I tried both converting through GPS-Babel to XML (.gpx) and importing in QGIS with the GPS-tool. Both times set to WGS 84, as is the setting in Ozi.

What should I do?


texts: Object-Oriented LULC Classification in Google Earth Engine

texts:
I am testing a object-oriented LULC classification approach in Google Earth Engine, using the code available at the following link: https://code.earthengine.google.com/4032fe102252b285ec249f395d0d6d0f.
The problem here is that the GEE reports this error message:

Layer error: Computation timed out.

I think that I exceeded the maximum number of pixels due to large AOI, because when I clip the AOI to a small area the code works out well.
Is there a way to fix /define the maximum number of input pixels or another way to avoid the Computation timed out?


texts:
I need to open a shapefile from ArcMap in R to use it for further geostatistical analysis. I've converted it into ASCII text file, but in R it is recognized as data.frame. Coordinates function doesn't work as soon as x and y are recognized as non-numeric.

Could you help to deal with it?


texts: QGIS Python - Convert the CRS of an in memory vector layer

texts: In a python script of qgis I load layers from different datasources using the QgsVectorLayer() and  the QgsRasterLayer() commands. The coordinate reference systems of the data sources differs. Since I have for some sources only the permission to read, I cannot transfrom the CRS of the datasource. Thus, i have to change the CRS from old CRS into new CRS in the in memory layer. What is the way to do it, e.g. to get an extent() in the new coordinates?


texts: Making selection color transparent when in editing mode in QGIS

texts:
in QGIS, how do I make the selection color transparent?
I need to split a polygon and i want to see the layer underneath it. But I have to select the polygon in order to cut it, and when i do that it turns yellow!


texts: TimeManager: "Have layers, but animation not possible"

texts:
I have used TimeManager before and really like however I have some problem now with QGIS 2.18.4 (can't pinpoint the problem to 2.18.4)

However, I use a PostgreSQL/PostGIS layer as input

I've tried to cast the time to time zone with/without time zone. QGIS recognises the type:

TimeManager manage to get the time as well:

But I got the error "Have layers, but animation not possible".
I've tried with the field "time_" instead, I've tried with the timeformat: %Y-%m-%d %H:%M:%S, %Y-%m-%dT%H:%M:%S, %Y-%m-%d %H:%M:%S+01

Any suggestion what I shall change?


texts: I have an eTrex 10. I'd like my GPS to remain in UTC (or GMT) time zone, but it does not offer this option as far as I can tell. I think "Western European" and "Ignore DST" is equivalent to GMT. Is this the best or only way to have data on this device recorded as close to UTC as possible?


texts: Spatial join is outputting seemingly incorrect result in ArcMap

texts:
I am attempting what I thought was a straightforward operation in ArcMap. I have a polygon layer where some polygons overlap. I also have a polyline layer. I would like to spatially join the polygon layer with the polyline layer using a "WITHIN" relationship.

The objective is to make sure every polyline which is either fully or partially within a polygon layer takes on the attribute from the (partially containing polygon). This is what the spatial join dialog looks like:

What I do see, however, is that many polylines do not obtain a containing polygon attribute and get a 0 instead. This occurs mostly in polygon border crossings and when a significant portion of the segment lies outside the polygon.

Which is very confusing because it implies that "WITHIN" and "COMPLETELY WITHIN" have the same behavior? I thought that "WITHIN" implied that even target features partially within the join features get matched...

What am I doing wrong here? Is there another way?


texts: How to use Boolean input in Arcgis with arpy? (Desktop edition)

texts:
So, my code goes something like this:

haveopenvel = arcpy.GetParameterAsText(10)
arcpy.AddMessage(haveopenvel)
if (haveopenvel == True):#error
    openvelftsec = Raster(openvelftsec)
    arcpy.MakeRasterLayer_management(openvelftsec,"openvelftsec","#",ext,"#")
else:
    ##Optional Model that can be

However, the loop always goes into the else loop even when I check mark the input in the tool. COuld anyone tell me what is going wrong here?


texts: Geocoding addresses to coordinates in QGIS

texts:
I've got an Excel file which includes known locations with their exact address (street, house number, city, state). This Excel file does not have a geographical column.
Also, I have a .shp file inside those boundaries of buildings.

I want to relate every address in the Excel file to a building from the .shp file in order to represent to the .xls file on a map and for further analysis.

Couldn't find a way to do it in QGIS. Any ideas?


texts:
I have QGIS 2.10.1 Pisa and GDAL 1.11.2, and I can't open MrSID (.sid) images en QGIS, I have lost three days trying to make it work, but I can't make it.

My Linux is 14.04 trusty tahr.


texts: Ajax request on WFS with BBOX doesn't change my number of features

texts:
I need to know how many features are displayed on my screen (depending on the BBOX). I read that you can do it with a WFS request and resulttype=hits but my number of features still stay the same, whatever I put in the BBOX.
Here's my code :

map.on('moveend', function(evt) {
            var extent = map.getView().calculateExtent(map.getSize());
            console.log(extent);
            var url = 'http://localhost:8080/geoserver/wfs?service=wfs&version=1.1.0&request=getfeature&typename=aisdata:ais_vessels& outputFormat=json&format_options=callback:getJson&resulttype=hits&bbox=' + extent.join(',');
            $.ajax({
                url: url,
                success: function(response) {
                    debugger;
                }
            });
        });

And on the following picture you can see that my BBOX is changing :


texts: Will core functionality in QGIS 3.0 exploit multi-threading?

texts:
As mentioned in a recent QGIS blog by Tim Sutton, the major changes to QGIS 3.0 are:

Updating Qt4 to Qt5
Updating PyQt4 to PyQt5
Updating Python 2.7 to Python 3
Improving the QGIS API

Also, a couple of years ago, Boundless seemed to have did some work on multithreading the execution of various processing algorithms back in 2015:

  Multithreading execution. Multithreading execution of Processing algorithms was implemented as part of the Google Summer of Code this year. We collaborated in mentoring that work, but still some work has to be done before the code is ready to be merged. We have already started reviewing the code and adapting it with the goal of having this important functionality ready before QGIS 3.0 is released,, and we hope to have it ready for testing soon.

Will core functionality in QGIS 3.0 finally exploit multi-threading?


texts:
I often work with MXDs that have a large number of group layers, and would like to find a specific layer by name (either TOC given name or an original file name) in the Table of Contents.

I tried Ctrl+F and search options, but that's not what I'm after.

Is there a shortcut or an option for this in ArcMap?


texts: Why does SDE.ST_NUMPOINTS() add extra vertices for lines with arcs?

texts:
Oracle 18c/10.7.1 EGDB:
I have a polyline FC with a feature that has an arc.

I want to count the number of vertices in the line using the SDE.ST_NUMPOINTS() function.
select
    sde.st_numpoints(shape) as num_vertices
from
    a_test_fc

NUM_VERTICES
------------
         302

The line only has 2 vertices. However, the result from SDE.ST_NUMPOINTS() was 302 vertices.
Question:
Why does SDE.ST_NUMPOINTS() add extra vertices to a line with arcs?
I'm guessing that it is densifying the arc. But, I don't understand why it would be doing that since it produces an incorrect result (defeats the purpose of the function).


texts: What are the units of the calculated area from a raster in decimal degree

texts:
I'm trying to calculate area for a raster file which contains 2 layers (0 for unsuitable habitat and 1 for suitable). It's been calculated successfully but I don't know what is the unit of the calculated values.

Calculated Area for 43873037 counts:365608.662141
               and for  3241800 counts: 27015.001513

Spatial reference: WGS84_world_Mercator,
Linear unit: Meter,
Angular unit: Degree,
Cell size:0.0083333338,
GRID format


texts: How to find depressions (sinks) and hills based on contour line features

texts:
I have contour lines (interval 0.1m), derived from a DEM (1m).

I used the thinning ratio to find only features, which are like circles.

I've imported these features into a Spatialite database.

The image below shows two groups of contour features (there are thousands of them). The left one shows the contour features of a small depression, the right one of a small hill.

My purpose is to find these feature groups, union them (or find the outer ring) and classify them to depressions and hills.

Sample of data:

ID;ELEV;thin_ratio
1;200;0.8

My first idea is to use the ST_Contains function to check if there is a feature with smaller or higher elevation then the feature which overlays and wrap it into a loop.

This query counts the contour features which are within the outer ring of a group. Excluding the outer ring: a.id!=b.id)

SELECT a.id, count(b.id) AS count FROM circle_test a, circle_test b
  WHERE ST_Contains(ST_BuildArea(a.geom),ST_BuildArea(b.geom))
    AND a.id!=b.id
      GROUP BY a.id

The other idea is to intersect them with a Depth in Sink raster (shown as blue in the image), processed with WhiteBox. The raster fills up sinks to a plane surface.
The problem is that there are big sinks with both depressions and hills within.

Solutions with Spatialite or Python are favored.


texts: OpenLayers how to maintain layer ordering while using SelectFeature control

texts:
I'd like my map to display a few layers in a particular order.  Some of these layers are Vectors, others are XYZs, and some of my Vector layers display attribution.

Initial layer ordering works fine until I attempt to provide attribution for some of my vector layers.  The problem comes when I use the SelectFeature control to handle vector attribution. Those layers that are added to this control are then rendered on top of the map's other layers (preventing the layers from being ordering in a meaningful manner).

I'm seeking a solution that allows the attribution of some vector layers while also allowing control over the layer rendering order.

I've filed an OpenLayers bug report: https://github.com/openlayers/openlayers/issues/661

Sample code can be found here:  https://gist.github.com/3537881

And an example of that same code can be viewed here:  http://labs.ecotrust.org/selectfeature/selectfeatureanomaly.html

This example demonstrates the re-ordering that occurs when the SelectFeature control is activated and deactivated (click on the Activate Control link under the map).

Again, I'm looking for a solution that allows me to maintain control over layer ordering while providing attribution for a subset of the rendered vector layers.

-Scott


texts: I'm looking for street/roads data in California with maximum speed limit attributed. TIGER data doesn't have the speed limit data, and open street map data is only attributed at about 25%. I'm looking for either open source or proprietary data. Is there any available?


texts: Is there any OSM based (or open source) online tool to draw circles on a map?

texts:
I'm looking for online tools to draw circles on a map (with accurate center and radius). By accurate edit, I mean:

edit center GPS coordinates
edit radius (in meters)

I tested https://www.mapdevelopers.com/draw-circle-tool.php. It does the trick, as I can draw circles with accurate center and radius, but I find the UI quite confusing.  Also, I'd prefer an open source project (so I can participate) based on OpenStreetMap.

I had the idea of using GeoJSON (on
http://geojson.io for example), but there's no circle feature yet in GeoJSON
- I tested https://geoman.io; it's OSM based but I can't edit circles center and radius

Is there any other tool to draw circles over a map (OSM preferred) and edit center and radius ?

Visual example (from https://www.mapdevelopers.com/draw-circle-tool.php):


texts:
I have hypsometric model of relief, and I need to visualize shadows. I was suggested to use GRASS command r.sunmask, so I tried to set time values and run command, but it gives me an error

"ERROR: Neither sun position nor east/north, date/time/timezone definition are complete"

So then I tried to put values in Position tab, which I founded by SOLPOS Calculator (http://www.nrel.gov/midc/solpos/solpos.html)

r.sunmask elevation=model@PLANINA output=model_sunmask year=2011 month=1 day=15 hour=9 timezone=1 east=45.909361 north=15.965538

But it didn't work out, too, it says "Calculating shadows from DEM..." but nothing's going on for hours of waiting. So I ask if that command works at all in GRASS, or I'm doing something wrong?

You can suggest some other application for same thing, except ArcGIS applications. Thank you.


texts: Can I use Dynamic Layer MapService on the client side using ArcGIS Online hosted weblayer?

texts:
I'm developing a native android application using ArcGIS Runtime SDK for Android. All layers are hosted in AGOL organizational account as feature layers and web layers.

Is is possible to use Dynamic Layer Map Service to render layers on client side ?


texts:
I am producing radar maps from the NEXRAD site near Omaha (OAX) using GrADS.  I'd like to make a basemap in QGIS.

I put a geotiff version of the GrADS output in QGIS and had it 'zoom to layer.'

In Composer I set the shape parameters to match GrADS 16:9 (in GrADS all lengths are 11", so 11x6.187).

The attached image superimposes my state and city overlay over the states and counties from the radar.  They don't match.

I'd like to make a bunch of basemaps in QGIS... if I can make them comport.  Can someone please help?


texts: Calculating speed along point feature class (.gpx) in ArcMap?

texts:
I have imported a .gpx file into ArcMap as a point feature class.
Each feature has a Point z as well as a DateTime.
I want to create a graph that shows speed as a function of distance. How do I do this?


texts: Convert a shapefile with variable values to a raster and distribute the values spatially in each grid cell

texts:
I want to convert a shapefile (polygons in WGS84) which contains 3 variables (float values) for each municipality of an area. In fact it contains 1035 municipalities. I want to convert it to a raster file  and distribute the values of municipalities to a specific grid with LCC projection .

After the suggestions described below in the answers/comments (to my initial question which I have edited) I reprojected the original shapefile with the municipalities to the final desired grid , I made the calculations with the area and then I rasterised the file as following:

ogr2ogr -t_srs "$outproj" emissions_municip_reproj.shp emissions_municip.shp

ogrinfo -dialect SQLite -sql "SELECT Esye_id, CO,NOx,SOx ST_Area(geometry) as AREA FROM emissions_municip_reproj" emissions_municip_reproj.shp

INFO: Open of $Database_outdir/emissions_municip_reproj.shp

# Add field, calculate and rasterize the shapefile
ogrinfo emissions_municip_reproj.shp -sql "ALTER TABLE emissions_gr_Cday_municip_reproj ADD COLUMN AREA float"

ogrinfo emissions_municip_reproj.shp -dialect SQLite -sql "UPDATE emissions_municip_reproj SET AREA =ST_Area(geometry)"
ogrinfo emissions_municip_reproj.shp -dialect SQLite -sql "UPDATE emissions_municip_reproj SET COarea = CO/ST_Area(geometry) * $RESX * $RESY"

gdal_rasterize -a COarea -ot Float32 -tr $RESX $RESY -te $xmin $ymin $xmax $ymax -a_nodata 0 -l emissions_municip_reproj emissions_municip_reproj.shp emissionsCO.tif

The result is almost correct. The problem is that it seems that values have not been distributed to all cells. In the picture below the municipality is shown with red line. It seems that over some part of the municipality there are no grid cells (no data).Any idea for this?


texts: How to uncompress a file geodatabase with open-source tools?

texts:
At 3.10.4 I tried to load a FGDB polygon feature class, but it failed, with an "Invalid Data Source" error.  I tried the load with both the Vector > Directory > ESRI FileGDB and Vector > Directory > OpenFileGDB methods:

It turns out that the FGDB had been previously compressed using the Arc Desktop Compress File Geodatabase Data tool.  Fortunately I have access to Arc Desktop and was able to uncompress the FGDB using the Uncompress File Geodatabase Data tool.

Following the uncompress, QGIS successfully loaded the feature class.

This experience leads me to believe that QGIS cannot load compressed FGDBs.  An internet search did not turn up anything.

Note that I'm not asking about FGDBs that have been zipped.

I am now curious; are there any open-source tools that will uncompress a FGDB?  Or allow QGIS to load a compressed FGDB feature class?


texts: Changing raster cells values inside of a polygon to the closest cells just outside the polygon!

texts:
I have been wrestling with this for a while and couldn't find out how I can do this. I want each cell of a raster inside of a polygon boundary to have the value of the closest cell just outside of the polygon.

I think the closest answer was the one at this question, but I couldn't apply it!

My plan was:

Use "Polygon to Raster" to change the polygon to a raster. The cells
inside the polygon become 1 and the points outside of the
polygon are nodata.
Use Con(IsNull(Raster),0,NULL) and assign it to Raster2
Add Raster1 and Raster2 to get nodata inside the polygon
It should be easy from here with nodata inside the polygon. We just use "Euclidean Allocation" and it does the job.

I'm stuck at step 2 and get an error:

  [SELECT * FROM VAT_Extent_RAS WHERE IsNull( "Value" )] Failed to
  execute (Con).

Do you have any thoughts? Is there an easier way to do this?

Update:

I'll try to summarize my progress here. I want to change all raster cells inside the polygon to "nodata" while the raster cells outside of the polygon remain the same. I used polygon to raster function first. This gives me a raster with 0 cell inside the polygon and "nodata" outside of the polygon. What do I do next?


texts:
I'm trying to import a .csv file to ArcMap which has more than one feature related to a single x, y coordinate. It doesn't seem to like it - can anyone help with the formatting?

e.g.

ID  name    x   y   veg name
1   swamp 1 208229  920604  green sprouts
2   swamp 1 208229  920604  red onions
3   swamp 1 208229  920604  cabbage
4   swamp 2 257589  935259  green sprouts
5   swamp 2 257589  935259  lentils


texts: Copying DEM attributes to vector point using QGIS?

texts:
I have a point that I drew on a tiff DEM file in QGIS. The point doesn't have any information in the attribute table except an ID of 1.

How can I make it copy the raster information based on the location it is on the DEM file?

For example, take/copy the height attribute from the location it is based on the raster/DEM file.

The closest solution seems to be the Point Sampling Tool plugin in QGIS. It extracted the values as needed.


texts: Replacing existing layer's data with new geojson in Carto?

texts:
I try to use CartoDB (https://carto.com/).
I created a layer with geojson file, but if the file is updated, how can I replace existing layer's data by new geojson?

I know, creating new layer and load new geojson, copying all styles and template settings from old layer to new and deleting old layer, is exactly same with replacing data.
But, "copying all styles and template settings from old layer to new" is very complicated task, so I want to replace just data only...

Do anyone know about how to do this?


texts:
I try to develop an app using QGIS Api. My code looks like this (it's an app with a gui) :

QgsApplication.setPrefixPath("C:\OSGeo4W64\apps\qgis", True)

qgs = QgsApplication([], True)

qgs.initQgis()

vlayer = QgsVectorLayer("C:\projects\files\Test_polygone.shp", "Test Polygones", "ogr")

if not vlayer.isValid():
    print("Layer failed to load!")

Each paths are good. No problem for there. I can't debug it because when i try to read vlayer.error().message(), it's an empty string.

Have someone an idea why it doesn't work ?


texts: How to use Collect Values with Cell Statistics in iterative model in ModelBuilder?

texts:
I'm trying to run multiple viewsheds in ArcGIS 10, reclassify these rasters and finally add them together. This is for a wind farm project where we need to repeat multiple viewsheds from various turbines and therefore the process lends itself particularly well to Modelbuilder.

My problem is that I am unable to get the Collect Values and Cell Statistics tools to work together to locate the reclassified rasters and add them together (see lower part of my model below):

There are a few threads dotted around that have been useful (here and here), but none have show their final models.

Can anyone show me how to correctly use the Collect Values tool to point to the reclassified rasters and use the Cell Statistics tool to sum these up?

Update with working model

Here is the revised (working) model for cumulative viewsheds. The user simple needs to point the model to a list of turbine point files for it to run. Thanks all.


texts: Neither "Photo2Shape", nor "Geotag and import photos" works

texts:
System:

Windows 7
QGIS 2.6 Brighton and QGIS 2.8 Wien
Photo2Shape 0.5.0
EXIFread 2.0.1
Geotag and import photos 1.1.0
EXIFtool 9.90

Problem:

EXIFtool and EXIFread are installed correctly. Both plugins are installed without any errors. But when I try to geotag some photos, nothing happens. I keeps working without progress. QGIS doesn't crash. I'm able to save and restart it.
I really tried and read all!
Anbody who can help?


texts:
When opening a data source in EPSG:25832 in a project with crs set to EPSG:31468 I am prompted to select a datum transformation (e.g. +nadgrids=BETA2007.gsb)

The nationwide standardized gridfile for the transformation of geotopography data BeTA2007 is already available and though implemented in QGIS.
By jan, 2019 the Freistaat Bayern / Germany switches coordinate reference system from EPSG:31468 to EPSG:25832 and for this purpose is going to provide an NTv2 gridfile "BY-KanU" („Bayern-Kataster nach UTM“). Various gridfiles are available in advance: https://www.ldbv.bayern.de/vermessung/utm_umstellung/trans_geofach.html from where I obtained a gridfile BY-KanU_Mittelfranken_Testdatei_vorlaeufig.gsb
How do I use these gridfiles in QGIS (applies for 2 and 3 as well)?
Under Settings | Options... | CRS it seems that I can define my own transformations, but I do not know how and where to select the desired gridfile.


texts: How can I fix the issue in the TIGER loader that causes it to fetch all files when loading Kansas?

texts:
I've discovered a bug in the TIGER loading script provided by Postgis 2.1.2.  When loading Kansas (state code 20) I get FTP statements like the following:

wget ftp://ftp2.census.gov/geo/tiger/TIGER2010/ZCTA5/2010/*_20* --no-parent --relative --recursive --level=2 --accept=zip --mirror --reject=html
wget ftp://ftp2.census.gov/geo/tiger/TIGER2013/FACES/*_20* --no-parent --relative --recursive --level=2 --accept=zip --mirror --reject=html
wget ftp://ftp2.census.gov/geo/tiger/TIGER2013/FEATNAMES/*_20* --no-parent --relative --recursive --level=2 --accept=zip --mirror --reject=html
wget ftp://ftp2.census.gov/geo/tiger/TIGER2013/EDGES/*_20* --no-parent --relative --recursive --level=2 --accept=zip --mirror --reject=html
wget ftp://ftp2.census.gov/geo/tiger/TIGER2013/ADDR/*_20* --no-parent --relative --recursive --level=2 --accept=zip --mirror --reject=html

in the loader script.

The problem with these statements is that all of the files in the FTP directories contain a '_2010' or '_2013', which means that these statements attempt to pull all of the files from the FTP source directory, not just the Kansas specific ones.

The problem can be fixed by replacing the '*_20*' with a 'tl_*_20*', but frankly the loader script generation functions are complex, and I'm not sure how to make this change without breaking anything.

The Postgis Github repo doesn't seem to accept issues which aren't pull requests, so I'm posting this question here in the hopes that one of the Postgis maintainers will see it and address the issue.

I can fix this manually for scripts I generate, but I imagine it would be more desirable to fix the root problem.  What's the best way to get this fixed?


texts: Problem in multiple dataframe visibility in arcgis 10

texts: I have added four dataframes on singe canvas. Within Arc environment everything works perfectly but when I export the map to .eps or .pdf or .jpeg two of them are working properly while the other two show a white background which destroys the look of the map. properties of the dataframes are same. kindly give ideas.


texts: Geoserver WMS character representation error (UTF-8)

texts:
I created a layer that is connected to an Oracle DB table. I can access the layer as WMS and as WFS. The character representation is correct in the WFS but has errors in the WMS service. I used Qgis and MapInfo for testing and I have the same problem with both.

Screenshot MapInfo info tool (WMS):

Screenshot MapInfo info tool (WFS):

I recognized that I have the same problem with all my WMS layers where german umlauts are used. In the Geoserver global settings the charset is set to UTF-8. How can I fix this?

EDIT: I tested the same WMS with ArcGIS Explorer and Gaia. Both can display the character correctly. The problem seems to be the software that is used to access the WMS... I will let you know when the ticket is answered.


texts: How to remove a layer from the Layerswitcher in Openlayers

texts:
I have an init-Function, there some layers are added with the option displayInLayerSwitcher = false; Thats ok, but in some situations i need to set this later in the application to true. So i want to reset the option. The ways i tried to figure out:

1:my_layer.displayInLayerSwitcher = true;
my_layers_switcher.redraw();
// Result is, that my_layer is resetted and in the LayerSwitcher nothing changed.

2: my_layer.displayInLayerSwitcher = true;
my_layers_switcher.destroy();
map.removeControll(my_layers_switcher);
my_layers_switcher = new OpenLayers.Control.LayerSwitcher();
map.addControll(my_layers_switcher);
// Same Result as by way 1

Anyone, who have an idea for this?


texts: I am using QGIS3 on a Mac and the Vector tool menu is missing entirely. The processing plugin is enabled and online it suggested deleting .qgis2, which I don't have. Can anyone help?


texts: How to query OSM server for supported zoomlevels?

texts: OpenStreetMap tile servers usually have 20 zoomlevels, but some have a different number. Is there a standard way to query the number of zoomlevels a tile server supports?


texts:
I am trying to delete field names where the field "FULL_NAME" = 1.

What am I doing wrong here?

I am not getting an error message, but nothing is being deleted from my file.  I am not sure why the delete field isn't working.

import arcpy
roads = "C:/hwy"
cursor = arcpy.da.UpdateCursor(roads, ["FULL_NAME"])
for row in cursor:
    row[0] == "1":
        cursor.deleteRow()
del row
del cursor


texts: How can I get the map to repeat to fill the screen?

texts:
I'm programming my first ArcGIS application. I've created a simple map. When I zoom out to a world view, the map doesn't fill the width of the screen. Ideally I'd like it to repeat until the available space is filled. How can I do this?

Here's the code I have so far:

var map = new esri.Map("map", {

basemap: "topo",
center: [0, 0],
zoom: 2
});

And here's the result:

How can I get the map to repeat to fill the white space on either side?

(Crossposted: http://forums.arcgis.com/threads/79757-How-can-I-get-the-map-to-repeat-to-fill-the-screen)

Update:

I've set the spatial reference to (I think) Web Mercator, and wrapAround180 to true, but this has no effect. New constructor:

var map = new esri.Map("map", {
            basemap: "topo",
            center: [0, 0],
            zoom: 2,
            sliderPosition: "bottom-left",
            spatialRefence: new esri.SpatialReference({ "wkid": 3857 }),
            wrapAround180: true
        });

What am I doing wrong?


texts: Accessing Data on Remote Server without Setting Data Sources everytime?

texts:
I'm in one of two offices within my company. The other office has a large library of GIS maps they've created on a drive they've named the M:\ drive. Of course, when I open the maps all of the links are broken since they're referencing data on their server, and not mine. As it stands I use the "Repair Data Source" feature and navigate to the folders, but it's often excruciatingly slow.

I'd like to open their maps without any broken links. I've mapped the drive through Windows, and used Catalog to create a folder connection, but obviously neither have worked.


texts: Leaflet: click feature to get external diagram

texts:
I am working on a Leaflet webmap which receives a layer withs polygons from ESRI REST API.

I want to implement this:
When the user clicks on a polygon a diagram shall pop up below the map.
The diagram shall use the attributes of the polygon (in my case measurements of different years).

ID   year1990   year1991   year1992   year1993    ...
1    -0.3       +2.2       -0.2       +0.5        ...
2    -4.3       -5         +6.8       +0.4        ...

I thought about transferring the polygon's attributes into a JSON and the visualize it with D3.js .

Does this makes sense to you? How would you solve this?


texts: Difference between MODIS sinusoidal and integerized sinusoidal grid

texts:
What are the differences between the MODIS sinusoidal and MODIS integerized sinusoidal grid? The sinusoidal grid seems to be currently used, whereas the integerized sinusoidal grid is marked obsolete (click "obsolete grids" at the very end of the page):

http://modis-land.gsfc.nasa.gov/MODLAND_grid.html

I'm also interested into why it is obsolete now, albeit it was invented later(?).

What advantages and disadvantages instoring and retreiving data do they have?


texts: Determining the Longest Segment of a Non-Continuous Polyline

texts:
I have a polyline that is non-continuous (please see picture, polyline is in red). I am looking to separate the polyline based on this geographic discontinuity so that I can calculate the length of each portion individually. Unfortunately I do not know how to separate the polyline in this manner, i have tried multiple things with no success. My best thought so far was to split the line at the vertices (each time the line changes angle) and then unsplit them, but that does not work for some reason. Thoughts? I am using ArcGis 10.5.1.

Note: the the dataset is large and contains quite a bit of variation. The Lines will often converge forming what looks to be a single line, though they each will independently follow the same path. These polylines are the result of a Least Cost Path analysis.

The Red numbers indicate the four different segments I am looking to individually calculate the length of. The Blue dotted line adjacent to segment 4 indicates a partial overlap from a different polyline that separates roughly half way through segment 4.

The overall goal of this analysis is to determine the longest polyline segment.


texts: OpenStreetMap to pgRouting missing attributes?

texts:
We are trying routing with pgrouting for our project but we noticed that NOT all the attributes we had in the row OSM data are in our database after we imported the data. we have checked in the original OSM file that the attribute tages are there but not in the database after using osm2pgrouting tool.  we followed the tutorial in http://workshop.pgrouting.org/chapters/osm2pgrouting.html . we needed the attributes tags like traffic_signals, lit.

have you ever experienced this problem?

could it be from osm2pgrouting conversion tool?


texts:
I've installed geopandas on ubuntu 17.10 and can't import geopandas. I can't find this particular error anywhere.

from pyproj import _proj
ImportError: cannot import name '_proj'

I'm using a virtual environment with Python 3.6.3 GDAL 2.2.1
I ran pip install geopandas which also installed dependencies.
Here is output of pip list:

bleach (2.1.2)
click (6.7)
click-plugins (1.0.3)
cligj (0.4.0)
cycler (0.10.0)
decorator (4.1.2)
descartes (1.1.0)
entrypoints (0.2.3)
Fiona (1.6.0)
geopandas (0.3.0)
html5lib (1.0.1)
ipykernel (4.7.0)
ipython (6.2.1)
ipython-genutils (0.2.0)
ipywidgets (7.1.0)
jedi (0.11.1)
Jinja2 (2.10)
jsonschema (2.6.0)
jupyter (1.0.0)
jupyter-client (5.2.1)
jupyter-console (5.2.0)
jupyter-core (4.4.0)
MarkupSafe (1.0)
matplotlib (2.1.1)
mistune (0.8.3)
munch (2.2.0)
nbconvert (5.3.1)
nbformat (4.4.0)
notebook (5.2.2)
numpy (1.14.0)
pandas (0.22.0)
pandocfilters (1.4.2)
parso (0.1.1)
pexpect (4.3.1)
pickleshare (0.7.4)
pip (9.0.1)
pkg-resources (0.0.0)
prompt-toolkit (1.0.15)
ptyprocess (0.5.2)
Pygments (2.2.0)
pyparsing (2.2.0)
pyproj (1.9.5.1)
python-dateutil (2.6.1)
pytz (2017.3)
pyzmq (16.0.3)
qtconsole (4.3.1)
setuptools (38.4.0)
Shapely (1.6.3)
simplegeneric (0.8.1)
six (1.11.0)
terminado (0.8.1)
testpath (0.3.1)
tornado (4.5.3)
traitlets (4.3.2)
wcwidth (0.1.7)
webencodings (0.5.1)
widgetsnbextension (3.1.0)


texts: Chinese encoding problem in ogr2ogr from DXF file to GeoJSON?

texts:
I have a dxf format file named test.dxf which has some Chinese characters inside. I tried to transform it to GeoJSON file using ogr2ogr -f GeoJSON test.geojson test.dxf and returned GeoJSON showing unrecognizable characters like this:

comp:甯岀背绉戯紙鍖椾含锛夊尰鑽鎶€鏈夐檺鍏徃"

I also tried with ogr2ogr -f GeoJSON test.geojson test.dxf ENCODING='UTF-8' but it doesn't work out.

How can I get the result correctly?

Related reference:

Converting shapefile from Shift_JIS to UTF-8 when the usual methods fail


texts: Clip raster by mask layer not working in QGIS 3.4.2?

texts:
I've been successfully using QGIS 2.18.21 to "clip raster by mask layer" and then run v.rast.stats. However recently v.rast.stats has stopped working in QGIS 2.18.1. The log does not show an error, it just causes QGIS to stop responding.

To work around this I downloaded QGIS 3.4.2. However, now when I try to run "clip raster by mask layer" it does not work and does not give an error in the log. The same thing happens when I try to run v.rast.stats. I've checked that all the layers are in the same CRS.

Has anyone else experienced this?


texts: GeoDjango: AttributeError: function 'GDALAllRegister' not found

texts:
Im trying to get started with GeoDjango. I am using this https://docs.djangoproject.com/en/1.8/ref/contrib/gis/tutorial/ tutorial.

When I run "python manage.py migrate"
This error occur: "AttributeError: function 'GDALAllRegister' not found"

I have installed GDAL through the OSGeo4W.

Im on a 32bit Windows7 machine with PostgreSQL 9.3 and postgis 2.1.

What should i do to solve this problem?

Unfortunately, I didn't find a resolution. But instead, I found a solution for a mac os x Yosemite machine.


texts: Fixing thousands of small overlapping polygons?

texts:
I'm using ArcMap 10.3.

I have a layer of different land use types. I found out that many polygons overlap each other and as a result the same polygon have two different land use types.

Using @radouxju's answer, I built a topology for the land use layer (with the mustn't overlap rule). Then, I used merge to fix the large overlapping polygons.

However, I have thousands of small overlapping polygons. I can fix all of them manually which will be time taking.

Is there a faster way to fix all these errors?

I have updated the map to show all land use types. Polygons are overlapping each others; not slivers.


texts: Adding Layer Files (.lyr, .lyrx) to QGIS

texts:
At college we use ArcGIS Desktop 10, but that does not function on my Mac so I have QGIS there instead.

I want to do my exercises with QGIS, but in ArcGIS Desktop we use .lyr files with ArcMap and .lyrx files with ArcGIS Pro, and they are not supported by QGIS.

Does anyone know, how I can connect the *.lyr files to my *.shp files to use them on my Mac with QGIS?


texts: Error in Built Up area Classification using Landsat 8 in google earth engine ( SVM, RF and Cart)

texts:
I am using Earth Engine to Classify an area using supervised classification with a combination of reflectance bands and indices. The built-up area gets mixed with fallow land (bare soil) and I have tried many different combinations and somehow I am unable to resolve it.
I have used built-up index, tillage index, ndvi, mndwi and evi if different combinations. I have referred to many journal papers and still unable to resolve it.
I have tried increasing and decreasing the training points too.
Could anyone please help me with this.
Here is the code for the classification. https://code.earthengine.google.com/20bc92c707702e1fdb32bd5325ff19af
This image shows the results where the cart and random forest classification mix with other classes and whereas SVM results fail to classify the built-up areas properly with same training points.
It would be great if the code is improved to get confusion matrix (kappa and overall accuracy).


texts: Export of reprojected layer from GEE fails

texts:
I am trying to export a resampled image from Google Earth Engine, and keep getting errors such as -

Error: Reprojection output too large (23734x8520 pixels).

As shown in the following screenshot -

The error can be reproduced using the sample script from https://developers.google.com/earth-engine/resample#reduce-resolution as follows -

// Load a MODIS EVI image.
var modis = ee.Image(ee.ImageCollection('MODIS/006/MOD13A1').first())
    .select('EVI');

// Display the EVI image near La Honda, California.
Map.setCenter(-122.3616, 37.5331, 12);
Map.addLayer(modis, {min: 2000, max: 5000}, 'MODIS EVI');

// Get information about the MODIS projection.
var modisProjection = modis.projection();
print('MODIS projection:', modisProjection);

// Load and display forest cover data at 30 meters resolution.
var forest = ee.Image('UMD/hansen/global_forest_change_2015')
    .select('treecover2000');
Map.addLayer(forest, {max: 80}, 'forest cover 30 m');

// Get the forest cover data at MODIS scale and projection.
var forestMean = forest
    // Force the next reprojection to aggregate instead of resampling.
    .reduceResolution({
      reducer: ee.Reducer.mean(),
      maxPixels: 1024
    })
    // Request the data at the scale and projection of the MODIS image.
    .reproject({
      crs: modisProjection
    });

// Display the aggregated, reprojected forest cover data.
Map.addLayer(forestMean, {max: 80}, 'forest cover at MODIS scale');

The script works fine and the resampled image forestMean is shown on the map.

However, trying to export the image with the following expressions leads to the above-mentioned error -

// Continental USA polygon
var pol = ee.Geometry.Polygon([[[-126.21093749999999,23.241346102386135],[-67.8515625,23.241346102386135],[-67.8515625,48.922499263758255],[-126.21093749999999,48.922499263758255],[-126.21093749999999,23.241346102386135]]]);

// Export
Export.image.toDrive({
  image: forestMean,
  description: "forest_example",
  scale: 926.6254330555,
  region: pol
});

It is important to note that the same error appears even when using a very small export region.

Update:

As Philipp Gärtner discovered, this seems to be an issue of export area size. Using a very small area indeed works.

For example, the following very small export area works -

While this one doesn't work -


texts: Is it possible to open *.dmp file and find out what cause ArcMap to crash?

texts:
Is it possible to find out what cause ArcMap to crash from *.dmp file? This files are located here:

  C:\Users\user\AppData\Local\ESRI\ErrorReports

I have made some Python add-in to ArcMap wich sometimes works and sometimes ArcMap crases. But I can't find out what part of my add-in causes this event. Is it possible to open *.dmp file and find something there? Or is there any other option how debug already compiled ESRI add-in?


texts: Problem in VersionEdit4.Reconcile4 method

texts:
I had an standalone application in which i am getting all the available versions in the database. Based on user filter he selects multiple versions and reconcile it.

 if (pVersionEdit4.Reconcile4(pVesrion.VersionInfo.Parent.VersionName, true, false, false, false) == true)

I am iterating user selected versions and generating conflict report.
But the previous one is Row level conflict check. But i want to check column level. For that my code changed to

if (pVersionEdit4.Reconcile4(pVesrion.VersionInfo.Parent.VersionName, true, false, false, true) == true)

But i am getting exception message Failed Security. Is any where i am missing.


texts:
In OpenJump there is a tool called ClosedVertexFinder which you can use to find vertices which are too close to each other.

https://github.com/ssinger/roadmatcher/blob/master/src/com/vividsolutions/jcs/qa/CloseVertexFinder.java

I would like to build a tool like this in QGIS.

Is there something that can already be used in order not to build the tool from scratch or will I have to calculate the distance between all vertices of a layer?


texts: Seeking QGIS equivalent to Path Distance of ArcMap

texts:
I am a university student and, in my University, the professors use and explain ArcGIS Desktop but I want to also learn QGIS.
I do "Path Distance" in ArcMap with the following parameters:

Input raster or feature source data: Network (a layer with 6 points)
Output distance raster: Network_dist.tif
Input surface raster (optional): Bythometry.tif
Max distance (optional): 9999999999999999999

How can I do this in QGIS?
The case is that:
We have 6 points and we want a raster with minimum distance from any of that points.


texts:
I have recently received my OS OpenData discs from Ordnance Survey, and for the last two evening tried to get this data into MapServer for ultimate use as a WMS for OpenLayers.

Safe to say I have struggled with how to get the data (.tif and .tfw) into a MapServer instance.

I have set Mapserver up correctly, in that I a directory structure such as:

MapServer
  data
  map
    mymapfile.map
  scripts
    gdal
    ...
  tmp

I am able to browse to my localhost (Windows) and retrieve a sample graphic from MapServer, so I am fairly confident I have MapServer setup correctly.

http://mapserver.localhost.com/scripts/mapserv.exe?map=C:\MapServer\map\mapfiletest.map&mode=map&FORMAT=image/png&WIDTH=256&HEIGHT=256

That is a sample of the Url I am using.

I have obviously Googled a fair bit over the last few evenings, but I dont seem to be finding any specific information to take me forward with getting my OS OpenData (streetview as an example) available through MapServer.

Does anyone have any links or specific information to help me out?


texts:
i have my Vector Data saved in a PostgreSQL PostGIS enabled Database.
The Data im displaying is very detailed, i want to simplify it through the simplify() method of Postgis.

I duplicated one of my tables (for test purposes to not mess up the real thing) and then used the simplify func"select ... simplify(...)..."

I got a result and inserted into the geom part of my new table.
After i created a new layer with this table as spatial information, whenever i click on Layer preview, im not greeted with a map, but with a download prompt for an xml which tells me i have an error.

Am i using the function wrong ?? or is my approach inside of geoserver wrong ??

Could i possibly save my current data as a shapefile, then simplify this and then let it be processed with pg2sql or something like that ??

Whenever i add the Layer to my Geoserver and press compute Bounding Boxes, geoserver comes up with 0, -1, 0, 1 .

Sorry if i have alot of questions but im really really stuck.


texts: I have a Shapefile that includes country divisions (Polygon). I have another Shapfile  that relates to the country's roadmap (Path). I also have another Shapfile  that relates to the village points (Points). How can I separate a part of the  Shapfile  for geographic divisions, which includes path and Points?


texts:
I've tried out the Hot Spot Map tutorial for ArcMap, which is working just fine, and I'm now starting to make some modifications so that it fits my data sets.
Instead of aggregating the points as in the tutorial I added a column for the weighting.

The result, using the Hot Spot Analysis (Getis-Ord Gi*) (Spatial Statistics) tool, can be seen in this picture:

What I'm trying to do now is to repoduce these results via MATLAB, just to get an idea of how the algorithm, which is described here, works. If I'm getting it right, the GI-score is calculated from the distance between one feature and all the others and their weighting.

Using these formulas I came to the following result:

As you can see, it looks like all the hot spots and cold spots are located in stripes.

Does ArcMap use a different algorithm from what is stated in the documentation or am I missing out on something else?

FYI, I used the default settings for the Hot Spot Analysis tool.


texts: How to set up a QGIS project with SQL Server which can be opened on another computer?

texts:
I am using QGIS with SQL server 2008.
I need to save project on one computer and then copy project to another 20 clients!

After I save project and try to open it on another computer I am getting error.

How I can setup project that can be used on another computer?


texts:
I have a set of GeoPackages created with QGIS. They have those layer_styles in it, which contain -- among other columns -- the styleQML column providing styling information.

I am reading and manipulating such GeoPackages with GeoPandas, and one of the actions I'm applying is renaming: renaming layers and columns in some layers -- to apply a naming convention of our project. Eventually, the new geopackage(s) will be saved to a new .gpkg.

My problem is with QGIS' styleQML content we found in layer_styles: QML is a XML document. I want to know what is the schema/model styleQML follows so that I can write a modifier (for whenever I change the name of layer/column, changes are properly applied to styleQML).


texts: How to get the file name from docked browser in QGIS 3x

texts:
Within QGIS 3x, there is a browser docked on the left hand side.

I would like to use this in a python plug-in, how can I use the highlighted file  and path within a Python Script? For instance to put into label, as in the above figure.


texts: Transforming coordinate in Swiss reference system LV95 into LatLong/WGS84?

texts:
I would like to transform coordinates expressed in the Swiss reference system LV95 into LatLong/WGS84. On the website of the Swiss  topography centre, there is an application to do so, however altitude data are required and I don't have them.

In Plotting species distribution on map in R, the user also wants to convert from Lv95 into LatLong. I was wondering if the solution he gets (reprojecting his shapefile into lat/long) could also work for me, as I don't have a shapefile but just separate coordinates.


texts: A accurate tool for digitization in qgis

texts:
I have some raster maps that  need to be digitized.
Me problem is that some of the shapes are very hard to be digitized.
Take for example the circles that need to be digitized in this map

Now matter how accurate I will be the result will never a perfect shape.

Is there some kind of plugin that provide better digitized tools for circles  or squares?


texts: Moving PostGIS database for Social Media to Cassandra?

texts:
I currently have a large and growing database of geo-social media. I'm using Postgres with postgis, and so far it's been quite good.

My use case thus far has been to query the database for all posts within specific areas, such as every tweet within each post/zip-code in the US/Canada (or other boundaries -- the posts are not geo-coded to the post/zip code. selection is through intersection of coordinates).

In the past this was relatively quick and easy to do, however, now that my data table has grown to over 350m posts, my queries take forever to complete. Each post code can take up to a day to retrieve the tweets from within its geometry.

I am currently debating what to do next. I know of a few solutions, but I'm curious what you think of them and if I'm missing anything.

throw more memory at postgres, and reduce disk searching. I'm currently working on doing this, but I wonder if it is a good long term solution.
find a better way to partition data into more tables. Currently, all of my posts are in the same table, with each column indexed. Given that I need to search through all of the posts to find the ones in the given areas, is this the wrong way to be doing it?
create a hybrid of Cassandra-PostGIS. Bulk store my data into Cassandra, and then pull out roughly the data I need (say, one state/province at a time) to a staging data table in PostGIS, then query individual post/zip codes using PostGIS. Is this overkill/horribly inefficient? I have looked into Postgres-XL, but it has been kiboshed by others in this project for cost issues.


texts: How to increase height of image header in Infowindow?

texts:
The CartoDB Infowindow default image header width is max 400, with height always set to a default, which looks to be about 250 so the image is really wide.

Is there a way to set the image to be more "square" let say 400x400 on the header or break the 400 max?


texts: Tilemill Styling Multiple feature Overlap some data

texts:
Please see the following Image.

I am using a single feature layer for road. Now I using the following code

    #roads {
  [ROAD_TYPE = 1] {line-width:8; line-color:#194536;}
  [ROAD_TYPE = 2] {line-width:8; line-color:#456321;}
  [ROAD_TYPE = 3] {line-width:8; line-color:#785632;}
  [ROAD_TYPE = 4] {line-width:8; line-color:#236485;}
  [ROAD_TYPE = 5] {line-width:8; line-color:#FFFFFF;}
  [ROAD_TYPE = 6] {line-width:8; line-color:#194536;}
  line-join: round;
  line-cap: round;
  line-smooth: 0.1;
  }

The problem is in the Black circles the white goes under the green road layer. But in the red circle the white layer overlap the green and other major layers. How can I solve it and make it a very nice style as Google or OSM map?


texts: Adding QTimer creates local variable error

texts:
In QGIS v2.18 I am writing a processing script combining two of the default scripts found in the Processing Toolbox, ‘Random Extract’ and ‘Differential Privacy’.  I believe I need to add a timer to allow enough time for Random Extract to create a temporary layer before Differential privacy begins working on it.  Through my studying I have come up with the following code:

QTimer.singleShot(1000, refreshLayers)

def refreshLayers():
    self.iface.mapCanvas().refresh()

However, I get the following error:
“local variable 'refreshLayers' referenced before assignment See log for more details”.

I take it the message is referring to the Log found in the GUI window... which is blank other than the error message.  Not sure if this helps but the variable name I am using is 'OUTPUT' and the actual temp layer name is called 'Extracted Points'.  Can someone lead me in the right direction?


texts: Why are ArcGIS Server SOC machines removed from host after Windows Update?

texts: Whenever Windows service pack updated on our servers, SOC machines were automatically removed from HOSTS. Is this a known behaviour? Are there any ESRI tech support articles on this? Machines are Win2008 R2 with ArcGIS 10.


texts: Seeking (Open Source) Mobile GIS app for survey?

texts:
I'm going to do an on-site survey for a business area development research.
Therefore I'll visit all of the biggest business district from a medium-scale city (100,000 population) to check present companies, economic sector, building type, etc. I prepared a survey sheet (like a questionnaire) that I would fill into Excel manually after to proceed with my work in QGIS.

Is there any open source mobile GIS application for smart phones that I could use instead?

I could already prepare or import the survey design in the app and do the survey with my phone and data is directly stored digital in a format to continue work (evaluation, analysis) on the computer after.

It would speed up the whole process a lot and could be used for further researches as well.


texts: Color ramp with legend based on attribute calculation in QGIS

texts:
Using the attribute table below, I want calculate a new field using the formula (field3+field4+field5+field6)/field7. Field 7 will be always larger than the addition of fields 3 through 6, so the new field will be a number between 0 and 1. Lastly, each feature should be symbolized differently using a color ramp.

For example, leinster will have a dark blue color with a value of 0.7 and connacht will have a light blue color with a value of 0.3

Below is what I want the map to look like. The map needs to include a legend as well.


texts:
My question is in the same vein as this one.

I'm attempting to convert precipitation lines (Precipitación Media Anual) to a polygon for Mexico in QGIS using the Polygonizer plugin.

After running the Polygonizer using the New Method, it returns the polygon but with a table full of NULL values after clicking Copy table field from line layer, the same as this problem.

The Old Method crashes at around 19% completion.

Here's a few related questions to this problem:

Is clipping the lines to smaller sizes the easiest option?
extracting each line class, polygonize, and then union

Though I'm unfamiliar with GRASS, is there a way to easily do this?

EDIT: I'm going to see if GRASS or PostGIS will do the trick

I'm open to any other ways to do this.

Thanks again for all the help!


texts: Figuring out coordinates of a line in a shapefile

texts:
I am working with a shapefile in R that is made up of 48 polygons. It is a map of Sweden with a line plotted to denote a biological boundary between north and south. I want to figure out the latitude of the line (see plot) but I cannot figure out how to extract that information from the shapefile.

I am wondering how to go about this. Can I plot each polygon separately to see if one/multiple polygons make up this line?


texts: Distance from Tract Centroid to all County Border Pairs

texts: I am starting with a shapefile of census tracts and a shapefile of counties. I am trying to calculate the shortest distance from the centroid of every census tract to every county bordering the county containing the census tract. For example, say County A borders 4 other counties (B, C, D, and E), and suppose census tract X is in County A. The output table will have four rows for census tract X. One row will have the shortest distance from the centroid of tract X to the border of county B. There will be three additional rows for X to C, X to D, and X to E. This will be the case for every census tract in the US. I will need generate the same file for census blocks as well. Does anyone have a python script that could be useful for this or advice for a GIS newcomer with this complex task?


texts: Automatically delete overlapping polygons

texts:
I have 200 000 polygons that needs to remain like circle and they can't be on top of each other. I need to find a way to correct the 8000 errors I've got automatically.

A picture's worth a thousand words :


texts: Pop up appearing blank after JSON parse?

texts:
My aim is to retrieve the attribute list of a WMS layer. The concept is to fetch it as JSON data, then parse it and display it in a pop up. I am using OpenLayers3 with local GeoServer and JavaScript.

When I enter the WMS request URL in browser, I get response as given below.

{"type":"FeatureCollection","features":[{"type":"Feature","id":"OSM_Grid_25000_include_hazard_line.60","geometry":{"type":"MultiPolygon","coordinates":[[[[623207.0023,1133268.8385],[609515.9426,1133223.6554],[609472.75,1147046.3421],[623158.4075,1147092.0522],[623207.0023,1133268.8385]]]]},"geometry_name":"the_geom","properties":{"OBJECTID_1":162,"OBJECTID_2":162,"OBJECTID":25925,"OSM_25K_IN":"C 43 K 3/SW","Map_No":39,"Map_KL_No":"KL 39"}}],"totalFeatures":"unknown","numberReturned":1,"timeStamp":"2019-09-19T12:06:17.336Z","crs":{"type":"name","properties":{"name":"urn:ogc:def:crs:EPSG::32643"}}}

Here, I need the value of "Map_No" to be displayed inside the pop up. Hence I wrote the code given below.

map.on('singleclick', function(evt) {
            var view = map.getView();
            var viewResolution = /** @type {number} */ (map.getView().getResolution());
            var url = tiledOSM.getSource().getGetFeatureInfoUrl(evt.coordinate, view.getResolution(), view.getProjection(), {
                'INFO_FORMAT': 'application/json',
                'PIXELRADIUS': 5,
                'FEATURE_COUNT': 10
            });
if (url) {
        var parser = new ol.format.GeoJSON();
        $.ajax({
          url: url,
          dataType: 'jsonp',
          jsonpCallback: 'parseResponse'
        }).then(function(response) {
          var result = parser.readFeatures(response);
          if (result.length) {
            var info = [];
            for (var i = 0, ii = result.length; i < ii; ++i) {
              info.push(result[i].get('Map_No'));
            }
            content.innerHTML = info.join(",");

          }
        });
      }

      overlay.setPosition(evt.coordinate);
    });

When I run this code, the pop up is appearing blank. Since I am a beginner, I am not able to understand it.

Should I add any JavaScript reference link to make this code work?


texts: How to highlight a road segment in my map?

texts: I have a table in QGIS that has 2 columns that I use to create road segments, StreetID and SegmntID. This creates a unique ID for each road segment. For example, for a road named Maine St, which I associated a 1 to identify, there might be 3 segments. So each segment will have a unique ID of 1-1, 1-2 and 1-3. I am going to print one sheet for each road segment. In Composer Manager, I have each sheet set up where you see that specific street segment in the middle of the map but you also see the ends of other segments attached to it in the same map view. I would like to be able to highlight that specific segment for that sheet so that it is obvious that that is the segment for that sheet. How would I set that up in the Composer Manager?


texts: How to adjust Python code from PyQGIS 2 to 3 to create a buffer on a shapefile

texts:
I am starting to learn PyQGIS and following tutorials online. However, I have QGIS 3.0 and most of the tutorials use QGIS version 2. I am trying to do a buffer on a point shapefile, but only have the code for Qgis 2. I have looked at the backwards incompatible changes (https://qgis.org/api/api_break.html) but cannot work out how to adjust the code. The 'GeometryAnalyzer' class, which contained the buffer function in Qgis 2 no longer exists, so I am not sure where to find it now. I am sure there is a simple solution..

The code for Qgis 2 as far as I know is as follows:

from qgis.analysis import QgsGeometryAnalzer
layer = iface.activeLayer()
QgsGeometryAnalyzer().buffer(layer, 'OutputFilePath/FileName.shp', 0.5, False, False, -1)

How can this be adjusted to Qgis 3?


texts:
How to labeling this rows and column of points so we could identified each point belong to what rows or what column depending on the polygon or starting point without depending by X or Y coordinate because the variations of points sometimes doesn't straight but askew by some degree?

I'am using arcmap btw, so arcpy or model builder or just simply labeling proccess would be big help.

sample created manually.


texts: Why does onDemand setting remove my buttons from my menu?

texts:
When I set the onDemand="false" attribute, my buttons disappear from my menu!

I'm setting the attribute to false so that I control its behaviour using the OnUpdate override in code. Am I doing anything wrong?

<ArcMap>  <Commands>    <Button id="SayHi_Btn" onDemand="false" class="SayHiButton" caption="Say Hi!" image="Images\SayHi_Btn.png" />    <Button id="SayHello_Btn" onDemand="false" class="SayHelloButton" caption="Say Hello!" image="Images\SayHello_Btn.png" />    <Button id="SayAloha_Btn" onDemand="false" class="SayAlohaButton" caption="Say Aloha!" image="Images\SayAloha_Btn.png" />  </Commands>  <Menus>    <Menu id="ChooseGreeting_RootMenu" caption="Choose Greeting" isRootMenu="true">      <Items>        <Button refID="SayHello_Btn" separator="true" />        <Button refID="SayAloha_Btn"  />        <Button refID="SayHi_Btn" />      </Items>    </Menu>  </Menus></ArcMap>


texts: MSSQL - Cast Lat/Long As A Native Geometry Type Into a New Column

texts:
I have two separate Columns, one for Lat and one for Long. What I need to do is to cast both into a third column that I have called Shape and is of type Geometry.

I was able to find the below as to how to cast into geometry but am unaware if this is the best way. Is also like suggestions as to how to drop it into the third column.

SELECT  *,
          geography::STGeomFromText('POINT(' +
                CAST([longitude] AS VARCHAR(20)) + ' ' +
                CAST([latitude] AS VARCHAR(20)) + ')', 4326) as GEOM,

          geography::Point([latitude], [longitude], 4326) as SAME_GEOM

  FROM **table**
  WHERE ('latitude' <> 0) AND ('longitude' <> 0) AND ('latitude' IS NOT NULL) AND ('longitude' IS NOT NULL)


texts: How do you sort a Standalone table that's already visible within ArcMap?

texts:
I'm using ITableWindow2.FindViaStandaloneTable--that works.  I can get the ITableWindow for the stand alone table and the table from that.

Then I try to sort using an ITableSort.  No exceptions, but the table doesn't sort.  Is ITableSort the right approach or is there a different approach to sorting within a Table Window?

ITableSort sort= new TableSortClass();
sort.Fields = "MDATE";
sort.Ascending["MDATE"] = true;
sort.Table = standAloneTable.Table;
sort.Sort(null);

The goal is to programmatically sort a field. The way an end-user would manually right click a column on an attribute table and sort ascending or descending.


texts: Problems with layer-selector code on joomla article

texts:
I have three visualizations I have created,you can see them separately on the following links
http://www.locattivitaecon.altervista.org/index.php/mineral-deposits-map
http://www.locattivitaecon.altervista.org/index.php/roads-port-airport-population-map
http://www.locattivitaecon.altervista.org/index.php/countries-regions-provinces-map
Now I want to create only one map with all the layers of the previous maps and to be able to select the maps separately. I have written the following code searching for different examples but I don't know what to write on the part where I put *. Any ideas?

All

Places

Cables


texts: Serve a subset of OpenStreetMap as a WMS from GeoServer?

texts:
I am trying to develop an application which should serve a small subset of OSM as a WMS for only my study area. I need it to be able to load higher-resolution tiles as the user zooms in.

I'm feeling a little lost right now. I can publish a WMS of a GeoTIFF in GeoServer. My first attempt was to save a QGIS project as a GeoTIFF, but obviously this is of course just a static image, so the resolution degrades when you zoom into it.

I download the "QTiles" plugin, which does a good job of generating a directory structure full of tiles, and so I am wondering, can I can simply somehow load this directory structure as a data source to GeoServer, and have it handle serving these tiles automatically?

My last-resort is to write a python script which will sew together the highest-resolution tiles into one contiguous image and just use this as the WMS data source, but I feel like this will probably not run optimally.

GeoServer version is 2.7.1.1


texts:
I'm just beginning with GeoDjango and just came across with term SRID.
What I've understood so far is that I'd need to know a little something about the map projection to for example calculate distancese between two locations. Threrefore:
What is the correct Nokia HERE maps SRID? And where to find this information?


texts: Smoothing grid of cells inSrcGIS for Desktop?

texts:
I have a grid of cells that outline sharp transitions between different neighbour values:

I would like to "smooth" the transition between the cells so that I obtain some sort of gradients between the different values.

How can I do it with ArcGIS?


texts: Using findAndReplaceWorkspacePaths to change datasource for ArcSDE geodatabase?

texts:
I manage to change my datasouce by right clicking on an mxd file and choosing a different sde database connection file as shown in the following image.

However When I try to do this using arcpy script the datasource it is not changed, but the new mxd is created (having the old datasource)

import arcpy
mxd = arcpy.mapping.MapDocument(r"C:\Users\SMALIS\Documents\ArcGIS\pythonSwitchSource.mxd")
mxd.findAndReplaceWorkspacePaths(r"Database Connections\sdeConnection1.sde",r"Database Connections\sdeConnection2.sde")
mxd.saveACopy(r"C:\Users\SMALIS\Documents\ArcGIS\pythonSwitchSourceNEW.mxd")
del mxd

Can somebody explain me why is that happening?

I use ArcGIS 10.2 for Desktop.
Oracle 11g


texts:
is there a tool in QGIS which is able to perform a spatial join with the possibility of define a snapping radius? For better understanding:

I have 2 point shapes with attributes I need to join by location. The problem is that I want to define the radius (max distance between related points). Is there a solution for my problem or a workaround?


texts:
I want to display a heatmap (png) on a globe in an android app.
The heatmap is generated programmatically. The heatmap consists of colored circles.
I upload the png to geosserver and display it as wms layer in my app.
I use EPSG 4326 as projection to display it on the globe.
Unfortunatly, the heatmap is distorted at the poles. It's shown in the picture.
How can I fix that?


texts: I am trying to calculate the optimal route from one archaeological site to multiple sites. I already have the r.cost raster. My question is if I can use r.drain to multiple destination points.


texts: Drop M values from MULTIPOLYGON with PostGIS?

texts: The title is self-explanatory. I wonder if there is a simple function like ST_Force2D for the M values?


texts: Scratch workspace for ArcGIS Desktop: file name conflict for multiple users

texts: I have installed ArcGIS Desktop on a server (it is not ArcGIS Server installation) where I want to give access to users to use toolbar (developed using Python-Addins) that I have developed. All intermediate files are written to a folder "C:/Temp" (I am not specifically writing those but only have specified scratch folder to be "C:/Temp"). If multiple users access it simultaneously, will there be a file naming conflict inside "C:/Temp"? If yes, is there any solution to it (how do I make file names unique even if multiple users are doing the same process or running same function at a moment)?


texts: Importing TIF raster in QGIS, default color band

texts:
Where can I change the default color bands for importing raster images?

When I import a TIF image the green color band is not set causing a purplish image.

When I set it in the properties it views correctly, but I have a lot of images to import and would like it set by default.

A week ago I was not having this problem, but I can't find how to set it back.


texts: Interpreting number on x-axis of semivariogram in ArcGIS Desktop?

texts:
Using ArcGIS 10.3.1

The attached image shows a semivariogram. The lag size is 65000 (meters), which makes sense based on my data and map units. My understanding of the x-axis is that it is the distance between the pair points. That is, points on the left are closer together, points on the right are farther apart. Why does my x-axis show such small numbers? 6.5 x 10 -5?  What are the units?

I realized by modifying it more that the end of the x-axis is the lag size * # of lags, or in this case, 65000*10, or 650,000. But again, the units should be 6.5 x 10^5, not 10^-5.


texts:
I am getting an Error while trying to execute "table to a domain".

The domain is "Nednaam":

The table is a .xls file, named "NedFlora", also a tekst file and looks like this:

My input:

My output:

What did go wrong? There are no double value's in the table.


texts: How is a database related to a .mxd file?

texts: What exactly is an mxd file, and how is it related to a database? Can a single database have multiple mxd files?


texts: Find duplicate vertices in SQL Server geometry line (ArcSDE)

texts:
I have a line in a ZM polyline featureclass that has a invalid geometry. My suspicion is that the line doubles back on itself somewhere, which I've found SQL Server doesn't like. Anyone know of a quick SQL method or query that could help me identify the suspect bad points that are jacking up my geometry? The string representation looks like this:

1835815.86 12887142.42 0 0, 1835816.72 12887142.68 170 170, 1835817.53 12887142.76 349.99 350, 1835817.52 12887142.76 559.99 560, 1835817.78 12887142.76 659.99 660, ....

Also, I'm wondering if I could use a regular expression and a look ahead and/or look behind to find duplicate numbers??


texts: GeoTools class VectorToRasterFactory is no more

texts:
I'm very new to GIS and GeoTools too. I'm upgrading org.geotools dependencies of a Maven project from 2.7.0 to 15.2. I got a compilation error about VectorToRasterFactory, It looks like there is no VectorToRasterFactory in the newest libraries. I didn't found any docs about how to substitute this class or how to get the same features.

UPDATE
I upgrade the libs but I got an issue during runtime. Here is the code:

...
GridCoverage2D cov = null;
try {
    org.geotools.process.Process p = Processors.createProcess(new NameImpl("vec", "VectorToRaster"));
    Map<String, Object> map = new HashMap<String, Object>();
    map.put(AbstractFeatureCollectionProcessFactory.FEATURES.key, srcCollection);
    map.put("attribute", "value");
    map.put("rasterWidth", dimension.width);
    map.put("rasterHeight", dimension.height);
    map.put("bounds", bounds);
    Map<String, Object> result = p.execute(map, monitor);
    cov = (GridCoverage2D) result.get("result");
    return (PlanarImage) cov.getRenderedImage();
} finally {
    if (cov != null)
    cov.dispose(true);
}
...

I can't get a Process instance, Processors.createProces... always returns null.


texts:
I am trying to fill sinks but only to a certain level.  For example - if I want to fill them by 5 meters - this means that a sink that is 15 meters deep - will fill only to 10 meters and a sink that is 3 meters deep would fill completely.

Using ArcGIS 10.0 - service pack 4.

Any help of this would be appreciated.

I thought I created a script to do what I wanted, but my client clarified something different than what I though.  Attached is an image of what I have and what I want for a custom fill tool.


texts: Create a parametrized SQL view based layer using geo.Geoserver Python API

texts:
I created "by hand" (using GeoServer HTML configuration interface) a layer based on the following SQL view:
SELECT p1.suid, p1.id, p1.geometry, p1.classname, p2.date, COALESCE(p2.percentage, 0) - COALESCE(p1.percentage, 0) AS percentage
FROM get_300km_percentages('%classname%', '%enddate%'::date, '%prevdate%'::date) p1
LEFT OUTER JOIN get_300km_percentages('%classname%', '%startdate%'::date, '%enddate%'::date) p2
ON p1.suid = p2.suid
ORDER BY
    percentage %order%
LIMIT
    %limit%

In the picture, we can see that all parameters, including %order% and %limit%" are ok. They were correctly guessed by the API.

How can I create a layer like that using geo.Geoserver or another Python-based GeoServer config library?


texts: How to update the Layer Tree View after changing the color ramp in PyQGIS?

texts:
After changing a raster style/color ramp (wich is working properly to define the class colors and values) using:

layer = processing.getObjectFromName('NDVIrec')
provider = layer.dataProvider()
myColorRamp = QgsColorRampShader()
myColorRamp.setColorRampType(QgsColorRampShader.INTERPOLATED)

CumulativeLower = 0.02
CumulativeUpper = 0.98
min, max = layer.dataProvider().cumulativeCut(1, CumulativeLower,
CumulativeUpper)
range = max - min

add = range / 11
int2 = min + add
int3 = min + 2 * add
...
int10 = min + 9 * add

colDic = {'1': '#a50026', ..., '11': '#006837'}

valueList = [min, int2, ... int10, max]

lst = [QgsColorRampShader.ColorRampItem(valueList[0], QColor(colDic['1'])),
       QgsColorRampShader.ColorRampItem(valueList[1], QColor(colDic['2'])),
       ...
       QgsColorRampShader.ColorRampItem(valueList[10], QColor(colDic['11']))]

myColorRamp.setColorRampItemList(lst)
myRasterShader = QgsRasterShader()
myRasterShader.setRasterShaderFunction(myColorRamp)
myPseudoRenderer = QgsSingleBandPseudoColorRenderer(provider,
                                                    1,
                                                    myRasterShader)

layer.setRenderer(myPseudoRenderer)
layer.triggerRepaint()
iface.layerTreeView().refreshLayerSymbology(layer.id())

I get the expected result in the map canvas, but the Layer Tree View does not update properly, as seen here (the class values wont show up):

The only way i found for showing the values is to open the layer properties editor and hit the "apply" or "ok" button, wich is not sufficient for my application.

Am i doing something wrong with the code for updating the Layer Tree View (bellow)?

iface.layerTreeView().refreshLayerSymbology(layer.id())