texts:
The "Styler" plugin for Geoserver promised flexible styling of your data, using the GeoExt library:

The user guide tells you to first install the REST plugin, and then extract the styler zip file into the /data/www Geoserver directory. However, the download page for Geoserver 2.2 does not include the REST plugin. Installing older versions (I'm using Geoserver 2.2) failed for me, with one error saying "This web site needs a different Google Maps API key" (from /geoserver/www/styler/ - I have a Google Maps API).

REST seems to have disappeared recently: versions up to 2.0.3 contained the option to download the REST zip package, whereas after that there is no mention of REST, and the installation guide remains.

I did find some evidence that REST may have died, however:

  This feature is no longer updated as an independent GeoServer download, but is available as part of the OpenGeo Suite. For more information on its development see OpenGeo project page

If Geoserver's "styler" plugin really is dead, then I should probably change the question to something more general like "how do I make OpenGeo's Styler work on Geoserver?" Just wanted to check with others who've had similar experience.

Context
I've set-up Geoserver on an AWS EC2 instance (testing map can be seen here).
What I'm trying to do is to make those boring grey polygons colourful and interesting, based on their attributes. I've tried exporting the sld that discribes the style below into Geoserver, but, at present, this does not work.


texts: Finding differences between shapefiles using ArcGIS for Desktop?

texts:
Every month I have to QA/QC some streets data for an entire county (over 47,000 records). The first step is to compare the two files. Currently, I'm using a tool that was created in house via VBA approx 5 years ago. It works compares the two files and reports back the differences based on a few factors (added records, deleted records, address changes, city changes, segment length/vertices changes). I then export these to individual shapefiles. Additionally, the tool works in ArcGIS 9.3 but will no longer work in ArcGIS 10.1

The tool works great, but it takes approximately an hour to complete.

I've started to work on hopefully improving this tool to speed up the performance. I've looked into a variety of methods including the Compare Feature Tool and Select by Location.

Both of these tools work significantly faster, however, they don't seem to capture everything that I need. The big issue I have with the Compare Feature Tool is that I cannot choose which field to compare on. It defaults to ObjectID when I need to compare by Segment_ID. There isn't an option to select that.

Does anybody have a good solution or ideas on how to achieve a full compare like I've described?

For some reason, ArcGIS isn't recognizing the join that I'm attempting per RyanDalton's answer.

But

Taking the advice of a join I have figured out how to find the added records and deleted records differences in the two files by:

Joining the Previous Month shapefile to the Current Month shapefile via the Segment_ID field
Performing a definition query where Segment_ID IS NULL
Exporting that to a new shapefile (Added Features)

I then reverse the joins and find the Deleted Features

Still working on a way to find the other differences listed and then string them all together in a model.

I can find all records that have changed by performing a select by location records from one of the datasets that are identical to the other dataset. I then switch the selection and it gives me all the records that are not identical.

I'll put it all in a model and hopefully it works just as well.


texts: Predicting the whole image through trained classifier and visualize it as classified landuse and landcover image(LULC)

texts:
I am working on LULC classification of sentinel-2 image(4-bands) using Machine Learning algorithms. I split the whole image into 64*64 patches and assigned class label (builtup, barren land, water and vegetation) to selected patches. For each class patch spectral features (mean, variance values of each band, mean normalized difference vegetation index(NDVI) for each class patch, texture features (grey-level co-occurrence matrix (GLCM) contrast, homogeneity, dissimilarity and correlation for NIR band of each class patch were calculated. I split the features into training and testing dataset and trained random forest classifier on training dataset. The problem is I don't know how to classify the whole image using trained classifier and visualize it as a classified image.
'''
import numpy as np
import pandas as pd
from sklearn.metrics import accuracy_score
from sklearn.ensemble import RandomForestClassifier
from sklearn.metrics import accuracy_score
from sklearn.model_selection import StratifiedKFold
from sklearn.model_selection import cross_val_score
from sklearn.metrics import confusion_matrix
train = pd.read_excel(r'D:\band_composite\traning_data\training.xls')#import traing dataset
test=pd.read_excel(r"D:\band_composite\test_data\test.xls")#import testing daatset
# (Xtrain,Ytrain for training dataset)
X_train = train.drop(['class'], axis=1)
y_train = pd.DataFrame(train['class'].values)
X_test = test.drop(['class'], axis=1)
y_test = test['class']
# Random forest model training
clf = RandomForestClassifier(max_depth=2, random_state=0)
model=clf.fit(X_train, y_train)

'''
# X_train shape
    (141, 15)
# y_train.shape
    (141, 1)

#X_training dataset excel image
  [ist five rows of X_train dataset][1]

  [1]: https://i.stack.imgur.com/uYjTy.png


texts: cumulative atlas in QGIS 2.14 and above?

texts:
I have a question about something which I think is possible but just can't get my head round it, and haven't found any answers online. I'm using $id= @atlas_featureid in QGIS 2.14 to create an set of images from the Atlas tool, and it works fine as usual. However, instead of having each map iterate through the dataset and show one row of data at a time, I'd like it to show one row of data and then add each row of data successively in a cumulative way. I think I can use a rule to achieve this, but not sure how.

Example image: https://2.bp.blogspot.com/-326CYkmBxH0/WOUFwYzYv6I/AAAAAAAACIk/bzgu-84DcX8QR23zFItNF2lZXnUGE7SqQCLcB/s1600/gun_homicides_cumulative_2015_100ms.gif

Individual Atlas extract - one day = one map.

To help explain, see this gif I created from a series of maps on US gun crime. Instead of one day's worth of data appearing at a time, I'd like it to build up over time so that on the Jan 31 frame there is 31 day's worth of dots on screen and then it's just a mess of dots by Dec 31 - I know it will probably look bad but I'd like to use it in other projects as well.

Does anyone know how this could be done?


texts: How to select polygons by distance to a single polygon in another layer?

texts:
I have a shapefile A with a group of polygons, lakes.
Now I want to select the polygons(lakes) that are within a specific distance to a polygon in a second shapefile, B, with the polygon a city.

In other words, I am looking for something like ArcGIS "Select By Location" and "are within a distance of the source layer feature"

I tried Select by Location, and also Select by Radius, but didn't manage to get any of them to work.


texts: I have a file of Tennessee zip codes and another excel file including 600,000 origin and destination addresses with zip codes and trip purposes. Since I cannot Geocode all 600,000 addresses, I Geocoded zip codes, and now want to create a flow map with joining/relating these two files together. Is it possible to do so?


texts: Create a Select by Attributes subquery from a locational selection in ArcGIS with multiple variables

texts: I have a large shapefile of GPS points. Each point has a unique point ID number, but they are all part of different clusters, where each cluster has a unique cluster ID number (there is no overlap between clusters). For example, points 1, 2, and 3 are part of Cluster 100 and points 4, 5, and 6 are part of Cluster 200. I have done a Select by Location query of these points based on their proximity to a road feature class, and now have about 1400 points selected. I would now like to select any additional points that share a common Cluster ID to the already selected points. For example, say I captured points 1 and 2 of Cluster 100 with the locational selection query, how would I capture point 3 of Cluster 100 as well, in an efficient query? Or is there a way to have a subquery in the Select by Location query that would look at the cluster ID attribute?


texts: Mimicking ArcGIS Combine function in PostGIS

texts:
I am trying to use PostGIS to create something similar to the attribute table produced when two raster files are combined using the ArcGIS Spatial Analyst function "Combine". This shows the count of unique combinations of values associated with each pixel across each raster.

I have a functioning query for a very small (8 pixel by 8 pixel) multiband raster which replicates this function where merged_rasters is my PostGIS raster table with 2 bands.

SELECT b1val, b2val, COUNT(*)
FROM(
    SELECT ST_Value(rast, 1, x, y) As b1val, ST_Value(rast, 2, x, y) As b2val
    FROM merged_rasters CROSS JOIN
    generate_series(1, 10) As x CROSS JOIN generate_series(1, 10) As y
    WHERE x <= ST_Width(rast) AND y <= ST_Height(rast)
    ) as merged_values
GROUP BY b1val, b2val;

However, when I try using the ST_Value on larger rasters, even tiled rasters with a tile size of 256x256, this becomes very slow.

How would I improve performance?


texts:
I am very new to GIS and am trying to use QGIS 2.8 to make a very large display map of Southern Ontario, roughly 4 feet by 8 feet. I would like to use aerial imagery, and am thinking of using the Bing aerial base map for that. The problem is, I would need to make an image file which when blown up to that size, doesn't look blurry.

I've vaguely heard of georeferencing being used to "stitch together" high-resolution images of small areas, which I thought might be helpful.

Are there any other efficient ways I could go about this?


texts:
Possible Duplicate:
  How to add a simple WFS layer from GeoServer to OpenLayers map?

I'm trying to simply add a WFS layer in OpenLayers. But completely stuck! Could someone please take a look at my code and suggest what I'm doing wrong?

I have a fresh install of GeoServer on a Linux (Ubuntu) VM. I'm trying to connect to one of the default layers (medford:hospitals).

I'm fairly confident I have all the parameters correct. The layer previews fine in GeoServer. My webserver is on a seperate server to GeoServer so hopefully no proxy issues?

Starting to wonder if this is some weird server issue?

Anyhow - the code:

            //Add WFS Layer
        var wfs = new OpenLayers.Layer.Vector("Medford Hospitals",
                {
                    strategies: [new OpenLayers.Strategy.Fixed()],
                    projection: new OpenLayers.Projection("EPSG:4326"),
                    protocol: new OpenLayers.Protocol.WFS({
                        version: "1.1.0",
                        url: "http://squirrel.bis.local:8080/geoserver/wfs",
                        featurePrefix: 'medford', //geoserver worspace name
                        featureType: "medford:hospitals", //geoserver Layer Name
                        featureNS: "http://medford.opengeo.org", // Edit Workspace Namespace URI
                        geometryName: "the_geom" // field in Feature Type details with type "Geometry"
                    })
                }
        );
        map.addLayer(wfs);

Any help and advice gratefully accepted!

Thanks,
Steve


texts:
I created a route layer that uses hatches to help determine addresses along a street centerline. About 90% of the time, the hatches go up from zero in the direction that the line was digitized but about 10% of the time (+/-) the hatches go in (what I would say is) the "wrong direction". I can symbolize this layer with arrows showing the direction of the line which shows me the instances where the hatches are wrong.

What determines the direction of the hatches and can I do a better job of controlling this?

Wow hatches or hatching or centerlines are new tags ? Wish me luck with this one!


texts: How do I embed a KML map using the "new" Google Maps?

texts:
I have a KML link that I can use to generate an embedded Google Map using "old" Google Maps and would like to use it instead to create a "new" Google Map embedding.

To generate the an "old" embedding, all I need to do is click on the link icon on the page for my KML, which generates an iframe I can use for embedding:

<iframe width="425" height="350" frameborder="0" scrolling="no" marginheight="0" marginwidth="0" src="https://maps.google.com/?q=http:%2F%2Fshare.abvio.com%2F9d18%2Fde58%2F4d38%2Fa340%2FCyclemeter-Cycle-20140629-0905.kml&amp;ie=UTF8&amp;t=m&amp;ll=42.429738,-71.31901&amp;spn=0.135997,0.21338&amp;output=embed"></iframe>

But this is an "old style" embedding, and I would like to generate a "new style" one from the same KML. Is there a simple mapping from one to the other or from my KML URL directly to the form of the "new iframe?


texts: Getting points-touching-a-line attributes back to lines layer using QGIS?

texts:
I'm using QGIS 2.14.4-Essen. I have two layers:

points.shp that contains points with a YEAR column
lines.shp that contains lines that exactly connect the points from points.shp

I would like to get the YEAR attribute from points.shp back to lines.shp. Each line is a single segment with one point at each end (see picture below). What I would like to get is YEAR from first point and YEAR from second point back into each line attributes.

For instance: line 1 is touching a first point with YEAR=2010 and a second point with YEAR=2011. I would like to get back something like '2010-2011' into line 1 attributes. The result should look like this:

id  | points
----+-----------------
1   | 2010-2011
2   | 2011-2012
3   | 2012-2016

Is there a way to achieve this using QGIS?


texts: Alternative to Rasterio Mask module for clipping a raster using masked geojson

texts:
I am clipping a raster (.tif file) using a mask vector (.geojson) file.

with rasterio.open(url) as src:
    out_image, out_transform = rasterio.mask.mask(src, shapes, crop=True)

But, this rasterio.mask module is taking some time to create this NumPy array out_image.

I want to create the clipped image as fast as possible.

Is there any alternative way (for eg. from rasterio.io import MemoryFile module) to clip a raster with a mask vector layer using Rasterio (only using Rasterio and not using any other libraries) ?


texts: Unable to get processing extent to match desired raster output

texts:
I have a few thousand ortho mosaic tiles that I've resampled to a lower GSD and am trying to cut into larger tiles. I developed a script that has worked perfectly for thousands and thousands of tiles over the course of many projects.

I am having a problem now where at the upper right corner of my current project area the tiles are not cutting into the proper image extent, but only for two of the five tiles I'm attempting to cut.

I'm using ArcGIS 10.2 for the clipping process. The script utilizes arcpy.Clip_management. The images are TIFFs. I've played with dozens of environment settings, including processing extent, and I get the same thing every time. The input data has RGB values of 0,0,0 outside of the bad extent it is clipping to. The script itself clips to a shapefile defining the area of interest, but has always matched the input extent regardless of whether or not there was actually image there.

The process clips to these bounds:

When I need it to clip to these bounds (and yes, I need/want all of those black pixels so image size will match all other tiles):


texts: What is the best database for an online Android location app?

texts:
I am working on Android app for locating food places in my city.
Using Google Map v2 api.

App shows user his current location, and the closest food spots around him.

I want to store location coordinates in an online database, and App should retrieve data
in JSON format from db to phone.
I am familiar with PHP and MySql.

Which database should i use to store locations ?
Is MySql good for that purpose ?
How will i store location data in it?


texts: Deleting SQL Server database in order to delete ArcSDE geodatabase?

texts:
I'm new to this ArcSDE stuff and would like to know the best way to cleanly delete an SQL Server database in order to delete the ArcSDE geodatabase. Is it as easy as right clicking on the database in SQL Server and hitting delete? Or is there a more preferred way to do this?

Also, is there a way to rename the SQL server database that the ArcSDE geodatabase is using? Or is it best to just start fresh with a new database?


texts: How to make a topographic map with a Digital Terrain Model shapefile?

texts:
I'm 2 weeks in self-learning QGIS so I don't know if my question is the same with the suggested similar ones.

I'm making a map like this. A map of a mountain with a "mountain-like" texture.

By u/coolrivers from r/qgis

I have a Digital Terrain Model shapefile of our region and how do I have a texture like that?

I tried the Batch Hillshader and its results into error everytime I run it.

This is what the shapefile looks like.

Edit. It's a digital terrain model, not a digital elevation model


texts: Combine unbounded multiple point features into 1 point feature

texts:
Is there any way that I can combine multiple point features into 1 point feature?

I tried 'Union' and 'merge', but neither of them works?


texts: ArcMap 10.5 Network Analyst Routing Unchecked Start Time

texts:
I'm doing a network analysis with ESRI's ArcMap 10.5.  As part of this analysis I'm being asked to compare routes at different days and times to determine the impact of traffic congestion on routing.  I'm curious, for a Route analysis if you leave the "Use Start Time" radio button unchecked under the "Analysis Settings" tab how does Network Analysis calculate travel time?  Is it simply using free flow, i.e. posted speed limit, conditions?  Is it using an average of travel times throughout the week?  I know my network data set includes traffic information and I get different results when I pick different days and times.  In fact, the results are exactly what you'd expect.  If you pick rush hour travel time is worse than on the weekends.  However, I have no idea what the method is when no start time is chosen.

I can't seem to find a definitive answer in the ESRI documentation.


texts: Making routable graph from Point and Line shapefile using Networkx?

texts:
I have one shapefile which consists of all the streets in the study area, and another shapefile which has points representing origins and destinations.

Is it possible to find shortest path between all origins and destinations, using the links in the link shapefile?

I primarily use Python, Networkx and QGIS for my work, so it would be great to use these tools for this task.


texts: Why is my CRS not projecting data points where they are meant to be using QGIS?

texts:
I am trying to plot pipelines which are expressed as multiple data points. For each data point there is a Northing and Easting which have been projected using the ED50 30N projection.

I need my data to be in WGS84 EPSG:4326. I have loaded a new QGIS with the basemap projection of WGS84. I have then uploaded the data as a delimited text file. When uploaded, I thought the points were not appearing. However, from closer invesigation they are just being projected far, far off the base map as seen below (bear in mind 'Peterhead' is actually situated within the world map, as seen in the 3rd figure ).

The figure above identifies two seperate pipelines which have been input. However the scale below is 0:1 (Which is as large as QGIS can go).

This figure can show how large the map is, as peterhead is situated within the world map.

Can anyone solve this issue. I have taken the first Easting and Northing idicated within the first figure and tested it within the map. It projects the data to the west coast of scotland seen in the bottom figure (nevertheless this is still wrong as it should be on the east coast of scotland). However, even if my coodinates are wrong my projections shouldnt be as far away as they are. does anyone know why?

I have tried exporting feature and saving it as a new shapefile with the WGS 84 projection, but it still occurs in the same place.

However, even if my coodinates are wrong my projections shouldnt be as far away as they are. does anyone know why?


texts: Google Earth Engine: Get a Geometry object based on user's drawing

texts:
In order to publish my code as a Google Earth Engine App, I need to interactively set the Area of Interest (AOI) of the analysis  based on the users's drawing (Rectangle in this case). I'm trying to use the .onDraw() method as described in the ui.Map.DrawingTools and ui.Map.GeometryLayer documentation examples but I cannot seem to make it work. I'm trying the following:

// Don't make imports that correspond to the drawn rectangles.
Map.drawingTools().setLinked(false);
// Limit the draw modes to rectangles.
Map.drawingTools().setDrawModes(['rectangle']);
// Add an empty layer to hold the drawn rectangle.
Map.drawingTools().addLayer([]);
// Set the geometry type to be rectangle.
Map.drawingTools().setShape('rectangle');
// Enter drawing mode.
var getAOI = Map.drawingTools().draw();

var getAOI = function(){
  Map.drawingTools().layers().get(0).toGeometry();
};

var AOI = Map.drawingTools().onDraw(getAOI);

Here is a link to the full script: https://code.earthengine.google.com/11a77302259a6d99c3a561a2a723fca0


texts: Get attribute information from a polygon feature using command line tool

texts:
I am creating some type of command line tool that will take a coordinate point(X and Y) and return an attribute from a polygon (shapefile) that surrounds the coordinate.  I found this

To get attribute data of polygon layer by giving Lat/Lon

but it deals with a SQL statement and does not return any information.

I am using shapefile but I can convert it to another format.

I havent tried anything but I was hoping someone would give me a starting point.


texts:
I'm trying to use GDAL polygonize to run a batch convert but I cannot for the life of my figure out a way to get it to read a text file as input that contains all of the commands.

Do you know the correct command to run this conversion in one step so I don't have to do it many times over?

Here's an example of a single line of my input txt file:

gdal_polygonize X:\Source-Imagery\Satellite\WaterMask\Collection_5\2000\ALL\MOD44W_Water_2000_ML1920.tif\MOD44W_Water_2000_ML1920.tif -f "ESRI Shapefile" X:\Source-Imagery\Satellite\WaterMask\Collection_5\2000\ALL\MOD44W_Water_2000_ML1920.tif\MOD44W_Water_2000_ML1920.shp

The txt file contains 176 lines to convert.


texts: after to make join between 2 tables a field is not recognize

texts: I made join between two attribute tables with a common field and then I tried to use python to read a field and now field is not recognized on the new table.. Is it normal? How can I solve this problem?


texts: Post-processing GPS Satellite Measurements

texts:
I am working on a system that reports its location to a server using very little power. In order to achieve an extremely low power usage, I am avoiding having to download both the almanac and the ephemeris data from the satellites. Downloading these from another source does not wield an acceptable power usage either. Instead, I am trying to post-process the on the server.

I came across the following question with a similar situation, where the user provided an answer himself:

  Storing only the satellite ID, its timestamp and the receiver's timestamp (when the Navigation Message arrived) is enough to calculate position. As hinted in Kirk Kuykendall's comment, this is equivalent to the satellite ID and pseudorange data stored in a RINEX file.

Then I found the following answer to another question about post-processing:

  To perform post-processing you typically need a GPS that can record raw pseudorange and carrier cycle data. Very few consumer level GPS support this functionality, and those that do require software modification.

I believe that these modules are pretty common nowadays. I am planning on using a Maestro A2235-H GPS Module, which is using the SiRFstarIV architecture. The GPS module supports both the NMEA protocol and the SiRF Binary Protocol (with the One Socket Protocol extension).

According to this One Socket Protocol Interface Control Document, the One Socket Protocol contains a Navigation Library Measurement Data message (Message ID 28, page 287). This message contains what I believe to be the necessary data to reconstruct the receiver's accurate position on the server. Each measurement, it reports:

Time Tag
Satellite ID
Estimated Satellite Time (GPS Software Time)
Raw Pseudorange (Without ionospheric, tropospheric or clock corrections)
Carrier Frequency
Carrier Phase

Now I'm looking for software that can post-process the receiver's GPS position using these satellite measurements. I collected the following libraries from the two questions above:

RTKLIB
GeoKettle
GPSTk

Unfortunately, this is where I got lost. I'm not sure which library supports my needs. I simply need a tool that converts the satellite measurements into a lon/lat position which I can incorporate into my code on the server. I have therefore the following questions:

I could generate a RINEX file for every measurement, but do I have to feed the RINEX file to these libraries? Or can I feed the satellite measurement data directly?

Do these libraries collect the required ephemeris data automatically?

Which library do you recommend for my situation?


texts: Check if point in GIS polygon fails to produce correct output

texts:
I've tried many, many many, ways to find out if a point is in a GIS polygon based on many great posts, but I can't get any to work properly.  Specifically I created a grid of points covering the bounding box of Japan, and now I want to keep only the ones that are actually over Japan.

Here is one way I tried to get a GIS shapefile for Japan:

world = gp.read_file(gp.datasets.get_path('naturalearth_lowres'))
JapanShapes =  world[world.name == "Japan"]
geoJapanShapes = gp.GeoDataFrame(JapanShapes2)
geoJapanShape = geoJapanShapes['geometry'].unary_union

I have millions of point data in a geodataframe, but here's just one point (Tokyo station).

allCoords = pd.DataFrame([[35.6812,139.7671]], columns=["lat","lon"])
allCoords['geometry'] = allCoords.apply(lambda row: "POINT ("+str(row.lat)+" "+str(row.lon)+")", axis=1)
allCoords = gp.GeoDataFrame(allCoords)
allCoords['geometry'] = allCoords['geometry'].apply(wkt.loads)

Now, if test whether Tokyo station is in Japan:

print(allCoords['geometry'][0].within(geoJapanShape))
>>>> False

However, the last time I checked Tokyo station was within Japan.  And if I do:

coordsWithinJapan = allCoords[allCoords.apply(lambda row: row['geometry'].within(geoJapanShapes), axis=1)]
>>>> AttributeError: ("'GeoDataFrame' object has no attribute '_geom'", 'occurred at index 0')

Based on several other posts I also tried:

from geopandas.tools import sjoin
gridWithinJapan = sjoin(geoHexData, geoJapanShapes, how='left', op='within')

Which just returns every point in the original data.  If I try with the inner option I instead get an empty dataframe.

So both methods (kinda) run on my data, but neither method is working to accurately report whether a point is within the boundary. Other versions gave me crs errors, but I don't how to deal with that considering my points are hand-made.  Another had a problem because maybe/apparently that shapefile is a multipolygon, which it couldn't handle.  No matter what I try, I seem to be doing something wrong, but I can't find the error.

  How do I actually determine whether my points are in my shape?


texts: Creating new variable conditional on another feature using Python Parser of ArcMap Field Calculator

texts:
I don't usually work in ArcGIS Desktop or Python.
I have a shapefile of road data and the roads have a 'class' feature. I want to categorize the 'Class' feature into just 5 groups and call that variable 'Type'. It's just a basic if/then but I'm screwing it up- here's an example of the data:

I looked at this link and tried to model my code after this:
Parser:
Python

Expression:
Reclass(!WELL_YIELD!)

Code Block:
def Reclass(WellYield):
    if (WellYield >= 0 and WellYield <= 10):
        return 1
    elif (WellYield > 10 and WellYield <= 20):
        return 2

I have no idea what I'm doing and I tried this
def Reclass(class):
    if (class== 'motorway' OR 'primary' OR 'trunk'):
        return 'Primary'
   el if (class== 'secondary'):
        return 'Secondary'


texts: Deploy GeoServer inside Apache Tomcat in Google Cloud Platform

texts:
I am currently working on a little GIS web app. I have deployed Apache Tomcat on Google Cloud Platform

and I want to install GeoServer on it. I have downloaded the web archive of GeoServer 2.16 but when I want to deploy it on Tomcat, I have an error as you can see below

From this question, it seems to be a file size error, as said "default tomcat only accepts 50MB in the WAR files to upload, and the web archive of GeoServer is 98.1 MB. So I need to change this"

Since my Tomcat is on the cloud, please how to change this constrains?


texts: find all circles intersecting another in a table

texts:
I am using postGreSQL (9.3) and postGIS (2.1).
I have a table containing names and circles (point + radius).

Given a name and a circle (point + radius), I would like to find all the database lines that contains (circles that intersect my circle) AND (share the same name). Is it possible using postGIS ?

Here is what I have so far, taken from other questions on this site.  Please note that I am just beginning to learn postGreSQL :

--create table

CREATE TABLE circles (
    city VARCHAR(50),
    geog geography
);

-- insert sample values

INSERT INTO circles (city, geog)
VALUES ('Lannion', geometry(ST_Buffer(geography(ST_GeomFromText('POINT(48.732084 -3.459144)')), 1000)));

INSERT INTO circles (city, geog)
VALUES ('Lannion', geometry(ST_Buffer(geography(ST_GeomFromText('POINT(48.732084 -3.459144)')), 2000)));

INSERT INTO circles (city, geog)
VALUES ('Lannion', geometry(ST_Buffer(geography(ST_GeomFromText('POINT(48.732084 -3.459144)')), 3000)));

INSERT INTO circles (city, geog)
VALUES ('Lannion', geometry(ST_Buffer(geography(ST_GeomFromText('POINT(48.732084 -3.459144)')), 4000)));

INSERT INTO circles (city, geog)
VALUES ('anotherName', geometry(ST_Buffer(geography(ST_GeomFromText('POINT(48.732084 -3.459144)')), 1000))); -- not selected because the name is different

-- search (need help for this one)
SELECT ST_Buffer(geometry(ST_MakePoint(48.732084, -3.459144)), 2500) FROM circles;

I can't make sense of the output of the SELECT query so far. Also, it doesn't have any constraint on the city name. I would like to have as output the first two inserts in this case.


texts:
I am making a python code for Arcgis9.3. In this code I have an input feature and I create a Output feature. I want that the output feature is the same geometry type as the input feature. At the moment my code to create the output feature is :

gp.CreateFeatureclass("C:/Users/Desktop/Test/", NameOutFC, "POLYLINE")

But I would like to get the input feature geometry type and use it to create my output feature :

gp.CreateFeatureclass("C:/Users/Desktop/Test/", NameOutFC, GeometryLikeInput)

I found the code with ArcPy :

desc = arcpy.Describe('C:/shapefile.shp')
geometryType = desc.shapeType

if geometryType == 'polygon':
print 'Polygon'
elif geometryType == 'point'

But I can't use it with arcgis9.3.


texts: Map on canvas renders after image saving

texts:
I'm a beginner to Qgis and Python;
i need a script (don't need it to be a plugin) to get maps of selected features from a layer satisfying certain conditions;
Here is my code:

from qgis.gui import *
from PyQt4.QtCore import *
from PyQt4.QtGui import *
import time
fnam ='C:/Users/User/.qgis2/Test/Output/'
canvas = qgis.utils.iface.mapCanvas()
#print canvas.size()
Nom =QInputDialog.getText(None,"Value", "Enter text",0)
fnam2 = fnam + Nom[0] + '.txt'
#print fnam2
outfile = open (fnam2, 'w')
layer = iface.activeLayer()
layer.setSelectedFeatures([])
#provider =layer.dataProvider()
selection=[]
for f in layer.getFeatures():
    n1 = f['Comune_1']
    n2 = f['Comune_2']
    n3 = f['Comune_3']
    n4 = f['Comune_4']
    n5 = f['Comune_5']
    if (n1 == Nom[0]) or (n2 == Nom[0]) or (n3 == Nom[0]) or (n4 == Nom[0]) or (n5 == Nom[0]) :
        selection.append(f.id())
        line = '%s\n' % (f['link'])
        unicode_line = line.encode('utf-8')
        outfile.write(unicode_line)
layer.setSelectedFeatures(selection)
#provider.select()
outfile.close()
mnam ='C:/Users/User/.qgis2/Test/Output/' + Nom[0] + '.tif'
mnam2 ='C:/Users/User/.qgis2/Test/Output/' + Nom[0] + '.png'
#print mnam
canvas.zoomToSelected()
canvas.refresh()
canvas.update()
canvas.saveAsImage(mnam,None,'tif')
#time.sleep(3)
canvas.saveAsImage(mnam2,None,"PNG")
#print 'Finished'

My problem is that saved image(s) show the map at the moment i press "run" button in Python console, not the actual map based on the last selection;
i tried to insert a time.sleep() call to delay saving, hoping it could help waiting rendering to complete before saving, but it doesn't work;
i also tried recursing the whole process with a for loop, neither working;
any tips on what am i doing wrong?

I'll try to explain in a better way:
in this project i have two layers/shapefiles:
one shows villages in my region with their boundaries;
the second one shows the arrangment of orthophotos in a 1:5000 scale of the territory (img1)

now professionals or village administrators may ask for orthophotos of a specific village;
running my script i wish:
- have a list of orthophotos of that specific village (this works fine with txt file i create);
- give a map showing the arrangement of requested orthophotos as a complementary aid.
Now, if i run the script selecting, say, "Potenza", i eventually can see in Qgis the updated map with desired selection, but saved images show the whole region:(I mean "Potenza.png" is identical to above image while it should be like the one below);
only now i noticed there are selections, so the problem seems to be in zoomToSelected command.

if i rerun the script and i select another village, say "Maratea", this is what i'll find as Maratea.png

(tif images are identical to png);

(When i'll have a working script selected "rectangles" will show image name, and i won't use inputdialog but i'll get village names from the first shapefile and run it only once for all villages).


texts: Extract the edge between two raster cells with different values

texts:
Is there a way in R to create lines (black lines in the figure below) along the edge between two raster cells that have different values ?

Update

I'm trying the functions clump() and rasterToPolygons() as in the post [https://stackoverflow.com/questions/28859181/how-to-get-contour-lines-around-the-grids-in-r-raster].

First, I reclassified my raster "r" to have cells with two different values:

new.values <- cbind(c(0, 1, 2, 3, 4, 5, 6, 7, 8), c(NA, NA, 1, NA, 2, NA, NA, NA, NA))
new.r <- reclassify(r, rcl=new.values)

Then, I used the functions rasterToPolygons and clump:

test.edge <- rasterToPolygons(clump(new.r), dissolve=TRUE)
plot(test.edge, add=TRUE)

And I obtained:

How can I extract the edges that have not contour lines (for example, edges near red arrows) because they are located between two different values?


texts: Getting the distance traveled across road segments

texts:
so I have a query that takes an encoded polyline, segmentizes that polyline into smaller pieces, and matches each point on that segmentized line to a "road segment" from a separate table in my database. This query works well and efficiently, matching each point to its closest road segment. What I want to be able to do is then find the total distance traveled within each returned matched segment, but I'm struggling with that part.

Here is the query:

WITH
  points AS (
  SELECT (ST_DumpPoints(
         ST_Segmentize(ml, .000032))
       ).*
  FROM ST_LineFromEncodedPolyline('xxx') AS ml
)
SELECT *
FROM (
  SELECT
     seg.xdsegid,
     pts.path,
     ST_Distance(pts.geom::geography, seg.geom::geography) AS distance
FROM points AS pts
JOIN LATERAL (
  SELECT st.xdsegid,
         st.geom
  FROM alabama AS st
  ORDER BY pts.geom <-> st.geom
  LIMIT 1
) AS seg
  ON true
  ORDER BY seg.xdsegid, path
) q
ORDER BY path;

And it currently returns rows that look like this:

I know I probably need to turn each segment into a line and get the distance that way, but I can't for the life of me figure out how to best go about grouping and measuring each one of the returned segments. Also, sometimes the segments will return in such a way that a small portion of the points will momentarily be matched to another segment before returning to the original. Not sure if this is problematic.


texts: Prevent parameter defaults in service definition file

texts:
I have a Python toolbox with a single tool designed for 10.2. The tool is intended to be a geoprocessing service on a server. So I run the tool, and after it succeeds, I right click the result and choose Share As --> Geoprocessing Service. In the wizard, I choose "Save a service definition file" and move through the wizard until it switches to the dialog where I can set service properties. I Analyze, and there are no errors or warnings.

At this point, I notice that my parameters are being given default values. In particular, they get the values that I used when I ran the tool to create the results file. ArcGIS's dialog will not let me modify the default parameters; the text boxes are read only.

I want to prevent this. Ideally, my parameters would have no default values, but empty string default values would also be acceptable. All my input parameters happen to be strings. How can I prevent ArcGIS from using the result's parameters as the defaults?


texts: My grids aren't displaying in degrees in QGIS

texts:
I have a added grid coordianates to my map but it is giving wrong/meaningless coordinates.
How do I make it display the correct coordinates for my map?


texts: Integrate leaflet module in a google maps application

texts:
Is it possible to use a leaflet application over a google maps application?

What I want is to integrate a leaflet module / application over a google maps application just by passing the google maps DOM element to my leaflet module and not just change the leaflet module baselayer.

Is it a way to do that?

I was thinking to create 2 distinct maps and "share the events" like the mouse interactions but I don't know if it could work.


texts: Get UTM regions overlapping with a geojson polygon

texts:
how can I get the list of UTM grids overlapping with a geojson polygon in web mercator? I am trying to download data off of Sentinel, and I am seeing if I can go directly from an input polygon to downloading a tile off of the Sentinel bucket on AWS/GCloud.
My current workflow is to use sentinelsat to get the scenes, but I want to download the scenes directly from the bucket without any intermediary libraries, and knowing the UTM grids will help me get there.


texts: Differences in first-order Polygon Contiguity in ArcGIS - Contiguity Edges Corners (CEC) vs. Contiguity Edges Only (CEO)

texts:
I am working on spatial autocorrelation, and using the spatial modeling relationship methods available in ArcMap.
In ArcMap, there are the following methods for polygon contiguity:
Contiguity Edges Only (CEO), which is using only coincident boundaries in order to create a spatial relationship (rooks) while Contiguity Edges Corners (CEC) uses coincident edges and corners of the polygon (queens). Link follows: https://pro.arcgis.com/en/pro-app/latest/tool-reference/spatial-statistics/modeling-spatial-relationships.htm
Now, I'm trying to figure out how exactly these methods are being used, and the pros/cons of each. I've been trying to look up the algorithms and spatial statistics related to these methods, unfortunately, to no avail. My guess is also that the Contiguity Edges Corners case would be much more computationally intensive, since it needs to explore more cases. Also, I am not quite sure how to interpret any spatial statistics generated from these 2 methods, and how to compare and contrast them in why different clusters may be highlighted (sometimes drastically differently) in ArcMap.
How does ArcMap (ArcGIS Desktop) implement this?
Another factor is weighting (by inverse distance, inverse distance squared, custom weight matrices, etc.) but I would like to save this for another time.


texts: Where can I find the latest WMX (JTX) Object Model Diagrams for ArcGIS 10?

texts:
I'm looking for a version that includes ESRI.ArcGIS.JTXUI--I'm hoping someone knows of an updated and online version of C:\Program Files\WMX\Developer Kit\Diagrams\ObjectModel.pdf.

For reference, I'm on ArcGIS 10 SP2:

ObjectModel.pdf is dated: 10/27/2010
SDK.chm is dated 6/4/2010

So these weren't updated by either service pack.


texts: Managing Leaders for Extent Indicators of Insets using ArcPy with Data Driven Pages?

texts:
I am working in data driven pages that need insets on some maps (not all). I found and have successfully managed to use this code below. However, I need to have leaders. Any ideas on how to manage extent indicator leaders in arcpy?

import arcpy, os
# Create an output directory variable
#
outDir = r"C:\TXT"

# Create a new, empty pdf document in the specified output directory
# This will be your final product
finalpdf_filename = outDir + r"\FinalMapBook5.pdf"
if os.path.exists(finalpdf_filename): # Check to see if file already exists, delete if it does
  os.remove(finalpdf_filename)
finalPdf = arcpy.mapping.PDFDocumentCreate(finalpdf_filename)

# Create a Data Driven Pages object from the mxd you wish to export
#
mxdPath = r"O:\Data\project_specific\cedar_milford\milfordtoCedar\Test_MXD\05858_hb_PHI_DE_cedar_to_milford_cut_sheets_JDS.MXD"
tempMap = arcpy.mapping.MapDocument(mxdPath)
tempDDP = tempMap.dataDrivenPages

# Create objects for the layout elements that will be moving, e.g., inset data frame, scale text
dataFrame = arcpy.mapping.ListDataFrames(tempMap, "Test")[0]
narrow = arcpy.mapping.ListLayoutElements(tempMap, "MAPSURROUND_ELEMENT", "North Arrow")[0]

# Instead of exporting all pages at once, you will need to use a loop to export one at a time
# This allows you to check each index and execute code to add inset maps to the correct pages
#
for pgIndex in range(1, 2 + 1, 1):

  # Create a name for the pdf file you will create for each page
  temp_filename = r"C:\TXT" + \
                            str(pgIndex) + ".pdf"
  if os.path.exists(temp_filename):
    os.remove(temp_filename)

  # The following if statements check the current page index against given values
  # If the current page index matches, it will execute code to set up that page
  # If not, the page remains as is
  # Note: If you created a text file containing this information, this is where
  # you would paste in that code

  # Code for setting up the inset map on the first page #
  if (pgIndex == 1):
    dataFrame.elementPositionX = 0.6644
    dataFrame.elementPositionY = 7.5622
    dataFrame.elementHeight = 3.025
    dataFrame.elementWidth = 5.8834
    insetExtent_1 = arcpy.Extent(614813.318397, 499413.131878, 615191.030877, 499607.335952)
    dataFrame.extent = insetExtent_1

  # Code to export current page and add it to mapbook
  tempDDP.exportToPDF(temp_filename, "RANGE", pgIndex)
  finalPdf.appendPages(temp_filename)

  # Clean up your page layout by moving the data frame and resetting its extent after exporting the page
  # This will reset the page to the basic layout before exporting the next page
  #

  dataFrame.elementPositionX = 1000 # Move inset data frame off the page
  narrow.elementPositionX = 1000
  dataFrame.scale = 350000000 # Change scale so extent indicator no longer visible in main data frame
  arcpy.RefreshActiveView()

# Clean up
#
del tempMap

# Update the properties of the final pdf
#
finalPdf.updateDocProperties(pdf_open_view="USE_THUMBS",
                             pdf_layout="SINGLE_PAGE")

# Save your result
#
finalPdf.saveAndClose()


texts: Buffer road line, taking into consideration elevation and upslope/downslope

texts:
I have a roads polyline layer. I have vector-based contour lines. I have converted those contours to a TIN, then created a raster based on that TIN. Now I have my roads and my elevation raster. Now I need a 200ft buffer on the upslope and a 150ft buffer on the downslope. There are many posts regarding this...but nothing simple. I need step by step directions.

I am using ArcGIS 10.4 Advanced with all the bells and whistles.


texts: Assign_vertex_id function does not exist

texts:
I would like to add unique start/end nodes for each feature in my road network.  I ran across some posts referencing the pgrouting assign_vertex_id function, which appears to do exactly what I want, but when I try to use it, I get an error saying the function does not exist.

PGRouting is installed and enabled on my database.  Using

SELECT pgr_version();

I can see that I am running 2.1.0.  Was this function removed at some point, and if so, is there a good alternative in place?

Update, here is the code I am trying to run, where tucson is a table of roads and the geometry is a Linestring

ALTER TABLE tucson ADD COLUMN source integer;
ALTER TABLE tucson ADD COLUMN target integer;
SELECT assign_vertex_id('tucson', 0, 'geom', 'gid');

And then I get the following error:

ERROR:  function assign_vertex_id(unknown, integer, unknown, unknown) does not exist
LINE 3: SELECT assign_vertex_id('tucson', 0, 'geom', 'gid');
           ^
HINT:  No function matches the given name and argument types. You might need to add explicit type casts.


texts: minidump error when importing dxf, causing QGIS to crash

texts:
QGIS (2.18) crashes when I use dxf importer to import a simple dxf file (Autodesk 2000 dxf format).

May I ask how to solve this problem?

error message:

  minidump written to
  C:\Users\XXX\AppData\Local\Temp\qgis-20170214-140259-18376-13556-0332f5a.dmp


texts: Prevent ArcGIS from maintaining extents during resize

texts:
Background:

I have a custom ArcGIS application written using ArcGIS Engine 10.0 that at its simplest form is just a map control inside of a resizeable form.

The Problem:

During resizing of the form, the performance is absolutely awful and in some cases the entire application crashes (depends on quality of machine).  This is because the MapControl constantly refreshes in order to maintain the same extents of the map.

Question:

Is there a way to prevent the MapControl from refreshing to maintain its extents?


texts:
I used the Bearing Distance to Line command (Data Management) in arc to convert UTM coords (i.e. of a plot center), bearing and distance to the figure shown below.  The end of each line represents a tree location.  How can I convert these lines to points in order to represent tree locations as a point shapefile?  Thanks for the help.

Edit:  The final model and final product based on help from others:


texts: NoData Value = 0, problem when creating NDVI, Erdas

texts:
I try to create an NDVI from Landsat and ASTER satellite images using Erdas Imagine. In Erdas, the satellite images have a nodata value of 0.

After creating an NDVI, the nodata value is still 0, but obviously there are also file pixels with the value 0 now (NDVI ranges from -1 to +1). So if I reset the NoData value to 0 (it is not done automatically), some pixels in my raster are missing since they are now classified as "nodata"

Any suggestions?


texts:
I've installed QGIS 3 today. I found the xyz tiles option and have added e.g. google, openstreet, and ESRI Topo. These all work and show the map. Now I'm looking for a topo map, so ESRI topo is a nice starting point. Unfortunately I find it a bit blurry on the font, and the lines a bit to light.

Does someone know working XYZ topo entries? I can find some XYZ tile servers for topo, but they are often not downloading tiles


texts: DLLs or examples to support writing a shp, shx, dbf file as memory stream for uses on Azure website?

texts:
I'm trying to read/write shapefiles for a web service for an azure website and I do not have the luxury of accessing the local file storage (which means that reading/writing from/to disk is out of the question as many of the DLLs read and write via local disk). I've come across libaries such as DotSpatial, NetTopologySuite, ShapeLib, GDAL's Ogr2Ogr which does very nice things for desktop when it comes to read/write shapefiles, but since I don't have access to local disk, I have to come up with something different.

What I have so working so far is reading the dbf and shp file from a blob storage and parsing the binaries to get the geometry and property data and send that down via json or text.

I'm currently working exporting my data (edit: the data is stored in SQL server table on Azure, this is not a geom table, the lat/lng are in floats and properties columns of the table) to shp, dbf, and shx files in webapi and send that down as a zip for download. What I come up with is building the binary array, which is very tedious. I'm hoping that there's a tool out there that can facilitate this process.


texts: How to calculate Night Flight time for a given flight between 2 locations?

texts:
I am currently working on a project (web app) at university to be able to calculate the amount of time spent during night flying. Currently, I have,  Departure Lat and Long and Departure Time in UTC  Arrival Lat and Long and Arrival Time in UTC Sunset and  Sunrise Times (Based on dates) for both Arrival and Departure airfields in UTC
An example:  Depart UK Edinburgh airport on the 23 April 2019 at 11:10 UTC  Arrive Saudi Arabia Riyadh Airport at 18:10 UTC Total Time 7hr  day flight = ?  night flight = ?
I know many places online do this type of calculation for you, so I am wondering if there is a formula or method available that allows me to do this calculation also taking into account the direction of travel etc (Which I can then convert to code)?


texts: Stopping Kernel Density from cutting off symbology?

texts: I need help configuring Kernel Density so that is does not cut off symbology on hot spots. (See pic below) I'm assuming it has something to do with the extent but I could not get it to display the way I want it.


texts: ArcGIS Server returning all feature attachments as text/plain

texts:
Using the REST endpoint on a feature service, it appears that all attachments are returned with a context type of "text/plain" rather than what is displayed for the attachment.  In my case, I have a pdf which has a content type of application/pdf, but when I click the link to view/download it (from the .../FeatureServer/0/1234/attachments page), it just dumps the encoded file as a text stream.

Does anyone know how to make this work properly?  My only thought it to just write up a custom servlet to return the file with the correct content type set, but this seems like something that should work properly through the REST endpoint, particularly since ESRI goes to the trouble of saving the actual content type in the attachments table.


texts:
I have one polygon file containing all German municipalities as well as a point file. For each municipality, I would like to compute the average of a certain attribute of all points within the municipality's borders. I use the spatial join command for this. However, I am getting the error message that CRS do not match and thus I do not receive any results.

The polygon file seems to be EPSG:25832, the following shows the .prj content

PROJCS["ETRS_1989_UTM_Zone_32N",
    GEOGCS["GCS_ETRS_1989",
        DATUM["D_ETRS_1989",
            SPHEROID["GRS_1980",6378137.0,298.257222101]
            ],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]
    ],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",500000.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",9.0],
    PARAMETER["Scale_Factor",0.9996],
    PARAMETER["Latitude_Of_Origin",0.0],
    UNIT["Meter",1.0]
]

The point file is in the format EPSG:31467 (Gauß-Krüger-projection in the 3. meridial zone, Ellipsoid Bessel, date Potsdam (central point Rauenberg)).

First, I add the polygon vector file. Then I import the text delimited file. Polygons and points overlap, but I cannot conduct a spatial join because of non-matching CRS (obviously).
However, when reprojecting the latter file to the former, polygons and points do not overlap anymore and I do not get any results of the join (the alert that layers have different CRS does NOT pop up anymore). The way how I reproject is by right clicking on the layer - save as ... - choose corresponding CRS.

What am I doing wrong?

In reference to @AndreJ's comment, here are the Extents for both unchanged layers.

Polygon Layer:

In layer spatial reference system units
xMin,yMin 280371.06,5235855.98 : xMax,yMax 921292.37,6101443.71
Layer Spatial Reference System
+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Point Layer

In layer spatial reference system units
xMin,yMin 3281000.00,5238000.00 : xMax,yMax 3920000.00,6102000.00
Layer Spatial Reference System
+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs


texts:
Can someone explain how the ogr2ogr simplify works?

I am trying to use raster data which has a spectral resolution of 30m. The vector file that I am trying to import into Google Fusion is large in size and so I was recommended to use ogr2ogr with tolerance of 0.00001.

Can someone explain the unit of the tolerance?

I assumed that it would be decimal degrees but the resulting KML still has around 10 decimal points.


texts: ArcGIS Online Print Template Legend Override

texts:
I have question that some of us in the forum raised in Stack Exchange in regards to the appearance word override in the legend (map). It seems the the map issued the legend (ArcGIS Online Basic Viewer) for Felda correctly (divided into different groups).

However, symbology for my feature hosted layers such as BN winning margin and 2013 Coalition winners, the word override appears everywhere in the legend. This issue also appeared when I built a WebApp builder configured for printing.

The solutions that was presented on the stack exchange that requires coding is mind boggling. Can anyone explain in simple GIS terms where possible on how this legend truly shows the values that ought to be?


texts: Create new field based on the sum of occurrences of a feature attribute

texts:
I have an Owner_Name field. I want to add a virtual field containing the number of properties each owner owns. (Num_Properties).
To figure this I believe I need the expression for the Num_Properties field to count the number of occurrences of the same owner name within all the features/rows.
I'm not sure what my expression needs to be.  I've tried so many!  I'm assuming that I should be using count().
count( "Owner_Name" , filter:= "Owner_Name" = attribute($currentfeature,'Owner_Name')

That doesn't work, just an example of what I've tried.  Apparently that elicits the same value on both sides of the "=" for the filter.


texts: Hierarchy not working in network analyst

texts:
I have created a simple network to test hierarchy. But apparently I am missing something.
I have created an attribute for my roads called hierarchy and given it values of 1 to 3. The network is set to use Hierarchy but it is apparently going through the shortest route not the hierarchy.

I am including a photo.

The red and yellow lines : hierarchy 1
The green lines: hierarchy level 2
the orange line: Hierarchy level 3

So in my photo I would imagine to go from point 1 to 2 , it should go through the red or yellow path but it goes through the orange with hierarchy 3. Am I missing something?

!   enter image description here


texts:
I am trying to calculate road network properties while using custom distance metric in sDNA – in my case time which differs in segments depending on availability of public transport.

Initially I have imported a table from Excel (where for different groups of polylines the different formulas of “time” were used) and have tried various number formats in attributes table: short, long, single, double, numeric, etc.

For integral analysis I choose custom distance metric and the attribute column as a “custom metric data field”, but the calculations are always unsuccessful with the following message: “Calculation expected numeric field “…”” not found in data”.

Calculations for the same network with angular distances and Euclidean radiuses work well. Symbology displays custom data field.

Any idea what could wrong in the data?


texts: Creating ArcSDE feature class in MS SQL Server?

texts:
Unfortunately I am totally dumb when it comes to programming, and very new to databases (mainly because I studied geography), so :

I have a Microsoft SQL Database Server, and a ArcSDE geodatabase. I made a view in SQL Server Manager, which I can see it in my geodatabase (in ArcCatalog), as a table. Which is great. But I also need the geographic information not just the data from the tables.

So how can I create a feature class in my geodatabase from a view that I created in SQL Server Manager. I did select the "shape" field, but sadly it's not that simple:P

I use: Microsoft SQL Server 2008, ArcGIS 10.0, ArcGIS Server 10.


texts:
What is the best practice to manage temporary blocked roads (works on the road) with pgRouting?

Is there a dedicated feature adding exception?

If not, i give an heavy weight to the edge?


texts: Where is the legend for graduated styling?

texts:
I created a map using graduated styling on a layer and there is a field for legend format and a column that says legend, but I can't find the legend anywhere.
Does anyone know how to access it?
(The only thing in the legends tab within the layer properties window is something that says transparency slider)?
Ideally I would like to save it as a png.


texts: Specifying Bayesian Spatially lagged dependent (SAR) variable regression in winBUGS or Stan

texts:
I have been looking for code example on how to specify a SAR spatially lagged dependent variable regression in winBUGS (geoBUGS) or, if possible, in Stan.

I mean, something of the format:

Y = beta_0 + beta_1*Y*W + ... + e

Where W is the weight matrix identifying neighborhood of geographical areas.

Notice that I am not looking for a CAR model, i.e. when the errors are spatially auto correlated, but specifically for a SAR lag model.

Can anyone help with some code example or any source suggestion?


texts: Data interchange between Civil3D 2016 and QGIS 2.18

texts:
Does anyone have examples of good practice for dynamic interchanging of data between QGIS and Civil3D?  Specifically, importing vector shapes from QGIS into Civil3D (and occasionally vice versa).

I am working on an engineering project which has a drafting team working in Civil3D, with our wider technical team using QGIS as a 'spatial sketchpad'. Often the design team will sketch out something e.g. an overall building footprint in QGIS and then pass this through to the drafters for further development.  We have been using the shapefile import facility in Civil3D but this has to be reimported every time a change is made.  Is there some format or database (GDB? PostGIS?) that could sit between the two software and dynamically update spatial data?


texts: How to populate the geometry_columns table

texts:
I have a SQL Server table of vector data (containing  latitude, longitude, UComp and VComp columns), that I am trying to add as a layer to QGIS.

Here is a sample of the data in the SQL Server data table:

And the table design:

I have successfully connected to my database from QGIS and also created the geometry columns table.

What data should I be adding to the [f_geometry_column], [coord_dimension] and [geometry_type] columns in the geometry_columns table?


texts:
One of my biggest frustrations with OSM is adding non-forest/wood areas to an otherwise completely tree-covered area of land.

I'd like to know if there's any way to invert the selection of say, grassland, to have everything around it be wood. Or conversely, if it's possible to place grassland on a wooded area that would cut it out, instead of just overlapping it?


texts: QgsVectorLayer setRenderer problem with lingering reference to old renderer

texts:
I have a stand-alone pyqgis application that can display vector data. I have implemented a routine that sets layer's renderer as below:

def set_sym(layer):
    symbol = QgsMarkerSymbolV2.createSimple({})
    symbol_layer = QgsSimpleMarkerSymbolLayerV2()
    symbol_layer.setColor(QColor(0, 0, 0, 255))
    symbol_layer.setName('circle')
    symbol_layer.setSize(1.0)
    symbol.appendSymbolLayer(symbol_layer)
    renderer = QgsSingleSymbolRendererV2(symbol)
    layer.setRendererV2(renderer)

It works. However, if I call it a second time on a layer, then immediately go to exit the program, it hangs & I have to kill the background python process to get out of it. I've already checked that the newly created renderer and the layer.rendererV2() are different and their symbol() are also different (I mean the memory addresses are different). So, it makes me think that the setRendererV2() routine is not replacing the renderer cleanly, hence hanging the program upon quitting right after the second call on the same layer. Please let me know if there is something I missed that caused this.


texts: QGIS RGB raster layer: min/max-values defined by user not working?

texts:
I tried several times to define the min/max-values for the RGB-display under properties>style: Benutzer min/max-Werte (user min/max-values) but it is not working. By chance it worked once but now it is not working any more. The other functions as to invert the color representation or the selection of the channel are working.

The raster file are created by the following process:

Scanned as an PDF
imported in Photoshop CS5 (modus RGB)
saved as (in Photoshop) tif (color profile not embedded, respectively sRGB IEC61966-2.1 not activated).

Does anybody have an idea where the problem is?


texts: When i import NetCDF file (.nc) to QGIS then I get a picture that consists of quite big pixels. Is it somehow possible to soften those pixels to get something like this:


texts: How to export layer with metadata in QGIS?

texts:
I have set some metadata for my layer that I would like to export. I export the layer using the following:
_writer = QgsVectorFileWriter.writeAsVectorFormat(layer,"path/layer.shp",'utf-8',driverName='ESRI Shapefile')

However the saved layer does not contain the metadata. How can I export the layer and keep the metadata?


texts: How can I fill polygon with strokes in openlayers 3?

texts:
I want to create a polygon and have it filled with strokes as in the image :

Currently my style function looks like this :

var tempObjectsStyleFunction = function(feature, resolution) {
return [
    new ol.style.Style({
        fill: new ol.style.Fill({
            color: feature.get('color')
        }),
        stroke: new ol.style.Stroke({
            color: feature.get('color'),
            width: 5
        }),
        image: new ol.style.Circle({
            radius: 20,
            fill: new ol.style.Fill({
                color: feature.get('color')
            })
        })
    })
];

};

What do I need to modify in order to support strokes in the fill ?


texts: Leaflet: Change displayed column on click

texts:
I have a GeoJSON layer in a Leaflet map showing buildings. The columns of the buildings layer contain electricity usage info by date (eg., the column 12052015 contains electricity usage on Dec 5, 2015).  I'm trying to make a button that will style the layer based on a different column.

Edit: here's where I define the GeoJSON layer from an AJAX request:

   var mylayer = L.geoJson(data, {style:style, onEachFeature: onEachFeature}).addTo(map, mapStyle);

So far, this works:

    $('#clicker').click(function () {
            mylayer.setStyle(teststyle());
    });

    function teststyle(feature) {

        return {
            fillColor: getColor(feature.properties["12052015"]),
            weight: 2,
            opacity: 1,
            color: 'white',
            dashArray: '3',
            fillOpacity: 0.7
        };
    }

But what I really want to do is pass the name of the new column as an argument, something similar to this example:

        $('#clicker').click(function () {
            mylayer.setStyle(teststyle(mylayer, "01052016"));
    });

    function teststyle(feature, colname) {

        return {
            fillColor: getColor(feature.properties[colname]),
            weight: 2,
            opacity: 1,
            color: 'white',
            dashArray: '3',
            fillOpacity: 0.7
        };
    }

When I do this,  however, I get this console error: "Uncaught TypeError: Cannot read property '12/10/2015' of undefined".

Anyone know why I'm having trouble passing in the feature argument?


texts: QGIS2 to QGIS3 plugin migration issue AttributeError: type object 'QObject' has no attribute 'connect'

texts:
I am trying to migrate a plugin. I have these two lines:

QtCore.QObject.connect(self.buttonBox, QtCore.SIGNAL(_fromUtf8("accepted()")), Configuration.accept)
QtCore.QObject.connect(self.buttonBox, QtCore.SIGNAL(_fromUtf8("rejected()")), Configuration.reject)

I thought after looking into some examples that in PyQt5 the lines should read:

self.QObject.connect(self.buttonBox, QtCore.SIGNAL(_fromUtf8("accepted()")), Configuration.accept)
self.QObject.connect(self.buttonBox, QtCore.SIGNAL(_fromUtf8("rejected()")), Configuration.reject)

but it just throws the error

AttributeError: 'ConfigurationDialog' object has no attribute 'QObject'

This is confusing because ConfigurationDialog is not found in this py script. ConfigurationDialog is a class in another script but nowhere in it does it make reference to QObject. Can I get some clarification on what where I am going wrong?


texts: Finding a Conceptual Site Modeling Software

texts:
I've been asked by my company to generate a conceptual site model for a project and I have no clue where to start. The purpose is to show groundwater flow through a few layers of rock either in 2D or 3D. I don't need a program like MODFLOW, just something to visualize it.

My search always ends at a website that describes how to make a CSM or I read articles with great pictures/figures with no explanation on who/what created it. What I want is to download (doesn't have to be free) software that will help my model look professional.

Does anyone have a suggestion for me?


texts: Setting SAGA folder for QGIS with OS X for macports install?

texts:
I am trying to get the toolbox to work in QGIS 2.6, but I keep getting an error message that SAGA is not properly configured. After much digging, I finally found I can set the SAGA folder in Processing > Options > Providers > SAGA. I installed SAGA with macports. The saga_cmd executable is in my /opt/local/bin folder, but setting this as my folder does not fix the problem. I continue to get the error message.

Does anyone know how I can configure QGIS to find my saga installation?


texts: Is there alternative to tools like MIPT and Discover in QGIS?

texts:
I am getting fond of QGIS but I am unable to find tools for frequency plannig and other stuff. Are there any alternative to MIPT and Discover of Mapinfo in QGIS?

MIPT is used for frequency planning in telecom. It creates cell sector and is useful for neighbour correlation. and discover has link here. and more information


texts: Measuring the ground coverage in % of a transparent getiff along a path of a vectorline in QGIS

texts:
I am searching for a way to measure the ground coverage in QGIS. I have an high resoluted orhtomosaic of an agricultural plot trial where i have extracted the organic material on the surface:

The result of the extraction is a transparent orthomosaic with the organic material on the surface:

Now, i am looking for a way to measure the percentage of organic material of a determined area. The area is demarced by a buffered vectorline.

In my case thera are 14 vectors.

Can someone show me how to do this?


texts: Improve Smartphone GPS with external RTK - will a simple antenna work?

texts:
Regarding RTK technique, if you have appropriate equipment for full RTK you will have up to cm accuracy. I wondered if we can use very simple external to device antenna (I'm not speaking about Garmin GPS device style build in antennas like we all had in past mobile phones) communicating with the smartphone via wifi/bluetooth/rfid + software layer running on the smartphone, to somehow take advantage on RTK technique in order to improve current GPS accuracy which is ~3 meters to let's say even something like ~1 meter?

Will we compromise accuracy by using very primitive external hardware - So is a simple antenna relevant for RTK?

Speaking about this article

So they said "one base station and one on a vehicle, connected via radio link" but not explaining + $900 is too much. The question is very simple: Can we have the equivalent to "base station" as the smartphone, and equivalent to "one on a vehicle" as simple piece of metal-e.g. antena, all together costs nothing comparing to smartphone price, in order to get 1 meter accuracy compromising the 1 cm accuracy because still we modified here cheaper constellation?


texts: Handling "bad layers" in QGIS version 1.9

texts:
I'm having a problem every time I try and open any project in qgis 1.9. It will load part of a few layers and then a yellow bar comes up across the top of the map area that says "Handle Bad layers: 23 of 23 bad layers were not not fixable."

(I was going to post a screenshot but apparently I don't have enough "reputation" to do so.)

It doesn't tell you any thing about what's wrong with any of these layers or allow you to do anything other than close the warning!

Does any one have any ideas? Please help!!!

***Sorry about that, this is from a project that I made. I'm using pgAdmin III and the layers are postGIS and it was working fine yesterday when I shut my machine down. I have the log message if that would help?

When you bring up the log messages, there are four tabs: General, Plugins, Python and PostGIS. This log was under the Python tab.

warning:/usr/lib/python2.7/dist-packages/osgeo/gdal.py:81: DeprecationWarning: gdal.py was placed in a namespace, it is now available as osgeo.gdal
  DeprecationWarning)

traceback:  File "", line 1, in
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 212, in startPlugin
    plugins[packageName].initGui()
  File "/usr/share/qgis/python/plugins/plugin_installer/installer_plugin.py", line 69, in initGui
    plugins.getAllInstalled()
  File "/usr/share/qgis/python/plugins/plugin_installer/installer_data.py", line 634, in getAllInstalled
    plugin = self.getInstalledPlugin(key, False, testLoad)
  File "/usr/share/qgis/python/plugins/plugin_installer/installer_data.py", line 531, in getInstalledPlugin
    exec("import %s" % key)
  File "", line 1, in
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 452, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/home/localadmin/.qgis2/python/plugins/image_cutter/__init__.py", line 26, in
    from image_cutter import Image_Cutter
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 452, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/home/localadmin/.qgis2/python/plugins/image_cutter/image_cutter.py", line 33, in
    from dlgImageCutter import ImageCutterDlg
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 452, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/home/localadmin/.qgis2/python/plugins/image_cutter/dlgImageCutter.py", line 30, in
    import gdal
  File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 452, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/usr/lib/python2.7/dist-packages/gdal.py", line 3, in
    deprecation_warn('gdal')

Under the PostGIS tab it just says: invalid PostgreSQL layer

I tried manually re-adding the layers that didn't load earlier today and it gave this under the PostGIS tab:

Erroneous query: SELECT styleQML FROM layer_styles WHERE f_table_catalog='howell' AND f_table_schema='original' AND f_table_name='oh_secondary' AND f_geometry_column='the_geom' ORDER BY CASE WHEN useAsDefault THEN 1 ELSE 2 END,update_time DESC LIMIT 1 returned 7 [ERROR: relation "layer_styles" does not exist LINE 1: SELECT styleQML FROM layer_styles WHERE f_table_catalog='how...

I have no idea what styleQML is...or layer_styles for that matter. lol.


texts: How to choose pixel size while doing interpolation in QGIS?

texts: How do you choose pixel size X and Y while doing TIN interpolation? What determines the size?


texts: Splitting Polygons with Lines in ArcGIS Pro?

texts:
I'm trying to use the split polygons with lines tool found here in ArcGIS Pro to subdivide a city's corporate boundary (polygon) with a street centerlines layer (lines) in order to subdivide the city into its blocks.

However, when I run it, the results are just the city's corporate boundary, the line geometry hasn't had any effect:

Any idea what I can do to get the results I'm looking for?


texts: Openlayers tiles source with poor resolution, how zoom out?

texts:
I am using Openlayers to display a map from Vienna Open Data (Capabilities) in an app for mobiles, which have higher DPI screens.
The resolution of the tiles is not so good (I suppose because they are JPEG instead of PNG).

Is it possible to use the tiles from the upper zoom level as if they where smaller tiles for the current level, faking the resolution?


texts: Changing shapefile longitude from -180-180 to 0-360 in ArcGIS?

texts:
I have a world shapefile. The shapefile which I have has Greenwich as the prime meridian. I got the world shapefile from this link.

I have to convert this shapefile into a similar world shapefile as shown below. The longitudes of output shapefile should range from 0 to 359 where the latitude remains the same.

Is there some way to do this thing in ArcGIS?

My output shapefile should look something like as shown above.


texts: How to view layers by data source in QGIS 3?

texts:
Is there a tool or panel in QGIS 3 that allows a user to view the map layers by data source? Similar to the List By Source view in the ArcMap Table of Contents?

I thought I had come across it, but now I can't find it. Perhaps I dreamt it up or perhaps I'm searching with the wrong vocabulary.

Does this capability exist in the core QGIS software? Or a plugin?


texts: Problem running i.gensigset in Grass. Any ideas?

texts: It just returns ERROR: Unable to create signature file.  It doesn't say why... Anyone had this before and found a solution?


texts: Difference between JOIN ON and WHERE in Virtual Layers in QGIS

texts:
There is already a question on this topic Difference between SpatialJoin with “…where…” and “join…on…”. However, it is dedicated to PostGIS.

Since in the Virtual Layer query window, I cannot simply proceed with EXPLAIN ANALYZE command before my SELECT, therefore cannot see the difference between two queries and moreover it uses SQLite and Spatialite to operate.

So, what is the difference between those queries in terms of Virtual Layers and QGIS?

SELECT d.id, COUNT(p.id) PoInPol
FROM districts AS d, points AS p
WHERE st_within(p.geometry, d.geometry)
GROUP BY d.id

Versus

SELECT d.id, COUNT(p.id) PoInPol
FROM districts AS d
JOIN points AS p ON st_within(p.geometry, d.geometry)
GROUP BY d.id

Those queries simply give a number of points inside each feature of districts layer.


texts: Is there any way to set polygon features fill color and outline programmatically?

texts:
By some process I am creating a polygon feature layer. Before adding the layer in the map, I want to set polygon's fill color to no color and set polygon's outline color and width. I am working with C#.

Is there a way to set polygon feature layer's fill color or outline in ArcObjects?


texts: Uploading files to local Geoprocessing service

texts:
I am trying to upload csv file to geoprocessing service, everything works fine if I replaced the URL and uploaded the file to any GP service published by ESRI but when I try to upload to my local GP service it fails. I am using ArcGIS Server 10.4.1 but I'm not sure what did I miss to configure in my service.

var UploadURL="https://MyServer/arcgis/rest/services/Density/GPServer/uploads/upload";
//var UploadURL="http://sampleserver6.arcgisonline.com/arcgis/rest/services/911CallsHotspot/GPServer/uploads/upload";
esriConfig.defaults.io.corsEnabledServers.push(UploadURL);


texts: Intersection of rasters in ArcGIS Desktop?

texts:
I have two set of raster layers, one for buildings (static) and one for depth of surface water due to pluvial flooding from simulation (dynamic, changes with each time step).

My aim is to generate an animation which shows buildings changing their color into red,yellow etc. as soon as water depth reach a certain value. For this purpose, I have to intersect building raster and water depth raster from different time instants.

I have tried it with CON() tool in spatial analyst but the problem with CON() tool is that only the cells where water depth is reached or exceeded in the building raster are changed and not the whole building cells.

Is there any way with which I can manage to change the color of the whole building (and not just the cells where water depth is exceeded) when even a single cell of that building reaches or exceeds a certain water depth?


texts:
I am working on a large modelling and mapping project that will require the processing of significant numbers of 1d models through ISIS Mapper.

Has anyone had any experience of scripting tools in ISIS mapper to batch process the extraction of model results or the interpolation of 1D results?

A google search and hunt through the information provided by Halcrow doesn't provide any suggestions.


texts: SAGA landform classification crashing on MAC, Windows, standalone or with QGIS

texts:
I am trying to apply the "TPI-based landform classification" algorithm from SAGA.
Originally in QGIS 2.14, on a MAC, then on windows and with SAGA standalone, but it quickly saturates memory/CPU and ends up crashing after a few minutes.

I have tried using a smaller elevation map, even with 10 x 10 pixels, as well as changing the projection of the data, but the issue remains the same.

Here is the log from QGIS:

Algorithm Tpi based landform classification starting...
io_gdal 0 -TRANSFORM 1 -INTERPOL 0 -GRIDS "/var/folders/c3/b7snydx51hl3ytjckcs6bq580000gn/T/processing698bae80d9d24acab1a1ba0e08e8edcd/f96bf62b13bb4c3183ece82e146a4dba/DTMAlvaresASTERSAGAproj.sgrd" -FILES "/Users/yannick/Documents/Alvares/Base dados/ASTER DTM/DTM_Alvares_ASTER_SAGAproj.tif"
ta_morphometry "TPI Based Landform Classification" -DEM "/var/folders/c3/b7snydx51hl3ytjckcs6bq580000gn/T/processing698bae80d9d24acab1a1ba0e08e8edcd/f96bf62b13bb4c3183ece82e146a4dba/DTMAlvaresASTERSAGAproj.sgrd" -RADIUS_A_MIN 0 -RADIUS_A_MAX 100 -RADIUS_B_MIN 0 -RADIUS_B_MAX 1000 -DW_WEIGHTING 0 -DW_IDW_POWER 1 -DW_IDW_OFFSET true -DW_BANDWIDTH 75 -LANDFORMS "/var/folders/c3/b7snydx51hl3ytjckcs6bq580000gn/T/processing698bae80d9d24acab1a1ba0e08e8edcd/1ed5f09736e4466a9521db2c1bee40fc/LANDFORMS.sdat"

## #####

SAGA Version: 2.2.3

library path: /Applications/QGIS.app/Contents/MacOS/lib/saga/
library name: libio_gdal
library : GDAL/OGR
tool : Import Raster
author : O.Conrad (c) 2007 (A.Ringeler)

Parameters
Grids: No objects
Files: "/Users/yannick/Documents/Alvares/Base dados/ASTER DTM/DTM_Alvares_ASTER_SAGAproj.tif"
Transformation: yes
Interpolation: Nearest Neighbor

loading: /Users/yannick/Documents/Alvares/Base dados/ASTER DTM/DTM_Alvares_ASTER_SAGAproj.tif

Driver: GTiff
Bands: 1
Rows: 470
Columns: 484

Transformation:
x' = -8.161063 + x * 0.000278 + y * 0.000000
y' = 40.094790 + x * 0.000000 + y * -0.000278

loading band [1/1]

/Users/yannick/.qgis2//processing/saga_batch_job.sh: line 4: 37503 Killed: 9 saga_cmd ta_morphometry "TPI Based Landform Classification" -DEM "/var/folders/c3/b7snydx51hl3ytjckcs6bq580000gn/T/processing698bae80d9d24acab1a1ba0e08e8edcd/f96bf62b13bb4c3183ece82e146a4dba/DTMAlvaresASTERSAGAproj.sgrd" -RADIUS_A_MIN 0 -RADIUS_A_MAX 100 -RADIUS_B_MIN 0 -RADIUS_B_MAX 1000 -DW_WEIGHTING 0 -DW_IDW_POWER 1 -DW_IDW_OFFSET true -DW_BANDWIDTH 75 -LANDFORMS "/var/folders/c3/b7snydx51hl3ytjckcs6bq580000gn/T/processing698bae80d9d24acab1a1ba0e08e8edcd/1ed5f09736e4466a9521db2c1bee40fc/LANDFORMS.sdat"
Converting outputs
Loading resulting layers

The following layers were not correctly generated.
Landforms
You can check the log messages to find more information about the execution of the algorithm


texts: How can I create Mapbox vector tiles with zoom extent 0-4?

texts:
when I use Mapbox and load a vector polygon.
after upload is done the file is converted to vector tiles.
it is then assigned a zoom extent:

How can I tweak this ? I would like the tiles created from this polygon to be available between zoom 0 and 4 and not 4 to 10.

It doesn't work when I try to override the zoom extent visibility from the style.


texts: How to visualize car park occupancy over time?

texts:
I am using QGIS Brighton 2.6.1.

I have a csv table with the exact entrance datetime and exact exit datetime of every car in an specific parking. Example of table:

Nº   DateTimeEntrance         DateTimeExit
--   ----------------         ------------
1    2014.01.01 20:15:16      2014.01.01 21:08:16
2    2014.01.01 20:15:17      2014.01.01 22:11:27
3    2014.01.01 20:15:17      2014.01.01 20:28:09
...      ....                     .....

I have the polygon of the parking and I know that the max capacity of the parking is of 300 cars.

I've been asked to make a heatmap or something similar to paint the polygon redder or lighter taking into account the occupation of the parking (number of cars inside the parking).

I first thought of using the Time Manager plugin but i can't display rasters (heatmaps) in that plugin, although i am starting to think that Time Manager is not the tool i need.

If I try to solve this with a vectorial layer, i don't exactly know how to store in the layer the current capacity of the parking in every single minute for later assign that layer to the Time Manager.

I need some light in this case.


texts: Error inserting new column in ArcSDE table (ERROR 999999)

texts:
I am using a python script to populate a table with values I get from rasters files. I use arcpy.AddField_management to add a new column and arcpy.CalculateField_management to insert the values for the new column.
The table has more than 450 columns and after executing arcpy.CalculateField_management I get the following error:

ExecuteError: ERROR 999999: Error executing function.
Underlying DBMS error [[Microsoft][SQL Server Native Client 11.0][SQL Server]Cannot create a row of size 8074 which is greater than the allowable maximum row size of 8060.]
Failed to execute (CalculateField).

I can avoid this error executing in SQL Server:

ALTER TABLE REBUILD

Is there anyway to execute this SQL server instruction from the script or another equivalent way to avoid the error?


texts: Setting up Hot Key to zoom to Map Scale?

texts:
I'm wondering if I can set up a hot key to zoom to 1:4500 whenever I use it.

Also if I can set up a hot key to zoom in to 1:1250.

I have these two extents set up in the custom extent area and when I use my mouse scroll they come up but I want to be able to hit Ctrl + Q or something along those lines and then my map will zoom to 1:4500.

Is this possible or can you only set up the hot key to turn on zoom where you zoom in by drawing a box on the area?


texts: Grass-gis grid and Shapefile grid: what's the difference?

texts:
I have a raster map with trees in grass (1 and NULL value) for 6m resolution. I would like to aggregate this data in a grid with a resolution of 250m.

For that I have created a grid with v.mkgrid and using v.rast.stats I what count pixel for each grid cell.

It works perfectly in 50 min for a grid of 81323 cells.

To compare with MODIS data I have used R and create a grid with the exact extend and resolution of my MODIS value

grid <- raster(extent(one.ym))
res(grid) <- res(one.ym)[1]
proj4string(grid)<-proj4string(one.ym)
gridpolygon <- rasterToPolygons(grid)
writeOGR(obj = gridpolygon,dsn="/home/user/",layer="mygrid_modis.shp", driver = "ESRI Shapefile", overwrite_layer = T

Once my grid imported in GRASS-GIS I run v.rast.stats but I spend 6h and it is still running.

What's the difference between shapefile grid and Grid from v.mkgrid? How can I improve my shapefile grid time?

By the way, I have created my grid in grass72 using r.to.vect on my landsat raster. It doesn't make a difference in CPU time. In both case I have defined the region size and resolution on my raster.


texts: Where can I find the Sentinel-2 HTML folder in Google Cloud Storage?

texts:
I'm using sen2cor 2.8 to make the atmospheric correction to Sentinel-2 L1C imagery downloaded from google cloud storage. I'm getting an error from sen2cor because it uses a file contained in the HTML directory which seems to be missing in the GCP bucket. I used this Downloading Sentinel-2 data from Google Cloud Storage with Python? to build the .SAFE folder.

How can I download the HTML folder? or How can I use sen2cor without this folder?


texts: Leaflet with EPSG:3765; GridLayer.js:661 Uncaught Error: Attempted to load an infinite number of tiles

texts:
I am not sure about Leaflet 1.7 possibility to show layers in EPSG:3765. Not sure because from Leaflet docs it is not possible but I found something like Proj4Leaflet and it looks like it can... I get some error:
GridLayer.js:661 Uncaught Error: Attempted to load an infinite number of tiles

My code:
//define crs
var crs= new L.Proj.CRS('EPSG:3765',
        '+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs',
        {
        resolutions: [128, 64, 32, 16,
        8, 4, 2, 1, 0.5,
        0.25, 0.125, 0.0625],
        origin: [477174.25,4882262.63],
        bounds:[208311.05,4614890.75,
        724721.78,5159767.36]
        }
        );

can someone explain what is represented with resolutions (is that zoom level), origin and bounds?
//map
var mymap = L.map('map',{
            zoomControl:false,
            crs:crs
            }).setView([477174.25,4882262.63], 15);

//WMS
var lyrDOF = new L.TileLayer.WMS('https://geoportal.dgu.hr/services/dof/wms', {
                layers: 'DOF5_2011',
                format: 'image/png',
                transparent: true,
                version: '1.3.0',
                crs: crs
            }).addTo(mymap);

Error is:
GridLayer.js:661 Uncaught Error: Attempted to load an infinite number of tiles.


texts: Is it possible to combine 2 or more images defined on single variables?

texts:
I first reclassified an existing land use classification and wrote the corresponding new classes on a new variable. Therefore I would like to know whether it is possible in GEE to combine these variables and their images / data again.

var aoi = point.buffer(200000);
var igbpLandCoverVis = {
  min: 1.0,
  max: 17.0,
  palette: [
    '05450a', '086a10', '54a708', '78d203', '009900', 'c6b044', 'dcd159',
    'dade48', 'fbff13', 'b6ff05', '27ff87', 'c24f44', 'a5a5a5', 'ff6d4c',
    '69fff8', 'f9ffa4', '1c0dff'
  ],
};

var LC2001 = ee.Image('MODIS/006/MCD12Q1/2001_01_01')
                .select(['LC_Type1'])
                .clip(aoi);

print(LC2001)

Map.addLayer(LC2001,igbpLandCoverVis, "LandCover2001")

//Reklassifizierung Grassland
var Grassland2001 = LC2001
    .eq(10)
    .or(LC2001.eq(6))
    .or(LC2001.eq(7))
    .or(LC2001.eq(8))
    .or(LC2001.eq(9));
print(Grassland2001)

Map.addLayer(Grassland2001, [],"Grassland2001", false);

//Reklassifizierung Forest
var Forest2001 = LC2001
    .eq(1)
    .or(LC2001.eq(2))
    .or(LC2001.eq(3))
    .or(LC2001.eq(4))
    .or(LC2001.eq(5));
print(Forest2001)

Map.addLayer(Forest2001, {},"Forest2001", false);

//Reklassifizierung Wetlands
var Wetlands2001 = LC2001
    .eq(11)
    .or(LC2001.eq(11));

print(Wetlands2001)

Map.addLayer(Wetlands2001, {},"Wetlands2001", false);

//Reklassifizierung Crop
var crop2001 = LC2001
    .eq(12)
    .or(LC2001.eq(14));

print(crop2001)

Map.addLayer(crop2001, {},"Crop2001", false);

//Reklassifizierung Settlements
var Settlements2001 = LC2001
    .eq(13);

print(Settlements2001)

Map.addLayer(Settlements2001, {},"Settlements2001", false);

//Reklassifizierung Desert
var Desert2001 = LC2001
    .eq(15)
    .or(LC2001.eq(16));
print(Desert2001)

Map.addLayer(Desert2001, {},"Desert2001", false);

//Reklassifizierung Water
var Water2001 = LC2001
    .eq(17);

print(Water2001)

Map.addLayer(Water2001, {},"Water2001", false);`


texts: How to draw a map of postal codes of a country?

texts:
I have data:

postal codevalue
00-002          50
00-003          70
...
99-440          510

where postal code is a postal code in Poland, and value is a number (it can be number of people, density, temperature etc.).
I want to create a map for this data. So I need a map with areas coloured depending of a value. Something like this here. So such data exists somewhere.

How can I get data that I can use to draw the map of Poland divided into areas with their postal code? Where can I find it?
What software can I use to do this?
If this is not a proper place for such question, where can I ask?


texts: Save raster (.tif) with R writeRaster function

texts:
I'm trying to modify a raster through an R script but I find an error when I try to save in the ".tif" format, while in the ".asc" format I do not find any problem.
The script I use is as follows:

library(raster)

path <- "/path/to/input/file.tif"
path2 <- "/path/to/output/file.tif"
open_file <- raster(path)
rasterzero <- open_file * 0
writeRaster(rasterzero, filename=path2, overwrite=TRUE)

The error I get is the following:

 *** caught segfault ***
address 0x565300000009, cause 'memory not mapped'

Traceback:
 1: .Call("RGDAL_SetProject", object, proj4string, PACKAGE = "rgdal")
 2: .gd_SetProject(object, ...)
 3: rgdal::GDALcall(transient, "SetProject", as.character(projection(r)))
 4: .getGDALtransient(x, filename = filename, options = options,     ...)
 5: .startGDALwriting(y, filename, options, setStatistics = setStatistics,     ...)
 6: .writeGDALall(x, filename = filename, format = filetype, ...)
 7: .local(x, filename, ...)
 8: writeRaster(rasterzero, filename = path2, overwrite = TRUE)
 9: writeRaster(rasterzero, filename = path2, overwrite = TRUE)

Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace

Any suggestions on the possible cause?
The same script on different computers runs without errors and the installation procedures of R and its packages are always the same.


texts: Azimuthal Equidistant Projection - Rotate projection

texts:
I am using Azimuthal Equidistant Projection to plot points in a map being drawn in an Android application.

The origin is the centre of the map.

I have the formulas from here and here. I am getting expected results.

Simplified version of the code being used by me in JavaScript:

// Add functions to Math
if (!Math.toRadians)
    Math.toRadians = function(angledeg) {
        return angledeg * Math.PI / 180;
    };
if (!Math.toDegrees)
    Math.toDegrees = function(anglerad) {
        return anglerad * 180 / Math.PI;
    };

/* Mean radius of Earth (in metres)*/
const WGS_1984_R = 6370997;

/* Function to convert geodetic coordinates to Azimuthal Equidistant Projection projection (cartesian metres)*/
function convertGeodeticToAEP(longitude, latitude, longitudeOrigin, latitudeOrigin) {
    var deltaLongitudes = Math.toRadians(longitude - longitudeOrigin),
            cosDeltaLongitudes = Math.cos(deltaLongitudes),
            cosLatitude = Math.cos(Math.toRadians(latitude)),
            sinLatitude = Math.sin(Math.toRadians(latitude)),
            cosLatitudeOrigin = Math.cos(Math.toRadians(latitudeOrigin)),
            sinLatitudeOrigin = Math.sin(Math.toRadians(latitudeOrigin)),
            // Angular distance from the centre
            c = Math.acos(sinLatitudeOrigin * sinLatitude + cosLatitudeOrigin * cosLatitude * cosDeltaLongitudes),
            k = c == 0 ? 0 : c / Math.sin(c),
            x = k * cosLatitude * Math.sin(deltaLongitudes) * WGS_1984_R,
            y = k * (cosLatitudeOrigin * sinLatitude - sinLatitudeOrigin * cosLatitude * cosDeltaLongitudes) * WGS_1984_R;
    return {x: x, y: y};
}

What would be the formulas to rotate the projection by a given angle, from the origin? Say modify the function to include an additional parameter for angle of rotation and then modify the output accordingly...


texts:
I have a set of points in ArcMap which I create a buffer for each point. Then I use intersect to extract the overlapping buffers. The result is something like this:

The lighter color are the overlapped layer. Weird thing is that on the left, the overlapped buffer is rather incorrect, is it? Any idea what is going on?


texts: Using ModelBuilder to count number of times pixel has value higher than zero for more than three consecutive rasters?

texts:
I am quite new to the iterators within ModelBuilder. I have 2500 rasters representing a temperature in a region of Norway on a given day at 1 km resolution. I am interested in when the temperature is above freezing for three consecutive days (as then a layer of ice forms on top of the snow).

The raster files are in the format TM_YYYY_MM_DD. I have been trying to use Calculate field to add a value of either 1 or 2 onto the DD value, however this doesn't work since I don't have the files saved in a julian day format.

Is there a way where I can iterative through the folder and create one raster for each set of three temperature rasters where the pixel values have a value of 1 if the temperature was above zero for three days, and a value of 0 if it was below zero on one or more of these days?

Following on from this, is there a way I can count the number of these events per pixel?


texts:
If seen, that the Tasks is a good way to configure workflows for customers.

Is there a way to put a button to a group of a tab on the ribbon that starts a previously built task?


texts: QGIS Join issue on tables with generated unique value index field

texts:
I'm trying to make a simple join in QGIS 3.10 between a table with a generated ID (using "Add unique value index field") named POINTID and a table created by the Distance Matrix algorithm with an InputID field (whose values come from the POINTID).

But this simply doesn't work as I get an error message during the processing saying that all my features could not be matched…

First, it looked like a duplicate of Joining attributes by field value not working but I double checked the fields type and it is the same: int/integer(0,0) as you can see in both the layers properties and in the logs.
I also checked that IDs are common to both tables so that the join can be done.

Strange thing but maybe irrelevant: the length of the generated field is 0. But the documentation of this function doesn't say anything about setting a field length (see enter link description here). The way the values are sorted in the table from Distance matrix is strange too (1,10,11…,2,20,21,…).
Also, it looks like the join works well on fields of type string, but I didn't find a way to generate unique ID fields of type string (I'm using the graphical modeler).
This gets me confused, any ideas what I'm doing wrong ?
[EDIT] : here is some simple csv data to check the process (SRID is 2154)
Table A
STR_ID;X;Y;Z
P_0;311794;6644698;-0.885
P_1;311803;6644708;-0.689
P_2;311813;6644718;-0.516
P_3;311823;6644728;-0.337

From table A, I generate a POINTID field using X column as class for algorithm "Add unique value index field".
Table B (one point from which to compute distance)
STR_ID;X;Y;Z
P_REF;312556;6646432;-0.885

Table which is the result of Distance Matrix (distance values may be wrong but it doesn't matter):
InputID;TargetID;Distance
0;0;392.1
1;0;378.3
2;0;365.8
3;0;358.4


texts: Running ArcPy script on unix/linux server?

texts:
I've written some python scrips that use the ArcPy module to do some pretty intensive processing. I'd like to be able to run these on the fast unix server we have in our department so that I don't have to tie up my machine for hours running these.

Firstly: is this possible?

If it is, then do I need to install some part of ArcGIS on the machine (I can't seem to find ArcGIS Desktop for linux/unix), do I need ArcGIS server on there, or do I just need to put the python modules on there, plus something to let it find the license?

I've tried to research this but have ended up getting very confused about the various parts of ArcGIS and what does what!


texts: Joining (merge, union) non-overlapping vector polygons with QGIS

texts:
I look for a solution how I can join several polygons that do not overlap with each other and have no common sides, by creating a new polygon. That is, it should be an output polygon  covering the "outer" vertices of several input polygons, but in general repeating their shape (see attached figure, red line).
I tried several features in QGIS, in particular Convex Hull and Minimum bounding geometry, but it didn't work out what I needed.


texts: I am trying to access heatmap plugin from Raster -> HeatMap -> HeatMap using PyQGIS.If my custom button is clicked, the heatmap plugin dialog box should opens.I know the heatmap is a c++ plugin.Is it possible to trigger the heatmap plugin using Python?