texts:
I am trying to install ArcSDE 10 for PostgreSQL in CentOS. I tried doing that but am stuck with same errors repeatedly and am thinking of giving one more shot from the scratch. I just want to make sure the versions/releases before I start over again. I am thinking of using CentOS 5.4 (for i686 architecture), PostgreSQL 8.3.8/8.4.1, PostGIS 1.4.0 and ArcSDE 10. Please let me know if I have to change any version/releases of any software.

Thanks!


texts: Any possible way to convert cyrillic encoding to latin in QGIS?

texts: I'm struggling with QGIS showing certain layers for cities and other places in cyrillic writing (when the encoding is set to UTF-8) inside the attribute table. I'm wondering whether there's an encoding (or a plugin maybe) available to show that cyrillic writing directly into latin writing?


texts: Show ValueMap values based on value of another field in QGIS widgets

texts:
I'm applying a Value Map widget to two text fields in my layer.
The two fields are as follows:

I have a hierarchy consisting of Landcover, and under that Landuse, where the color is code for the available option. For example, Residential is only a Built-up landcover, and Open only has Grassland and Woodland to chose from.
I want to apply a widget to the Lower field, where based upon the selected value of Upper, it display only the available choices for Lower
How can this be achieved with QGIS Attribute Form widgets, or other Layer Properties configurations?


texts:
I want to draw a polygon on the world map in different map projections, for example, in mercator proj or stereographic one. They must correspond to ICAO ares G (https://aviationweather.gov/satellite?gis=off). But I got a trouble with geo coordinates. I see that my geo coordinates are differ than they should be. For example, in EPSG:3395 or EPSG:5588 my geocoords coincides in one point on the map, but in EPSG:4326 they are out of bounds of the world map. How can I draw geometry accordingly to geo coordinates? I am using GEOTOOLS.
Look at the point (-155.27,67.3623) -  it is on the north west of North America, but after reprojection this point appears somewhere in the Pacific Ocean


texts: Numeric digitise icon not in digitise bar

texts: I am using QGIS 2.4 on Mac Os.  I have the Numeric digitise plungin installed but it does not appear in the Digitise Tool Bar.  I have tried reinstalling it but it makes no difference.  I have used it in the past but not often so I have no idea when it disappeared.


texts: Markov chains, cellular automata, and raster algebra

texts:
I need to build a probabilistic cellular automaton for simulating land-use changes and am wondering how to implement this in a GIS. I understand that Markov chains are a way of doing this, and IDRISI has some sort of module for Markov chains, but I have no access to IDRISI. Someone said earlier that, with raster data, this can, in theory, be done with any GIS that can do raster algebra. Presumably this means that it can even be done in ArcGIS, and I think I've figured out how to build cellular automata in ModelBuilder (once I know with certainty how to make iteration work there).

Still, I can find next to no practical advice on this. Once I have two datasets, one showing the land use at a certain point in time (time A) and the other at a later time (time B), I imagine doing the following (correct me if I'm wrong):

If there are (for example) 3 types of land use, then there are 9 possible things that can happen to a cell from time A to time B; thus, reclassify each dataset such that each land use is represented by a different number (for example, 9, 11, and 14 at time A and 3, 4, and 5 at time B);
Divide one raster by the other, obtaining a raster with 9 possible cell values.
Determine the percentage of cells in the new raster that has each of the 9 values (I forgot, how does one do this in ArcGIS?); those will be the probabilities of each type of change (or lack thereof) in land use.

What I cannot figure out, though, is how to apply these probabilities during the next iteration. Finally, I wonder if there is a way of making an external factor (say, proximity to roads) influence the spatial distribution (not the overall likelihood) of change in my model.

Can anyone recommend where to find information on how to do this in ArcGIS?


texts: Error when using Raster Calculator in ArcGIS to compute NDVI, Process output is binary, limited to pixels that are either 0 or -1

texts:
I am using WorldView-3 satellite images that have been orthorectified and corrected to sensor reflectance values to compute a set of normalized difference indexes, including NDVI in ArcGIS version 10.31. While I am aware that it is possible to automatically compute NDVI in Image Analysis, I want to try the operation in Raster Calculator since NDVI is not the only index that I wish to compute.

However, whenever I run the operation (with any of the indexes that I am interested in using), the result is always the same. I obtain a binary raster output, in which pixels are either 0 or -1. The exact same pixels, no matter the type of index and bands that I choose. This is weird since when I manually calculate these indexes for the reflectance values in a set of control points, the result shows that the image should lead to an output that as expected ranges from 1 to -1.

The formula that I have been inputting into the raster calculator is the following.

((NIR2 - Red)/(NIR2 + Red))

To add more context, the site that I am working at is a peninsula in the north of Madagascar. The 0 values are always restricted to the land and the -1 values are restricted to the surrounding ocean and a small lagoon. The WorldView dataset is stored as a 16-bit TIFF file, containing the 8 multispectral bands that the satellite's sensor pics up. To run the raster calculation, I extracted the bands using ArcCatalog and opened them independently in ArcGIS.

Has anyone ever run into a problem like this?

Or does anyone have an idea of what is causing the problem?


texts: Polygons merge/union using two colums values

texts:
I'm using HydroSHEDS watersheds shapefiles for a proyect. Watershed are delinated applying the rules of Pfafstetter codification.

I need to merge/union the different polygons in a same shapefile. I need to select values for 2 different colums (values are the same, one is the HYBAS_ID and the other is the NEXT_DOWN) in the table of attributes, that way I can contruct the upstream basin. Also I need to keep some attributes and also sum others.

I'm using ArcGis 10.3.

Any idea in how can I do this?


texts: mapping number of days a site is flooded from daily inundated areas maps

texts:
I have a set of daily inundated area polygons and would like combine them into one map that can give me a map showing the number of days a site remained inundated.

I have ArcGIS Desktop.


texts: Speeding-up las files processing in a chunk with buffer

texts:
My aim is to do tree segmentation from lidar data. I have 15 las files and I want to merge them before segmentation. I realized that using chunk with buffer helps to get a proper segmentation.

When I use segmentation for each las tile, 1 meter gaps are occurring between the tiles. Therefore, I believe that before running segmentation codes, it is better to use catalog and chunk buffer. The problem is that my 2 laptops could not finish segmentation. It seems that it takes too much time and my laptops are not so powerful.

Is there anyway to merge las tiles and run the segmentation without any gap?

Here are the codes that I used in R:

ctg = readLAScatalog("path/to/catalog")
opt_chunk_buffer(ctg) <- 0
opt_chunk_size(ctg) <- 500
opt_output_files(ctg) <- "path/to/new/catalog/Forest_{ID}"
newctg = catalog_retile(ctg)
las <- readLAS(newctg, filter="-keep_class 2 5")
dtm <- grid_terrain(las, algorithm = knnidw(k = 8, p = 2))
las_normalized <- lasnormalize(las, dtm)
lasfilternoise = function(las, sensitivity)
{
  p95 <- grid_metrics(las, ~quantile(Z, probs = 0.95), 10)
  las <- lasmergespatial(las, p95, "p95")
  las <- lasfilter(las, Z < p95*sensitivity)
  las$p95 <- NULL
  return(las)
}
las_denoised <- lasfilternoise(las_normalized, sensitivity = 1.2)
chm <- grid_canopy(las_denoised, 0.5, pitfree(c(0,2,5,10,15), c(3,1.5), subcircle = 0.2))
ker <- matrix(1,5,5)
chm_s <- focal(chm, w = ker, fun = median)
algo <- watershed(chm_s, th = 4)
las_watershed  <- lastrees(las_denoised, algo)
trees <- lasfilter(las_watershed, !is.na(treeID))
hulls  <- tree_hulls(trees, type = "concave", concavity = 2, func = .stdmetrics)


texts:
Is there anyone who has some experience or knowledge on 3D visualization and/or GIS softwares capable of displaying 3D bathymetric data and 3D ship model together? Nearly all GIS softwares work with 3D data in grid format but I especially would like to learn whether it is also possible to put a 3D ship model on bathymetric model and dynamically visualize a ship's navigation by displaying its underwater sections.

Thanks,


texts: Transform coordinates in MapInfo with GSB data

texts:
Does anybody know if MapInfo (or MapBasic) is capable to do coordinate transformation using the NTv2 (National Transformation version 2) method?

Currently Switzerland changes their national coordinate reference system from CH1903 to CH1903+. For transformation the authoritiy provides a CHENYX06a.gsb file in NTv2 Format. Can I use this in MapInfo?

When I check the MapInfo Program folder it already contains a couple of .gsb files (but not the swiss one, CHENYX06a.gsb) and a file called NTv2.xml where you can see some reference data. So, in principle MapInfo should be capable to perform such transformation but I don't know how to trigger their usage.

When I try the normal conversion from CH1903 to CH1903+ the coordinates are shifted but the errors are up to 0.8 meters. Just for comparison, when I do the conversion in PROJ.4 using the CHENYX06a.gsb file, then the errors are in the range of Millimeters!

Appendix:

Looks like solution from Peter Horsbøll Møller does not work. Here the result of my test:

Easting     Northing    New Easting MI   New Northing MI   New Easting correct   New Northing correct    Delta meters
602030.68   191775.03   2602030.4        1191775.2         2602030.74            1191775.03              0.380
617306.30   268507.30   2617306.6        1268507.9         2617306.92            1268507.87              0.321
776668.10   265372.68   2776669.3        1265372.2         2776668.59            1265372.25              0.712
497313.29   145625.44   2497312.1        1145626.0         2497312.65            1145626.14              0.568
722758.81   87649.67    2722758.5        1087648.6         2722759.06            1087648.19              0.694

"correct" values are taken from Swisstopo sample data. In case somebody likes to reproduce the issue, here you can download the GSB-File: CHENyx06 dataset

Maybe one problem is, both coordinates are using the Datum (as far as I understand). In PROJ4 I use this string for transformation.

cs2cs -f "%.3f" +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 +x_0=600000 +y_0=200000 +ellps=bessel +nadgrids=CHENYX06a.gsb +units=m
    +to +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 +x_0=2600000 +y_0=1200000 +ellps=bessel +nadgrids=@null +units=m

The results are matching in range of Millimeters to the officaly values provided by Swisstopo.


texts: WMS - get point data in computer readable form (other than map tile)

texts:
is it possible to get a data from some map services on questions like

is there a water on [X,Y] (i.e. blue colour on map)
is there a dessert on [X,Y] (i.e. yellow colour on map)
similarly for any available layer

In language of computers, to know, where on some location (point,radius) is some layer present and more ideally, what are its properties (river=1,lake=0,waterfall=0,...)

I get get these data visually just by looking on some map tile presented on screen, but this is hard for computer calculations. I have read WMS reference http://docs.geoserver.org/stable/en/user/services/wms/reference.html#getfeatureinfo but I am not any more clever from that.

Thanks


texts: Intersect, Spatial Join?: One Line crossing Multiple Polygons, assigned wrong attribute

texts:
I am looking for any possible solutions to a problem I am encountering with ArcGIS.

I am using census block groups (CBG) to assign streets an income level value based on their location within the CBGs. I have tried using the intersect method which alleviates the problem but creates another in that it now splits the streets at the boundary of census block groups rather than at intersections as needed (need to keep them whole).

If I use spatial join, it keeps the street lines in their original form but in many instances (data with over 40,000 streets) they are being assigned the wrong attributes.

For example, in the images attached, a small portion of the street is located in a higher income CBG (yellow) with the majority being located a the lower income CBG (green). Using spatial join (match option: intersect), the street is assuming the attributes of the higher income (yellow).

Is there any method for running a spatial join (or other tool) that will assign attributes of a single polygon category (ex. Low income) to lines (streets) based on having the majority of their length in a single polygon when it is contained in two or more different polygons with varying attributes?

TLDR: (using image example): spatial join is assigning yellow polygon attributes to street though only a small portion of the line is within it. I am looking for a way to have attributes of any polygons joined to any lines which have their greatest majority within them (ex. green would be assigned to street not yellow) I could do this manually, however I need to ensure this will be done accurately for all 40,000 streets.


texts:
I'm trying to match small segments with a larger segment they are the most probably related to: relatively close, similar bearing, and facing each other.

Here is a typical example of the data I have:

Here I'd need to match segment 652 to 198969, while having 711 and 707 not matching anything.

I've looked for different methods, in particular the Hausdorff distance (based on the answers here). I computed it using PostGIS but I'm getting odd results: the shortest distance I get is between 707 and 198985, and 652 has a greater distance to 198969 than to 198985 for example (I can add the query and results if needed).

Is Hausdorff actually the correct method to solve this? Are there other approaches? I thought of simply creating a set of checks on the parameters I mentioned (distance, bearing, etc.) but I'm afraid of having to add a whole bunch of conditions to handle edge cases or things like thresholding on how much they're facing each other.

Update: I found a method which seems like an acceptable compromise:

I first find the 10 nearest black segments from the blue one I'm trying to match (using the PostGIS <-> operator) which are less than 10 meters away.
I then create a new segment by finding the nearest points to the ends of the blue segment on each of the black ones (using ST_ClosestPoint) and filter out the results whose length is less than 90% of the blue one  (meaning the segments aren't facing, or that the bearing difference is more than ~20°)
Then I get the first result sorted by distance and Hausdorff distance, if any.

There might be some fine tuning to do but it seems to do an acceptable job for now. Still looking for any other methods or additional checks to run if I missed some edge cases.


texts: Is ArcMap's Fixed Zoom In/Out algorithm documented anywhere?

texts:
I only occasionally use the Fixed Zoom In/Out buttons on the Tools toolbar of ArcMap and had always assumed that they zoomed in/out by a set scale factor.

However, using ArcGIS for Desktop 10.2, I just set my scale to be 1:10,000,000 and then hit the Fixed Zoom In button 10 times while recording the subsequent scale each time and was surprised to see:

10,000,000
8,000,000
6,000,000
5,000,000
4,000,000
3,000,000
2,500,000
2,000,000
1,500,000
1,250,000
1,000,000

Dividing each scale denominator by its successor shows a sequence of:
1.25, 1.33, 1.2, 1.25, 1.33, 1.2, 1.25, 1.33, 1.2

Curious about what happened in between the scales above I experimented a little to discover that the first issuance of the Fixed Zoom In button simply takes the scale to the first of those listed above and then starts following the same pattern.  Try starting with 1:8,000,001 to see ...

8,000,001
8,000,000
6,000,000
5,000,000
4,000,000
3,000,000

This seems like quite a neat solution, but not one that I had realised was in place - does anyone know if it has been documented anywhere?

I've looked in the Online Help without success.

In answer to @DevdataTengshe, the scales currently in my available list (see below) seem to bear no relation to those used by Fixed Zoom In/Out.

However, when I add 1:7,500,000 into my available list it does seem to insert that as an additional scale step ...

10,000,000
8,000,000
7,500,000
6,000,000
5,000,000
4,000,000

This makes the algorithm being used even more mysterious to me.


texts: How to get a list of features from Oracle via FME?

texts:
I have a table in Oracle containing only the names of other tables of varying feature types, which I'd like to process (i.e.export to csv). this task is to be scheduled via FME Server. the feature types to work on are only known at runtime. Therefore a test like Mapperz suggested would not work in this situation.

There are several examples how to do this based on names in text files or using cx_oracle (see http://fmepedia.safe.com/articles/How_To/Pass-a-List-of-Tables-or-Layers-to-an-FME-Reader-using-Python-Scripted-Parameters).

My situation is that cx_oracle isn't availabe and a text file won't do it as entries will change frequently (but the latter approach via scripted parameter is actually working fine).

The link above states in the beginning that: (...) For many situations, using an FME reader to read the table list and then the FeatureReader transformer to read the data may be a better approach than the one described below. (..)

I tried to follow this approach, which is also described by Mark Irland FME Features to read parameter, set up the reader for the table with the names of the other tables and connected it to a FeatureReader. It is configured to query the feature types specified in the attribute below (which is set to a concatenated attribute as described at "Feature Type List" in http://blog.safe.com/2011/05/fmeevangelist78/ or in another test case to the relevant column of the table ).
All ports are then connected to a shapefile-writer.

The result is that the entries of the table with the names are output via the initiator port and the features (with geometries) themselves to the 'other' port.
As result the geometries are shown in the inspector but aren't contained in the shapefile.

I also tried exposing fme_feature_type on the reader then taking it on the FeatureReader, but this makes only appear the feature names in the .dbf file with no geometries at all.

the hole thing is quite related to How to get a list of features from FME? Python? and Mark Ireland answer, and might be only causing troubles because Oracle is involved.

I also checked:
http://fmepedia.safe.com/articles/How_To/How-to-Read-All-Feature-Classes-from-Multiple-ESRI-Geodatabases / http://screencast.com/t/5OjFTrkiZL:
this approach won't work with oracle as you have to specify credentials.

so how could I get the FeatureReader to query the table names and put each in different shapefiles?

Unfortunatly the workbench code I used cannot be pasted here because of size limitation.

Added Image shows the workbench used: note in the lower left corner the concatenated result based on which the FeatureReader is supposed to query the database

Version information:
    #! LAST_SAVE_BUILD="FME(R) 2013 SP2 (20130605 - Build 13499 - WIN64)"


texts: Workaround for grouped layers to stay collapsed when opening project

texts:
Using QGIS 2.6.1 on Windows 7, whenever I open a project containing layers in various groups, they appear expanded. This is regardless if the project is saved when the layers are collapsed. I'm guessing this is a little bug but is there a workaround?

I am open to python, not open to a newer version of QGIS although it would be useful to know if this is solved.


texts: Defunct online mapping service with abstract or simplified route map

texts:
Back when it was common for people to print paper driving directions from a website like MapQuest (say, around 2000) I remember one of the mapping services offering an alternative format for the route map, which simplified each road to a straight line, in much the same way that subway maps often simplify the geometry of the routes.
In terms of the visual appearance on the page, this style of map always reminded me of a sentence diagram:

Does anybody remember what mapping service offered this style of map, and what they called it? Does anyone still have any examples, or know of a service that still offers this style of map?


texts:
I'd like to implement routing on my map using OSRM and Openlayers 3 but there don't seem to be any examples or tutorials anywhere. Everyone seems to be utilising Leaflet for this as there is a plugin available.

I understand that I would need to add a layer to hold the geometry and then make a request to viaRoute on the routing server but I'm not sure what form this request should take. I know it requires parameters to be passed in the request in this format:

http://router.project-osrm.org/viaroute?loc=52.503033,13.420526&loc=52.516582,13.429290&instructions=true

My javascript is still a little patchy so I'm not sure how to construct the request and response to be passed to my geometry layer.

Can anyone provide an example of this or perhaps give some guidance as to how it can be achieved?


texts: Can't create PNG file with "st_png" function PostGIS and Java

texts:
I have this code in Java:

conn = DriverManager.getConnection ("jdbc: postgresql: // localhost: 5432 / postgis_23_sample" "postgres", "admin");

conn.setAutoCommit (false);

String s = "SELECT ST_AsRaster (ST_Buffer (ST_GeomFromText ('LINESTRING (50 50,150 150,150 50)'), 10), 200,200, ARRAY ['8BUI' '8BUI' '8BUI'] ARRAY [118,154,118] ARRAY [0 , 0,0]) ";

PreparedStatement sGetImg = conn.prepareStatement (s);
SGetImg.executeQuery ResultSet rs = ();
FileOutputStream fout;

try
{

    rs.next ();
    fout = new FileOutputStream (new File ("d: \\ files \\ save1x.png"));
    fout.write (rs.getBytes (1));
    fout.close ();

}

catch (Exception e)
{

    System.out.println ("Can not create file");
    e.printStackTrace ();

}

rs.close ();
sGetImg.close ();
conn.close ();

catch (SQLException se)
{

    System.out.println ("Could not connect: print out a stack trace and exit.");
    se.printStackTrace ();
    System.exit (1);

}

....

I created the file but apparently the PNG file is corrupt, you can not open. That I can be doing wrong.

I reviewed the bytes and the query returns the array of bytes me, but does not write the file apparently well.

Postgres 9.4 y PosGis 2.2


texts: How to figure out the unit of measurement from table fields?

texts: I have a couple of fields with predefined values. In particular I have an 'area' field that I need to figure out its unit of measurement (e.g. km, miles, acres). What's the easiest way to do this?


texts: Change coordinate system in multiple MXDs ArcGIS 9.3 using Python

texts:
I have to loop through some folders which may contain *.mxd files and change their coordinate system. I already found this piece of Python which should do the trick for ArcGIS 10.x (from here):

prjFile = prjPath + "/" + "NAD 1983 UTM Zone 10N.prj"
sr = arcpy.SpatialReference(prjFile)
for DF in MAP.ListDataFrames(MXD):
    DF.spatialReference = sr

However, it looks like there is nothing similar to the arcpy.mapping module in 10 available in 9.3. I need to access the MXD file and then the dataframe(s) it contains, in order to change the coordinate system.

I could not find anything useful so far, although I assume that other people may have similar issues.

My question therefore, is there any (simple) workaround to this using Python (arcgisscripting) in ArcGIS 9.3?


texts: ArcGIS Pro line symbology: adding gap in line where a shape marker is

texts:
I have a created a line symbology template with shape markers at regular intervals along the line, as shown in the picture below.

I want to make it so that where there is a shape marker there is a gap in the line, like this:

Is there a way to do this?


texts: Geocoding addresses based on own reference data

texts:
I have a historical addresses dataset in form:

|Name   |Description  |

Where "Name" is a single field with house number, street name, city and country.  I have digitised two georeferenced historical maps of the location to cover the period and account for street name or geometry changes and for the addition of or removal of streets compared to present day.

This gives me point data for individual addresses and attributed polylines with geometry for the streets in the city:

AddressName: "59 Great Strand Street, Nulltown, GISland"
DigitisedMapName: "Great Strand Street, Nulltown, GISland"

I now wish to join these up matching based on the street name, giving the street geometry to the address data and plotting the address a certain % along the street based on the house number.

I have tried the ArcMap 10.2.2 address locator but get 0% matched leaving only a manual selection (not feasible given the number of records).

Is it possible to complete this task in PostgreSQL instead to strip the mapped geometry from that table and give it to the AddressName record?  Or are the wheels going to fall off in the joining accuracy?


texts:
What is the best way to look at and query data (e.g. group by) in a file geodatabase?

I'm thinking something similar to opening up a Personal Geodatabase in Access, creating queries, adding columns, etc.

The data I'm working with is too large and way too slow to query in Personal Geodatabases, but I still want to be able to view the data outside of ArcGIS and have the updates relayed back to GIS when I work with it there.

Am I better migrating the data to a proper database?

I'd rather not have to do this.


texts: Using ModelBuilder to replicate series of tasks without deleting output features first?

texts:
I am calling a number of files that I want to clip etc in ModelBuilder which clips data and sends it to a specific feature class.

I want to re-run that same process for which I have converted the model to python and I am calling from Scheduled task in Windows 7 to update data in an application that staff use.

What I would like to know is if the feature classes I am clipping have changed, will it automatically do the clip, or do I have to delete the feature classes that have already been created first, and then re-do the Clip?


texts:
Is there a way to print an atlas in a pdf/A file for each page, being pdf/A an international standard?

I'm using QGIS 2.18.


texts:
I am trying to reclassify the CORINE map for Hungary to 4 classes in R. but after reclassification, when I import the asc file to QGIS it loads it empty. Am I doing something wrong?

x<-scan("C:/Users/HP/Documents/Corine_Hungary/Corine_Hungary_2006.asc", sep=" ")
row=length(x)
l <- 1:row
for (i in seq(along=l)){
      if (is.na(x[i])){
      x[i]=0
    } else if (x[i]==0){
      x[i] = 0
    }else if (x[i]<4){
      x[i]=1
    }else if (x[i]==5){
      x[i]=1
    }else if (x[i]==4){
      x[i]=2
    }else if (x[i]<10){
      x[i]=2
    }else if (x[i]<30){
      x[i]=3
    }else if (x[i]<35){
      x[i]=2
    }else if (x[i]<45){
      x[i]=4
    }else if (x[i]>45){
      x[i]= "NODATA"
    }
  }

x=t(x)
write.table(x, "c://Users/HP/Documents/Corine_Hungary


texts:
I want to create a template_postgis for connecting postgis via qgis

So I follow the step on following page

how to create a new postgis database for usage with QGIS under windows

But when I follow the steps there's a problem happened

createdb: database creation failed: ERROR:  template database "template1" does not exist

So I google it for answer but no one can help me to create a template1 database!

Any suggestion?


texts: Hibernate Spatial - 'Invalid endian flag value encountered' Exception

texts:
I'm trying to run a simple query in Hibernate Spatial 4.0 on PostgreSQL 9.3.  I have a number of objects in a table with latitude/longitude values, and I'm trying to query objects that fall within a given radius of a particular latitude/longitude.  The geometry values seem to be persisted without any problem, and are defined like this in my entity class:

@Column(columnDefinition = "Geometry", nullable = true)
@Type(type = "org.hibernate.spatial.GeometryType")
private Point coordinates = null;

I don't have any errors when persisting objects with the coordinates value set.  However, when I run a query, I see the following Exception:

javax.servlet.ServletException: javax.servlet.ServletException: javax.ejb.EJBTransactionRolledbackException: org.hibernate.exception.GenericJDBCException: could not extract ResultSet

<snip />

org.postgresql.util.PSQLException: ERROR: Invalid endian flag value encountered.

The query I'm running looks like this:

Query query = this.entityManager
        .createQuery(
                "SELECT v FROM MyEntity v WHERE within(v.coordinates, :filter) = true",
                MyEntity.class);
query.setParameter("filter", point);

I came across this question on StackOverflow that mentions the same error, but I'm declaring my coordinates property the same way that is mentioned in the answer and still get this error.

For the record, I've tried using PostGIS 1.5.2 and PostGIS 2.1.0.  I've also tried different version of my PostgreSQL JDBC driver, from 8.4 to 9.3.  Regardless of the library version I still encounter this issue.

Can anyone shed any light on what could be happening here?  Is my query wrong?  Is my property not defined correctly?  Is there something else I should be trying?  I'm completely stuck and at a loss for what could be the issue here.  Any advice would be GREATLY appreciated.


texts:
After invoking:

shp2pgsql -W "latin1" theshape.shp > thesql.sql

My console responds with:

Shapefile type: PolygonZ
Postgis type: MULTIPOLYGON[2]

Note the shapefile type is PolygonZ. The result is, I have to manually edit the thesql.sql file and wrap the geometry column in ST_Force_2D() before executing the SQL, otherwise, Postgres complains I'm trying to import a 3D shape into a 2D column.

How can I force dimensionality to 2 with shp2pgsql?
Alternatively, can I force it in QGIS?


texts: Whether NOAA AVHRR NDVI V5 already remove cloud?

texts:
We want to use NOAA/CDR/AVHRR/NDVI/V5 datasets for analysis, but we are not sure if the NDVI already remove the cloud and cloud shadow or not?
Need we use 'QA' to remove the cloud and cloud shadow further?


texts: How do I edit an individual label in QGIS that has data from two columns on two lines, but also needs a word wrap?

texts:
So I have a complex question. Initially I needed to know how to do a word wrap on a single label in QGIS because I have a county name that is too long for the county it belongs to. I was told how to do a word wrap function.

However, my labels also have data from a second column in the attribute table that goes beneath the county name.

My question is how do I combine both functions in the expression so that the one county with its long name word wraps, but all the labels also get the data from that second column? For reference, here is the link to the first question I asked: How do I edit a single label in QGIS 3.4 (Madeira)?


texts: Accuracy assessment of RFclassificatiton in dzetsaka plugin QGIS

texts:
I tried random forest classification on dzetsaka plugin in QGIS to discriminate different plant species using sentinell(10m) and planet(3m) image. Classification is performed but i am confused about its accuracy. How can i access the accuracy of this classification?
Thanks


texts: r.viewshed.cva error: No points found, nothing to be exported

texts:
I need to do an intensive viewshed and I am going to use our university SCC (Linux).

I am running a test viewshed with only 13 observers (just to build the code - eventually the script will be used on an OP shapefile with 700k points).

I've imported them into GRASS, and they are overlaying properly on the DEM. The input is simple (1.65m offset of the observer, no targets, 0.13 refraction. I've used these same parameters in ArcMap and it works fine.

Here is the code:

r.viewshed.cva -c input=pelop1arcm34b@PERMANENT vector=PI_14@PERMANENT output=test observer_elevation=1.65 memory=1000 refraction_coeff=0.13

For the output, I get this error:

WARNING: No points found, nothing to be exported
Note that the routine is constrained to points in the current geographic region.
Calculating "Cumulative Viewshed" map
ERROR: Raster map not found
Traceback (most recent call last):
  File "C:\Users\nsusm\AppData\Roaming\GRASS7\addons/scripts
/r.viewshed.cva.py", line 181, in <module>
    main()
  File "C:\Users\nsusm\AppData\Roaming\GRASS7\addons/scripts
/r.viewshed.cva.py", line 169, in main
    input=(",").join(vshed_list), output=out,
method="count")
  File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.4.0\etc\py
thon\grass\script\core.py", line 421, in run_command
    return handle_errors(returncode, returncode, args,
kwargs)
  File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.4.0\etc\py
thon\grass\script\core.py", line 332, in handle_errors
    returncode=returncode)
grass.exceptions.CalledModuleError: Module run None
['r.series', '--q', 'input=', 'method=count', 'output=test']
ended with error
Process ended with non-zero return code 1. See errors in the
(error) output.


texts: How to remove black box around georeferenced grayscale image in QGIS?

texts:
i'am trying to georeference a grayscale image (*.tif) in QGIS 2.2.0 Valmeria. It is working but QGIS draws a black box around the georeferenced raster image. If i use the option set 0 transparent, not only the black box is transparent but also the black pixels in my gray scale image. But setting black or white transparent isn't an option for me, because it is a gray scale image.

Thanks a lot for you help!

Rob


texts: Unable to connect to postgres server in Database Manager

texts:
I try to connect to my postgres server, hence I got this failure

And following error:

2018-04-26T18:02:05 WARNING Traceback (most recent call last):
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_model.py", line 441, in rowCount
            self._refreshIndex(parent, True)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_model.py", line 490, in _refreshIndex
            if item.populate():
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_model.py", line 167, in populate
            if not connection.connect():
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\postgis\plugin.py", line 103, in connect
            return self.connectToUri(uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\plugin.py", line 109, in connectToUri
            self.db = self.databasesFactory(self, uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\postgis\plugin.py", line 70, in databasesFactory
            return PGDatabase(connection, uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\postgis\plugin.py", line 111, in __init__
            Database.__init__(self, connection, uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\plugin.py", line 218, in __init__
            self.connector = self.connectorsFactory(uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\postgis\plugin.py", line 114, in connectorsFactory
            return PostGisDBConnector(uri)
          File "C:/OSGEO4~1/apps/qgis/./python/plugins\db_manager\db_plugins\postgis\connector.py", line 128, in __init__
            os.remove(sslCAFile)
        PermissionError: [WinError 5] Access is denied: 'C:/Users/Matrix/AppData/Local/Temp/tmp_basic_{453d7688-e4d9-44e5-b2c5-51860a063061}.pem'

Any suggestions?


texts: PyQGIS: Properly trigger actionPan() tool when feature is added

texts:
I wrote a plugin which creates a custom map tool that adds features. I would like to cleanly connect to the pan tool after a feature is added. I've done it in the following way and it works, but I get the QMessageBox an increasing amount of times as I use the plugin. Any suggestions on a cleaner/better way to do this?

def initGui(self):
    self.action.triggered.connect(self.run)
    self.moreAction.triggered.connect(self.runWithMoreStuff)
    self.clickTool.connect(self.addToMyLayer)

def setPanTool(self):
    self.iface.actionPan().trigger()
    QMessageBox.information(None, 'debug', 'pantriggered')

def addToMyLayer(self):
    self.myLayer.addFeature().trigger()

def addToAnotherLayer(self):
    self.anotherLayer.addFeature().trigger()

def run(self):
    self.canvas.setMapTool(self.clickTool)
    self.myLayer.featureAdded.connect(self.setPanTool)

def runWithMoreStuff(self):
    self.canvas.setMapTool(self.clickTool)
    self.myLayer.featureAdded.connect(self.addToMyLayer)
    self.anotherLayer.featureAdded.connect(self.setPanTool)

The idea is if one button is pushed, i'll add just the one feature then set the pan tool. If the other button is pushed, i'll add a feature, then another feature on another layer, then set the pan tool.


texts: I cannot choose the preference layer after I add my GDB file in QGIS

texts:
After I add my Geodatabase file into QGIS. Usually, I can select the preference layer that I want to use. (See the attached picture). [This box usually appears after I add my GDB file where all feature classes in the gdb will be listed, number of features, IDs, and geometry type]. However, there is no box anymore. So, I could not select the preference layer from the GDB file.

I use QGIS 3.10. Can anyone help me to fix this?


texts:
I wish to have a Online 3D Map, which would show my Data, i.e both Raster and Vector, draped on a DEM that I have.

Due to the requirements of using my own DEM & Raster data, Google Earth API is pretty much out of the running.

Are there any projects which can do this?


texts:
I use ArcGIS 10.1 & arcobjects.net.

I have the following code:

        ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
        gp.OverwriteOutput = true;

        ESRI.ArcGIS.ConversionTools.FeatureClassToShapefile tool = new ESRI.ArcGIS.ConversionTools.FeatureClassToShapefile();

        tool.Input_Features = @"F:\test buffer of JAMNAGAR SHP FILES ex3\FOOTPATH.dwg";
        tool.Output_Folder = @"C:\Temp1";
        gp.Execute(tool, null);

but after executing this code I got an error

  "Error HRESULT E_FAIL has been returned from a call to a COM
  component. "

If anyone have any suggestion or any sample for this conversion or any link then send me please


texts: QGIS 1.6.0 problem during printing & export to PDF

texts: When trying to print or export my map to pdf something happens to my document.  When trying to print, the outcome always results in the printer producing a solid grey rectangle, and when I try to export to pdf the file becomes damaged and will not open.  I'm trying to print a 36x24 in map if the size makes any difference.  Has anyone had a similar problem?


texts:
I am new of pgrouting.

I have created a new network using osm2po with both edges and nodes.
Everything works fine, and I have tested the system using the Pgrouting Layer plugin for QGIS but I have noted that I am not able to use a table of POI that are not intersecting my network.

I am not able to figure out how to add these POI as nodes, is there any easy way to solve this? I have already checked the workshop manuals and other discussions, but it seems to me that they don't explain how to do this operation, or, better, I understand that the POI should be on the edge in oder to split the lines before to create the table with the nodes, but if I have a table with hundreds of POI (e.g. all the bar/pubs in a city), and all these points are not in the edge, is there a way to fix it automatically?

Using QGIS 2.99 I see that there are some interesting tools for network analysis, if I try the shortest path algoritm using my table of POI I have this result

I would like to reach the same result using Pgrouting.

Is it possible?


texts: Selecting contour lines with higher elevations on either side

texts:
I want to select contour lines that have lower elevation values from neighboring contour lines.  I'm using ArcGIS but do not have the hydrology tool set.  I can use open source if needed.  I thought about cutting the contour lines out of a road buffer, and selecting from the cut contour layer using sql query by tping something to the effect of: elevation value (n) is less than elevation value
(n plus or minus one.) I however, don't think this type of query can be done, and I don't know how to make the contour lines part of a network to perform spatial sequential analysis. If I did I would create the network using the roads layer.  I hope that makes sense.


texts:
I am working wit a win10 64 bits 16mb of ram. I want to install Lyon on my computer and not get rid of 2.83. My 2.83 is a stand allone install. I thought I could install Lyon usion OSGEO 4w installer. Nevertheless, even though I tried several times, I can not get acces to the grass plugin. I can not found it on the plugin manager nor I can found a toolbar. I tried to uninstall the OSgeo4w install several times, ask for to isntall Lyon, grass and Saga and no luck at all. I found at http://mappinggis.com/2015/08/novedades-del-plugin-grass-en-qgis-2-10-pisa/ (in spanish) that Lyon handles grass differently than 2.83 and that in order to have access to Grass I must install it using the stand allone version.

How do I a Install Lyon with Grass support using the OSGEO  installer and not the standalone?


texts: Is the PostGIS Raster primary key required to be named rid?

texts:
In all examples from the PostGIS Raster reference the primary key in tables with a raster column is named rid. I wonder if that naming is expected by any functionality of PostGIS itself or any client software connecting to PostGIS (such as mapserver for instance).

In other words:

Is the name rid for the primary key in raster tables a requirement, or can I use id as a primary key without running in to trouble?


texts: CartoDB: Info window from multiple layers, one geo, one supporting data

texts:
I want to have one table (readonly_table) with supporting non-geo data (e.g. project code, name, manager, etc.) and another table (write_table) with geo point, project code, perhaps a few other user entered properties e.g. start date.

And then I want to have these tables joined by project_code as key in a visualization so the info window displays both tables' fields on refresh (assuming readonly_table has a matching project_code).

Example: readonly_table has project codes 1, 2, 3, along with their other field values)  The write_table is empty.  User goes into the map, clicks "add feature" on the write_table layer.  Adds a point.  Clicks the "edit data" option.  Types in a project code (2). Hits save and close.

Now when the user clicks on the same point on the map, he sees one info window with the merged results of write_table fields (project code, start date) and readonly_table fields (name, manager, etc.)

Is this doable within the CartoDB dashboard?  If not, can you show a very basic example of an info window displaying from two tables based on a common field join (and geo data from one of the tables)?


texts: add description to a new polygon and save -> postgres + geoserver + openlayers -> wfs-t insert

texts:
actually it is possible to draw a new polygon with openlayers and save it over WFS-T-insert (geoserver wfs-t 1.1.0) in my postgres-db:

var save = new OpenLayers.Control.Button({
title: "Save Changes",
trigger: function() {
    if(edit.feature) {
    edit.selectControl.unselectAll();
    }

var win = new Ext.Window({
    title: 'Beschreibung des Fehlerreports',
    html:'<br>',
    html:'Bitte beschreiben Sie den aufgezeigten Fehler in der Darstellung:',

    layout:'fit',
    width:500,
    height: 300,
    modal: true,
    autoScroll: true,
    shadow: false,
    closable: true,
    buttons: [
    {
        scope: this,
        text: 'Speichern',
        handler: function() { win.hide(); saveStrategy.save(); }
    },{
        scope: this,
        text: 'Abbrechen',
        handler: function() { win.hide(); }
    }]
});
win.show();
//alert('saved');
},
displayClass: "olControlSaveFeatures"
});

In the saveing-window i need 3 forms to store some information (character) together with the polygon in the DB. This is the structure of the DB and the the published layer from geoserver:

DB - Postgres:
http://dl.dropbox.com/u/19066498/db.jpg

Geoserver - public Layer:
http://dl.dropbox.com/u/19066498/geoserver.jpg

Any idea what is a possible way to store these informations together with the wfs-t insert post??


texts: Getting StatPlanet to work with custom shapefile and data?

texts:
I'm assigned to get StatPlanet Plus working with some 3rd-party data to display a thematic map, on a Mac.

I have the custom regions and attribute data in a map.shp/dbf. And I've configured the settings.csv and data.csv according to the manual instructions in the User Guide. (There is some doubt about whether the Excel macros will work properly, especially on a Mac, so I've elected to do the configuration manually in order to eliminate unknowns.)

That guide has good instructions on what to do, but it's hard to find any complete, simple examples of the data.csv, settings.csv, map.shp and map.dbf files (other than the complex world example in StatPlanet_Shapefile_(ESRI)). Most of the maps at http://www.statsilk.com/maps/download-free-shapefile-maps do not include shapefiles. I eventually found that the map for Bhutan did. I used that one as a template for my own data, making only minimal changes to get it to work with my shapefile -- but alas, it still doesn't work.

What works for me:

The StatPlanet_Shapefile_(ESRI) map out of the box, with its many indicators in several categories, shows the countries of the world, colored according to the selected indicator.
The StatPlanet_Shapefile_(ESRI) map, with the shp/dbf/settings/data for Bhutan that I downloaded.
The StatPlanet_Mac app for shapefiles, with my own map.shp and map.dbf dropped in, but without modifying data.csv/settings.csv: This displays the regions from my shp, but of course doesn't color them according to our data (I wouldn't expect it to).

What doesn't work:

The StatPlanet_Mac app for shapefiles, with my own map.shp and map.dbf, and with my modified data.csv/settings.csv: This displays the StatPlanet splash screen, and nothing further. No error messages are given (even when I run the app from the command line). I imagine there's something wrong with the data.csv and the settings.csv, but as far as I can tell, I've done everything according to the instructions.

Here are the parts of those files that I've modified from the ones that came in the box:

settings.csv

Set the ID column name next to DBF-ID to GNS_ADM1
Put in the names and IDs of provinces of Suriname at the bottom under NAMES, e.g. "Nickerie   NS14".

data.csv

Deleted all rows below the first indicator (row 4)
Put in my own numbers in row 4, columns L through U
Put in ID's in row 1, columns L through U (these are values from the GNS_ADM1 column, e.g. NS19)

If you want to look at my file contents in detail, here they are:

settings.csv
data.csv
map.shp
map.dbf

And in case anyone asks, yes, I copied the map.dbf and map.shp into both the map/ and web/map/ folders.

So my question boils down to...

Can anybody see anything wrong with how I've set up my settings.csv or data.csv?
Does anyone know how to get StatPlanet Plus to output error messages?

So far I've only been using the desktop version of StatPlanet rather than the web version, since it wasn't clear how to use the web version with your own region vector data (and the web also introduces another layer of complexity unnecessary for this assignment).


texts: Geometry converts from Polygon to Multipolygon in PostGIS

texts:
I'm having difficulties with the Geometry types in PostGIS and/or ArcPy and hope some of you can shed a light on what I'm missing or how I can fix it.
I have a .gdb containing Polygons. The Table in ArcGIS lists the features as Polygons and using ArcPy using the following code I get the Output "Polygon" as well.
desc = arcpy.Describe("Path to Featureclass")
geometryType = desc.shapeType
print(geometryType)

Now when I upload the featureclass using ogr2ogr:
ogr2ogr -f "PostgreSQL" “PG:host=[host] port=[port] dbname=[dbName] user=[user]" “[PathtoFileGDB]” -overwrite -skipfailures –progress

And then access it via psycopg2:
import psycopg2
from shapely import wkb

conn = psycopg2.connect(host...)
conn.autocommit = True
cursor = conn.cursor()

cursor.execute('''select * from testfeatures;''')
record = cursor.fetchall()

conn.commit()
conn.close()

I'll get Multipolygons as a geometry. I already ran the "Repair Geometry" Tool in ArcGIS Pro before uploading the .gdb and the features display fine at every turn...
This sums up my problem, if you need any more information don't hesitate to ask and I'll try to provide it.


texts: Different values between a windowed read and a normal read in rasterio

texts:
I am struggling to understand why I obtain different results using the same coordinates in this example:

import rasterio
#data link: https://earthexplorer.usgs.gov/download/external/options/SENTINEL_2A/3022286/INVSVC/

pathToImgFolder = ""
pathData = pathToImgFolder + "L1C_T21LXC_A001666_20170701T140052\S2B_MSIL1C_20170701T140049_N0205_R067_T21LXC_20170701T140052.SAFE\GRANULE\L1C_T21LXC_A001666_20170701T140052\IMG_DATA\T21LXC_20170701T140049_B01.jp2"
boudingBoxCoordinates = (606859.0750363453, 8241169.269917269, 607219.0750363453, 8241529.269917269)
points = [(607014.0750363453,8241374.26991727),(607024.0750363444,8241374.26991727),\
          (607034.0750363453,8241374.26991727),(607044.0750363453,8241374.26991727),\
          (607054.0750363442,8241374.26991727)]

with rasterio.open(pathData) as img:
    bandWindow = rasterio.windows.from_bounds(*boudingBoxCoordinates, img.transform)
    winTransform = rasterio.windows.transform(bandWindow,img.transform)
    bandData = img.read(1, window = bandWindow)
    allData = img.read(1)
    for point in points:
        rBand,cBand = rasterio.transform.rowcol(winTransform,point[0],point[1])
        rFull,cFull = img.index(point[0],point[1])
        bandVal = bandData[rBand,cBand]
        fullVal = allData[rFull,cFull]
        print(f"{fullVal} {'=' if bandVal == fullVal else '!='} {bandVal}")

When I execute this example i get this result:
1202 = 1202
1330 != 1202
1330 != 1202
1330 = 1330
1330 = 1330

I think it is related to the spatial resolution of the band(60m), but I am not sure and I do not know what to do.
The img.profile is:
{'driver': 'JP2OpenJPEG', 'dtype': 'uint16', 'nodata': None, 'width': 1830, 'height': 1830, 'count': 1, 'crs': CRS.from_epsg(32721), 'transform': Affine(60.0, 0.0, 600000.0, 0.0, -60.0, 8300020.0), 'blockxsize': 192, 'blockysize': 192, 'tiled': True}


texts:
I have a shapefile with polygons of beaver ponds and a DEM of the area of river watershed where ponds are located.

I want to check if a river gradient and a river valley shape effect area of a pond.

I've created graphs of a river vertical profile and river cross sections in pond locations using 3D analyst of ArcGIS but I want to see my ponds on those graphs.

Is there a way to insert the pond data into the created graphs or some other way to show on one graph a pond and river profile?


texts: QGIS Graphic Modeler issues with Raster Calculator

texts:
A colleague and I have downloaded QGIS Dufour and are trying to use the Graphical Modeler to carry out a raster subtraction with a DSM and a DTM. This fails using both the Grid Difference and Raster Calculator tools with a message box saying the output file could not be opened.

If the same operation is attempted using the Raster Calculator outside of the Graphical Modeler it works fine and a good result is achieved.

Has anyone got any ideas or a solution?

Thanks,

Ed


texts: Map algebra from files in multiple directories

texts:
I'm a rookie in python programming.
I want to do map algebra for multiple rasters but these rasters are in different folders. Here's the code ive created so far:

#Importing modules
import arcpy
from arcpy import env
from arcpy.sa import *

#Files name
name1 = "BR_CRU_PRE"
name2 = "BR_CRU_ETP"

#Workspace
path = r"C:\Users\Ibytee\Documents\FUNCEME\ParaGiullian\Dados_ArcGIS\PRE_ETP"
env.workspace = path

#Checking for extension
arcpy.CheckOutExtension("Spatial")

#Creating raster lists
folders = ["PRE","ETP"]
env.workspace = path +"\\"+ folders[0]
a = arcpy.ListRasters(name1+"*", "TIF")
env.workspace = path +"\\"+ folders[1]
b = arcpy.ListRasters(name2+"*", "TIF")

#Executing the loop
for p, etp in zip(a, b):
     print p
     outRaster = Raster(p) - Raster(etp)
     outRaster.save(path+"\\PRE_ETP_SUB\\"+"sub_"+name1+"_%s" %etp)
     print "sub_"+name1+"_%s" %etp

This script seems to be accessing the right folders and creating the raster list, but when it comes the time to do the map algebra, looks like the script searches the rasters in the wrong folder. Here's the error I get:

C:\Python27\ArcGIS10.1\python.exe C:/Users/Ibytee/Documents/eclipse_workspace/Scripts/teste.py
BR_CRU_PRE_ANO1971.tif
Traceback (most recent call last):
  File "C:/Users/Ibytee/Documents/eclipse_workspace/Scripts/teste.py", line 24, in <module>
    outRaster = Raster(p) - Raster(etp)
RuntimeError: ERROR 000732: Input Raster: Dataset BR_CRU_PRE_ANO1971.tif does not exist or is not supported

Process finished with exit code 1


texts: Generating raster (spec. GeoTiff) out of mesh file (.ply)?

texts:
I created a mesh file using OpenDroneMap program out of 100 geotagged jpgs, and I would like to know how to generate a raster image out of the resultant .ply file?

Has anyone done this?

Meshlab has the option to make a screenshot but it returns very poor quality image. I saw that Agisoft Photoscan has this ability as well, but are there any free software packages that can do this?


texts: Develope ArcGIS Desktop Add-Ins with ArcGIS Runtime SDK

texts:
I need to learn develop add-in in ArcGIS 10.3. I have this book :

http://www.amazon.com/Beginning-ArcGIS-Desktop-Development-using/dp/1118442547

but he was using visual studio 2010 and ArcGIS Runtime SDK 10.1.
I have Visual Studio 2013 and ArcGIS Runtime SDK 10.2.5

He used ArcGIS Add-Ins Wizard. but i can't find Add-Ins Wizard in visual studio 2013. How can i develop Add-In in SDK 10.2.5 ?


texts: QGIS Problem with extract and warp modules in modeler?

texts:
I have problems using the gdal modules extract and warp in the modeler. They run correctly when started without the modeler, but when I integrate them in the modeler they won't work.
The only input I do set is the raster layer.

The process ends with the known remark:

  "Oooops! The following output layers could not be open
  Output Layer: C:\Users\sven\AppData\Local\Temp\processing\d2999fed6cd643a280f809eb0bee2373\OUTPUTALG0.tif
  The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm
  Checking the log information might help you see why those layers were not created as expected"

The Log states everything is OK and ends with:

  "Executed 1 algorithm total.
  Converting outputs.
  Loading resulting layers"

(If more information is needed ask for it)

The ogr-protocol states that the "data source is invalid"

As mentioned everything works fine, if I dont use the modeler.


texts: Export stacked layers of images as a thumbnail URL in GEE

texts:
I am computing different normalized index like NDVI, NDWI, NDBI, UI and adding them as layers stacked on top of Sentinel 2-1C Image patch.

I would like to download the generated stacked image as a thumbnail (jpg or png) file.

I can do image.getThumbnailURL on each layer and extract them as separate thumbnail images, but I want to export the entire stacked-layer as a single image.

Code Snippet

var LAT = 18.7332794
var LNG = 82.811036
var geometry = ee.Geometry.Point([LNG, LAT]).buffer(10000).bounds()

function mask2clouds(image) {
  var qa = image.select('QA60');

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

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

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

var collection = ee.ImageCollection('COPERNICUS/S2')
              .filterDate('2015-06-01', '2015-12-31')
              .filterBounds(geometry)
              .sort('CLOUDY_PIXEL_PERCENTAGE').limit(5)
              .map(mask2clouds)
              .select(['B4', 'B3', 'B2', 'B8', 'B11', 'B12'])

var median = collection.median()
var roi = median.clip(geometry)

// NDWI: NIR - GREEN
var ndwi = roi.normalizedDifference(['B3', 'B8'])
var ndwiMasked = ndwi.updateMask(ndwi.gte(0.2))

// NDVI: RED - NIR
var ndvi = roi.normalizedDifference(['B8', 'B4'])
var ndviMasked = ndvi.updateMask(ndvi.gte(0.2))

// NDBI
var ndbi = roi.normalizedDifference(['B11', 'B8'])
var ndbiMasked = ndbi.updateMask(ndbi.gte(0.05))

// BARE
var bare = ndvi.lt(0.2).and(ndwi.lt(0))
var bareMasked = bare.updateMask(bare)

Map.setCenter(LNG, LAT, 12)

Map.addLayer(roi, {bands: ['B4', 'B3', 'B2'], min: 0.0, max: 0.3}, 'roi')
Map.addLayer(ndwiMasked, {'min': 0.5, 'max': 1, 'palette': ['cf92c6', '67001f']}, 'ndwi')
Map.addLayer(ndviMasked, {'min': 0.5, 'max': 1, 'palette': ['98d7c7', '00441b']}, 'ndvi')
Map.addLayer(ndbiMasked, {'min': 0.5, 'max': 1, 'palette': ['fec253', '662506']}, 'ndbi')
Map.addLayer(bareMasked, {}, 'bare')

I am attaching the link to code editor for reference to what I am trying to achieve.

Edit

I am able to concatenate the ndvi, ndwi, ndbi normalized images and extract them as a thumbnail.

Code Snippet

// NDWI: NIR - GREEN
var ndwi = roi.normalizedDifference(['B3', 'B8'])
// NDVI: RED - NIR
var ndvi = roi.normalizedDifference(['B8', 'B4'])
// NDBI
var ndbi = roi.normalizedDifference(['B11', 'B8'])
// BARE
var bare = ndvi.lt(0.2).and(ndwi.lt(0))

var multiband = ee.Image.cat([ndvi, ndwi, ndbi, bare])
var url = multiband.getThumbURL({'min':-1, 'max': 1, 'dimensions': 512, 'format': 'jpeg'})

This however adds in a lot more color than masking and adding, which might be complicated for the model I am trying to build.

Here, is the link to the updated code.


texts: Measure the Length of a Line on a Map after a MOUSE_MOVE Event

texts:
One of the actions a user can perform using the ESRI Flex Viewer application is to draw a shape on a map (e.g. a line). Once the user is finished drawing the line, they will double-click the mouse which fires a DrawEvent.DRAW_END event.  When this is fired, the line is measured using a call to an API function geometryService.project([geometryObject], spatialReference);.

I want to be able to call this function after a MouseEvent.MOUSE_MOVE event in order to measure the line on the fly, instead of only after the DRAW_END event.  Unfortunately, the only two DrawEvents that get dispatched are the DrawEvent.DRAW_START and DrawEvent.DRAW_END.  These ESRI classes are compiled classes, so I cannot make any changes to them.

Is there a way to add an eventListener to a MOUSE_MOVE event that can grab the current geometry of the line that is being drawn and call my API function to measure on the fly?

Thanks in advance.


texts: Launch Attribute Table from Sublayer Menu

texts:
I've added

layerControlLayerInfos: {
             menu: [{
                  label: 'Query Layer...',
                  iconClass: 'fa fa-search fa-fw',
                  topic: 'queryLayer'
              }, {
                  label: 'Open Attribute Table',
                  topic: 'openTable',
                  iconClass: 'fa fa-table fa-fw'
              }]
            },

What are the additional steps need to launch a attribute table from a sublayer menu?

How do I subscribe other widgets to the topic defined above?


texts: Two cartodb visualization on the same html page

texts:
I need to display 2 maps separate by vertical line on the middle of the html, it means on diferents  and make some buttoms or aplicate the same zoom on both at the same time.

I managed to put one script code map on the first div but when I do the same to the other on the second div it invalidates the previous map and only shows the second map.

Im newbie on javascript. How can I solve this issue?

I know that with Iframe can do it but I need to make functionalities "outside" dinamically with 2 maps at the same time like zoom or filter's or a timelapse for the two at a time. Any ideas for this?


texts: How does Nearest Neighbor Image Resampling work in ArcGIS?

texts:
I am trying to understand how nearest neighbor for the resampling of Image datasets works in ArcGIS.

The output raster cell value is the value of the closest cell value in the input raster:

In this case the center of each output cell is the middle cell of each 3x3 input cells.

what happens if they are all at the same distance?
if the output has half the dimension of the input, the center of the output will have the same distance to the 4 closest neigh input cells ?

InRas1=6x6
OutRas=3x3

He gets the majority of cell values? No

Or am I missing something here?


texts: OpenLayers 3: Vector Layer loses styling when map pans or changes zoom

texts:
I'm working on switching from OpenLayers 2 to OpenLayers 3 and supplementing my reading of The Book of OpenLayers 3 by trying some things out on my own. There's a project that I'm building here that should be fairly simple: take form inputs and use them as parameters for building a vector layer. The interface works as intended, save for one problem.

Whenever I pan the map or zoom it, the vector layer seems to disappear. I can't for the life of me imagine what's going on. Anyone able to help me figure out what I can do to keep the layer in place?

To make things easier, I have created a JSFiddle, please feel free to use at will.


texts: PostGIS: Is it possible to identify overlapping polygons in the same layer using their geomteries?

texts:
I have a single table containing the geometries of ~3 million polygons and I am looking to calculate the number of polygons that overlap with another.

I know that you are able to do this in QGIS using the Topology Checker Plugin in Plugin Manager, I have used this when comparing much smaller numbers of polygons however I don't see how this would handle such a large amount of polygons. Also, from my research you are unable to use the results i.e. get a list of id's/display a map

After deciding the Topology Checker Plugin wouldnt work I found this:

http://postgis.17.x6.nabble.com/Method-to-remove-overlaps-in-a-layer-td4621793.html

I followed the steps and created the topology for all of the polygons successfully however my relation table wasn't populated and therefore the final step didn't work:

SELECT r.element_id as face, array_agg(r.topogeo_id) as geoms
FROM topo3.relation r
WHERE r.element_type = 3 -- 3 == face
GROUP BY face
HAVING count(r.topogeo_id) > 1;

I am running out of ideas and was wondering if anyone could possibly help


texts:
layer = iface.activeLayer()

# get unique values for 'severity' field
fni = layer.fieldNameIndex('severity')
unique_vals = layer.dataProvider().uniqueValues(fni)

# define categories to use in symbology
categories = []

for val in unique_vals:
    # initialise the default symbol for this geometry type
    symbol = QgsSymbolV2.defaultSymbol(layer.geometryType())

    # configure a symbol layer
    layer_style = {}
    layer_style['color'] = '%d, %d, %d' % (randrange(0,256),randrange(0,256),randrange(0,256))

    layer_style['outline'] = '#000000'
    symbol_layer = QgsSimpleFillSymbolLayerV2.create(layer_style)

    # replace default symbol layer with the configured one
    if symbol_layer is not None:
        symbol.changeSymbolLayer(0, symbol_layer)
    else:
        print "success"

    # create renderer object
    category = QgsRendererCategoryV2(val, symbol, str(val))
    # entry for the list of category items
    categories.append(category)

 # create renderer object
 renderer = QgsCategorizedSymbolRendererV2('severity', categories)

# assign the created renderer to the layer
if renderer is not None:
    layer.setRendererV2(renderer)

layer.triggerRepaint()

I have used this post, and copied the code above in to the python console. It appears to work fine, as I end up with the correct number of categories, i.e. 1 for every unique value in my chosen attribute, however after categories.append(category) it prints out False for every iteration of the loop. End result is fine though, so would anyone be able to explain why it's printing out False?

I thought it might have been if symbol_layer is not None: etc. but I added an else: print "success" and it's not doing that so doesn't look like that's the issue. Like I said, it does work but I'm curious

perhaps worth mentioning the layer I'm using is from a postgres database


texts: Determining distance to nearest location along stream network by each cell using GRASS?

texts:
I have a river network in raster and locations of many dams on the river.

In the image above, the BLUE point represent one random cell. Dam 1 is downstream to the cell.

I am looking to calculate the distance to the closest dam for each cell. And the distance should be calculate along the river network.

In the above picture, the distance from the blue point to Dam2 should be returned to the blue cell.

Does anyone know how to solve the problem using GRASS? I need the result as a raster file.

Some of my guess are given below.

I'm guessing some tools like r.cost?


texts: Saving a raster with writeRaster and getting an error

texts:
I am trying to save a Rumble Index derived from a LAS point cloud catalog. The code follows an example from and it does provide a graphic output with the plot() command. However, I wish to save it to my hard drive with the writeRaster command and I get an error. The code I am running is below.
source("D:/TAGB/Scripts/0000_setup.R")

rumple_index_surface = function(cluster, res)
{
  las = readLAS(cluster)
  if (is.empty(las)) return(NULL)
  las <- filter_surfacepoints(las, 1)
  rumple <- grid_metrics(las, ~rumple_index(X,Y,Z), res)
  bbox <- raster::extent(cluster)
  rumple <- raster::crop(rumple, bbox)
  return(rumple)
}

project <- readLAScatalog("D:/TAGB/2018/LAZ_N2/")

opt_chunk_buffer(project) <- 1
opt_select(project) <- "xyz"
opt <- list(raster_alignment = 25, automerge = TRUE)
output <- catalog_apply(project, rumple_index_surface, res = 25, .options = opt)

plot(output, col = height.colors(50))

rumple <- raster(output)

writeRaster(rumple, name = "D:/TAGB/2018/Metrics/rumpleIndex_25m.img", overwrite=TRUE)
#> Error in (function (classes, fdef, mtable)  :
#>  unable to find an inherited method for function ‘writeRaster’ for signature ‘"RasterLayer", "missing"


texts: Pixel-based image classification (Orthoimage or Raw image?

texts:
I am studying image classification and wanted to produce some land cover maps.
However, i can not seem to find resources that states whether you need to geometrically correct (orthorectification) your raw satellite  image first before you conduct your image classification or its ideal to use the raw satellite image data to retain the original pixels for a more effective image classification.
Or would it depend on the type of classification whether pixel-based or object-based classification?
Thank you for your help!


texts: Transforming from WGS84 UTM30N to British National Grid (BNG)?

texts:
I am using ArcGIS 10 and trying to work in British National Grid.

My shapefile is in WGS84 UTM30N.

I passed the shapefile across to a colleague and we both transformed into BNG.

When we compared the data there was a significant difference in our results.

I have since cross-checked with QGIS, Global mapper and a website or two and got a whole range of results but in general the most common is the result my colleague got.

What have I done wrong?
Why is there such a variation (~150m) in values? (1 or 2m error is
probably ok for what I am looking at.)
Which value is "right"?


texts: Transforming TM75(29903) squares to ETRS89(3035)?

texts:
A bit of a long winded explanation here so please bear with me. I have received grid square data in ETRS89(3035) which I need to cross reference with out own grid squares which are in TM75(29903). However I am having issues with the transformation. I have created some images in QGIS to demonstrate.

The following shows a 10km fishnet created in 3035 bounding Ireland. The pink squares are the 3035 data which I received.

However if I simply re-project my 29903 squares to 3035 you end up with the following issue. You can see that the green squares are tilted. The layer is re-projected but since the squares were drawn using the parallel lines on 29903 they appear tilted in 3035. I realised the issue was likely the same as explained here Projection Transformation from ETRS89 LAEA to WGS1984?. The re-projection is done but the squares need to be essentially redrawn in 3035.

I surmised that the best way to approach this would be to get the centroid of each of my 29903 squares, re-project that point to 3035 and then rebuild my 10km square buffer. I am using shapely to rebuild the buffer. This works..sort of. But as the following image shows the squares are essentially rotated individually and so are slightly out of position or overlapping.

At this point, the way forward which I am considering is iterating through my re-projected centroid points and finding the 3035 square from the 3035 fishnet which contains said centroid and then setting my updated geometry to that square. Does that seem like a valid approach?

The goal is to reproject my 29903 squares to 3035 and have them be on the same grid square scheme as the 3035 data which I need to cross reference. The same applies vice versa in fact. Included below is the script which I am testing to do the reprojection. Note I have not implemented my fishnet approach yet.

import ogr, osr
import os
from shapely.geometry import Point

driver = ogr.GetDriverByName('ESRI Shapefile')

# input SpatialReference
inSpatialRef = osr.SpatialReference()
inSpatialRef.ImportFromEPSG(29903)

# output SpatialReference
outSpatialRef = osr.SpatialReference()
outSpatialRef.ImportFromEPSG(3035)

# create the CoordinateTransformation
coordTrans = osr.CoordinateTransformation(inSpatialRef, outSpatialRef)

# get the input layer
inDataSet = driver.Open(r'I:\Code\Python\input\test_shape\test_shape.shp')
inLayer = inDataSet.GetLayer()

# create the output layer
outputShapefile = r'I:\Code\Python\output\test_shape\test_shape.shp'
if os.path.exists(outputShapefile):
    driver.DeleteDataSource(outputShapefile)
outDataSet = driver.CreateDataSource(outputShapefile)
outLayer = outDataSet.CreateLayer("basemap", outSpatialRef, geom_type=ogr.wkbMultiPolygon)

# add fields
inLayerDefn = inLayer.GetLayerDefn()
for i in range(0, inLayerDefn.GetFieldCount()):
    fieldDefn = inLayerDefn.GetFieldDefn(i)
    outLayer.CreateField(fieldDefn)

# get the output layer's feature definition
outLayerDefn = outLayer.GetLayerDefn()

# loop through the input features
inFeature = inLayer.GetNextFeature()
while inFeature:
    # get the input geometry
    geom = inFeature.GetGeometryRef()
    print(geom.Centroid().GetX())
    centroid = geom.Centroid()
    centroid.Transform(coordTrans)
    geom = Point(centroid.GetX(), centroid.GetY())
    geom = geom.buffer(5000, cap_style=3)
    geom = ogr.CreateGeometryFromWkt(geom.wkt)
    geom.AssignSpatialReference(outSpatialRef)
    print(geom.Centroid().GetX())

    # create a new feature
    outFeature = ogr.Feature(outLayerDefn)
    # set the geometry and attribute
    outFeature.SetGeometry(geom)
    for i in range(0, outLayerDefn.GetFieldCount()):
        outFeature.SetField(outLayerDefn.GetFieldDefn(i).GetNameRef(), inFeature.GetField(i))
    # add the feature to the shapefile
    outLayer.CreateFeature(outFeature)
    # dereference the features and get the next input feature
    outFeature = None
    inFeature = inLayer.GetNextFeature()

# Save and close the shapefiles
inDataSet = None
outDataSet = None


texts:
When I select a object in a particular layer using qgis in order to get information about it, this object don't get resalted in red and in the qgis log a get the following error:

not found space object -1456845214

This issue occurs in some objects, but not all.


texts: Sunlight based on GPS coordinates and UTC+date

texts:
I'm trying to create a sunlight map. I found the following question: Sunlight layer into CartoDB

However, my interest is in trying to calculate the amount of sunlight manually. Can you please point me to the formulas/documentation on how I could compute the amount of sunlight (sunrise/sunset) based on GPS coordinates and UTC/date?


texts:
I do not have a great experience in ArcGIS, but recently I received 3 datasets in ArcGIS Geodatabase format, which I would like to work on. From these 3 datasets I can open only 1, but not the other 2, because after the successful layer selection ArcGIS simply crashes without any error message.

I am using ArcGIS Desktop 10.3. I also tried to open them with the most recent version of the free ArcGIS Explorer on another computer, with exactly the same result. I contacted the office from whom I received the data and we verified in their office that they can open the data on their computers from the very same external hard disk I was using before, therefore the data is valid. They were using ArcGIS 10.1 and 10.3 on their computers.

I am totally out of idea what is causing this issue, and it is also not helping that there are no error messages at all. If you have any advise please share.

UPDATE:
For request, the progress of facing the issue through screenshots is shown below:

Step 1: adding data

Step 2: selecting layer

Step 3: crashing

Update #2: crash info by Windows

Problem signature:
  Problem Event Name:   BEX
  Application Name: E3.exe
  Application Version:  2.0.0.3400
  Application Timestamp:    55b913ce
  Fault Module Name:    FileGDB.dll
  Fault Module Version: 2.0.0.3400
  Fault Module Timestamp:   55b948ff
  Exception Offset: 000c82f9
  Exception Code:   c0000409
  Exception Data:   00000000
  OS Version:   6.1.7601.2.1.0.256.48
  Locale ID:    1038
  Additional Information 1: 9f85
  Additional Information 2: 9f8512db2514abc8411174b85319aee2
  Additional Information 3: 9270
  Additional Information 4: 92708c4c8f00576bcf7d8ffce69d82f8


texts: identify percentage of raster cells within a polygon

texts:
As per the title I want to determine the population inside multiple polygons according to different land use zoning (High Residential, Low Residential, National Parks etc) that's in raster format.

To do this I'm thinking I should calculate the percentage of raster cells inside the polygon and then join the tables so I can calculate the population.

So far I've reclassified these land use zones to make 3 different binary rasters.
e.g. For High Res:- 1 for high Res and 0 for the other two (National Parks and Low Res); For Low Res:- 1 for Low Res and 0 for High Res and National Parks.

The next step would be to use "Zonal Statistics by Table" and join the table to the polygon table that has population information. However it seems I'm just adding up polygon population essentially coming out with the same number.


texts: Reduced resolution when exporing raster with layout view using ArcGIS

texts:
I have an ArcGIS project with raster and vector layers. I use the function Export map ... to export layout view to a tiff file with different DPI values (150 and 300). The resulting image is very different from what is displayed in the ArcGIS. The resolution of the raster layer is highly distorted regardless of the DPI option.

How to fix this issue?

Pic 1. ArcGIS layout view (natural raster layer resolution).

Pic 2. Export result - DPI 150 (raster layer resolution reduced).

Pic 3. Export result - DPI 300 (raster layer resolution reduced same way as with DPI 150).


texts: NBI Coordinate Data Missing place-holding zeroes for DMS format

texts:
When working with National Bridge Inventory (NBI) coordinates I have found what appears to be an error in the raw data provided.  I couldn't find any other mention of this specific error when searching for a solution to my problem.  There seems to be dropped zeroes in the seconds portion of the DMS format used.  According to the NBI Coding Guide the format of their coordinates is DDMMSSss.  After I had converted this data to decimal degrees, plotted the points and compared it to known bridge locations I realized that several of these coordinates were incorrect.  After trying a few different possibilities I discovered that there were missing zeroes in their data.  For example:

Raw Data
Lat 34212979 Long 84565142
Lat 34172838 Long 84584821
Lat 34171338 Long 84452538
Lat 34069870 Long 84461200
Corrected Data
Lat 34212979  Long 845605142
Lat 341702838 Long 84584821
Lat 341701338 Long 844502538
Lat 340609870 Long 84460001200

Where the coordinate was 84 degrees, 56 minutes, 5.142 seconds the NBI data would record this without the place-holding zero in front of the seconds so that, according to their formatting, it would read 51.42 degrees instead.  Even when the seconds are 0.012 the raw data would simply omit these place-holding zeroes and record 1200 in the seconds place.  Out of 216 records there were 64 that were missing these place-holding zeroes.  Had I not had pre-existing bridge locations to compare these to then 1/4 of my geocoding would have been incorrect.

My question is, How can I identify this incorrect data more easily in the future so that I don't have to go point by point to correct the formatting error?  So far, the only solution I found was to double check the coordinates of all the points that I know to be incorrect.  Also, has anyone else noticed such a problem with their raw data from the NBI? I didn't see any explanation for this in the NBI Coding Guide.


texts: Spatial Analysis Grouping/Clustering Question

texts:
If this question is incorrectly titled, I apologize; I wasn't sure exactly how to best phrase it. I'm working with my city's housing department to develop zones for building inspectors. One of the attributes of each building footprint/polygon is total building square footage. The department is requesting that I find a way to group all the buildings into zones that have roughly equivalent total square footage. I briefly looked into Grouping Analysis and Location-Allocation, but I'm not sure either of those is what I am looking for. IIRC from graduate school, ESRI does have a tool for this, but it's been almost 10 years since graduate school.

Is there such a tool, or am I remembering incorrectly?

If not, could anyone suggest a good way to do this?

Spatial statistics and analysis are not really my strong suit.


texts:
I am trying to run Solar analysis in ArcGIS Desktop and it is an enormous DEM of a very mountainous region.     I am going to have to split the DEM.

Are there any issues with doing this?

For example, if you analyze a section separately from an adjacent mountain that would have shaded this section.


texts: Using Attribute Values of Rasters in ArcGIS Desktop?

texts:
I have 4 Shapefiles which basicly have 7 important fields and are basicly looking like:

------------------------------------------------------------------
| Name | Connection 1 | Con 2  | Con 3  | Con 4 | Con 5  | Value |
------------------------------------------------------------------
|Tab1_A|    Tab2_A    | Tab2_B | Tab3_F | -     | Tab4_x |  10   |
------------------------------------------------------------------

So the first collumn contains the name of the polygon, the 2nd to 6th collumn contain Connections this polygon has with polygons from the other Shape-Files and the last collumn contains a value.

My Task is, to create a Raster, which contains the maximum sum of all values from conntected shapes in one cell.

My current approach is:

Convert the 4 Shapes into 7 Rasters each (one per collumn)
4 Raster Calculator: Tab1: Value + If(Con1 == Tab2.Name;Tab2.Value;0)+...
Another Raster Calc: MAX(Raster1,...,Raster4)

My Problem is, that I cannot compare the Name-Raster and the Connection-Rasters, because the Name is not the Raster-Value, but an Attribute.

Is there any way to compare an attribute in the Raster Calculator, or do you have a better idea, how to solve this problem?


texts: How do v.distance 'upload' parameters actually work?

texts:
I'm trying to use v.distance (using QGIS 3.4) to calculate the distance in meters between a layer of water shutoff points and the boundary of the nearest building footprint polygon. Even better would be the distance to building footprint with a matching "ASSOC_SC" value (associated service connection).

I have tried many, many different combinations to try and satisfy the 'upload' section of the function parameters and I cannot for the life of me understand how exactly it works.

Do I need to create an empty column of the right data type in my 'from' layer for the chosen measurement result to go and then choose that in the "Column name of the nearest..." section? (e.g. to match cat for identifying nearest, distance for listing the actual measurement, etc.)  When I do that, it does not  work.

If I'm trying to measure based on a matching ID attribute, why does it ask me to select "to_attr)" in the "'upload: Values..." section and then not list any of the columns in the "Column name of the nearest..." section?

What does that "Column name of the nearest..." language actually mean?

It allows me to choose an ID column from my 'to' layer in the dropdown below but nowhere allows me to select from my 'from' layer what to match...

Has anyone made use of this function from the QGIS processing toolbox like this?

All layers are in the same CRS, etc.


texts: What happens in the DB when a user updates a record in a versioned table?

texts:
ArcGIS Desktop 10.7.1/Oracle 18c.
In a table that is versioned (with the option to move edits to base):
What happens in the database when a user updates a record?
It would be helpful to know what exact steps are performed in the database, and in what order.


texts: How does Mobile GPS work in Windows Tablets?

texts:
Can anyone tell me how GPS chips generally works in Windows8 Tablets? As far as I can tell Tablets do not have a separate GPS chip like those in handheld GPS receivers - It seems to be combined into the Network chips. With my current tablet, gps doesn´t start without first logging into the network with a mobile SIM card. They seem to first get a fix according to the location of the last successful IP Address login and then refine the position using triangulation from satelites.

I would be grateful if anyone with experience with such devices could clear this up for me.


texts: Efficiently selecting related records using ArcPy?

texts:
Below is the code I'm using to replicate the "related tables" button in ArcMap.  In ArcMap that button selects features in one feature class or table based on the selection of features in another related feature class or table.

In ArcMap I can use that button to "push" my selection to the related table in a matter of seconds.  I was unable to find anything built in to arcpy that replicates the button so I used some nested loops to do the same task.

The code below loops through a table of "treatments".  For each treatment, it loops through a list of "trees".  When a match is found between the ID fields of treatment and trees, a selection occurs in the tree layer.  Once a match is found for a treatment, the code does not continue to search the tree layer for additional matches.  It goes back to the treatment table, selects the next treatment and again searches through the tree feature class.

The code itself works fine, but it is agonizingly slow.  The "treatment table" in this case has 16,000 records.  The "tree" feature class has 60,000 records.

Is there another more efficient way to recreate what ESRI is doing when it pushes the selection from one table to another? Should I be creating an index for the tables?  NOTE: This data is stored in an SDE.

 # Create search cursor to loop through the treatments
treatments = arcpy.SearchCursor(treatment_tv)
treatment_field = "Facility_ID"

for treatment in treatments:

    #Get ID of treatment
    treatment_ID = treatment.getValue(treatment_field)

    # Create search cursor for looping through the trees
    trees = arcpy.SearchCursor(tree_fl)
    tree_field = "FACILITYID"

    for tree in trees:

        # Get FID of tree
        tree_FID = tree.getValue(tree_field)

        if tree_FID == treatment_FID:
            query = "FACILITYID = " + str(tree_FID)
            arcpy.SelectLayerByAttribute_management(tree_fl, "REMOVE_FROM_SELECTION", query)
            break


texts:
I am learning how to use Arc Hydro and attempting to change the Z Units because I get a warning message saying the Spatial Reference Doesn't Match.

The DEM was in Meters, and I adjusted the Z Unit to 100.

My DEM is in a GCS and when I project it to the same PCS my vector data is in, it puts them in two very different places, even though they are in NC State Plane 3200 FIPS (feet).

How do I fix this?


texts:
I want to contribute to GeoServer development. I am well experienced in programming and since I am now completely working on GIS projects, I think I will find enough time for these tasks.

Where should I start? Are there to do lists or tasks for developers? Are there basic jobs for beginners?

What would you suggest?


texts: Upload of geojson to Mapbox studio showing strange artifacts

texts:
I want to upload Geojson data to Mapbox studio. The upload works correctly, but the data in Mapbox show artifacts, which do not exist in the original data.
It seems that some vertex has a wrong position.
Please see a screenshot of the original data (in Qgis):

And this is the same section of the data in Mapbox Studio:

I tried using the Mapbox UI as well as the Mapbox cli: both datasets have artifacts, but different ones.
Original data: https://we.tl/t-vWHnewZ2TI


texts: Processing Modeler QGIS Issues With Interpolation Step

texts: I am trying to interpolate Sampled Points with IDW Interpolation in the QGIS processing modeler.  For some reason I am not able to select the Sampled Points with the IDW Interpolation.  Anyone know why that would happen?


texts: Filter Sentinel2 Image Collection With IRRADIANCE Metadata

texts:
I am trying to filter an Image collection in google earth engine, but it seems that my metadata filter based on irradiance has no discernable effect on the output.  Values above the threshold are still kept in the ImageCollection, but if I reduce the value too low, no images make it in (even though some should)   What am I doing wrong?  The other metadata filter seems to work (Cloudy pixel percentage).

var sentinel2 = ee.ImageCollection('COPERNICUS/S2')
.filterMetadata("SOLAR_IRRADIANCE_B4","less_than", 2000)
.filterDate('2015-06-22','2018-04-09')
.filterMetadata('CLOUDY_PIXEL_PERCENTAGE','less_than',20)
.filterBounds(Caselton)
.select(['B4','B3','B2'])
.sort('system:time_start',true);

// Create an image time series chart.
var chart = ui.Chart.image.series({
  imageCollection: sentinel2,
  region: Caselton,
  reducer: ee.Reducer.mean(),
  scale: 200
});

// Add the chart to the map.
chart.style().set({
  position: 'bottom-right',
  width: '500px',
  height: '300px'
});
Map.add(chart);

// Outline and center Either Pond
var caselton = ui.Map.Layer(Caselton, {color: 'FF0000'}, 'CMT');
Map.layers().add(caselton);
Map.setCenter(-114.49052095413208,37.903421736152076, 15);

// Create a label on the map.
var label = ui.Label('Click a point on the chart to show the image for that
date.');
Map.add(label);

// When the chart is clicked, update the map and label.
chart.onClick(function(xValue, yValue, seriesName) {
  if (!xValue) return;  // Selection was cleared.

  // Show the image for the clicked date.
  var equalDate = ee.Filter.equals('system:time_start', xValue);
  var image = ee.Image(sentinel2.filter(equalDate).first());
  var s2Layer = ui.Map.Layer(image, {
    gamma: 1.3,
    min: 0,
    max: 2500,
    bands: ['B4','B3','B2']
  });
  Map.layers().reset([s2Layer, caselton]);

  // Show a label with the date on the map.
  label.setValue((new Date(xValue)).toUTCString());
});


texts:
I'm using Weighted Overlay from Spatial Analyst Tools in ArcMap to calculate the weights of two raster layers while calculating, the process just works in an area, not the whole area I'm trying to do.

The weights are just calculated in the area of RN Layer, not in the LandCover area.


texts: From everywhere to everywhere in ArcGIS Network Analyst?

texts:
ArcGIS Network Analyst question:

I have a network dataset representing an archaeological site with ~500 buildings. Each doorway of these buildings represents a terminal node. I would like to simulate aggregate foot traffic patterns at the site.

One way to do this would be to generate least cost network paths from every node to every node. That is, generate 500 x 500 route matrix (250,000 routes). That is the same as doing a closest facility analysis for every building at the site. Of course, manually doing a closest facility analysis 500 times would be ridiculously time consuming.

Can I do it using ModelBuilder?


texts: Change SVG symbol size based on a field value range

texts:
I have been experimenting with QGIS for about three weeks now.

I am using a SVG symbol that I created and successfully loaded into the directory for QGIS 3.10.  I can make it appear on the map.

It represents the size of a colony and I would like the size to reflect the number of inhabitants, based on ranges of numbers.

I have created a layer to represent these colonies, with a field called "number"

What I want to do is have the size of the SVG symbol for each colony change based on the value in the field "number"

Ideally, I want to use ranges of numbers, so a colony between 1-10 inhabitants would be the smallest size symbol, 11-50 the next size and so on.

I have tried using the data defined override box, but clearly have no idea of correct syntax to either identify that I want the size to change based on a number range taken from the "Number" field or to specify that it is the symbol size I want to change.

If anyone can point me to an appropriate tutorial, I have been unsuccessful finding one on my own!

If you have a process that I am missing that would be great

If there is a way to do data defined, I would be interested in that too.

My programming background is an ancient use of OpenScript (Proprietary for the old Toolbook program) and a minuscule amount of whatever you use to make excel do tricks!


texts: "Import Textfile as Layer" generates only 1 feature

texts:
I have a text file containing adresses with coordinates wich looks like this:

[...]
N;DEBYv00018000101;A;09;4;75;120;0001;00000;9;;4496255,330;5571332,790;Brunnenstraße;95182;Döhlau;, Kr Hof, Saale;Döhlau;;
N;DEBYv00018000106;A;09;4;75;120;0001;00000;10;;4496153,010;5571475,360;Dahlienweg;95182;Döhlau;, Kr Hof, Saale;Döhlau;;
N;DEBYv00018000111;A;09;4;75;120;0001;00000;11;;4496086,420;5571470,390;Dahlienweg;95182;Döhlau;, Kr Hof, Saale;Döhlau;;
N;DEBYv00018000117;A;09;4;75;120;0001;00000;13;;4496066,170;5571481,620;Dahlienweg;95182;Döhlau;, Kr Hof, Saale;Döhlau;;
[...]

No field names are given in a first line and fields are separated by ;. The file has 1504 lines and what I get when opening this file is not the expected 1504 features with 19 attributes each but only one feature with 28575 attributes as shown in the screenshot:

The settings seem correct to me, so what am I probably missing?


texts:
I am trying to teach my self Django and its spatial extension geoDjango.
There is one thing, which is still a bit foggy in my mind.
When using the trio: Postgis, geoserver, openlayers is perfectly clear for me how I will display the data on the map.
But if geodjango doesn't use geoserver, how do I display the data on the map using Openlayers? Should I export the data as geojson or kml and then use openlayers functions to render them? That would be very limiting.

What's the proper way to do this?


texts: Efficient route planning between a set of points and with a set of people

texts:
I have the following use case:

A company needs to make winter service at a set of addresses in a city. To be as fast as possible covering all addresses an efficient route calculation is required. However, it shouldn't be one single route between all addresses but a set of routes depending on the number of workers. So if the company has 4 workers the algorithm should find the most efficient SET of 4 routes between all addresses in which no address is covered twice. The start and end of each route is same (being the building of the company) and each route should cover a similar amount of addresses

Is there an algorithm for that particular use case? I am using QGIS 2.18.


texts: Adding hyperlink in a popup which comes from a polygon?

texts:
I need to add a hyperlink in a popup which comes from a polygon. The following piece of script shows the text which needs to be a link.

var polygon = L.polygon([
[51.569, -4.144],
[51.568, -4.145],
[51.57, -4.146],
[51.57, -4.145]
])

.addTo(map);

polygon.bindPopup(" RTK Survey Area ");


texts: Extrude Shapefile to KML based on field value

texts:
I would like to extrude shapefile polygons to a KML based on a value stored in the attribute table.  I am hoping to do this using open source software (QGIS, FWTools, GDAL, etc.).  I looked at ogr2ogr using the -zfield option, but the polygons were not extruded as a KML, they were just an outline.

Any help would be appreciated,

Thanks...


texts: Using arcpy.da.UpdateCursor with dictionary?

texts:
I'm working through some da.UpdateCursor scripts to make updates to a fgdb. Although the script below works successfully, I want to use the dictionary that I have created rather than the manual [ ] system I have for demo purposes.

Can anyone provide some guidance please as I have researched this quite significantly and can't get an answer after some hours of effort and no error message. (Using Win 7, ArcGIS 10.2, Python 2.7.5)

The code below updates 'book_ref' where it finds the matching 'book_id' in 'seq = .......'

    seq = [[7643, 7625], [9644, 2289]]
    Att_Dict = {}
    for pair in seq:
        Att_Dict[pair[0]] = pair[1]

    with arcpy.da.UpdateCursor(in_fc1, ['book_id', 'book_ref']) as cursor:
        for row in cursor:
            pozzy = row[0]
            if pozzy in Att_Dict:
                row[1] = Att_Dict[pozzy]
                cursor.updateRow(row)

How do I get it to use the values from my dictionary?

seq = {7643: 7625,
       9644: 2289,
       4406: 4443,
       7588: 9681,
       2252: 7947}


texts: How does CartoDB do international geocoding?

texts: I work on projects overseas so I wanted to see how it does international geocoding -- can or does it use addresses, town names, neighborhood names, or what? And what areas does that cover?


texts:
We set up a QGIS server and are using QGIS to connect to a WFS layer. Everything is working fine within the desktop application. Then we use qgis2leaf to generate a html/js web map. That all works fine, except when I use the WFS layer. For some reason I am getting the Unexpected token < error, as can be seen here (http://qgis.webmuni.com/export).

So basically what I think is happening is that the call that is made is (http://qgis.webmuni.com/qgis_mapserv.fcgi?SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=zones2&SRSNAME=EPSG:4326&outputFormat=text/javascript&format_options=callback:processJSON) and it expects it to return JSON, but the server is returning XML.

I don't know why the server wont return JSON(P) and I am not sure if there is something I am doing wrong or if we need to configure the server in some way to enable JSON capabilities.


texts: Postgis/SpatiaLite SQL query to connect MySQL, Access or Oracle database

texts:
In SpatiaLite or Postgis, how can one create a SQL query to read a table from a  MySql or Oracle database?

I don't know **the correct SQL syntax to connect to a table in a MySql or Oracle database **.


texts: Query Geoserver using lat/long with UTM store

texts: I have been using a Geoserver for years with my developers querying it using lat/long to get US parcel boundaries from a SQL Server store. I got a dataset with UK parcel boundaries but it's in EPSG:27700 (UTM). I can query it fine using UTM values but I need to use lat/long and I'm not sure what the trick is. The Geoserver lists the native SRS as EPSG:4326 for some odd reason. I think it should say EPSG:27700 for native and I want EPSG:4326 for declared, with either "reproject native to declared" or "force declared." Any guesses as to what I'm doing wrong?


texts: Filling attribute tables with values from the same attribute table

texts:
Just wondering  if there's an easy way to fill null values in one column in an attribute table with values in another column of the same attribute table by just writing a query "just the "null values" to be filled with details I already have in another column".Found this query
 field 1= sec volt, field 2: Lv connection

def updateValue( !SecVolt!, !LVConnection! ):
if !SecVolt! is None: return !LVConnection!
else: return '0'

but in my case I don't want to replace the value with "0" but with values in another column.


texts: Coordinate Point Drops on timeseries in QGIS

texts: I have timestamp data for every minute of a day (1140)for multiple objects by id, with a corresponding lat long. I want to get the time in which an id disappears off the map or off the csv. What would be the best tool, and is QGIS the best route?


texts: how to decrease Raster size without altering cell size?

texts:
Is there a way to decrease the size of a raster image without decreasing cell size or quality?

I´m using ArcInfo but also have Saga-Gis or FME.

I could abviously change an 8bit raster to a 1bit raster and thereby lose the colour information, but I was thinking that there is maybe other "secret" or more lesser known methods.