texts:
I have to regularly create KML polygons of map area outlines (white outline, QGis width 1.754, Google Earth width 5). I tend to have map image source files. I know how to georeference a source image to a raster, and would like to create area outlines from raster sections. Example:

This is the existing map of SW Pacific swordfish stocks. I'm looking to update this.

This is the management area image which will be the source file. There are two possible routes for me to create maps from this:

Just do the subarea boxes (1N, 1C etc).

1.1 Do this by hand in google earth. Time consuming but gets it done.

1.2 Make a csv of WKT points, import into QGis with QuickWKT. Firstly: lots of work to manually enter it all. Secondly: since subarea 2 boxes cross over the +/- 180 degrees line, this is handled wrong by QGis, splaying all the way around the world rather than across the 180 line.

1.3 georeference the image then trace the subarea borders with the polygon creation tool, save as individual polygons, etc. But might as well do in Google Earth in that case.

Do the high-seas areas only. This would be ideal. These are the bits outside of the 200 nautical mile buffer zones around land masses. Does anyone know if it's possible to create a polygon outline from an area of georeferenced raster, e.g. if I first colour the map image in Gimp/Paint/etc, and flood fill the high seas areas in 6 different colours (1 per subarea), then georeference, could I create polygon by colour, convert polygon to outline, save outline as KML? Coloured areas example below, ignore the green dots.

Will this all be confounded by the +/- 180 degree line problem?

I've found this plugin which claims to convert coloured areas to polygons.... but it doesn't work well.


texts: Openlayers doesn't want to show map throught Mapserver

texts:
I have a configured Mapserver with WMS service, I know Its working because I tried in QGIS and I can see the map. The problem is that openlayers doesnt want to show the map via mapserver service. Yes, the projection is same...

Index.html

<!doctype html>

  <meta content="text/html; charset=utf-8" http-equiv="Content-Type">
  <title>MapServer WMS OpenLayersben</title>

MapServer WMS OpenLayersben

mapserv file:

MAP
  NAME "MAPSERVER QUICKSTART"
  EXTENT 19.64 46.12 20.77 46.80
  UNITS dd
  SHAPEPATH "/home/"
  SIZE 1024 768

  IMAGETYPE PNG24

  PROJECTION
    "init=epsg:4326"
  END

  WEB
    METADATA
      ows_enable_request "*"
    END
  END

  LAYER
    NAME "Csongradutak"
    STATUS ON
    TYPE LINE
    DATA "csong_utak"
    CLASS
      STYLE
        COLOR 0 0 0
        OUTLINECOLOR 0 0 0
      END
    END
  END


texts: Adding legend for World Ocean Base map in ArcGIS Desktop?

texts:
I am using World Ocean Base as a background for a map I created in ArcGIS Desktop 10.6. I have been asked by a reviewer to include the meaning of the shading in the legend.

I know the shading on the base map relates to bathymetry, but is there a way to add this information to the legend as it does not automatically appear as with other layers?


texts: Integrating sea level curve data in a bathymetric grid

texts:
I am attempting to apply these data recording where sea-level was in the Gulf of Mexico (in meters below mean sea level) over the past 20,000 calendar years at a 500-year resolution in ArcMap:

What would be the best way to integrate this into a bathymetric raster so that the resulting map could be visualized with the various "shorelines" over time? Ideally, a .GRD file could be exported with an added column to the existing x, y, and z data, which may be tricky for visualization because it's the vertical shoreline progression is not always in a constant direction with time.


texts: How to change a Esri Scale Bar's font color

texts:
I'm working on my first ESRI web map project.  In the Default.Aspx page the scale bar is initiated as shown below :

dojo.require("esri.dijit.Scalebar");
map.on('load', function () {
     //resize the map when the browser resizes
     var scalebar = new esri.dijit.Scalebar({ map: map, scalebarUnit: 'english' });
}

Which produces the scale bar as shown below :

How do I change the default grey color for the bar to red?

Console Output :

I'm using ArcGIS API for JavaScript version 3.x.


texts: Does ArcGIS Pro utilize response paging in their new wfs feature, and, if not, are there work-arounds?

texts:
ArcGIS Pro now has the ability to add a feature layer by WFS - however, it appears the functionality is limited in some respects based on the size of data requested.  Does anyone have knowledge of whether the WFS functionality in Pro uses WFS response paging, and, if not, if there are other known work-arounds for loading many features from an OGC WFS layer?

In ArcGIS Pro, there are two 'modes' for loading data, On Demand, and Snapshot.  On Demand mode loads only what is in the boundary, so you can scan around the map and see different features, and is not related to this question (unless someone knows how to easily build a layer off of such panning around).

Snapshot mode will load whatever you have set as the maximum features to return.  As stated in the documentation:

  "To avoid performance issues, the default value is set to 3000 features."

What I have found so far is you can load tens of thousands of features, or fewer than 3k, of course depending on the layer, specifically how much data is sent for the features (eg how many attributes).  ArcGIS 'fails silently' when some kind of limit is hit, and I am not able to determine what defines that limit:  it just won't display anything, and won't display an error message.  If you then lower the limit, it will display features again after a pause, if it will in fact display them.

Best for me would be to get response paging working, and that starts by knowing if Pro can even utilize response paging, and if so what arc-specific settings if any must be used there are to get it working.  Failing that, does anyone have work-arounds, or know of settings that can be changed to allow larger loads, or some way to see more information about the reason when a WFS fails to load or doesn't display data?


texts: Why are there changes in location when printing inserted points?

texts:
Añadí una capa de texto delimitado por comas y luego agregue una imagen satelital de google earth. Realizo un diseño de impresión y al imprimir cambia la ubicación de los puntos en varios metros. Cual es la solucion?

(Edit) Google Translation:

I added a layer of comma delimited text then added a satellite image from Google Earth. I do a print design and during printing the locations of the points change by several meters.  What is the solution?


texts: Multiple layers from one WMTS-capabitilies source in Openlayers

texts:
I want to use multiple layers from one WMTS-source. I started from the WMTS Layer from capabilities example. I would like to collect the WMTS-layers in a layer-group baselayers, that I'm adding later to the map together with other overlay-layers.
What I've tryed so far is the following code, but I always get an error (TypeError: a is undefined), pointing to the code: var map = new ol.Map({ ...
I cannot find the error in my script, but maybe there is a cleaner way to code, what I intend to do?

Html-file

<!DOCTYPE html>

    <meta charset="utf-8">
    <title>Multiple WMTS Layers from Capabilities</title>
    <link rel="stylesheet" href="https://openlayers.org/en/v4.3.1/css/ol.css" type="text/css">

    <link href="./css/layerswitchercontrol.css" rel="stylesheet">

The javascript code:

    function initialize() {
      var parser = new ol.format.WMTSCapabilities();
      var baselayers;

    var toplayers = new ol.layer.Group({
        title: 'Regions',
        layers: []
    });

      fetch('http://wmts1.geoportail.lu/opendata/wmts/1.0.0/WMTSCapabilities.xml').then(function(response) {
        return response.text();
      }).then(function(text) {
        var result = parser.read(text);

        var hybrid_s = ol.source.WMTS.optionsFromCapabilities(result, {
          layer: 'hybrid',
          matrixSet: 'EPSG:3857'
        });
        var ortho_2016_s = ol.source.WMTS.optionsFromCapabilities(result, {
          layer: 'ortho_2016',
          matrixSet: 'EPSG:3857'
        });

    baselayers = new ol.layer.Group({
        title: 'Background-Layers',
        layers: [
            new ol.layer.Tile({
              source: new ol.source.OSM(),
              opacity: 1.0,
              title: 'OpenStreetMap',
              baseLayer: true,
              visible: false,
            }),
            new ol.layer.Tile({
              opacity: 1.0,
              title: 'Hybrid',
              baseLayer: true,
              visible: true,
              source: new ol.source.WMTS(/** @type {!olx.source.WMTSOptions} */ (hybrid_s))
            }),
            new ol.layer.Tile({
              opacity: 1.0,
              title: 'Ortho 2016',
              baseLayer: true,
              visible: false,
              source: new ol.source.WMTS(/** @type {!olx.source.WMTSOptions} */ (ortho_2016_s))
            })
          ]
        });
      });

    var map = new ol.Map({
        target: 'map',
        view: new ol.View({
            center: [685000,6380000],
            zoom: 11
        }),
        controls: ol.control.defaults().extend([
            new ol.control.LayerSwitcher() // http://viglino.github.io/ol3-ext/examples/map.switcher.html
        ]),
        layers: [baselayers, toplayers],
        interactions: ol.interaction.defaults().extend([
            new ol.interaction.Select({
                layers: [baselayers, toplayers]
            })
        ])
    });
};


texts:
I would want to use arcpy in my python code. I am able to import it through the ARCGIS desktop python console. But I am not able to import it into Python IDLE. I get the following error

Traceback (most recent call last):

File "C:\Python26\lib\site-packages\win32com\client\dynamic.py", line 246, in
__getitem__raise TypeError("This object does not support enumeration")

I mentioned the following folders as a part of my python path.

C:\Program Files\ArcGIS\Desktop10.0\arcpy;
C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy;
C:\Python26\ArcGIS10.0\Tools\Scripts

I found the same issue being mentioned.

import arcpy yields "TypeError: This object does not support enumeration"

The user says it is due to the following:

  I wrote for Arcgis 9 and 10 in the same directory. Something about the initialization logic in import arcpy finds the custom arcgisscripting.py used for running python 2.6 with arcgis 9.3 when they are located together.

I have only ARCGIS 10 and i am not aware of how arcpy looks up for its modules.


texts: How do I control the behavior of the QMessageBox

texts:
I'm working on my first plug-in.  In an attempt to validate the user selections on the dialog, I am using a message box to return messages to the user if an error is found.  The messgae box displays properly, however when the user selects the ok, the dialog closes along with the message box.  How do I control the behavior of the message box so only it will close and not the dialog too?

    def validate(self):

    message = ''

    if len(self.myListWidget.selectedItems()) == 0:
        message = 'Choose at least 1 Soil Survey Area'

        QMessageBox.critical(self, 'Execution Error', message)


texts: Is ArcGIS API for Python or ArcREST faster than arcpy.da for editing hosted features?

texts:
I have a large hosted feature layer on ArcGIS Online of 52,000 points. I have been building the feature layer over time by simply using the Append tool in ArcGIS Pro. Now I need to update some text and date values in these records, and I am finding that the computer takes over an hour to do so using arcpy.da.UpdateCursor. Is there a better way to do this? Do ArcGIS API for Python or ArcREST offer improved performance over an UpdateCursor? I have not used these libraries yet.
Here is an example of the simplest UpdateCursor that took over an hour - but I have many longer ones I need to do updating 10 or more fields for tens of thousands of records.
with arcpy.da.UpdateCursor(Layer,['date_approved', 'date_cleared']) as UC:
    for row in UC:
        if row[0] is not None:
            row[1] = row[0]
            UC.updateRow(row)


texts:
I have installed python 2.5.1 and FWTools2.4.7 in my laptop(Windows 7)
when i run my pythonscript in FWTools i get following error
Import error:No module named _gdal
can any one suggest how to solve this problem?


texts: Calculate volume of a building using LiDAR data

texts:
I would like to use 3D spatial metrics that introduced in the paper titled "Lidar and VHRS data for assessing living quality in cities - An approach based on 3D spatial indices". I have a lidar dataset including 14 las files and I have a shapefile including 416 buildings.
I wanted to calculate cubic volume of each building using LiDAR data. As I understand,I need to clip las catalog for each object in the shapefile and then run voxels on every clipped las file. I tried to do that but I could not understand how to get the volume calculation for each building. Is there any way to write volume value into the shapefile?


texts: MSSQL : How to build a geography with holes using the SqlGeographyBuilder

texts:
I'm trying to save a polygon with holes in it in a spatial layer using the SqlGeographyBuilder.

Using the below code, it doesn't take into account any holes.

How do i tell the SqlGeographyBuilder to start drawing a hole in the polygon?

case "polygon":
    {
        SqlGeographyBuilder gb = new SqlGeographyBuilder();
        gb.SetSrid(4326);
        gb.BeginGeography(OpenGisGeographyType.Polygon);
        gb.BeginFigure(vm.Coordinates[0].Latitude, vm.Coordinates[0].Longitude);
        foreach (var cooodinate in vm.Coordinates.Skip(1))
        {
            gb.AddLine(cooodinate.Latitude, cooodinate.Longitude);
        }
        gb.AddLine(vm.Coordinates[0].Latitude, vm.Coordinates[0].Longitude);
        gb.EndFigure();
        gb.EndGeography();

        DbGeography geog = DbGeography.FromText((gb.ConstructedGeography.STArea() < gb.ConstructedGeography.ReorientObject().STArea() ? gb.ConstructedGeography.ToString() : gb.ConstructedGeography.ReorientObject().ToString()), 4326);
        feature.Geog = geog;
        _unitOfWork.Save();
        break;
    }


texts:
I am an extreme neophyte with QGIS, so I maybe doing this all wrong. I have loaded the OpenLayers plugin as described in several online posts (and youtube video), and it is shown as listed in my plugins as installed (with the X mark). But, as shown in the youtube video, the plugin is not listed under the Plugins menu item, nor are any tools shown.

I am using QGIS Chugiak v2.4.0 on a 64-bit Windows machine.


texts: What is getting nearest points between two different but nearby (almost intersecting) datasets called?

texts:
I need some terminology help as I am still learning GIS. This is related to, but not the same as, What is getting nearest points between two different but overlaping datasets called?

I have two datasets.  One with the address of every fire and it's latitude/longitude of small census tract.  One with the location of every parcel of land (property) within that census tract.  I need to join or match the two datasets so I end up with a map of every property that has ever had a fire and a list of the properties without ever having a fire.

To help answer this question let me show two different datasets with this nearby situation:

Here is the original practice data from the City of Pittsburg that is over an 8 year period:

In it, every red point is a building fire, while every brown point is a building.  Each green point is the centriod of a relevant building. To clarify the green points are the centers of the brown regions (ie. I checked), its just that the green ones are of actual places people live).   So every fire building data point (red) don't exactly match up to the exact buildings.
The problem is because some of Pittsburg fire data (red) was anonymized and had some data collection issues. From my estimate about 20% of the Pittsburg fire data needs to be joined by the nearest polygon or green center point to it.

