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

InfoWindows in Nested Polygons (Google Maps Api 3)


I have a large polygon with thick borders inside which there are several smaller polygons something like multiple city polygons inside a state border. The problem with this is that the info-windows for some of the polygons don’t show up when i click the polygon.(i have set up a click listener for each smaller polygon)

Toolvalidator script is not working in ArcGIS Server


I have made a small GP tool in ArcMap, which is also having some validation code in tool validator class. And it is working fine as GP tool in ArcGIS desktop(Arcmap). But when i published the same tool as Geo-Processing service added as toolbox in ArcMap, it is not considering any validations.
Ex: Disable/Enable some textboxes.

What should i do to integrate these validations in GP service also?

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