QGIS How to multiply (repeat n times) a point by field


I have a point layer in QGIS that has a field named “Amount” that stores an integer value for every point in the layer.

I want to create N points in the exact place of the feature that has N as value in the amount “field”. I mean, if a point has in the amount field 16, I need to create 15 more points in that exact lat long (1+15 so that there are 16 in total) of that point.

After this process, the sum of “Amount” values in the starting layer should be equal to the sum of points in the new layer.

I need it to create a heatmap with the new layer. What is the easyest way of achiving this?

Regards

Exporting feature class into multiple feature classes based on field values?


I have a feature class with over 2,000 features, and I need to make them all individual feature classes based on a field.

I know there has to be a way to do this, but I just can’t figure it out!

Any takers?

Mapinfo : How to cut a raster with a vector layer?


With mapinfo 12, i would like cut my raster layer .tab with one of my vector layer.
How to do this basic task with mapinfo?
Is it possible with mapinfo to convert .tab raster table toward any standard formats (ecw, tif, etc)

Thank you.

Is there a geoprocessing tool that I can use to find the closest road to a point on a map?


I have a Network Dataset layer that I use to perform routing. Even if a starting or ending point is not near a road in the Network Dataset layer, I still want to perform routing using those points.
Is there a way that I can get the routing to start at the nearest place on the ND, to the starting point, or end on the route at the closest place near the ending point?

ArcGIS 10.1 – Error No spatial reference exists


I’m encountering a new error on some geoprocessing. It originally showed up while running commands in arcpy but now it crops up whenever I run some tools. I’ve been using the Extract by Mask on some rasters, and then using the Focal Statistics function on the resulting clipped grid. The main error that’s showing up is:

ERROR 999999: Error executing function
Unexpected operation
No spatial reference exists

A screen grab of the error messages is attached. The other two errors (010302 and 010067 ) have no helpful information. I’ve been running similar processes for months and haven’t encountered this before. Both the input layers (a raster and polygon feature) are in WGS84 with their coordinate system defined. Both of these input layers are in a File GDB and I did export them to a shapefile and GRID and the process executed successfully. I could run this series of processes outside of the GDB but it’d be nice not to.

Any thoughts on why this is happening?error message

Using valueAsText, giving me unicode data instead of a list


I am learning Python to create toolbox right now. When I’m trying to get the parameter, using valueAsText, the date returned is unicode instead of a list that I was expecting.

This is the code that I was using:

inFeatures = parameters[0].valueAsText

I was hoping the input data would be stroed in a list in this way, put when I use for loop, this is what I have:

    messages.addMessage(type(inFeatures))
    for inFeature in inFeatures:
        messages.addMessage(inFeature)

enter image description here

It is dividing the path into single characters.

Any way I could get the result in a list?

Total area of ArcGIS Intersect output is NOT equal to the smaller input?


I have two polygon inputs that share the same projected coordinate system – one is land cover data, the second is parcel data (of smaller total area). When I execute a standard geoprocessing Intersect, and then recalculate the geometry of the Output’s polygons, the resulting total area of the Output polygons is not equal to the total area of parcel data total area (Output is ~200 acres < the parcel data input). I might anticipate slivers in the original land cover data, but the Output shows that all area is classifiable (no null Land Cover classes). Is this a normal occurrence? Or is there something I can/should do to get these areas to match exactly?

Also, not sure if order matters, but I listed the Land Cover data first, and the Parcel Data second, when adding Intersect inputs.

I’m using ArcGIS 10.x for Desktop.

How to obtain data to create average precipitation in 2012?


I am currently doing a final project, which I have to create a map of 2012 average precipitation of lower 48 states by using natural neighbour interpolation. The problem is that how can I get the data? I already have the data of 2012 from NOAA website in .CSV format. There are more than 80000 geographical locations in the file (http://www1.ncdc.noaa.gov/pub/orders/cdo/434898.csv) Would it be possible for me to use them to create the map and if anyone has any idea for
me about getting a data please help.

Clipping MrSID files in ArcGIS Desktop?


Does anyone have experience subsetting/clipping a large NAIP imagery file in MrSID format? Some of the NAIP files we have are by county and they have really large file sizes.

The files are from the FSA NAIP Imagery.

How does it work if you subset a MrSID file? Does ArcGIS has to uncompress it first before it can proceed to clip or subset an area of interest?

ArcGIS Linear referencing tools duplicating records and erroneous measures


I am using LocateFeaturesAlongRoutes_lr, MakeRouteEventLayer_lr, and MakeRouteEventLayer_lr to combine three different overlays into one based on this question I asked month ago. The code is as follow:

try:    
    print "ntttLocating features along roads",
    arcpy.LocateFeaturesAlongRoutes_lr(FTR_TO_RS_TEMP, ORIG_RD_FTR,
                                       "CONTROL_SE", "0 Meters",
                                       FTR_TO_RS_DB_TEMP,
                                       "RCSKEY LINE FMEAS TMEAS", "FIRST",
                                       "DISTANCE", "ZERO", "FIELDS",
                                       "M_DIRECTON")

    print "ttt[OK]nttcombining features and PDT schedule",
    arcpy.OverlayRouteEvents_lr(FTR_TO_RS_DB_TEMP,
                                "RCSKEY LINE FMEAS TMEAS",
                                LP_SCHEDULE, "RCSKEY LINE START END",
                                "UNION", ROAD_LP_JOINED_TEMP,
                                "RCSKEY LINE FMEAS2 TMEAS2",
                                "NO_ZERO", "FIELDS", "INDEX")

    print "ttt[OK]ntttNew event route layer",
    arcpy.MakeRouteEventLayer_lr(ORIG_RD_FTR, "CONTROL_SE",
                                 ROAD_LP_JOINED_TEMP,
                                 "RCSKEY LINE FMEAS2 TMEAS2",
                                 LP_EVENT_LAYER_TEMP)

except Exception, e:
    tb = sys.exc_info()[2]
    print "Line %i" % tb.tb_lineno
    print e.message

For some reason, the routine is causing the following two issues:

  1. Duplicate records of overlayed route events. For example, when ArcGIS creates a new record based on a union of two others, it would duplicate it 5 to 10 times.
  2. New fmeas and tmeas that are -0.000023 to 0 and 0 to 2.999977 instead of 0 to 3 as in the original layers

I am quiet certain that the two issues are related and was wondering if anyone can shed a light

Question and Answer is proudly powered by WordPress.
Theme "The Fundamentals of Graphic Design" by Arjuna
Icons by FamFamFam