What exactly is this type of matching or joining between nearest lat/long of different datasets called in GIS terminology?

My concern is that with my GA fire sample data I might have about a similar problem when I receive the rest of the 5-8 year data.  So I am trying to figure out how to deal with that 20% possibility (even if its minor in the data I have right now)


texts:
Is there a (Windows) application that takes multiple, consecutive tracks, and merge them into a single GPX file?

I tried deadware "GPS Track Editor", but neither Merge nor Join did the trick:


texts: Creating columns in PostGIS using osm2pgsql hstore tags

texts:
While I feel this may come down to a database question, it feels GIS-related enough that I feel comfortable posting it here. I'll be happy to take it elsewhere if it belongs there, though!

I've started the arduous journey of mangling OSM data to work in an RDBMS format that is useful to me. One of the major problems I am running into is the hstore data type. While it is quite useful for being able to hold an extremely large dictionary of values, they are hard to "get at", in my experience. Specifically, the osm2pgsql tool creates addr:flats, addr:housenumber, and addr:interpolation fields, but leaves out some that are quite obvious to me, such as city, postcode, and street. I would like to parse this data to populate new columns with it. What methods would you recommend for accessing such data?

Thanks


texts: Extracting Netcdf climate data for multiple locations at different times using R?

texts:
I have daily sea surface temperature and air temperature data from 1980 to 2018. How do I extract values for multiple (up to 100) specific locations at different dates in one go to csv. Right now I can only extract one location at a time.

How can it be done in bulk in R?


texts:
Is there a function of arcpy that will set the layer's hyperlinks (as set in layer properties display tab) to a given field? Or another method I might set this for many layers at a time.

The field used will be the same for each layer.


texts: longitude and latitude values calculated using gdal differs from values displayed in various tools like (photoscan , global mapper etc..)

texts:
I have used gdal to read geotiff file. I tried to calculate longitude and latitude value for each pixel and line. Here is part of my code.

   for (int l_idx = 0;l_idx < line;l_idx++ )
{
    for (int p_idx = 0;p_idx < pixel;p_idx++)
    {
        if (elevation_data[(l_idx * pixel + p_idx)] > 0 && written_count < count)
        {
            x[(l_idx * pixel + p_idx)] = geo_transform[0] + geo_transform[1] * (p_idx) + geo_transform[2] * (l_idx);
            y[(l_idx * pixel + p_idx)] = geo_transform[3] + geo_transform[4] * (p_idx) + geo_transform[5] * (l_idx);
            x[(l_idx * pixel + p_idx)] += geo_transform[1] / 2.0;
            y[(l_idx * pixel + p_idx)] += geo_transform[5] / 2.0;

            fprintf_s(fp, "pixel and line %f,%f,", x[(l_idx * pixel + p_idx)], y[(l_idx * pixel + p_idx)]);
            coordinates_info->Transform(1, &x[(l_idx * pixel + p_idx)], &y[(l_idx * pixel + p_idx)]);
            fprintf_s(fp, "long and lat %f,%f, ", x[(l_idx * pixel + p_idx)], y[(l_idx * pixel + p_idx)]);
            fprintf_s(fp, "elevation %f\n", elevation_data[(l_idx * pixel + p_idx)]);

        }

    }
}

I saved my output into the file. I have imported same geotiff file in photo scan ,global mapper tool.I didn't(using above code) get any  longitude and latitude  or Easting and Northing value displayed in those tools. May I know the reason. Is there any mistake in my code or DO i need to any conversion or code?? or How to check the correctness of my output(correct long & lat)?


texts: Determine when a route enters and leaves a polygon

texts:
I have all of the segments of a route from White Horse, YK to Miami, FL. Now I want to determine how many miles I travel in each State / Province.

I have Polygons for each State and have been able to use Ray Casting to test whether or not each Segment is in a State but is there a faster way to determine when a route enters and leaves each State without testing ever single segment? Doing each segment would require over one hundred thousand lookups.


texts: Looking for a list of tools/add-ins for ArcGIS Explorer

texts:
I'm hoping to compile a list of add-ins & resources for ArcGIS Explorer.  I intend to install it at home, and apparently, there are a number of server feeds out there that can significantly augment the capabilities of this free ESRI software (IMS/WMS servers, toolbox items, etc, etc).

I thought many of you might be interested in having access to such a compilation as well, especially since it's a free software.  I'm quite anxious to see what's out there!!

Try your best to give at least a brief description of the add-in/server feed you post in your answers.

Thanks everyone


texts:
I have a bunch of geotiffs stored in a PostGIS DB that i am mounting on my OpenLayers map.

They all have the geoTransform information in them. I just want to know how i can mount them on my OpenLayers map with the correct affine transformation without running a gdalwarp on all the images. Though i wouldn't mind calling gdalwarp on the fly if possible. I just don't want to store those warped images in the DB (they take up 8x more space).

Can OpenLayers rotate my image based on the GeoTransform info in the GeoTiff? If so how?:

This is the code i'm using to overlay my tiff:

var bounds = MyOrg.UI.Map.Results.Layer.features[array_index].geometry.bounds;

var options = {
    numZoomLevels: 15,
    isBaseLayer: false,
    resolutions: MyOrg.UI.Map.getMap().layers[0].resolutions,
    maxResolution: "auto",
    projection: MyOrg.UI.Map.Projections.Mercator,
    strategies: [new OpenLayers.Strategy.Fixed()]
};

var graphic = new OpenLayers.Layer.Image(
    'test',
    'GranuleImage.ashx?granuleID=' + MyOrg.UI.Map.Results.Data.arrGranules[array_index].id + '&thumb=0',
    bounds,
    bounds.getSize(),
    options
);

MyOrg.UI.Map.getMap().addLayer(graphic);


texts:
I have a PostGIS query that throws an error somewhere down the line:

 SELECT ST_Disjoint(ST_Transform(location, 26910),ST_Buffer(ST_Transform(ST_SetSRID(ST_MakePoint(-121.3991535, 37.7458),4326), 26910), 100))

ERROR:  transform: couldn't project point (19.9597 50.0466 0): latitude or longitude exceeded limits (-14)

There may have been some error in point collection which is causing this error because there are points in location outside of the projection limits .

How can I skip records/points in my query that can't be projected??  In other words, how can I skip records/points that are outside of the projection limits?


texts: Trouble with creating new attributes on large shapefile?

texts:
I have a point shapefile with ~2.5 million points.  I am creating a new attribute converting ft to m, and QGIS (2.18.15) is having trouble working with the datasets.  I have seen suggestions to only do this through the main window, which sped up the creation of the attribute, but now the hangup seems to be saving the changes.  This has been occupying 30 of the write capacity of my hard drive for over 30 minutes now.

Do I just need to be more patient or is this an issue with QGIS?  I did this previously when I had access to ArcGIS and it went pretty quickly.


texts: Issues updating field values programmatically for particular field data type on 64bit Windows computer using QGIS 2.8.7

texts:
I have an issue with a plugin for QGIS 2.8 I'm developing when using it on a Windows 64bit machine (as opposed to 32bit).

I am trying to update values in a particular field in my plugin, and the values are not updated when running the plugin on the 64bit machine.

I notice that the Generic Type and Identified Type for this field are different when opening the project on each type of machine. The field has a Generic Type "int" and Identified type "Integer" on the 32 bit machine, while on the 64 bit machine the same field has a Generic Type "QString" and an Identified type "Integer64". The code below does work when updating a field of Generic Type "double" (Identified Type "real" on the 64 bit machine) and Generic Type QString (with identified Type "string"), the values are updated. However, I am wondering if I can also get it working for the former case.

UPDATE Something strange I am noticing on the 64bit computer, is that I am unable to keep a field in the Generic Type "double" nor Generic Type "int". Once this field is created in int or real--when the project is saved, closed, and reopened, the field Generic Type becomes QString and Integer64 once again (the data type that I am unable to update via my plugin).

UPDATE #2 Further testing has shown that the values are updated and the field type stays "int" when editing on a 64bit machine using QGIS 2.8.3. The issues I had were with another 64bit machine running QGIS 2.8.7. Could this be an version issue instead of a system issue?

