topic Re: how to create own surface profile service? in Geoprocessing Questions

how to create own surface profile service?

Divyaprasad — Fri, 30 Apr 2010 08:53:11 GMT

Hi

The GPservice used to get surface profile in ESRI samples is " http://sampleserver2.arcgisonline.com/ArcGIS/rest/services/Elevation/ESRI_Elevation_World/GPServer/ProfileService"

How do i create my own service to get the same result with my data than ESRI sample data..

Or is it possible to get surface profile from my data using esri's profileservice? cos when i tried i ended up getting errors. im confused whether its my data or service.. im using DEM.

could somebody please give solution for this please. Im totaly messed up.

Re: how to create own surface profile service?

AndrewRettig — Tue, 04 May 2010 15:47:38 GMT

I would like to create my own service as well. Any word on this?

Andrew

Re: how to create own surface profile service?

Divyaprasad — Thu, 20 May 2010 16:45:05 GMT

why is it that whenever a question is posted on surface elevation profile service creation for oneself, there is no answer !!! i mean is it something which is not possible ?? kindly acknowledge..

Re: how to create own surface profile service?

StefanP__Jung — Fri, 12 Nov 2010 09:15:25 GMT

I am also very interested to create my own surface profile GPService.

Did anyone get some more information here?

Re: how to create own surface profile service?

GeraldLee__GISP — Tue, 01 Mar 2011 15:22:31 GMT

i just have to chime in as well. this would be nice to know.

Re: how to create own surface profile service?

dcopple — Wed, 02 Mar 2011 14:21:25 GMT

I would like to know how to use ESRI's profile service using javascript. I am actually trying to create a slope calculator. However to try to answer other questions. I would think you could create a geometry service (which usually accounts for x,y coordinates and then be able to calculate a z value at that point based on your data) In a profile your line would have several x,y coordinates and calculate the z values at each one of those x,y points. The result would be your profile.....

Now anyone know how to write javascript to use ESRI's profile service using Javascript API?

Re: how to create own surface profile service?

SeanCook — Wed, 09 Nov 2011 20:54:55 GMT

TTT!

I would also like to see an answer to this. Is ESRI actively ignoring these?

Re: how to create own surface profile service?

SertaçTAKIL — Wed, 28 Dec 2011 13:49:40 GMT

I also would like to know how it's done?
Seems no one from ESRI interested in this subject.

Re: how to create own surface profile service?

IainCampion — Wed, 29 Feb 2012 05:02:58 GMT

ditto,

