texts:
I had created an shaplefile on QGIS. That shapefile came from a georeferenced raster.

I have created that with the CRS EPSG: 3857 - WGS84/PseudoMercator, than i converted to EPSG:4326, WGS 84.

I have open this new project and set the WGS 84 coordinate system to project. After add this new vector to QGis it is not georeferenced.

Can you help me with this issue?


texts:
I need to develop a shopping mall GIS. I have no idea to develop it. The shopping mall GIS should function like GPS (guides user to desire shops). Do I need to set up some GIS device at each corner of shopping mall, in order to keep track location of user? What are those devices?

Does anyone involve in GIS development before? I need your help.

Through  an extensive research, WiFi Positioning System (WPS) is the only feasible solution.

One of the features of my system:
- Example, you connect to a shopping mall WiFi using your mobile phone (iPhone, Windows Phone 7, Android Phone, etc). It will redirect to the main page. You can find desired location by clicking "search" button. It will guide you the way to the specific shop (like google map). A red dot in the map represents your current location. The red dot will move if you are moving. WIFI Positioning System is the way to achieve this. How to setup WIFI Positioning System? What kind of device I need? Does anyone have experience and knowledge in this kind of development? I am totally new in this kind of development. Please share your opinion and guidance.


texts: MacOS QGIS3.8.2, Opening QGIS prompts '....save the current project.'

texts:
When opening QGIS Zanzibar the following prompt message occurs. If the splash is marked on, in the settings, the prompt remains hidden but awaits for user prompt at startup.

Help to get rid of message.

This was not present in QGIS 3.4.


texts:
Is there an aplication to add tree menu for the map showed with openlayers?

I mean show a map on my site exactly how is showing in google earth for example, my layers organized in folder/subfolder, with checkboxes.


texts: can anyone suggest me for changing mapset from current location using command prompt.I have compiled the grass gis 6.4.4 and when ever i want to use the commands every time i am opening gui and setting the mapset and location. without starting grass gui, i want to run my module. is it possible? And i have one more query regarding the mapset.my module is every time it is looking for same mapset and i have to load that only, but i want to test it for other data of different region. even if i change mapset and location in gui , libgrassgis-dll missing it is showing and not working.without grass gui i want to run it.may i know where the problem is and is?


texts: How can we export shapefile to SVG format using PyQGIS?

texts:
I am doing a project where I have to export file from shapefile to SVG format using PyQGIS.

I have read some articles regarding SimpleSvg plugin but I don't know how can we achieve the same task using PyQGIS.


texts: Intranet GIS app with JS web app builder?

texts: Is there a work-around for building an intranet GIS app with ESRI JS web app builder (without sharing data to the public)?  Does it mean that I have to create a web map in ArcGIS com or is there an alternative?


texts: Saving (converting) raster objects in R?

texts: I want to save my raster plot that produced with R to ArcGIS raster (.img, .tif and ASCII) file.


texts: How to connect polygons to make a circle?

texts:
I created polygons from a raster on which I ran an edge detection filter. I turned the edges into a polygon layer (multipart splitted into single parts). The intention was to have circular polygons at the end, but now there many single polygons forming one circle but are not touching each other, so I cannot connect these. The single yellow polygon in the figure below is an example.

Can anybody tell me how to connect all polygons in order to form a circle?
The polygon is available as a download here. I work in SAGA, GRASS, or QGIS, or in R.


texts:
I am trying to show layer from Geoserver as WFS with OpenLayers, but when I open http://localhost:8080/geoserver/www/Geoinformatika2/index.html, it gives me only OSM map as a base, but there is no polygon that is defined in layer

var osnova = new ol.layer.Tile({
          title: 'OSM',
          type: 'base',
          visible: true,
          source: new ol.source.OSM({
            layer: 'OSM'
        })});
var logori = new ol.layer.Vector({
          url: 'http://localhost:8080/geoserver/mydb/ows?service=WFS&version=2.0.0&request=GetFeature&typeName=mydb:geoserver&maxFeatures&outputFormat=application/json',

format: new ol.format.GeoJSON()
        })
});

var map = new ol.Map({
  target: 'map',
  controls: ol.control.defaults().extend([
    new ol.control.ScaleLine(),
    new ol.control.ZoomSlider()
  ]),
  layers: [
    new ol.layer.Group({
      'title':'Osnova',
      layers:[ osnova ]}),

    new ol.layer.Group({
      title:'Adrse',
      layers:[logori]})
    ],
  view: new ol.View({
    center: ol.proj.fromLonLat([17.86249816417694, 44.60232779122055]),
    zoom: 8
  })
});

as the response I get

    {"type":"FeatureCollection","totalFeatures":1,"features":[{"type":"Feature","id":"geoserver.2","geometry"
:{"type":"Polygon","coordinates":[[[6414273.98269,5015670.1361],[6414000,5015000.1361],[6416273.98269
,5017670.1361],[6414273.98269,5015670.1361]]]},"geometry_name":"koordinate","properties":{"stratisteid"
:2,"naziv":"s","tip":"logor","zrtve_broj":6000,"tipgeometrije":"Polygon"}}],"crs":{"type":"name"
,"properties":{"name":"urn:ogc:def:crs:EPSG::31276"}}}

If I try to get it like this

var osnova = new ol.layer.Tile({
          title: 'OSM',
          type: 'base',
          visible: true,
          source: new ol.source.OSM({
            layer: 'OSM'
        })});
var logori = new ol.layer.Vector({
        source: new ol.source.Vector({
            format:new ol.format.WFS({
                version:"1.1.0",
                url: "http://localhost:8080/geoserver/mydb/wfs",
                featurePrefix: "mydb",
                featureType: "mydb:geoserver",
                featureNS: "http://localhost:8080/geoserver/mydb",
                geometryName: "koordinate",
                srsName: "EPSG:4326"
            })
        })
});

var map = new ol.Map({
  target: 'map',
  controls: ol.control.defaults().extend([
    new ol.control.ScaleLine(),
    new ol.control.ZoomSlider()
  ]),
  layers: [
    new ol.layer.Group({
      'title':'Osnova',
      layers:[ osnova ]}),

    new ol.layer.Group({
      title:'Adrse',
      layers:[logori]})
    ],
  view: new ol.View({
    center: ol.proj.fromLonLat([17.86249816417694, 44.60232779122055]),
    zoom: 10
  })
});

it also gives me the only OSM map, but I can't see GET request for WFS layer. I can't figure what is wrong
URI of workspace is http://localhost:8080/geoserver/mydb
layer name is mydb:geoserver
if I try to open layer from the link: http://demo.opengeo.org/geoserver/topp/ows?service=wfs&version=2.0.0&request=GetFeature&typeNames=topp:states&maxFeatures&outputFormat=application/json
it shows it with no problem, so it looks like openlayers is not getting layer from geoserver


texts: Conversion of a GeoJSON into EWKB format?

texts:
I see libraries such as shape from shapely.geometry that support conversion to WKB and WKT binaries. What I am unable to find is guidance on converting to EWKB specifically.

I suspect this is because libraries that support WKB conversion also support EWKB? How do I know if I am creating an EWKB vs WKB? Are all WKB also EWKB (but not the other way around)?

If EWKB is an extension of WKB, how would I check that, say, Shapely, supports EWKB format?

I am trying to understand if EWKB is standard or something that was just built for PostGIS (where I have seen it used) and how I can make sure that data I am producing for a function that expects EWKB can be correctly formatted.

Because of @shongololo's comment, I wanted to add a snippet of what I am going presently to generate the WKB of the shape from GeoJSON.

from shapely.geometry import shape
import binascii
# convert geom to wkb format
geom_str = json.dumps(route['geometry'])
geom_geo = geojson.loads(geom_str)
geom_shp = shape(geom_geo)
geom_wkb = geom_shp.to_wkb()
geom_cln = binascii.hexlify(geom_wkb).decode('ascii')


texts: Finding River Lots in Alberta, Saskatchewan and Manitoba using QGIS?

texts:
I have to border farmlands in Canada (Alberta, Saskatchewan and Manitoba). Usually they give input in the form of DLS like SE-9-56-22-W4.So I use https://locator.scadalink.com/dls2latlng.html to find location and border farmlands.(like in Screenshot below)

Each Section has 640 Acres,Half 320 Acres and Quarter 160 Acres.

Currently I am facing issues, I recieved input as RiverLots like RL-8-55-22-W4. I don't know how to find these RiverLot location.

I have tried Scadalink Locator and LegalLandConverter but no use.

How do I find and border RiverLots using QGIS?


texts: How to access geoserver web services on a backend server through openlayers web application hosted on another server via IIS?

texts:
I am struggling with how to access geoserver web services in my openlayers application. I setup a test machine that hosted both geoserver and the openlayers app as well as a postgres database to publish data from. I was able to just add the services as layers to my application by specifying crossorigin: anonymous. Now I need to move to a different environment where geoserver resides on a separate server, and the postgres database also resides on its own server, and these two servers are not public facing, but are on the same network as a public facing server which will host the openlayers app. When I try to access the application internally I get a CORS error and the services won't load; externally the server name included in the service url just can't be resolved. The web server is running IIS on port 80, and geoserver is running on port 8085 on its own backend machine. What do I need to do be able to see my services in my web application in this 3 tier environment?

UPDATE:

I tested the solution below, and it works (whether I set the proxy based url or not) when i am logged into a virtual machine on our network. However when I try to login to the application via a browser on a machine outside of the network, which I can do by supplying my domain credentials, the services still don't load. I get this error for each of my service URLs:  Failed to load resource:net::ERR_NAME_NOT_RESOLVED.


texts: Resize Polygon just by outer edge using OpenLayers 2?

texts:
I am writing a program in openlayers2 of polygon resizing. So resizing function take two parameters 1.Scale and 2nd Origin. So my problem is I don't want to resize by taking origin. I want to change size by other methods, like By taking outer edge etc. So can someone guide any other method to resize polygon ?

function resizeFeatures(scale) {
        polygonFeature.geometry.resize(scale, origin);
        vectorLayer.redraw();
    }


texts: How do you change the Window Icons in QGIS 1.7

texts:
I am a beginner to QGIS 1.7, and I want to change the icon files in QGIS, (wherever we see the 'Q' icon I have to change that). I tried and I got success in the case of the mainWindow, and also somewhere it was reflecting but the problem was with menuitem named 'Layer'.

When ever I trigger 'add postgis layer' a widget opens but the icon 'containing Q' was still there...

So could you suggest a method through which I can change all the window Icons in the  QGIS with the icons that I want?


texts: PyQGIS toggle group visibility (and recursively subsub...group visibility)

texts:
I would like to set visibility of a layer group and all of its children (and grand-children and so on...) using PyQGIS (QGIS 2.14.1).

Assume the groups name is g0 with subgroups g0[1...n] and subsub...groups g0[1...n][1...n][...].

QgsProject.instance().layerTreeRoot().findGroup('g0').setVisible(False)

sets the group and ALL what comes beneath invisible, so far so good:

But vice versa

QgsProject.instance().layerTreeRoot().findGroup('g0').setVisible(True)

sets the group visible and leaves all its children untouched, wich results in nothing visible at all:

QgsProject.instance().layerTreeRoot().findGroup('g0').updateChildVisibility()