Below is the code I used to update the values:

    self.Nature_idx = self.resLayer.fieldNameIndex('MATERIAU')
    #This field below is the problematic one, with inconsistent data types on 32 vs 64 bit computers
    self.Diam_idx = self.resLayer.fieldNameIndex('DIAMETRE')
    self.resMain_idx = self.resLayer.fieldNameIndex('MAINT')
    self.fuitesNature_idx = self.fuitesLayer.fieldNameIndex('NaturCanal')
    #This field below is the one carrying the value I would like to carry over and update in the previous field. It has a generic type "Double" and Identified Type "real"
    self.fuitesDiam_idx = self.fuitesLayer.fieldNameIndex('DiamCanal')
    self.fuitesOrigin_idx = self.fuitesLayer.fieldNameIndex('OriginInfo')

    self.fuitesSelectList = []
    self.fuitesNatureList = []
    self.fuitesDiamList = []
    self.fuitesOriginList = []

    for feature in self.fuitesLayer.getFeatures():
        if feature.id() == fid:
            self.fuitesSelectList.append(feature.id())
            self.fuitesNatureList.append(feature.attributes()[self.fuitesNature_idx])
            self.fuitesDiamList.append(feature.attributes()[self.fuitesDiam_idx])
            self.fuitesOriginList.append(feature.attributes()[self.fuitesOrigin_idx])

    caps = self.resLayer.dataProvider().capabilities()
    fid = self.selectList[0]
    newNature = self.fuitesNatureList[0]
    newDiam = self.fuitesDiamList[0]

    #FYI: I get integers here
    #QMessageBox.information(None, 'DEBUG', newDiam)

    self.resMaint_idx = self.resLayer.fieldNameIndex('MAINT')
    if caps & QgsVectorDataProvider.ChangeAttributeValues:
        attrs = {self.Nature_idx : newNature, self.Diam_idx : newDiam, self.resMaint_idx : 'fuite_reparee'}
        self.resLayer.dataProvider().changeAttributeValues({fid : attrs})


texts:
I'm trying to reproject a .tiff file from a wgs84 projection to a south america albers equal area conic projection. To do this I have created the following script (it will later be modified for a high performance computing task hence the need for a new function):
# Packages
library(raster)

# Define the projection- South America Albers Equal Area Conic
saaeac <- '+proj=aea +lat_1=-5 +lat_2=-42 +lat_0=-32 +lon_0=-60 +x_0=0 +y_0=0 +ellps=aust_SA +units=m +no_defs'

# Define the reprojection function
reproject <- function(input, output){
  mb <- raster(input)
  mb_p <- projectRaster(mb, crs = saaeac)
  writeRaster(mb_p, output)
}

reproject(input="/rds/general/user/jaj20/home/jaj20/MapBiomas_1997.tif", output="/rds/general/user/jaj20/home/jaj20/MapBiomas_1997p.tif")

When I run the script however, the following errors are returned:
Error in projectRaster(mb, crs = "+proj=aea +lat_1=-5 +lat_2=-42 +lat_0=-32 +lon_0=-60 +x_0=0 +y_0=0 +ellps=aust_SA +units=m no_defs") :
  input projection is NA
Calls: reproject -> projectRaster
In addition: Warning message:
In .newCRS(value) :
  +proj=longlat +datum=WGS84 +no_defs is not a valid PROJ.4 CRS string
Execution halted
mv: cannot stat 'MapBiomas_????.tif': No such file or directory

Why is this the case? I copied the projection string from epsg so I know it to be correct.


texts: Obtain channel/stream network from a DEM in QGIS

texts: This might sound like a basic question, but running SAGA's Strahler order won't do the job for me. The lengths of the streams look way too short and this is quite a mountainous area. The result makes sense, but I need to see the the longer networks. If I can get large any size of basin with r.watershed, why can't I see the corresponding channel network? These only correspond to very small basins.


texts:
I use QGIS and don't know why there is a North indicator (see the figure).

Anyone know how did it come and how to remove it?


texts: How to set scale-rule based renderer in QGIS using Python

texts:
I'm trying to create a QGIS app using Python, where some layers are only visible at certain scales.

My code for the renderer is:

    ROAD_RULES = (('Freeway', "CLASS_CODE = 0", 'dark green', 0.2, Qt.SolidLine, None),
        ('Main', "CLASS_CODE IN (1, 2, 3)", 'black', 0.2, Qt.SolidLine, None),
        ('Minor', "CLASS_CODE IN (4, 5)", 'dark gray', 0.15, Qt.SolidLine, (0, 2000000000)),
        ('Track', "CLASS_CODE IN (6, 7, 11, 12)", 'dark gray', 0.15, Qt.DashLine, (0, 2000000000)),
        ('Proposed', "CLASS_CODE = 9", 'gray', 0.15, Qt.DotLine, (0, 2000000000)))

    # roads
    # create new rule-based renderer and get root rule
    symbol = QgsSymbolV2.defaultSymbol(self.roadLayer.geometryType())
    renderer = QgsRuleBasedRendererV2(symbol)
    root_rule = renderer.rootRule()

    # based on rules, create clones and set their details
    for label, expression, color_name, width, style, scale in ROAD_RULES:
        rule = root_rule.children()[0].clone()
        rule.setLabel(label)
        rule.setFilterExpression(expression)
        rule.symbol().deleteSymbolLayer(0)
        symbol_layer = QgsSimpleLineSymbolLayerV2()
        symbol_layer.setColor(QColor(color_name))
        symbol_layer.setWidth(width)
        symbol_layer.setPenStyle(style)
        rule.symbol().appendSymbolLayer(symbol_layer)
        if scale is not None:
            rule.setScaleMinDenom(scale[0])
            rule.setScaleMaxDenom(scale[1])
        root_rule.appendChild(rule)

    # delete default rule and apply renderer
    root_rule.removeChildAt(0)
    self.roadLayer.setRendererV2(renderer)

This works OK as long as I only zoom in and out of my map canvas using self.mapCanvas.zoomIn() and self.mapCanvas.zoomOut() methods.  But if I zoom in and out using the mouse wheel, or if I pan, the scale changes and the maximum and minimum scales no longer work as intended.

When I say 'the scale changes', I mean that even when the map later looks the same on screen, mapCanvas.scale() will report a different scale figure (sometimes orders of magnitude different).

Is there a way to stop this, so that scale-dependent layers operate consistently in my app as they would in QGIS itself?


texts: Calculate average of the temperature value from 4 surrounded points in PostGIS

texts: There are two-point tables. One includes temperature value and the other does not have this information. I would like to calculate the average of the temperature for the 2nd table based on the four closest points around the first ones. I wonder how can I do it in PostGIS.


texts: How to access ExtendedData/SimpleData from KMLfile using openlayers 2.x?

texts:
Is is possible to access  SimpleData tag using Openlayers 2.x library? I have a KML file exported from QGis and it stores description data SimpleData tag "<>". i know that there is a way getting the name and the description tag like this :feature.attributes.name , feature.attributes.description. the kml attribute data are like this:SimpleData> name="NAME_0">Syria /SimpleData> so i suppose that i should call it by the name of the attribute.

Any help ?


texts: How to easily toggle outlines and fills for polygon layers (in QGIS 3.x)?

texts:
So I want to be able to easily toggle outlines and fills for polygon layers in the Layers List.
I managed to make two rules (both with the rule 'true') and then they will show up in the Layers List:

But... I want to make sure that all of the symbol options stay in the legend after unchecking them. What happens now is that when I uncheck the fill, it disappears from the layers list, and then I need to double click to go back into the symbology settings and check it from there. But... that kinda defeats the point of having a convenient way to toggle the outlines and fill styles independently.
Are there any easy solutions for this that I couldn't find in the settings? Or shall I open a feature request?


texts: nodes geographic coordinates with pgRouting Query

texts:
Can we get the coordinate of nodes belongs to shortest path???

I give brief description, I calculate the shortest path with A* algorithm using Pgrouting, but what I'm looking for coordinates (long, Lat) of all the nodes that belong to path???


texts:
I have created the layout directory
C:\Program Files (x86)\GeoServer 2.15.1\data_dir\layouts

In the directory I have created an xml file with the following:

<layout>
    <decoration type="image" affinity="bottom,right" offset="6,6" size="80,31">
        <option name="C:\Program Files (x86)\GeoServer 2.15.1\data_dir\layouts" value="Uis_Extensions.png"/>
    </decoration>

    <decoration type="scaleline" affinity="bottom,left" offset="36,6"/>

    <decoration type="legend" affinity="top,left" offset="6,6" size="auto"/>
</layout>

The expectation is that when you go to layer preview the decoration should show overlayed on the WMS. I am getting none of that. What am I doing wrong?


texts: Add a nodata mask or alpha band to read-only image?

texts:
How to add a nodata mask or alpha channel to a source image that is read-only or should be treated as such?
For example I have this satellite image that uses jpeg compression in a geotiff. Since the source already uses a lossy compression we don't want to run it through any process that will resample the valid pixels again, thus losing more information.  We also don't want to pay the storage penalty of converting to another file using lossless compression because a) we're duplicating the source, and b) the lossless-file will be several times larger than the source.

Yellow is current nodata (actually 0,0,0). Black jagged perimeter is what we want to turn into nodata.
The gdal nearblack utility is a very convenient method of removing this unwanted collar in one step but we could also use a polygon mask or some other method. Unfortunately all of the methods I'm aware of create redundant copies of the data and/or resample our value pixels which we want to leave unchanged.


texts: Error message while running Strahler order for river delineation

texts:
I wanted to make river network using Strahler order in SAGA with filled DEM as the input, but there is an error message keep appearing:
QGIS version: 3.14.0-Pi
QGIS code revision: 9f7028fd23
Qt version: 5.11.2
GDAL version: 3.0.4
GEOS version: 3.8.1-CAPI-1.13.3
PROJ version: Rel. 6.3.2, May 1st, 2020
Processing algorithm…
Algorithm 'Strahler order' starting…
Input parameters:
{ 'DEM' : 'C:/Users/PERSONAL/Desktop/QGIS_results_HDP/fill_34.sdat', 'STRAHLER' : 'C:/Users/PERSONAL/Desktop/QGIS_results_HDP/river_34.sdat' }

io_gdal 0 -TRANSFORM 1 -RESAMPLING 3 -GRIDS "C:/Users/PERSONAL/AppData/Local/Temp/processing_KZPmeS/bbc1c29883054780acfc3d030361a759/fill34.sgrd" -FILES "C:/Users/PERSONAL/Desktop/QGIS_results_HDP/fill_34.sdat"
ta_channels "Strahler Order" -DEM "C:/Users/PERSONAL/AppData/Local/Temp/processing_KZPmeS/bbc1c29883054780acfc3d030361a759/fill34.sgrd" -STRAHLER "C:/Users/PERSONAL/Desktop/QGIS_results_HDP/river_34.sdat"

C:\Users\PERSONAL\Documents>set SAGA=C:/PROGRA~1/QGIS3~1.14/apps\saga-ltr
C:\Users\PERSONAL\Documents>set SAGA_MLB=C:/PROGRA~1/QGIS3~1.14/apps\saga-ltr\modules
C:\Users\PERSONAL\Documents>PATH=C:\PROGRA~1\QGIS3~1.14\apps\qgis\bin;C:\PROGRA~1\QGIS3~1.14\apps\Python37;C:\PROGRA~1\QGIS3~1.14\apps\Python37\Scripts;C:\PROGRA~1\QGIS3~1.14\apps\qt5\bin;C:\PROGRA~1\QGIS3~1.14\bin;C:\Windows\system32;C:\Windows;C:\Windows\system32\WBem;C:\PROGRA~1\QGIS3~1.14\apps\Python37\lib\site-packages\pywin32_system32;C:\PROGRA~1\QGIS3~1.14\apps\Python37\lib\site-packages\numpy\.libs;C:/PROGRA~1/QGIS3~1.14/apps\saga-ltr;C:/PROGRA~1/QGIS3~1.14/apps\saga-ltr\modules
C:\Users\PERSONAL\Documents>call saga_cmd io_gdal 0 -TRANSFORM 1 -RESAMPLING 3 -GRIDS "C:/Users/PERSONAL/AppData/Local/Temp/processing_KZPmeS/bbc1c29883054780acfc3d030361a759/fill34.sgrd" -FILES "C:/Users/PERSONAL/Desktop/QGIS_results_HDP/fill_34.sdat"
____________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
____________________________

