POST
|
THANKYOU!!! this worked like a charm, sorry late reply, i forgot i posted this!
... View more
03-27-2020
03:54 PM
|
1
|
0
|
1025
|
POST
|
So here is the scoop, I am creating a Mod map for A video game called DAYZ, However, the Mapping software it uses only uses UTM ZONE 31N. I need to use existing polygon datasets to populate things like trees but the polygons display in a different location than my terrain, my data is from UTM zone 11 N. Canada. So basically what i need is a algorithm to move and scale the shape files to line up with my height map or DEM in 31 N.
... View more
12-12-2019
08:05 AM
|
0
|
2
|
1162
|
POST
|
I have 2 sets of street data. One set is divided into small sections (essentially there is a small line every block), this data set is for asset management, IE things like manholes and sewer lines are also given a street id which is the the unique value for each line segment. The other data set is more so for operational value, It is just continuous by street name so instead of the streets being split at every block to identify its street id it is one big line for every street/avenue. What i am trying to do is add the individual street id's from the split data set to the continuous line data set without splitting any lines or having duplicate features. thanks for the help!
... View more
05-11-2018
01:43 PM
|
0
|
1
|
876
|
POST
|
THIS WORKS! import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("3D")
arcpy.env.overwriteOutput = True
# Set the current workspace
arcpy.env.workspace = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\soil"
Dir = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\ClippedMaps"
GRID="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\GRID"
GRIDlyr=arcpy.MakeFeatureLayer_management(GRID, "GRIDlyr")
# Get and print a list of TIFs from the workspace
rasters = arcpy.ListRasters("*")
for raster in rasters:
print raster
name = Dir + "\\" +"clp_"+ str(raster)
print name
arcpy.RasterDomain_3d(raster, "in_memory/temp", "POLYGON")
# Process: Central Feature
arcpy.MeanCenter_stats("in_memory/temp", "in_memory/temp1")
# Process: Select Layer By Location
arcpy.SelectLayerByLocation_management(GRIDlyr, "intersect", "in_memory/temp1", "", "NEW_SELECTION")
# Clip Raster
arcpy.Clip_management(raster,"#", name,GRIDlyr, "#", "ClippingGeometry")
# Delete in_memory
arcpy.Delete_management("in_memory")
print "processing " + raster + " complete..."
print "All processing is now finished"
... View more
02-07-2017
12:06 PM
|
1
|
0
|
136
|
POST
|
so far i have this but no return on my clip import arcpy # Set the current workspace arcpy.env.workspace = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\soil" #variables R1 = "102I16" R2 = "102I09" R3 = "102I08" R4 = "092L13" R5 = "092L12" R6 = "092L11" R7 = "092L10" R8 = "092L08" R9 = "092L07" R10 = "092L06" R11 = "092L05" R12 = "092L04" R13 = "092L03" R14 = "092L02" R15 = "092L01" R16 = "092E16" R17 = "092E15" R18 = "092E14" R19 = "092E10" R20 = "092E09" R21 = "092E08" R22 = "092E07" # Local variables: Output_Raster ="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\gridlist.tif" gridlist= R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R14,R15,R16,R17,R18,R19,R20,R21,R22 grid="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\GRID" arcpy.MakeFeatureLayer_management(grid, "gridlyr") arcpy.SelectLayerByAttribute_management("gridlyr", "NEW_SELECTION", " NTS_SNRC = 'gridlist' " ) print gridlist # Get and print a list of TIFs from the workspace rasters = arcpy.ListRasters("*", "TIF") for raster in rasters: print(raster) while raster==gridlist: arcpy.Clip_management("raster", Output_Raster, "gridlist", "", "", "NONE", "NO_MAINTAIN_EXTENT")
... View more
02-06-2017
12:53 PM
|
0
|
1
|
697
|
POST
|
I like this idea! So after i could split aprt the mosiac using split raster?
... View more
02-06-2017
10:52 AM
|
0
|
0
|
697
|
POST
|
I did selct a grid portion and created a layer form selcted data then proceeded to do the clip. Your way sounds easier. Is there a way to script this for I have a lot of maps to do.
... View more
02-03-2017
11:36 AM
|
0
|
8
|
697
|
POST
|
I have multiple geo referenced scanned maps from the 1950's, Like 60 of them. I need to clip some white space off the edges so I can apply a fishnet. So the main question is I have a 50k nts grid where the maps sit perfectly within each grid tile, When i try to clip the map to the grid it does nothing because it is reading the outside border of the nts grid. how would i get the clip to read each individual grid tile for a clip? thanks
... View more
02-03-2017
11:17 AM
|
0
|
11
|
1775
|
POST
|
SOLVED, silly mistake on my part, I had to right click in the aux.xml web page and "save as". thanks for all your help!
... View more
02-03-2017
11:10 AM
|
0
|
0
|
2347
|
POST
|
so basically it is not reading any coordinate system, ive checked off the use world file box, still not working.
... View more
01-29-2017
09:40 PM
|
0
|
1
|
2347
|
POST
|
thanks for the reply dan, recieved an ERROR -1: unidentified coordinate system, I think it is because the coordinate system info is stored outside the tif in one of the additional files.
... View more
01-29-2017
02:42 PM
|
0
|
1
|
2347
|
POST
|
I have a raster dataset at: http://www.env.gov.bc.ca/esd/distdata/ecosystems/TEI_Scanned_Maps/B01/B01-4779/ I need to geo reference the tif file but I have no clue what to do with the assosiated files as I am new to aux.xml, tfwx, etc. Some clarity on how to properly project these would be awesome! thanks a million.
... View more
01-29-2017
02:22 PM
|
0
|
7
|
4346
|
Title | Kudos | Posted |
---|---|---|
1 | 03-27-2020 03:54 PM | |
1 | 02-07-2017 12:06 PM |
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:25 AM
|