only touches the direct children and not recursively all the descendants:

I wonder why manually setting the root groups (g0) visibility does result in recursively setting the visibility of all beneath it at all.

Question: Do I have to implement some recursive function to set the visibility of a group and its whole big family, or is there a more convenient way to accomplish this? Or is here something buggy?


texts: Making automatic updated date and time field in QGIS?

texts:
I want to make a date and time updated field as I edit an attribute field in this format yyyy-mm-dd 00:00:00 (= 2016-05-08 11:04:00) - my local date/time. The data is in a QGIS 2.14 shapefile and vector lines as screenshot below.
When a feature as NULL in the 'name' field, I want to enter a name for the road and the 'mod' field must receive de modified date/time like '2016-05-08 15:16:00'.
I already tried some options and failed.
I need help to start from beginning.

I create a new field, date/time type, but stays with NULL value after I edited, in the "name" field, or the feature itself.
I used too, the expression now(), but stays with NULL value after I edited.

Note: the date seen there, 2016-04-14 is the first value created with "Field Calculator".

I am using QGIS 2.14.


texts:
I'm currently stuck in a project because whenever i try to open GRASS console (shell) it doesn't appear anything, and i cant find the reason why.

Im using QGIS Desktop 3.4.3 with GRASS 7.4.2

What i can do?


texts: Error -Collection query aborted after accumulating over 5000 elements

texts:
I have created a code that suppose to take sentinel images and calculate the NDVI for a speficif polygon I have. I need it to run it for 2 years, but I have realized that for some reason the dataset is way too big and I get wrong numbers, for example, if I ask for the size of the dataset between 01/01/2018 to 03/01/2018 I get that the size of the dataset is 5099 images which doesn't make any sense.

this is the code I have up to now:

/**
 * Function to mask clouds using the Sentinel-2 QA band
 * @param {ee.Image} image Sentinel-2 image
 * @return {ee.Image} cloud masked Sentinel-2 image
 */
function maskS2clouds(image) {
  var qa = image.select('QA60');

  // Bits 10 and 11 are clouds and cirrus, respectively.
  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;

  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

  return image.updateMask(mask).divide(10000);
}

// Map the function over one year of data and take the median.
// Load Sentinel-2 TOA reflectance data.
var dataset = ee.ImageCollection('COPERNICUS/S2')
                  .filterDate('2018-01-01', '2018-01-03')
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
                  .select('B2','B3','B4','B8','QA60')
                  .map(maskS2clouds);

var rgbVis = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};

var clippedCol=dataset.map(function(im){
   return im.clip(geometry);
});

//test if clipping the image collection worked
Map.centerObject(geometry);
Map.addLayer(clippedCol.median(), rgbVis, 'RGB');

// Get the number of images.
var count = dataset.size();
print('Count: ',count);
print(clippedCol);//here I get the error messege "collection query aborted after accumulation over 5000 elements
print(dataset,'dataset');//the same error here

//function to calculate NDVI
var addNDVI = function(image) {
  var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI');
  return image.addBands(ndvi);
};

//NDVI to the ckipped image collection
var withNDVI = clippedCol.map(addNDVI);

// Test the addNDVI function on a single image.
var ndvi1 = withNDVI.select('NDVI').mean();

var colorizedVis = {
  min: 0.0,
  max: 1.0,
  palette: [
    'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
    '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
    '012E01', '011D01', '011301'
  ],
};

Map.addLayer(ndvi1,colorizedVis,'test');
Map.centerObject(geometry);

my goal is in the end to see the NDVI mean value for every image.


texts: Can I update the service on ArcGIS for Server programmatically?

texts:
There is Windows Service. This service adds user raster data to the Raster Dataset in sde database.

The Raster Dataset registred on the ArcGis Server. This Raster Dataset also was publishing on ArcGis Server. When the user data has been added to this Raster Dataset, the data is not updated on ArcGis Service.

How update service on the ArcGis Server from my Windows Service? (I know IP address, password, and service name)


texts: Pyqgis: How to take output from Qgs.ZonalStatistics and join by ID with point layer

texts:
What I'm trying to do:

Extract median by a single polygon from multiple GeoTIFF rasters stored in a folder (Solved)
Add each median value to a point shapefile (Match/join by: id column for each point in shape = raster filename)

Basically the project has three sources:

One folder containing 558 GeoTIFF raster's which are viewshed results from different observer points. All raster's have the same extent and are named 1.tiff, 2.tiff, 3.tiff .... 558.tiff
One vector point layer with the above observer points. The ID field is numbered from 1 to 558 (which corresponds to the raster filenames)
One vector polygon layer which represents a planned new development area (much smaller extent then the raster's). This is the polygon that raster median is calculated from.

What I got so far:

By using Python console from within Qgis and Qgs.ZonalStatistics I have, with help from @Joseph, been able to extract median values by polygon (3) from all TIFF's in the folder (1) with this code:

import glob, os
from qgis.analysis import QgsZonalStatistics
vectorlayer = iface.activeLayer()
rasterfolder = 'D:/effect_test/'
os.chdir(rasterfolder)
for lyr in glob.glob("*.tiff"):
    QgsZonalStatistics(vectorlayer, lyr, attributePrefix=lyr, rasterBand=1, stats=QgsZonalStatistics.Statistics(QgsZonalStatistics.Median)).calculateStatistics(None)

Right now the statistics are added to a new column in the new planned development polygon layer I use to query the raster's with (3). The columns are named by filename with "median" added without speces behind the .tiff extention, see example below:

But this is not optimal since I now have to cut the data out of the attribute table, rearrange it in excel and then join it with the observer points vector layer, and on top of that I have to do it in two batches since 255 is the column limit in a shapefile.

How can I re-write the above code to take the median output I get from Qgs.ZonalStatistics and append it (join by id/filename) with the vector point layer containing all the observer points (2)?


texts: Which arguments are incorrect when running ArcGIS example for arcpy.mapping.ListBookmarks and exporting to PDF arcpy.mapping.ExportToPDF?

texts:
I want to use the bookmarks in ArcMap 10.2.2 and 10.3 MXDs to loop through the bookmarks, zoom to the extent of said bookmark and then export a PDF by renaming the PDF with the name of the bookmark. It seems like a simple task especially when Esri offers a greate example for this (from http://resources.arcgis.com/en/help/main/10.1/index.html#//00s300000060000000)

import arcpy
mxd = arcpy.mapping.MapDocument(r"C:\Users\Kim\Documents\Work_related\NursingSchool\P20_OnlineWellStudy\PBH_interpolationSurfaces.mxd")
df = arcpy.mapping.ListDataFrames(mxd, "Dataframe")[0]
for bkmk in arcpy.mapping.ListBookmarks(mxd, data_frame=df):
     print bkmk
    df.extent = bkmk.extent
    outFile = r"C:\Users\Kim\Documents\Work_related\NursingSchool\P20_OnlineWellStudy\PBH_interpolationSurfaces\\" + bkmk.name + ".pdf"
    arcpy.mapping.ExportToPDF(mxd, outFile, df)
del mxd

BUT, I've run this code and variances of it, with the same error:

Runtime error
Traceback (most recent call last):
  File "<string>", line 8, in <module>
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\utils.py", line 181, in fn_
    return fn(*args, **kw)
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\mapping.py", line 1139, in ExportToPDF
    layout.exportToPDF(*args)
AttributeError: Invalid destination path

What does this error mean? I don't know what I need to change in the parsing error for the ExportToJPEC. I've changed this to ExportToPDF and get the same error. I thought it might be triggered from the bookmark "printBookmark" being inaccessible, but I can access that by simple using a print statement:

Bookmark(name=u'printBookmark', extent=<Extent object at 0x23227670[0x22c36ad0]>).

Is this related to the difference between ArcGIS versions (10.1. vs. 10.2.2?)


texts: How to turn the QGIS "Coordinate capture" plugin on?

texts: I have enabled the Coordinate Capture Plugin via the Plugins menu. But I cannot seem to "turn on" or activate the Plugin and get it running. How do I do that?


texts: Floor height creation from Aerial Lidar/Imagery only

texts:
I was wondering if it is possible to generate the floor heights of buildings from aerial lidar? We don't have terrestrial lidar or stereo imagery. I don't need the ground DEM level or the roof (DSM) level but the height above the ground where the ground floor is.

The main current requirement is -A number of buildings are on stilts and built to allow for seasonal flooding so when flood models are brought in they shouldn't be classified as high risk since they are designed to allow water to flow under them. Other building are not on stilts but raised off the ground so a max flood height of 845 for example won't affect the building as it's min floor level is 845.5.

Lots of photogrammetry software using Stereo imagery/terrestrial lidar data will allow us to automate the creation of 3D building footprints (with floor height) and build 3D models for visualizing the effect of flooding (an example workflow) BUT what can we do if the only input is aerial 4ppm Lidar and 125mm orthophotos? I looked at coupling it with google streetview but many houses are not covered as they are on private roads away from public view.

See Building extraction with LiDAR data (Improve Process) and Simulate flood to assess Asset inundation footprints for background.

We have arcgis 10.1 and ERDAS Imagine/ER Mapper and I am also testing a number of lidar analysis tools like LP360, Lastools, Lidar Analyst and E3De for this.

Thanks,


texts: Are there any Objective-C frameworks for shapefile drawing?

texts: Can someone point me to library for this? Which can read shapefiles and render them?


texts: QGIS report composer to show related tables

texts:
Is it possible to include in the QGIS report composer a fields group section with tables that are related on a 1 to many relationship?
For example a point layer "HOUSES", with a KEY that relates to a table "FAM" that has a list of family members data (name, sex, age , etc). But I can´t get it to work in qgis report composer.

UPDATE:


texts: ArcGIS or FME - Line feature to linestring or multistring

texts:
There are Line feature classes as ShapeFile. The aim is to convert the data set to a CSV using linestring or multilinestring.

Is there a way in ArcGIS or FME?


texts: Group nearby points and find center point for each group

texts:
I'm currently working with a point shapefile derived from OSM that displays intersecting vertices of a road network. This shapefile has had duplicate geometries removed, resulting in a single point per link intersection. Although OSM is rather close to actual road topology, this leaves a junction with divided roads using 4 points at 1 road junction. I need to merge these 4 points into 1 to create a single point at a junction. There are occurrences where the junction will have 3 points, these will also need to be merged to 1 center point as well. A graphic displaying what needs to be done is attached. The blue points are the original junction points, while the green center point is what is needed. The shapefile consists of over 20,000 unique points, therefore manually completing this task is not efficient.

I have used Arcmap 10.6 and ArcGIS Pro (similar tools) to try converting the roads into polygons and getting the center line from the polygon, where I can then split at intersections and get single points at intersections, however this didn't work for extracting the correct points. I have also tried snapping points within a distance, however this didn't create the new center point I'm looking for. I have also tried dissolving the point shapefile based on common attributes, but the attributes are not consistent enough and this also doesn't create a new center point. I have spent hours researching this topic, but always came up empty handed.

ArcMap 10.6 and ArcGIS Pro with advanced license. My scripting knowledge is almost non-existent.


texts:
I am creating a single polygon by merging different polygons in ArcMap. When I merge them, slivers are created. So I have to delete those slivers manually. I don't have a license for ArcInfo so tools like Eliminate and others are not working.

Are there any other options for removing slivers all at once in ArcMap?

I am working on ArcGIS 10.


texts: Displaying Modified and Size in Contents panel of ArcMap's Catalog Window?

texts:
In ArcCatalog 10.0 SP2 it is possible to display the Modified and Size column values, and to sort by them.

However, in the Catalog window of ArcMap 10.0 SP2, it does not seem to be possible to do this.

The reason I am looking for this behaviour is that I have a large number of MXD files in a folder and want to be able to use the Catalog window to locate, and switch between, any of the last 10-12 that have been modified.  Without being able to sort them based on the Modified column (because it is not visible there), it is very cumbersome to scroll through them.

If this functionality is already present can somebody point me at where to find it, please?

Please note that using the Recent maps list is of no use for this because the maps I want to access have been created by a Python script and not previously opened in ArcMap.


texts: How to create one .mbtiles from thousands of GeoTIFFs?

texts:
I have about 7000 of Geotiffs (3600x2400px each).

I want to make a single .mbtiles file from them.

Can you tell me, how can I do it?


texts: Function "PropertyExists" in GeoServer style returns NullPointerException

texts:
I am using GeoServer 2.15.2.  With this as the "ogc:Filter" in my "FeatureTypeStyle" element ...

<ogc:PropertyIsEqualTo>
  <ogc:Function name="PropertyExists">
    <ogc:Literal>foo</ogc:Literal>
  </ogc:Function>
  <ogc:Literal>true</ogc:Literal>
</ogc:PropertyIsEqualTo>

... and GeoServer throws a NullPointerException (NPE):

2019-09-10 14:37:16,615 [http-apr-81-exec-9] ERROR - null [org.geotools.renderer.lite]
java.lang.NullPointerException: null
    at org.geotools.filter.function.PropertyExistsFunction.evaluate(PropertyExistsFunction.java:94)
    at blah, blah, blah

The error occurs when attempting to view the layer via "Layer Preview" in the GeoServer web app.

Has anyone successfully used this function (lately) and if so can you provide a sample?

I tried a couple different spatial data sources, one being a layer in a GeoPackage and the other a ShapeFile just for kicks.  Both result in an NPE.

Update - 2019-09-11:

This is the extent of the stack trace while still in the "geotools" package.  The remaining 125 lines are outside of "geotools" in the geoserver, tomcat and spring packages (and I don't think that portion is helpful).

2019-09-10 14:37:16,615 [http-apr-81-exec-9] ERROR - null [org.geotools.renderer.lite]
java.lang.NullPointerException: null
    at org.geotools.filter.function.PropertyExistsFunction.evaluate(PropertyExistsFunction.java:94)
    at org.geotools.filter.visitor.SimplifyingFilterVisitor.visit(SimplifyingFilterVisitor.java:477)
    at org.geotools.filter.FunctionExpressionImpl.accept(FunctionExpressionImpl.java:143)
    at org.geotools.filter.visitor.DuplicatingFilterVisitor.visit(DuplicatingFilterVisitor.java:125)
    at org.geotools.filter.visitor.DuplicatingFilterVisitor.visit(DuplicatingFilterVisitor.java:170)
    at org.geotools.filter.visitor.SimplifyingFilterVisitor.visit(SimplifyingFilterVisitor.java:557)
    at org.geotools.filter.IsEqualsToImpl.accept(IsEqualsToImpl.java:154)
    at org.geotools.filter.visitor.SimplifyingFilterVisitor.collect(SimplifyingFilterVisitor.java:228)
    at org.geotools.filter.visitor.SimplifyingFilterVisitor.visit(SimplifyingFilterVisitor.java:160)
    at org.geotools.filter.AndImpl.accept(AndImpl.java:45)
    at org.geotools.renderer.lite.StreamingRenderer.getStyleQuery(StreamingRenderer.java:1302)
    at org.geotools.renderer.lite.StreamingRenderer.getFeatures(StreamingRenderer.java:2192)
    at org.geotools.renderer.lite.StreamingRenderer.processStylers(StreamingRenderer.java:2156)
    at org.geotools.renderer.lite.StreamingRenderer.paint(StreamingRenderer.java:892)

Example feature the produces error:

wkt_geom    fid
LineString (-78.20535245444592931 38.73629880547370163, -78.20466261356823168 38.73735892913973089) 1

Full SLD that produces error:

<?xml version="1.0" encoding="UTF-8"?>
<StyledLayerDescriptor version="1.0.0"
 xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"
 xmlns="http://www.opengis.net/sld"
 xmlns:ogc="http://www.opengis.net/ogc"
 xmlns:xlink="http://www.w3.org/1999/xlink"
 xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
  <NamedLayer>
    <Name>test line</Name>
    <UserStyle>
      <FeatureTypeStyle>
        <Rule>
          <Name>rule1</Name>
            <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
              <ogc:PropertyIsEqualTo>
                <ogc:Function name="PropertyExists">
                  <ogc:Literal>foo</ogc:Literal>
                </ogc:Function>
                <ogc:Literal>true</ogc:Literal>
              </ogc:PropertyIsEqualTo>
            </ogc:Filter>
            <LineSymbolizer>
              <Stroke>
                <CssParameter name="stroke">#0000FF</CssParameter>
              </Stroke>
            </LineSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor>


texts: Import multiple xyz files and rasterize them

texts:
How can I load multiple XYZ files in QGIS?

I want to load about 100 files to create from each one a raster so we can overlay our network on them.

Manually it's too much work. I tried to create a model where I load all the files with the same parameters and then rasterize them. My thought is to run the model as a batch when its finished.

I didn't find an algorithm that loads all the xyz files with the same parameters. My file format is custom delimiters (space) and the point coordinates are the fields containing x , y values. If do this over layer > add a delimited layer it's fine. But I want to import all the files automatically from my folder without setting every time manually the parameters. The second step is to rasterize this layers then one by one. I thought the model builder was a good idea. Maybe there is another way for my doing.

How can I do this an easy way?


texts: Default transparency with GeoServer layer groups with overlapping raster layers

texts:
I have deployed two pan-sharpened raster from Pleiades imagery covering different areas with an overlap. Individually they are opaque (see here and there) but when inside a layer group (with default settings), they turn approx 50% transparent and it seems there is nothing you can change (see here). If alone in another layer group the raster are displayed without transparency (see there). It seems this is a because there is a overlap. Is there a way to solve this in GeoServer or is it mandatory to cut the imagery in order to have no overlap at all?

EDIT : adding a screenshot of one of the layers settings:

and another screenshot regarding the layer group settings:


texts: Making a pivot table in ArcGIS using Standard license

texts:
Having issued with using the pivot table tool in ArcGIS Pro due to the license level, I want to know if any geographer here knows how to make the pivot table of pivot table-like output using a set of different geoprocessing tools.

Pivot table tool is only accessible through Advanced level of license while mine is only Standard level.


texts: How to Embed KML files (Google Earth) into a website without the google gadget?

texts: I am trying to embed a KML file onto our webpage and the way it was previously done, using the KML embedding tool from Google no longer works, it seems that google has taken the gadget down or something.  Is there another way of doing this that isn't overly complex? Any Suggestions and or advice is greatly appreciated


texts: Grass i.landsat.acca says image is cloud free but it is not

texts:
I'm working on an area in Amazonia with high cloud persistance which makes global data on land cover quite useless at small scales, so I wanted to learn to process my own Landsat images, with GRASS functions.

I downloaded one Landsat 5 image from http://earthexplorer.usgs.gov/,  corresponding to 7 raster files, one for each band;
- Imported them into GRASS;
- Applied i.landsat.toar function to get top of atmosphere reflectance : The full code line is i.landsat.toar --o input=LT52270571991288CUB00_B output=LT52270571991288CUB00 metfile=/home/LT52270571991288CUB00_MTL.txt.
- Finally, I applied i.landsat.acca function to the outputs of the previous function.

What I get is the following:
"Processing first pass...

Preliminary scene analysis:
* Desert index: 0.00
* Snow cover: -nan %
* Cloud cover: -nan %
* Temperature of clouds:
** Maximum: 0.00 K
** Mean (cold cloud): -nan K
** Minimum: 10000.00 K
Result: Scene cloud free
Removing ambiguous pixels..."

This suggests that the image is cloud free, which is difficult to see looking at the image, and checking the metadata I see that the cloud cover is 29%.

Is there a step that I missed somewhere? I tried with a neighbouring Landsat scene, using exactly the same code, and in that case the correct cloud cover is found (same as in the meta data file).


texts: Showing split polygons in different layers in QGIS

texts:
I used the split feature to split a polygon into two parts.

Is it possible to show the split parts in sub-layers in the layers panel, so that I can hide them individually?

I am using QGIS 3.


texts: Overlapping polygons causing issues with Tabulate Area (Spatial Analyst) Tool

texts:
I am attempting to create a table that summarizes habitat information for a given survey area. My habitat information is in the form of a local raster file, with codes for specific habitat classes. My survey areas are circular buffers saved as a local polygon shapefile. Many of these polygons are overlapping. When I run the Tabulate Area (Spatial Analyst) tool, the values for overlapping polygons totals to less than what they should be. It appears that for overlapping polygons, only the area not overlapped by another polygon is being tabulated. If I select a single polygon and run the Tabulate Area tool on it, the values that are output are correct.

Is there a way to work around this without running the tool for each individual polygon?

Or if that is the only way, is there any suggested Python code I could use?

I am using ArcMap (ArcGIS 10.3.1).


texts: Converting Polygon to Raster in ArcMap, with a specified spatial extent

texts:
I have a vector/polygon map of UK soil %sand values which is using the British National Grid coordinate system.

I would like to convert this to a raster map with 50m x 50m resolution, with the bottom left hand corner starting at national grid 0,0 (or any number divisible by 50).

When I use the Polygon to Raster tool I am not given an option to change the extent. I have tried changing the extent of my dataframe, but the resulting raster map does not have the same extent as the one I specify.

Does anyone know how to set the extent of exported raster maps?


texts: Make an optional input parameter in QGIS Processing Algorithm Script

texts:
My version of QGIS is 2.8.3

I have a tool which takes parameters like this:

##Input_Vector_Layer=vector
##Input_Raster_Layer=raster
##Distance=number 1.0
##Height=field Input_Vector_Layer

# Do Something...

Is it possible to make the fourth parameter optional?


texts: extra records when joining tables in ArcGIS

texts:
I've been having a problem in ArcGIS 10.4.1 when performing larger joins (either attribute joins or spatial joins)--that is, with around 2,000 to 6,000 points. The join seems to work fine, but then when I export the joined tables out to a new table (all this is in one geodatabase), some of the records are repeated.  That is, I may start with 3400 points in one table, join another smaller table, and then export the result (also 3400 records), but end up with 3412 points in the new table (the 12 extra records seem to be duplicates).

Any idea why this is happening, and how to handle it?  Is it some kind of internal memory issue?

So far the only "workaround" I've found is to select and then copy all the records in the joined attribute tables out to Excel, then import them back to ArcGIS and re-plot the points, which is fairly tedious.  Splitting the datasets into 2 or 3 smaller sets also seems to resolve the issue, but again, is annoying.  Any advice would be most welcome.


texts: Making bus routes in ArcGIS Network Analyst?

texts:
I am a novice.

Aim is to make bus routes in ArcGIS Network Analyst, so specifying two points on the map, I could see what buses to ride with minimal layovers and time.

Please tell me where it is written in detail.


texts:
I want to export my point layer as CSV with the point coordinates as "xy", QGIS does the job quite well. One problem arises though, and that is that the Easting UTM coordinate is lacking the first two numbers (defining the UTM zone)

So for example, I m working in UTM zone 32(N), the output coordinates look like:

x         | y
538790.26 | 5289086.88

I need the leading "32" here ("x"), e.g. 32538790.26, which is unfortunately missing. Is there a way to export the full UTM coordinates to the .csv column?


texts:
I'm trying to merge my DEM files into 2 separate files, one will have 4 files and the other will have 3. I've been reading every post from when someone asked this before for over an hour and I'm still having problems. Here's a look at what I'm trying to do:

I'm having problems using the DEM to raster tool in ArcGIS 10.2 for Desktop (I keep on getting a 999999 error and haven't been able to locate anything that I'm doing wrong) I'm trying to do that and then use mosaic to new raster (something I found in a previous discussion). Any suggestions or a basic step by step on something better to do would be nice.


texts: Solar Terminator (day/night boundary line) generator

texts:
I'm working on a map with an orthographic projection that requires to show the solar terminator, but I've found very difficult to draw the geometry of the day/night boundary in a specific day of the year. I'm looking for a way (or a tool) to generate (or import) this line in QGIS, and being able to determine its path for each day of the year.

The only reference I could find on the forum is this question, which has a different focus, though, and doesn't provide an answer: Anybody know of a public WMS out there that serves a solar terminator?

Below a screenshot of my draft output, where the terminator has been added manually in Photoshop.


texts: Ogr2Ogr - PostGIS and using existing column for Feature ID

texts:
I am importing existing GeoJSON files into a PostgreSQL database using Ogr2Ogr. Some of these files have existing columns we want to use as unique identifiers, and some do not.

I've been unable to find a way to set a given column as a primary key or use an existing column to seed the primary key using variations of this:

.\ogr2ogr.exe -f "PostgreSQL" PG:"{connectionString}" "test.geojson" -nln test -preserve_fid

Is this something Ogr2Ogr supports? Or do I just need to hack it through SQL?


texts: Why is shapefile non recognised in the dataset?

texts:
Do CartoDB datasets accept .shp files for the geometry?

I should use them to build a map but I can't do it. Can you help me?


texts: Executing Python script from .NET AddIn and also write its standard output and geoprocessor messages to ArcMap Python Window?

texts:
Can a Python script be executed from a .NET AddIn and also write its standard output and geoprocessor message to the ArcMap Python Window?

I have no problem call python scripts with Process.Start or GeoProcessorClass.Execute.  But I have an existing AddIn (10.1) and I want to add a new button that will call the python script and have it behave the same way a right-click Load... behaves in the ArcMap python window.

Is there a way to do this?

I'm currently using this procedure.

And it does run the script--I'm just after the same out to the python window.


texts: Calculating areas of different classes in Google Earth Engine?

texts:
I am attempting to calculate the area of 5 different classes in my classified map. GEE tutorials showed me a method (https://developers.google.com/earth-engine/tutorial_forest_03#calculating-pixel-areas), however I struggled to apply it because I dont understand how to select a specific class before using .reduceRegion to calculate area because the example is quite different to my map.

Can someone guide me through how to calculate area, whether using the same method, or another?

Here's the link to my map: https://code.earthengine.google.com/1171bdeecd776800b334a49f60d5d443

Here is my code

var image=ee.Image('LANDSAT/LC08/C01/T1_SR/LC08_169062_20180205');

//define parameters

var green = image.select('B3');
var swi = image.select('B6');
var mndwi = green.subtract(swi).divide(green.add(swi)).rename('MNDWI');

var mndwipara = {min: 0, max: 0.7, palette: ['white', 'blue']};

Map.addLayer(mndwi, mndwipara, 'MNDWI');

//apply threshold to select only positive pixels
var lake_mask= mndwi.gt(0);
Map.addLayer(lake_mask);

//classification attempt - scum
var scum=image.expression(
'(nir)*0 + (nir < 0.018)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt - sediment
var sediment=image.expression(
'(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green>
nir/red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
 });

//classification attempt - lowbiomass
var lowbiomass=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir < red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
 });

//classification attempt-scum
var scum=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red > 4.07)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
 });