SAGA Version: 2.3.2 (64 bit)
____________________________
library path: C:\PROGRA~1\QGIS3~1.14\apps\saga-ltr\modules\
library name: io_gdal
library : GDAL/OGR
tool : Import Raster
author : O.Conrad (c) 2007 (A.Ringeler)
processors : 8 [8]
____________________________

Parameters

Grids: No objects
Files: "C:/Users/PERSONAL/Desktop/QGIS_results_HDP/fill_34.sdat"
Select from Multiple Bands:
Alphanumeric Sorting: yes
Transformation: yes
Resampling: B-Spline Interpolation

loading: C:/Users/PERSONAL/Desktop/QGIS_results_HDP/fill_34.sdat

Driver: SAGA
Bands: 1
Rows: 19998
Columns: 96657
loading: fill_34

C:\Users\PERSONAL\Documents>call saga_cmd ta_channels "Strahler Order" -DEM "C:/Users/PERSONAL/AppData/Local/Temp/processing_KZPmeS/bbc1c29883054780acfc3d030361a759/fill34.sgrd" -STRAHLER "C:/Users/PERSONAL/Desktop/QGIS_results_HDP/river_34.sdat"

C:\Users\PERSONAL\Documents>exit
Execution completed in 691.15 seconds
Results:
{'STRAHLER': 'C:/Users/PERSONAL/Desktop/QGIS_results_HDP/river_34.sdat'}

Loading resulting layers
**The following layers were not correctly generated.
• C:/Users/PERSONAL/Desktop/QGIS_results_HDP/river_34.sdat
You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.**

I have tried to change the target directory to the C:/ or D:/ to prevent the long directory name, if that is the problem, but the message keep appearing and I still cannot run the Strahler.
PS. I am using QGIS 3.14 and 3.16.


texts: First time using i.segment.uspo in GRASS

texts:
I am using the module i.segment.uspo in GRASS for the first time and I wanted to know if anyone here haws experience with that add-on in GRASS so I can get some hints on how to use it more efficiently, because I could use i.segment but I feel like it wouldn't be the same because, for once, I don't know how to properly set up meanshift. So, can anyone help me understand better how i.segment.uspo works or are you as new to this as I am?

P.S If someone could point me to some examples, or even share, other than the ones that are in the documentation file, that would be great.


texts: Selection of lines that lays mostly in defined zone (majority overlap)

texts:
I have lines of surface runoff and I need to select the ones that are lying mostly in internal buffer of polygon. I'm trying to make some script in python, but it's very important to make the script simpliest as it's possible, because there are 30 milions of this lines (and some about 650 000 polygons) and this is only a part of whole problem. Has somebody any suggestion?
I'm using ArcGIS 10.2.

There are the raw data.


texts: Mismatch between two methods projecting hdf5 using GDAL

texts:
I want to convert subdataset from h5 to tiff format and project the data with WGS84.
The product is https://www.geo.vu.nl/~gwerf/GFED/GFED4/ name is GFED4.1s_2010.hdf5 and the subdataset is sd39 without coordinate information stored.
gdalinfo GFED4.1s_2010.hdf5 -sd 39

Driver: HDF5Image/HDF5 Dataset
Files: GFED4.1s_2010.hdf5
Size is 1440, 720
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  720.0)
Upper Right ( 1440.0,    0.0)
Lower Right ( 1440.0,  720.0)
Center      (  720.0,  360.0)
Band 1 Block=180x45 Type=Float32, ColorInterp=Undefined
  Metadata:
    burned_area_01_burned_fraction_long_name=GFED4s burned fraction. Note that this INCLUDES an experimental "small fire" estimate and is thus different from the Giglio et al. (2013) paper
    burned_area_01_burned_fraction_units=Fraction of grid cell

I got some code from other answers listed in the community and tried 2 methods to project the data into WGS84 coordinate system.
Method 1
I look through the study area from https://daac.ornl.gov/VEGETATION/guides/fire_emissions_v4_R1.html and use gcp function.
gdal_translate -of VRT -a_srs EPSG:4326 -gcp 0 0 -180 90 -gcp 1440 0 180 90 -gcp 0 720 -180 -90 -gcp 1440 720 180 -90 HDF5:" GFED4.1s_2010.hdf5"://burned_area/01/burned_fraction neu.vrt
gdalwarp -of GTIFF -tps neu.vrt test2.tif

Method 2
Because I found there are two datasets containing lon and lat information, I create vrt from subdatasets.
  SUBDATASET_704_NAME=HDF5:"GFED4.1s_2010.hdf5"://lat
  SUBDATASET_704_DESC=[720x1440] //lat (32-bit floating-point)
  SUBDATASET_705_NAME=HDF5:"GFED4.1s_2010.hdf5"://lon
  SUBDATASET_705_DESC=[720x1440] //lon (32-bit floating-point)

The code is
gdal_translate -of VRT HDF5:"GFED4.1s_2010.hdf5"://lat lat.vrt
gdal_translate -of VRT HDF5:"GFED4.1s_2010.hdf5"://lon lon.vrt
gdal_translate -of VRT HDF5:"GFED4.1s_2010.hdf5"://burned_area/01/burned_fraction burn.vrt
gdalwarp -geoloc -t_srs EPSG:4326 -overwrite burn.vrt test4.tif

And finally, when I export the two figure into Arcmap, they have a mismatch with each other, like this

Because I am new in dealing with satellite data and I guess maybe the method 2 is more accurate? I wonder why the mismatch occurs and which should be used for my following analysis?


texts: Can targets placed on the facades of a building be used to increase the accuracy of a point cloud of a building generated from Drone/UAV/UAS imagery?

texts: Can targets (used for registering point clouds) mounted on a building be used to increase the accuracy of a point cloud generated from UAV imagery?


texts: How to create MPK using ArcGIS Runtime WPF?

texts:
I would like to create a Map PacKage with ArcGIS Runtime WPF.

Have you an example ?


texts: How can I add OpenLayers layer from the Python console?

texts: I am new to the Python console. How can I add a layer from the OpenLayers plugin from the Python console?


texts: How to use Esri custom geographic transformation in FME Desktop 2015 Workbench?

texts:
I am trying to reproject DWG files from one projection system to other using a custom transformation created by me with ArcGIS.

When I execute it reprojects the data but does not use the custom transformation parameters.


texts: Clip raster by mask layer gdalwarp error in QGIS

texts:
I've been trying to find why my gdalwarp isn't allowing me to clip a raster to the shape of a shapefile.

GDAL command:
gdalwarp -ot Byte -of GTiff -cutline /Users/johnmorris/Desktop/Geodata/StudyArea.shp -crop_to_cutline -dstnodata 0.0 /Users/johnmorris/Desktop/Geodata/TotalRevenueDraft. /var/folders/z4/bpntn2_964d4hvl5l37wfny80000gn/T/processing_cce36f0eaee64b1aa98a58af4c157ad2/9ca388a9371f44deb2fa4eab49ba0026/OUTPUT.tif
GDAL command output:
/bin/sh: gdalwarp: command not found

Execution completed in 0.08 seconds
Results:
{'OUTPUT': }

Loading resulting layers
The following layers were not correctly generated.
/var/folders/z4/bpntn2_964d4hvl5l37wfny80000gn/T/processing_cce36f0eaee64b1aa98a58af4c157ad2/9ca388a9371f44deb2fa4eab49ba0026/OUTPUT.tif
You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.

It keeps saying it doesn't understand the command gdalwarp. I have GDAL 2.3.1 installed and am using QGIS 3.0.3

Please help!


texts: Adding short CRS description to composer template?

texts:
How can I add CRS description as automated text for composer template? so anytime someone open a project, it uses the project CRS to print description of CRS like: WGS84 / Mercator 1

@project_crs gives me ESPG code, @project_crs_definition gives me more than I want in a different format I want.

All I want is simply WGS84 / Mercator 1 or UTM Zone.


texts: How to animate daily data in QGIS time manager plugin by changing color of polygon instead of plotting markers

texts:
I 'm a beginner in QGIS in general and geo spatial data visualization.
Given a vector type map (.tab file) of districts geographically spanning my business area and csv file containing weekly sales volume per district. I tried to animate the weekly sales per district using QGIS time manager plugin.
How to animate data from the csv by changing the vector map color , instead of plotting separate colored markers ?

I expected my result would be like this :
https://underdark.files.wordpress.com/2011/09/thl_population.gif


texts: Pure python library to handle raster data (not GDAL, rasterio)

texts:
I would like to do some basic operations over some downloaded raster files  (tiff) and include the process in a package, but I don't want the users of my package to install GDAL, neither rasterio.

In the process I basically download every band of a multiband raster separately, and need to merge all into one single raster tiff. I know it can be done with the mentioned packages, but I need a pure python solution for this simple task. I was thinking in scikit-image but I am not sure how it will be for handling geographic information, mainly CRS.


texts: Upgrading from Qgis 2.4 to 2.12 on Ubuntu 15.04

texts: This should be a simple question. On Ubuntu 15.04 how can I get Qgis upgraded from 2.4 to 2.12?


texts: Get field attributes from Esri JavaScript API query

texts:
I'm trying to get the geometry in an array from a queryTask using the Esri JavaScript API v4.7. I'm sure I'm missing something but I dont know what.  I want to take geometry from this to then create a new record in another service.

require([
"esri/tasks/QueryTask",
"esri/tasks/support/Query"
],
function(
QueryTask,
Query
){
  var LayerUrl = "https://services8.arcgis.com/R1PAHHLiCM5plgJt/arcgis/rest/services/MICC_Planning/FeatureServer/0";
  var queryTask = new QueryTask({
    url: LayerUrl
  });
  var lotp="15MPH7941";
  var q="Lotplan LIKE '"+lotp+"'";
  var query = new Query();
  query.returnGeometry = true;
  query.outFields = ["*"];
  query.where = q;
  queryTask.execute(query).then(function(results){
    console.log(4);
    window.alert(results.features.geometry);
  });

  // When resolved, returns a count of the features that satisfy the query.
  queryTask.executeForCount(query).then(function(results){
    console.log(results);
  });

});


texts: How to remove line overlap hierachically in PostGIS with ST_Difference?

texts:
I try to remove the overlap of line geometries in one table  based on a class attribute. My table consists of two columns containing line geometries and a number "class" from 1 to 3.
CREATE TABLE lines(class INT, geom geometry );

INSERT INTO  lines VALUES
    (1, 'SRID=4326;LINESTRING(1 1, 1 4)'),
    (2, 'SRID=4326;LINESTRING(1 3, 1 6)'),
    (3, 'SRID=4326;LINESTRING(1 5, 1 9)');

Now I'd like to extract the part of geometry 2 wich is not intersecting with geometry 1 and the part of geometry 3 wich is not intersecting with geometry 2. I did this already in Python:
class_1 = shp[(shp.Class == 1)]
class_2 = shp[(shp.Class == 2)]
class_3 = shp[(shp.Class == 3)]

part_1 = gpd.overlay(class_2,class_1, how='difference')
part_2 = gpd.overlay(class_3,class_2, how='difference')
part_3 = class_1