come on esri .. chime in for f(**( sake ... even if the answer is no.

Re: how to create own surface profile service?

KyleShannon — Sat, 11 Dec 2021 09:47:36 GMT

I managed to hack this together at one point, but found it was too slow for a web service. I rewrote the script using gdal for raster access and it was hundreds of times faster. The code for the arcpy script is below, with an example output graph. The commented code was for grabbing points from a feature set (for a toolbox or gp service). This ugly, but it works. disclaimer: This is hacked, old and ugly. sorry.

from polyline import *
import matplotlib.pyplot as plt

def fetch_cell_value(x, y):
    xy_string = "{0} {1}".format(x, y)
    band = "1"
    
    arcpy.GetCellValue_management(r, xy_string, band)
    results = arcpy.GetMessages()
    lines = results.split('\n')
    value = float(lines[2])
    return value

#for arc script
#fs = arcpy.GetParameter(0)
#rows = arcpy.SearchCursor(fs)

# Assume feature set has only one feature and one shape.
#row = rows.next()
#shape = row.getValue("SHAPE")
#arcpoints = shape.getPart(0)

#points = list()
#for p in arcpoints:
#    points.append([p.X, p.Y])

raster = "c:/data/someraster.tif"
r = arcpy.Raster(raster)

#test points
points = [[290626.474, 4875792.8594],
          [282358.29,4863026.666],
          [297836.3843, 4870699.62]]

p = poly_line(points)

n_points = 100
#if p.length() / n_points < cell_size:
#    skip = cell_size
#else:
#    skip = int(p.length() / n_points)
skip = p.length() / n_points

dists = list()
i = 0
while i * skip < p.length():
    dists.append(i * skip)
    i += 1

values = list()
for d in dists:
    x,y = p.d_line(d)
    v = fetch_cell_value(x,y)
    values.append(v)
print("Time:{0}".format(time.clock() - start))

plt.plot(dists, values)
plt.xlabel("Distance (m)")
plt.ylabel("Elevation (m)")
plt.title("Elevation Profile")
plt.grid(True)
plt.show()

The poly_line is a little class I wrote to handle densifying lines. It assumes euclidean geometry (planar) so doesn't handle earth curvature. You could probably use some type of densify function for an arc type polyline and use the points in the shape. I didn't try it.

poly_line:

import math

##
# Class representing a segment of a polyline
#
class segment:
    p1 = [0,0]
    p2 = [0,0]
    ##
    # @brief Construct a segment
    #
    def __init__(self, p1, p2):
        self.p1 = p1
        self.p2 = p2
    ##
    # @brief Get the length of the segment
    #
    # @todo Test for /0
    #
    def length(self):
        a = self.p2[0] - self.p1[0]
        b = self.p2[1] - self.p1[1]
        l = math.hypot(a,b)
        return l
    ##
    # @brief Get the slope of the segment
    #
    # @todo Test for /0
    #
    def m(self):
        a = self.p2[1] - self.p1[1]
        b = self.p2[0] - self.p1[0]
        m = a / b
        return m
    ##
    # @brief Get the y intercept of the segment
    #
    # @return y intercept
    def b(self):
        b = self.p2[1] - (self.m() * self.p2[0])
        return b
    ##
    # @brief Get the coordinate of a point some distance from point 1
    #
    # @param d distance from point 1
    # @return x and y coordinate
    def d_line(self, d):
        if d == 0:
            return self.p1[0], self.p1[1]
        if d > self.length():
            return None
        m = self.m()
        u = d / math.sqrt(m * m + 1)
        if self.p2[0] < self.p1[0]:
            u = -u
        v = (m * d) / math.sqrt(m * m + 1)
        if self.p2[1] < self.p1[1]:
            v = -v
        return (self.p1[0] + u, self.p1[1] + v)
##
# Class representing a polyline composed of segments
#
# @see segment
#
class poly_line:
    segments = list()
    def __init__(self, points):
        for i in range(len(points) - 1):
            s = segment(points, points[i+1])
            self.segments.append(s)
    ##
    # @brief Get the length of all the segments
    #
    # @return total length
    #
    def length(self):
        l = 0.0
        for s in self.segments:
            l += s.length()
        return l
    ##
    # @brief Get the coordinate of a point some distance from the polyline 
    # origin
    #
    # @param d distance from the origin
    # @return x and y coordinate for the point
    #
    def d_line(self, d):
        if d == 0:
            return self.segments[0].d_line(0)
        l = 0.0
        for s in self.segments:
            if l >= d:
                break
            dd = d - l
            l += s.length()
            seg = s
        return seg.d_line(dd)

I just ran it and wrote the attached output, but no guarentees it will work for anyone. I think that the service ESRI wrote used lower level raster access, because of it's speed. I got comparable speed with the gdal version.

Sorry if I didn't help, but it is possible. I stopped working on it, but it may help others. I will watch the thread and try to answer questions.

kss

Re: how to create own surface profile service?

KarenRobine — Tue, 08 May 2012 17:00:11 GMT

I just found this:

http://www.forumsig.org/archive/index.php/t-28565.html

related to the initial question of using your own data. I haven't delved in yet, but will be in the future...

Karen

Re: how to create own surface profile service?

KevinHibma — Tue, 08 May 2012 17:41:31 GMT

Im only seeing this thread now - there are a lot of solutions to this problem.

Of course everyone has seen the SOE which is available as a sample from the Applications Prototype lab:
Configurable SOE: http://www.arcgis.com/home/item.html?id=3b65b141041044ca97d83eb8264da08c
Basic Viewer to go with the SOE: http://www.arcgis.com/apps/OnePane/basicviewer/profile.html

kss has a nice script to perform this analysis above. (also using matplotlib to draw the output). This should work real well at 10.0 but be sure to install matplotlib on your server machine. At 10.1 matplotlib comes with ArcGIS (Desktop and Server)

New at 10.1 is a tool "Stack Profile". This tool produces a table of X,Y,Z and optional graph output. Please note that the graphing component is not supported in Server, so you can use the tool to produce your table of X,Y,Z.
http://resourcesbeta.arcgis.com/en/help/main/10.1/index.html#//00q9000000mm000000
I've attached a HTML file which consumes the "Stack Profile" service.
Basically it allows you to draw a line on a map, and looks for the result table. It builds up some arrays and plots the chart. (Makes use of JavaScript API 2.7)
Make use to update the references to your server.

Note, the attached web page and the one referenced above are pretty similar to each other and one could extend/modify either page to suit individual needs.
The difference is how you "draw" the profile : does the service create an image, or does the client render an image with X,Y,Z information?

Re: how to create own surface profile service?

KarenRobine — Tue, 08 May 2012 18:32:34 GMT

Hi Kevin:
Thanks so much for the response. That 10.1 tool looks interesting.

But back to the question at hand:
We are looking for a way to replace the data utilized in the SOE, or the ProfileService (at http://samplesserver2.arcgisonline.com/arcgis/services Elevation/ESRI_Elevation_World)
I also checked out:
http://sampleserver4.arcgisonline.com/ArcGIS/rest/services/Elevation/ESRI_Elevation_World/MapServer
and it looks like the data it uses for this is a Global Elevations layer, which doesn't look to be quite precise enough for my purposes.
The SOE that's at
http://resources.arcgis.com/gallery/file/arcobjects-net-api/details?entryID=87BEC705-1422-2418-34B5-308930DE323A

may handle this, but the source provided is in Java (and I haven't researched the code enough yet to see if it supports what I'm asking, or if its only a client app). So am wondering if that's the only solution for the problem? (ps. would be nice to get the model associated with the ProfileService?)..

Thanks again....

Re: how to create own surface profile service?

KevinHibma — Tue, 08 May 2012 18:43:30 GMT

Karen,

In the instructions to the SOE setup it explains you attach the SOE to a map service with a single band raster layer which will act as your elevation service.
So yes, you can and must use your own data to set it up. The SOE gets tied to that map service.

SETUP
1) Create a Map Document (mxd) with at least one elevation layer; the elevation layer must be a single band raster layer. Other types of layers are ok to use because the SOE will only expose a resource for the single band raster layers.

2) Publish the Map Document as a Map Service and then stop the map service after it has been successfully published.

3) Right-click on the map service and select �??Service Properties�??. In the �??Map Service Properties�?? dialog, select the �??Capabilities�?? tab.