//classifcation attempt bleached scum
var bleached_scum=image.expression(
'(nir)*0 + (nir > 0.018 && red > green && green/nir > 1.3 && green > 0.065)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt - microphytobenthos
var microphytobenthos=image.expression(
'(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green< nir/red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
 });

//classification attempt - highbiomass
var highbiomass=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red < 4.07)*1', {
 'nir': image.select('B5'),
 'red': image.select ('B4'),
'green': image.select('B3'),
 });

/*
lake = 1
scum = 2
bleached_scum = 3
microphytobenthos = 4
highbiomass = 5
sediment

*/

var lake_mask_ = ee.Image(1).mask(lake_mask).toInt();
var scum_ = ee.Image(2).mask(lake_mask.mask(scum)).toInt();
var bleached_scum_ =
ee.Image(3).mask(lake_mask.mask(bleached_scum)).toInt();
var microphytobenthos_ =
ee.Image(4).mask(lake_mask.mask(microphytobenthos)).toInt();
var biomass_ = ee.Image(5).mask(lake_mask.mask(highbiomass)).toInt();

 var test = ee.ImageCollection([lake_mask_,scum_,bleached_scum_,microphytobenthos_,biomass_,]);

 var test = test.reduce(ee.Reducer.max());

 print(test);
 Map.addLayer(test,{min: 1, max: 5, palette: ['0000FF',
 '1BCFFF','66ff33','F4FF0B','FA0007']},'Classification');


texts: Layer labels not showing properly on the composer

texts:
I'm currently working on a map for the different brazilian biomes.

 I have to label them and show it on my final map, but when I'm working on the composer, this happens:

That's the exact size of my map. If you take a close look, you'll see that the "CERRADO" word appears cut in the middle on the upper- left side. I've looked around and found a pretty similar situation here but the solution proposed doesn't apply. It still displays like the second picture. I need all three labels to be inside the map window I created on the composer. Any thoughts?


texts: Understanding strange measure distance results from QGIS?

texts:
My colleague just told me about incorrect measure distance tool results in QGIS. For reference we used this pond.

While our local map provider Mapy.cz

and Google Maps:

both show about 970 meters distance, QGIS shows me about 1500 meters:

I first though it is OTF reprojection issue so I switched it off. But problem remains. The results are same for various Openlayers map services (OSM, Bing etc.) but also for OSM derived offline data.

I have my data in EPSG:3857 - WGS84 Web Mercator in order to be able to let the beginner users switch between offline and online data without slowing down the map loading and avoiding the need of OTF reprojection.

Is there something wrong with this particular projection? I would not expect such a huge differences for this small area.


texts:
I am trying to plot raster map using Pyhton. Image is not being saved, though it appears.
Moreover, pixels are not appearing on the exact place. I am getting image as shown below (which is WRONG, however, lat-lon labeling is correct.)

Below is the correct orientation, which I have exported from ArcGIS.

There is some mistake in placing the pixel at right location, but I ma not getting it. How to correct it?

Here is my script.

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from osgeo import gdal
from numpy import linspace
from numpy import meshgrid
import numpy as

#################################################################################

ds = gdal.Open('Path\\To\\Raster.tif')
data = ds.ReadAsArray()

gt = ds.GetGeoTransform()
proj = ds.GetProjection()

#################################################################################

xres = gt[1]
yres = gt[5]

xmin = gt[0] + xres * 0.5
xmax = gt[0] + (xres * ds.RasterXSize) - xres * 0.5
ymin = gt[3] + (yres * ds.RasterYSize) + yres * 0.5
ymax = gt[3] - yres * 0.5

x_center=(xmin+xmax)/2
y_center=(ymin+ymax)/2
#################################################################################

fig = plt.figure(figsize=(18, 10))
ax = fig.add_subplot(111, axisbg='w', frame_on=True)

m = Basemap(llcrnrlon=xmin,llcrnrlat=ymin,urcrnrlon=xmax,urcrnrlat=ymax,
             projection='tmerc', lat_0 = y_center, lon_0 = x_center)

parallels = np.arange(15.,30.,0.25)
m.drawparallels(parallels,labels=[1,0,0,0],fontsize=12, linewidth=0.4)

meridians = np.arange(70.,90.,0.25)
m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=12, linewidth=0.4)

x = linspace(0, m.urcrnrx, data.shape[1])
y = linspace(0, m.urcrnry, data.shape[0])

xx, yy = meshgrid(x, y)

m.pcolormesh(xx, yy, data, cmap=plt.cm.jet)

plt.show()
plt.savefig('Path\\To\\Save_Image.png', bbox_inches='tight', pad_inches=.2, dpi=600)


texts:
I'm making an evaluation to replace the TomTom commercial Database (ex TeleAtlas) with Open Street Map data.

What I need to do is use the street map as a base layer for my WebGIS applications. I use MapServer as application server and I made styles for features of the TomTom DB based on attributes values. I want to do the same with OSM data.

Looking on the web I found many references to MapBox, TileMill, preconfigured style schemas and so on... but I have some questions/doubts about all that.

1) If I need to render this data through MapServer, can I avail of all these instruments?

2) Is there a way to convert those styles to MapServer .MAP files?

3) I imported OSM data from SHP files downloaded from a web site that regularly creates them from updated OSM data. I read something about OSM2PGSQL, OSM2PO, OSMOSIS and so on... but I don't really understand what they are and how they work and if I can use them with my actual database. Do I really need to import OSM data with such tools?


texts: Vector tile layer over Google Maps or MapKit in iOS?

texts: I'm looking for a way to add a vector tile layer instead of the built in raster layer for Google Maps or MapKit in iOS but I can't seem to find anything that isn't a complete replacement like MapBox.  Does this exist or is there a way to render the vector to UIImage (raster) on device.


texts:
I am trying to use the "Tabulate Area 2" tool in the following update:

http://blogs.esri.com/esri/arcgis/2013/11/26/new-spatial-analyst-supplemental-tools-v1-3/

I have overlapping polygons in which I'd like to calculate the area of lake coverage, hence why the original "Tabulate Area" will not work.

Here is a picture of the output table when I run "Tabulate Area 2":

The feature zone is working fine, using the ID column I am inputting, but the feature class is using ObjectID and not using any other attribute, see the example below:

The areas that it is calculating are not in metres, I think they are in decimal degrees.

How can I change the class field the "Tabulate Area 2" tool is using for my feature class? and How can I make this tool calculate area in metres?

I am using ArcGIS 10.1 with an ArcView license.


texts: QGIS 3.2.0 on Mac Sierra 10.12.6 not loading Processing / Vector / other toolbars

texts:
I've already changed the Environment to enable "Use Custom Variables" and added a "Prepend" "PATH" ""/Library/Frameworks/GDAL.framework/Programs:/Library/Frameworks/Python.framework/Versions/3.6/bin:""

I'm behind a corporate firewall, but I've updated my proxy settings to make sure plugins are working correctly. The "Processing" plugin shows up under "installed."

