Latest Contributions by JoshuaAhmed

‎09-16-2016

Thanks, Stephen. Is there a way to see these dates without having access to the account (i.e., can you see them on the published story map?). Josh

‎09-16-2016

Hi, I'm looking to use story maps as an academic tool but need to have a way of monitoring when the latest updates were made to the story. Is there a way to add this to the map automatically when its published? Any help would be much appreciated Josh

‎10-02-2014

Hi all, I have a similar script where I am going through a file and converting all .txt files into events layers, and finally into points files. The points are X and Y coords. I have tried changing all manner of things yet I still get a Runtime error: object: error in executing tool message. Any ideas where the script is going wrong. I know it's something in the MakeXYEventLayer_management but cannot figure out where... Any help will be welcomed. ------------------------------------------------------------------- import arcpy, sys, os from arcpy import env from arcpy.sa import * arcpy.CheckOutExtension("Spatial") arcpy.env.overwriteOutput = True #Set the script to be able to overwrite files Path_table = "J:/Annual_time_series/NDVI_images/Beni1/Centrelines/" #path for root, dirs, files in os.walk(Path_table): for name in files: beni3_reclass = os.path.join(root,name) #name if beni3_reclass.endswith('_CL.txt') == True: tilename = name[0:-4] #changed to keep my filenames print tilename + " being processed." table = Path_table + name in_x_field = "cl_x" in_y_field = "cl_y" out_layer = Path_table + tilename + "_layer" #saved_layer = Path_table + tilename + ".lyr" #Convert txt file table to an events layer arcpy.MakeXYEventLayer_management (table, in_x_field, in_y_field, out_layer, {"J:/Annual_time_series/NDVI_images/coordinate_system/WGS 1984 UTM Zone 19N.prj"}, {}) print arcpy.GetMessages() #save to featureclass - this makes it a permenant shape points file Path_points = "J:/Annual_time_series/NDVI_images/Beni1/Centrelines/CL_coords" in_features = out_layer out_featureclass = Path_points + tilename + "_CL_pts.shp" arcpy.CopyFeatures_management (in_features, out_feature_class, {}, {}, {}, {}) ------------------------------------------------------------------------- The error message is: Traceback (most recent call last): File "E:/CLtxt_to_points.py", line 22, in <module> arcpy.MakeXYEventLayer_management (table, in_x_field, in_y_field, out_layer, {"J:/Annual_time_series/NDVI_images/coordinate_system/WGS 1984 UTM Zone 19N.prj"}, {}) File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 6348, in MakeXYEventLayer raise e RuntimeError: Object: Error in executing tool Thanks, Josh

‎04-11-2014

Hi all, I have been creating composite banded images from Landsat imagery using to produce an RGB true-colour image, and sometimes pan-sharpening these images. It all appears to work ok and I can view the image in Arcmap once it's done. However, when I transfer the image to another computer and open it in QGIS it appears as the image but with a blue or red cover on the top (I assume from one of the bands). I also have a problem viewing the image in Coreldraw, where the image doesn't appear as a composite. I hope this can be fixed and look forward to hearing any suggestions.

‎02-06-2014

Hi Josh, Try passing a list of the rasters to the Composite Bands tool. Try the following: import arcpy, sys, os from arcpy import env from arcpy.sa import * arcpy.env.overwriteOutput = True path_images = "J:/Proj/Annual_time_series/Landsat Images/1987/1_001_070/" rasterList = [] for root, dirs, files in os.walk(path_images): for name in files: beni87 = os.path.join(root,name) if beni87.endswith(('B10.TIF', 'B20.TIF', 'B30.TIF', 'B40.TIF', 'B50.TIF', 'B60.TIF', 'B70.TIF')): rasterList.append(beni87) tilename = name[0:-7] print tilename + "being processed." out_raster = path_images + "\\" + "composite.tif" arcpy.CompositeBands_management(rasterList, out_raster) Also, when posting code, be sure to enclose it with CODE tags using the # symbol above. Thanks, that was perfect! Been bugging me all day. Josh

‎02-05-2014

Hi, I have a folder containing 7 tif rasters that I would like to combine into one composite tif. I have to do this for many images so I have written a python script to speed up to process. However, I'm having some trouble collecting all of the bands from the folder and creating a composite. Any help would be much appreciated. import arcpy, sys, os from arcpy import env from arcpy.sa import * arcpy.env.overwriteOutput = True #set the script to overwrite files #the path to where the images are stored - will need to change for each set of images as the end folders are different path_images = "J:/Proj/Annual_time_series/Landsat Images/1987/1_001_070/" for root, dirs, files in os.walk(path_images): for name in files: beni87 = os.path.join(root,name) if beni87.endswith[('B10.TIF', 'B20.TIF', 'B30.TIF', 'B40.TIF', 'B50.TIF', 'B60.TIF', 'B70.TIF')] == True: #collect band 1,2 and 3 image tilename = name[0:-7] #remove the band number and .tif from the end of each image to be renamed print tilename + "being processed." #the list of banded images need to be put here, or the file containing the list needs to be put here in_rasters = path_images + name #input the bands from the file out_raster = path_images + tilename + "composite.tif" #insert the name of the image between the quotation marks #start composite bands processing arcpy.CompositeBands_management(in_rasters, out_raster)

Online Status	Offline
Date Last Visited	‎11-11-2020 02:23 AM

My Ideas

Latest Contributions by JoshuaAhmed

Re: Is there a way to automatically produce a last updated tag to a story map?

Is there a way to automatically produce a last updated tag to a story map?

Re: Create XY Event Layer / Save to Layer file / Copy features to shapefile

Make composite bands image, problems outside of arcmap

Re: Arcpy create composite images from several single band tifs

Arcpy create composite images from several single band tifs