September 10, 2013, 12:31 am
Hi,
I have a python geoprocessing that takes a lot of time. The geoprocessing calculates geometries using points, polygons, buffers, and cursors.
Do you have any hint of how to optimize the geoprocessing?
Would it be shorter in time if I used the Java API?
Thanks.
↧
September 10, 2013, 3:33 am
Hi
I have 10 geoprocessing services on my server. But I have problem with these GP Services. When I run GP Services with large area the GP service dont work. Is there any limitation on ArcGIS Server by default?
Regards
A.Ponker
↧
↧
September 10, 2013, 9:21 am
Having issues with the titles (X-Axis) and (Y-Axis) in my graphs not working properly. I have it set in a Python script to make it "KWH_USAGE" but when the graph is created - it shows "kwh_usage"... Why is it only showing everything as lower case?
Also - I want to include the data from a specific field in a feature class to show in the title. So if I have 14884512 in a field called "TWACs_Number" - I want the title to include that. How would I write that?
Thank you
↧
September 10, 2013, 10:11 am
I am running a geoprocessing tool in Modelbuilder that saves a temporary feature class which is normally added into the table of contents at the completion of the tool. I then review it for completeness and append it to a permanent feature class. Sometimes, even though the tool works properly, the temporary feature class is not added into the TOC. Is there some way to connect to and add data from the In_memory location?
↧
September 10, 2013, 2:40 pm
Hi,
I have a large raster file, and I am trying to get the weighted average values of a number of pixels located either entirely or partially inside of a shapefile boundary, but I'm not sure how to do this. I was thinking something along the lines of cutting the raster file based on the shapefile boundaries, but I'm not really sure. Any advice would be appreciated.
Thanks,
Scott
↧
↧
September 10, 2013, 2:59 pm
Hi,
If you generate the Semivariogram using Geostatistical wizard. it will show you the Semivariogram map at the bottom of the window. (which looks like a raster map) (following image)
Attachment 27332
when you move the mouse pointer on the map, two numbers will pop up for each pixel
one is the Value of the cell and the other one is the Weight.
I understand that the "Value" should be the average of moment of inertia of the pairs in that specific bin
But I don't know what is the "Weight"? I could not find anything regarding to this number in the manual.
I would appreciate if someone could shed some light on this for me.
Regards
E
↧
September 10, 2013, 8:29 pm
I have a tool in a geoprocessing service that executes from Toolbox, but not from the REST endpoint on 10.1. In AGS 10.0, it executes fine from either. The problem is "No module named PParams". This is just a module in the same folder as the Tool (PParams.py and PParams.pyc).
the code to import it is just
Code:
from PParams import pparams
I did a "print sys.path" to check and that folder is in sys.path.
Any ideas why it cannot find the module when run through the REST endpoint? At the moment, my guesses I am looking into are
1) That perhaps when 10.1 tries to package it up (have not been able to find any docs on exactly what it does), that it is leaving out some items, though that is only conjecture.
2) Some change in how the environment is handled in 10.1 that it could not find another .py file in the same folder as the tool script.
THANKS!!!!!!
↧
September 11, 2013, 1:58 am
Hello,
Here is a follow-up from another question:
http://forums.arcgis.com/threads/918...es-with-Python
As a workaround, I first created a python list with strings pointing to individual raster bands ("raster1.tif/Band_1", "raster1.tif/Band_2", "raster1.tif/Band_3",etc).
Then I used it in this small script:
Code:
def raster_overlay(input_files):
for raster in input_files:
ExtractMultiValuesToPoints(in_shp, raster, "BILINEAR")
in_shp = "temp/shapefiles/points.shp"
raster_overlay(raster_files)
The shapefile contains around 50 points, and each raster is around 15MB.
If the list of rasters is small (more or less <20), everything works as expected. However, I have around 50 raster images. Using the complete list, Python crashes ("python.exe has stopped working") after processing 21 images or so.
The values of the first 20 rasters are correctly appended to the shapefile (one column for each raster).
I monitored python.exe in the Windows 7 Resource Monitor. Memory usage of Python starts around 15 MB and then gradually increases to around 700 MB (and then Python crashes). So, memory leakage might be the reason for the crashes (although the PC has 8 GB of RAM).
Is this a known bug?
Is there a workaround?
By the way, I posted the same question earlier here:
http://gis.stackexchange.com/questio...crashes-python
↧
September 11, 2013, 6:47 am
I have a geoprocessing serviced that takes two parameters: netCDF (file) and Date (str). I'm working off of the demo found
here.
The issue I have is getting the netCDF file to the geoprocessing service. By creating an upload feature I can see how I would be able to provide the needed parameter but I already have all the necessary files server side.
How can I get my web application to read local netCDF data into the geoprocessing service?
Any suggestions would be greatly appreciated.
Thanks!
↧
↧
September 11, 2013, 9:25 am
Hi all,
I have WRF climate model output data that I am trying to rasterize in Arc. The WRF cells are 3km x 3km in this scenario, and WRF gives me the latitude and longitude of the centroid of each cell. Upon trying to rasterize this point data set, I cannot manage to create the same grid WRF is using because it appears ArcGIS auto populates the cell with no manual override.
This grid also has to be not running North-South by angled at about 330 degrees. This second capability I have not been able to find in Arc. I have tried using Point to Raster as well as creating my own raster with the known bounds of the grid, but I run into one or bothe of the problems described. You can see what I am talking about in the image below.
Attachment 27354
Has anyone here created a raster with an angled (non North-South) orientation? How is it done?
I am using ArcInfo 10.0.3 build 2414, and have all possible extension rights.
Thanks for any help!
Ona Strikas
PhD Student
Dept of Geology
UNC- Chapel Hill
↧
September 11, 2013, 11:06 am
I know this is a very elementary question so please excuse what might sound so simple. I am trying to calculate the total area of a raster dataset that has a specified Cell Size of 1 m. I used IDW to interpolate carbon lbs/per area and I’m confused as to how to understand the Raster properties. I can see how many Columns and Rows there are and I know that the inputs are in UTM and I also see the No Data values but I’m unclear as to how to understand all of this. I’ve read and re-read the Help articles and remain just as confused. I see the Extents but I think these numbers correspond the geographic location of the data, not the size of the area. I also used the Zonal Stats tool to Sum the amount of Carbon in the study area and the table shows Count, Area and Sum but I don’t know how to interpret these.
Any help is greatly appreciated.
Thankyou!
Natasha
Attachment 27364
↧
September 12, 2013, 5:39 am
I am trying to make a model that essentially duplicates the Distributed Geodatabase/Extract Data Wizard toolbar that you can use in ArcMap. The Extract Data wizard works fine in ArcMap with my GeoData service, but I want to make it easier for users to extract the data by using a model. I cannot find this specific tool in the Toolbox. Anyone have any ideas?
Thanks
↧
September 12, 2013, 8:23 am
Hi all,
I wonder if anyone have some ideas about how to solve this geoprocessing situation:
1: I have a raster layer called mixlyr_apr (representing depth values)
2: I have a series of raster layers called: vv_10_apr; vv_20_apr; vv_30_apr; vv_40.2_apr; vv_55.5_apr; vv_78.5_apr; vv_113_apr and vv_242.6_apr (representing vertical velocity for a particular depth)
All rasters have the same extent, cell size, etc.
Say, for example, the cell value of mixlyr_apr is 60. In that case I need to produce a raster layer with the average values (per cell) of the vertical velocity layers that represent the depths below the value of the mixlyr (60 in this case). That is to say, the average values (in this case) of vv_10_apr; vv_20_apr; vv_30_apr; vv_40.2_apr and vv_55.5_apr
Hope doesn't sound too confusing and some of you can suggest some ideas.
Many thanks
↧
↧
September 12, 2013, 3:51 pm
Hi,
I have a number of rasters showing precipitation values for different months. I'd like to create a time series to visually show how this changes over time. I've tried googling for answers on how to do this, but have come up short. Any advice?
Thanks,
Scott
↧
September 13, 2013, 7:11 am
Thanks in advance.
The title says it all. I need to find the total lengths in each cardinal direction (N, NE, E, SE, S, SW, W, NW) for all rivers per county in multiple counties.
I had thought about using flow direction (spatial analyst) but this does not provide what I need from what I can tell.
Thanks for any and all help!
↧
September 13, 2013, 8:22 am
When I attempt to tie the output of "XY to Line" to the next step, "Add Field", model builder simply stops after "XY to Line", thinking it's all done. It does not execute subsequent steps. Is there something wrong with doing this? It seems to work manually...
Attachment 27438
This is excerpted from a larger flow. Two methods of geocoding produce two sets of X,Y points which I make into lines, then get the distance in miles for visualization and comparison purposes.
↧
September 13, 2013, 9:56 am
Hello All,
I am not able to find a tool (or function) that would add geographic coordinates say DD or DMS to a shape file that is in UTM (or other projected systems). When I add a field and calculate using the !shape.extent! in the window, I get only the UTM coordinates. When I use the ADDXY tool, the same.
How can I add additional geographic coordinates to the shape file without re-projecting the datasets?
Thanks,
Samuel Rajasekhar
Code:
import arcpy
import os
from arcpy import env
arcpy.env.overwriteOutput = False
arcpy.env.workspace = r'\\GISARL01\AerialStaging\BPSDelivery_080813\BPS8-Delivery\WO_80_Early_1980s_CIR_NHAP_WI_IL_IN_OH_GA_2nd\Index'
# Set the input feature (Shape file, geo-databse polygon feature, etc of Index)
inputFeature = "Foot_Prints_WGS84_NoCollar.dbf"
arcpy.AddField_management(inputFeature, "MinLat", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"MinLat Field added"
arcpy.AddField_management(inputFeature, "MinLon", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"MinLon Field added"
arcpy.AddField_management(inputFeature, "MaxLat", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"MaxLat Field added"
arcpy.AddField_management(inputFeature, "MaxLon", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"MaxLon Field added"
arcpy.AddField_management(inputFeature, "CenterLat", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"CenterLat Field added"
arcpy.AddField_management(inputFeature, "CenterLon", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")
print"CenterLon filed added"
arcpy.CalculateField_management(inputFeature, "MinLat", "!SHAPE.EXTENT.YMIN!", "PYTHON_9.3", "")
print"MinLat Field calculated"
arcpy.CalculateField_management(inputFeature, "MinLon", "!SHAPE.EXTENT.XMIN!", "PYTHON_9.3", "")
print"MinLon Field calculated"
arcpy.CalculateField_management(inputFeature, "MaxLat", "!SHAPE.EXTENT.YMAX!", "PYTHON_9.3", "")
print"MaxLat Field calculated"
arcpy.CalculateField_management(inputFeature, "MaxLon", "!SHAPE.EXTENT.XMAX!", "PYTHON_9.3", "")
print"MaxLon Field calculated"
arcpy.CalculateField_management(inputFeature, "CenterLat", "!SHAPE.CENTROID.Y!", "PYTHON_9.3", "")
print"CenterLat Field calculated"
arcpy.CalculateField_management(inputFeature, "CenterLon", "!SHAPE.CENTROID.X!", "PYTHON_9.3", "")
print"CenterLon Field calculated"
↧
↧
September 13, 2013, 11:48 am
This is a confusing one. I have two sets of data, my address points and the addresses for the phone numbers in my county. I want to compare the two to make sure I have all the phone numbers with an address in my address points and see what address points don’t have phone numbers. What makes this interesting is that not all of the address points I have will have a phone number to match and some of the addresses in the phone data are listed multiple times. The phone company data also does not list any road types only the name of the road. Businesses having more than one phone line but would only have one address point. I also have more than one address point with the same number and road name. 101 Main St for example could be an address point in all four of the cities in our county.
What phone number address does the phone company have and I don’t have?
What address points do I have the phone company doesn’t?
Thanks for any help.
↧
September 14, 2013, 9:46 am
I'm doing a viewshed analysis, using LiDAR data as my elevation dataset. My observation points are located on the top of buildings which range from 20-50 ft., but I need the observation points to have an elevation value below the top of the building (i.e @ 10 ft.). I can't set this using the SPOT field because SPOT requires values greater then the datasets elevation value. How to I get around this?
↧
September 15, 2013, 8:15 pm
Hi All,
I need to overlay a shapefile on webmap using geoprocessing service. Can anyone help me build a model which will accept full path of the shapefile and return me the feturelayer/featurecolection.
I am trying using make feature layer, but don't just understand how to create a model parapeter to pass the path of a shapefile.
Similary i have to overlay geotiff as well, using the same machanism. Please advice.
Any hwlp would be highly appreciated.
Rgds,
Pankaj Roy
↧