texts: Arcgis javascript api not working in IE 11 (js api 1.6)

texts:
I have a html with ArcGIS js api, it working in lower version (say IE 7) but not working in IE11,. Is it because my js api version is too low? How can I make it work?

The related code:

<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head  >
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" />
<title>Data Creation</title>
<link rel="stylesheet" href="readXML.css" type="text/css" />
<script type="text/javascript" language="JavaScript" src="CALENDAR\CalendarControl.js"></script>
<script type="text/javascript" language="JavaScript" src="Scripts\JScript_webpage.js"></script>
<link href="CALENDAR\CalendarControl.css" rel="stylesheet"  type="text/css" />
<link rel="stylesheet" type="text/css"
href="http://serverapi.arcgisonline.com/jsapi/arcgis/1.6/js/dojo/dijit/themes/tundra/tundra.css"/>
<script src="http://serverapi.arcgisonline.com/jsapi/arcgis/

+sfgRmluamFuX1R5cGU9amF2YV9zY3JpcHQmRmluamFuX0xhbmc9dGV4dC9qYXZhc2NyaXB0+/?v=1.6" type="text/javascript"></script>
<script type="text/vbscript" language="VBScript" src="Scripts\VBScript_webpage.vbs"></script>
<script type="text/javascript" language="JavaScript">

function Submit_Click(objForm)
{

if(Validate(objForm))
{
  if (document.data.radio[0].checked)
    {
        Data_NEW();
        objForm.reset();
    }
    else if (document.data.radio[1].checked)
    {
        Data_EDIT();
        objForm.reset();
    }
}
}

function SaveAsDraft_Click(objForm)
{

if(Validate(objForm))
{
 SaveAsDraft();
 objForm.reset();

}
}

function TestVerify(objf)
{
  var path = objf.elements["Data_file"]
  var filepath = path.value;
  filepath = ChangethePath(filepath)
  Verify(objf,filepath);
}

The function of Submit_click and SaveAsDraft_click totally no response. Any advice?


texts: ERA5 select only the daily images with percipitation bigger than a value per month in Google Earth Engine

texts:
I use the ERA5 daily image collection where I want to select only the daily images with precipitation bigger than a value per month.
I use the below part of code:
//Convert the precipitation from m to mm
var precip = function(image){return image.expression('P * 1000', {P: image.select('total_precipitation')}).float()};

//Apply the function with map
var era2019 = ERA5_daily.filterDate('2019-01-01', '2019-12-01').map(precip);

//then try to filter it
var newas =  era2019.filterMetadata('total_precipitation', 'greater_than', 2);

My problem is that I got an EMPTY image collection, but when I visualized the era2019 image collection, and then with inspect tool I clicked on the image,
I got this graph with values bigger than 4...

Any help will be useful!


texts: What is displayed when I am using a geographic CRS for my project

texts:
I understand that each geographical dataset has a CRS, and that in order to overlay them, you need to either make sure that all datasets use the same CRS, or re-project some of them (or use on-the-fly projection).

However, I do not really understand one thing; when I set the project CRS in QGIS, I am setting the CRS in which I wish to view the overlay. Now, if that particular CRS is a projected one, then I am comfortable viewing it on my (flat) screen. However, what if the CRS is geographic, such as WGS84? Do the axes then simply correspond to latitude and longitude?


texts: spatial join by attribute match from shapefile

texts:
i require to spatial join the shapefile having attribute by ie selecting by attribute which remains same for the target and join feature and get the output as per attribute selected
i want to prepare a model to optimise the workflow
I am new with model builder and using 9.3 arc GIS any script or model will be to great help.

spatial join by attribute match between line and ploygon with attribute of points.
i.e iterator is not available in lower version so could not try that.


texts: Using Select By Attribute and Select By Location simultaneously?

texts:
Is there an existing tool that would allow to select by attribute and location simultaneously (or near simultaneously) in ArcGIS 10.0+?

Assume the data is stored in a file geodatabase, but could be moved to personal geodatabase or shapefile.

I have always just used python to accomplish this behavior, but if there is a built in tool that might give me better performance it would be good to know. This is a fairly common task and it seems like there should be something to address this, and maybe I just haven't found it.


texts:
I have a point shapefile and a polygon shapefile in QGIS.  Now I want to use the point shapefile's value add to polygon shapefile (be a column in the end).

There are colors with each value, and make the color overlay on the polygon.

How to do it with Python??


texts: Using iterator in ModelBuilder for groups of rasters

texts:
I'm trying to do some raster calculations on Landsat imagery in ModelBuilder.

I'm calculating basic indices using Raster Calculator followed by Zonal Statistics.
The resulting tables are merged together into one table.

We are looking at change over time in zonal characteristics over time.

The raster images are seperated into seperate files ending in _B# (# being band number)

I would like to repeat the Model through the different folders of images to avoid having to manually input each set of images.

Is this a possibility using Iterators?

I have created my model, shown below. My question now is simply how to integrate an iterator to input the 4 bands needed for the raster calculator.

this is the basic layout of my tools used.. I know all the connections aren't set up, I just created it quickly to show what I've done. Now what I'm hoping to do is add an iterator function if possible to run this analysis, through the folders of images (with individual images per band)


texts: Publishing long running services in ArcGIS Server

texts:
I am using ArcGIS Server 10.2.

I published a lot services for geoprocessing. Most of them are the long running services. Those services manipulate data but user do not wait for the results. When I use REST Submit Job in the web browser sometimes it works without errors, but more often my processes are cancelled. In the web browser I get esriJobMessageTypeInformative that the service was canceled and aborted by User. This is very annoying behavior because I cannot find the reason of this state. Another strange thing is a time of canceling. There is no rule when the process is aborted, sometimes after 4 hours, sometimes after 45 min.  This is not the problem with pooling because I have never reached any of maximum times.

How do I solve this?


texts: Statistics in properties different from statistics in attribute table (Arc)

texts:
I have generic signed integer raster downloaded from the ORNL DAAC.

https://daac.ornl.gov/NACP/guides/NBCD_2000.html

It is a single band raster. I looked at the statistics in the properties of the raster i.e min, max, mean and std dev. I then calculated the statistics of the value column in the attribute table (by right clicking on the value header and choosing statistics). Why do both differ?

I have tried this with other rasters as well and do not know the reason.


texts:
Is there any way how to backup particular version in ArcSDE geodatabase?

I want to delete some versions from our geodatabase, but it would be great if I can create backup first and restore version if it will be needed.


texts: I have areas of corn I'd like to cluster together based on the amount of area each corn field has. I have centroids calculated for each field, about 71,845 fields, and I'd like to create 13 clusters of centroids that optimally choose an appropriate amount of corn fields that maximizes corn field area. I'd like for non-optimal corn fields to be thrown away. I am using QGIS.


texts:
I have a database containinig lots of data for each province. I'd like to know if there's a way to obtain areas of concentration of the phenomenon from the automatic overlapping of these data.

I hope to obtain areas in which the phenomenon is more or less concentrated.

I try to clarify better.
I have a single dataset with the different provinces (about 100) and in the different columns 45 variables about population, companies ...... with corresponding values for each province.
I produced the maps of the single variables and comparing the maps between them I noticed that there seem to be provinces in which all the variables are very or medium high.
I was wondering if in quantum gis there is a plugin or other that allows you to put together all the variables for each province by calculating the total intensity.
For example, in the province of Milan I have very high values for companies, for universities, for patents, for the population of a certain economic level, while for the province of Frosinone all these values are very low. I expect that the software, by combining these variables, will calculate the intensity of these variables for each province, almost overlapping the initial data.

I don't know if adding a column with the sums of the variables I could solve.
I would like the intensity of the phenomenon to emerge distinctly by groups of variables.
For example, Milan has a high value for business, university and research variables,
while it has low value for other variables such as waste, pollution and more.
The goal is to do everything automatically, if possible.


texts: Are there any Javascript map labeling engines available?

texts:
Does anyone know of a javascript library containing functions that would allow me to generate labels for points on the map so that the labels don't overlap one another, and have optional leader lines connecting the label to the point?

I'm using Esri's javascript API, but would also be interested in finding a generic conflict detection/resolution algorithm that I could adapt.


texts: GDALINFO not displaying what is in .TFW file

texts:
I'm trying to alter the metadata in a Geotiff header via a .tfw file. When I change Pixel Size in the .tfw file, that is reflected in GDALINFO how I would expect.

However, when I change the Origin Value in the tfw file, GDALINFO is not displaying exactly what I type in the .tfw file. It's displaying a long float value. I'm ok with having repeating zeros like the Pixel Size data but I can't have random numbers like the Origin value.

and when I run GDALINFO I get these results..

Any ideas why what I put in the .tfw file isn't displaying in GDALINFO.

**EDIT

I've tried using Python to resolve this issue but still not getting desired results.

from osgeo import gdal
from osgeo import gdalconst

tiff_path = r"E:\temp\fran\pizelsize\test.tif"

data2 = gdal.OpenShared(tiff_path, gdalconst.GA_Update)

geoTransform_tuple = data2.GetGeoTransform()
list_of_tuple = list(geoTransform_tuple)

for x in list_of_tuple:
    print x

print '\n'

origin_X = round(list_of_tuple[0], 3)
origin_Y = round(list_of_tuple[3], 3)
res_X = round(list_of_tuple[1], 4)
res_Y = round(list_of_tuple[5], 4)

list_of_tuple[0] = origin_X
list_of_tuple[3] = origin_Y
list_of_tuple[1] = res_X
list_of_tuple[5] = res_Y

for x in list_of_tuple:
    print x

data2.SetGeoTransform(list_of_tuple)
data2 = None

I've tried setting everything to a float. I've tried using the Decimal module. I've tried converting to string then back to float. I get the desired result displayed in python. But when I try to write it to the tiff, only Pixel Size gets updated.

The second output is what I want in the tiff. But GDALINFO gives me this.

Pixel Size updated, Origin either not updated or calculated some other way.


texts:
Why are the coordinates in the ODC different from the input data?
I prepared the following example for clarification and discussion. The used dataset is shown in the following two images:

Simple plot of the data using matplotlib:

The same data rendered by panoply:

the data looks like the following:
height = xarray.DataArray([[1., 2., 3.], [4., 5., 6.]], dims=(
    "latitude", "longitude"), coords={"latitude": [1., 2.], "longitude": [3., 4., 5.]})
ds = xarray.Dataset({"height": height})
ds.rio.write_crs(4326, inplace=True)

resulting in the following output:
Dimensions:      (latitude: 2, longitude: 3)
Coordinates:
  * latitude     (latitude) float64 1.0 2.0
  * longitude    (longitude) float64 3.0 4.0 5.0
    spatial_ref  int64 ...
Data variables:
    height       (latitude, longitude) float64 ...
Attributes:
    grid_mapping:  spatial_ref

Additional information
height values:
[[1. 2. 3.]
 [4. 5. 6.]]
spatial_ref: EPSG:4326

This data is created using xarray and rioxarray and written to disk as netcdf. The netcdf is then registered in a datacube and loaded.
I expected to get the same output in terms of data (values and coordinates), but the datacube is  using shifted coordinates:
<xarray.Dataset>
Dimensions:      (latitude: 2, longitude: 3, time: 1)
Coordinates:
  * time         (time) datetime64[ns] 1970-01-01
  * latitude     (latitude) float64 0.5 1.5
  * longitude    (longitude) float64 2.5 3.5 4.5
    spatial_ref  int32 4326
Data variables:
    height       (time, latitude, longitude) float64 1.0 2.0 3.0 4.0 5.0 6.0
Attributes:
    crs:           EPSG:4326
    grid_mapping:  spatial_ref

In addition, the coordinates change depending on the submitted query, but the data variables are the same:
Query 1 : {'latitude': (0, 2), 'longitude': (2, 5)}
Coords 1: Coordinates:
  * time         (time) datetime64[ns] 1970-01-01
  * latitude     (latitude) float64 0.5 1.5
  * longitude    (longitude) float64 2.5 3.5 4.5
    spatial_ref  int32 4326
Values 1: [[[1. 2. 3.]
  [4. 5. 6.]]]

Query 2 : {'latitude': (1, 3), 'longitude': (3, 6)}
Coords 2: Coordinates:
  * time         (time) datetime64[ns] 1970-01-01
  * latitude     (latitude) float64 1.5 2.5
  * longitude    (longitude) float64 3.5 4.5 5.5
    spatial_ref  int32 4326
Values 2: [[[1. 2. 3.]
  [4. 5. 6.]]]

I could not find anything regarding this topic in the documentation, the issues, and gis.stackexchange.
Do you need any additional information?
The script to create these outputs is published as gist. It requires a running datacube environment, psql installed, and the delete_odc_product.sql file. The script performs the following steps:

Create some files:

netcdf file
datacube product specification
datacube dataset specification
datacube.conf

Init datacube database if not done already
Add the product if not done already
Add the dataset if not done already
Loads the dataset
Print metadata and data
Remove product from datacube using psql and the above mentioned delete_odc_product.sql
Removes all created files

Versions

os: Ubuntu 20.04.2 LTS
python: 3.8.5
datacube: Open Data Cube core, version 1.8.3


texts: Problem measuring area with gps coordinates

texts:
I am trying to measure area from gps coordinates but fail to get the correct result (7200 square meters). I used the formula from this post but I get way too high number as result (should be in square meters). Here is my current work in jsfiddle. Any idea what is wrong?

Math.radians = function(degrees) {
  return degrees * Math.PI / 180;
};

area = Math.radians(x2 - x1) * (2 + Math.sin(Math.radians(y1)) + Math.sin(Math.radians(y2))) + Math.radians(x3 - x2) * (2 + Math.sin(Math.radians(y2)) + Math.sin(Math.radians(y3))) + Math.radians(x4 - x3) * (2 + Math.sin(Math.radians(y3)) + Math.sin(Math.radians(y4))) + Math.radians(x5 - x4) * (2 + Math.sin(Math.radians(y4)) + Math.sin(Math.radians(y5)));

area = Math.abs(area * 6378137.0 * 6378137.0 / 2.0)
alert(area);


texts: How to sync Airship legend component data with CARTO VL legend data

texts:
I am trying to sync Airship category legend data as a VueJS component but it looks like the data property is not digesting the values properly:

<legend orientation="vertical" data="[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object]"></legend>

This is my pseudo-code. I have a legend VueJS component:

Vue.component('Legend', {
    'template':
    `
    <as-legend-category
        ref="widget"
        :bind="$attrs"
    >
    </as-legend-category>
    `,
    props: {
        data:{ type: Array }
    },
    data: function() {
        return {
            widget: null,
            isLoading: true
        }
    },
    watch: {
        data: function(newValue){
            this.isLoading = false;
            this.widget.data = newValue;
        }
    },
    mounted: function() {
        this.widget = this.$refs.widget
        const widget = this.widget;
    }
});

This is how I load it on my HTML:

<legend orientation="vertical" :data="computedLegendData"></legend>

I got the legend data from CARTO VL when the layer is updated:

layer.on('updated' , function () {
  vm.legendData = layer.viz.color.getLegendData().data;
});

And this is my vue-config.js:

// ignore airship tags
Vue.config.ignoredElements = [/as-\w+/];