But "processing" is still missing from the toolbar, from the list under System > Options. So I'm not even able to reset that to the default.

I did not install this first in another language. This was directly downloaded from: https://qgis.org/en/site/forusers/download.html


texts: Seeking smallest device able to interface with GPS unit?

texts:
For GPS units that talk over USB, what's the smallest device you could plug the GPS in to in the field, in order to upload from the device to the GPS and download from the GPS to the device?  Netbook?  Smartphone?  Ipad? Kindle Fire? other?

I'm looking for more of a turnkey end-user out-of-the-box (insert buzzphrase here) solution that's on the mainstream market.

The application is uploading search segments boundaries and insertion routes for search and rescue teams, before the full IC computer setup arrives.


texts:
How do I load this REST service in QGIS?

https://tpwd.texas.gov/arcgis/rest/services/GIS/Ecosys1/MapServer

The ArcGIS REST API Connector returns the error:

[Errno 1] _ssl.c:504: error:14090086:SSL
routines:SSL3_GET_SERVER_CERTIFICATE:Certificate:certificate verify failure


texts: Applying Symbology From Layer to Include Label Properties, Feature Scale Range, etc

texts:
Using the Applying Symbology from Layer tool appears to not capture a layer's table of contents visibility, label properties, feature scale range, or label scale range. I have many feature classes that are repeatedly used in numerous MXDs for which I need to change these properties. I need a process that is similar to said tool because I can run it in a fairly short script and apply it to all my MXDs.

I could use a script similar to that below, but would be required to (in addition to adding lines for the other properties) alter the feature class in ListLayers and the applicable properties many times in order to edit all the feature classes as needed.

mxdList = ['MXD1', 'MXD2', '...']
fcList = ['FC1', 'FC2', '...']

for doc in mxdList:
    mxd = arcpy.mapping.MapDocument('Path\\' + doc)
    for item in fcList:
        for lyr in arcpy.mapping.ListLayers(mxd, item):
            # lyr filename matches FC ToC name less extension
            arcpy.ApplySymbologyFromLayer_management(lyr,
                r'path\lyrFiles\/' + item + '.lyr')
            if lyr.supports('VISIBLE'):
                        lyr.visible = False
            # continue with other properties...
arcpy.RefreshTOC()

Is there a quick way to change these other properties for many MXDs? This would not be an issue if lyr files captured the properties of table of contents visibility, label properties, feature scale range, or label scale range.

[Edit]

I just found this link, which leads me to think ArcObjects in Python could be a solution to my question. This looks like a can of worms, however, for someone who is fairly new to the use of Python. Are there other alternatives besides extending code above, using ArcObjects, or dragging feature classes from one ArcMap window to another?


texts: Cropping a point cloud with PDAL also strips it off

texts:
Description of the problem:

When cropping a .ply point cloud using PDAL pipeline:

    {
      "pipeline":[
        "input.ply",
        {
          "type":"filters.crop",
          "bounds":"([52435.72, 63951.79], [53477.36, 58371.35])"
        },
        {
          "type":"writers.ply",
          "filename":"output.ply"
        }
      ]
    }

it also "strips" the point cloud.
Here are some screenshots from CloudCompare:

Fig. 1 PDAL Cropped cloud view from above (xy plane).

Fig. 2 PDAL Cropped cloud view from front (xz plane).

Fig. 3 Original point cloud view from front (xz plane). Exact same viewing position and orientation as the second image above. The original point cloud is denser and more homogeneous.


texts:
I want to work with kml/kmz files in the Google Maps API, but I want to remap the color scales.

Is that possible?

Here are example kml/kmz files where I'd like to change the color representation of the data, but I don't see any method to do so.  I can see them fine using Google Earth or the Google Maps API.

world population density:
https://neo.sci.gsfc.nasa.gov/servlet/RenderData?si=875430&cs=rgb&format=KMZ&width=3600&height=1800
world topography:
https://neo.sci.gsfc.nasa.gov/servlet/RenderData?si=196466&cs=rgb&format=KMZ&width=3600&height=1800
US current snow depth:
http://www.nohrsc.noaa.gov/snow_model/GE/latest_nsm_depth.kmz


texts:
I have successfully configured Mapnik and generate_tiles.py on my computer so I can generate Mapnik tiles from the OpenStreetMap map database.

However, I would like to generate tiles similar to OpenCycleMap's ones. If I'm right tile configuration is in the osm.xml file. Is this configuration of OpenCycleMap available somewhere on the internet? If yes, where?


texts: Choosing most suitable GIS for archaeology?

texts:
I will start on project in archaeology where I'll have to start using GIS. I only tried using ArcGIS a while ago but for the rest i have no experience.

What platform would you suggest I start on?

I'm looking for something that has analysis tools I need for archaeology (least cost route, line of sight..) and it doesn't take ages to learn how to use it. I would prefer open source, but that isn't mandatory (university has licenses).

More important is how steep is the learning curve for the software?


texts: Showing same layer multiple times using ArcMap

texts:
I am using ArcMap.

I am working with a layer for which I have two separate values that I want to display using symbology. Ideally, I would like to duplicate the layer and place them side by side, with one image displaying symbology for the first variable and the other for the second variable. When I try to copy and paste the layer however, the labels for both variables show simultaneously on the same image.

How can I create two separate layers?


texts:
I am making a web application with OL3 and have integrated a 3d viewer with cesium. I have two data sets that are being loaded with GeoJSON as vector WFS. Both data sets have height information contained and are published as such in GeoServer (Elevation is set).

I am trying to give one WFS an extrusion (building locations) and the other a height (aircraft elevation). Everything I have found seems to be based off of cesium geometry. Is it possible to do this? Do I need to convert to primitives first?

** EDIT **
Here is the buildingSource script:

var geojsonFormat2 = new ol.format.GeoJSON();
var buildingSource = new ol.source.Vector({
    loader: function(extent, resolution, projection) {
    var url = 'http://local:58088/geoserver/wfs?service=WFS&' +
    'version=1.1.0&request=GetFeature&typename=Thunderbirds:Nevada&' +
    'outputFormat=text/javascript&format_options=callback:loadFeatures2' +
    '&srsname=EPSG:3857&bbox=' + extent.join(',') + ',EPSG:3857';

$.ajax({
    url: url,
    dataType: 'jsonp',
    jsonp: true
    }).done(function(response) {
        WFSformat = new ol.format.WFS(),
        sourceVector.addFeatures(WFSformat.readFeatures(response)) }); },
strategy: ol.loadingstrategy.tile(ol.tilegrid.createXYZ({
    maxZoom: 19 })) });


texts: How to create polyline of varying distance from origin point using azimuth?

texts:
I am looking for a Python script that does the following:

Uses a point (from a point shapefile) as the origin.
Creates about 3 km line from point
Based on an azimuth provided in the original point shapefile.
Creates new entry in another shapefile with the new line

I imagine this involves trig, but I can't decide the best way to go about it.


texts: How to build a mobile app from Flex web app

texts:
We have a Flex application based on ArcGIS flex API in the final stage of development and ready to move to testing/production. The application uses some map services and a complex GP tool. My boss wants a quick mobile app out of it. I have never worked with mobile app. I want to learn- what it takes to convert this application to mobile app- software, coding effort. Some useful links related to ESRI or Adobe will be helpful.

Thanks


texts: MapServer Integration with ZOO Project on Windows 10

texts:
I installed ZOO Project using the package installer (x64-72). To add MapServer, I installed ms4w using the zip file (ran apache-install.bat) and added the zoo project application. I also added the mapserverAddress pointing to the mapserv.exe on Apache in main.cfg of ZOO Project (running on IIS).

When I execute HelloPy or any other service with "useMapserver = true" in the services configuration file, I get the following error

  ZOO Kernel failed to process your request, receiving signal 11 = SIGSEGV

Also, in the dataPath directory, two files, a .data file and a .zip file are added. If I pass in a shapefile's zip for HelloPy, the .data file shows some unreadable content while the .zip contains a .cst file, which in turn only shows some character encoding (ISO ... ).

Can someone please help me with this issue?

EDIT: I found that a part of my issue was that I was not passing in a file that could be understood by GDAL. On making this change I can see that the .map file is generated along with the two files mentioned above. However, the SIGSEVG error still exists.

EDIT 2:
If I pass a WMS result to MapServer (say image/jpeg, image/geotiff) a mapfile is generated along with the other 2 files but the SIGSEGV error persists. If I pass in image/geotiff8 instead, then the result is a link for a getmap request. However, this link does not show any result.

If WFS results are passed in instead, like a kml, then the result is a vector map instead of the desired raster. What should I do?


texts: How to add custom REST controller with @RestController annotation?

texts:
I have this dilemma.

As from 2.14 version of GeoServer, there is a new way to write rest controller
using @RestController annotation and extending RestBaseController.class

But after starting the server I can't see the custom rest path.

How do I create custom REST services with the new way?


texts: gdal_merge to merge hundreds of png files

texts:
I need to merge about 760 png files into one and am looking at using gdal for this. This is related to the following question Getting a constant background color for transparency

P:\2012\183_TownPlanning_Symbology\Working\Raster_Layer_Creation\tests\BA\25k>gd
al_merge -o ba_25k_gdal -of png -init "255 255 255" *.png

should work but gives

  Format driver png does not support creation and piecewise writing.
  Please select a format that does, such as GTiff (the default) or HFA
  (Erdas Imag ine).

When I switch to HFA or GTiff format and try

P:\2012\183_TownPlanning_Symbology\Working\Raster_Layer_Creation\tests\BA\25k>gd
al_merge -o ba_25k_gdal -init "255 255 255" *.png

