<\/HEAD>
Ingrese la contraseña:b7b7b7b7b7b7b7b7
Conectado exitosamente a ArcGIS Enterprise A
Buscar en GIS los datos de entrada para el análisis
Índice Modificado por Humanos
En [2]:
# Buscar el elemento de capa de imágenes del Índice Modificado por Humanos por título
item_hmi = web_gis.content.search('title:Human Modified Index', 'Imagery Layer')[0]
item_hmi
">
">
Salida[2]:
Índice Modificado por Humanos Una medida del grado de modificación humana, el índice varía de 0.0 para una condición de paisaje virgen a 1.0 para las áreas más modificadas.

Capa de Imágenes por djohnsonRA
Última modificación: 06 de julio de 2017
0 comentarios, 2 vistas
En [3]:
# Buscar el elemento de capa de imágenes DEM por título
item_dem = web_gis.content.search('title:USGS NED 30m', 'Imagery Layer')[0]
item_dem
Salida[3]:
Delimitación y extensión del área de estudio
En [4]:
font-weight: bold; font-size: 25.998px; margin: 0.538em 0px 0px;" id="toc-hId-1356367382" id="toc-hId-1356367382">Mostrar los datos de análisis
In [9]:
# Obtener una referencia a la capa de imágenes del portal item
lyr_hmi = item_hmi.layers[0]
# Establecer la extensión de la capa a la extensión geográfica del área de estudio y mostrar los datos.
lyr_hmi.extent = extent_studyarealyr_hmi
In [10]:
# Obtener una referencia a la capa de imágenes del portal item
lyr_dem = item_dem.layers[# Extraer los datos del Índice Modificado por Humanos dentro de la geometría del área de estudio
hmi_clipped = clip(raster=lyr_hmi, geometry=geom_studyarea)
hmi_clipped
In [14]:
# Extraer los datos de Elevación dentro de la geometría del área de estudio
elev_clipped = clip(raster=lyr_dem, geometry=geom_studyarea)
elev_clipped
In [15]:<2Div3e<2Div3e<2Div3e<2Pre3e<2Span3e# Extraer los datos de Pendiente dentro de la geometría del área de estudio
<2Span3eslope_clipped<2Span3e <2Span3e=<2Span3e <2Span3eclip<2Span3e(<2Span3eraster<2Span3e=<2Span3elyr_slope<2Span3e,<2Span3e geometry<2Span3e=<2Span3egeom_studyarea<2Span3e)
<2Span3e# Aplicar la función Stretch para mejorar la visualización de la capa slope_clipped.
<2Span3eslope_clipped_stretch<2Span3e <2Span3e=<2Span3estretch<2Span3e(<2Span3eraster<2Span3e=<2Span3eslope_clipped<2Span3e,<2Span3estretch_type<2Span3e=<2Span3e'StdDev'<2Span3e,<
dra<2Span3e=<2Span3e'true'<2Span3e)
<2Span3eslope_clipped_stretch<2Span3e
<2FPre>
<2FDiv>
<2FDiv>
<2FDiv>
<2FDiv>
< div >< div
style=color:#000;padding:.5 em .5 em .5 em .4 em >
Paso 1 : NormalizaciónUtilice la función Remap para normalizar cada conjunto de datos de entrada a una escala común de 1 - 9 , donde 1 = menos adecuado y 9 = más adecuado.
< div >< div
stile=color:#303f9f;border-top :solid transparent
1 px ;padding :.4 em ; >En [16]: < / div > < div > < div stile = "background:#f7f7f7;border :solid #cfcfcf
1 px ; " > < div stile = "background :transparent;border:none;margin :.4 em ; " > < pre clase = "lia-code-sample line-numbers language-none " > < span stile = "color:#408080 ; " ># Crear un mapa de colores para mostrar los resultados del análisis con nueve colores que van desde rojo hasta amarillo y verde.
</ span > clrmap < span stile = "color:#666666 ; " > = </ span > < span > [[ < span stile = "color:#666666 ; " >1 </ span > , < span stile = "color:#666666 ; " >230 </ span > , < span stile = "color:#666666 ; " >0 </ span > , < span stile = "color:#666666 ; " >0 </ span > ] , [ < span stile = "color:#666666 ; " >2 </ span > , < span stile = "color:#666666 ; " >242 </ span > , < span stile = "color:#666666 ; " >85 </ span > , < span stile = "color:#666666 ; " >0 </ span > ] , [ < span stile = "color:#666666 ; " >3 </ span > , < span stile = "color:#666666 ; " >250 </ span > , < span stile = "color:#666666 ; " >142 </ span > , < span stile = "color:#666666 ; " >0 </ span > ] , [ < span stile = "color:#666666 ; " >4 </ span > , < span stile = "color:#666666 ; " >255 </ span > , < span stile = "color:#666666 ; " >195 </ span > , < span stile = "color:#666666 ; " >0 </ span > ] ,
[
5,
255,
255,
0], [
6,
197,
219,
0], [
139<\/SPAN>,<\/SPAN> 181<\/SPAN>,<\/SPAN> 0<\/SPAN>], [8<\/SPAN>,<\/SPAN> 86<\/SPAN>,<\/SPAN> 148<\/SPAN>,<\/SPAN> 0<\/SPAN>],<\/SPAN>
<\/SPAN><\/SPAN><\/SPAN><\/SPAN><\/SPAN><\/SPAN><\/SPAN> [<\/SPAN>9<\/SPAN>,<\/SPAN> 38<\/SPAN>,<\/SPAN> 115<\/SPAN>,<\/SPAN> 0<\/SPAN>]<\/SPAN>]<\/SPAN><\/PRE><\/DIV><\/DIV><\/DIV><\/DIV><\/DIV>En [17]:<\/DIV># Normalizar los datos de elevación<\/SPAN>
elev_normalized = remap(raster=elev_clipped,
input_ranges=[0,490, 490,980, 980,1470, 1470,1960, 1960,2450,
2450,2940, 2940,3430, 3430,3700, 3920,4100],
output_values=[9,8,7,6,5,4,3,2,1], astype='U8')
# Mostrar imagen con mapa de colores de los datos de elevación reclasificados<\/SPAN>
colormap(elev_normalized, colormap=clrmap)&nbps;<\/PRE><\/DIV><\/DIV><\/DIV><\/DIV>Salida[17]:<\/DIV>
<\/DIV><\/DIV><\/DIV><\/DIV><\/DIV>En [18]:<\/DIV># Normalizar los datos de pendiente
<\/SPAN>slope_normalized<\/>slope_normalized = remap(raster=slope_clipped,&nbps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsps;nbsp;sinput_ranges=[0,1,1,2,2,3,3,5,5,7,7,9,9,12,12,15,
                              15,100],
output_values=[9,8,7,6,5,4,3,2,1], astype='U8')
# Mostrar una imagen con mapa de colores de los datos de pendiente reclasificados
colormap(slope_normalized,colormap=clrmap)<\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/><\/> class="">],<\/SPAN>
output_values<\/SPAN>=<\/SPAN>[<\/SPAN>9<\/SPAN>,<\/SPAN>8<\/SPAN>,<\/SPAN>7<\/SPAN>,<\/SPAN>6<\/SPAN>,<\/SPAN>5<\/SPAN>,<\/SPAN>4<\/SPAN>,<\/SPAN>3<\/SPAN>,<\/SPAN>2<\/SPAN>,<\/SPAN>1<\/SPAN>],<\/SPAN> astype<\/SPAN>=<\/SPAN>'U8'<\/SPAN>)
<\/SPAN># Mostrar una imagen con mapa de colores de los datos HMI reclasificados
<\/SPAN>colormap<\/SPAN>(<\/SPAN>hmi_normalized<\/SPAN>,<\/SPAN> colormap<\/SPAN>=<\/SPAN>clrmap<\/SPAN>)<\/SPAN><\/PRE><\/DIV><\/DIV><\/DIV><\/DIV>Out[19]:<\/DIV>
<\/DIV><\/DIV><\/DIV><\/DIV><\/DIV>Paso 2: Ponderación<\/H3>Utilice el operador de multiplicación sobrecargado * para asignar un peso a cada conjunto de datos normalizado según su importancia relativa para el resultado final.<\/P><\/DIV><\/DIV><\/DIV>In [20]:<\/DIV># Aplicar pesos a los datos normalizados usando el operador de multiplicación sobrecargado
# "*".
<\/SPAN># - Índice Modificado Humano (Human Modified Index): 60%
<\/SPAN># - Pendiente (Slope): 25%
<\/SPAN># - Elevación (Elevation): 15%
<\/SPAN>hmi_weighted<\/SPAN> =<\/> hmi_normalized * 0.6
slope_weighted = slope_normalized * 0.25
elev_weighted = elev_normalized * 0.15Paso 3 : Suma
Sume los conjuntos de datos ponderados para producir un resultado final del análisis.< /p >
# Calcular la suma de los conjuntos de datos ponderados usando el operador de suma sobrecargado "+"
result_dynamic = colormap(hmi_weighted + slope_weighted + elev_weighted, colormap = clrmap, astype = 'U8')
result_dynamic
< divclass = " "style = " color:#d84315;padding :0.4em " > Out[21] : < / div > < divclass = " "style = " padding :0.4em " > < img__jive_id = "360677 "class = " image-13 jive-image "src = "
https://us.v-cdn.net/6038851/uploads/legacyfs/online/360677_pastedImage_13.jpg "style = " border :0px ; " /> < / div > < / div > < / div > < / div > < / div > < divclass = " "style = " color:#000000;background:-webkit-linear-gradient(left,transparent -40px,transparent 1px,transparent 1px,transparent 100%);border-width :1px;border-style :solid;border-color :transparent;padding :5px " > < div > < divclass = " "style = " color:#000000;padding :.5em .5em .5em .4em " > < h3id = "toc-hId--312933377 "id = "toc-hId--312933377 "style = " font-weight:bold;font-size :18.004px;margin :.777em 0px 0px;color:#000000;padding :.5em .5em .5em .4em " > El mismo análisis también se puede realizar en una sola operación </ h3 > </ div > </ div > </ div > < divclass = " "style = " color:#000000;background:-webkit-linear-gradient(left,transparent -40px,transparent 1px,transparent 1px,transparent 100%);border-width :1px;border-style :solid;border-color :transparent;padding :5px " > < div > < divclass = " "style = " color:#303f9f;border-top-width :1 px;border-top-style :solid;border-top-color :transparent;padding :.4 em ;font-family:"Courier New", Courier, monospace;font-size:.875rem;background-color:#f7f7f7;border-width :1 px;border-style :solid;border-color :#cfcfcf;margin :.4 em ;white-space :pre-wrap;text-align:left;color:#303f9f;font-weight:bold;font-style:none;text-decoration:none;text-transform:none;font-variant:none;font-feature-settings:none;-webkit-font-smoothing:auto;-moz-osx-font-smoothing:auto;-webkit-text-stroke-width:.01 px;-webkit-text-stroke-color:#303f9f;-webkit-text-fill-color:#303f9f;-moz-text-fill-color:#303f9f;-moz-text-stroke-width:.01 px;-moz-text-stroke-color:#303f9f;-moz-text-fill-color:#303f9f;-moz-text-stroke-width:.01 px;-moz-text-stroke-color:#303f9f;-moz-text-fill-color:#303f9f;-moz-text-stroke-width:.01 px;-moz-text-stroke-color:#303f9f;-moz-text-fill-color:#303f9f;-moz-text-stroke-width:.01 px;-moz-text-stroke-color:#303f9f;-moz-text-fill-color:#303f9f;-moz-text-stroke-width:.01 px;-moz-text-stroke-color:#303f9f;-moz-text-fill-color:#303f9f;-webkit-font-smoothing:auto;-moz-osx-font-smoothing:auto;font-family:"Courier New", Courier, monospace;font-size:.875rem;background-color:#fff;border-radius:.25rem;margin-bottom:.75rem;margin-top:.75rem;padding-left:.75rem;padding-right:.75rem;padding-top:.75rem;padding-bottom:.75rem;text-align:left;color:#333333;font-weight:bold;font-style:none;text-decoration:none;text-transform:none;font-variant:none;font-feature-settings:none;text-shadow:none;text-indent:none;text-overflow:none;text-align-last:start;text-anchor:start;text-decoration-line:none;text-decoration-style:none;text-decoration-thickness:auto;text-emphasis-style:none;text-emphasis-position:none;text-emphasis-color:none;text-combine-upright:none;text-orientation:mixed;text-underline-position:auto;text-underline-offset:auto;text-underline-thickness:auto;text-shadow:none;text-indent:none;text-overflow:none;text-align-last:start;text-anchor:start;text-decoration-line:none;text-decoration-style:none;text-decoration-thickness:auto;text-emphasis-style:none;text-emphasis-position:none;text-emphasis-color:none,text-combine-upright:none,text-orientation:mixed,text-underline-position:auto,text-underline-offset:auto,text-underline-thickness:auto,text-shadow:none,text-indent:none,text-overflow:none,text-align-last:start,text-anchor:start,text-decoration-line:none,text-decoration-style:none,text-decoration-thickness:auto,text-emphasis-style:none,text-emphasis-position:none,text-emphasis-color:none,text-combine-upright:none,text-orientation:mixed,text-underline-position:auto,text-underline-offset:auto,text-underline-thickness:auto">
result_dynamic_one_op = colormap(
raster=raster(
# Capa del índice modificado humano (Human modified index layer)
0.60 * remap(raster=clip(raster=lyr_hmi, geometry=geom_studyarea),
input_ranges=[0.0,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.5,
0.5,0.6,0.6,0.7,0.7 class="">,
0.8, 0.8,1.1], output_values=[9,8,7,6,5,4,3,2,1]) + # Capa de pendiente 0.25 * remap(raster=clip(raster=lyr_slope,
geometry=geom_studyarea),
input_ranges=
[0, 1, 1, 2, 2, 3, 3, 5, 5, 7, 7, 9, 9, 12, 12, 15,
15, 100],& nbsp; output_values=[9, 8, 7, 6, 5, 4, 3, 2, 1])& nbsp; +& nbsp;
# Capa de elevación 0.15 * remap(raster=clip(raster=lyr_dem , geometry=geom_studyarea),
input_ranges=[-90 ,250 ,250 ,500 ,500 ,750 ,750 ,1000 ,1000 ,1500 ,1500 ,2000 ,2000 ,2500 ,2500 ,3000 ,3000 ,5000],
output_values=[9 ,8 ,7 ,6 ,5 ,4 ,3 ,2 ,1])& nbsp;),& nbsp;
colormap=
clrmap , astype=
'U8')
result_dynamic_one_op
< div clas s= " "> < div clas s= " "> < div clas s= " "> < div clas s= " "> < div clas s= " "> En [25]: </ div > < div clas s= " "> < div clas s= " f7f7f7 ; border :1px solid cfcfcf ; "> < div clas s= " "> < pre clas s= " lia-code-sample line-numbers language-none "> < span styl e=