Revisar los datos de capacidad de desalojo

[edithvilla]

Procesar la base de datos de la red de desalojo: \MEGADAPT\mega_insumos\mayo2018\ykim\Sewee_network_october2017\Tables for Contingency Matrix\Final Data\Infraestructure_Decision_Criteria_AGEB.xlsx.

Infraestructure_Decision_Criteria_AGEB.xlsx. Hojas: dry_AHEB_DF y wet_AHEB_DF Actividades:

• [x] Promediar los datos que vienen en las dos hojas las columnas: E, F, G y H

• [x] Normalizar los datos de la columnas E, F, G y H

• Los datos de la columna F "capacity" y los datos de la columna I "weghited_total_degree_nm" son los datos que se van a ocupar en el modelo bayesiano

Quedo en espera de la indicación en qué carpeta se respaldará la información @serranoycandela

[edithvilla]

• [x] Promediar los datos que vienen en las dos hojas las columnas: E, F, G y H

Se han promediado los datos de las columnas E, F, G y H. Los datos están en shape: infraestructura_prom.shp

@serranoycandela en donde se están colocando los resultados?

[edithvilla]

Acuerdos de la reunión 9 de agosto, 2018 con Luis

• [x] 1. Revisar en el informe de Infraestructura de la red de desalojo (Yeowon) ¿cómo se hizo la intersección de la red con las AGEB?

• [x] 2. Documentar todo el procedimiento que se va hacer con la información de la topología de la red

• [x] 3. Identificar los NA verdaderos a partir de la red de desalojo simplificada (lluvias y secas)

• [x] 4. Identificar los NA que se puedan corregir a partir de la red de desalojo primaria que nos proporcionó sacmex (lluvias y secas)

• [x] 5. Identificar los NA que se encuentran en el centro de la ciudad que se tienen que corregir FALTA definir qué valor se le asignará (promedio de vecinos) (lluvias y secas)

[edithvilla]

1. Revisar el el informe de Infraestructura de la red de desalojo (Yeowon):

Para la capacidad ocupó la ecuación de Manning (m3/s) Q = Flow Rate (m3/s) A = Flow Area (m2), For the case of no free surface = pipe cross-sectional area kn = 1.0 n = Manning’s Roughness Coefficient, For concrete channel = 0.011 Rh = Hydraulic Radius (m) S = channel slope (m/m)

Densidad = El volumen del drenaje dividido por el área del AGEB = m3/m2

[edithvilla]

1. Revisar el informe de Infraestructura de la red de desalojo. Al no encontrar información sobre cómo se obtuvo la capacidad de la red de drenaje y cómo pasó los datos a nivel de AGEB; se le mandó un correo a Yeowon con copia a los coordinadores del proyecto, este es el mensaje:

We are trying to use your data of sewer drainage network and we have a couple of questions, we imagine that you first did the calculations on the tubes of the sewer and then you calculated the values for the AGEB polygons.

1. Do you have the intermediate files where you calculated the capacity for the sewer tubes?

2. What procedure did you follow to pass the values from the tubes to the AGEB polygons?

estamos en espera de su respuesta, mientras tanto se seguirá trabajando en corregir los datos con valores N/A

[edithvilla]

está ligado al issue #36

[edithvilla]

Se calculó la capacidad del drenaje con la ecuación de Manning @serranoycandela los nuevos cálculos se encuentran en el archivo:

capacidad_drenaje_dry.csv

capacidad_drenaje_wet.csv

en:

\MEGADAPT\mega_procesamiento\banco_datos\capacidad_drenaje\bd

[edithvilla]

En la reunión del 5 de agosto con Luis, se acordó que tengo que entregar 5 capas:

• [x] 1. Volumen de agua por AGEB

\SIG\desarrollo\megadapt\procesamiento\volumen_drenaje_ageb\volumen_cdmx_ageb.shp

• [x] 2. Número de alcantarillas por AGEB, @aleestrada nos proporcionará la capa de número de coladeras por AGEB con sus respectivos metadatos.

\SIG\desarrollo\megadapt\procesamiento\alcantarillado\alcantarillado_ageb.shp

• [ ] 3. Capacidad de desalojo en la red primaria y secundaria (Q50 y Q100)

• [x] 4. Capacidad de desalojo del drenaje profundo, este calculo es para las AGEB donde hay lumbreras)

[edithvilla]

Se necesita corregir:

• [x] 1. La columna de la pendiente (columna AB), falta multiplicar el promedio del 2 al millar por la longitud del segmento
• [x] 2. La columna de "slp_yeo" se necesita convertir a números absolutos, así que se crea la columna "AA"
• [x] 3. Se necesita corregir las fórmulas de las Q de yeo para poder indicar la nueva columna "slp_yeo_abs"
• [x] 4. Calcular el volumen del segmento

[edithvilla]

@aleestrada puedes colocar en esta carpeta la capa de alcantarillado de la CDMX

\MEGADAPT\mega_insumos\septiembre2018\aestrada

y por este medio colocas la fuente de información

[aleestrada]

@edithvilla ya cargué los archivos del alcantarillado de la CDMX, los datos se obtuvieron en el 2016.