4) From the list of capabilities, check on the new �??Get Elevations�?? capability and click the OK button.

5) Start the map service.

6) Go to the Services Directory admin page and clear the REST cache.

7) Navigate to the Services Directory page of the Map Service and you�??ll now see the Elevations SOE listed at the bottom of the page in the �??SUPPORTED EXTENSIONS�?? section. Follow the available links to view the available REST resources and operations for the Elevations SOE.

EDIT:
Note, the SOE is custom code: there is no "model" which produces the result. The code (I could only imagine) would be something similar to the code in the above Python script.

Re: how to create own surface profile service?

KarenRobine — Tue, 08 May 2012 18:48:56 GMT

Thanks Kevin.. That's perfect.. Should have read those direction first before responding (oops)...

Now I get to see if the Elevation layer can be a WCS type of layer (which is what ours will be)..

Re: how to create own surface profile service?

KarenRobine — Tue, 08 May 2012 20:02:27 GMT

Also, is there a .NET version of that same SOE out there anywhere? The only source for the SOE is a Java-based version. All the rest of the stuff in that .zip file are clients.

Re: how to create own surface profile service?

DrewDowling — Tue, 08 May 2012 21:05:27 GMT

I would like to second Karen's request for .NET source code for the elevation profile SOE if it exists.

Re: how to create own surface profile service?

nicogis — Wed, 06 Feb 2013 14:47:16 GMT

news: esri is coming out with the new version x 10.1 and the included source net and java
here you can see my preview beta http://sit.sistemigis.it/Samples/Elevations/
however you download esri version ...

Re: how to create own surface profile service?

nicogis — Tue, 09 Apr 2013 16:15:22 GMT

news: esri is coming out with the new version x 10.1 and the included source net and java
here you can see my preview beta http://sit.sistemigis.it/Samples/Elevations/
however you download esri version ...

Now you can download new version 10.1 soe elevations:
http://www.arcgis.com/home/item.html?id=60bf17f2d59441839135adc25654a7d9

Re: how to create own surface profile service?

PranaviGaddipati — Tue, 27 Aug 2013 05:08:44 GMT

Hi.
I'm working with ArcGIS 10.1 and I'm using the stack profile tool in the python script. It is working perfectly fine but the sampling distance used to get the values is by default the raster pixel size. For me, it is 25.

But I want it to be a configurable parameter while using this stack profile tool. The python script for this tool doesn't have an input option as such for the sampling distance. The current inputs I'm using are the polyline i'm drawing and the raster.

Please let me know if there is any alternative for this through the python script?