Research article Special Issues

Optimal wind farm sitting using high-resolution digital elevation models and randomized optimization

  • We investigate the problem of wind farm design in isolated mountainous areas. We first describe a remote sensing approach for the terrain reconstruction of complex terrains. We then employ a well--known evolutionary optimization algorithm to find the optimal wind farm layout. Although the algorithm has been efficiently used for off--shore or smooth on--shore areas, we show that its performance is significantly affected by the complex topography. Moreover, we illustrate how a priori information can be exploited to improve both the computational time and efficiency of the optimization algorithm.

    Citation: John Koutroumpas, Konstantinos Koutroumpas. Optimal wind farm sitting using high-resolution digital elevation models and randomized optimization[J]. AIMS Energy, 2015, 3(4): 505-524. doi: 10.3934/energy.2015.4.505

    Related Papers:

    [1] Jianzhong Shi, Ying Song . Mathematical analysis of a simplified general type-2 fuzzy PID controller. Mathematical Biosciences and Engineering, 2020, 17(6): 7994-8036. doi: 10.3934/mbe.2020406
    [2] Baohua Hu, Yong Wang, Jingsong Mu . A new fractional fuzzy dispersion entropy and its application in muscle fatigue detection. Mathematical Biosciences and Engineering, 2024, 21(1): 144-169. doi: 10.3934/mbe.2024007
    [3] Li Yang, Kai Zou, Kai Gao, Zhiyi Jiang . A fuzzy DRBFNN-based information security risk assessment method in improving the efficiency of urban development. Mathematical Biosciences and Engineering, 2022, 19(12): 14232-14250. doi: 10.3934/mbe.2022662
    [4] Qiwen Wang, Guibao Song, Xiuxia Yang . Mixed-attitude three-way decision model for aerial targets: Threat assessment based on IF-VIKOR-GRA method. Mathematical Biosciences and Engineering, 2023, 20(12): 21514-21536. doi: 10.3934/mbe.2023952
    [5] Yuting Zhu, Wenyu Zhang, Junjie Hou, Hainan Wang, Tingting Wang, Haining Wang . The large-scale group consensus multi-attribute decision-making method based on probabilistic dual hesitant fuzzy sets. Mathematical Biosciences and Engineering, 2024, 21(3): 3944-3966. doi: 10.3934/mbe.2024175
    [6] Han Wu, Junwu Wang, Sen Liu, Tingyou Yang . Research on decision-making of emergency plan for waterlogging disaster in subway station project based on linguistic intuitionistic fuzzy set and TOPSIS. Mathematical Biosciences and Engineering, 2020, 17(5): 4825-4851. doi: 10.3934/mbe.2020263
    [7] Muhammad Akram, Ahmad N. Al-Kenani, Anam Luqman . Degree based models of granular computing under fuzzy indiscernibility relations. Mathematical Biosciences and Engineering, 2021, 18(6): 8415-8443. doi: 10.3934/mbe.2021417
    [8] Bo Sun, Ming Wei, Wei Wu, Binbin Jing . A novel group decision making method for airport operational risk management. Mathematical Biosciences and Engineering, 2020, 17(3): 2402-2417. doi: 10.3934/mbe.2020130
    [9] Yevgeniy Bodyanskiy, Olha Chala, Natalia Kasatkina, Iryna Pliss . Modified generalized neo-fuzzy system with combined online fast learning in medical diagnostic task for situations of information deficit. Mathematical Biosciences and Engineering, 2022, 19(8): 8003-8018. doi: 10.3934/mbe.2022374
    [10] Muhammad Akram, Saba Siddique, Majed G. Alharbi . Clustering algorithm with strength of connectedness for $ m $-polar fuzzy network models. Mathematical Biosciences and Engineering, 2022, 19(1): 420-455. doi: 10.3934/mbe.2022021
  • We investigate the problem of wind farm design in isolated mountainous areas. We first describe a remote sensing approach for the terrain reconstruction of complex terrains. We then employ a well--known evolutionary optimization algorithm to find the optimal wind farm layout. Although the algorithm has been efficiently used for off--shore or smooth on--shore areas, we show that its performance is significantly affected by the complex topography. Moreover, we illustrate how a priori information can be exploited to improve both the computational time and efficiency of the optimization algorithm.


    As an effective and low-cost method to achieve high precision deformation information of large areas, Interferometry Synthetic Aperture Radar (InSAR) technology is widely used in topographic survey, earthquake monitoring, ground surface subsidence monitoring and many other applications [1,2,3]. However, conventional InSAR is affected by the temporal decorrelation and atmospheric effects, which seriously restricts the accuracy and reliability in deformation monitoring [4]. Based on plenty of relatively stable, high coherent points from time-series SAR images and implementing differential Interferometry Synthetic Aperture Radar (DInSAR) on them, Permanent Scatterer InSAR (PSInSAR) can mitigate the spatial and temporal decorrelation and atmospheric delay effect to get relatively high accuracy and reliable information of the earth surface deformation [5,6,7,8,9,10,11,12].

    One of the key steps for PSInSAR technology is to analyze and process those stable scatterers with the high coherent characteristic. The phase stability of PS determines that the deformation results of these scatterers is reliable and can reflect the ground deformation if the density of PS is sufficient [13,14]. So how to identify and select PS automatically and effectively from SAR images is one of the most critical procedures in PSInSAR technology.

    Ferretti et al. firstly proposed the algorithms of PS identification and selection, including amplitude dispersion thresholding, temporal correlation thresholding and phase deviation thresholding [15,16], and some related experiments were implemented correspondingly. Among them, the amplitude dispersion thresholding needs a large number of images for the statistical analysis. However, in the amplitude dispersion thresholding method, only the varying of amplitude is used without considering the corresponding amplitude value. In addition, its reliability is related to the quality of the image preprocessing and the setting of the threshold [17,18]. The results based on temporal correlation thresholding are actually some stable regions rather than individual pixels. Their reliability is related to the window size and the selection of the threshold directly. Phase dispersion thresholding only considers phase information, so the results are easily affected by decorrelation, system thermal noise and so on [19]. Kampes et al. proposed a new PS selection algorithm by setting an amplitude threshold to identify high coherent points, which can detect high amplitude points easily [20,21]. Considering the phase stability, Hooper presented a new PS identification method based on the spatial correlation analysis of interferometric phase and implemented some related deformation extraction experiments [22,23].

    To automatically identify PS, Chen et al. put forward a Double-Thresholding algorithm (DTA) by combining both correlation coefficient threshold and amplitude dispersion index threshold [18]. Luo et al. proposed a method named amplitude information dual-threshold approach (AIDTA). With this method, those points with high amplitude values are firstly selected as permanent scatterer candidates (PSCs), then the eventual PS are taken out of PSCs by setting a suitable amplitude dispersion index (ADI) threshold [24]. Azadnejad et al. presented a new method which is based on minimizing Amplitude Dispersion Index (ADI) by finding the optimum scattering mechanism to increase the number of PSC and PS pixels [25]. Similarly, taking into account the characteristics of strong reflectivity and stable backscattering of PS to radar waves, Hu et al. presented a new method by seriously setting tri-threshold with time-series correlation coefficient, amplitude dispersion index and amplitude index [19]. Based on the characteristics and flaws of all kinds of PS selection methods, such as amplitude stability, phase stability, temporal sequence correlation and amplitude dispersion index threshold, Long et al. proposed a new approach considering the union and intersection of temporal sequence correlation coefficient threshold and amplitude dispersion index threshold [26]. Tao et al. presented a new method for identifying and selecting PS by introducing wavelet phase analysis and taking into account the stability of amplitude and phase [27].

    All the above mentioned methods enriched the PS identification technology and achieved superior effects and they are suitable for those land cover with clear characters (e.g. unstable plants and stable man-made buildings) areas. However, these methods identify PS by setting one or several thresholds under the “black or white” rule without taking the fuzzy characteristics of the PS into account and can easily fall into a logical paradox especially in the area covered by objects possessing ambiguous boundary for interested features. For example, in an mountain area covered with various reflective land covers (there is no clear boundary of the reflection intensity), considering both the amplitude and amplitude dispersion index thresholds, it might be obviously unreasonable that if a pixel, with amplitude equal or greater than 64 and amplitude dispersion equal or less than 0.4, is selected as PS, while another pixel with the amplitude 63.99999 and amplitude dispersion 0.2 is excluded from PS. In fact, this can exclude many valuable PS, which usually leads to unreliable subsidence analysis results because of sparse density of PS. Based on this, this paper tries to find a suitable PS selection method to deal with those areas with fuzzy features. The fuzzy set (FS) theory was proposed in 1965 [28], since then it is widely used in various fields [29,30,31]. It adopts quantitative method to study and deal with the fuzzy phenomenon and always makes the solution in conformity with the practical situation. In other words, FS method can achieve more accurate result by a fuzzy point of view. In order to obtain more reasonable PS in those areas covered by objects possessing ambiguous boundary for interested features, a new fuzzy PS (FPS) concept is first proposed based on fuzzy set theory in this paper. Taking the two indexes amplitude and amplitude dispersion in AIDTA as an example, the selection method is presented. Two groups of real data are performed to validate the proposed method. The results indicate that, compared with traditional method, FPS can not only eliminate the negative effect of the “black or white” rule because of setting single threshold on each PS characteristic, but also gain a unified quantitative index to guide the selection and quality classification of PS by fusing different PS characters comprehensively.

    The paper is organized as follows. Section 2 introduces the detail description of amplitude information dual-threshold approach (AIDTA). The fuzzy PS identification method is presented in Section 3. Test results and data analysis are given in Section 4. The last part of the paper is the conclusions attained through this study.

    Based on the characteristics of PS with strong and stable Radar reflectivity, a method named amplitude information dual-threshold approach (AIDTA) was proposed [24]. Firstly, the pixels with high amplitude values are selected as PS candidates (PSC) using amplitude thresholding in view of their strong reflectivity. Then, PS are selected from PSC according to an amplitude dispersion index thresholding for their stable reflectivity.

    1) Consider $N$ time-series images covering the same area. The amplitude threshold is computed as:

    ${T_A} = \min \left\{ {\frac{1}{{m \times n}}\sum\limits_{i = 1}^m {\sum\limits_{j = 1}^n {{A_{{k_{i, j}}}}\;\left| {k = 1, 2, \cdots , N} \right.} } } \right\}$ (1)

    where $m$ and $n$ represent the number of rows and columns of the SAR images; ${A_k}$ is the calibrated amplitude array of the ${{\rm{k}}^{th}}$ SAR image [32]. From Eq (1), the minimum of the amplitude mean of each image is set as the amplitude threshold.

    For a fixed ${T_A}$, PSC are those pixels satisfying the following formula:

    $\min \left\{ {{A_{{k_{i, j}}}}\left| {k = 1, \ldots , N} \right.} \right\} \geqslant {T_A}$ (2)

    The above amplitude threshold considers only the strong reflectivity without regarding to the scattering stability of PS, making the result not reliable enough. In the following step, PS are further selected from the PSC using amplitude dispersion index threshold for their stable reflectivity.

    2) Amplitude dispersion index ${D_A}$ and phase dispersion index ${D_\varphi }$ are used to quantitatively express the stability of the PS scattering. The distribution of the amplitude values $f(A)$ is given by the Rice distribution [16]:

    $f(A) = \frac{A}{{\sigma _n^2}}{I_0}\left[ {\frac{{Ag}}{{\sigma _n^2}}} \right]{e^{ - {{({A^2} + {g^2})} \mathord{\left/ {\vphantom {{({A^2} + {g^2})} {(2\sigma _n^2)}}} \right. } {(2\sigma _n^2)}}}}$ (3)

    where ${I_0}$ is the modified Bessel function; $\sigma _n^2$ is the power of complex circular Gaussian distributed noise; $g$ is complex reflectivity and $A > 0$.

    The shape of the Rice Distribution depends on the SNR (i.e., the ratio ${g \mathord{\left/ {\vphantom {g {{\sigma _n}}}} \right. } {{\sigma _n}}}$). For low SNR, the probability density function (PDF) tends to a Rayleigh distribution. For high SNR (${g \mathord{\left/ {\vphantom {g {{\sigma _n}}}} \right. } {{\sigma _n}}} > 4$), $f(A)$ approaches a Gauss distribution, and the following equations hold:

    $\left\{ σAσnR=σnIDφσnIgσAmA=DA
    \right.$
    (4)

    where ${\sigma _{nR}}$ and ${\sigma _{nI}}$ are the noise standard deviation of the real and imaginary part; ${m_A}$ and ${\sigma _A}$ are the mean and the standard deviation of the amplitude values ${A_k}$; ${D_A}$ is called amplitude dispersion index (ADI). According to the above equation, the amplitude dispersion index is a measure of phase stability, at lease for high SNR pixels [16].

    Based on the above analysis, the selected PSC can meet the high SNR, and the points are selected as PS if the amplitude dispersion is below a threshold in the second step which means that the final PS are with stable reflectivity. Therefore, PS can be finalized through further analysis of the amplitude dispersion index on PSC.

    The basic idea of fuzzy set is to make the absolute membership degree in common sets flexible, specifically, by expanding the membership degree of element from the only 0 or 1 to arbitrary value in [0, 1]. Some relevant concepts are briefly described as follows [33]:

    Let $U$ be the universal set. A fuzzy set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} $ of $U$ can be mathematically described as:

    $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} = \left\{ {\left( {u, {\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} }}\left( u \right)} \right)\left| {u \in U} \right.} \right\}$ (5)

    where ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} }}\left( u \right):U \to \left[ {0, 1} \right]$ represents the membership degree of an element belonging to $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} $.

    $\lambda {\rm{ - }}$ cut of a fuzzy set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} $ in $U$ can be defined as

    ${(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} )_\lambda } = \left\{ {u{\rm{ | }}{\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} }}\left( u \right) \geqslant \lambda , u \in U} \right\}$ (6)

    where $\lambda $ is called a threshold. Obviously, ${(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} )_\lambda }$ is a common set and ${(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} )_{{\lambda _1}}} \supset {(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{F} )_{{\lambda _2}}}$ if ${\lambda _1} < {\lambda _2}$.

    T-norms, T-conorms and negated function are the basic operations on fuzzy sets to compute the membership degree of their union, intersection and complement, respectively. Minimum, Algebra product and bounded product are common T-norms operators.

    PS are with the characteristics of steadier scattering and stronger radar echo signal. While “steadier” and “stronger” are fuzzy descriptions of PS properties. Based on previous research PS and non-PS are classified by setting threshold respectively for each characteristic. One pixel will be excluded from PS once it does not meet any threshold rule. This in turn will lead to an “either this or that” logic paradox especially in the area covered by objects possessing ambiguous boundary for interested features. Because it is very difficult to set a precise and proper threshold in practice. For example, if pixels whose amplitude greater than or equal to 367.54 and amplitude dispersion less than or equal to 0.35 are identified as PS, it is more reasonable to classify the pixel with amplitude 367.30 and amplitude dispersion 0.27 to PS, while the previous method may make a rejection choice. Based on this, the concept of fuzzy PS is proposed and a fuzzy identification method is designed accordingly.

    On the base of fuzzy set theory, each point is no longer classified into either “PS” set or “no-PS” set absolutely, but has a membership degree to describe its level belonging to “PS”.

    Let $P$ represent all pixels in a study area. Each pixel $p \in P$ has $m$ feature indexes ${{\rm{C}}_1}, \ldots , {C_m}$ which are used to describe the characteristics of “PS”. Accordingly, there will be a fuzzy vector $V(p) = \left( {{u_1}(p), \ldots , {u_m}(p)} \right)$, where ${u_i}(i = 1, \ldots , m)$ is the value of $p$ in ${C_i}$. Define $m$ fuzzy sets ${{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A}}_{1}}, \ldots , {\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} _m}$, where ${\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} _i} = \left\{ {\left( {u, {\mu _{{{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }_i}}}\left( u \right)} \right)\left| {u \in {U_i}} \right.} \right\}$ is a fuzzy set describing feature ${C_i}$, ${U_i}$ is a universal set containing all values of $p$ in ${C_i}$. Since all the fuzzy sets are to describe the feature of PS, fuzzy PS set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} $ can be given by $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} = \bigcap\limits_{i = 1}^m {{{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }_i}} $ and the membership degree can be calculated by the following formula

    ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}\left( p \right) = T_{i = 1}^m{\mu _{{{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }_i}}}({u_i}(p))$ (7)

    where $T$ represents T-norms and Algebra product operation $ \bullet $ is selected in this paper. Thus, the fuzzy PS set can be written as $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} = \left\{ {(p, {\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}\left( p \right))\left| {p \in P} \right.} \right\}$. Accordingly, the $\lambda {\rm{ - }}$ cut of $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} $ can be symbolized as ${\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} _\lambda } = \left\{ {p{\rm{ | }}{\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}\left( p \right) \geqslant \lambda , p \in P} \right\}$, which means all the pixels whose membership degrees are greater than or equal to $\lambda $.

    According to the above theory, each pixel has a corresponding membership degree, which can be considered as the level it belongs to “PS”. A larger membership degree means the higher possibility that a certain pixel will be regarded as PS. It must be pointed out that, how to select a proper membership function in practice will directly affect the identification results. Since there is not a universal law to select a membership function, the reasonable and proper formulas for membership functions are more or less formulated with subjectivity or by experience.

    The new PS identification algorithm based on fuzzy set (fuzzy PS identification for short) are as follows:

    1) Select the study area. A matrix is used to describe the whole pixels which can be expressed by position.

    2) Choose $m$ feature indexes related to properties of PS, written as ${{\rm{C}}_1}, \ldots {C_m}$, such as “amplitude”, “amplitude dispersion” or “time-series correlation coefficient”.

    3) Define fuzzy set ${\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} _i}$ for ${C_i}$, such as “high” amplitude, “low” amplitude dispersion and so on. They are usually constructed based on some empirical models.

    4) Determine fuzzy PS set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} = \bigcap\limits_{i = 1}^m {{{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }_i}} $, ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}\left( p \right) = T_{i = 1}^m{\mu _{{{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }_i}}}({u_i}(p))$. The operator $T$ is decided on practical problem.

    5) Set membership degree threshold $\lambda $ to select the points satisfying ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}\left( p \right) \geqslant \lambda $. This is usually related to the threshold of traditional method.

    6) Determine practical PS considering the traditional method and the fuzzy method.

    The proposed method is a new characteristic fusion tool to identify PS. By fusing two or more characteristics, the practical PS can be selected on one hand, and the quality value of each PS can be obtained on the other hand. The former is conducted by controlling the membership degree threshold instead of setting single threshold for each characteristic. From Eq (7), it can be seen that the membership degree in fuzzy PS set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} $ is a fusion result taking the entire feature into account. Therefore, it is more reasonable than the traditional method by setting single threshold one by one.

    Twenty-seven TerraSAR-X images acquired from April 2012 to January 2013 are selected as the first test data set. The cropped images are 500 × 500 pixels, which cover part of 18,207 working face in Tunlan coal mine and its geographical location is about from east longitude 111 to 112 degree, and north latitude 37 to 38 degree.

    Taking the amplitude information dual-threshold approach (AIDTA) as an example, the traditional method (AIDTA) and fuzzy amplitude information dual-threshold approach (FAIDTA) are compared to show their relationship and difference. For simplicity, PS selected by AIDTA are shortly written as PS-AIDTA, PS selected by FAIDTA are simply described as PS-FAIDTA, and PS selected by FAIDTA but excluded from AIDTA is expressed as PS-FAIDTA/AIDTA.

    The calculated amplitude threshold for this test data was ${T_A} = 104.11$ based on Eq (1). Meanwhile, the amplitude dispersion threshold was set as ${T_B} = 0.32$ to make sure the selected PS candidates with high amplitude stability. Then, the calculated number of PS-AIDTA is 709. The FAIDTA method can be described as follows:

    Define that $ U = {({p}_{i, j})}_{I\times J},(i = 1, \dots , I, j = 1, \dots , J)$ represents all the pixels in study area, ${p_{i, j}}$ is the position of the $i$ row, $j$ column pixel, $I = 500$ and $J = 500$ are the number of rows and columns. The minimum amplitude value of time-series (amplitude value for short) ${c_{i, j}}$ and amplitude dispersion value ${q_{i, j}}$ for pixel ${p_{i, j}}$ were selected to describe the features of PS. Correspondingly, two fuzzy sets were defined.

    The first fuzzy set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} $ means that the amplitude value $c$ is high. The second fuzzy set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{B} $ means that the amplitude dispersion value $q$ is low. As a kind of methods for constructing fuzzy membership function, the assignment method is used in this paper [34,35]. The membership functions of $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} $ and $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{B} $ are shown as follows, which are Cauchy distributions

    ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }}(c) = \left\{ 00<cTA/211+(cTA/2TA/10)2.5TA/2<cmax{cij}
    \right.$
    (8)
    ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{B} }}(q) = \left\{ 02TBqmax{qi,j}11+(2TBq2TB/5)40<q<2TB
    \right.$
    (9)

    Figure 1 shows the membership functions of amplitude value and the amplitude dispersion value. It can be found that the former is monotone increasing and the latter is monotone decreasing. The membership degrees of the amplitude threshold (${T_A} = 104.11$) and the amplitude dispersion threshold (${T_B} = 0.32$) are greater than 0.9, and the membership function of the latter is more sensitive than the former, which conforms to the real expectation.

    Figure 1.  Membership Degree Functions of amplitude and amplitude dispersion. Note that $c$ is the amplitude value (unit: dB) and $q$ is the amplitude dispersion value (Dimensionless). The membership degree $\mu $ is a dimensionless quantity with a range [0, 1].

    According to fuzzy set theory, the membership degree of any pixel belonging to “PS” is

    ${\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} }}(p) = {\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{A} }}(c(p)) \bullet {\mu _{\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{B} }}(q(p))$ (10)

    The calculated minimum membership degree of PS-AIDTA is 0.92384.

    Figure 2 shows three overlay sub-graphs. The green ‘*’ means the PS-AIDTA, and the red ‘o’ represents the PS-FAIDTA. PS-AIDTA are all the same in the three sub-graphs because the thresholds of amplitude and amplitude dispersion are fixed, while PS-FAIDTA are different under three membership degree thresholds, and both of them are on top of the mean intensity image in study area. From all the three sub-graphs, it can be found that the PS-AIDTA and the PS-FAIDTA are mostly distributed in the strong reflected intensity area. With the decrease of membership degree threshold, the number of PS-FAIDTA grows. When the membership degree threshold was set as T = 0.92384 (the minimum membership degree of PS-AIDTA), the PS-FAIDTA include all the PS-AIDTA and the former is much more than the latter. This is consistent with the fuzzy PS theory—the $\lambda $ -cut of a fuzzy set is monotone decreasing with $\lambda $.

    Figure 2.  PS distribution under different membership degree thresholds in the mean-intensity map for Case 1. (a) with membership degree threshold T = 0.94; (b) with membership degree threshold T = 0.93; (c) with membership degree threshold T = 0.92384.

    Table 1 shows the statistic information under different membership degree thresholds. There are some abbreviations: TMD is the threshold of membership degree; n-FAIDTA, n-FAIDTA/AIDTA, and n-intersection stand for the number of PS-FAIDTA, PS-FAIDTA/AIDTA, and the intersection of FAIDTA and AIDTA respectively. It can be found that with the decrease of membership degree threshold, the number of PS-FAIDTA, PS-FAIDTA/AIDTA and their intersection increase. When the membership degree threshold was set as T = 0.92384, the PS-FAIDTA include all the PS-AIDTA and the former is about 37 percent more than the later. Figure 2 and Table 1 indicate that the PS-FAIDTA method is feasible and it can select some different pixels based on FAIDTA method. Following, the quantitative characteristics of these PS-FAIDTA/AIDTA are showed and analyzed.

    Table 1.  The number of PS-FAIDTA, PS-FAIDTA/AIDTA and intersection under different membership degree thresholds.
    TMD n-FAIDTA n-FAIDTA/AIDTA n-intersection
    0.92384 974 265 709
    0.925 941 234 707
    0.928 891 190 701
    0.93 841 146 695
    0.94 628 13 615
    0.95 440 1 489
    0.96 277 0 277

     | Show Table
    DownLoad: CSV

    Figure 3 shows the characters of the PS-FAIDTA/AIDTA under different membership degree thresholds. Sub-graph (a) offers the amplitude values of these scatterers, (b) shows their amplitude dispersion values and ${\rm{n}}$ is the number of selected pixels. It is clearly shown that all of them are not far from the “PS” characters. And as the membership degree threshold increases, the number of PS-FAIDTA/AIDTA decreases, while the amplitude values of these scatterers are closer to the threshold of AIDTA and the amplitude dispersion values get smaller and smaller, which conforms to the PS characteristics better. This can be confirmed from the two subplots on the right-most side of Figure 3 especially. According to the traditional method (AIDTA), these scatters should be belonged to “no-PS” because they do not simultaneously meet the requirements that the minimum amplitude value of the pixel sequence should be greater than a fixed threshold 104.11 and the amplitude dispersion should be smaller than 0.32. However, it is not a logical judgment to some extent. For our study, the first requirement is very strong. In practice, there are some pixels whose amplitude values are slightly smaller than the given threshold while their amplitude dispersion values are much lesser than the given threshold. Obviously, they conform to “PS” better. Similarly, there is the same problem for amplitude dispersion index threshold in traditional method. Actually, these problems are embodied in several traditional methods beside AIDTA because the fuzzy properties of PS are not considered. This demonstrates the significance of the FAIDTA method. In addition, it should be worth noting that some PS-FAIDTA/AIDTA may not be selected as practical PS because their characteristics are far from PS somewhat. However, with the increase of the membership degree threshold, the reliability of the PS-FAIDTA/AIDTA is improving. Therefore, we can select those with high membership degree as practical PS.

    Figure 3.  Characters of the PS-FAIDTA/AIDTA under different membership degree thresholds, (a) is about amplitude. (b) is about the amplitude dispersion results.

    Figure 4 shows the distribution of membership degrees of those PS-FAIDTA/AIDTA under different membership degree thresholds. We find that the membership degrees are relatively greater compared to the minimum membership degree of PS-AIDTA 0.92384 on the whole and the lesser the membership degree threshold is, the more PS-FAIDTA/AIDTA can be selected as. In practice, it is improper for them to be all excluded from PS. Take the sixth point in the third sub-graph for example, its amplitude is 100.31 and amplitude dispersion is 0.1888. Relatively, the former is slightly less than the threshold 104.11, while the latter is much less than 0.32. By our subjective judgment, it is more consistent with the characteristics of “PS”. Its membership degree, greater than 0.955 from the Figure 4, confirms the judgment, which reflects the rationality of the new method.

    Figure 4.  The membership degrees of the PS-FAIDTA/AIDTA under different membership degree thresholds.

    Table 2 lists eight typical pixels whose membership degrees are greater than 0.94, amplitude are greater than 100 and amplitude dispersion are less than 0.32. Clearly, the amplitude of all the pixels are slightly smaller than the threshold 104.11 relatively. That is why these points are excluded by AIDTA. In fact, it is reasonable for them to be selected as PS based on their high amplitude stability seen from amplitude dispersion and negligible undesirable amplitude information. Taking the eighth pixel of Table 2 for example, its amplitude 102.42 is only a litter below the threshold 104.11 a little relatively, while its amplitude dispersion 0.2711 is much smaller than 0.32. Comprehensively considering, it’s a better choice to classify it as a PS.

    Table 2.  The amplitude (A) and amplitude dispersion (AD) values of some pixels in PS-FAIDTA/AIDTA.
    Index Values
    AD 0.2463 0.2629 0.2694 0.2708 0.1888 0.2390 0.2710 0.2710
    A 101.38 101.96 101.52 101.75 100.31 101.52 101.29 102.42

     | Show Table
    DownLoad: CSV

    Another merit of FAIDTA is that quality results of all the interested pixels can be achieved based on their membership degrees in study area. Figure 5 shows two typical quality maps about PS-FAIDTA under membership degree threshold T = 0.8 and T = 0.9 respectively (mathematically named 0.8-T cut and 0.9-T cut of fuzzy set $\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle\thicksim}$}}{P} $). In Figure 5, the brighter the pixel is, the greater the membership degree it has. Thus it can be seen that the FAIDTA method can not only recognize the effective PS, but also obtain the quality result of these pixels. In other words, FAIDTA can also obtain the quantitative description for the selected PS, which may be used for evaluating the reliability of them.

    Figure 5.  Quality maps of PS-FAIDTA under membership degree threshold T = 0.8 and 0.9. (a) with membership degree threshold T = 0.8; (b) with membership degree threshold T = 0.9.

    In addition, the proposed method is also compared with the traditional algorithm (AIDTA) from the numerical computational costs. To compare the efficiency of different algorithms objectively, the same hardware and software were adopted. The hardware was a PC with a Core i7-8565U CPU and 16 Gb RAM. Matlab software (version R2017b) running on the Win10 system was used to design the different algorithms. For the test data in this case, the new method costs about 0.5 second and AIDTA costs about 0.05 seconds. Because the membership degrees of the candidate points need more computational costs. It should be pointed that the cost of the new method is closely related to the parameter setting in selecting the candidate points.

    Fifty-six ERS images acquired from April 1992 to January 2003 over Luxembourg are selected as the second test data set. The cropped images are 150 × 650 pixels and cover about 39 km2. The same methods and processing strategy are used. The amplitude dispersion threshold is set as 0.35. The calculated amplitude threshold and minimum membership degree of PS-AIDTA PS are 367 and 0.9263 respectively in this test. For keeping the paper concise, parts of the results are shown below.

    Figure 6 shows three corresponding overlay sub-graphs. The green ‘*’ and red ‘o’ mean the same as Figure 2. Similar regularity it demonstrates as Figure 1.

    Figure 6.  PS distribution under different membership degree thresholds in the mean-intensity map for Case 2. (a) with membership degree threshold T = 0.94; (b) with membership degree threshold T = 0.93; (c) with membership degree threshold T = 0.9263.

    Table 3 shows the statistic information about the PS-AIDTA and PS-FAIDTA under different membership degree thresholds. The abbreviations are the same as Table 1. When the membership degree threshold was set as T = 0.9263, the PS-FAIDTA include all the PS-AIDTA. Other analyses are similar to Figure 6.

    Table 3.  The number of PS-FAIDTA, PS-FAIDTA/AIDTA and intersection under different membership degree thresholds.
    TMD n-FAIDTA n-FAIDTA/AIDTA n-intersection
    0.9263 890 141 749
    0.93 844 99 745
    0.94 729 26 703

     | Show Table
    DownLoad: CSV

    Table 4 lists nine typical pixels with membership degrees greater than 0.94, amplitude greater than 360 and amplitude dispersion less than 0.35. Clearly, the amplitude values of all the pixels are slightly smaller than the threshold 367.54 calculated by AIDTA method relatively. That is why these points are excluded by AIDTA. In fact, it is reasonable for them to be selected as PS based on their high amplitude stability seen from amplitude dispersion and negligible undesirable amplitude information. Taking the sixth pixel of Table 4 for example, its amplitude value 367.30 is slightly below the threshold 367.54, while its amplitude dispersion 0.2757 is much lesser than 0.35. Comprehensively considering, it’s a better choice to classify it as a PS.

    Table 4.  The amplitude (A) and amplitude dispersion (AD) values of some pixels in PS-FAIDTA/AIDTA.
    Index Values
    AD 0.2751 0.2729 0.2867 0.2652 0.2757 0.2465 0.2514 0.2961
    A 365.87 360.84 365.40 361.51 367.30 366.04 363.89 363.83

     | Show Table
    DownLoad: CSV

    Figure 7 shows the distribution of membership degrees of those PS-FAIDTA/AIDTA under different membership degree thresholds. It can be found that the membership degrees are relatively greater compared to the minimum membership degree of PS-AIDTA 0.9263 on the whole and the lesser the membership degree threshold is, the more PS-FAIDTA/AIDTA can be selected as. The twenty-second point, whose amplitude value is 366.04, amplitude dispersion value is 0.2465 and membership degree is greater than 0.95, indicates the illogical judgment of excluding all of them by AIDTA.

    Figure 7.  The membership degrees of the PS-FAIDTA/AIDTA under different membership degree thresholds.

    From Cases 1 and 2, it can be found that AIDTA method always lose some key pixels which could be identified as PS for application while FAIDTA can make up this shortcoming. In other words, combing the two methods can usually lead to a nicer result for PS identification. Firstly, traditional method AIDTA is used to get some PS with clear characteristics. Then, some supplement pixels with high membership degree selected by FAIDTA method are obtained. As a consequence, it can increase the density of practical PS. Sub-graphs (a), (b) in Figures 2 and 6 show the selection results for the combining method. Moreover, from Tables 1 and 3, we can find that the new method plays a greater role in case 1 (the percent of n-FAIDTA/AIDTA is higher). In other words, the extent of the increase is different for different study areas. This depends on the land cover types. The vaguer the difference between land covers, the greater the extent of increase and the more sensitive the selection result. For instance, if a study area is covered with various reflective land covers, more reasonable PS can be selected based on the FAIDTA method and a slight change of the threshold will lead to great changes of the selection result. However, an opposite result can be shown in areas covered only by crops and some man-made corner reflectors. Therefore, the FAIDTA method can further improve the PS technology. In addition, it is clearly that the quality results of the selected PS can be also obtained, which can be used to evaluate the degree of the reliability for the selected PS. This is a remarkable advantage compared with traditional methods.

    Based on the fuzzy sets theory and features of PS pixels, a new concept of fuzzy PS was proposed and a fuzzy identification method was designed. Taking the amplitude information dual-threshold approach (AIDTA) as an example, the characters of PS were quantitatively described by membership degrees of amplitude and amplitude dispersion. Accordingly, a fuzzy AIDTA (FAIDTA) method was formed. Two groups of real data indicated that FAIDTA could recognize more effective and reliable permanent scatterers so as to compensate for traditional method AIDTA. In addition, quality results of the selected PS also could be obtained from the membership degrees, which could be used for quantitative evaluation of the selected points in time-series analyzing step of PSInSAR technology. That is to say, the new method not only can provide more practical PS to promote the time series analysis reliability for ground subsidence, but also can present quality value of each PS which can be used to evaluate the result of deformation based on PSInSAR technology. Generally, the proposed method can be extended based on other traditional PS identification methods, which is an effective feature fusion tool for PS identification and can improve the PSInSAR technology. Furthermore, how to apply the selected PS to subsidence monitoring is our future work. Additionally, heuristic and PSO algorithms can also be used to solve this problem, which are our other further research.

    This work was supported by Science Foundation of Jiangsu Province (No: BK20170248) and National Natural Science Foundation of China (No: 51804297).

    The authors declare no conflict of interest.

    [1] S. CHOWDHURY, J. Z HANG, A. M ESSAC, AND L. CASTILLO, Exploring key factors influencing optimal farm design using mixed-discrete particle swarm optimization, in 13th AIAA/ISSMO Multidisciplinary Analysis Optimization Conference, 2010.
    [2] S. CHOWDHURY, J. Z HANG, A. M ESSAC, AND L. CASTILLO, Unrestricted wind farm layout optimization (uwflo): Investigating key factors influencing the maximum power generation, Renewable Energy, 38 (2012), pp. 16–30.
    [3] C.-R. CHU AND P.-H. C HIANG, Turbulence effects on the wake flow and power production of a horizontal-axis wind turbine, Journal of Wind Engineering and Industrial Aerodynamics, 124(2014), pp. 82–89.
    [4] J. F ENG AND W. Z. SHEN, Wind farm layout optimization in complex terrain: A preliminary study on a gaussian hill, in Journal of Physics: Conference Series, vol. 524, IOP Publishing, 2014, p. 012146.
    [5] J. S. G ONZÁLEZ, M. B URGOS PAYAN, AND J. M. R IQUELME -S ANTOS, Optimization of wind farm turbine layout including decision making under risk, Systems Journal, IEEE, 6 (2012), pp. 94–102.
    [6] J. S. G ONZÁLEZ, A. G. G. R ODRIGUEZ, J. C. M ORA, J. R. S ANTOS, AND M. B. PAYAN, Optimization of wind farm turbines layout using an evolutive algorithm, Renewable Energy, 35 (2010), pp. 1671–1681.
    [7] S. GRADY, M. H USSAINI, AND M. M. ABDULLAH, Placement of wind turbines using genetic algorithms, Renewable energy, 30 (2005), pp. 259–270.
    [8] H. GU AND J. WANG, Irregular-shape wind farm micro-siting optimization, Energy, 57 (2013), pp. 535–544.
    [9] N. HANSEN, The cma evolution strategy: a comparing review, in Towards a new evolutionary computation, Springer, 2006, pp. 75–102.
    [10] H.-S. HUANG, Distributed genetic algorithm for optimization of wind farm annual profits, in Intelligent Systems Applications to Power Systems, 2007. ISAP 2007. International Conference on, IEEE, 2007, pp. 1–6.
    [11] T. JÄGER, R. M CKENNA, AND W. FICHTNER, Onshore wind energy in baden-württemberg: a bottom-up economic assessment of the socio-technical potential.
    [12] N. O. JENSEN, A note on wind generator interaction, 1983.
    [13] I. K ATIC, J. H ØJSTRUP , AND N. JENSEN, A simple model for cluster efficiency, in European Wind Energy Association Conference and Exhibition, 1986, pp. 407–410.
    [14] K. KOUTROUMPAS, E. C INQUEMANI , P. KOURETAS, AND J. LYGEROS, Parameter identification for stochastic hybrid systems using randomized optimization: A case study on subtilin production by bacillus subtilis, Nonlinear Analysis: Hybrid Systems, 2 (2008), pp. 786–802.
    [15] K. KOUTROUMPAS, E. C INQUEMANI , AND J. LYGEROS, Randomized optimization methods in parameter identification for biochemical network models, Proceedings of FOSBE 2007, (2007).
    [16] X. LI , J. WANG, AND X. ZHANG, Equilateral-triangle mesh for optimal micrositing of wind farms, in Proceedings of the 14th WSEAS international conference on Computers, Corfu Island, Greece, 2010, pp. 23–25.
    [17] J. E. G. M ARTÍNEZ, Layout optimisation of offshore wind farms with realistic constraints and options, Master’s thesis, Delft University of Technology, 2014. Available:http://www.lr.tudelft.nl/fileadmin/Faculteit/LR/Organisatie/Afdelingen_en_Leerstoelen/Afdeling_AEWE/Wind_Energy/Education/Masters_Projects/Finished_Master_projects/doc/Julian_Gonzalez_r.pdf.
    [18] J. C. M ORA, J. M. C. B ARÓN, J. M. R. S ANTOS, AND M. B. PAYÁN, An evolutive algorithm for wind farm optimal design, Neurocomputing, 70 (2007), pp. 2651–2658.
    [19] G. MOSETTI, C. P OLONI, AND B. DIVIACCO, Optimization of wind turbine positioning in large windfarms by means of a genetic algorithm, Journal of Wind Engineering and Industrial Aerodynamics, 51 (1994), pp. 105–116.
    [20] B. PÉREZ, R. M ÍNGUEZ , AND R. GUANCHE, Offshore wind farm layout optimization using mathematical programming techniques, Renewable Energy, 53 (2013), pp. 389–399.
    [21] E. S. POLITIS, J. P ROSPATHOPOULOS, D. C ABEZON, K. S. H ANSEN, P. C HAVIAROPOULOS, AND R. J. BARTHELMIE, Modeling wake effects in large wind farms in complex terrain: the problem, the methods and the issues, Wind Energy, 15 (2012), pp. 161–182.
    [22] P.-E. R ÉTHORÉ, P. F UGLSANG, G. C. L ARSEN, T. B UHL, T. J. L ARSEN, H. A. M ADSEN, ET AL., Topfarm: Multi-fidelity optimization of offshore wind farm, Proceedings of the Twenty-first, (2011), pp. 516–524.
    [23] R. A. RIVAS , J. C LAUSEN, K. S. H ANSEN, AND L. E. JENSEN, Solving the turbine positioning problem for large offshore wind farms by simulated annealing, Wind Engineering, 33 (2009), pp. 287–297.
    [24] S. M. F. RODRIGUES, P. B AUER, AND J. P IERIK, Modular approach for the optimal wind turbine micro siting problem through cma-es algorithm, in Proceedings of the 15th annual conference companion on Genetic and evolutionary computation, ACM, 2013, pp. 1561–1568.
    [25] B. SAAVEDRA-M ORENO, S. S ALCEDO-S ANZ, A. PANIAGUA-T INEO, L. P RIETO, AND A. PORTILLAFIGUERAS, Seeding evolutionary algorithms with heuristics for optimal wind turbines positioning in wind farms, Renewable Energy, 36 (2011), pp. 2838–2844.
    [26] M. SAMORANI, The wind farm layout optimization problem, in Handbook of Wind Power Systems, Springer, 2013, pp. 21–38.
    [27] S. ¸SI ¸SBOT , Ö. T URGUT, M. T UNÇ, AND Ü. ÇAMDALI, Optimal positioning of wind turbines on gökçeada using multi-objective genetic algorithm, Wind Energy, 13 (2010), pp. 297–306.
    [28] M. SONG, K. C HEN, Z. H E, AND X. ZHANG, Wake flow model of wind turbine using particle simulation, Renewable energy, 41 (2012), pp. 185–190.
    [29] M. SONG, K. C HEN, Z. H E, AND X. ZHANG, Bionic optimization for micro-siting of wind farm on complex terrain, Renewable Energy, 50 (2013), pp. 551–557.
    [30] M. SONG, K. C HEN, Z. H E, AND X. ZHANG, Optimization of wind farm micro-siting for complex terrain using greedy algorithm, Energy, 67 (2014), pp. 454–459.
    [31] J. T ZANOS, K. M ARGELLOS, AND J. LYGEROS, Optimal wind turbine placement via randomized optimization techniques, in Proceedings of the 17th Power Systems Computation Conference, Stockholm, Sweden, 2011, pp. 22–26.
    [32] M. WAGNER, K. V EERAMACHANENI, F. N EUMANN, AND U.-M. OÂÂZREILLY, Optimizing the layout of 1000 wind turbines, European Wind Energy Association Annual Event, (2011), pp. 205–209.
    [33] C. A. WALFORD, Wind turbine reliability: understanding and minimizing wind turbine operation and maintenance costs, United States. Department of Energy, 2006.
    [34] C. WAN, J. WANG, G. YANG, X. L I , AND X. ZHANG, Optimal micro-siting of wind turbines by genetic algorithms based on improved wind and turbine models, in Decision and Control, 2009 held jointly with the 2009 28th Chinese Control Conference. CDC/CCC 2009. Proceedings of the 48th IEEE Conference on, IEEE, 2009, pp. 5092–5096.
    [35] C. WAN, J. WANG, G. YANG, AND X. ZHANG, Optimal micro-siting of wind farms by particle swarm optimization, in Advances in swarm intelligence, Springer, 2010, pp. 198–205.
    [36] J. WANG, X. L I , AND X. ZHANG, Genetic optimal micrositing of wind farms by equilateraltriangle mesh, INTECH Open Access Publisher, 2011.
    [37] V. YAKHOT, S. O RSZAG, S. T HANGAM, T. G ATSKI, AND C. SPEZIALE, Development of turbulence models for shear flows by a double expansion technique, Physics of Fluids A: Fluid Dynamics (1989-1993), 4 (1992), pp. 1510–1520.
    [38] M. H. ZHANG, Wind Resource Assessment and Micro-siting: Science and Engineering, Wiley, 2015.
  • This article has been cited by:

    1. Yongjiu Feng, Yilun Zhou, Yanling Chen, Pengshuo Li, Mengrong Xi, Xiaohua Tong, Automatic selection of permanent scatterers-based GCPs for refinement and reflattening in InSAR DEM generation, 2022, 15, 1753-8947, 954, 10.1080/17538947.2022.2070678
    2. Weitao Liu, Qiushuang Zheng, Lifu Pang, Weimeng Dou, Xiangxi Meng, Study of roof water inrush forecasting based on EM-FAHP two-factor model, 2021, 18, 1551-0018, 4987, 10.3934/mbe.2021254
    3. Dan Wang, Yier Lin, Liang Hong, Ce Zhang, Yajie Bai, Zhen Zhen Bi, Naeem Jan, Detection of the Driver’s Mental Workload Level in Smart and Autonomous Systems Using Physiological Signals, 2022, 2022, 1563-5147, 1, 10.1155/2022/5233257
    4. Yun Wang, Songpeng Wang, Ram Bilas Pachori, Lei Chen, 2024, Surface deformation monitoring and analysis using PS-inSAR in Bortala Mongol Autonomous Prefecture, 9781510680012, 127, 10.1117/12.3029813
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5755) PDF downloads(1154) Cited by(0)

Figures and Tables

Figures(10)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog