|
POST
|
This is a problem I come across regularly and feel there's room to save a step or two. 1.) Field data comes back from GPS or field notes in a table containing UTM Zone, Easting, and Northing fields. There may be multiple UTM Zones (e.g. Zone A and B). 2.) Display XY data, using Zone A as CRS. 3.) Export all data to a single feature class, using the chosen UTM Zone (A). 4.) Set Data Frame CRS to Zone B. 5.) Select and export records in Zone B to new feature class, using Data Frame CRS. 6.) Start editing and delete records in Zone B from the feature class in Zone A. 7.) Merge Zone A and Zone B feature classes into a single feature class with common CRS. As you can see, this is somewhat convoluted. How do you handle tabular data in multiple CRSes?
... View more
01-11-2016
11:04 AM
|
1
|
27
|
12379
|
|
BLOG
|
The "Copy Text from Picture" function is a great tip. Thanks.
... View more
01-11-2016
10:09 AM
|
0
|
0
|
333
|
|
POST
|
There's no attachment. Have you tried this as a custom CRS?
... View more
01-08-2016
03:19 PM
|
2
|
1
|
4506
|
|
POST
|
I'm not an expert, but I think it depends on where you live. Or at least it definitely did 20 years ago when I lived in an area without 911. It may be the case now that everyone has access to 911 somewhere, but they'll simply redirect your call to the appropriate local emergency service provider. I will go out on a limb and say that Canadian 911 services infrastructure lags behind that in America. I'm sure someone in the emergency field would know - we can all rest assured that is not me! edit: it looks like about 98% of Canadians are now covered by 911 services. But, as I said before, I'm sure a good portion of those services are acting remotely and would simply relay the call locally.
... View more
01-08-2016
02:22 PM
|
1
|
3
|
4465
|
|
POST
|
Just because I've been beckoned by Chris Donohue, GISP, I'll add my largely uninformed two cents. I can't imagine that such a database exists for the public. If it did, any business could immediately call every resident. Also, I'm not sure if your plan is to automatically call or send texts, but for texts, at least, you do need prior consent from the recipient to conform with anti-spam laws, which means building your own database. edit: I suppose one way to get your phone numbers through somewhat icky backdoor methods would be to screen scrape a website like White Pages for Canada - 411.ca with your list of addresses, but make sure it doesn't violate their terms of service (which I'm sure it would).
... View more
01-08-2016
01:02 PM
|
1
|
5
|
4465
|
|
POST
|
Here's a start that gets me close. I'm not 100% sure that it's even possible to fit this dataset exactly, in all parts of the raster at the same time. The pertinent information can be found here (I'll note that they indicate that you can use the lat/long file [txt and xlsx attached], but you'd have to do a fair bit of manipulation to tie the points to the ASCII grid values): Grid Specification: North Polar Stereographic ni (# rows) = 706, nj (# cols) = 706, bottom left grid point = (1,1) Grid rotation -10° Pole position = (353, 353) Grid resolution at 60°N = 23812.5 m So, your header should look like: NCOLS 706
NROWS 706
XLLCENTER 1
YLLCENTER 1
CELLSIZE 23812.5
NODATA_VALUE 0 [or another value like 999 if you want to include values of 0] You also need to define a custom coordinate reference system to make it line up with the world, once you've converted ASCII To Raster: false_easting and false_northing = 8405812.5 = cellsize * pole position central_meridian = 10 The final result is:
... View more
01-08-2016
11:59 AM
|
1
|
1
|
2998
|
|
POST
|
'Gewaesser Rechtswert ohne 32 >= 494225.671350 AND Gewaesser Hochwert in UTM Koordinaten >= 6028859.486114 AND Gewaesser Rechtswert ohne 32 <= 502301.705257 AND Gewaesser Hochwert in UTM Koordinaten <= 6036832.306366'. Remove the spaces in your field names.
... View more
01-08-2016
09:45 AM
|
1
|
0
|
3431
|
|
POST
|
Can you provide the line where you specify the workspace? On a side note, why does your code flipflop between 'arcpy' and 'gp'? Are they both pointing to arcpy? Have you tried this with the superior data access module SearchCursor? Here's how to post in code tags (it's a bit convoluted): Posting Code blocks in the new GeoNet rows = gp.SearchCursor(workspace + "\\parcels_ut_town_list.dbf")
for row in rows:
town = row.town
# Create Output Directory
towndir = town.strip()
# Create Output Directory
gp.AddMessage("Processing" + " " + town)
print "Processing" + " " + town
gp.AddMessage("Creating Directory For" + " " + town)
print "creating outfolder for "+town
outfolder = gp.CreateFolder_management(workspace, towndir) ##Sometimes quitting here
expression = "TOWN = " + "'" + town.strip() + "'"
parcels = outfolder + "\\parcels"+towndir+".kmz"
# Process: Make Feature Layer
print "making feature layer" ##-And sometimes quitting around the next line-
arcpy.MakeFeatureLayer_management(parcelsProject_shp, parcelsProject_Layer, expression, "", "FID FID VISIBLE NONE;TOWN TOWN VISIBLE NONE;COUNTY COUNTY VISIBLE NONE;GEOCODE GEOCODE VISIBLE NONE;STATE_ID STATE_ID VISIBLE NONE;MAP_BK_LOT MAP_BK_LOT VISIBLE NONE;PARENT PARENT VISIBLE NONE;PROP_LOC PROP_LOC VISIBLE NONE;PROPLOCNUM PROPLOCNUM VISIBLE NONE;TYPE TYPE VISIBLE NONE;FMUPDORG FMUPDORG VISIBLE NONE;FMUPDAT FMUPDAT VISIBLE NONE;FMSRCORG FMSRCORG VISIBLE NONE;Shape_STAr Shape_STAr VISIBLE NONE;Shape_STLe Shape_STLe VISIBLE NONE;Shape Shape VISIBLE NONE;Shape_ST_1 Shape_ST_1 HIDDEN NONE;Shape_ST_2 Shape_ST_2 HIDDEN NONE")
# Process: Layer To KML
gp.AddMessage("Converting town to KML")
print "Converting town to KML"
arcpy.LayerToKML_conversion(parcelsProject_Layer, parcels, "0", "false", parcelsProject_Layer, "1024", "96", "CLAMPED_TO_GROUND")
#Copy Metadata
gp.AddMessage("Copying parcels Metadata For" + " " + town)
print "Copying parcels Metadata For" + " " + town
shutil.copy(metat, outfolder + "\\" + "parcels.txt")
gp.AddMessage("Done Processing parcels For" + " " + town)
print "Done Processing parcels For" + " " + town
#Create The Zipfile
#gp.AddMessage(output + "\\Parcels_" + towndir)
zpfile = zipfile.ZipFile(output + "\\Parcels_" + towndir + ".zip", 'w')
zpfile.write(outfolder + "\\parcels"+towndir+".kmz", "Parcels_" + towndir + ".kmz", zipfile.ZIP_DEFLATED)
zpfile.write(outfolder + "\\parcels.txt", "Parcels_" + towndir + ".txt", zipfile.ZIP_DEFLATED)
zpfile.close()
#Delete Uncompressed Files
shutil.rmtree(outfolder)
gp.AddMessage("All Done Processing" + " " + town)
del rows
... View more
01-08-2016
09:22 AM
|
1
|
7
|
2988
|
|
POST
|
Can you post more of your code (in code tags) for context?
... View more
01-07-2016
02:57 PM
|
0
|
10
|
2988
|
|
POST
|
If you're going to use projectAs, it needs to be on a polygon-by-polygon basis, not the list: >>> inSr = arcpy.SpatialReference(4326)
... outSr = arcpy.SpatialReference(3857)
... features = []
... features2 = []
... feature = [[1,2],[2,2],[2,3],[1,2]]
... polygon = arcpy.Polygon(arcpy.Array([arcpy.Point(*coords) for coords in feature]), inSr)
... polygon2 = polygon.projectAs(outSr)
... features.append(polygon)
... features2.append(polygon2)
... arcpy.CopyFeatures_management(features,r'in_memory\poly')
... arcpy.CopyFeatures_management(features2,r'in_memory\poly2')
... View more
01-07-2016
02:35 PM
|
2
|
2
|
4489
|
|
POST
|
Try: elseif InStr([Modified], vbNewLine) then Score =7
... View more
01-07-2016
11:30 AM
|
1
|
2
|
1569
|
|
POST
|
I want to get into the sustainable development, natural resource management, and environmental services industries, companies that work with the government or private sector, and eventually work my way up to a natural sciences manager but it's really hard to know what these big companies want. I would start finding the website for any small to mid-size environmental consulting companies in your area, see if they offer GIS services, and show up resume-in-hand. This way, you can usually speak to someone (sometimes a manager, not a HR person), and see if they have an opening - if not, they know you're around and they know you're motivated. If no one can see you, no harm done, move onto the next. A lot of these companies may have a need for casual or seasonal help, especially in spring/summer when field work is ramping up/ongoing. I would still apply for large company postings, but you'll likely get screened out immediately with no direct experience. Other than that, staying active by volunteering in some GIS capacity, like you're doing, is great. Unfortunately, the first job is the hardest (obviously) so just prepare for rejection and legwork. PS - I'm speaking as an environmental consultant, and know that people get hired where there isn't a job posting all the time.
... View more
01-07-2016
10:36 AM
|
4
|
0
|
12417
|
|
POST
|
The Feature to Polygon tool is where it's always been, AFAIK. Data Management -> Features Are you sure the difference isn't that you used to have ArcInfo (Advanced) licensing, and now do not?
... View more
01-07-2016
10:11 AM
|
2
|
1
|
1606
|
|
POST
|
The data type you use depends on what format you reasonably expect and require from your user. If you can rely on your user to type in or (copy and paste) a full polygon array, like: [[[6,5],[7,8],[8,7],[6,5]]] ...then you could simply use a string. Similarly, you could input a table of coordinates, or type and parse an input table like Luke mentions. If you go with the above array option, the following slightly modified script should work: import arcpy, ast
array = arcpy.GetParameterAsText(0) # a simple string
# A list of features and coordinate pairs
feature_info = ast.literal_eval(array)
# A list that will hold each of the Polygon objects
features = []
for feature in feature_info:
# Create a Polygon object based on the array of points
# Append to the list of Polygon objects
features.append(
arcpy.Polygon(
arcpy.Array([arcpy.Point(*coords) for coords in feature])))
# Persist a copy of the Polyline objects using CopyFeatures
arcpy.CopyFeatures_management(features, r"C:\junk\poly.shp")
... View more
01-06-2016
09:42 AM
|
1
|
0
|
4591
|
|
POST
|
Feature templates store default symbology and attributes, but not geometry. As Joe indicates, copying and pasting in an edit session would be the fastest way to duplicate geometry.
... View more
01-05-2016
01:14 PM
|
1
|
0
|
2955
|
| Title | Kudos | Posted |
|---|---|---|
| 1 | 11-25-2015 01:51 PM | |
| 1 | 08-30-2013 02:22 PM | |
| 1 | 04-12-2011 11:19 AM | |
| 1 | 09-17-2021 09:43 AM | |
| 1 | 04-04-2012 12:05 PM |
| Online Status |
Offline
|
| Date Last Visited |
07-15-2023
12:11 AM
|