I get

  ERROR 4: `*.png' does not exist in the file system, and is not
  recognised as a supported dataset name.

  Traceback (most recent call last):   File
  "C:\OSGeo4W\bin\gdal_merge.py", line 509, in
      sys.exit(main())   File "C:\OSGeo4W\bin\gdal_merge.py", line 392, in main
      ulx = file_infos[0].ulx IndexError: list index out of range

but the format seems correct according to http://www.gdal.org/gdal_merge.html

How can I tell it to include all png files in the directory? I can pipe *.png > png_list.txt but can I get gdal to read a list from a text file?
So does this mean merging png's in gdal is not possible?

Are there other ways of achieving this? I have tried Global Mapper but that estimates a day to complete it and I am just running it. Erdas crashed. I am trying ArcMap 10.1 import into mosaic dataset and export as raster -still running but estimates a 35gb png file.

This is from 161 tiles (25k resolution) over 65,000 sq km. But the source png's are only 73Mb.


texts: How long does geoserver store the cookie it generates from authentication via j_spring_security_check?

texts: I have a website which connects to a geoserver, it fetches sessions from geoserver via j_spring_security_check. How long is the cookie given by geoserver valid and can I edit its expiration?


texts: Details on 'pixel' and 'line' variables for gdal_translate gcp inputs

texts:
I am attempting to use gdal_translate and gdal_warp to georeference an image similar to the plugin for QGIS "Georeferencer". When I use QGIS, georeferencer produces this GDAL script:

gdal_translate -of GTiff -gcp 1.47803e+06 -5.10486e+06 1.47803e+06 5.10486e+06 -gcp 1.47862e+06 -5.10514e+06 1.47862e+06 5.10515e+06...................

gdalwarp -r near -order 1 -co COMPRESS=NONE  ...........

What are the 'pixel' and 'line' (the first two variables after -gcp )?

I had assumed they were the x and y coordinates of the pixel to place the GCP, but my tiff isn't 1.5e+6 pixels wide, and I dont understand why the line would have a negative value.

here is the gdalinfo response for the tiff I am trying to georeference:

Driver: GTiff/GeoTIFF
Files:     C:\Users\PC_2\Downloads\uav_processing_pipeline\testing_z_field\mosaics%2F2018-    11-02_with_laz%2F8-fid12-park_rd_rcd30_2018-11-02_green_nogcp.tif
Size is 19206, 16391
Coordinate System is:
PROJCS["NZGD2000 / New Zealand Transverse Mercator 2000",
GEOGCS["NZGD2000",
    DATUM["New_Zealand_Geodetic_Datum_2000",
        SPHEROID["GRS 1980",6378137,298.257222101,
            AUTHORITY["EPSG","7019"]],
        TOWGS84[0,0,0,0,0,0,0],
        AUTHORITY["EPSG","6167"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4167"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",173],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",1600000],
PARAMETER["false_northing",10000000],
UNIT["metre",1,
    AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","2193"]]
Origin = (1477791.815500000000000,5105524.866330000600000)
Pixel Size = (0.074660000000000,-0.074660000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_SOFTWARE=pix4dmapper
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 1477791.816, 5105524.866) (171d28'14.69"E, 44d11'35.73"S)
Lower Left  ( 1477791.816, 5104301.114) (171d28'13.67"E, 44d12'15.38"S)
Upper Right ( 1479225.735, 5105524.866) (171d29'19.27"E, 44d11'36.59"S)
Lower Right ( 1479225.735, 5104301.114) (171d29'18.25"E, 44d12'16.24"S)
Center      ( 1478508.775, 5104912.990) (171d28'46.47"E, 44d11'55.99"S)
Band 1 Block=19206x1 Type=Float32, ColorInterp=Gray
  NoData Value=-10000


texts: ArcMap attribute joins generating incorrect values

texts:
Background:
At 10.8.1, ArcMap attribute joins run without complaint, yet generate incorrect results.
As shown in the screenshot below, the tables component (parent) and chorizon (child) are joined 1:1, each via the field cokey.
Both are file geodatabase attribute tables, stored in the same FGDB.
I have conducted the join using two different techniques:

From the attribute table window dropdown > Joins and Relates > Join... (both with the Keep all Records and Keep only Matching Records options, always using the Validate Join button)
Using the Add Join tool (with Keep All Target Features unchecked)

Problem:
When reviewing the output, I notice that some of the child table values are incorrect following the join.  For example, the screenshot shows two selected rows where the joined sandtotal_r values are different than the source table (6.8 vs. 63.7 and 13.5 vs. 29.5).  In fact, 6.8 and 13.5 do not exist anywhere in the source attribute table!
The incorrect values always appear, regardless of the join method used.
When I conduct the join using QGIS 3.18.2, the output is correct.  Thus it appears that the problem is with ArcMap, not the data.
Any thoughts?
@FelixIP suggested I try the Join Field tool, which worked correctly.  However, the original problem still exists.  I suspect a bug.


texts: How can I insert the result of the dijkstra algorithm as a feature in a table with PyQGIS?

texts:
I ask for your help with the PyQGIS API.
Being a new user of GIS I would like to carry out an operation.

It is for me to look for the shortest path in a graph using the DIJKSTRA algorithm.

The code of the QGIS doc after some small adaptations works well except that I do not recover anything at the end if it is not a graphic plot. After the execution I would like to retrieve the shortest path found (the result) as an entity that I insert in a table and so on for possibly other shorter paths that I will calculate.

from PyQt4.QtCore import *
from PyQt4.QtGui import *
from qgis.core import *
from qgis.gui import *
from qgis.networkanalysis import *

vl = qgis.utils.iface.mapCanvas().currentLayer()
director = QgsLineVectorLayerDirector(vl, -1, '', '', '', 3)
properter = QgsDistanceArcProperter()
director.addProperter(properter)
crs = qgis.utils.iface.mapCanvas().mapRenderer().destinationCrs()
builder = QgsGraphBuilder(crs)

pStart = QgsPoint(-0.835953, 0.15679)
pStop = QgsPoint(-1.1027, 0.699986)
tiedPoints = director.makeGraph(builder, [pStart, pStop])
graph = builder.graph()
tStart = tiedPoints[0]
tStop = tiedPoints[1]
idStart = graph.findVertex(tStart)
idStop = graph.findVertex(tStop)
(tree, cost) = QgsGraphAnalyzer.dijkstra(graph, idStart, 0)

if tree[idStop] == -1:
  print "Path not found"
else:
  p = []
  curPos = idStop
  while curPos != idStart:
    p.append(graph.vertex(graph.arc(tree[curPos]).inVertex()).point())
    curPos = graph.arc(tree[curPos]).outVertex();
  p.append(tStart)
  rb = QgsRubberBand(qgis.utils.iface.mapCanvas())
  rb.setColor(Qt.red)
  for pnt in p:
    rb.addPoint(pnt)

May somebody help?


texts: When trying to load a shapefile into QGIS, a little photograph icon comes up instead of a line.

texts: QGIS does not seem to recognize the geometry or extents of the file, even though they must be present, as the site where I got the file includes a preview of the line shapefile over Bing Maps. The shapefile is from PASDA, and all other files I have downloaded from this website have worked without problems. The link to the website is as follows http://www.pasda.psu.edu/uci/DataSummary.aspx?dataset=888


texts: gdal_translate or gdalwarp to define projection

texts:
Which of these is the best approach to define (AKA. set) a projection for a tiff?

1) gdal_translate -a_srs EPSG:32182 in.tif out.tif

2) gdalwarp -s_srs EPSG:32182 -t_srs EPSG:32182 in.tif out.tif


texts: Applying Different Label to Data Driven Pages Current Target?

texts:
I'm exporting the pages to png files. The pages are divided by lots. Each page has a label with the lot number.

In the case of the sample included, 18-385 is the lot being exported. Pieces of several other neighboring lots are also printed.

What I would like to do if possible is format the label of the lot being exported different from the neighboring lots.

The first picture is what I have so far, and the second picture is what I would like to do.

I am very new to ArcGIS.
I'm doing the exports in the Geoprocessing->Python box. I don't know yet how to get to a full IDE where I can save code.
I'm using ArcGIS 10.1

If it helps, here's the code I'm using.

mxd = arcpy.mapping.MapDocument("CURRENT")
ddp = mxd.dataDrivenPages
basePath = "C:\\Downloads\\Dev\\ArcGIS\\Fulton\\"
for pageNum in ddp.selectedPages:
   mxd.dataDrivenPages.currentPageID = pageNum
   distNum = mxd.dataDrivenPages.pageRow.DistNum
   savePath = basePath + distNum
   fileName = savePath + "\\" + mxd.dataDrivenPages.pageRow.DistPage + ".png"
   print fileName
   if not arcpy.os.path.exists(savePath):
      arcpy.CreateFolder_management(baseDir,dirName)
   arcpy.mapping.ExportToPNG(mxd, fileName,resolution=300)
del mxd

What I have:

What I want:


texts: Overlay with the feature collection properties


texts: How can I use the OSM SRTM plugin for Osmosis program?

texts:
I am using Osmosis and I have a test .osm map. I want to install a plugin from the following Github repo:

https://github.com/locked-fg/osmosis-srtm-plugin/

I have followed the steps to "install" it by creating a 'plugins' folder and placing the binaries of the plugin in the floder as mentioned in there. Unfortunately, it shows the following error on running this command:

osmosis.bat --read-xml myMap.osm --write-srtm --write-xml mySrtmMap.osm

This gives the following command-line errors:

Errors screenshot

Can anybody please explain if I am missing something?


texts: Getting longitude/latitude coordinates from Addresses

texts:
I got a CSV file with more than 15.000 address, and I need to show them on a map of the South-West of France.
I've found the MMQGIS Plugin on QGIS, but it does so without the right SCR.
Do you have any idea to help me convert my addresses? I searched on the Internet and didn't find any adequate solution.


texts: Geoserver to PostGres database using a hostname

texts:
I am trying to connect GeoServer to a PostGres Database with PostGIS extensions installed, but I need to use the host name becauuse the IP address may change.  I am running geoserver and postgis in linked docker containers (see https://docs.docker.com/userguide/dockerlinks/), and restarting a container may result in a change of IP address.  However, I have full control over the hostname.  My /etc/hosts contains the entry:

172.17.0.244    postgis

If I configure GeoServer to use host 172.17.0.244 it connects fine, but if I specify postgis if doesn't connect.  Is this just not possible with GeoServer, or am I doing something wrong?


texts:
I have these spatial points:

long <- c(-1.769717, -1.767766, -1.769170, -1.768462, -1.768389)
lat <- c(54.92214, 54.92207, 54.92113, 54.92167, 54.92165)
points <- data.frame(long, lat)
coordinates(points) <- c('long', 'lat')
proj4string(points) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
points <- spTransform(points, CRS("+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"))##

I can calculate if any two points fall within 20m of each of other like this:

points_matrix <- gWithinDistance(points, dist = 20, byid = TRUE)
diag(points_matrix) <- NA

      1     2     3     4     5
1    NA FALSE FALSE FALSE FALSE
2 FALSE    NA FALSE FALSE FALSE
3 FALSE FALSE    NA FALSE FALSE
4 FALSE FALSE FALSE    NA  TRUE
5 FALSE FALSE FALSE  TRUE    NA

If any two points lie within 20m of each, I need to remove one of these points. The points_matrix shows that points 4 and 5 lie within 20m of each other. How can I remove either 4 or 5 (doesn't matter which one) from points and keep all other points.


texts: Using ogc:filter and ogc:and gives me an empty result

texts:
I want to render my points using the point stacker transformation. I am creating a style  so that when zoomed in far enough the individual points are in 4 different categories using 4 different colors. To do this I'm trying to use 4 different rules and in each rule I  try to filter the points based on two things.
The first where 'count' Is Equal To '1' where the attribute count has been generated by the point stacker.
Second where 'Attribute' Is Equal To 'Attribute value' which is already stored in the data and consists from string values.

My problem is that when I am trying to to filter the points using 'ogc:And' my output is empty.

My code using the filters is:

    <Rule>
      <Name>rule1</Name>
      <Title>MyTitle</Title>
      <ogc:Filter>
       <ogc:And>
        <ogc:PropertyIsLessThanOrEqualTo>
          <ogc:PropertyName>count</ogc:PropertyName>
         <ogc:Literal>1</ogc:Literal>
        </ogc:PropertyIsLessThanOrEqualTo>
        <ogc:PropertyIsEqualTo>
          <ogc:PropertyName>AttributeColumn</ogc:PropertyName>
         <ogc:Literal>MyValue</ogc:Literal>
        </ogc:PropertyIsEqualTo>
       <ogc:And>
      </ogc:Filter>
      <PointSymbolizer>
        <Graphic>
          <Mark>
            <WellKnownName>circle</WellKnownName>
            <Fill>
              <CssParameter name="fill">#FF0000</CssParameter>
            </Fill>
          </Mark>
          <Size>8</Size>
        </Graphic>
      </PointSymbolizer>
    </Rule>
    <Rule>

When I am using the ogc:And statement the output of this rule is empty. When I use ogc:Or there is an output, which makes sense. However I can only use a rule with ogc:Or once because it is always true and geoserver ends up displaying all points in the same color regardless of the other rules set. So I need to use ogc:And. I am very new to geoserver and sld and I do not understand why my statement is giving me an empty output. When I validate my statement it gives me no errors and the rest of the code works fine.

After these 4 rules I have 2 more rules that define the clustering of the points, I basically copied that part from the OpenGeoSuite website and it works fine.

geoserver version is 2.8.2

Any thoughts why I get an empty output or perhaps a workaround this problem?


texts: PostGIS: Duplicated Rows when Drawing State Boundaries

texts:
I have two shapefiles. One is the DMA and another of states. What I'm trying to do is to draw a state boundary on DMAs (which can cross state boundaries).

I've asked this question and I recieved an answer, which solve the problem but then I noticed there are duplicates in the polygons. Below is the query that creates the state line boundaries.

CREATE SEQUENCE polyseq_1;

CREATE TABLE boundaries AS
SELECT
  nextval('polyseq_1') AS id,
  b.name as state_name,
  a.dma_1 as dma_1,
  CASE
     WHEN ST_Within(a.geom,b.geom)
     THEN a.geom
     ELSE ST_Multi(ST_Intersection(a.geom,b.geom))
  END AS geom
FROM tl_2015_us_state_4326 b
JOIN dma_boundary a
ON ST_Intersects(a.geom, b.geom);

CREATE TABLE boundaries_polygons AS
SELECT *
FROM (SELECT id, state_name, dma_1, (ST_Dump(b.geom)).geom
from boundaries b) as temp
WHERE ST_GeometryType(temp.geom) = 'ST_Polygon'
;

Below is the results, you can see that there are duplicates of the same DMA in the same state. Ideally, I'd like to merge all these polygons together into one by their name and state. For instance, combine all the Washington-DC DMA for Washignton, then all the Washington-DC DMAs for Virgina, and all the Washington-DC DMAs for Maryland into 3 distinct parts.


texts: Passing conditional statements into arcpy.da.UpdateCursor

texts:
How do you pass conditional statements from a list into a data access module update cursor?

I am using ArcGIS v10.3.1.

a = r"...\Default.gdb\New_Shapefile"
ConList = ["r[1] >= 1"]
with arcpy.da.UpdateCursor(a, ["Id", "apples"]) as uCur:
    for r in uCur:
        if ConList[0]:
            r[0] = 99
            uCur.updateRow(r)

Result is shown, highlighted rows should hold the pre-existing value of zero.


texts: Limitations of iterate row selection in ArcGIS ModelBuilder?

texts:
inb4 using real script ;), I'm looking for some advice on arc's model builder.

The iterate row selection in Arc's model builder allows you to make multiple table views based on field groupings of a table. These table views then can be used for other processes in the model.

I'm using the table views to calculate multiple variables which are then populated into a different feature dataset. More or less, I'm using the iterate row selection to run processes on sets of related records (MANY) in one table and then using the results to populate a second table (ONE).

One of the fields in the second table requires a variable that is calculated in part by a group of records that is not the group of records identified by the iterate row selection. In other words, thus far each field only requires the set of records made into a table view by the iterator. But as mentioned above, this new field requires a table view not based on the selections made possible in the iterate row selection--a different relationship is required.

There are some problems with this in model builder. My approach now is to take the inline variable output of the iterate row selection and use it in an sql expression of a seperate Make Table View within the model. This Make Table View is based on the original table (used by the iterator) with the variable in an sql expression such that a different column is used to base the make table view on for each iteration. example: "c%Value%" = 1. (in the table there are a number of column headers that begin with the letter c and then a number. |c1|c2|c3|c4|...c24|. My reasoning is, with each iteration, a new %Value% is produced (1-24), and with each new %Value% the separate Make Table View "should" produce a table view based on the inline variable in the expression, which then can used to make the necessary calculations which then can be calculated into the second table's problematic field.

After running the model, I find that the fields thus far calculate correctly, but the problematic field has values that are quite wrong. Running the make table view independently produces the correct variables as well calculating those values into the target dataset, but when ran with the iterator, it does not work.


texts: How to clean up a polygon by removing vertices and lines within?

texts:
After merging a set of shapefiles I've ended up with polygons that have vertices and lines within them. Is there a simple tool to remove these vertices and clean the polygon in ArcMap 10.2?

I realise the same question has been asked here, but it does not state the GIS-program and I cannot find the tool mentioned in the answer (delete ring in the advanced digitising toolbar).

EDIT: when zooming in very far, the lines and dots appear to be tiny triangles. They are not separate features, but holes in the polygons that arose from merging polygons that had not been lined up nicely. I am now looking for a way to repair there holes.

I reckon there is a simple tool to do this, but I cannot seem to find it.


texts:
I am using QGIS and qgis2web to make a simple historical map. I have many problems because I am an amateur, but I am solving them through a search.

But as you can see in the picture, I have a problem with my font now!

China is output properly set to 'points' but other countries are set to 'map units' and are not output properly!

I made the font type 'Points' at first, and it worked well on qgis2web. Second, I used 'Meters at scale' but it didn't work normally. I searched and found that qgis2web does not yet support Metres at Scale, and instead, Map units will work.

So I tried setting it as 'Map units' for the last time, and it doesn't work as well. I want to make the font small when zoom out and big when zoom in. So I used Map units... :(

What should I do?

Oh, I'd like to ask you another question. I set up a font separately, but qgis2web wiki explains 'The export only specifies which font a label uses - it does not export the font itself'! I tried to put the .woff2 fonts file in the webfonts folder, but it doesn't work properly.

What should I do?


texts: Distance between two layers along the road network

texts:
I have four layers. They are:

Crime
School
Census tract
Street

I want to calculate the distance between one school to each crime incident along with the road network. Moreover, Those school and crime incidents should be within one census tract. For example, the attached census tract has two schools (Green color-A, B).
How can I calculate the distance between A to each crime within the particular census tract and B to each crime incident in that census tract along with the road network?


texts: Project only required shapefiles in folder - update

texts:
I have updated the script from my previous question;

Project only required shapefiles in folder

I can't get the if fcspatialRef != spatialRef: to work. The script won't filter the shapefiles in ArcGIS, any advice/suggestions?

import arcpy

# Set Workspace
arcpy.env.workspace = arcpy.GetParameterAsText(0)

#Access shapefiles
fcList = arcpy.ListFeatureClasses()

try:
# Projection template
  template = arcpy.GetParameterAsText(1)
# Output folder
  outFolder = arcpy.GetParameterAsText(2)

# Get the spatial reference
 spatialRef = arcpy.Describe(template).SpatialReference

 #Loop through shapfiles in folder and reproject
 for fc in fcList:
    fcspatialRef = arcpy.Describe(fc).SpatialReference
    if fcspatialRef != spatialRef:
        arcpy.Project_management(fc, outFolder + "\\" + fc, spatialRef)

# Return any errors
except:
    arcpy.AddMessage(arcpy.GetMessages ())

# Print shapefile list
for fc in fcList:
    print fc


texts: SAGA GIS getting 90 degree slope in the whole terrain

texts:
I´m studying a jungle area in Peru, and when I try to calculate the slope with the "Slope, Aspect, Curvature" tool from SAGA v7.7.0, I get that the whole terrain has a 90-degree slope. In the contrary, when I use the ArcGIS slope tool, I get normal values. My DEM is in meters, same units for the x, y, and z axis. Currently trying with the 9 parameter 2nd order polynom (Zevenbergen & Thorne 1987) method in SAGA.
How can I do this?


texts: How to determine left and right side of crossing lines?

texts:
Given two lines of which we know the starting XY points and intersecting point, I need to find which line is left and which line is right. The lines are always crossing and can be turned in any direction. Here are 4 possible scenarios as an image where A is always on the left line:

Below are the same four scenarios as data, where every scenario has two crossing lines. Each line has a start point (x, y) and an end (x, y) point. In the scenarios below, the starting point of A is always on the left of the starting point B and therefore lineA is on the left of lineB:
line_A1 = ((1, 1), (0, 0))
line_B1 = ((0, 1), (1, 0))

line_A2 = ((0, 1), (1, 0))
line_B2 = ((0, 0), (1, 1))

line_A3 = ((0, 0), (1, 1))
line_B3 = ((1, 0), (0, 1))

line_A4 = ((1, 0), (0, 1))
line_B4 = ((1, 1), (0, 0))

middlepoint_C = (0.5, 0.5) #The middle point in all 4 scenarios is (0.5, 0.5)

How can we determine left and right side of the crossing lines in the four scenarios above?
Preference for solution would be with built in Python modules or shapely, but all solutions are welcome since we are kinda stuck.


texts: Showing labels as they will be printed/exported in my print layout?

texts:
I set up a print layout for a map I want to print in QGIS. Some of my layers have labels on them. When I zoom in and out in the screen, without changing anything else on the print layout, some labels shift position relative to the page and other features. (Sometimes they even appear or disappear.) I am not zooming inside the map, I am zooming the view of the entire layout.

The thing is, when I print the layout to paper, the labels won't move when I move my eyes closer or further from the page. I'd like to see in the QGIS print layout what the labels will look like when printed. How can I do this? I'm fine with QGIS choosing where to place the labels. I just want it to be consistent in the Print Layout with what will actually be printed.

I couldn't find any options that sounded like they would change this on QGIS 3.10.1. I tried locking the layers and/or the map item, and that only did what it'd be expected to do; it did not change the behavior of the labels. I looked around for a "symbology scale" option and found none. Am I missing something?


texts: Unable to open hgt downloaded from dds.cr.usgs.gov

texts: I downloaded a SRTM hgt file from https://dds.cr.usgs.gov/srtm/version1/SRTM30/w100n40/w100n40.dem.zip, after unzipping I am not able to open it in QGIS or even gdalinfo. I downloaded other hgt files from this other site http://dwtkns.com/srtm30m/ and can open them without trouble; I read the file names have to do with the CRS of the file, so I played a little with the one that doesn't open, trying to follow the names of the one that does open, but not even this way I could read it; the one that opens starts with latitude (N19W097.hgt) whereas the one that does not open starts with the longitude (w100n40.hgt).


texts:
The question is related to interpretation of grids representing soil types separately, like salty soils:

TAXNWRB_Gleyic.Solonetz

TAXNWRB_Calcic.Solonetz

TAXNWRB_Haplic.Solonetz

TAXNWRB_Mollic.Solonetz

TAXNWRB_Gypsic.Solonchaks

TAXNWRB_Haplic.Solonchaks..Sodic

TAXNWRB_Haplic.Solonchaks

It is written that the units are "Predicted probability in percent" but the pixel values are in range between 0-10 (in one case even 12?).

Are those values representing percentage divided by 10? Or it is something else behind them?


texts: Exporting Raster Layer using "Export Raster Data" vs Copy Raster

texts:
Is it possible to access the "Export Raster Data" tool in arcpy?

Some background on the steps I am trying to process:
I have a raster loaded into an mxd.  I have converted the raster to a temporary layer ("make raster layer" tool) file so that tabular data (in the form of a .csv file) can be joined (via Add Join tool) based on a common field. I then export the raster layer with the joined tabular data to a new raster.  My issue is with the last step: export to new raster.

The "Export Raster Data" tool that is accessed by right clicking on the raster layer to bring up the Data/Export Data menus produces a new raster in about 30mins.  The file is exported to a file geodatabase with the same spatial reference and non compression type.

The "Copy Raster" tool provides the same output as "Export Raster Data" tool but takes 14hours to run on the same dataset.  The file is exported to a file geodatabase and optional values are left blank.

I have hundreds of rasters that I need to export so I have been building a python script to do the process for me (convert to raster layer, join tabular data, and export to new raster). But since the "Copy Raster" tool is taking too long to run, I want to find a way to use the "Export Raster Data" tool in the script.

I have also tried the "Copy" tool and the "Raster to Other Format Conversion" tool but they do not handle raster layer files.

I am running ArcGIS 10.2 Basic with Spatial Analyst Extension. Below are the python snippets extracted from ArcMap once the tools were run in the mxd.  The last tool (copy raster) took 14hrs to run in arcmap so I'm looking for an alternative (such as the export raster data tool).

print "     make raster layer"
arcpy.MakeRasterLayer_management("D:/Alces_Projects/AO_Alberta_NEW/Data/Indicators/Climate/ClimateChange_Processing_CAAB_20160411.gdb/blank_canvas_Tier3Tile5_v1f_20160411","blank_canvas_Tier3Tile5_v1f_20160411_layer","#","D:/Alces_Projects/AO_Alberta_NEW/Data/Indicators/Climate/ClimateChange_Processing_CAAB_20160411.gdb/blank_canvas_Tier3Tile5_v1f_20160411","#")
print "     add join to csv data"
arcpy.AddJoin_management("blank_canvas_Tier3Tile5_v1f_20160411_layer","Value","D:/Alces_Projects/AO_Alberta_NEW/Data/Indicators/Climate/ClimateChange_Software_Output_Files/DecadeScenarios/Tier3Tile5_Decade_1921_1930.csv","ID1","KEEP_ALL")
print "     copy joined raster to new raster"
arcpy.CopyRaster_management("blank_canvas_Tier3Tile5_v1f_20160411_layer","Tier3Tile5_Decade_1921_1930","#","#","2147483647","NONE","NONE","#","NONE","NONE")

I thought I would edit this post to include the solution that I found to the problem.  Turns out the issue was with the join to the csv file.  I converted the csv to DBASE file and the copy raster process ran in minutes instead of hours.


texts:
I'm trying to get Mapbender3 to work with Apache2.4 and PHP5.6. I have Apache and PHP talking to each other and installed correctly. I have followed the instruction to install Mapbender3 found here:
https://doc.mapbender3.org/en/book/installation/installation_windows.html

The purpose is to set up a server in order to test Geoserver and Mapbender3's abilities to talk to each other and produce a webmap.

NEW ISSUE: I go to the symphony config page. http://localhost/mapbender3/config.php I am greeted with the text below. I make the changes to the php.ini in the Apache folder and re-check the configuration. The issues don't disappear. Do I need a php.ini in the Mapbender folder as well?

This script will guide you through the basic configuration of your project. You can also do the same by editing the ‘app/config/parameters.yml’ file directly.

MAJOR PROBLEMS
Major problems have been detected and must be fixed before continuing:
date.timezone setting must be set
Set the "date.timezone" setting in php.ini* (like Europe/Paris).

RECOMMENDATIONS
Additionally, to enhance your Symfony experience, it’s recommended that you fix the following:

intl extension should be available Install and enable the intl extension (used for validators).

a PHP accelerator should be installed Install and/or enable a PHP accelerator (highly recommended).

realpath_cache_size should be at least 5M in php.ini Setting "realpath_cache_size" to e.g. "5242880" or "5M" in php.ini* may improve performance on Windows significantly in some cases.

short_open_tag should be disabled in php.ini Set short_open_tag to off in php.ini*.

PDO should have some drivers installed (currently available: none) Install PDO drivers (mandatory for Doctrine).

* To change settings, create a "php.ini".

When I go to my Apache mapbender page http://localhost/mapbender3/app.php I get a fatal error. How do I fix this?

Fatal error:
Uncaught exception 'Exception' with message 'DateTime::__construct(): It is not safe to rely on the system's timezone settings.
You are *required* to use the date.timezone setting or the date_default_timezone_set() function.
In case you used any of those methods and you are still getting this warning, you most likely misspelled the timezone identifier.
We selected the timezone 'UTC' for now, but please set date.timezone to select your timezone.'
in C:\Mapbender3\mapbender\src\Mapbender\CoreBundle\DependencyInjection\MapbenderCoreExtension.php:30

Stack trace:
#0 C:\Mapbender3\mapbender\src\Mapbender\CoreBundle\DependencyInjection\MapbenderCoreExtension.php(30): DateTime->__construct('now')
#1 C:\Mapbender3\vendor\symfony\symfony\src\Symfony\Component\DependencyInjection\Compiler\MergeExtensionConfigurationPass.php(59): Mapbender\CoreBundle\DependencyInjection\MapbenderCoreExtension->load(Array, Object(Symfony\Component\DependencyInjection\ContainerBuilder))
#2 C:\Mapbender3\vendor\symfony\symfony\src\Symfony\C in C:\Mapbender3\mapbender\src\Mapbender\CoreBundle\DependencyInjection\MapbenderCoreExtension.php on line 30


texts:
In my project, I need to calculate each village's distance along a
    route to a seaport.

    However, as you can see in the picture(now deleted for space), not all
    the villages are connected to roads already, so, first of all, I want
    to ask how to create the shortest path from each village to the road.
    Since there are more than 400 such villages, I'm looking for an
    automatic way to do this job.

    My second question is how to finding the shortest path and distance
    base on the network. I'm looking for an automatic way like running
    certain python functions to do this too.

  I have done the first one by using the "Near" tool in "Analysis" and
  the "XY to Line." "Near" generates the closest coordinates on the road
  and "XY to Line" can draw lines from villages to the nearest
  coordinates.

  I'm still frustrated by the second part.

  I have tried to use the way suggested by user3338197 in the first
  answer. I merged existing roads and the roads I created into one
  shapefile and built a network based on that and remained the default
  setting unchanged (I'm using 10.3). However, it always shows

    Warning: No "Destinations" found for "Location 1" in "Origins"

  or

    Warning: No "Facilities" found for "Location 1" in "Incidents"

  I tried to figure out what is the problem, so I ran a "New Route"
  problem, and I got a result like this:

  [![screenshot2][2]][2]

  [1]: http://i.stack.imgur.com/mg1uF.jpg   [2]:
  http://i.stack.imgur.com/XTl7y.png

  As you can see, the first and the second location are connected by
  roads, but the network analysis failed to generate a route between
  them. The most important thing I'm trying to get is the distance from
  each village to the port.

  Hope someone can help me fix this problem.

In the end, it's because some junctions are not connected as they seems to be. Solved by "Integrate."


texts:
So v.out.svg works perfectly fine for my purpose, except that it doesn't export polygon IDs, so I have no way to identify them in the output. I want to use one of the attributes (specifically the Name attribute) as the ID of each polygon.

This dialog even allows me to specify "Attribute(s) to include in output SVG", but selecting my attribute in this option doesn't change the output in any way.

I need the attribute to be in SVG content (as an XML attribute of the corresponding  element.

Is there another way of doing what I'm after?


texts:
I have polygons represented in Leaflet.
Then I added some markers over these Polygons.

var = customMarker = L.Marker.extend({
    options:{
        Layer: ''
    }});

I extended marker to store the layer ID of my polygon.
How can I know which what is the Layer ID of my Marker?


texts: Clarification for WMS Service in MapServer

texts: I am trying to understand the WMS service in MapServer. What I've done so far is making a complete map with a roads network some buildings and parks. Then I create an HTML template and by calling the url I get my map. My question is: is this map considered a WMS service?  This confuses me because I came across this http://mapserver.org/ogc/wms_server.html that explains how to set up a WMS server and so I set it up and now I can call fot the GetCapabilities and GetMap requests. So, this confuses me. Various resources on the internet refer to a simple map as a WMS. On the other hand, the link refers to the WMS as something a bit more different. Eventually what is the difference?


texts: ogr2ogr CSV conversion with two points into column names

texts:
I have a CSV with ":" in all column names. I am trying to convert it with ogr2ogr to another format. I don't need all fields

example for formating:

#0:imgID,1:camSN,2:lat,3:lon,...,57:fields
000000001,201002422,4.8192835254,-75.7043506327,...
000000003,201002422,4.8192743299,-75.7042594643,...
000000004,201002422,4.8192692652,-75.7042137224,...

where I use:

ogrinfo Outfile.csv -dialect SQLite -sql "SELECT substr(1:camSN , 2) as serial FROM Outfile"

I get:

INFO: Open of `Outfile.csv'
      using driver `CSV' successful.
ERROR 1: In ExecuteSQL(): sqlite3_prepare_v2(SELECT substr(1:camSN , 2) as serial FROM Outfile):
  near ":camSN": syntax error

Where I use:

ogrinfo Outfile.csv -dialect SQLite -sql "SELECT 1:camSN,2:lat,3:lon FROM Outfile"

I get:

INFO: Open of `Outfile.csv'
      using driver `CSV' successful.
ERROR 1: In ExecuteSQL(): sqlite3_prepare_v2(SELECT 1:camSN,2:lat,3:lon FROM Outfile):
  near ":camSN": syntax error

where I use:

ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:4326 -oo X_POSSIBLE_NAMES=lon* -oo Y_POSSIBLE_NAMES=lat* -f "ESRI Shapefile" Informe.shp Outfile.csv -dialect SQLite -sql "SELECT camSN,lat,lon,height,roll FROM Outfile"

I get only Informe.dbf and Informe.prj NOT shapefile (Informe.shp)


texts: Does anyone know how to export a style file like QGIS qml or SLD format from Global Mapper? Gm is very good at importing files like .tab with styles preserved, but I cannot find a way to write these out in a format QGIS can understand.


texts: Reproducing raster color rendering from ArcMap in QGIS?

texts:
I have a 16-bit satellite geotif. These images can be visualised in a number of ways, depending on which stretch you apply.

The default stretch applied by ArcGIS to these images produces a "natural" result, which is good.

However, when I view a raster (satellite geotif file) in QGIS, the color is very different than in ArcMap. The result is over-exposed and less "natural".

Can I maintain/replicate ArcGIS's color rendering / stretch for users viewing the satellite tifs in QGIS?


texts: Trying to rescale a raster with color table but resampling goes wrong

texts:
I'm trying to rescale a palletized raster using QGIS (2.14.21), and GDAL (2.2.3), using a bilinear resampling filter, but the resulting colors in the new raster are weird...

Here is the command line:

gdalwarp -r bilinear -tr 0.00006 0.00004 -dstnodata 0 -of GTiff C:/source.tif C:/out.tif

The source raster (source.tif):

The resulting raster (out.tif):

Any idea what could cause this?

Nota:

I don't have this problem with the "near" filter. But it is the only filter working properly in this case.
The magenta line in the bottom is the no-data value.