Hello everyone, am trying to calculate the flooded areas of a coastal area, I am using sea level rise and storm surge. The storm surges given are small stretches of the coastal line and I am stuck on how to go about this using arcpy. Below is the table of my storm surge with S10 meaning 10 years and S100 meaning 100 years with each row representing a stretch on the coastal area. I tried to clip each row with the input raster data but i am getting errors. fc = arcpy . GetParameterAsText ( 2 ) inputraster = arcpy . GetParameterAsText ( 3 ) with arcpy . da . SearchCursor ( fc , [ "SHAPE@" ] ) as cursor1 : for rows in cursor1 : geom = rows [ 0 ] del cursor1 cursor = arcpy . da . SearchCursor ( fc , [ 'S10' , 'S100' ] ) row = cursor . next ( ) while row : while row [ 0 ] : outputclip = outputfolder + str ( row [ 0 ] ) + "_C" clipgeo = geom extent = str ( clipgeo . Extent ) arcpy . Clip_management ( inputraster , extent , outputclip , " #","ClippingGeometry","MAINTAIN_EXTENT") Any help is really appreciated. Thanks. izmircoastdiva1.lpk.zip JenniferAcunto_2-1759942506381.png

I tried following all the steps but I am getting incorrect output . Also i can not get around how to write a python script for the process

Thanks, let me try working with this first and see what happens.

You may be interested in this post too: https://community.esri.com/people/GTaylor-esristaff/blog/2017/06/16/3d-flood-tools

nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

How to calculate flooded areas for the shore line

stashavint

Hello everyone, am trying to calculate the flooded areas of a coastal area, I am using sea level rise and storm surge. The storm surges given are small stretches of the coastal line and I am stuck on how to go about this using arcpy. Below is the table of my storm surge with S10 meaning 10 years and S100 meaning 100 years with each row representing a stretch on the coastal area.

I tried to clip each row with the input raster data but i am getting errors.

fc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">)</SPAN>
inputraster <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> cursor1<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> rows <SPAN class="keyword token">in</SPAN> cursor1<SPAN class="punctuation token">:</SPAN> 
&nbsp;&nbsp;&nbsp;&nbsp; geom <SPAN class="operator token">=</SPAN> rows<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> 
<SPAN class="keyword token">del</SPAN> cursor1
cursor <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'S10'</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">'S100'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
row <SPAN class="operator token">=</SPAN> cursor<SPAN class="punctuation token">.</SPAN>next<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">while</SPAN> row<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp; <SPAN class="keyword token">while</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; outputclip <SPAN class="operator token">=</SPAN> outputfolder <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">+</SPAN> <SPAN class="string token">"_C"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; clipgeo <SPAN class="operator token">=</SPAN> geom
&nbsp;&nbsp;&nbsp;&nbsp; extent <SPAN class="operator token">=</SPAN> str<SPAN class="punctuation token">(</SPAN>clipgeo<SPAN class="punctuation token">.</SPAN>Extent<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Clip_management<SPAN class="punctuation token">(</SPAN>inputraster<SPAN class="punctuation token">,</SPAN>extent<SPAN class="punctuation token">,</SPAN>outputclip<SPAN class="punctuation token">,</SPAN>"<SPAN class="comment token">#","ClippingGeometry","MAINTAIN_EXTENT")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Any help is really appreciated. Thanks.

izmircoastdiva1.lpk.zip

JenniferAcunto_2-1759942506381.png

Finde mehr Beiträge mit dem Tag

Accepted answers

All comments

DanPatterson_Retired

Code Formatting... the basics++ to format your code so syntax and structure can be checked

DanPatterson_Retired

Ok... what are the errors. since the first 7 lines of code essentially do nothing and line 5 is indented improperly and even if it wasn't, you delete the cursor1 object so it is gone. and where is getparameterastext numbers 0 and 1?

stashavint

that was just a part of the code i was writing, the error i get says :

ExecuteError<SPAN class="punctuation token">:</SPAN> ERROR <SPAN class="number token">000622</SPAN><SPAN class="punctuation token">:</SPAN> Failed to execute <SPAN class="punctuation token">(</SPAN>Clip<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN> Parameters are <SPAN class="operator token">not</SPAN> valid<SPAN class="punctuation token">.</SPAN>
ERROR <SPAN class="number token">000628</SPAN><SPAN class="punctuation token">:</SPAN> Cannot set input into parameter rectangle<SPAN class="punctuation token">.</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

repost a fixed up version of your code. The error states that one of the parameters is incorrect

http://desktop.arcgis.com/en/arcmap/latest/tools/analysis-toolbox/clip.htm

and the code in its current state makes it hard to figure out what specifically this is without guessing

XanderBakker

I notice that you are clipping a raster and provide a string version of the rectangle. See the code below:

<SPAN class="keyword token">import</SPAN> arcpy
ext <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Extent<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1952602</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">294196</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1953546</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">296176</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> str<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">print</SPAN> <SPAN class="string token">" "</SPAN><SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>str<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">" "</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">4</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>

extent <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{0} {1} {2} {3}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMin<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMin<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>XMax<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMax<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> extent
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

The first print statement will give you the extent as you generate it (simply converting the extent object to string), the second what you could do to get a valid syntax of the extent and the third a more readable version.

This will yield:

1952602 294196 1953546 296176 NaN NaN NaN NaN
1952602 294196 1953546 296176
1952602.0 294196.0 1953546.0 296176.0<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

To be able to see what is going wrong, you should provide the entire code since there are several things missing (parameters 0 and 1, how is the path defined, etc) that could be causing errors. And if it is not the code, it might be data related. So if possible provide the data too.

Nimm an der Diskussion teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.