OpenCV es una biblioteca libre de visión artificial originalmente desarrollada por Intel. Se ha utilizado en infinidad de aplicaciones que abarcan desde sistemas de seguridad con detección de movimiento, aplicaciones de control de procesos donde se requiere reconocimiento de objetos y, por qué no, análisis de imágenes en el campo SIG.
En un post anterior se usó ndimage, del módulo python scipy, para acceder a subarrays de bloques 3x3. Sin embargo, con el módulo cv2 de OpenCV y numpy el procedimiento es más directo.El código aquí delineado se basó en ésta referencia, con la ayuda de esta otra, y se probó con el ráster tif, 9x9, de la imagen siguiente:
El código a continuación emplea una ventana de 5x5 para muestrear todos los subarrays y, para soslayar los efectos de borde, se crea un anillo de valores -999 (padding) en la imagen original para que los subarrays de los píxeles de los bordes puedan ser extraidos con el arreglo 5x5 en su ventana de convolución al igual que los píxeles interiores.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | import numpy as np import cv2 img = cv2.imread('/home/zeito/pyqgis_data/dem_part_raster.tif', cv2.IMREAD_UNCHANGED) img = img.astype(np.float32) kernel = np.ones((5,5),np.float32) (iH, iW) = img.shape[:2] (kH, kW) = kernel.shape[:2] # "pad" the borders of the input image so the spatial # size (i.e., width and height) are not reduced pad = (kW - 1) / 2 print img img = cv2.copyMakeBorder(img, pad, pad, pad, pad, cv2.BORDER_CONSTANT, value =-999) # loop over the input image, "sliding" the kernel across # each (x, y)-coordinate from left-to-right and top to # bottom for y in np.arange(pad, iH + pad): for x in np.arange(pad, iW + pad): #extracting subarrays according kernel roi = img[y - pad:y + pad + 1, x - pad:x + pad + 1] print roi |
Cuando el código anterior se ejecutó en la Python Console de QGIS se puede corroborar que los 81 subarrays se imprimen con el arreglo esperado; tal como se observa a continuación:
[[ 2703. 2781. 2836. 2894. 2928. 2966. 3002. 3037. 3078.] [ 2682. 2767. 2819. 2863. 2887. 2916. 2950. 3004. 3033.] [ 2652. 2729. 2776. 2809. 2840. 2880. 2932. 2963. 2983.] [ 2671. 2703. 2738. 2775. 2820. 2858. 2878. 2919. 2952.] [ 2671. 2703. 2738. 2775. 2820. 2858. 2878. 2919. 2952.] [ 2644. 2670. 2704. 2747. 2784. 2812. 2857. 2901. 2935.] [ 2616. 2654. 2687. 2717. 2751. 2788. 2829. 2873. 2907.] [ 2633. 2672. 2696. 2725. 2751. 2780. 2811. 2855. 2896.] [ 2636. 2672. 2697. 2723. 2749. 2782. 2812. 2843. 2866.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ -999. -999. 2703. 2781. 2836.] [ -999. -999. 2682. 2767. 2819.] [ -999. -999. 2652. 2729. 2776.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ -999. 2703. 2781. 2836. 2894.] [ -999. 2682. 2767. 2819. 2863.] [ -999. 2652. 2729. 2776. 2809.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2703. 2781. 2836. 2894. 2928.] [ 2682. 2767. 2819. 2863. 2887.] [ 2652. 2729. 2776. 2809. 2840.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2781. 2836. 2894. 2928. 2966.] [ 2767. 2819. 2863. 2887. 2916.] [ 2729. 2776. 2809. 2840. 2880.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2836. 2894. 2928. 2966. 3002.] [ 2819. 2863. 2887. 2916. 2950.] [ 2776. 2809. 2840. 2880. 2932.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2894. 2928. 2966. 3002. 3037.] [ 2863. 2887. 2916. 2950. 3004.] [ 2809. 2840. 2880. 2932. 2963.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2928. 2966. 3002. 3037. 3078.] [ 2887. 2916. 2950. 3004. 3033.] [ 2840. 2880. 2932. 2963. 2983.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 2966. 3002. 3037. 3078. -999.] [ 2916. 2950. 3004. 3033. -999.] [ 2880. 2932. 2963. 2983. -999.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.] [ 3002. 3037. 3078. -999. -999.] [ 2950. 3004. 3033. -999. -999.] [ 2932. 2963. 2983. -999. -999.]] [[ -999. -999. -999. -999. -999.] [ -999. -999. 2703. 2781. 2836.] [ -999. -999. 2682. 2767. 2819.] [ -999. -999. 2652. 2729. 2776.] [ -999. -999. 2671. 2703. 2738.]] [[ -999. -999. -999. -999. -999.] [ -999. 2703. 2781. 2836. 2894.] [ -999. 2682. 2767. 2819. 2863.] [ -999. 2652. 2729. 2776. 2809.] [ -999. 2671. 2703. 2738. 2775.]] [[ -999. -999. -999. -999. -999.] [ 2703. 2781. 2836. 2894. 2928.] [ 2682. 2767. 2819. 2863. 2887.] [ 2652. 2729. 2776. 2809. 2840.] [ 2671. 2703. 2738. 2775. 2820.]] [[ -999. -999. -999. -999. -999.] [ 2781. 2836. 2894. 2928. 2966.] [ 2767. 2819. 2863. 2887. 2916.] [ 2729. 2776. 2809. 2840. 2880.] [ 2703. 2738. 2775. 2820. 2858.]] [[ -999. -999. -999. -999. -999.] [ 2836. 2894. 2928. 2966. 3002.] [ 2819. 2863. 2887. 2916. 2950.] [ 2776. 2809. 2840. 2880. 2932.] [ 2738. 2775. 2820. 2858. 2878.]] [[ -999. -999. -999. -999. -999.] [ 2894. 2928. 2966. 3002. 3037.] [ 2863. 2887. 2916. 2950. 3004.] [ 2809. 2840. 2880. 2932. 2963.] [ 2775. 2820. 2858. 2878. 2919.]] [[ -999. -999. -999. -999. -999.] [ 2928. 2966. 3002. 3037. 3078.] [ 2887. 2916. 2950. 3004. 3033.] [ 2840. 2880. 2932. 2963. 2983.] [ 2820. 2858. 2878. 2919. 2952.]] [[ -999. -999. -999. -999. -999.] [ 2966. 3002. 3037. 3078. -999.] [ 2916. 2950. 3004. 3033. -999.] [ 2880. 2932. 2963. 2983. -999.] [ 2858. 2878. 2919. 2952. -999.]] [[ -999. -999. -999. -999. -999.] [ 3002. 3037. 3078. -999. -999.] [ 2950. 3004. 3033. -999. -999.] [ 2932. 2963. 2983. -999. -999.] [ 2878. 2919. 2952. -999. -999.]] [[ -999. -999. 2703. 2781. 2836.] [ -999. -999. 2682. 2767. 2819.] [ -999. -999. 2652. 2729. 2776.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2671. 2703. 2738.]] [[ -999. 2703. 2781. 2836. 2894.] [ -999. 2682. 2767. 2819. 2863.] [ -999. 2652. 2729. 2776. 2809.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2671. 2703. 2738. 2775.]] [[ 2703. 2781. 2836. 2894. 2928.] [ 2682. 2767. 2819. 2863. 2887.] [ 2652. 2729. 2776. 2809. 2840.] [ 2671. 2703. 2738. 2775. 2820.] [ 2671. 2703. 2738. 2775. 2820.]] [[ 2781. 2836. 2894. 2928. 2966.] [ 2767. 2819. 2863. 2887. 2916.] [ 2729. 2776. 2809. 2840. 2880.] [ 2703. 2738. 2775. 2820. 2858.] [ 2703. 2738. 2775. 2820. 2858.]] [[ 2836. 2894. 2928. 2966. 3002.] [ 2819. 2863. 2887. 2916. 2950.] [ 2776. 2809. 2840. 2880. 2932.] [ 2738. 2775. 2820. 2858. 2878.] [ 2738. 2775. 2820. 2858. 2878.]] [[ 2894. 2928. 2966. 3002. 3037.] [ 2863. 2887. 2916. 2950. 3004.] [ 2809. 2840. 2880. 2932. 2963.] [ 2775. 2820. 2858. 2878. 2919.] [ 2775. 2820. 2858. 2878. 2919.]] [[ 2928. 2966. 3002. 3037. 3078.] [ 2887. 2916. 2950. 3004. 3033.] [ 2840. 2880. 2932. 2963. 2983.] [ 2820. 2858. 2878. 2919. 2952.] [ 2820. 2858. 2878. 2919. 2952.]] [[ 2966. 3002. 3037. 3078. -999.] [ 2916. 2950. 3004. 3033. -999.] [ 2880. 2932. 2963. 2983. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2858. 2878. 2919. 2952. -999.]] [[ 3002. 3037. 3078. -999. -999.] [ 2950. 3004. 3033. -999. -999.] [ 2932. 2963. 2983. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2878. 2919. 2952. -999. -999.]] [[ -999. -999. 2682. 2767. 2819.] [ -999. -999. 2652. 2729. 2776.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2644. 2670. 2704.]] [[ -999. 2682. 2767. 2819. 2863.] [ -999. 2652. 2729. 2776. 2809.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2644. 2670. 2704. 2747.]] [[ 2682. 2767. 2819. 2863. 2887.] [ 2652. 2729. 2776. 2809. 2840.] [ 2671. 2703. 2738. 2775. 2820.] [ 2671. 2703. 2738. 2775. 2820.] [ 2644. 2670. 2704. 2747. 2784.]] [[ 2767. 2819. 2863. 2887. 2916.] [ 2729. 2776. 2809. 2840. 2880.] [ 2703. 2738. 2775. 2820. 2858.] [ 2703. 2738. 2775. 2820. 2858.] [ 2670. 2704. 2747. 2784. 2812.]] [[ 2819. 2863. 2887. 2916. 2950.] [ 2776. 2809. 2840. 2880. 2932.] [ 2738. 2775. 2820. 2858. 2878.] [ 2738. 2775. 2820. 2858. 2878.] [ 2704. 2747. 2784. 2812. 2857.]] [[ 2863. 2887. 2916. 2950. 3004.] [ 2809. 2840. 2880. 2932. 2963.] [ 2775. 2820. 2858. 2878. 2919.] [ 2775. 2820. 2858. 2878. 2919.] [ 2747. 2784. 2812. 2857. 2901.]] [[ 2887. 2916. 2950. 3004. 3033.] [ 2840. 2880. 2932. 2963. 2983.] [ 2820. 2858. 2878. 2919. 2952.] [ 2820. 2858. 2878. 2919. 2952.] [ 2784. 2812. 2857. 2901. 2935.]] [[ 2916. 2950. 3004. 3033. -999.] [ 2880. 2932. 2963. 2983. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2812. 2857. 2901. 2935. -999.]] [[ 2950. 3004. 3033. -999. -999.] [ 2932. 2963. 2983. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2857. 2901. 2935. -999. -999.]] [[ -999. -999. 2652. 2729. 2776.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2644. 2670. 2704.] [ -999. -999. 2616. 2654. 2687.]] [[ -999. 2652. 2729. 2776. 2809.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2644. 2670. 2704. 2747.] [ -999. 2616. 2654. 2687. 2717.]] [[ 2652. 2729. 2776. 2809. 2840.] [ 2671. 2703. 2738. 2775. 2820.] [ 2671. 2703. 2738. 2775. 2820.] [ 2644. 2670. 2704. 2747. 2784.] [ 2616. 2654. 2687. 2717. 2751.]] [[ 2729. 2776. 2809. 2840. 2880.] [ 2703. 2738. 2775. 2820. 2858.] [ 2703. 2738. 2775. 2820. 2858.] [ 2670. 2704. 2747. 2784. 2812.] [ 2654. 2687. 2717. 2751. 2788.]] [[ 2776. 2809. 2840. 2880. 2932.] [ 2738. 2775. 2820. 2858. 2878.] [ 2738. 2775. 2820. 2858. 2878.] [ 2704. 2747. 2784. 2812. 2857.] [ 2687. 2717. 2751. 2788. 2829.]] [[ 2809. 2840. 2880. 2932. 2963.] [ 2775. 2820. 2858. 2878. 2919.] [ 2775. 2820. 2858. 2878. 2919.] [ 2747. 2784. 2812. 2857. 2901.] [ 2717. 2751. 2788. 2829. 2873.]] [[ 2840. 2880. 2932. 2963. 2983.] [ 2820. 2858. 2878. 2919. 2952.] [ 2820. 2858. 2878. 2919. 2952.] [ 2784. 2812. 2857. 2901. 2935.] [ 2751. 2788. 2829. 2873. 2907.]] [[ 2880. 2932. 2963. 2983. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2812. 2857. 2901. 2935. -999.] [ 2788. 2829. 2873. 2907. -999.]] [[ 2932. 2963. 2983. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2857. 2901. 2935. -999. -999.] [ 2829. 2873. 2907. -999. -999.]] [[ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2644. 2670. 2704.] [ -999. -999. 2616. 2654. 2687.] [ -999. -999. 2633. 2672. 2696.]] [[ -999. 2671. 2703. 2738. 2775.] [ -999. 2671. 2703. 2738. 2775.] [ -999. 2644. 2670. 2704. 2747.] [ -999. 2616. 2654. 2687. 2717.] [ -999. 2633. 2672. 2696. 2725.]] [[ 2671. 2703. 2738. 2775. 2820.] [ 2671. 2703. 2738. 2775. 2820.] [ 2644. 2670. 2704. 2747. 2784.] [ 2616. 2654. 2687. 2717. 2751.] [ 2633. 2672. 2696. 2725. 2751.]] [[ 2703. 2738. 2775. 2820. 2858.] [ 2703. 2738. 2775. 2820. 2858.] [ 2670. 2704. 2747. 2784. 2812.] [ 2654. 2687. 2717. 2751. 2788.] [ 2672. 2696. 2725. 2751. 2780.]] [[ 2738. 2775. 2820. 2858. 2878.] [ 2738. 2775. 2820. 2858. 2878.] [ 2704. 2747. 2784. 2812. 2857.] [ 2687. 2717. 2751. 2788. 2829.] [ 2696. 2725. 2751. 2780. 2811.]] [[ 2775. 2820. 2858. 2878. 2919.] [ 2775. 2820. 2858. 2878. 2919.] [ 2747. 2784. 2812. 2857. 2901.] [ 2717. 2751. 2788. 2829. 2873.] [ 2725. 2751. 2780. 2811. 2855.]] [[ 2820. 2858. 2878. 2919. 2952.] [ 2820. 2858. 2878. 2919. 2952.] [ 2784. 2812. 2857. 2901. 2935.] [ 2751. 2788. 2829. 2873. 2907.] [ 2751. 2780. 2811. 2855. 2896.]] [[ 2858. 2878. 2919. 2952. -999.] [ 2858. 2878. 2919. 2952. -999.] [ 2812. 2857. 2901. 2935. -999.] [ 2788. 2829. 2873. 2907. -999.] [ 2780. 2811. 2855. 2896. -999.]] [[ 2878. 2919. 2952. -999. -999.] [ 2878. 2919. 2952. -999. -999.] [ 2857. 2901. 2935. -999. -999.] [ 2829. 2873. 2907. -999. -999.] [ 2811. 2855. 2896. -999. -999.]] [[ -999. -999. 2671. 2703. 2738.] [ -999. -999. 2644. 2670. 2704.] [ -999. -999. 2616. 2654. 2687.] [ -999. -999. 2633. 2672. 2696.] [ -999. -999. 2636. 2672. 2697.]] [[ -999. 2671. 2703. 2738. 2775.] [ -999. 2644. 2670. 2704. 2747.] [ -999. 2616. 2654. 2687. 2717.] [ -999. 2633. 2672. 2696. 2725.] [ -999. 2636. 2672. 2697. 2723.]] [[ 2671. 2703. 2738. 2775. 2820.] [ 2644. 2670. 2704. 2747. 2784.] [ 2616. 2654. 2687. 2717. 2751.] [ 2633. 2672. 2696. 2725. 2751.] [ 2636. 2672. 2697. 2723. 2749.]] [[ 2703. 2738. 2775. 2820. 2858.] [ 2670. 2704. 2747. 2784. 2812.] [ 2654. 2687. 2717. 2751. 2788.] [ 2672. 2696. 2725. 2751. 2780.] [ 2672. 2697. 2723. 2749. 2782.]] [[ 2738. 2775. 2820. 2858. 2878.] [ 2704. 2747. 2784. 2812. 2857.] [ 2687. 2717. 2751. 2788. 2829.] [ 2696. 2725. 2751. 2780. 2811.] [ 2697. 2723. 2749. 2782. 2812.]] [[ 2775. 2820. 2858. 2878. 2919.] [ 2747. 2784. 2812. 2857. 2901.] [ 2717. 2751. 2788. 2829. 2873.] [ 2725. 2751. 2780. 2811. 2855.] [ 2723. 2749. 2782. 2812. 2843.]] [[ 2820. 2858. 2878. 2919. 2952.] [ 2784. 2812. 2857. 2901. 2935.] [ 2751. 2788. 2829. 2873. 2907.] [ 2751. 2780. 2811. 2855. 2896.] [ 2749. 2782. 2812. 2843. 2866.]] [[ 2858. 2878. 2919. 2952. -999.] [ 2812. 2857. 2901. 2935. -999.] [ 2788. 2829. 2873. 2907. -999.] [ 2780. 2811. 2855. 2896. -999.] [ 2782. 2812. 2843. 2866. -999.]] [[ 2878. 2919. 2952. -999. -999.] [ 2857. 2901. 2935. -999. -999.] [ 2829. 2873. 2907. -999. -999.] [ 2811. 2855. 2896. -999. -999.] [ 2812. 2843. 2866. -999. -999.]] [[ -999. -999. 2644. 2670. 2704.] [ -999. -999. 2616. 2654. 2687.] [ -999. -999. 2633. 2672. 2696.] [ -999. -999. 2636. 2672. 2697.] [ -999. -999. -999. -999. -999.]] [[ -999. 2644. 2670. 2704. 2747.] [ -999. 2616. 2654. 2687. 2717.] [ -999. 2633. 2672. 2696. 2725.] [ -999. 2636. 2672. 2697. 2723.] [ -999. -999. -999. -999. -999.]] [[ 2644. 2670. 2704. 2747. 2784.] [ 2616. 2654. 2687. 2717. 2751.] [ 2633. 2672. 2696. 2725. 2751.] [ 2636. 2672. 2697. 2723. 2749.] [ -999. -999. -999. -999. -999.]] [[ 2670. 2704. 2747. 2784. 2812.] [ 2654. 2687. 2717. 2751. 2788.] [ 2672. 2696. 2725. 2751. 2780.] [ 2672. 2697. 2723. 2749. 2782.] [ -999. -999. -999. -999. -999.]] [[ 2704. 2747. 2784. 2812. 2857.] [ 2687. 2717. 2751. 2788. 2829.] [ 2696. 2725. 2751. 2780. 2811.] [ 2697. 2723. 2749. 2782. 2812.] [ -999. -999. -999. -999. -999.]] [[ 2747. 2784. 2812. 2857. 2901.] [ 2717. 2751. 2788. 2829. 2873.] [ 2725. 2751. 2780. 2811. 2855.] [ 2723. 2749. 2782. 2812. 2843.] [ -999. -999. -999. -999. -999.]] [[ 2784. 2812. 2857. 2901. 2935.] [ 2751. 2788. 2829. 2873. 2907.] [ 2751. 2780. 2811. 2855. 2896.] [ 2749. 2782. 2812. 2843. 2866.] [ -999. -999. -999. -999. -999.]] [[ 2812. 2857. 2901. 2935. -999.] [ 2788. 2829. 2873. 2907. -999.] [ 2780. 2811. 2855. 2896. -999.] [ 2782. 2812. 2843. 2866. -999.] [ -999. -999. -999. -999. -999.]] [[ 2857. 2901. 2935. -999. -999.] [ 2829. 2873. 2907. -999. -999.] [ 2811. 2855. 2896. -999. -999.] [ 2812. 2843. 2866. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ -999. -999. 2616. 2654. 2687.] [ -999. -999. 2633. 2672. 2696.] [ -999. -999. 2636. 2672. 2697.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ -999. 2616. 2654. 2687. 2717.] [ -999. 2633. 2672. 2696. 2725.] [ -999. 2636. 2672. 2697. 2723.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2616. 2654. 2687. 2717. 2751.] [ 2633. 2672. 2696. 2725. 2751.] [ 2636. 2672. 2697. 2723. 2749.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2654. 2687. 2717. 2751. 2788.] [ 2672. 2696. 2725. 2751. 2780.] [ 2672. 2697. 2723. 2749. 2782.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2687. 2717. 2751. 2788. 2829.] [ 2696. 2725. 2751. 2780. 2811.] [ 2697. 2723. 2749. 2782. 2812.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2717. 2751. 2788. 2829. 2873.] [ 2725. 2751. 2780. 2811. 2855.] [ 2723. 2749. 2782. 2812. 2843.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2751. 2788. 2829. 2873. 2907.] [ 2751. 2780. 2811. 2855. 2896.] [ 2749. 2782. 2812. 2843. 2866.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2788. 2829. 2873. 2907. -999.] [ 2780. 2811. 2855. 2896. -999.] [ 2782. 2812. 2843. 2866. -999.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]] [[ 2829. 2873. 2907. -999. -999.] [ 2811. 2855. 2896. -999. -999.] [ 2812. 2843. 2866. -999. -999.] [ -999. -999. -999. -999. -999.] [ -999. -999. -999. -999. -999.]]
No hay comentarios:
Publicar un comentario