const vm = new Vue({
    el: '#app',
    data: {
        legendData: [],
        ...
    },
    computed: {
        computedLegendData: function() {
            const data = [];
            for (let index = 0; index < this.legendData.length; index++) {
                const element = this.legendData[index];
                data.push({
                    color: this.rgbToHex(element['value']),
                    label: element['key'],
                    type: 'point',
                    width: 16
                })
            }
            return data;
        },
        ...
    },
    methods: {
        rgbToHex: function (color) {
            return "#" + ((1 << 24) + (color.r << 16) + (color.g << 8) + color.b).toString(16).slice(1);
        }
    },
    ...
}

As you have noticed I have to calculated a new legend data object because what CARTO returns is not enough and colors need to be transformed to hex values.

You can check the code here.


texts: Geospatial One-Stop no longer accessible?

texts:
I was wondering what happened to the "Geospatial One-Stop" site?

The link:

http://www.geodata.gov

is no longer accessible.

Is the portal superseded by another government website, or is it abandoned?


texts: PostGIS ST_Distance in meters discrepancy

texts:
Could someone please explain why there is such a significant difference between the results of this query:

SELECT
    ST_Distance(ST_Transform(previous_geom,3857), ST_Transform(current_geom,3857)) AS distance_transform_exp,
    ST_Distance(previous_geom::geography, current_geom::geography) AS distance_geometry_exp
FROM
(
    SELECT
        ST_SetSRID(ST_MakePoint(-111.9096893, 40.7411742),4326) AS previous_geom,
        ST_SetSRID(ST_MakePoint(-111.9092079, 40.74135181),4326) AS current_geom
) points_table

Results:

 distance_transform_exp | distance_geometry_exp
+-----------------------+-----------------------
       59.6050787601899 |           45.19182216
(1 row)

I understand there would be a difference between geometric (Cartesian) calculation results and geodetic (spheroid) alas for the example points which are relatively close, the discrepancy would be expected to be negligible.


texts: ATS77 MTM 5 to NAD83 CSRS UTM 20N Coordinate Converter?

texts: I've been looking for a simple way to check the accuracy of my transformation using a set of coordinates for a single point but have been unable to find an online tool or piece of desktop software that can do the ATS77 MTM 5 to NAD83 CSRS UTM 20N conversion since Franson CoordTrans and GeoSuite Calc are unable to do this. Any suggestions?


texts: How can a truck network be built / extracted using osm2po?

texts: Is it possible to adjust the config of osm2po to build a truck / lorry (in German: LKW) network?


texts: Mosaicking landsat imagery before LULC classification

texts:
I have to perform change detection of several years for Central Indian Landscape. I have corrected the atmospheric aberrations of the Landsat 8 imageries using the semi-automatic classification plugin of QGIS and then mosaicked them using the merge tool. One problem I am seeing is that the imageries have differences in spectral signatures for the same land feature (say, for example, agriculture area).

Would it be better if I do the land use land cover classification first and then merge the classified layers?

Please note I will be doing supervised classification in ERDAS Imagine 9.2.


texts: I am looking for an alternative solution to openlayers for satellite imagery in QGIS?

texts:
Does anyone have any other suggestions for sources of satellite imagery on QGIS? Openlayers is becoming extremely temperamental in QGIS 2.6 (Google imagery not bothering to work anymore, Can't zoom in below 1:2000 on all layers, Doesn't work at all with print composer). There must be another solution to having a decent satellite image?

Can anyone point me to any alternative sources or solutions?


texts: OSM data convert to old microstation dgn format

texts:
I'd like to convert OSM data into dgn file (version J), with my projection. I'd like you it as a reference file. Is this possible?

I try to save it in QGIS. Is this solve with Microstation V8 version program?

When I open it in Microstation J, the dgn content show on the zero coordinates.

What is my mistake?

In QGIS, I save file in Microstation DGN, after that I open it Microstation J but the projection is not good.


texts:
I created a cost map from an elevation map an a (degree) slope map using r.walk (version released on 2015-03-12) in QGIS 2.14.3. According to the documentation the resulting map should display the time needed to reach each raster cell from my starting point. As far as I understood, the unit of time should be seconds. I wanted to mask the area reachable in 3600 seconds, but it is only about 800 m wide on more or less flat terrain.

Could the unit of time be something else than seconds? Or is there some other error in my thinking.


texts:
We have been using Google Maps V3 during the development of a web service. We've been struggling with getting the maps integrate good with the design, since the maps look kind of grey. Now when i look at Apple Maps, they look so much nicer (in standard mode that is). The colors are brighter, and the maps doesn't look "grey and boring" like Google Maps do.

My questions are,
1. Is there any way to get Google Maps to look like Apple Maps do?
2. If not, is there any (preferably free) maps we could use that looks something like Apple Maps?


texts: Arcpy Update Cursor (arcpy.da.UpdateCursor()) not saving changes to 'SHAPE@LENGTH'

texts:
Looked around and found other Cursor questions, but nothing dealing with updating 'SHAPE@LENGTH'. I am running ArcGIS 10.2.1 and have a multipoint shapefile that I am grouping into polygons based on a shared attribute. However, sometimes only one point has the attribute and so a very tiny polygon is drawn, as described by the minimum bounding geometry documentation (http://resources.arcgis.com/en/help/main/10.2/index.html#//00170000003q000000).

I wanted to loop through the new polygon shapefile and use Update Cursor to access the tiny polygons' Geometry ('SHAPE@LENGTH') and bump it up to a visible level. I am unfortunately not understanding something because the changes to the length are not saved. I figured it would be harder than just reassigning the area value, but so far I can't understand what is going wrong:

import os
import arcpy

# Point data directory
polyFiles = 'C:\Users\...polyFiles'

# Polygon file made from the points
boxes = os.path.join(polyFiles, 'OE_AK_poly.shp')

# If length is less than 1000m bump it up so it's visible on the map
with arcpy.da.UpdateCursor(boxes, ["SHAPE@LENGTH"]) as cursor:
    for row in cursor:
        if row[0] < 1000:
            row[0] = 1000
        cursor.updateRow(row)


texts: Is there a way to Collect Lat-Long in a Google Docs Form?

texts:
I'm helping out a small non profit. They are collecting some simple data using a Google Docs Form. The need to collect exact latlong (upto 5 meters accuracy is fine).

Ideally the end user would click on a map and that lat-long would be submitted along with the rest of the form.

I have seen the mapbox blog talk about this here. But that is using the mapbox API.

I was wondering if there exists a ready made template to collect the point on top of Google Maps. I don't want to spend half a day converting the mapbox sample to Google Maps API, and then find out that something already existed.

So Is anyone aware of template where a point was collected from Google Maps in Google Forms?


texts:
Possible Duplicate:
  Access-Control-Allow-Origin error in web map.

I am using Drupal and geoserver for my project.
Drupal is running on localhost:80 and geoserver on localhost:8080. I want to send data from drupal to geoserver, and get it back to drupal to display it on a map. When I send the data to geoserver through the GetCapabilities URL I am not able to see the points in Layer->preview and the GML file is empty.

When I view the error in firebug its giving:

  XMLHttpRequest cannot load
  http://localhost:8080/geoserver/wfs.
  Origin http://localhost is not allowed
  by Access-Control-Allow-Origin.

can anyone tell me what this means and how to fix it, thanks.


texts: Normalize and Results show as Percentage

texts:
I've Normalized my data and everything looks okay. However, its appearing (on my map legend and in the Table of Contents) as a ratio (Ex: 0.000-0.005) but I would rather it appear as a percentage (Ex: 0%-5%).
As an attempt at a solution, when I Right click on the layer property and Open Attribute table, there's no field in the table as NORMALISED DATA. Therefore, I'm having problems right clicking on that field, going to Numeric and changing it to percentage.

Note:
The process I used to generate my Normalised data (which is been shown as ratio) is this:

I went to symbology>Quantities>Graduated Colors and selected my first field for Value and my other field for Normalization. My Classes is 5.

I am using ArcGIS Desktop 10.3.1.


texts: Transferring Z values from point layer to line layer using ArcMap?

texts:
I have two shapefiles, point and polyline, and I am trying to transfer Z values from point layer to polyline layer to calculate the pipe slope (by using manhole's depth and z value).  In the attachment you can see the attribute table of each layer.

What procedure should I use?


texts: Retrieve value from another layer by field calculator

texts:
I have two vector layers, one point layer named "element_ep", and another layer line named "reseau_ep".

In the layer "element_ep" there is a field 'nu_affaire' which contains the business number of each entity.

I want to retrieve this business number (from the field calculator) by creating a 'nu_affaire' field in the 'reseau_ep' layer.

The Qgis #refFunctions # plugin is well suited for this operation, but I want to get the code expression to get there.

Here are the layers and attribute tables.


texts: How to open an output set that is a list of double values in a geoprocessing service

texts:
we are trying to reduce the amount of data we are returning from the server by eliminating returning a polyline and instead returning a collection of values of type GPDouble.

When the transaction is kicked off inside our GeoProcessing Service via Execute, the first thing we do is create the parameters we are either going to process or return.  for example to get the target list we execute this code:

        IFeatureClass targetClass;
        filter = null;
        _gputilities.DecodeFeatureLayer(GPUtilities.GetParameterValueByName(parameterValues, "TargetList"), out targetClass, out filter);
        if (targetClass == null || targetClass.FeatureCount(null) == 0)
        {
            message.AddError(2, "Could not open Target input dataset.");
            return;
        }
        message.AddMessage("Target List has been defined.");

I changed my OutList parameter definition as follows (old definition commented out)

            //IGPParameterEdit3 outputPoints =
            //    GPUtilities.CreatePolylineRecordsetParameterEdit("OutputList", "Output List", "MsramSDE.DBO.Standard_Polyline",
            //        esriGPParameterDirection.esriGPParameterDirectionOutput,
            //        esriGPParameterType.esriGPParameterTypeDerived, new GPFeatureRecordSetLayerTypeClass(), true);
            IGPParameterEdit3 outputPoints = GPUtilities.CreateParameterEdit("OutputList", "Output List",
                esriGPParameterDirection.esriGPParameterDirectionOutput,
                esriGPParameterType.esriGPParameterTypeDerived, new GPDoubleTypeClass(), true, "0.00");
            parameters.Add(outputPoints);

None of our other GeoProcessing Services return anything but a FeatureLayer.  So I lack any example and could not find one via Google for opening up an Output Item that gives me a way of writing out the double values.  Anyone ever done this and how do I do it?  I've looked through the GPUtilities calls to see if there was an obvious DecodeBaseType or some such call but could not find one.  Should this be defined as a record set instead?


texts: Selecting rows where values of two columns are same in ArcGIS for Desktop?

texts: I would to have your guidance on a type of data selection in ArcGIS. I would like to select the rows of a feature "where the values of two columns are the same". as a standard approach there is no such a thing in ArcGIS and I am wondering how it is possible to do it.


texts: Clipping polygon layer according to DEM elevation using QGIS?

texts:
I have a polygon layer on top of a DEM.  Is there a way to clip the polygon layer according to elevation, ie the clipped result being only those areas of the polygon layer between x and y altitude?

The result needs to be a vector layer.


texts: Tuning PostGIS for production environments?

texts:
What PostGIS documentation have you come across that has been helpful when tuning PostGIS for a production environment?

I would like to get together with my DBA to setup a Postgresql/PostGIS installation that is fit for production. I have read that there is some tweaking involved to achieve this, and I was hoping to find the answer on the refractions web site.

So far I have found some documents on the OpenGeo site helpful, like this one.

And this old forum post is the kind of information I have found helpful, this is probably just basic DB stuff but to me its good stuff.

I would be interested to see what resources have helped others in achieving a stable production installation of PostGIS.


texts:
I have two lines

I want to disole where the two lines intersect to become like this

How to do that in QGIs 3.16?

I digitizing two streets. And I modify the width of the lines.


texts: I can't find my content layers in my web scene on ArcGIS Online

texts:
I am new on ArcGIS Online. I started uploading geographic layers; in my webmap I can load them when I search for layers in my contents, like shown below

But in my webscene I can't find any layers?


texts:
We have an existing postgis database. Our application is developed in Java, but queries database through a jni dll, which appears to be not very effective for our purpose.

I've managed to migrate to jdbc, and it's working well, but the next step would be to rely on a framework such as hibernate-spatial, which i don't know at all.

How mature is it? I've seen some people here using it, what do you think of it?

It looks that it supports only cartesian coordinates and it may be a concern for us for distance computation.


texts: For each cell with value 1 in a raster, create a raster with only that value

texts:
I have a binary raster in R. It has only 55 cells with value 1 and all the others have value 0 or NA. Now I want to write 55 rasters. For each cell with value 1, I want to write one raster. In these final rasters I only want to have one cell with value 1. My Idea is to loop through all the cells of the input-raster and in case it finds a cell with value 1 write a new raster with only that cell beeing of value 1 and all others of value 0 or NA.

The input raster looks like this:

class      : RasterLayer
dimensions : 1482, 1274, 1888068  (nrow, ncol, ncell)
resolution : 12.5, 12.5  (x, y)
extent     : 347279.6, 363204.6, 6255207, 6273732  (xmin, xmax, ymin, ymax)
crs        : +proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
source     : memory
names      : raster
values     : 0, 1  (min, max)

My idea would lood kind of like this:

for (i in 1:length(raster){
    if (raster[i] == 1){
    name = paste0("rastercell_", i, ".tif")
    ras = raster(...)
    writeRaster(ras, name)
    }
}

It's just that I don't know how to create the raster (ras = raster()) from scratch with only that single cell being of value 1 and all others 0 or NA.


texts: How is geometry loaded into QGIS Midvatten plugin obs_lines table?

texts:
I have created a csv file that contains the coordinates for a series of polyline vectors and I am trying to load it into Midvatten's obs_line table. I have made several attempts at this but the data gets loaded using all fields except the geometry field. I first tried the layout given in the tutorial example data set:

WKT;obsid;name;place;type;source
"LINESTRING (136518.11581965966615826 6972638.38746192120015621,136526.23101053971913643 6972634.87753117084503174,136527.8896709508553613 6972634.69279802031815052,136531.11793618529918604 6972632.84546670503914356)";BE1;BE1;Beeston;line;connect_obs_points

I used the import dialog box that is summoned by 'Import data to database' > 'Import data using general csv format' > load from file and chose UTF-8 and then file and then chose the destination as obs_lines. All 78 rows were loaded minus the geometry. I deleted the data and tried again changing the first item on the header from 'WKT' to 'wkt_geom'. I had seen 'wkt_geom' in the obs_point attribute table and thought it might be the correct form but the geometry still failed to load.

So finally I deleted the data once more and rearranged the csv table to match the column list using the column name 'geometry':

obsid;name;place;type;source;geometry
BE1;BE1;Beeston;line;connect_obs_points;"LINESTRING (136518.11581965966615826 6972638.38746192120015621,136526.23101053971913643 6972634.87753117084503174,136527.8896709508553613 6972634.69279802031815052,136531.11793618529918604 6972632.84546670503914356)"

This time I no longer got the warning that I had loaded data without geometry, so I thought it may have worked. But the attribute table did not appear to have geometry in it. I wondered if the geometry was hidden from view so I tried running a section plot and got the following python error:

if geom.wkbType()==QGis.wkbLineString:#...and that the activelayer is a line vector
Attribute Error: 'NoneType' object has no attribute 'wkbType'

I assume that this is telling me that there is no geometry.

How do I need to change the csv file to get the geometry to load? I guess I could try rounding the coodinates to 1 or 2 decimal places. The large number of decimal places resulted from converting from EPSG 4326 to EPSG 3857.

I have now rounded all the coordinates to 2 decimal places and retried loading with the three different header names as previously ('WKT','wkt_geom' and 'geometry') and did not manage to load the geometry. So have now ruled out rounding the coordinates as a possible solution.


texts: Why won't QGIS show markers outside the window?

texts:
I'm having an issue with the QGIS2.18.3: I have grid center points, where I want to put square shaped markers to create a grid. The problem is that if a center point is outside the window, the program won't draw the marker - leaving that part blank.

What would you suggest? Can I force somehow QGIS to draw everything despite it would not see it?


texts: Applying batch process with ModelBuilder in ArcMap?

texts:
I have done a process in a Model Builder in ArcMap 10.5 version and in order to execute it once over a lot of rasters I have used raster iteration.

It works well, but I cannot  specify differents parameters to each raster such as do a clip with differents shape.. So I want to know if can we use the batch option in ModelBuilder.

I add a picture with a part of my model, if it is another way to do the iteration you can suggest me ( the labels are in spanish but only is important the iteration -hexagon shape box-).


texts: Using Python Parser in field calculator to update column from another column

texts:
I am an experienced ArcMap user who is fairly new to Python. I am having trouble getting the script to look for a list of values and assign a new label in a new column. The script runs without errors but it does not fill in the new column.
Do the two columns need to be the same type?
I have provided the code that I use in the code block of the field calculator. In the other half of the field calculator I have
def reclass(Landing_Factor):
    if Landing_Factor >=1.0 and Landing_Factor <=1.99:
      return "A"
    elif Landing_Factor >= 2 and Landing_Factor <=2.99:
      return "B"
    elif Landing_Factor >= 3 and Landing_Factor <=3.99:
      return "C"
    elif Landing_Factor >= 4 and Landing_Factor <=4.99:
      return "D"
    elif Landing_Factor >= 5 and Landing_Factor <=5.99:
      return "E"
    elif Landing_Factor >= 6 and Landing_Factor <=6.99:
      return "F"
    elif Landing_Factor >= 7 and Landing_Factor <=7.99:
      return "G"
    elif Landing_Factor >= 8 and Landing_Factor <=8.99:
      return "H"
    elif Landing_Factor >= 9 and Landing_Factor <=9.99:
      return "I"
    else:
      return "N/A"

BENCH_1 = TextValue( !BENCH_1! )


texts: Does the psql prompt have a working directory in the file system?

texts:
I am learning postgis. This tutorial says "Navigate to where you downloaded the shapefile." Does that mean navigate to this location in the #psql prompt or navigate to this location with the $propmt? In other words: does the psql shell have a working directory in the file system that you can change with cd, like in *nix?

I am confused because the tutorial says to run shp2pgsql after navigating to that location -- but shp2pgsql is a postgres command, not a bashcommand. But the sql prompt lists tables, records etc -- which does not have to do with the file system. So does the sql prompt know about the file system or not?


texts: Count number of grids intersected by multiple polygons in layer

texts:
I am using QGIS. I would like to know the number of 10x10 grids that intersect with each polygon of my layer.

Layer 1: This layer contains a list of species. Species A, Species B, etc. Each species has its own polygon or multiple polygons within the layer depicting its geographic distribution across the globe.
Layer 2: I have a second layer, a grid of 10x10 cells that I created using:
Vector-> Research Tools-> Vector Grid.

I want to know the number of grids (Layer 2) that intersects with each species polygon(s) (Layer 1). I know that if Layer 1 had only one species, then I could use Vector->Data Management Tootls-> Join Attributes By Location.

However, I do not know how to loop through Layer 1 so that it does this for each species.

This link seems to be asking for something similar, but the solution uses ArcGIS, and I am using QGIS.


texts:
I'm thinking of using Openlayers for an interactive map of a building. I'll have custom tile layers for each floor of the building's interior overlaying regular world map tiles.

The building is on a bit of an angle, though, with respect to north, and the client wants it to be displayed with the walls vertical and horizontal. This means rotating north to be 10 degrees or so from the vertical. This is what ruled out the Google Maps API.

I know Openlayers can rotate a map (example at http://openlayers.org/en/v3.1.0/examples/rotation.html) but when rotated the tiles don't look so good since they no longer match the screen pixels 1:1.

For that reason, and also for ease of drawing, it'd be good if I could make the map tiles already rotated (horizontal and vertical walls) and then somehow make Openlayers aware of this.

Then when the rotation of the map is introduced the tiles would be in their native orientation once again.

Is this possible with Openlayers? How would it be done? Or must I make the tiles with north being up?

The alternative would be to abandon real world coordinates and orientation completely, which would reintroduce Google Maps as a possible technology.


texts: Empty fields after joining attributes by location to join point to polygon layer

texts: I'm trying to join a point layer to a polygon layer to get a specific field in the attribute table into the polygon layer. I want the polygon layer "Zone" to have the field "Improved" from the point layer "Permits". Is there a way to join based on location? All the points in "Permits" are located within the respective polygons in "zone". I've tried performing a Join Attributes by location" but the fields in the new layer end up being blank.


texts: I'm trying to create an introductory GIS course that will involve an ideally simple web mapping component. I was underwhelmed that ArcGIS online doesn't support rasters, and surprised that MangoMap and GeoCommons don't appear to either at the moment. Is anyone familiar with a webmapping site with a fairly simple UI that allows the uploading of raster data in addition to vector?


texts: Downloading multiple files from TIGER FTP site

texts:
Hello everyone! I've been downloading TIGER data per county from the FTP site one at a time, with some states having hundreds of files. I was curious if there was a faster way to download maybe all the files for each county for each state at once? Thanks for your time and help!

Attached is a picture of what i'm dealing with, would like to eventually download all the files for each state at once versus one at a time.


texts: Using polygons to apply unique value to all raster cells within those polygons?

texts:
In ArcGIS 9.3, I have a raster file that I would like to "mask" with polygons, giving a single value to all the raster cells that overlap the polygons.

How can I do this?

Here's a screen shot, where the black polygons are my polygon layer, and the rest is raster-based.


texts: Extracting data, not colors, from geotiffs in R

texts:
I'm an experienced R user new to GIS, but finding my way. I've recently run into an issue that I haven't been able to find any satisfactory documentation on (but could just be searching poorly). I'm attempting to extract data from a geotiff. However, when I load said file, all I can get are color values.  For example:

aMap <- raster("./tiff/mymap.tiff")
extract(aMap, areaCoords, buffer=2000)

[[1]]
[1] 247 247 240 247 240 247 247 247 247 247 240 240 240 247 247 247 240 240 247
[20] 247

where areaCoords is a spatial points data frame. I'm using the sp and raster libraries. Now, I can add a layer argument to raster, and get other channels (I'm assuming this goes in RGB order), but, the data is a single value in different units alltogether.

1) Is there something I'm missing in loading in a geotiff for getting data rather than just colors?

2) If all there is to extract are colors, any suggestions on converting from colors to a known scale (e.g., perhaps I can get the RGB of the low and high values of the scale - perhaps).


texts: Exporting results of ArcMap Data Reviewer to shapefile

texts:
I am new to the ArcMap 10.8.1 Data Reviewer extension and trying to use it to identify slivers.
I can produce a data reviewer table, but how do I download the results as a shapefile?
When I run the sliver check on polygons I get the message “Unable to save Reviewer geometry for RecordID 1195 as its spatial reference is incompatible with Reviewer workspace spatial reference. Datum of geometry is not equal to Reviewer Dataset datum” for all the records, but they are still written to the table.
These records do not show as selected in ArcMap and I can’t view the selected records which I need to do.


texts: How do you solve the map publishing problem with GeoServer 2.1.1 on Windows 7 x64?

texts:
Well, my problem is that I have GeoServer installed over a Windows 7 x64 system with both Java (JRE 6u29) x32 and x64 installed and I can't project my layers on the website. The preview on GeoServer works fine, but it doesn't work on my site (it's a localhost server, a testing server)

I have Apache + OpenLayers + GeoServer.

I have tried Tomcat + OpenLayers + GeoServer (war) with similar results.

The server status shows JAI libraries marked as false. I say this because I think the problem may be on the JAI library. I searched about and found that JAI works only on x32 systems so I don't know how to solve this. I need to work with x64 architectures.

I tried to install JAI over the x32 version of Java but my server stop working after that.

My local call to GeoServer through OpenLayers Instance is this:

<script type="text/javascript">
    var mapa, capa;

    function init(){
        mapa = new OpenLayers.Map('mapwms');
        capa = new OpenLayers.Layer.WMS(
            "Catastro",
            "http://localhost:8080/geoserver/Catastro/wms",
            {layers: 'Catastro:Predios'}
        );

        map.addLayers([capa]);
        map.zoomToMaxExtent();
        map.addControl(new OpenLayers.Control.LayerSwitcher());
    }
</script>

So I need to figure out why the GeoServer preview is working and not my code.

I'm new on this kind of software so I need some help

I put "http://localhost:8080/geoserver/Catastro/wms" because that is how appear on the preview source code.

I'll really appreciate your help please ;)

NOTE: I've changed the question title to a more clear one ;)

SOLVED: Thanks all for your help. I follow your advice and my server is now working.


texts:
I have this error in R:

  Error in GetProducts() : Application error

  In addition: Warning message: package ‘MODISTools’ was built under R version 3.2.3

Can you help me?


texts: disable onEachFeature at zoom level, leaflet

texts:
I have a leaflet map that zooms and fits on the bounds of the polygon that is clicked via onEachFeature.

Is there a way to disable this at a certain zoom level? I add new layers on the 'zoomend' event, I just can't figure out if disabling on 'click' function is possible.

function onEachFeature(feature, layer) {
        layer.on('mouseover', function() {
            this.setStyle({
                'fillColor': 'gray'
            });
        });
        layer.on('mouseout', function() {
            this.setStyle({
                'fillColor': '#fff'
            });
        });
        layer.on('click', function() {
            map.fitBounds(layer.getBounds())
        });

    }

var json_countiesJSON = new L.geoJson(json_counties, {
        style: doStylecounties,
        onEachFeature: onEachFeature
    }).addTo(map)

map.on('zoomend ', function(e) {
        if (map.getZoom() < 10) {
            map.removeLayer(clusterWells)
            json_countiesJSON.setStyle({
                'fillOpacity': '0.7'
            });
        } else if (map.getZoom() >= 10) {
            map.addLayer(clusterWells);
            json_countiesJSON.setStyle({
                'fillOpacity': '0.0'
            });
        }
    });


texts: Getting list of result from GeoPy Nominatim

texts:
How can I get list of possible result in Nominatim?
My code is
geolocator = Nominatim(user_agent="my_email@myserver.com")
location = geolocator.geocode("School")
location.raw

but unfortunately it only shows one result
geolocator = Nominatim(user_agent="my_email@myserver.com")

location = geolocator.geocode("School")

location.raw

{'place_id': 174509566,
 'licence': 'Data © OpenStreetMap contributors, ODbL 1.0. https://osm.org/copyright',
 'osm_type': 'way',
 'osm_id': 353533683,
 'boundingbox': ['35.5912514', '35.594524', '-78.8117411', '-78.8064663'],
 'lat': '35.5928655',
 'lon': '-78.80864240000005',
 'display_name': 'School, Fuquay-Varina, Wake County, North Carolina, United States',
 'class': 'place',
 'type': 'neighbourhood',
 'importance': 0.36}

But I wanted to get list of possible result.


texts: Extracting NDVI over time using Google Earth Engine

texts:
I am trying to extract NDVI data for provinces in Cambodia over time. The following codes work for a specific year (2002).

Can someone show me a better way to collect the data over time automatically on GEE?

var collection = ee.ImageCollection('LANDSAT/LE07/C01/T1_32DAY_NDVI')
.filterDate('2002-01-01', '2002-12-31').select('NDVI');

var mean = collection.reduce(ee.Reducer.mean());

var meanDictionary = mean.reduceRegions({
  collection: Cambodia,
  reducer: ee.Reducer.mean(),
  scale: 30,
});

Export.table.toDrive(meanDictionary,"NDVI","NDVI","NDVI_Cambodia");


texts: Spatial join in Geopandas is only finding itself

texts:
I ran my code successfully last week, cleared the kernel and now it does not work. It must be something simple I'm missing, but I cannot seem to puzzle it out. The purpose of the code is to find all of the polygons inside each of the polygons with a hole. In the data I provided here, there are two polygons with holes in the test_donuts_1.shp file and there are three polygons (including the two polygons with holes) in the test_poly1.shp file.

    import geopandas as gpd

    test_poly1 = gpd.read_file('test_poly1.shp')
    test_donuts_1 = gpd.read_file('test_donuts_1.shp')
    intersections_df = gpd.sjoin(test_poly1, test_donuts_1, how='inner', op='within')

This code will run without error for me, but the output is incorrect. It will say that they are only contained within themselves. But I know that if you union the three polygons in test_poly1 together, you get a single polygon. What I mean is that there are two polygons inside one of the test_donuts_1 polygons and there is one polygon inside of the other test_donuts_1 polygon, but this is not caught by sjoin. I have tried setting the crs for both geometries to be the same in my main file even though it is not reflected in the code above.

Can anyone see where I went wrong with this?


texts:
When I'm serving up maps and UtfGrids with tileStache I would like to also include some metadata with the layers so that my leaflet application can understand the context of the layer. For example: category, data, title, etc.

What is the best way to handle metadata with TileStache?


texts:
Specific question on Fill tool in ArcGis 10.2.2, I got a filling problem, I think holes are not fully filled when I use the fill tool, how do I change that ?

2 images of the same area for an example, in green, manually delimited sinkholes

Image 1: flow accumulation, you can see there is a concentration of flow in 1 point

Image 2: difference between fill tool and original DEM, you see in the same red circle there is no fill why ? how do I change that ?

I changed Z limit without any change


texts:
I've enabled the plugin qgis2web, after getting one WFS layer in to QGIS I started the gis2web plugin. Because it is a WFS layer and I want leaflet to use that WFS I don't check the Encode to JSON option.

The preview shows only the base layer. If I do enable the Encode to JSON option, my layer does show.

I've tried with different WFS-services in QGIS 2.10.0

What am I doing wrong?


texts: Looping through raster files to create attributes using dictionaries in ArcPy

texts:
I have a lot of output land cover raster files in .tif and want to only select the .tif files and change color map of all of them by taking a script that has attributes for Red, Green, Blue and LC Type. I have a list of raster values for colors and the associated land cover type; however, I know I am supposed to take the proportion of each RGB value out of 255 to properly assign the new color. I want to loop through each .tif file create the fields (Red, Green, Blue, and LC Type) and assign the appropriate value using a dictionary in ArcPy.

filepath = 'C:/directory'
env.workspace =filepath
for filename in os.listdir(filepath):
  if filename.endswith('.tif')

dictionary = {
    1: [255,0,0 'water'],
    2: [209, 255, 115 'Evergreen Needleleaf Forest'],
    3: [85,255,0 'Evergreen Broadleaf Forest'],
    4: [112,168,0 'Deciduous Needleleaf Forest'],
    5: [190,255,232 'Deciduous Broadleaf Forest'],
    6: [255,211,127 'Mixed Forest'],
    7: [0,255,0'Woodland'],
    8: [0,99,0 'Wooded Grassland'],
    9: [214,157,188 'Closed Shrubland'],
    10: [197, 0, 255'Open Shrubland'],
    11: [0,77,168 'Grassland'],
    12: [0,132,168 'Cropland'],
    13: [190,210,255 'Bare Ground'],
    14: [255,170,0 'Urban and Built-up'],
    17: [255,170,0 'Urban and Built-up']}
#create the four attributes, R,G,B, and LC
#create an update cursor
#loop through each .tif file create fields and assign values to records


texts: converting a python list to arcmap table

texts: I have values stored in a list that need processing in the costDistance tool in ArcMap 10.2 (desktop). For this they need to be converted into a table or file readable by arcmap's tools. Anyone with experience doing this?


texts: Are there any news about the release date of QGIS3

texts:
Since I will make a workshop on the new features of qgis 3, the members will ask the release date of QGIS3.
Having a look to
https://www.qgis.org/de/site/getinvolved/development/roadmap.html
only shows 2.18 until the next year, QGIS 3 has completly disapeard from the road map until 2019...?


texts: Integrate Google Streetview map into QGIS?

texts:
I'd like to show a the Google Maps layer that has the routes that streetview has taken in my QGIS 2.18.x project.  I'm not looking for the actual street view imagery.  I've seen multiple posts about including the imagery, but I'm interested in tiles. :)  I have QuickMapsServices plugin installed, and was hoping there was an easy way to integrate something.

Here's a screen shot of what I'm looking to duplicate.


texts: Counting points in polygon and extracting the max count of an attribute

texts:
My first query on here. I'm using QGIS.
I've created a gridded map by counting points in each grid:

This was generated by counting the number of observation points per grid(polygon).
The attribute table in the observation (polygon) layer looks like this:

But I now wish to create a similar gridded map, but instead of just counting the number of points within each grid, I wish to count the contribution of each author_id within each grid and add a column which shows the highest contributor within each grid.

Is this doable?


texts:
I have recently purchased a Phantom 4 pro, and am reliving my 35 year ago survey life, but am having problems with QGIS!

I have both jpeg and tiff ortho images saved by Agisoft Standard which appear to be OK! I have the openlayers plugin and georeferncer plugins installed in QGIS and both appear to be working.

I live in Western Australia so I have set both defaults to epsg.28350 94 zone 50.

I can load the google background and then can load the georeferencer image, both jpg and tiff. I do the gcp's and hit the transformation button with the button checked to install the image in QGIS. The image then appears to be installed in QGIS, but all I have is a blank white background with the red gcp' points on the geotiff image in QGIS. However when I look in my windows file manager I can see a distorted image

What am I doing wrong?


texts: Could not load file or assembly ESRI.ArcGIS.ADF.Local?

texts:
I just tried to register my extension for ArcGIS 10 in the 10.1 beta and it failed when I tried to register the dll:
"Could not load file or assembly 'ESRI.ArcGIS.ADF.Local' Version=10.0.0.0."

Do I have to re-compile it against 10.1 dll-s in order to make it work?
For the 9.x version all I had to do was to compile them against the base level (9.2) and it worked in 9.3 as well.


texts: Composition of features into attributes field using QGIS

texts:
For every feature of my map, I have to:

I localize the feature with its coordinate X-Y
I zoom into the canvas in order to center the feature and to obtain the scale wanted (depending of every category of feature). Every group of feature has an optimal scale to take the screenshot. This data is available in a attribute field named : "scale_screenshoot"
I export as image the map to capture the PNG, into a folder named "screenshoots" at the root of the QGIS file.
I resize the PNG (depending of the optimale size for this category of feature). This data is available in an attribute fields named: "size_height_screenshoot" and "size_width_screenshoot"
I write the path of the screenshoot into the attribute field named "path_screenshoot"

Option : it would be nice if every png produced get a worldfile.

There are several hundred of features, is there a way to automate this task?
One of my final goal is to to deliver an Excel sheet with map features listed and a "screenshoot" (path of the png) for every feature.


texts:
What is the difference between QGIS for Windows and QGIS for Android?

Is the functionality the same?


texts: GeoServer Web UI Tables - Default amount of items to show

texts:
As it stands, all of the HTML tables in GeoServer only load 25 items at a time. Is it possible to set this to a higher number so I don't have to go through numerous pages to find my items?
Aware that the search bar exists also - am hoping there's a quick bit of code I can change in the config.


texts: Performing change analysis between two time series images (2010-2000)

texts:
I performed the difference between two images, 2010-2000, with the command in Arcmap10: Image Analysis and raster difference,based on the pixel value difference in the two image tables (2010, 2000). The pixel value in each image represents the land cover aggregated class, with 6 class at all, and I need to represent the change value in the table of the output image as a combination of the single values in the two original images, ex: image 2000 pixel value=1; image 2010 pixel value=3, the pixel value of the output image should be: value= 13, which means that the land cover classes changed from 1 in 2000 to 3 in 2010.

How can I remap the pixel value in the output table like that?


texts: QGIS plugin: Problems importing resources (resources_rc) file - plugin doesn't load - PATH problems?

texts:
I'm building qgis plugin and I can't find solution for this error.

File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 478, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
ImportError: No module named resources_napoved_rc

For everybody that will be asking I've built python resources file:

pyrcc4 -o resources_napoved_rc.py resources_napoved.qrc

I still can't find a way to make it work. I always get same error.

On top of the script I have:

import resources_napoved_rc.py

I'm using ui file directly from qtbuilder.  Any ideas how to go forward ? I'm assuming this must be some kind of path problem or something similar.


texts:
Maybe it's a stupid question, but I don't know how to connect layer style with database in QGIS.

When I save the style to postgres database in QGIS, it automatically saves the style in the table layer_styles in public scheme.

But when I want to reconnect it to the layer, it shows only all other layer styles without required one and in the window 'styles related to the layer' there is a message "No styles found in the database". Do you know, where the problem could be?


texts: Projection issue when matching layer projected with SRC4326 with other layer on SRC3857

texts:
When I'm trying to visualize a multipolygon layer projected on SRC4326, on a base layer from Bing Maps (also tested with Google Maps), my polygons doesn't show up on the ocidental part of the map, as well on the antimeridian region, as shown in this image:

The project is set for a dynamic reprojection on SRC3857. Any hints on solving these two problems? This layer was originally a GeoJSON polygon collection, but now ingested on PostGIS.


texts: Create label from different fields in ArcGIS Desktop?

texts:
I want to create label like this one in ArcGIS Desktop 10.3.1.

Is there anyway to create label like this?

Each data come from different field.

I can create a label with box or a number with degree. But what I'm looking for is how to stack label like in the example with the automatic resize of the box depends on charactor width.

Can I create label like this in one label not from an annotation by typing and rotate manually?

Now I can create the label. By using the label expression and use split callout for text background to create the box but the problem is when I increase or reduce font size the box margin does not fit the charactor. Is there anyway better than this?

This what happen when I increase the font size.

when I reduce the font size, box doesn't fit the charactor

and my label expression is
[MFA] + vbnewline + [track01] + "° " + [IDENT] + " " + [track02] + "°" + vbnewline + [DIST]


texts: Generating a custom flat projection in pyproj

texts:
I have little experience with pyproj, but I want to achieve the following task:
I have a trajectory of a robot, given in lat/long coordinates which I want to convert to x/y Coordinates. And I want to do this using a "flat" cartesian coordinate system originating at a lat/long combination of my choice. I realize that this projection only works for a small area, but that is okay with me, as the area I am looking at is small (<10km in all directions.)
Here is what I achieved so far, which is a transverse mercator projection (or so I was told).

I have used the following code.

import pyproj as proj

crs_wgs = proj.Proj(init='epsg:4326')  # assuming you're using WGS84 geographic
    crs_bng = proj.Proj(init='epsg:27700')  # use a locally appropriate projected CRS
    # then cast your geographic coordinate pair to the projected system
    x, y = proj.transform(crs_wgs, crs_bng, s['gps_lng'], s['gps_lat'])

So this question is two-fold:

Which Projection do I use to achieve a local flat cartesian coordinate system
Howto define this projection at my given lat/long combination

I was told to ask this question in a new thread, see my old question. Also a similar question has been asked, but it seemed to be too unspecific to be answered.


texts: How to can I run the OpenGeo Suite Community Edition under Apache Tomcat on Windows server 2003?

texts:
I have installed opengeo suite community edition on ubunutu 10.10 and it is working with apache tomcat and it all works great. So I decide to do the same on a Windows 2003 server at work.

Note: On the windows 2003 server machine I already have apache tomcat installed and working correctly on port 8080

So what I did was download the opengeo suite community edition for windows from opengeo.org and run the installer. I change the default port to 8081 and everything works fine! However it installs using a different web server (JBOSS?? or J something?), and this server only seems to operate when I am logged onto the the machine which is undesirable. For example I can't access geoexplorer from a remote machine unless I am also logged onto the server.

Has anyone got the openGeo suite community edition running with Tomcat? and how did you do it?

Thanks
Ando


texts:
This problem is similar to retail presence optimization. I'd like to know which GIS's could help us do the following:

Given a set of point locations,

create a circle of radius n meters around each.
each circle has a unique score.
replace the circle's unique score with the sum of its score and the scores of all circles with which it intersects.
create a polygon(s)--the union of all circles with sum >= a threshold.
Layer in a new set of points and take the set of all these new points that are contained within the polygon(s).
Export as a table those points and their metadata.

I once upon a time was reasonably fluent with Manifold, and I think this could be done in Manifold, I think. We're interested in a GIS, for now, costing no more than Manifold.  Shortest learning curve super appreciated too.

Thanks in advance for any guidance you can provide,
Steve


texts:
I'm very excited about the new release of QGIS (2.0), and would like to install it on my Fedora 18 machine. I couldn't find any relevant RPMs, and the documentation on the QGIS main site for Fedora is pretty old (16/17).

Did anyone did it and can throw some tips?

Thanks!


texts: Selecting more features with "ILIKE" and "IN" operators using QGIS

texts:
I would like to select more features with expression in a specific field, which contains the name of the attribute (IN operator is useful). Moreover, the name could not be exact, so the ILIKE operator is necessary.
How can I match these 2 operators in one expression?
For example, I need to select 2 city called 'Rome' and 'New York'.
"city" in ('Rome','New York') -> selected both
"city" in ('rom','new yor') -> not selected

Any suggestion to select in the second case in QGIS?


texts: How to split a polygon into quadrants in R/Qgis?

texts:
I am solving the same problem as How to split a polygon into quadrants?, except that I am looking for a solution in R or QGis. So I have a polygon layer shapefile, which is unfortunatelly not exactly rectangular (due to some reprojections etc.). I would like to split the "squares" into 2x2, 4x4 or 8x8 smaller squares, as equal as possible.

Note that I am not looking for tool to generate a new grid, but to split the one I already have.


texts: Getting hillshade error "GDAL provider: Cannot open GDAL dataset" in QGIS?

texts:
I am using a DEM and using menu Raster > Terrain analysis > Hillshade. When I type a simple output file name "MyHillshade" and then click "OK", I get a two-part error message that says:

"GDAL provider: Cannot open GDAL dataset MyHillshade.tif: `MyHillshade.tif' does not exist in the file system, and is not recognised as a supported dataset name.

"Raster layer: Provider is not valid (provider: gdal, URI: MyHillshade.tif"

How can I repair this situation?


texts: Filling polygons with points or polygons using PyQGIS?

texts:
I am trying to find a method to, as efficiently as possible, fill irregular shape polygons with points (defined spacing) or alternatively, polygons of defined size ( which i can find the point centroid).

I am using QGis, Python


texts: Mapinfo, (or other?) Interpolation of lat/long and finding minimum in a region

texts:
I have a series of lat/long coordinates, each with a corresponding Z value. Each coordinate is spaced roughly 40km apart (this is stored as a MapInfo table of points).

What I'd like to do is to "refine" this bunch of coordinates so that I have interpolated values of Z every (say) 4km apart.

I have a separate table of regions, and I'd like to get the minimum value of my (refined) Z in each region. I hope this makes sense!

I have tried to:

create a thematic grid using my lat/long, but then what I want to do is find the minimum in each region using the thematic grid, and that doesn't work.
I cannot browse my thematic grid either (I assume because Mapinfo stores this as a raster?)

Am I going about this the right way? Can I do this in Mapinfo Prof.? Unfortunately, all my experience with MapInfo is doing a series of very repetitive tasks, and I'm a little outside my comfort zone here.

If I can't do this in MapInfo, what tool should I use?

Thanks,
JD


texts: Transforming Real-time Local position system coordinates into GPS/WGS84

texts:
I have a local coordinate system built with Ubisense sensors and tags (it is a real-time locating system used to automatically identify and track the location of objects or people). The local coordinate system is simply a cartesian coordinate system, which I can define the origin and axies at will except the rule that x-y panel must be horizontal, thus the z axis is vertical. After setting up the system, I can get coordinates of tags in the contained area where the sensors can get signals emitted from the tags and locate the tags. Now, I have the problem to convert the (x,y,z) coordinates into global system coordinates.

I made a model myself which took the origin of local coordinate system as a reference point and calculate the difference of other points in lat and lon to make the transformation. However, I guess there must be more standard and accurate way to deal with such problem. Since I am new in this field, I did google something and read some pages. However, there is rarely content relevant with this topic(much I found is from some standard to another, however is not a standard coordinate system, it is cartesian but not geocentric).

Can anyone give any hints or suggestions for this problem?


texts: How to create field that attributes what boundary a set of points reside in?

texts: I want to be able to assign a point layer a new attribute field with the polyline that the points are within. What method would work best in this scenario?


texts: Dissolving Points using QGIS Processing?

texts: Is it currently possible to create multipart features of points data based on one or more fields?


texts: gdalwrap incorrectly re-project wide angle

texts:
I'm trying to re-project Meteor-M 2 satellite image with gdalwrap. I have GCP points which I put into .vrt file:

<VRTDataset rasterXSize="1568" rasterYSize="480">
  <GCPList Projection="EPSG:4326">
    <GCP Id="0" Pixel="1568.5" Line="0.5" X="62.20476647354344" Y="42.09684997440528" Z="0.0" />
    ... lots of them
    </GCPList>
    <VRTRasterBand dataType="Byte" band="1">
      <SimpleSource>
        <SourceFilename relativeToVRT="1">output.png</SourceFilename>
        <SourceBand>1</SourceBand>
      </SimpleSource>
    </VRTRasterBand>
</VRTDataset>

Then I executed gdalwrap with the following arguments:

  gdalwarp -tps -overwrite  -t_srs epsg:3857 -of GTIFF output.vrt output10.tif

Then I cut this geotiff into tiles:

  gdal2tiles.py --profile=mercator -z 3-8 output10.tif tiles

So far so good. But when I draw both GCP and tiles onto open street map I'm getting:

Re-projected image and GCPs differ greatly.

gdalinfo:

Driver: GTiff/GeoTIFF
Files: output10.tif
Size is 1656, 583
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"],
    AUTHORITY["EPSG","3857"]]
Origin = (2926801.391389018390328,6009567.154649123549461)
Pixel Size = (2411.811113223014218,-2411.811113223014218)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 2926801.391, 6009567.155) ( 26d17'30.86"E, 47d24'42.83"N)
Lower Left  ( 2926801.391, 4603481.276) ( 26d17'30.86"E, 38d10'12.31"N)
Upper Right ( 6920760.595, 6009567.155) ( 62d10'12.90"E, 47d24'42.83"N)
Lower Right ( 6920760.595, 4603481.276) ( 62d10'12.90"E, 38d10'12.31"N)
Center      ( 4923780.993, 5306524.215) ( 44d13'51.88"E, 42d57'51.12"N)
Band 1 Block=1656x4 Type=Byte, ColorInterp=Gray

So the question is:

How to re-project wide-angle camera? Looks like gdal incorrectly process it. I.e. pixel size in the middle of the picture is ~1km. While gdalinfo makes it ~2km. Which might be ok on the edges.

Another question:

Why only "-tps" can handle non-linear distortions on the edges?


texts: Calibrating Landsat Level-1 Precision and Terrain (L1TP) corrected data?

texts:
My study will involve Landsat OLI, ETM+ and TM scenes for forest biomass change assessment (Time Series Analysis) using different vegetation indices and textures. According to the official site https://landsat.usgs.gov/landsat-collections " Landsat scenes with the highest available data quality are placed into Tier 1 and are considered suitable for time-series analysis. Tier 1 includes Level-1 Precision and Terrain (L1TP) corrected data that have well-characterized radiometry and are inter-calibrated across the different Landsat instruments".

Keeping the above in view I take that I can use them straightforwardly assuming that all three OLI, ETM+ and TM stand at the same line and can be used interchangeably since to me it does not make any sense to do any further calibration when these Surface Reflectance Products have already been pre-processed and made ready for respective use and carry well characterized radiometry across the different Landsat instruments.
Kindly make me clear please if I stand correct or still I would need any further correction or adjustment of my data.


texts: Convert Polygon Feature Centroid to Points using Python

texts: I would like to convert some polygon based shp files that have several polygon features into points for each feature that would essentially represent the centriod of each polygon feature.  I know in ArcGIS world I could use Feature To Point tool but I would like to keep this in a script that could be run on PC's that don't have arcpy on them so I'm looking for an open source alternative to that.  Is anyone aware of a library I could use for this along with some direction on how to leverage it accomplish this?


texts: Examples and physical materials that are helpful for teaching GIS principles

texts:
I am collecting some of the resources and materials that are helpful when teaching the basics of GIS to extend my current toolbox. Some of those I use myself when teaching introductory courses. I am looking for other hands-on materials (preferably, non-digital) that can help understand some of the GIS/geography related concepts such as projections, spatial analysis, and others.

Probably one of the most famous examples is a trick with an orange peel when a navel orange is being peeled to demonstrate the projection from a globe to a flat map. National Geographic has detailed instructions on this excellent tutorial along with other materials for constructing 3D models of Earth from paper (read more about Investigating Map Projections).

I also use an overhead projector for demonstrating with transparent plastic sheets how a route can be drawn (with a pen) while avoiding forbidden territories over multiple layers (conceptually, path analysis and geographic overlay illustrations). Great guidelines are available in the book The Really Practical Guide to Primary Geography.

What is great about these examples is that students can try this themselves - to see the holes in the orange peel and to see how hard it is to find a route taking into consideration the areas of multiple sheet layers.

What are other known hands-on materials that can help illustrate spatial/GIS concepts? Things that are related to computer science (such as database management) that have something related to geography would also fit.


texts: Writing custom data-source or plugin to load Google My Maps data into QGIS

texts:
My goal is quite straightforward: display maps created and mantained in Google My Maps as a layer in QGIS.

The way I thought about doing this is the following:

Write a Python script that acts as some sort of server;
Get this server script to fetch data from Google My Maps (as KML) and serve this data somehow to QGIS, most probably using WFS API or similar;
Add the server address (something like http://localhost:1234 I guess) as a readonly layer in QGIS.

Then, I could edit my maps using Google My Maps, but would be able to keep displaying their current state in QGIS.

So my questions are:

Does this scenario make sense? Is this a viable thing to do?
Is this "webserver" approach the best one? Should I be aiming at writing some sort of plugin instead?
Is it something that has already been done before?


texts: Debugger for ArcGIS Java Script API 3.x?

texts:
I am working on ArcGIS Java Script API version 3.2/3.3 and using IDE as aptana & notepad++(Currently testing all applications on Firefox and using Firebug while testing). Getting several errors while using JSAPI so is there any recommended Debugger available for ArcGIS JSAPI ?

I am also going through resource link.


texts: Where can I find Specifications for Canadian National Topographic Map series?

texts:
At one time the Natural Resources Canada had the specifications for the National Topographic Series (NTS) maps online at http://www.nrcan.gc.ca/earth-sciences/cartospecs. It included things like the exact CMYK colour numbers, size and width of every symbol in points and so on.

Web.archive.org has a record of it here (be patient, time travel is slow), but it only has the text. The image examples are missing. Before that I remember having a pdf document of the same, but I can't locate it now. I do have a hardcopy from 1988, but that is hard to search and a little long in the tooth...

Does anyone know of another place for this useful resource?

(And yes, I've already asked  NRCAN, but I've yet to get an answer other than "someone will get back to you". I'm sure I will eventually get some kind of answer, but I'm getting impatient P:-)


texts:
I have a raster file with classifcation codes, and want to reclassify valid codes to 1, and all others to 'NoData'.

RMAP2 = arcpy.sa.RemapValue([[3.2, 1],
                             [4.1, 1],
                             [4.2, 1]])

LCA = arcpy.sa.Reclassify(LCA, "LCCODE_USE", RMAP2, "NODATA")

The output however only contains the values 0 and 127, with the 0 value coinciding with just one of the missing values.

Any idea what the hell is going on?


texts: How to add customization to Configurable Map Viewer (CMV)?

texts:
Looking to learn more about CMV Viewer. For example, how would I add capability like:

map on load:

for learning how to configure CMV Viewer for my own needs, I'd like to start by learning to add "click" for "show coordinates"...

here's the example of the code I'd like to add onto my copy of CMV Viewer. Can someone teach me how to add on this type of customization?

Here's what I'd like to add:
http://developers.arcgis.com/javascript/sandbox/sandbox.html?sample=map_infowindow


texts:
I recently got an script tool in a tbx file to play with.

I'd like to see how the tool runs in a debugger to get a better idea of how it works. I've tried several paths to this, but all of them seem to have failed:

Open in visual studio, or some other IDE: the windows explorer, of course, does not recognize .tbx as a container, and just displays garbage.
Right click the tool and use the "Debug" option: tool runs as usual, but I have no control over stepping, breakpoints, etc. I can't see the code as it is running.
Right click the tool and use the "Edit" button: python script opens in NotePad. There doesn't seem to be a way to choose another editor. Of course, notepad doesn't have a build or compile option, so I can't modify the code even to print output.
Copying the code into Visual Studio and running it there: when trying to run debug from ArcMap, I run into this issue


texts: is ellipse inside ellipse in WGS84 ellipsoid

texts:
My question is quite clear:

Assume that there is an ellipse, namely EL1. Parameters of the ellipsoid are (focus11 [latitude, longitude], focus12[latitude, longitude], radius1).

And a second ellipse, namely EL2. Which is defined with the same parameter set (focus21 [latitude, longitude], focus22[latitude, longitude], radius2).

The reference ellipsoid is WGS84.

What is the method to find if EL1 contains EL2 ?

So far the best approach I come up with is to convert the EL1 and EL2 to polygons. But it is quite costly.

Is there any idea?
Any help is appriciated.
Waiting for your replys ...


texts: List items are replaced with new reference when adding in ArcObjects?

texts:
Based on the query i am getting some features in pFeatureCursor. I am adding that features into a List(IFeature). When iam iterating the pFeatureCursor and adding the values into pFeatureList all values in the list are replacing with the new reference pFeature

Please  refer my code below

//List is declared globally
List<IFeature> pFeatureList = new List<IFeature>();

//In a method i am trying to add features in a list
//based on queryfilter
pFeatureCursor = pFeatureClass.Search(pQueryFilter, true);
IFeature pFeature = pFeatureCursor.NextFeature();
while (pFeature != null)
{
  pFeatureList.Add(pFeature);
  //lstQuerValue is the listbox
  lstQueryValue.Items.Add(pFeature.get_Value(1).ToString());

  //When below line is execute whatever the recent pFeature
  //value is replacing all the items and duplicating in pFeatureList
  pFeature = pFeatureCursor.NextFeature();
 }

Please give some idea how to avoid it.


texts: Publishing MapService to ArcGIS Server problem (Slow and cannot publish data from SQLDatabase)

texts:
I am trying to publish a map to ArcGIS Server with one layer from SQL Server.

We have:

1. Data:

SQL Database with data
Most of Feature Classes are stored in Feature Sets

2. SQL Server (sqlsrv) with Data

3. Machine with ArcGIS Server agserver

This is a virtual machine
ArcGIS Server 10.3 installed
ArcGIS Server is running using domain account (agsacc)
agsacc account is an Administrator on agserver machine (Can access any folder)
Database (sqlsrv) is registered and verified (using Operational System Authentification)
I can see and use data from SQL sqlsrv in ArcMap on 'arcgisserver' machine

Problems:

When I am trying to publish a shapefile (150kb)

Service published successfully, but takes 1-2 hrs to publish this service.

When I am trying to publish anything from sqlsrv, it takes also about 2 hrs, but:

Publishing process failed with the error:

ERROR 001369: Failed to create the service.
Failed to execute (Publish Service Definition).
Failed at Thu Mar 24 18:46:18 2016 (Elapsed Time: 5 minutes 23 seconds)
Failed.
ERROR: the server's publishing job didn't succeed.
Failed to execute (UploadServiceDefinition).

I can see in ArcGIS Server Manager that my service is "stopping"
When I am trying to stop it and run again, I have the following error:

 Failed to construct instance of service 'test4.MapServer'. Can not open file C:\arcgisserver\directories\arcgissystem\arcgisinput\test4.MapServer\extracted\v101\test4.msd. The system cannot find the path specified. Probable cause: The file is inaccessible to Server.

There is nothing (0 files) in this service folder

What I have tried:

-Set log level to Verbose (couldn't find anything). I am not sure what should I be looking for.

-Repaired ArcGIS Server Installation

-Reinstalled ArcGIS Server (still doesn't work)

Update, 160328 1500

Just turned off Antivirus and shapefile published for a couple of seconds.
However, map with data from sqlserver

Got the following errors in the Log:

-Core server call to create service failed.

-Error executing tool.: ERROR 001369: Failed to create the service. Failed to execute (Publish Service Definition).

-Failed to create the service.: ERROR: Service failed to start, Internal server error.

-Failed to construct instance of service 'test20160328_2.MapServer'. Service startup timeout

Is there is a way to understand where is the problem?


texts:
This is the first time im using python/arcpy. I keep getting an error with import arcpy command.

Can anyone help me find out why?

import arcpy

  Runtime error  Traceback (most recent call last):   File "",
  line 1, in    File "c:\program files
  (x86)\arcgis\desktop10.3\arcpy\arcpy__init__.py", line 24, in

      from arcpy.toolbox import *   File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy\toolbox.py", line 359, in

      from .management import Graph, GraphTemplate   File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy\management.py", line 26, in

      import _management   File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy_management.py", line 14, in

      import _graph   File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy_graph.py", line 27, in

import numpy

  ImportError: No module named numpy


texts: Creating map of New Mexico where deaf people live in county from US Census database?

texts:
I would like to get some feedback from those who have been doing or perform on US Census Map.

Is there a way for me to get information on Deaf person living in New Mexico. Would US Census be helpful?

On the other hand, do I need to find more information from the state of New Mexico?

I wanted to show a map of where each county where a deaf person lives? Is that possible or is it the word "deaf" too broad to be included in the US Census database?


texts: QGIS to ArcMap file delivery via geopackage

texts:
When I'm saving a new layer made in QGIS (v 3.14), it gives me only an option of saving it as .gpkg file. I sent these to a person using ArcMap for further work and he's insisting that I send him all (.shp, .dbf, .shx, .prj, .sbn, .sbx). I have used ArcMap before, but within last year I have access only to QGIS.
Is it possible that ArcMap cannot open .gpkg file or there is something else that is the problem?


texts: Sharing labels and have multiple leader lines for like features in ArcMap?

texts:
I have many labels in ArcMap and want to use one label with multiple leaders to call out features whose labels would end up being the same. Is this possible?

As it is right now it is not possible for me to fit all labels without being able to do this.

Using arcmap 10.3


texts: Gdalinfo and QuantumGIS does not recognize Spatial Reference even though Python does

texts:
I apply the following Spatial Reference to a file:

spatialRef = osr.SpatialReference()
spatialRef.ImportFromEPSG(3411)
outraster.SetProjection(spatialRef.ExportToWkt() )

geotransform = (-3850000.0, 25000.0 ,0.0 ,5850000.0, 0.0, -25000.0)
outraster.SetGeoTransform(geotransform)

If I use "icechart.GetGeoTransform()" and "icechart.GetProjection()" in Python, the output looks correct.

>>>icechart = gdal.Open(r'C:\Users\max\Desktop\test\icechart_persistencemap.tif')
>>>icechart.GetProjection()
'PROJCS["NSIDC_Sea_Ice_Polar_Stereographic_North",GEOGCS["GCS_Unspecified datum based upon
the Hughes 1980 ellipsoid",DATUM["D_",SPHEROID["Hughes_1980",6378273,298.279411123064]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Stereographic_North_P
ole"],PARAMETER["standard_parallel_1",70],PARAMETER["central_meridian",-45],PARAMETER
["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]'

>>>icechart.GetGeoTransform()
(-3850000.0, 25000.0, 0.0, 5850000.0, 0.0, -25000.0)

However, gdalinfo tells me that there is no Spatial Reference defined, neither does QuantumGIS.

When I get the Projection and Geotransform from another file, rather than defining it manually, gdalinfo finds the Spatial Reference:

outraster.SetGeoTransform(icechart.GetGeoTransform())
outraster.SetProjection(icechart.GetProjection())

Also when I reproject the raster above with gdalwarp, gdalinfo and QGIS recognize the Spatial Reference:

os.system('gdalwarp -s_srs EPSG:3411 -tr 25000 -25000 -t_srs EPSG:3575 -of GTiff '
+ infile + ' ' + outfile)

So what is wrong with the first lines of code in this post? Any advice?


texts: How to populate my table with data in excel format?

texts: I have 160000 rows of data in Excel and I want to put them in my table in PostgreSQL. How can I import X,Y coordinates from my excel file to my table in postgresql?


texts: How to use GDAL to import lat/lon points from ESRI File Geodatabas into MySQL

texts:
I have a file geodatabase from here:
ftp://ftp.csc.noaa.gov/temp/MarineCadastre/AIS.SampleData.zip

Using ogrinfo, I can see that it contains a layer called "Broadcast_" that has lat/lon points:

C:\Users\Dirk Beer\Desktop\AIS.SampleData>ogrinfo Miami.gdb Broadcast_
INFO: Open of 'Miami.gdb'
      using driver 'FileGDB' successful.

Layer name: Broadcast_
Geometry: Point
Feature Count: 1365578
Extent: (-80.282093, 25.591130) - (-79.751260, 25.950708)
Layer SRS WKT:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4326"]]
FID Column = OBJECTID
Geometry Column = Shape
SOG: Integer (0.0)
COG: Integer (0.0)
Heading: Integer (0.0)
ROT: Integer (0.0)
BaseDateTime: DateTime (0.0)
Status: Integer (0.0)
VoyageID: Integer (0.0)
MMSI: Integer (0.0)
ReceiverType: String (0.0)
ReceiverID: String (0.0)
OGRFeature(Broadcast_):1
  SOG (Integer) = 0
  COG (Integer) = 77
  Heading (Integer) = 511
  ROT (Integer) = 128
  BaseDateTime (DateTime) = 2008/12/31 23:59:00
  Status (Integer) = 0
  VoyageID (Integer) = 11
  MMSI (Integer) = 366883280
  ReceiverType (String) = b
  ReceiverID (String) = 003669972
  POINT (-80.161767 25.773433)

I would like to import that table, with the points into a MySQL database. However, when I try, the points column does not show up in the database (the column "SHAPE" contains binary data, not lat/lon points):

ogr2ogr -f MySQL MySQL:ais2,user=ais,password=ais Miami.gdb -lco ENGINE=MyISAM

show columns in broadcast_;
OGR_FID
SHAPE
sog
cog
heading
rot
basedatetime
status
voyageid
mmsi
receivertype
receiverid

How can I import the lat/lon points?

Thanks for your help.


texts: I am trying to create a KML from a layer composed of a large number of polygons.  Because of the large size, it is quite slow to load the KML in Google Earth.  Is there any way to create the KML so that it is only the polygons that are in the display area are visible?


texts: QGIS calculate corners and create polygon

texts:
I would like to create a polygon in QGIS. I have photos with four points (X) I can see, which are not supposed to be the corners of the polygon, but are needed to measure these corners (Y).

The setup looks like this:

X   Y     Y   X

X   Y     Y   X

What I know is where the four points X are, and that when I connect the two upper X, on this line, at 1/6 and 5/6 are the upper points Y (the same applies for the lower points). I would like to load the photo into QGIS, tell it where the four points X are, and get a polygon with four corners at points Y.

I have hundreds of these photos, so a semi-automatic or quick way would help a lot.


texts: Attribute Tables with drop down menus from lookup tables with joins

texts:
I am attempting to create drop down lists when editing attributes in Arcgis Desktop in order to eliminate human errors. Through using domains I managed to carry this out successfully and the drop down lists work well.
However I have an application that caters for lookup tables and not domains, and in order to avoid having to modify both the lookup table and the domain I want to create drop down lists for attributes on Desktop using lookup tables with joins.
Is there a way of doing this as I can't seem to carry this out successfully?


texts: Converting clustered hexagonal grid units into polygons

texts:
I created a hexagonal grid using MMQGIS for my project and after some selection, I ended up with several clusters of units (as can be seen in the picture). I define as a cluster all the small hexagonal units that are connected.

What I wanted now is to convert these clusters into polygons so I can calculate the area of each polygon after. How can I do so?

QGIS version: 2.18.17


texts:
I have ~7,000 2 line dbf files (zonal stats output) that I would like merged into one table, all the same format. I have used the arcpy merge tool, as seen below.

    #merge tables together

    outfile = "C:\\Pivot\\hourly.dbf"
    arcpy.env.workspace = "C:\\Pivot\\dbfs"
    listTable = arcpy.ListTables()
    print("...")
    arcpy.Merge_management(listTable, outfile)
    print("All done!")

and it works but INCREDIBLY slowly. Adding about one file every 10 minutes. Is there a way to concatenate these files using a dbf module? Or is there a reason this is functioning so slowly? I've already processed all of the dbfs a couple times with other tools and they weren't nearly this slowly.


texts:
We are making a simple webmap with a Postgis database and Geoserver. So now we have the layers working in Geoserver and we can see the cartography and the the features using GeoExt, but we're not able to see the database related to them. Our commitment is to have a frame at botton of the page, or a popup, where the attribute table could be seen once a single feature is selected o maybe the entire layer if we select it.

Any idea on how to do this?


texts: Elevation interpolation along a (sewer) line

texts:
Here's what I want to do.  I have a vector sewer line feature class that contaings a field for the high invert elevation and one for the low invert elevation.  I can turn them into a raster no problem but what I can't figure out is how to interpolate between the two ends (high and low) and assign the elevations to the raster cells.

Ultimately what I need to do is compare the sewer line raster to an elevation raster to determine likely overflow points in the event of a backed up line.  My thought is it would be similar to creating a reservoir with a blockage being the dam.

I hope I articulated my question sufficently.


texts: Smoothing raster map using ArcGIS Desktop?

texts:
I have a raster map of US Midwest which is very sparse, i.e. the pixels of interest are few enough to be almost invisible when viewed at a scale where all states of US Midwest are visible. I would like to follow the approach outlined in this PNAS paper (http://www.pnas.org/content/110/10/4134.full) to create a better map, but not sure how to replicate it in ArcGIS Desktop.

The PNAS paper outlines the steps as follows:

  Because of the small sizes and scattered distribution of change areas,
  it was difficult to visualize regional patterns of LCLUC at the
  original 56-m spatial resolution. As a result, we used spatial
  smoothing techniques to create a regional change surface that
  highlighted local hotspots of change. Related approaches are used in
  fields such as spatial epidemiology to generate stable estimate of
  disease rates (48) but have not been broadly applied in the field of
  land change science. In our smoothing approach, change pixels at 56-m
  spatial resolution were first aggregated to the percentage of change
  at 560-m resolution. This was done by taking 10-by-10 blocks of 56-m
  pixels (i.e., 100 pixel blocks) and summing the binary change within
  each block (Fig. S4A). Next we used a 2D kernel smoother to compute a
  smoothed estimate of percent change for each of the 560-m resolution
  pixels (Fig. S4B). A quartic kernel function was used to calculate
  moving averages across the study area at a bandwidth of 10 km. The
  same quartic kernel function was used to smooth percent change from
  corn/soy in 2006 to grassland in 2011. Finally, we generated a
  smoothed map of grassland cover in 2006 by aggregating grassland
  presence at 56-m resolution to percent grassland cover at 560-m
  resolution, and then smoothing this aggregated cover layer by using
  the same 10-km quartic kernel. This smoothed grassland cover layer was
  subsequently used as the denominator in generating a map of relative
  rates of grassland conversion.

As far as I understand, this is the flowchart:

Use block statistics in ArcGIS to sum 10x10 pixels of 56-m raster to 560m raster
2D kernel smoother: not sure how to do this
Quartic kernel:  not sure how to do this

Not sure how to progress beyond step 1


texts: Identifying Parallel line segments at fixed distance in same shapefile layer using ArcGIS Desktop?

texts:
I have a single, polyline shapefile for a utility grid.  This shapefile displays overbuild lines by offsetting them from each other approximately 40 feet, and running the lines parallel to their true location.

I am trying to identify those areas where I have parallel lines within 40 feet, likely meaning they occupy the same corridor.

How do I do this using geoprocessing functions, as opposed to Python scripts?


texts: Can't fill the attribute table programmatically

texts:
I've got a shp-file, here is it's attribute table:

It's easy to fill the attribute table of some pFeature programmatically if I get the feature from cursor:

For example, "RECNO" Field(3) has a numeric format

Set pFeature = pFeatureCursor.NextFeature
Set pFeature.Value(3) =2
pFeature.store

Usually, it works fine, I shouldn't even start the Edit session for this.

But when I add a new feature and try to fill the attribute table in this code below, the error "Object required" appears.

Dim pMxDoc As IMxDocument
Dim pFeatureLayer As IFeatureLayer
Dim pFeature As IFeature
Dim pFeatureClass1 As IFeatureClass
Set pMxDoc = ThisDocument
Set pFeatureLayer = pMxDoc.FocusMap.Layer(0)
Set pFeatureClass1 = pFeatureLayer.FeatureClass

Set pFeature = pFeatureClass1.CreateFeature

' HERE IS THE PROBLEM
Set pFeature.Value(3) = 2


texts: Converting unknown type of coordinates to lat/lon

texts:
Our company is using Total Station (Surveying Machine using SOP) to calculate GPS coordinates and its results are something like this

Easting: 3250193.00, Northing: 1115573.00, Elevation: 656.00
Easting: 3249214.00, Northing: 1111220.00, Elevation: 735.00
...

I want to convert these points to decimal based lat/lon (so that can map those points to google maps). Can someone help me out in converting these points?

The above points are somewhere near to Muzaffarabad, Azad Kashmir. (Some where near to 34.3629, 73.4775)


texts: Creating pop-up box if condition is true in table using ArcPy with ArcMap

texts:
I'm creating a script that will grab layout data from a table for the current Building Number of the current building my map series is on. It works fine, but I need it to pop up a window if multiple entries exist in the table for that same Building Number. In this window I would need to be able to choose one of the entries identified by their BLDNAME attribute.
Is this possible with ArcPy?


texts: Openlayers 3 Export map with highlights and popups

texts:
I try to export map in openlayers 3 with this example (Export Map Example) and it works. But when I export the map I get this picture :

However I want to get this picture with highlights and popups the way I see them in the browser :

Is there any way to do this?