Los datos e obtuvieron de la página de Datos Abiertos de la Ciudad de México, en ella indican que la fuente de los datos es la Secretaría de Obras y Servicios.

http://www.datosabiertos.cdmx.gob.mx:5000/dataset/alcantarillado-en-la-ciudad-de-mexico

[edithvilla]

El archivo dbf del shape de Alcantarillado_DF_2015.shp viene mal: Se corrigió y solo se quedó un los atributos id (se le asignó un número consecutivo) y se calculó la columna "x" y "y"

Se guardó la información original como csv: "alcantarillado_cdmx_2015.csv" en \SIG\desarrollo\megadapt\procesamiento\alcantarillado\insumos

\MEGADAPT\mega_insumos\septiembre2018\aestrada

[edithvilla]

• [x] 4. Capacidad de desalojo del drenaje profundo, este calculo es para las AGEB donde hay lumbreras)

• [x] Capa de ageb con datos de lumbreras

[edithvilla]

@serranoycandela

Capa de capacidad de desalojo del drenaje profundo (Q100 y Q50) \SIG\desarrollo\megadapt\procesamiento\drenaje_profundo\capacidad_drenaje_profundo.shp

Capa de lumbreras por ageb \SIG\desarrollo\megadapt\procesamiento\drenaje_profundo\lumbreras_ageb.shp

[edithvilla]

Calcular el número y densidad de alcantarillas por AGEB

Responsable: Edith Villa

Fecha: 10 de septiembre de 2018

Capas de insumo: Alcantarillado_DF_2015.shp (EPSG: 4326) ageb_utm.shp (EPSG: 32614)

Ubicación de trabajo: C:\Dropbox (LANCIS)\SIG\desarrollo\megadapt\procesamiento\alcantarillado

Procesamiento:

1. Se cambió la proyección de la capa de Alcantarillado_DF_2015 (geográficas) a alcantarillado_utm.shp (proyección utm)

2. Se identificaron errores en la tabla

\SIG\desarrollo\megadapt\procesamiento\alcantarillado)

3. Los errores se corrigieron en QGIS, se crearon tres columnas en la tabla de atributos con los siguientes identificadores: "id", "x" y "y"

4. Para la columna "id" se le agregó un número consecutivo, para las columnas "x" y "y" se calculó mediante la Calculadora de campos/Geometría/$x y Calculadora de campos/Geometría/$y

5. Se eliminaron las tres columnas de origen que presentaron errores.

6. En QGIS mediante la herramienta Vectorial/Herramientas de análisis/Contar puntos de polígonos, se creó una capa nueva: alcantarillado_ageb.shp en esta capa se encuentran los campos ageb_id y num_alcant, en este segundo campo están el número de alcantarillas que se encuentran por AGEB.

[edithvilla]

Agregar a la capa de lumbreras el identificador del AGEB que le corresponde

Responsable: Edith Villa

Fecha: 11 de septiembre de 2018

Capas de insumo: lumbreras.shp ageb_cdmx.shp

Ubicación de trabajo: C:\Dropbox (LANCIS)\SIG\desarrollo\megadapt\procesamiento\drenaje_profundo

Procedimiento en QGIS:

Vectoria/Herramientas de gestión de datos/Unir atributos por localización

Capa vectorial objetivo: lumbreras.shp

Unir capa vectorial: ageb_cdmx.shp

Seleccionar: dentro

Estadísticas para resúmen: max

Se crea una nueva capa: lumbreras_ageb.shp a la cual se le agregó la columna del atributo "ageb_id"

Nota: se hicieron un par de correcciones en el campo de "MUID" con los identificadores de las lumbreras. Se corrigió RR_L1 y RR_L2 por IRRL1 y IRRL2, ya que corresponde a las lumbreras del Interceptor Río de los Remedios.

[edithvilla]

Cálculo de volumen de agua (m3) en las estructuras de desalojo por AGEB

Objetivo: calcular el volumen de agua en las estructuras de la red primaria y secundaria por AGEB.

Responsable: Edith Villa Fecha: 6 de septiembre de 2018 Ubicación de trabajo: C:\Dropbox (LANCIS)\SIG\desarrollo\megadapt\procesamiento\volumen_drenaje_ageb

1. Insumos

   red_desalojo_secas.shp (esta es una capa de líneas y contiene los nodos iniciales y finales del segmento, diámetro de las estructuras y nombre de las estructuras)

   ageb_cdmx.shp (esta es una capa de polígono, en su tabla de atributos solo cuenta con el campo identificador del ageb: ageb_id)

2. En el sistema de información geográfica QGIS se hicieron las siguientes operaciones:

Con la herramienta de Intersección en el menú Vectorial/Herramientas de geoproceso/ se hizo la intersección de la capa de la red_desalojo_secas y ageb_cdmx para poder obtener solo los tramos de la red que se encuentran dentro de cada AGEB. Se obtuvo la capa: tp_tramos_ageb.shp. Con la operación de la intersección le crea un campo nuevo: ageb_id para poder colocarle a cada tramo el identificador del ageb en donde se encuentra.

Sobre la capa de tp_tramos_ageb.shp se calculó:

• length, se calculó mediante la herramienta Vectorial/Herramientas de geometría/Exportar / Añadir columnas de geometría y se selecciona la capa de tp_tramos_ageb.shp y se creo una nueva capa: tp_tramos_ageb_v2.shp

  Sobre la capa tp_tramos_ageb_v2.shp se crearon los siguientes campos:

• diam_metro, se creó un campo nuevo de número decimal (real) para calcular el diametro en metros para cada tramo con la calculadora de QGIS: (DIAMETER/100)

• volumen, se creó un campo nuevo de número decimal (real) para calcular el volumen para cada tramo con la calculadora de QGIS: (pi() ("diam_metro" ^ 2)/4) "length" .

Mediante un código de python (volumen.py) se hizo el cálculo de la suma de los volumenes por cada tramo para cada AGEB y se le agregó esta columan de vol_dren al shape de ageb_cdmx.shp. Esta capa se renombró volumen_cdmx_ageb.shp

Codigo para obtener la suma de los volumenes de cada tramo de cada AGEB

Se abren capas vectoriales

tramos = QgsVectorLayer("C:/Dropbox (LANCIS)/SIG/desarrollo/megadapt/procesamiento/volumen_drenaje_ageb/tp_tramos_ageb_v2.shp", "tramos", "ogr") agebs = QgsVectorLayer("C:/Dropbox (LANCIS)/SIG/desarrollo/megadapt/procesamiento/volumen_drenaje_ageb/ageb_cdmx.shp", "ageb", "ogr")

Identifica todos los tramos para cada AGEB y hace la suma del volumen del conjunto de tramos

for ageb in agebs.getFeatures(): print "para el ageb",str(ageb['ageb_id']) suma_volumen = 0 for tramo in tramos.getFeatures(QgsFeatureRequest().setFilterExpression('"ageb_id" = ' + "'" + str(ageb['ageb_id']) + "'" )): print "el tramo",str(tramo['EDGE_NUM']),"tiene un volumen",str(tramo['volumen']) suma_volumen += tramo['volumen'] print suma_volumen

Se asocia la suma del volumen al campo vol_dren para cada AGEB

agebs.startEditing()
ageb['vol_dren'] = suma_volumen
agebs.updateFeature(ageb)
agebs.commitChanges()

[edithvilla]

Capacidad de desalojo de la Red de Drenaje Profundo

Responsable: Edith Villa Fecha: 14 de septiembre de 2018

Capas de insumo: drenaje_profundo.shp ageb_cdmx.shp

Ubicación de trabajo: \SIG\desarrollo\megadapt\procesamiento\drenaje_profundo\insumos

Se calculó la capacidad de desalojo de la red de drenaje profundo desde QGIS

Con la herramienta "Calculadora de campos" se hicieron los calculos para obtener las capacidades de desalojo para Q50 y Q100:

1. Se crea la columna radio

   "Diameter" / 2

2. Se crea la columna tirante100

   "Diameter"

3. Se crea la columna tirante50

   "Diameter" / 2

4. Se crea la columna ang_t_100

   2 * ACOS (1 - ("tirante100 / "radio"))

5. Se crea la columna ang_t_50

   2 * ACOS (1 - ("tirante50 / "radio"))

6. Se crea la columna a_hidr100

(( "Diameter" ^ 2)/8) * ("ang_t_100" - sin("ang_t_100" ))

7. Se crea la columna a_hidr50

   (( "Diameter" ^ 2)/8) * ("ang_t_50" - sin("ang_t_50" ))

8. Se crea la columna peri_h100

   "radio" * "ang_t_100"

9. Se crea la columna peri_h50

   "radio" * "ang_t_50"

10. Se crea la columna rh100

    "a_hidr100" / "peri_h_100"

11. Se crea la columna rh50

    "a_hidr50" / "peri_h_50"

12. Se crea la columna Q100

    "a_hidr100" (1 / 0.011) ("rh100" ^ (2 / 3)) * ( "slope" ^ (1 / 2))

13. Se crea la columna Q50

    "a_hidr50" (1 / 0.011) ("rh100" ^ (2 / 3)) * ( "slope" ^ (1 / 2))

14. Se guardaron los cálculos en la capa: drenaje_profundo_desalojo.shp en:

\SIG\desarrollo\megadapt\procesamiento\drenaje_profundo\insumos

15. Se guardó la capa solo con los campos: MUID, Descriptio, Q_100, Q_50

[aleestrada]

Capaciad de desalojo de la red de drenaje primario

Objetivo: Revisar los datos de la red de drenaje primario.

Como resultado de la reunión que se tuvo con Luis el día 12 de noviembre, se decidió considerar los datos de la pendiente del tubo como 2/1000 sin multiplicarlo por la longitud.

Se subieron los nuevos datos a la ruta: MEGADAPT\mega_procesamiento\banco_datos\capacidad_drenaje\bd c Con el nombre: capacidad_drenaje_v2

@edithvilla se encargará junto con @serranoycandela de realizar la nueva capa de capacidad de desalojo por AGEB para el drenaje primario, para el periodo de secas.

[edithvilla]

continua en el #74