but running following code, the query comes to no end:
SELECT
CASE WHEN lines.class = 2 THEN ST_Difference(lines.geom, (SELECT ST_Union(lines.geom) FROM lines WHERE class = 1 ))
     WHEN lines.class = 3 THEN ST_Difference(lines.geom, (SELECT ST_Union(lines.geom) FROM lines WHERE class = 2 ))
ELSE lines.geom
END
FROM lines;

How can I write the code most clever?


texts: Is the cell size dependent on the number of features for the Euclidean Distance tool?

texts:
I have attempted to run the Euclidean Distance tool, using ArcMap 10.5.1, on a feature class of 50 points (hospitals), but I received an invisible output. Using 'Identify' I get 'NoData' for the pixel values.

Is the cell size dependent on or limited by the number of features in the input FC?

What I think is more likely is the output cell size is causing the output issues. I want to do a Weighted Sum on my raster data, and I am trying to make the cell sizes in Euclidean Distance consistent for the raster math. I have grocery store and hospital FCs as separate inputs/outputs. I don't understand why this> is happening, but the default output cell sizes are drastically different for these two point FC inputs. One is 000s, and the other is 0.0000000 -E^3. I have tried to find an output cell size that works for both but haven't been successful (hospitals output raster is invisible/NoData).


texts: One-way ArcGIS replication copying __ATTACH tables instead of base tables

texts:
I'm trying to get one-way replication working in ESRI ArcGIS 10.  Here's a quick summary of the steps I followed:

Open map in ArcMap with the relevant data layers I want to replicate
Choose Create Replica from distribution toolbar
Choose One-Way Replica, Parent-to-Child
Use Simple Model (as opposed to Full Model)

This works for me with one problem.  Several of the featuresets that the data layers are based  off of have __ATTACH tables.  When I perform the replication, the __ATTACH tables get replicated but the base tables do not.

For example, there is a layer based on a featureset called Parcels that may have attachments stored in a SQL Server table called Parcels_ATTACH.  In the replicated copy of the geodatabase, there is a Parcels_ATTACH table but no Parcels table.

If it makes any difference, the Parent Geodatabase is stored in SQL Server 2008 R2 and the replicated child is stored in SQL Server Express 2008 R2.


texts: How to install PostGis on EC2 micro-instance

texts:
I tried to install PostGIS on my Amazon EC2 microinstance using sudo yum install postgis2 and got following error message:

$ sudo yum install postgis2
Loaded plugins: priorities, security, update-motd, upgrade-helper
amzn-main                                                | 2.1 kB     00:00
amzn-updates                                             | 2.3 kB     00:00
Setting up Install Process
No package postgis2 available.
Error: Nothing to do

How can I install PostGIS on an EC2 microinstance?


texts: Is there a way to add drop-down list of predefined variables to R-script menu in SEXTANTE QGIS?

texts: I would like to create a drop-down list of variables (similar to what it is done for layers input) for my R script menu in SEXTANTE. Is it possible?


texts: QGIS measure tool window has disappeared?

texts:
I'm using QGIS 2.18.14. I'm able to measure with the tool, but there is no dialog box/window with the actual results.

Where is it hiding?


texts: Combining polygon and raster using ArcGIS Desktop?

texts: How do I combine the raster data and county outline polygon so that only the elevation data within the county boundaries is displayed?


texts: Feature becomes invalid when reading through WFS/GML in GeoServer

texts:
I'm working with the TIGER shapefiles, specifically the Tract layer, and my I'm trying to read selected features via WFS using GeoServer, then save those to a new layer. However, when the shapes come out of GeoServer via WFS, they appear to get altered in a way that makes a significant number of them invalid.

In the below image, I have the exact same feature (Census Tract 41047010701). The green feature is it coming directly from PostGIS (using QGIS Database connection). The red is the exact same feature, just passed through GeoServer and output using it's WFS service. It's the exact same shape being read from the underlying PostGIS table, I haven't even copied the feature at this point.

I thought it could be related to reprojecting from TIGERS native epsg::4269 to epsg::4326, but I just turned that off in GeoServer and I get the same result.

This is just one example, but it happens to a ton. Below is an output of running ST_IsValidReason() on a table after features have been read through WFS and inserted into another table.


texts: GeoServer 2.13.3 unable to delete from Oracle 11gR2 over WFS-T

texts:
I'm attempting to upgrade GeoServer 2.9.2 to 2.13.3. GeoServer uses the Oracle Data Store plugin (updated to gt-jdbc-oracle-19.3 as part of the upgrade). Oracle is 11gR2.

One of the applications using GeoServer creates and deletes features in an Oracle-based layer over WFS-T (version 1.0.0). With 2.9.2 it can successfully create and delete. However with 2.13.3 it can still create but cannot delete.

GeoServer logs the following error:

Caused by: java.io.IOException: Error occured calculating bounds for LAYER_NAME
at org.geotools.jdbc.JDBCDataStore.getBounds(JDBCDataStore.java:1345)
at org.geotools.jdbc.JDBCFeatureSource.getBoundsInternal(JDBCFeatureSource.java:513)
... 118 more
Caused by: java.sql.SQLSyntaxErrorException: ORA-00933: SQL command not properly ended

Getting the exact SQL from Oracle requires me to jump through a few hoops within our organisation, so I don't expect to have that information until next week some time. Is this a known issue? Are there compatibility issues between this version of GeoServer and Oracle 11gR2?


texts:
We have some shapefile data and we tried to render them according to scope (boundingbox) from user.

For example, the client may send request like this:

layer=river&width=256&height=256&bbox=....

Then I will generate an image accordingly.

I know this is similar to a WMS server, however I do not want to use GeoSever or something else. They are heavy. I just have two shapefiles  with a fixed style.

And I found that GeoSever use GeoTools as the render library, so I tried to use GeoTools myself, however I found that GeoTools does not provide an acceptable performance.

So I wonder if there are any alternative solutions?


texts:
I am able to add and remove cartoDB layers over my google maps but after removing layer the legend stays . how to remove that legends ?
 here is code that i used

document.getElementById('school').onclick = function()
        {
            var layerUrl = 'http://astrosnehakul123.cartodb.com/api/v2/viz/4c13821a-e192-11e3-8f8f-0e73339ffa50/viz.json';
            if (school.checked)
            {
                        // add cartodb layer
                      cartodb.createLayer(map, layerUrl)
                      .addTo(map,0)
                      .on('done', function(layer)
                      {
                        layer1 = layer;
                        var sublayer = layer.getSubLayer(0);
                        sublayer.on( function(e, pos, latlng, data) {
                          cartodb.log.log(e, pos, latlng, data);

                    });

                    sublayer.on('error', function(err) {
                      cartodb.log.log('error: ' + err);
                    });

                  })
                  .on('error', function() {
                    cartodb.log.log("some error occurred");
                  });
        }
        else{

            layer1.remove();
            layer1.clear();
        }
    }


texts: QGIS Difference Algorithm: multilinestring won't work with multipolygon layer

texts:
I am working with two vector layers in QGIS (location: Indonesia). I want to 'subtract' the area of layer 1 ('Forest_buffer') of layer 2 (Scenario_1). I am using the Difference algorithm. Layer 2 is the input layer, and layer 1 is the overlay layer.
I've used this algorithm before, but when I use layer 1 ('Forest_buffer'), the algorithm always fails after about 7%. The error I get says something about that a multilinestring can't be combined with a multipolygon layer:
"Object couldn't be written to Difference_19c76041_d4e5_4481_97ed_a76ef3cf8d3d: Could not add object of geometrytype MultiLineString to layer of MultiPolygon."
How to fix this? I already used the 'Fix geometries' on both layers, and the buffer trick. Both don't work. Could I perhaps somehow convert multilinestrings in the layer to polygons?
Here is the wetransfer link to the two data-sets: https://we.tl/t-gshlI4TbCU
Hopefully someone can help this QGIS beginner out!


texts: Convert geodata from .rdb and .mpr files

texts:
I've got an COM-server and activeX client (looks like delphi app), that display some geodata.
Data stores in files with .rdb and .mpr extensions.

I find out that rdb is oracle database file format. There is no running oracle server on my machine, but application still works.

Is there any tools which could help me convert this data to another preferred format (.shp or .tab or something else)?


texts: Downloading pitem file from ArcGIS Online as shapefile using QGIS

texts:
I'm wanting to export the data connected to a "pitem" file (an ArgGIS "portal item" file).
Here in ArcGIS Online, you can click on "Open in ArcGIS Desktop", and it downloads the pitem file. I'm able to upload it into ArcGIS. I want to export the data as a shapefile.
Data:
https://www.arcgis.com/home/item.html?id=d0559466f6994dd4b8e79298018b2fb8
I'd really like to do this in QGIS.  Currently, I'm having to download ArcGIS just to obtain this data.
I found instructions on how to open it in QGIS at How to open ArcGIS online data using QGIS?, but when I went to Layer > Add Layer > Add ArcGIS FeatureServer Layer and copied and pasted the URL given on the bottom of the right side of the page, and clicked "connect", it said "The query returned no layers".


texts: How can I change start and end point of a linestring?

texts:
I have a LINESTRING from a pgRouting query. The start point and/or end poit of this line does not connect exactly the points I've clicked on the road and give me all the road segment of that edge (passing the point to next vertice) or stops in the nearest vertice (before the point).

I need to change the points to connect the line result (red lines on the images) near the points I've clicked on the map (blue circles).

My Kung-pgRouting-fu is like this tutorial http://workshop.pgrouting.org/chapters/wrapper.html#return-route-with-network-geometry (function pgr_fromAtoB()) and is working fine.

The function pgr_fromAtoB() in tutorial gets the nearest roads (edges) from the points (clicked on web interface an can be arbitrary locations) and connect its vertices using pgRouting.

EDIT: I'm testing this info now: How to find the nearest point projected on the road network?. I think this could help me.


texts: Best way to store polygon geospatial data in postgreSQL?

texts:
I have to store longitude and latitude coordinates for several places in a PostgreSQL db table. Each row in this table has a field with a set of pairs {latitude1, longitude1, ....latitudeN, longitudeN} wich represent an enclosing polygon.
I stumbled upong PostGIS, and I'm wondering if it is necessary to use Polygon datatype for this particular field in my table, like depicted in this example
I'm asking because the only thing I'm gonna do with those polygons is drawing them afterwards in a Google Maps, I don't really think I'm gonna do any reasoning on them (yet)...thus, would be PostGIS overkill for this? Besides all the time I have to spend in installing and learning PostGIS...
Any toughts?
Thanks!


texts: How to get all gps points on a road from osm data?

texts:
I'm trying to extract gps points on roads in a region. For each road, I need a gps point on that road every 10 meters.

I'm aware that Google Map API offer such way points, but due to the latency-sensitive nature of this application and amount of data needed, I cannot afford to do network communications to Google server.

I came across OSM this wonderful opensource map data. But as a newbie, it's very unclear how I can extract gps points on roads from OSM. I know roads are stored as "line string" with keyword "highway". Some geojson file includes start and end gps locations for each road. But how do I extract all gps points (Lat,Lng) along the road?


texts: GeoJSONTopoJSON roundtrip breaks leaflet on antimeridean


texts:
I'm implementing an application that uses GeoServer deployed on apache tomcat. I would like to know if there is a way to handle an error case, such as when GeoServer is hung. For example, I'd like to display a message to user or something like that on client side.

I've been told that :

If GeoServer hangs, you will receive no answer. You should set a read
time out for your GeoServer HTTP requests. After the timeout, you can
send back whatever you want.

