POST
|
I am using Pro to visualize MODIS Thermal Anomalies and Fire 8-day composite (MOD14) HDF's -- FireMask and MaxFRP. Last year, when I viewed the same data in an older version of Pro, I am quite sure I was able to see variation in at least one of these variables across the 8 days of the composite. Now, I am adding the hdf as a Multidimensional Raster, which I seem to have to view using the time slider. Every day looks identical in both variables (in Properties, Number of Days = 0). Alternatively, adding it as a Multiband Raster, every band has the same name in Extract Bands and looks identical on extraction. Is there any way to just unpack the hdf directly into 8 separate single-band rasters per variable?
... View more
09-30-2019
11:01 AM
|
0
|
0
|
455
|
POST
|
That *is* what I needed. Apparently you cannot assign it to a variable, which (mistakenly) led me to believe that approach was out of date. Really what I needed is: df = pandas.DataFrame(_layer.query().sdf) And now I have that -- thanks!
... View more
04-30-2019
06:03 AM
|
1
|
2
|
5440
|
POST
|
I understand that SpatialDataFrame has been replaced by a "spatially enabled dataframe." It looks like the former was a function in arcgis.features while the latter goes through pandas. I cannot get either to work! On advice from another thread, I rolled back my pandas to 0.23.4, but that has not helped (arcgis api version 1.6). Using tab-complete, pandas.DataFrame.spatial does not appear to have any methods. Apparently arcgis.features.GeoAccessor has a .from_layer method that yields a pandas dataframe, but I wasn't able to define a feature layer from a url. Old Way (SpatialDataFrame): New Way ("spatially enabled dataframe" via pandas): Third way (GeoAccessor): from arcgis.features import FeatureLayer, GeoAccessor
FL = FeatureLayer('http://[my institution].maps.arcgis.com/home/item.html?id=fbcae48fb3e64ad886ee84ddb225f39a')
df = GeoAccessor.from_layer(FL) That FL does not have properties and cannot be mapped, so it's no surprise from_layer fails as well. Any thoughts?
... View more
04-29-2019
03:31 PM
|
1
|
4
|
6325
|
POST
|
I downgraded to 0.23.4, but it did not solve the problem. pandas.DataFrame.spatial does not seem to have a method .from_layer, at least not according to tab-complete. (using arcgis 1.6)
... View more
04-29-2019
09:38 AM
|
0
|
0
|
765
|
POST
|
I have a rather complicated symbology I would like the user to adopt. I'd prefer not to save it as a template and have them import it, not having trust that they will use exactly the same field names for example.
... View more
11-01-2018
06:46 AM
|
0
|
0
|
229
|
POST
|
The process is quite straightforward in Raster Calculator, but if you have already committed to the Composite Bands/Mosaic Dataset/Image Analysis window, it is completely opaque. Raster Calculator doesn't work in that context, and while there are functions for "Apparent Reflectance" and -- inexplicably -- "Sentinel 1 Radiometric Calibration," it's not at all clear how one would put together such a function for Landsat.
... View more
10-05-2018
05:12 AM
|
1
|
0
|
717
|
POST
|
Does anyone know the scope of the Hands-On: Python Scripting for GIS People pre-conference workshop at UC 2017? Will it use the python API or is it arcpy via the desktop?
... View more
03-28-2017
06:24 AM
|
1
|
2
|
1227
|
POST
|
Much stripped down, this works: my_line_geom = {'paths' : [[[-78.8,38.8],[-78.8,38.2]]]} my_line = arcgis.geometry.Polyline(my_line_geom) my_line_feature = arcgis.features.Feature(my_line) my_line_feature_set = arcgis.features.FeatureSet([my_line_feature]) map1 = gis.map('Orkney Springs') map1.draw(my_line_feature_set) map1 I have not been able to draw both features of a two-feature feature set. Only the second one shows up. Nor have I had any success setting the symbology.
... View more
02-25-2017
02:56 PM
|
0
|
1
|
1059
|
POST
|
If you actually have a polyline geometry, try following it with lines = arcgis.features.Feature(fc_lines) lines1 = arcgis.features.FeatureSet([lines]) and then map1.draw(lines1) instead of add_layer. If not, to get from the text string to a geometry, you'd need fc_lines = arcgis.geometry.Polyline(fc_string_thing) before all that. Hope this works; I got it working for pandas-dataframe-to-polygon.
... View more
02-21-2017
07:44 AM
|
0
|
4
|
1059
|
POST
|
I have read the "Smart Mapping" Guide about how to set a renderer, but I cannot find a list of possible renderers, nor can I find how to set custom parameters like a polygon border with a chosen weight and color with no fill. Is there a detailed list somewhere? I don't see it in the docs. I've tried {"color":[208,32,144]} in map.add_layer( ) for example. That passes with no error, but the polygons show up blue (ignoring my suggested fuchsia-pink color). I'm inclined to believe it's possible, since line [18] of Using and updating GIS content | ArcGIS for Developers displays an example including such symbology parameters.
... View more
02-21-2017
07:31 AM
|
0
|
0
|
1021
|
POST
|
Specifically, is it possible to get percentile cutoffs and CoV? Can one define a custom statistic?
... View more
02-08-2017
07:52 AM
|
0
|
1
|
736
|
POST
|
Thanks Rohit. I had tried converting it to a Layer, but I overlooked FeatureSet. It works!
... View more
02-07-2017
06:38 AM
|
0
|
0
|
688
|
POST
|
I am able to produce: so to the best of my knowledge, I am authenticated. I have also been able to display maps. The specific line giving the no-spatial-analysis error is: PlotDiscs = create_buffers(PlotPoints, [11.2838], None, 'Meters', output_name='PlotDiscs400m') where PlotPoints is PlotPoints also displays as it ought to on the map. ps -- changing the url from the one they gave me when I signed up for a developer account to arcgis.com as your link instructs for an organization makes no difference. I am part of an organization, but I signed up as an individual, and we haven't yet connected me to them. Full error message: AMani-esristaff wrote: Jill Kelly could you paste your code and the error message that is returned? Also to connect to a GIS, you need to specify the URL, either to ArcGIS Online or a Portal for ArcGIS as shown here: https://developers.arcgis.com/python/guide/working-with-different-authentication-schemes/#ArcGIS-Online-using-an-Organization-account
... View more
02-01-2017
11:00 AM
|
0
|
0
|
1567
|
POST
|
Is this error more likely due to technical issues (that I'm trying to create a buffer in an unprojected/geographic crs, following the example in the california forest fires demo) or access (having failed to properly authenticate or have "credits" or do the equivalent of "check out an extension" or some such)? gis = GIS(username='mine',password='mine') What else would I need for the GIS to support spatial analysis?
... View more
02-01-2017
06:46 AM
|
0
|
6
|
2222
|
POST
|
"item" seems a little vague. I can, for example, add_layer a set of points: df1 = pandas.read_csv('CTpolygon.csv', header=0, names=['longitude', 'latitude']) rows = np.random.choice(df1.index.values, 499) df = df1.ix[rows] CTpoints = gis.content.import_data(df) map1 = gis.map('Connecticut') # focus on CT map1.basemap = 'topo' map1.add_layer(CTpoints) map1 How is a polygon feature any different in terms of being an "item"?
... View more
01-26-2017
12:03 PM
|
0
|
2
|
688
|
Title | Kudos | Posted |
---|---|---|
1 | 04-30-2019 06:03 AM | |
1 | 04-29-2019 03:31 PM | |
1 | 10-05-2018 05:12 AM | |
1 | 03-28-2017 06:24 AM | |
1 | 01-25-2017 10:50 AM |
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:24 AM
|