POST
|
Hi all, I have SPOT 5 imagery that I would like to display through a mosaic dataset with natural colors. My imagery doesn't have the normal structure and metadata as delivered from SPOT so I need to add the images as raster datasets. I found this thread http://forums.esri.com/Thread.asp?c=93&f=982&t=271727in the old forums and intend to use the algorithm suggested by Bill Huber in the second post. How can I best achieve this in a mosaic dataset? I've looked at the "Band Arithmetic" function but I can't figure out how to combine it with other functions to get what I want. Thanks! Nooski
... View more
01-20-2012
02:54 AM
|
0
|
0
|
2235
|
POST
|
Eric, thanks for your input! I had something like that in mind too. A potential problem with that though is that it runs the risk of creating abrupt changes in the surface around the points after �??burning�?� the dem. I will try that and also to do a focalmean on the resulting surface, probably iterating this a couple of times. As you can tell, the pixels under the input points are more important than the rest of the surface. I realize this might sound a little strange but that is the directives in this project. Do you have an idea for an interpolator that will give me the smallest difference between the input points interpolated value at that pixel? Thanks!
... View more
10-28-2011
11:29 AM
|
0
|
0
|
773
|
POST
|
Steve, Point to raster is worth a try, in combination with other tools. I do need a surface that covers the entire area. Why do the values have to be the same as the interpolated cell values? Well, it is a requirement of the project that's non negotiable:) . I also can't shift the position of the input locations. It seems like some combination of methods/tools might be needed. Basically I start with XYZ points and need to arrive at a surface, but what I do in between doesn't matter as long as the exact Z is kept in the pixels where the input points are.
... View more
10-28-2011
07:21 AM
|
0
|
0
|
773
|
POST
|
Thanks Steve, now I understand how it works. The question remains if there is a way to create a surface where the pixel "containing" the input point gets the exact value from the point. From what you have said it doesn't seem so...
... View more
10-28-2011
02:55 AM
|
0
|
0
|
773
|
POST
|
Try the Create Fishnet tool, specify Cell width and Cell height.
... View more
10-27-2011
12:56 PM
|
0
|
0
|
347
|
POST
|
Michael, The resulting surface needs to be modeled like this and also needs to be a raster. I thought about starting with a TIN but then I�??d need to convert it to raster which involves an interpolation anyway. Steve, No, my input points have no distinct pattern to them so the pixels created are not centered below the points. When you say �??The cell center is where the prediction is being made.�?�, to my understanding that means that the interpolator uses the point (which might be located at the edge of a pixel) to read the value, draws a curve to fit through the point and assigns the pixel the value that is represented by the curve at the center of the pixel? See my poor attempt at art-work in the attached image. Is this a correct understanding? And most importantly for this problem, is there any way to let the pixel get the exact value of the input point, in the image 123? Also, I missed to link to the thread in the GA forum, here it is: http://forums.arcgis.com/threads/29053-Spline-not-working-as-Exact-interpolator?p=144981#post144981 Thanks guys, I am really grateful for any help you can give here!
... View more
10-27-2011
12:47 PM
|
0
|
0
|
773
|
POST
|
Hi folks, I have to create a surface from points where the value in the surface at the location of each point is exactly the same as the z-value in the point. To me, it seems that the Spline interpolator (from Spatial Analyst) should do exactly this. Reading from the help it says: �??The Spline tool uses an interpolation method that estimates values using a mathematical function that minimizes overall surface curvature, resulting in a smooth surface that passes exactly through the input points.�?� To test if it really worked the way I expected, I used the Extract Values to points tool without the option to interpolate to draw the values from the pixels in the same location as my points. Looking at the resulting table, there is a difference of up to 20 percent in the values from my original points and the raster. I found this thread in the Geostats forum http://forums.arcgis.com/threads/29053-Spline-not-working-as-Exact-interpolator?highlight=spline, but it doesn't answer the question either. Can anyone give an explanation for why this is happening or perhaps a way to achieve the goal? Is there a setting in the Spline tool I am missing perhaps? Thanks, Nooski
... View more
10-26-2011
12:50 PM
|
0
|
9
|
3151
|
POST
|
You cannot edit multipatch vertices at 10. You can check to see that they are closed or not with the Is Closed 3d Tool, but editing vertices of a multipatch is not possible. You can move the multipatches around and up and down via the 3d Editor Toolbar. I am also including the help on the supported multipatch edits. I would recommend checking to see if the multipatches are closed, because they will error out if they are not closed. Jeff, I am facing a similar problem where I have multipatches created from TIN�??s that I want to intersect with another Multipatch with the Intersect 3D tool. My problem is that the Multipatches aren�??t closed and like you say there is no way to edit Multipatches in ArcGIS to make them closed. My workflow is: - Create TIN�??s from XYZ-points. - Extrude the TIN�??s with Extrude Between to create Multipatches representing volumes. (Create a cross-section Multipatch through digitizing a line, extruding and saving to multipatch by 3D layer to FC) - Intersect my volumes with the cross section Basically, I am following the workflow mentioned here: http://forums.arcgis.com/threads/6973-Creating-Cross-Sections-in-ArcScene?p=67522&viewfull=1#post67522 So my question is whether there is a recommended method out there that will let me create closed multipatches when I�??m going from my XYZ points to volumes? Or if there us a recommended tool for closing Multipatches through editing? Although I fear this might be a very time consuming task when the multipatches are created from extruded surfaces. Thanks, Nooski
... View more
08-24-2011
04:42 AM
|
0
|
0
|
327
|
POST
|
Hi Kevin, Use the Select Data tool, found in Data Management/General toolset. That will let you keep working on any output from the Network Analyst solvers. Cheers, Nooski
... View more
06-22-2011
01:21 AM
|
0
|
0
|
128
|
POST
|
All, Have anyone gotten the "Explode multipoints to points with LAS lidar attributes" script found here: http://arcscripts.esri.com/details.asp?dbid=16285 to work in ArcGIS 10? I know the Point to raster tool can use BLOB fields, but I need to handle the attribute-info on the individual points. Thanks, Nooski
... View more
06-09-2011
03:42 AM
|
0
|
1
|
843
|
POST
|
All, I can not get the LAS to Multipoint tool to store attributes of the LAS file. I have tried on several different datasets but all I get in the attribute table are the standard "Shape" and "PointCount" fields. When I try to store e.g. Intensity, there is no column in the attribute table. There isn't really any room for a mistake here as it is just one parameter to set, but who knows... See: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00q800000099000000.htm Has anyone else experienced this and have a solution to it? Thanks!
... View more
02-22-2011
06:46 AM
|
0
|
1
|
395
|
POST
|
Adding to my post above, there's an explanation in the help why you need to handle the result object. Have a look here: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z0000000n000000.htm In the note under the "Getting results from a tool" heading, it says: The result object's getOutput method returns a Unicode string from result objects that have output values. This is an important consideration when running a tool such as GetCount, which provides a count of records in a table, or CalculateDefaultClusterTolerance, which provides a cluster tolerance value. To convert the value to the expected type, it will have to be converted from a unicode string using built-in Python functions such as int() or float().
... View more
02-11-2011
05:46 AM
|
0
|
0
|
821
|
POST
|
Hi Hydronaut, What specific problem do you have with the code in the link? I believe that with some of the changes suggested in the comments in the link you sent, this should work. I was removing empty layers from mxd�??s with a similar method. What I found (in v.10) is that you need to make an integer of the results object from GetCount_management. So if you look at comment #5 in the link, which says: result = arcpy.GetCount_management(fc)
print result
if result == 0:
print fc + �?� is empty.�?�
arcpy.delete_management(fc)
Change that to result = arcpy.GetCount_management(lyr)
count = int(result.getOutput(0))
if count == 0:
arcpy.mapping.RemoveLayer(df, lyr)
print "Removing " + lyr
else:
print "Keeping " + lyr
In the code above, substitute arcpy.mapping.RemoveLayer for arcpy.Delete_management. Also, you might need to handle FC's inside feature datasets like the first comment suggests.
... View more
02-11-2011
12:37 AM
|
0
|
0
|
821
|
POST
|
Thanks anant! That works! What is the plan for 10, will this be fixed in the future? Nooski
... View more
01-27-2011
01:50 AM
|
0
|
0
|
224
|
POST
|
All, arcpy.ListDatasets and I can�??t seem to agree. I want to list all CAD datasets in a folder and according to the help I should be able to use arcpy.ListDatasets with the �??CAD�?� keyword. If I write this, I get all my CAD data but also everything else:
import arcpy
from arcpy import env
env.workspace = "C:\\data"
datalist = arcpy.ListDatasets("*", "All")
for data in datalist:
print data
If I write the following code, nothing gets listed:
import arcpy
from arcpy import env
env.workspace = "C:\\data"
cadlist = arcpy.ListDatasets("*", "CAD")
for cad in cadlist:
print cad
In the �??Listing data�?� section it lists CAD as one of the keywords. (http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z00000011000000.htm) But in the ListDatasets topic there is no mention of the CAD keyword. ( http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//000v0000001m000000.htm ) For now, I get by with this, which isn't optimal. cadlist = arcpy.ListDatasets("*.dwg", "All") Has anyone successfully used the CAD keyword and could you then help me out with what works? ESRI, I think the help section needs to be corrected. Thanks!
... View more
01-21-2011
12:28 AM
|
0
|
3
|
739
|
Title | Kudos | Posted |
---|---|---|
1 | 07-17-2020 09:21 AM | |
1 | 06-11-2019 12:17 AM |
Online Status |
Offline
|
Date Last Visited |
yesterday
|