but since I'm pretty beginner at these stuff, I've no idea how to do that. would anyone please explain me more how to set a readtime out for your GeoServer HTTP requests ??

EDIT:
Can I convert the configuration in proxyhost to message displays to the user ?


texts: Splitting second word of attribute using Python Parser of ArcMap Field Calculator?

texts: This is Field Calculator in ArcMap. I am trying to split the word after the first space if the words starts with "San" like San Francsico, San Diego and etc. It return an error. Below is the script.


texts: Using Conditional Statement in SLD GeoServer

texts:
Is there any way to use the conditional statement in SLD? somehow using if_then_else like this Pseudocode:

if (10 =< variable =< 100) {
    then Multiply variable by 2
} else {
    then Multiply it by 4
}

Or if it is also possible to have a multiple if_then_else

if (10 =< variable =< 100) {
    then Multiply variable by 2
} if (100 < variable =< 200) {
    then Multiply it by 3
} else {
    then Multiply it by 4
}

Example.sld

<ogc:Filter>
    <ogc:Function name="if_then_else">
        <ogc:PropertyName>variable</ogc:PropertyName>
        <!-- check if the "variable" is between 10 and 1000 -->
        <ogc:PropertyIsBetween>
            <ogc:LowerBoundary>
                <ogc:Literal>10</ogc:Literal>
            </ogc:LowerBoundary>
            <ogc:UpperBoundary>
                <ogc:Literal>100</ogc:Literal>
            </ogc:UpperBoundary>
        </ogc:PropertyIsBetween>

        <!-- if it is between this range then multiply it by 2 -->
        <ogc:Mul>
            <ogc:PropertyName>variable</ogc:PropertyName>
            <ogc:Literal>2</ogc:Literal>
        </ogc:Mul>

        <!-- else multiply it by 4 -->
        <ogc:Mul>
            <ogc:PropertyName>variable</ogc:PropertyName>
            <ogc:Literal>4</ogc:Literal>
        </ogc:Mul>
    </ogc:Function>
</ogc:Filter>


texts: Reducing "noise" in spatially regressed points

texts:
I'm working on spatial data on a grid. I have two variables that I think to be correlated and I would like to build a model (see Fig attached).

I did some exploratory analyses following the instructions in the spdep R package and I found (as expected) spatial autocorrelation in the residuals of the classical linear regression (OSL) between these two variables.

I used the Spatial Autoregressive Model to estimate the parameters of the spatial regression.

I suppose that the noise seen in the Fig attached is partially due to spatial autocorrelation in the data.
I would like to plot my data taking in account spatial autocorrelation of the residuals, in other words I would like to reduce the "noise" across observation, so that the values fit more to the regressed line.

Would be this possible? Does it make sense my reasoning? I'm quite new to spatial statistics and I would like to learn more.


texts:
I want to create an attribute-form that displays some fields in value-related dropdowns. The goal is to limit the options in the second dropdown dependent upon the selection in the first dropdown. In my case, I have the following two fields that should be related:

List "Schwachstelle" with attributes ID (PK) and Name
List "Schwachstelle_spez" with attributes ID (PK), Name and ID_Schwachstelle

The key-attributes between the two lists are the ID of "Schwachstelle" and ID_Schwachstelle of "Schwachstelle_spez".

The filter I wrote for the value-dependent dropdown ("Schwachstelle_spez")  is the following:

"ID" = current_value( 'Schwachstelle_ID' )

However, this does still show all possible values in the "related" dropdown.

