Latest Contributions by anTonialcaraz

Just to add: The script works with the FCs in Geographic Coordinates (degrees) using the global extent (below) I would expect the same in Mollweide Projected (meters) especially as you are measuring distances... arcpy.env.extent = arcpy.Extent(-180.0, -90.0, 180.0, 90.0) env.workspace = r"H:\OUTPUT.gdb" fcdist = arcpy.ListFeatureClasses("") for fc in fcdist:     output = EucDistance(fc, 5, 0.5, "")     output.save(env.workspace + "\\" + fc + "_250km")‍‍‍‍‍‍‍‍‍‍ ... View more

To be more precise when using the tool in ArcMap session:       It works when I use the "same as layer..." (in this case I use a Bathymetry raster layer) in the Processing Extent option within Environment Settings.       If I use " As specified Below" option and I enter the extent manually it does not work. I don't quite get it to be honest ... View more

As I previously mentioned, no problem at all when using the tool in ArcMap and setting the  global Extent ... View more

Thanks Dan. The Euclidean Distance calculation is the end part of a larger script. In the previous step I do project the FCs to Mollweide. It is just when using the Extent that does not work but I guess it may still be related with the problem? # reproject to Mollweide and create 250km buffer fcprj = arcpy.ListFeatureClasses() for fc in fcprj:     mollweide = arcpy.SpatialReference("Mollweide (world)")     arcpy.BatchProject_management(fc, OUTworkspace, mollweide, "", "") fcdist = arcpy.ListFeatureClasses("*1") for fc in fcdist:     output = EucDistance(fc, 250000, 10000, "")     output.save(OUTworkspace + "\\" + fc + "_250km")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Projected Coordinate System: World_Mollweide Projection: Mollweide False_Easting: 0.00000000 False_Northing: 0.00000000 Central_Meridian: 0.00000000 Linear Unit: Meter Geographic Coordinate System: GCS_WGS_1984 Datum: D_WGS_1984 Prime Meridian: Greenwich Angular Unit: Degree ... View more

Hi again, Thanks to everyone that have had a look into this. I guess it is a tricky one. Just to add on my question/query: I have a macro written in VB by someone at the company quite a while ago that does what I intend to do (improve) in Python. I have no knowledge of VB so not really sure how to translate this into Python. I assume that if this can be done in VB it should also by possible in Python? Please see the VB code below and if you can help with this it would be greatly appreciated. Const STYLE_FIELD = "CSID" Const LABEL_FIELD = "Description" ' -------------------------------------------- Public Sub LabelStylesFromAnotherField() Dim pMxDoc As IMxDocument Dim pGeoFtrLyr As IGeoFeatureLayer Dim pLegInfo As ILegendInfo Dim pLegGrp As ILegendGroup Dim pLegCls As ILegendClass Dim f As Long     Set pMxDoc = ThisDocument     Set pGeoFtrLyr = pMxDoc.FocusMap.Layer(0)     Set pLegInfo = pGeoFtrLyr.Renderer     Set pLegGrp = pLegInfo.LegendGroup(0)         For f = 0 To pLegGrp.ClassCount - 1         Set pLegCls = pLegGrp.Class(f)         pLegCls.Label = GetStyleLabel(pLegCls.Label, STYLE_FIELD, LABEL_FIELD, pGeoFtrLyr)     Next f         pMxDoc.UpdateContents         MsgBox "Updated Symbol Labels", vbInformation, ""     End Sub Private Function GetStyleLabel(sStyle As String, sStyleField As String, sLabelField As String, pGeoFtrLyr As IGeoFeatureLayer) As String Dim pFtrCls As IFeatureClass Dim lStyFldIdx As Long Dim lLblFldIdx As Long Dim lID As Long Dim pQryFltr As IQueryFilter     Set pFtrCls = pGeoFtrLyr.DisplayFeatureClass     lStyFldIdx = pFtrCls.FindField(sStyleField)     lLblFldIdx = pFtrCls.FindField(sLabelField)     If lStyFldIdx = -1 Or lLblFldIdx = -1 Then         GetStyleLabel = sStyle         Exit Function     End If         Set pQryFltr = New QueryFilter     pQryFltr.WhereClause = sStyleField + " = '" + sStyle + "'"         lID = pFtrCls.Select(pQryFltr, esriSelectionTypeIDSet, esriSelectionOptionOnlyOne, Nothing).IDs.Next         If lID = -1 Then         GetStyleLabel = sStyle         Exit Function     End If         GetStyleLabel = CStr(pFtrCls.GetFeature(lID).Value(lLblFldIdx))     End Function‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I see what you mean. Thanks so much Dan. rasterlist5 = arcpy.ListRasters(StageAge + "_*") for raster in rasterlist5:     print raster output_mean = CellStatistics(rasterlist5, "MEAN", "DATA") output_mean.save(env.workspace + "\\" + raster[:11] + "ann")‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I know it works like as below. Just thought there would be a more elegant and economic way to put it using ListRasters rasterlist5 = [StageAge + "_t_srf_apr", StageAge + "_t_srf_aug", StageAge + "_t_srf_dec", StageAge + "_t_srf_feb", StageAge + "_t_srf_jan", StageAge + "_t_srf_jul", StageAge + "_t_srf_jun", \ StageAge + "_t_srf_mar", StageAge + "_t_srf_may", StageAge + "_t_srf_nov", StageAge + "_t_srf_oct", StageAge + "_t_srf_sep"] output_mean = CellStatistics(rasterlist5, "MEAN", "DATA") output_mean.save(OUTWorkspace + "\\" + raster[:11] + "ann")‍‍‍‍‍‍‍‍‍‍ ... View more

Thank you Dan. yes, the print statement gives me the list of rasters I want to calculate the mean from. I thought I was giving the script the list of rasters with the variable rasterlist. I don't quite understand why it only processes the last raster in the list. ... View more

just to add that I'm using 10.2.1 at the moment. many thanks ... View more