(I took https://www.youtube.com/watch?v=8shIZrwlQt4 as an example)

This is what the layer, where I specify the Dropdown-values, looks like:

And this is the layer, on which the upper layer is dependent on:

Attribute table of "Schwachstelle":

Attribute table of "Schwachstelle_spez":

If the user were to choose "Nutzungskonflikte", "Umwege" would not apperar in the following dropdown.


texts: Is there a county-level maps dataset that includes independent cities in virginia?

texts:
The following code shows the dataset I am working with:

library(maps)
CountyMapData <- map_data("county")
head(CountyMapData)

       long      lat group order  region subregion
1 -86.50517 32.34920     1     1 alabama   autauga
2 -86.53382 32.35493     1     2 alabama   autauga
3 -86.54527 32.36639     1     3 alabama   autauga
4 -86.55673 32.37785     1     4 alabama   autauga
5 -86.57966 32.38357     1     5 alabama   autauga
6 -86.59111 32.37785     1     6 alabama   autauga

There are several cities in Virginia which are independent and not part of counties.  For example, Alexandria City and Bristol City. I would like to be able to have polygons for these too.

Is there a dataset such as the above that includes independent cities but is otherwise at the county-level?


texts: Ways of merging weather data from different timezones

texts:
Does anyone know of any simple means for merging (inverse distance weighting method) hourly weather files (precip, temp, etc.) from different timezones?

I am creating weather data for HUC8s, and am using the 5 closest NCDC weather stations. But some HUC8s cross over different timezones. Has anyone figured out the easiest method for aligning weather data?

My thought is to put everything into GMT - so that 10am is 10am everywhere, and track the offset in a new column in my data. I believe this would work, but am having trouble completely wrapping my head around it:

Focusing on a watershed that cross from CST into EST - if it was raining at 10am EST, and also raining at 9am CST (raining at the timezone division, with weather stations on either side of the timezone division), it's raining at the same time, we are just marking the 'time' differently in each location. However, if I transformed my date/time into GMT, with a new column tracking the offset (maybe to convert it back later), it would be the same GMT time in both locations that the rain was occurring? I effect, it is just shifting all the weather into the same time zone, and then merging based on IDW.

Does anyone have any simple solutions I am overlooking?


texts: What is the Sentinel-2 level-1B to level-1C interpolation method?

texts:
The Sentinel-2 level-1C raster data in cartographic geometry is resampled from level-1B data in sensor geometry, using B-spline interpolation. Which order of B-spline interpolation is used?

ESA Sentinel Online states just that:

  Radiometric Interpolation

  Interpolation estimates the radiance values of the target point, knowing the radiance of neighbouring pixels. A linear algorithm is used for this operation with B-spline functions.

and I haven't been able to find the information by googling.


texts:
I created a TMS using gdal2tiles command. The tiles where generated correctly and I can see them using the leaflet and openlayers html pages also created by gdal2tiles. Now, I'm trying to open the TMS in QGIS, using gdal and the TMS mini driver.

I created the xml file describing the TMS but when I load it in QGIS, I get a blank screen with some random colored pixels on the top part of the screen. Any Idea what could be wrong? I'm using QGIS 2.18.0, 64bis from OsGeo4W (windows 10).

Bellow is my XML file:

<GDAL_WMS>
<Service name="TMS">
    <ServerUrl>./${z}/${x}/${y}.png</ServerUrl>
    <ImageFormat>image/png</ImageFormat>
</Service>
<DataWindow>
    <UpperLeftX>-20037508.34</UpperLeftX>
    <UpperLeftY>20037508.34</UpperLeftY>
    <LowerRightX>20037508.34</LowerRightX>
    <LowerRightY>-20037508.34</LowerRightY>
    <TileLevel>14</TileLevel>
    <TileCountX>1</TileCountX>
    <TileCountY>1</TileCountY>
    <YOrigin>top</YOrigin>
</DataWindow>
<Projection>EPSG:3857</Projection>
<BlockSizeX>256</BlockSizeX>
<BlockSizeY>256</BlockSizeY>
<BandsCount>4</BandsCount>
<DataType>byte</DataType>
<ZeroBlockHttpCodes>204,303,400,404,500,501</ZeroBlockHttpCodes>
<ZeroBlockOnServerException>true</ZeroBlockOnServerException>
<Cache/>
</GDAL_WMS>


texts: Generate world file from PyQgis Print composer

texts:
I've managed to save Qgis print composer maps by looking at the example below:

http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/composer.html#output-to-a-raster-image

But how can I save a world file along with the raster image from PyQgis?


texts: Server searching for wrong PostGIS extension file version?

texts:
Our database is hosted on a shared server, spatial queries are now failing unexpectedly returning the following error:

ERROR:  could not access file "$libdir/postgis-2.0": No such file or directory
********** Error **********
ERROR: could not access file "$libdir/postgis-2.0": No such file or directory

Querying the extentions via pg_available_extensions shows the extension file versions to be 2.2.5.

I'm assuming the file it should be accessing is now $libdir/postgis-2.2.5.

How can I point the server to the correct file?


texts: How to preserve dgn cad file attributes when exporting from microstation and importing into Arcmap?

texts: I'm trying to import water utility lines from Microstation into ArcMap. The water utility lines are symbolized by size in Microstation, but when I try to view the attribute table in ArcMap the attributes seem to have not transfered over, the pipe size attribute is blank. Any advice or a direction to turn? I'm poor with microstation/CAD but have a solid background in ArcGIS. Thanks!


texts: How to convert coverage (.adf) to shapefile?

texts: Is there any tool or software except ArcGIS which I can use for converting coverage (.adf) files to shapefiles?


texts: Using pandas in Python window of ArcGIS for Desktop?

texts:
How can I import pandas module in ArcGIS at figure. I want to calculate poisson and negative binomial, so I think I need import pandas or other modules.

import pandas as pd
       Runtime error
       Traceback (most recent call last):
       File "<string>", line 1, in <module>
       ImportError: No module named pandas

import numpy as np
   from scipy.stats import nbinom
   def _ll_nb2(y, X, beta, alph):
   mu = np.exp(np.dot(X, beta))
   size = 1 / alph
   prob = size / (size + mu)
   ll = nbinom.logpmf(y, size, prob)
   return ll

 Runtime error
 Traceback (most recent call last):
 File "<string>", line 2, in <module>
 File "C:\Python27\ArcGIS10.3\lib\site-packages\scipy\stats\__init__.py", line 338, in     <module>
 from .stats import *
 File "C:\Python27\ArcGIS10.3\lib\site-packages\scipy\stats\stats.py", line 184, in    <module>
 import scipy.special as special
File "C:\Python27\ArcGIS10.3\lib\site-packages\scipy\special\__init__.py", line 593, in <module>
 from ._ellip_harm import ellip_harm, ellip_harm_2, ellip_normal
 File "C:\Python27\ArcGIS10.3\lib\site-packages\scipy\special\_ellip_harm.py", line 7, in <module>
  from ._ellip_harm_2 import _ellipsoid, _ellipsoid_norm
 File "_ellip_harm_2.pyx", line 2, in init scipy.special._ellip_harm_2    (scipy\special\_ellip_harm_2.c:7480)
File "C:\Users\erdogan\Anaconda\Lib\ctypes\__init__.py", line 10, in <module>
  from _ctypes import Union, Structure, Array
 ImportError: DLL load failed: %1 not proper for Win32 application.


texts: Calculating predominant value according to area within polygons using QGIS?

texts:
I am using QGIS 2.18.1 Las Palmas, in Windows 10, and that's the issue I want to clarify:

I've got 2 polygon shapes. 1- One with river basins, and 2- other one with different levels of dessertification risk:

In the attribute tables for each one, I've got names for each river basin in the first map, and in the second one, I've got 5 different levels (differenced by colors in the picture)

So, what I need is not only cut or intersect both of them. What I need is to know which level of dessertification risk is more predominant, within each river basin. In other words, I need a final shape similar as the river basins shape, but with an additional column saying the level of dessertification risk which occupies the biggest % of surface within each basin.

For example: I will zoom into one of the basins, setting the river basins without filling, so you can see the dessertification risk layer behind:

We see different colors within the polygon. Initially, it seems that the most predominant color is purple, or could be the light green.

What I need is to determine which one covers the biggest % of surface within this polygon, and set it as a new value for the river basins polygons.

I've tried different things but I don't find the solution.


texts: QGIS & GeoMedia .mdb .gws compatibility?

texts:
i cannot open Geomedia file directly from QGIS 1.8.0! Gdal driver >= 1.9.0

reference is made to below link

qgis-connecting-geomedia-mdb-warehouse

anyone have any experience in using GeoMedia file in QGIS???


texts:
I have a large number of Raster Aerial images which I have been building VRT files for.  I'm also building overview for them with the following code

gdaladdo -r average -ro Y:/vrt/sussex_2016_aerial.vrt 2 4 8 16 32 64

The orginal images are GEOTiffs and I would like to create compressed overviews but I'm struggling to understand the code required, which I need to add to the above line.  I would be grateful if someone can point me in the right direction?


texts: Adding WMTS service as Basemap in Mapbox?

texts:
I am trying out Mapbox and I do not see any ways to add my WMTS service as a basemap in Style in Mapbox.

Instead of using the Mapbox basemap, I will like to use a WMTS service as the basemap. There is not options to add WMTS in Tilesets. Styles only allow me to add a file or Tilesets/Datasets. I am trying to figure out how to add a WMTS to Style so that I can use it as the basemap before publishing it. I am also very new to Mapbox, trying to figure out how to use.


texts: QGIS: How to store the GPS coordinates in the attribute table?

texts:
I got vector points on a map. I already have the elevation information in the table for each vector point and the point-id. Now I want to get the GPS coordinates for each point stored in this attribute table. How can I do that?

Overall my aim is to get to know the gradient on this line which I divided into those points.


texts: Any way to export images from ArcGIS FeatureServer?

texts:
I'm trying to export images from an ArcGIS server.
Exporting MapServer layer images is just okay, but FeatureServer doesn't support exporting to images. Although it can query for features and the layerInfo contains drawingInfo.
Also the server can display the features in the map.
Current idea is to use selenium, but moving map will be tough one.
Anyone knows about a tool that would do the job?


texts:
In FME I have the following in the AttributeManager transformer to replace a value in a column

@ReplaceRegEx(@Value(Street_Address),"^([UNITunit]*.?\d+[\/-]?\d*\w).*",\1)

This should work as in the example in

https://regex101.com/r/CwKOSx/3

But instead it's doing the following

Truncated value `1 Keith Williams Drive CARDWELL QLD 4849' for column `Street_No' of MapInfo TAB file `\\sipv-gis01\Projects\Jobs\2019.0225_Property Update Automation\Working\FME_Scratch\Property_w2.tab' to the column's defined width of 26
Truncated value `Foxtail Avenue CARDWELL QLD 4849' for column `Street_No' of MapInfo TAB file `\\sipv-gis01\Projects\Jobs\2019.0225_Property Update Automation\Working\FME_Scratch\Property_w2.tab' to the column's defined width of 26
Truncated value `Unit 1/15 Keith Williams Drive CARDWELL QLD 4849' for column `Street_No' of MapInfo TAB file `\\sipv-gis01\Projects\Jobs\2019.0225_Property Update Automation\Working\FME_Scratch\Property_w2.tab' to the column's defined width of 26

Where I need the following for these

1
[no value]
Unit 1/15

I am not sure if this is a regex or fme question but how do I get it to show just the first match.

UPDATE

Using String Searcher with regex

^(UNIT\s)?\d+[-/]?\d*\w*

I get mostly correct output except

The model is as below


texts: How to retrieve answer from sendfile_osm_oauth_protector

texts:
I made a request by using the enter link description heresendfile_osm_oauth_protector
from subprocess import check_output

command = "curl -b 'text from cookie_output_file' https://osm-
    internal.download.geofabrik.de/asia/philippines-latest-internal.osm.pbf"
out = check_output(command)

how do I find the link to download in this output? it is like a complex text.


texts:
I want to access a feature class inside a geodatabase (gdb) from inside my shinyapp. I can access the fc class once I provide the gdb path:
 gdb_data = arcgisbinding::arc.open(gdb)

but how do I get the information of the path to my gdb as in input in my shiny app? fileInput doesn't work as a .gdb acts like a folder. I tried the package shinyFiles but that did not lead to success, as I could not navigate to any drive (I remained in Home, could not navigate to any depth in the files). I used this https://community.rstudio.com/t/shiny-directory-input/29160/2 as a template.
UI:
ui <- fluidPage(
 mainPanel(
            shinyDirButton("dir", "Input gdb", "Input gdb"),
            verbatimTextOutput("dir", placeholder = TRUE)
          )
        )

Server
server <- function(input, output, session) {

  shinyDirChoose(
    input,
    'dir',
    roots = c(home = '~'),
    filetypes = c('', 'txt', 'bigWig', "tsv", "csv", "bw")
  )

  global <- reactiveValues(datapath = getwd())

  dir <- reactive(input$dir)

  output$dir <- renderText({
    global$datapath
  })

  observeEvent(ignoreNULL = TRUE,
               eventExpr = {
                 input$dir
               },
               handlerExpr = {
                 if (!"path" %in% names(dir())) return()
                 home <- normalizePath("~")
                 global$datapath <-
                   file.path(home, paste(unlist(dir()$path[-1]), collapse = .Platform$file.sep))
               })
  }


texts:
I am trying to get the coordinates of each node of the polygon I create. I would like the attribute table to display the name of the polygon (ex. address) and then the coordinates of the nodes (ex. nodeA= Ax and Ay, NodeB Bx and By, NodeC Cx and Cy, NodeD Dx and Dy) is there a easy way to create all the needed polygons and then have the coordiends of the nodes displayed in the attribute table?

any help would be great!!

Thanks


texts: Delete an SSL certificate for ArcGIS Server

texts:
I added an SSL certificate to ArcGIS Server, and named it *.mydomain.net.

ArcGIS Server does not seem to like the asterisk though; when trying to open the REST interface (https://ags.mydomain.net/arcgis/admin/machines/AGS.MYDOMAIN.NET/sslcertificates/*.mydomain.net) the browser just shows Bad Request (Firefox) or Cannot find webpage (IE). Replacing * with %2A does not work either.

Obviously, I cannot delete this certificate through the web interface now.

Is there a way to do this manually? Maybe there's a directory, or registry enter or whatever where AGS stores certificates?

(Note: I already added the same certificate with a valid name too, and everything works, it's just that I want to clean up the mess)


texts: Using "Elseif" "OR" "And" statement to create label from two fields

texts:
I'm an admittedly novice when it comes to scripting but I have been able to create a basic 'else if' with 'OR' expression to create a label for a fuse feature class I have constructed. I am now trying to combine an additional field into this statement but have been at a lose on how to get there.

In the code shown below I am labeling fuses by assigning a number (1, 2,3) based on its phasing code (A, AB, ABC); respectively. I now also want to assign an "I"(installed) or a "P"(proposed), preferably in bold, that corresponds with my "Construction Status" field.

Function FindLabel ( [phasingcode] )
Dim a
If [phasingcode] = "A" OR [phasingcode] = "B"  OR  [phasingcode] = "C" then
    a = "(1)"
Elseif [phasingcode] = "AB" OR [phasingcode] = "BA" OR  [phasingcode] = "AC" OR [phasingcode] = "CA" OR [phasingcode] = "BC" OR  [phasingcode] = "CB" then
  a = "(2)"
Else
 a = "(3)"
end if
FindLabel = a
End Function


texts: PostGIS  operator does not return true distance

texts:
I want to find the closest polygon to a point. My idea was to use the <-> operator for it. According to the docs it says:

  For PostgreSQL below 9.5 only gives centroid distance of bounding
  boxes and for PostgreSQL 9.5+, does true KNN distance search giving
  true distance between geometries, and distance sphere for geographies.

I use PostGIS 2.2 and PostgreSQL 9.5.7. So, the <-> operator  should use the "true KNN distance", but it does not: it still orders by the centroid_distance.

SELECT id,
       ST_Distance(my_point.geom,
                   my_polygons.geom) AS real_distance,
       ST_Distance(my_point.geom,
                   ST_Centroid(my_polygons.geom)) AS centroid_distance

FROM (
    SELECT 'blue' AS id, ST_Polygon(ST_GeomFromText('LINESTRING(0 3, 2 3,  2 4, 0 4 ,0 3)'),4326) AS geom
    UNION
    SELECT 'green' AS id,  ST_Polygon(ST_GeomFromText('LINESTRING(2 0, 8 0 , 8 2, 2 2 ,2 0)'),4326) AS geom

) AS my_polygons,

(
    SELECT ST_SetSRID(ST_Point(1, 1),4326) AS geom
) AS my_point

ORDER BY my_point.geom <-> my_polygons.geom

Result:

|       | real_distance | centroid_distance |
---------------------------------------------
| blue  |       2       |        2.5        |
| green |       1       |        4          |

I expected that the <-> operator orders the polygons according to the real_distance and not the centroid_distance.

Am I doing something wrong?


texts: Is it possible to remove (Actions) and (Derived) from the QGIS Identify dialog?

texts:
When identifying a layer in QGIS 2.0.1 Dufour, the Identify Results dialog shows (Actions) and (Derived) keywords:

Is it possible to remove these?


texts: Where can I download NDVI or landsat images of south china that isn't in an .img file format

texts: Everything i seem to find is an hdf file and i can't open it in erdas. so is there a way to convert it to open in erdas or do i need to find somewhere else, earth explorer and glovis are the two places i've tried so far.


texts: Showing positive and negative values in legend of ArcMap?

texts:
I want to draw this Legend by ArcGIS 10.2. This legend starts with a big blue circle and decreases in size till 0, then in 6 its bigger and red. I know its related to symbology and multiple attributes but which part could do it for me I don't know.

I can't show both red and blue circles together.
It starts with a big size and goes to smallest. I can set it by hand but its not completely same in both negative circles and positives circle because sometimes the intervals are not the same.

.imgur.com/Rmckn.png


texts: In pyqgis, how to add a feature to a layer without opening any edit form?

texts:
I want to add a feature to a polygon vector layer in QGIS.

I have:

myLayer.featureAdded.connect(doSomething)

And when I click on a button:

iface.setActiveLayer(myLayer)
myLayer.startEditing()
iface.actionAddFeature().trigger()

So when the button is clicked I can draw a polygon on the canvas. When I end the drawing the usual edit form appears.

Is there a way to prevent it from opening or to close it immediately? (I want to edit the attributes programmatically in doSomething).


texts: How to configure QGIS (SAGA, GRASS algorithms) in Ubuntu

texts:
I have installed GRASS GIS following this install manual.

I have installed SAGA GIS following this install manual.

and finally I have installed QGIS 2.18.1  following the official details for install in Ubuntu Debian/Ubuntu.

The problem is I can't use algorithms from GRASS, SAGA and TAUDEM in the QGIS
interface.

I have toolbox processing plugin and I have configure the options providers for all software (SAGA, GRASS) and not working.

Only GDAL tools are working in the plugins and QGIS algorithms. Standalone out of QGIS, SAGA and GRASS GIS work fine.

How to configure QGIS to work like work in Windows?

System details QGIS 2.18.1 Ubuntu 16.04.

GRASS ERRORS:

SAGA GIS ERROR

SAGA GIS ERROR

TOOLBOX OPTIONS SETTINGS:

Finally I think so this is solution for my another question Why can't I use processing algorithms in a standalone PyQGIS script?


texts:
i do some experiment place two sensor node between them 4 meter as distance.

then i get these GPS information:

node1:  Latitude: 32, 33.255  , Longitude: 36, 0.6632

node2:  Latitude: 32, 33.228  , Longitude: 36, 0.6588

then to try to find the distance between them i use this Link, but i surprised that the distance calculated is 50 meters !!!!!

how can i go over this problem !


texts: Migrating ArcPy scripts from ArcGIS 9.3 to 10?

texts:
What changes are needed to migrate scripts from the ArcGis 9.x environment/gp-object to the ArcGis 10/arcpy environment?

Can anyone provide a good overview what the changes are?

Can someone share pitfalls etc from their own experience?


texts: How can I tell who owns a particular piece of property on a map?

texts:
I already saw mentions of maponics and reportall but they both seem to be address specific one by one sort of tools.

I basically would like to know who owns various pieces of land...ideally if I can see that data with a simple map overlay. I read there was a way to do this in google maps but I can't get that sorted.

It doesn't have to be free, it could be an online service I can use myself or it could be someone that I pay it doesn't really matter. I just want to be able to send a link or screenshot to a particular section of google maps and get the owner/availability/price so that I can make contact and purchase the land.


texts:
I'm using this code as a python script tool in ArcGIS 10.1:

import arcpy, os
arcpy.env.workspace = r'W:\S&P\s&p techs\Emily\Errors.gdb'

#Looping through dissolved feature classes, adding 'Name' field and writing
#feature class name in the added field.
fcs = arcpy.ListFeatureClasses()

for fc in fcs:
    arcpy.AddField_management(fc, "Name", "TEXT", field_length = 50)
    with arcpy.da.UpdateCursor(fc, "Name") as cursor:
      for row in cursor:
        row[0] = fc
        cursor.updateRow(row)

All it's supposed to do is add a field to some feature classes in a database. I double checked that the path is right, so I don't think that's the problem. This code has also worked before, but for some reason it's not working anymore. I get an error that says this when I run it:


texts: Using ST_* functions with different origin results in: Function not found

texts:
I've built a query which uses ST_asGeoJSON(); and If I run it on my server through pgadmin with credentials of UserX it runs like a charm.

But when connecting through an webapp with other credentialsof UserY it returns me the following error:

  SQLSTATE[42883]: Undefined function: 7 ERROR: function st_asgeojson(public.geometry) does not exist
  LINE 17: ,ST_AsGeoJSON(vbo.geopunt) AS vbo_geopunt
  ^
  HINT: No function matches the given name and argument types. You might need to add explicit type casts.

So I tried to explicitly cast to ::point but then It still does not work.

If I connect my webapp with the working credentials of UserX it also produces the same error.

Any idea what I'm missing?


texts: How do I create a shapefile of self-defined zones, using clusters of zip codes (ArcGIS)?

texts:
I have 6 zones in the United States that I defined by clustering all US zip codes. I now need a shapefile that depicts each of these zones in a different color. Here is an illustration of what I would like the map to look like:

What I am working with: A dataset containing every US zip code, their approx Long and Lat, and their zone designation. I have access to ArcGis 10. My knowledge of the software is limited, but I am a quick learner!

Ideally, what I would end up with is a single shapefile with each of the zones in a different color. The point of this is to use that shapefile in an app, which will project points onto a map of the US. That part I can do, I just need to user to be able to see which zone the point is in.


texts: In Leaflet I want to add a legend when I turn on an overlay, then when I turn off that overlay I want the legend to disappear and when I turn on a second overlay I want a second legend to appear. I have 2 png legends that I am using. I'm using 'overlayadd' & 'overlayremove' but can't get it working


texts: Is there a way to prevent a WMX user from adding an embedded attachment?

texts:
I would like to allow Linked and URL attachments but prevent embedded attachments.

The ManageAttachments    privilege seems to be all or nothing.  Does anyone else require this particular configuration? How did you accomplish it?


texts:
I am trying to get Solar Radiation data from WorldClim version 2.0 in R,
However, I am using the library raster in R  but it seems to not be possible to download this data. (since there is only functions for bioclimatic data, Temperature, altitude, and precipitation)

Does anyone know how to download this data from the website using R?

I have specific coordinates. And for my bioclimatic data I have used this code, which worked:

library(raster)
library(sp)

r <- getData("worldclim",var="bio",res=10)
r <- r[[c(1,12)]]

names(r) <- c("Temp","Prec")

lats <- c(48.45435)
lons <- c(-2.04734)

coords <- data.frame(x=lons,y=lats)
points <- SpatialPoints(coords, proj4string = r@crs)
values <- extract(r,points)
df <- cbind.data.frame(coordinates(points),values)
df


texts: Splitting a bounding box into n amount of bounding boxes

texts:
Firstly, I've never worked in the GIS area before nor do I have much knowledge of the mathematics behind some of the stuff I am trying to implement. Perhaps i'm getting in over my head.

My Situation

I am creating a geospatial map viewer in Qt. So far I have a MapTile class which stores its center location (lat and long) and the width of the box in kilometers, which may or may not be a good idea. From this information I can construct a bounding box from the center point and the width by using the algorithm found here.

The MapTiles are all square, no irregular shapes or rectangles. Is this a bad implementation?

My Problem

Now I am stuck with implementing my MapTileMatrix class which stores a multidimensional array of MapTiles[rows][cols]. I figured I could tell the MapTileMatrix what its bounding box is, and then have it generate the MapTiles position and width from this data. I haven't been able to find any algorithm online or mention of how to split a bounding box down into n number of sub-boundingboxes.

Is there a method of dividing a bounding box into say 4 equal parts? What about 9, 16 or 25? If this is not possible perhaps I should reconsider my design.

On a side note:

Are there any guides / tutorials to implementing a map viewer application? I have read over the Web Map Tile Service and am using this as a guide. This has been my introduction to GIS programming and I have only been working on it for a few days so have little to no knowledge of the domain.


texts: How do I draw transparent polygons with ArcObjects?

texts: How do I draw transparent polygons in ArcObjects? I mean is it possible to set alpha channel?


texts: Enabling C++ standards when compiling GEOS using cmake/ccmake?

texts:
I read from RFC # 7 of GEOS that the old autoconf based build system may be removed in the future, in favor of using cmake.

I sometimes need to customize geos and add functions that uses the newest C++ features (e.g. c++17). I know how to enable such features on the command line with autoconf, but

How to enable e.g., c++17 in the cmake system?

-- EDIT --

I have looked in SO, and the answers there suggest adding things to CMAKE files. However, I use ccmake, and do not want to edit any cmake files.

Also, to clarify, one reason I ask here instead of SO, is that GEOS is a common GIS package which may or may not be intended for C++17. The adoption of C+11 happened only recently.

That's why I think it's better to ask about future GEOS options here, rather than SO.


texts: How to show a specific geometry on the map and hide the others?

texts: I have a FeatureLayer, i used the method IFeatureLayer.Search() and the interface IQueryFilter to get a specific Feature, and i want to show only this Feature on the Map, how can i do that ?


texts:
I am having trouble with an WMS layer that calls the Ordnance Survey/Digimap OpenStream data, but however I call the layer it times out without ever displaying a map.

I have registered with Digimap and have my own API key. I have successfully created a WMS connection that allows me to call an OpenStream layer, but once I have clicked 'Add' to place the layer in the map, nothing appears until I receive a notice saying 'connection timed out'. This problem occurs whatever parameters I set for the OpenStream entry (username and password, CRS options, tile size... ). I can successfully create and load other WMS connections and add a layer, such as the following: http://ows.terrestris.de/osm/service

I am not seeing firewall reports that suggest that my firewall is blocking incoming or outgoing connections, and QGIS is an allowed application (allow TCP/IP and UDP in and out).

System is QGIS 2.14.15 ltr on Windows 7 Ult 64 bit.

Is this a temporary fault with Edina (could others test, please), or have I got something basic wrong?


texts: Expression Builder in QGIS in editing Attribute Table

texts:
I am working with a simple attribute table that has 4 columns;
'Xcoord'
'Ycoord'
'id'
'id2'

'id2' is a new field that i added to the attribute table.
'id' is a field showing land use classes with values ranging from 0-210.

I wish to make 'id2' be filled up with new values ranging from 1-6, based on their value in the 'id' field.

For example, I wish to write an expression for;

when 'id' are values 10,20,30 and 40, i want them to have a value of 1 in 'id2'
and 'id' values of 50,60,70 as 'id2' value of 2
and so on and so forth till an 'id2' value of 6.

Basically like a hierarchical aggregation of land cover classes.

I am trying to self-learn the QGIS expression builder but I am encountering many difficulties trying to code the expression properly.

Here is an image of the


texts: Create Area or Polygon based upon a Series of points

texts:
i have a large list of points that i would like to merge to create an Polygon based upon an attribute.
Essentially an area that encompasses all points.

x   |   y   |   attr
____________________
1   |   1   |   One
2   |   2   |   One
1   |   2   |   One
2   |   1   |   One
1.5 |   1.5 |   One

in the above example i would expect to see a box around the fifth point with an attribute of One

my data is held in SQL Server and i could either export this into qGIS to resolve or solve within SQL.

is there any way of re-solving this, or does the issue of generating a Link between the 5th and 1st point prevent this?


texts: How do I link two-digit Dutch postcodes to latitude and longitude for mapping data to a map of the Netherlands?

texts:
I want to show data on a map of the Netherlands. To ensure my Dutch participants' privacy, I only asked them to specify the first two numbers of their postcode.

I found latitude and longitude for all Dutch postcodes at http://postcodedata.nl, but these are for the full postcode, and I don't know how to map these more 'spatially specific' latitudes and longitudes to my less fine-grained data.


texts: Shortest path creation between two lines

texts:
I've a raster of temperature of a river, I'd like to extract values horizontally, so from bank to bank, shortest path. so basically extract pixel values along the red lines in the schema I did. What is the easiest way to generate these red lines ?

I've Arcgis 10.4, Qgis, matlab. Tried a lot of stuff already but nothing convincing, I extracted the banks and transformed the lines to points, numbered them and planned to joint #1 from left bank to #1 of right bank, but that doesn't creates shortest path. Then I wanted to create perpendicular lines to every point, taking a line between precedent and successive point as a reference, but my coding skills in matlab are insufficient.

I want use these extracted data to calculate lateral STD and different stats.


texts: Creating custom CRS from existing CRS EPSG:22780?

texts:
I am currently trying to figure out how to create a custom CRS in QGIS based on Deir ez Zor / Levant Stereographic (EPSG:22780). For this I am working with a map which shows both the EPSG:22780 coordinates and the coordinates of the custom CRS.
I have to work with the following criteria:

The custom CRS should be specific to a rather small and limited area (see map).
The negative coordinates of EPSG:22780 should be replaced by positive and "simple" coordinates (three-digits in pre-decimal position; see map).

The following topo-map, which I had to redact as it is not published yet, shows how the old and new coordinates should correlate:

I tried to figure out how to change the proj.4 parameters of EPSG:22780*, but as I am doing this for the first time and without any guidance I am kind of lost atm.

Here they are:

+proj=sterea +lat_0=34.2 +lon_0=39.15 +k=0.9995341 +x_0=0 +y_0=0 a=6378249.2 +b=6356515 +towgs84=-190.421,8.532,238.69,0,0,0,0 +units=m +no_defs

*(I got these from the Proj.4 library in QGIS and noticed, that epsg.io provides other values for +towgs84 (+towgs84=-187.5,14.1,237.6,0,0,0,0) - in case this is somehow relevant. The other values are the same.)